The Feni Island group is located at the southeastern end of the NW trending Ta-bar-Lihir-Tanga-Feni (TLTF) island chain, northeastern Papua New Guinea. This island chain is renowned for hosting alkaline magmas, geothermal activity, copper-gold mineralization and mining. The purpose of this paper is to present the geological observations of Feni Island within the context of the regional tectonic setting and evolution of the TLTF. The TLTF arc is a younger arc situated within the Greater Melanesian Arc, specifically within the New Ireland Basin bounded by New Ireland to the west and the extinct Kilinailau Trench in the east. The formation of the TLTF is attributed to post-subduction extensional tectonics relating to the opening of the Manus Basin and the magma is sourced from subduction-modified upper mantle. Normal, ex-tensional faulting is observed throughout the TLTF.

This study is describing airborne geophysical survey, carried out in Northeastern Russia (Archangelsk region) in order to aimed at geological mapping of Grib kimberlite pipe, which is known to have strong presence of remanent magnetization. Airborne data were inverted using two techniques, the magnetization vector inversion (MVI) and a QDIK inversion, a proprietary algorithm developed at the Irkutsk State University (ISTU), which allows extracting four parameters from airborne magnetic data (magnetic susceptibility, inclination, declination and Koeniksberger ratio). The inversion results were further analyzed and verified against the known geology of the Grib kimberlite. The study allowed to better map geology of Grib kimberlite using only uninvasive (remote sensing) technology

The acquisition of regional, 1 metre resolution LIDAR in the PNG Highlands combined with 3D modelling in MoveTM has revolutionized our understanding of the evolution of the Porgera gold mine and areas of new potential. The new 3D model demonstrates active pull-apart basins along a regional transfer during ongoing fold and thrust deformation. When overlain on regional aeromagnetic data, new potential is revealed. The Porgera gold mine is one of the richest in the world and lies in a wide valley at an elevation of 2800m surrounded by mountains up to 4000m high in the middle of the PNG Highlands. It lies along a major lineament, the Porgera Transfer Zone (PTZ), cutting across the orogen that is associated with a 50 km offset of the ophiolite belt. To the NW of the mine there is an extensive belt of low- to high-grade metasediments that formed in deep water during Jurassic rifting, but were metamorphosed during Eocene subduction to the north. Subduction ceased due to accretion of the Sepik Terrane which caused inversion and mild erosion of the Porgera area. The Late Oligocene onset of wrenching in the Mobile Belt to the north placed that area into extension, emplacing metamorphic core complexes, and led to regional subsidence in the Early Miocene. Collision of the margin with the Melanesian Arc in the Middle Miocene caused Late Miocene to Pliocene orogenesis creating the broad mountain belt that we see today, that is still active, as shown by the 2018 MW 7.5 earthquake. Regional mapping of the area using high resolution LIDAR in association with limited field mapping, analysis of gravity and magnetics data and drilling of 300-500m deep core-holes has allowed development of a detailed 3D structural model. The Porgera valley is constrained laterally by the Eastern Boundary and Western Boundary dextral strike-slip faults that lie along the PTZ cutting across the orogeny. However, whilst the position and nature of the faults is clear, due to bending, fracturing and offset of major synclines and anticlines, the faults do not crop-out as significant through-going features. Rather, the dextral offset of basement is soft-linked to the Tertiary limestones at surface due to intervening thick, ductile Cretaceous shales. Fold and thrust structures are ubiquitous in the Tertiary limestones, but the youngest structural features are large extensional faults, particularly around the Porgera valley. The intrusive underlying the Porgera ore-body was emplaced at 6.0+-0.3 Ma, the time of maximum compression during orogenesis. At this time the Eastern Boundary fault and Western Boundary fault were both active allowing a pull-apart basin to form locally along the PTZ at the south-eastern boundary of the Jurassic metasediments. This enabled emplacement of the intrusive. The gold-bearing fluids from the intrusive and metasediments were brought up through the Mesozoic muds to the contact with the Paleogene carbonates where the Porgera ore body was emplaced.

The study of fractures in the subsurface is very important in unconventional reservoirs since they are the main conduits for hydrocarbon flow. For this reason, a variety of equivalent medium theories have been proposed for the estimation of fracture and fluid properties within reservoir rocks. We experimentally investigated the feasibility of applying the Galvin model in fractured tight stones. For this proposal, three artificial fractured tight sandstone samples with the same background porosity (11.7% ± 1.2%) but different fracture densities of 0.00, 0.0312, and 0.0624 were manufactured. The fracture thickness was 0.06 mm and the fracture diameter was 3 mm in all the fractured samples. Ultrasonic P- and S-wave velocities were measured at 0.5 MHz in a laboratory setting in dry and water-saturated conditions in directions at 0°, 45°, and 90° to the fracture normal. The results were compared with theoretical predictions of the Galvin model. The comparison showed that model predictions significantly underestimated P- and S- wave velocities as well as P-wave anisotropy in water-saturated conditions, but overestimated P-wave anisotropy in dry conditions. By analyzing the differences between the measured results and theoretical predictions, we modified the Galvin model by adding the squirt flow mechanism to it and used the Thomsen model to obtain the elastic moduli in high- and low-frequency limits. The modified model predictions showed good fits with the measured results.

In this study, UAV (Unmanned Aerial Vehicle) photogrammetry was used to monitor the ground displacement on the slope below a coal waste-rock stockyard, and to investigate the role of rainfall on this displacement. The study area is a mountainous region in Korea, where coal mining continues, and coal waste-rock is stored on the slopes of the mountain. In this area, ma-terial removal work was undertaken to reduce the steepness of the slope and driving forces in order to prevent disasters, but the strategy requires continuous monitoring to confirm the stabilization of the slope. For slope monitoring, a total of six UAV photogrammetry campaigns were conducted between April 2019 and August 2020. As a result of data processing, an orthoimage and DSM (Digital Surface Model) were generated. The ground displacement was estimated through four ways: slope extraction, displacement area evaluation, horizontal displacement, and vertical displacement analysis. During the study period, the maximum vertical displacement was 3.3 m and the maxi-mum horizontal displacement was 3.5 m. The horizontal displacement was effectively evaluated through the measurement of the movement of the drainage system using orthoimages that captured with a periodic survey. The effects of rainfall on ground displacement were also investigated. A very clear linear relationship between rainfall intensity higher than 20mm/d and ground displacement was identified. Accumulated amount of rainfall also showed good correlation with slope movement, but the frequency of rainfall intensity below 20 mm/d showed relatively low correlation with ground dis-placements.

This paper deals with the quantitative analysis of measured fracture-induced electromagnetic radiation (FEMR) in the vicinity of the Dead Sea Transform using the Angel-M1 instrument, which enables the recording of FEMR signals in a 3D manner. Results show both the possibility of estimating the sizes of micro-fractures that are the sources of the radiation, and of assessing the direction of their location to the measuring device and the range of magnitude Mw of the impending "events" (EQ's) associated with the present FEMR measurements. Moreover, the relation between the measured FEMR activity (the number of FEMR hits per unit time) and the FEMR event magnitudes shows consistency with the Gutenberg-Richter relationship for the region. Such measurements could therefore constitute a preliminary 'field reinforcement' towards an EMR valid method for a real earthquake forecast, which would provide much earlier warnings than the seismic ones. Obviously, the present observed FEMR measurements can only be used as an assessment of the stress concentration and micro-fracturing in the region since they relate to the very initial nucleation phase of a "virtual" (extremely low Richter scale value) earthquake, but they do provide the necessary feasibility test for a prediction method since all lab-measured FEMR features were confirmed also in the field.

The determination and classification of seafloor sediment types are crucial for the exploitation of marine resources, construction of marine engineering, and maintenance of marine ecological en-vironments. Automatic classification of seafloor sediment based on acoustic telemetry data is an important method to quickly understand the type of a large range of sediment. Currently, most studies on sediment classification are based on multi-beam backscattering intensity data, which is a relatively single data type. Besides, the low-dimensional details of standard U-Net gradually weaken in the propagation process, limiting the accuracy of sediment classification. Therefore, this study proposes an automatic classification method of seafloor sediment types based on an im-proved U-Net and K-means clustering algorithm, using multi-beam water depth, sub-bottom profile, and sample test data in the Northern Slope of the South China Sea. Six sediment types, including gravelly muddy sand, sand, silty sand, less muddy silt, muddy silt, and silty mud, were identified in the study area. Additionally, the study area was divided into four sedimentary en-vironment zones, including the shelf sedimentary area, the upper shelf slope break sedimentary area, the lower shelf slope break mixed sedimentary area, and the slope deposition area, based on the results of sediment classification and geological background. The sedimentary environment zones were found to be distributed along the trend of the shelf slope break line. The results of this study not only provide an important supplement to the existing classification methods of seafloor sediments but also contribute to the understanding of the sedimentary environment and process of the Northern Slope of the South China Sea.

Comparison with conventional reservoirs, the prediction of pressure in ultra-deep tight sandstone formations is difficult. The prediction of seismic pressure is more challenging than well logging pressure prediction. The main methods for seismic pressure prediction include equivalent depth method, Eaton method, Fillippone formula, and its modified versions. Among them, the Eaton method is widely used and has good effectiveness. However, this method relies on difficult-to-obtain normal compaction trend line, which leads to low prediction accuracy in space. To address this issue, a method combining the Eaton method and Collocated cokriging is proposed. The Eaton formula is used to predict formation pressure at the well, with compressional wave velocity as the covariate for predicting the main variable, formation pressure. By simulating the shear wave velocity based on seismic compressional wave velocity, the influence of various parameters on the prediction results is analyzed, and the accuracy of this method is verified by comparing it with other methods. The proposed method is then applied to predict formation pressure in the ultra-deep formations of the Junggar Basin. The simulation results show that the Collocated cokriging method achieves higher planar accuracy and better matches the experimental expectations in terms of prediction results. The application results also demonstrate the scientific effectiveness of the combined method, which has achieved good results in practical production applications.

The sinking of cavities in urban areas is a very important geological risk to which many towns in Puglia are exposed. The study, mapping, geological and speleological description, classification and cataloging of forms and types of cavities is of fundamental importance for each urban center because it is linked to past local anthropic and natural processes, unique to each site. All this could lead to an enhancement of existing underground cavities in urban areas through conservation and continuous monitoring. Unfortunately in many cases these underground cavities have been used as landfills and subsequently abandoned. One of these cavities at the end of March 2007 collapsed in the inhabited center of Gallipoli causing damage to the structures and fortunately not to human lives. In the area relating to the collapsed cavity, a series of geophysical investigations were undertaken with the use of ground penetrating radar methodology in an attempt to delimit the area of collapse and understand the possible types of interventions for restoring the area. In the same area the measures were repeated 16 years later (December 2022) due to another collapse. The results show a much more widespread danger situation.

Abstract: This paper summarises an early and successful piece of citizen science, performed within The University of Queensland Seismograph Stations (UQSS) observatory, in cooperation with colleagues at CSIRO. It was designed to en-courage young STEM students from Brisbane high schools to engage in “real” re-search, back in 1995. Having completed the project report, their (analog) results sat in a cupboard until the report was dusted off and the project was re-analysed in 2022 by an honors student, considering timely climate change applications for the study. This is a time when science is changing considerably from analog to digital medium and operational methods. The original project was called Earthquake generated T phases on BRS Seismograph (Brisbane, Q’ld)- a predictor for Tasman Sea Tsunamis? [1] Fortunately, seismology is a very collaborative field. The research question has since changed. There is a lot of data analysis involved in the science of recording earthquake signals, with auxiliary definitive catalogues, observers logbooks, housing of the recordings themselves (analog and digital) and the software mediums that change over time. In other words, a lot of headaches can be encountered in longitudinal data collection study such as this. The citizen science students used a pre-prepared decadal collection (1980-90) derived from the BRS observatory data catalogue. BRS is part of the global World-Wide Seismograph Station Network (WWSSN). Currently in Australia, university earth science observatories have diminished, and in their place, public seismic networks (PSN) have evolved, either in back-yard sheds or school science labs. The level of expertise required fits the role of advancing citizen science for a real science advantage. This is a topical citizen disaster preparedness action area for today’s climate emergency, which is threatening the globe and all lifeforms. Citizen engagement or mobilisation is an essential preparedness strategy. Keywords: citizen science 1; seismology 2; earth observation 3; climate change 4; education for sustainable development 5; disaster preparedness 6; disaster alerts 7

In order to study the surface settlement law of subway tunnel engineering site overlying soft soil layer under the action of earthquake, ABAQUS finite element analysis software was used to study the surface settlement law of subway tunnel construction site overlying soft soil layer under the action of three different seismic waves, including El-Centro wave, Kobe wave, and artificial seismic wave with exceedance probability of 10%, based on the Dalian Metro Line 5 with overlying soft soil layer as the engineering background, and the law was analyzed theoretically. Then, the surface settlement law of tunnel top was studied based on El-Centro wave and soft soil parameters, and the sensitivity analysis of soft soil parameters on surface settlement of tunnel top was carried out by orthogonal test method. The results show that under earthquake action, the settlement of the strata within a certain range of the tunnel roof is significantly greater than that of the surrounding strata, forming a settlement trough with a width of about 8 to 20 meters. The width of the settlement trough under the El-Centro wave is the largest, about 19.6 meters, surpassing the artificially synthesized seismic waves with a probability of 10%, which is about 15.6 meters, while the width of the settlement trough under the Kobe wave is the smallest, about 8.5 meters. The ground surface within a range of about 20 meters above the tunnel roof is most strongly affected by the seismic waves and the special lithology of the overlying soft soil layer, and the settlement is the largest. The settlement law of the settlement trough in the overlying strata of the tunnel conforms to the Peck formula. Increasing the elastic modulus of the silty soil can reduce the settlement of the ground surface above the tunnel roof, increasing the Poisson's ratio of the silty soil will increase the settlement of the ground surface above the tunnel roof, increasing the cohesive force of the silty soil to 20 kPa will basically stabilize the settlement of the ground surface above the tunnel roof, and increasing the internal friction angle of the silty soil will basically not change the settlement of the ground surface above the tunnel roof. The sensitivity of the soft soil parameters to the settlement of the ground surface above the tunnel roof is in the order of Poisson's ratio, elastic modulus, cohesive force, and internal friction angle.

Shallow landslides pose serious threats to human existence and economic development, especially in the Himalayan areas. Landslide susceptibility mapping (LSM) is a proven way for minimizing the hazard and risk of landslides. Modeling as an essential step, various algorithms have been applied to LSM. In this study, information value (IV) and logistic regression (LR) were selected as representatives of the conventional algorithms, categorical boosting (CatBoost) and conventional neural networks (CNN) as the advanced algorithms, for LSM in Yadong county, and their performance was compared. To begin with, 496 historical landslide events were compiled into a landslide inventory map, followed by a list of 11 conditioning factors, forming a data set. Secondly, the data set was randomly divided into two parts, 80% of which was used for modeling and 20% for validation. Finally, the area under the curve (AUC) and statistical metrics were applied to validate and compare the performance of the models. The results showed that the CNN model performed the best (AUC 0.974 and accuracy=93.3%), while the LR model performed the worst (AUC 0.974 and accuracy=93.3%) and CatBoost model performed better (AUC 0.974 and accuracy=93.3%). Besides, the LSM constructed by the CNN model did a more reasonable prediction of the distribution of susceptible areas. As for feature selection, did a more detailed analysis of conditioning factors but the results were uncertain. The result analyzed by GI may be more reliable but fluctuates with the amount of data. The conclusion reveals that the accuracy of LSM can be further improved with the advancement of algorithms, by determining more representative features, which serve as a more effective guide for land use planning in the study area or other highlands where landslides are frequent.

Future climate change is expected to affect components of the hydrological cycle. However geologic factors can sometimes be more important than the climate with regarding to change hydrological processes. The Keriya River is a large glacier-fed river in the western Tibet Plateau, its upper reaches flow across active fault zones, which have experienced four medium-sized earthquakes within four months 1974. These active faults have made difficult to identify normal meltwater changes using a snow-ice melt model. Because of the seismic events, the river has experienced substantial changes in their drainage systems and complex hydrological processes. Both the Mann-Kendall trend test and a correlation analysis were applied to detect the monthly patterns of stream flows, the change points and the correlations among both summer discharge and air temperature before and after seismic events. We found significant correlations between monthly air temperature and meltwater flow in the summer June during the pre-earthquake of 1958-1974, but none in July and August in the post-earthquake of 1975-1980, and we found no correlation following the earthquake. The correlations were the inverse by losing the former and setting up the latter alternately since the multiple earthquakes. Clearly, the earthquakes and changes in the faults can greatly induce hydrological processes.

In China, the majority of mountainous regions are characterized by complex topography and a delicate, sensitive geological environment. Coupled with a generally underdeveloped infrastruc-ture and numerous unreasonable human engineering activities, these regions are often highly susceptible to geological disasters. Geological hazards can cause significant damage to human lives and property, impeding the development of mountainous areas. Consequently, researching the assessment of geological hazard vulnerability is crucial for disaster prevention, emergency management, and economic development in these regions. This study focuses on Xuanwei City and selects eight factors for evaluation, including elevation, gradient, slope aspect, normalized vegetation index, stratigraphic lithology, distance from faults, distance from rivers, and distance from roads. These factors are chosen based on a comprehensive analysis of the spatial and tem-poral distribution of geological hazards and disaster incubation conditions. Two paired models, the deterministic coefficient model + logistic regression model (CF+LR) and the information quan-tity model + logistic regression model (I+LR), were employed to quantitatively assess the study ar-ea. The accuracy of these models was evaluated using ROC curves and AUC values. The results indicate that: (1) The AUC values for the CF+LR and I+LR coupled models are 0.799 and 0.772, respectively, demonstrating that both models can objectively and reliably assess the vulnerability to geological hazards in the study area; (2) Based on the CF+LR model calculations, the geological hazard susceptibility of Xuanwei City can be categorized into four zones: extremely high suscepti-bility (6.09%), high susceptibility (31.08%), medium susceptibility (32.26%), and low susceptibility (30.57%); (3) The CF+LR model more accurately represents the evaluation results and offers a strong reference value.

(1) Background: The intermediate-depth seismicity in the Vrancea region (Romania) is characterized by localized and persistent earthquake activity that culminates about two or three times in a century with the occurrence of a large event (M 6.5). Here we have revisited some important seismicity characteristics, using earthquake catalog data spanning two different time periods: 2005 – 2013 and 1960 – 2000. (2) Methods: we have determined the b-value of the frequency-magnitude distribution of earthquakes, using a maximum likelihood procedure, and estimated the parameter to quantify anomalous seismicity rate decreases and increases. (3) Results: by using data from the first period, we have confirmed the existence of a decreased b-value in the deepest part of the seismogenic zone; by using data from the second period, we have statistically confirmed the seismic quiescence that preceded the occurrence of the 1977 M7.4 Vrancea earthquake. (4) Conclusions: the decreased b-value has been interpreted either in terms of an increased lithostatic stress with depth or as an indicator of the depth range where the next major Vrancea earthquake may occur. The time variation of the seismicity parameter may reveal anomalous seismic quiescence and increased earthquake rates that may precede the occurrence of large earthquakes.

Nepal is one of the most seismically active regions in the world, as highlighted by the recent devastating 2015, Mw7.8 Gorkha earthquake and a robust assessment of seismic hazard is paramount for the design of earthquake-resistant structures. In this study we present a new probabilistic seismic hazard assessment (PSHA) model for Nepal. We considered data and findings from recent scientific publications, which allowed us to develop a unified homogenized seismicity catalogue, propose alternative Seismic Source Characterization (SSC) models including up-to-date parameters of major thrust faults like Main Frontal Thrust – MFT and Main Boundary Thrust – MBT, while also considering existing SSC models and various hazard modelling strategies within a logic tree framework. The sensitivity analyses show the hazard levels are generally higher for SSC models integrating the major thrust faults, followed by homogenous volume sources and smoothed seismicity approach. The hazard maps covering entire Nepal is presented as well as the Uniform Hazard Spectra (UHS) for 5 selected locations (Kathmandu, Pokhara, Biratnagar, Nepalganj and Dipayal) at return periods of 475- and 2,475-years considering Vs,30=760 m/s. The results obtained are generally consistent with most recent studies. However, a notable variability in hazard levels and several discrepancies with respect to the Nepal Building Code NBC105: 2020 [1] and Global Hazard Model, GEM [2] are noted and possible causes are discussed.

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.

According to the structural features, sedimentary evolution, types of sedimentary facies, source rocks features, diagenesis evolution, reservoir features and hydrocarbon formation mechanism, exploration status and hydrocargon resource potential, this article analyzed the major control fac-tors and hydrocarbon distribution rules of the lithologic reservoirs in continental fault basin. The research shows that three major control factors and one coupled factor (fractures act as a tie) in-fluence hydrocarbon formation of lithologic traps in Nantun Formation of Tanan Sag, which in-clude sand body types, effective source rocks and effective reservoirs. With the increase in depth, enough hydrocarbon generate in source rock under thermal evolution. Hydrocarbon generation and expulsion are more intensive when it comes to the depth threshold and critical conditions of hydrocarbon supplying are met. Traps that surrounded or contacted by source rock, or communi-cated by fault are qualified to form reservoirs. As buried depth increases, hydrocarbon genera-tion-expulsion intensity grows and the trap is more petroliferous. Hydrocarbon accumulation and reservoir formation are also controlled by sandbody accumulation conditions. When it meets the critical conditions of hydrocarbon generation and concrete oil and gas are charged in, better phys-ical properties of sandbody always indicate more hydrocarbon accumulation in the trap. Alloca-tion of sand type, effective source rock and effective reservoir is optimized under the effect of fractures, and coupled hydrocarbon reservoir with these three elements.

Several coastal barriers experienced significant erosion and change in shape throughout the Mediterranean coasts over the past decades, and the issue has become of increasing concern for scientists and policymakers. With reference to a case study and by meteorological and geomorphological investigations, this note aims to infer recent shoreline dynamics and define coastal phenomena shaping the beach-dune barrier of Cesine Lagoon (South Italy), to finally point out the weather-driven processes. A challenge for this research was to obtain significant results from publicly available sources and simple and inexpensive methods. The retreat of the dune toe, the accretion of washover fans, and the formation of gravel beaches have been related to the analyzed wind-wave conditions. Net longshore transport was found in accordance with the direction of the more intense winds. The role of extreme events in the shaping of coastal landforms is yet to be established even if they greatly increase the vulnerability to flooding of the study area. The results achieved so far are starting points for further data collection and analysis in the perspective of assessing the impact of climate changes on the lagoon barrier.

The correlations of seismic noise, recorded during 1997-2022 at 229 stations around the world are investigated. The property of waveforms, known as the Donoho-Johnston threshold, is analyzed, which separates the absolute values of the orthogonal wavelet coefficients into "small" and "large". Network of 50 reference points is considered. For each reference point, a time series with a time step of 1 day is calculated as the median of the values of the DJ index from the 5 nearest stations. For each reference point, the correlations between the values of the DJ index at other reference points are calculated in a sliding time window of 365 days. Average values of the DJ index decrease, which is interpreted as indicators of global seismic hazard increasing. The maximum distances between reference points with strong correlations show an explosive increase after 2011. Bursts in the mean coherences between the length of day (LOD) and the values of the DJ index at the reference points precede the release of seismic energy with a delay of about 530 days. High amplitude of the DJ index response to LOD for the times 2020-2022 was identified, which may forecast a strong earthquake in the second half of 2023.

As semi-airborne mineral exploration has limited budgets, it is critical to design experimental procedures that generate data maximizing the desired information We investigate the effects of transmitter-receiver geometries for a variety of anomalies and semi-airborne layouts. Our simulations indicate that 200 m flight line spacing and 100 m point distance are the optimal trade-off between coverage and survey progress for various targets. Based on the target size and distance between the transmitter and the target, the transmitter length should be at least equal to the length of the target and for more than 1 km, at least two or three times the target size. Likewise important are the location and direction of the transmitter cables, which can have great impact on the result of inversion and should be parallel to the target strike. By using more than one transmitter, better results are obtained. If the strike of the target is known, transmitters should be parallel to each other and if not, it is better to use perpendicular transmitters. Results showed that the optimal distance between transmitters is 3 km. Our simulations show that it is even possible to recover targets just below the transmitter in corresponding areas of masked data.

Northeastern margin is a natural experimental field for studying crustal extrusion and expansion mechanism. Accurate crustal deformation pattern is a key point in the analysis of regional deformation mechanism and seismic hazard research and judgment. In this paper, the present-day GPS velocity field on the northeastern margin of the Tibetan Plateau was obtained from encrypted GPS observations around the Haiyuan-Liupanshan fault zone, combined with GPS observations on the northeastern margin of the Tibetan Plateau from 2010 to 2020. Firstly, we divided the study area into three relatively independent blocks: ORDOS block, Alxa block and Lanzhou block; secondly, the accurate fault distribution of the Haiyuan-Liupanshan fault zone was taken into account to obtain the optimal inversion model; finally, using the block and fault back-slip dislocation model, the inversion got the slip rate distribution, locking depth and slip deficit rate of each fault. The results indicate that the Laohushan Fault and Haiyuan Fault are dominated by left-lateral strike-slip, while the Liupanshan Fault is dominated by thrust dip-slip, and the Guguan-Baoji Fault has both left-lateral strike-slip and thrust dip-slip components. The maximum locking depths of the Laohushan Fault, Haiyuan Fault, Liupanshan Fault and Guguan-Baoji Fault are 5 km, 13 km, 15 km and 10 km respectively, and the locking of Haiyuan Fault is strong in the middle section and weak in the eastern and western section. The Haiyuan Fault is still in the post-earthquake stress adjustment stage. The slip deficit rate decays from 3.6 mm/yr to 1.8 mm/yr from west to east along the fault zone. Combined with geological and historical seismic data, the results suggest that the mid-long term seismic risk in the Liupanshan Fault is high.

Coastal flooding exacerbated by climate change is recognised as a major global threat which is expected to impact more than a quarter of people currently residing in Pacific Islands countries. While most research in the last decade has focused on understanding the dynamics and impacts of future coastal flooding from extreme sea levels, relative sea level rise (RSLR) effects on tsunami hazards are not well understood. Far-field or distant-sourced tsunamis tend to have relatively lower impacts in Pacific Island states compared with locally sourced events, but there is limited understanding on how the impact of far-field tsunamis changes over time due to RSLR. Using the hydrodynamics software BG-Flood, we modelled the Tōhoku tsunami from propagation to inundation in Samoa under incremental SLR to examine the effects that RSLR has on changing the exposure of the built environment (e.g., buildings) to a far-field tsunami. Outputs of maximum tsunami inundation and flow depth intensities which incorporate incremental SLR were then combined with digital representations of buildings and depth-damage functions in the RiskScape multi-hazard risk modelling software to assess the changes in building exposure over time. Results indicate that present day buildings exposure in Samoa to a Tōhoku-oki type far-field tsunami will increase by approx. 600% with 1 m RSLR by 2080–2130, and approx. 2,350% with 2 m RSLR by 2130–2140. These findings provide a useful baseline for tsunami hazard risk assessment under changing sea level conditions in analogous island environments.

Geological structures such as joints and faults in rock mass have significant influence on open pit mining. Hence, it is critical to develop a understanding of dynamic joint behaviour under blasting loading. This in turn can provide both theoretical and practical guidance to improve blasting rock fragmentation and associated bucket excavating efficiency. In this paper, delayed blasting on the highwall bench at Baiyunebo open-pit mine was used as an example, a nonlinear joint blasting model was also constructed. By simplifying the blasting wave propagation velocity, the P-wave normal incidence to the joint was obtained. The peak vibration velocity was 0.33m/s at 3.0s. The S-wave reflected by the joint interface was first reflected backward and then forward, which the peak vibration velocity was 0.027 m/s. By combining the relevant stress and displacement theories of type I and II cracks, the equipotential diagrams of the stress and displacement field with the vibration velocity of the particle were obtained. σx was positive in the direction of 0~330° and subjected to tensile stress, whereas σy was positive in the direction of 0~180° and under tensile stress. The longitudinal σy along joint was low in compressive stress distribution area and did not affect surrounding rock at the time. The stress concentration appeared in the lower right corner. Based on the continuous behavior of stress wave in joints, the asymmetry and continuous changes reflected in the whole process could not be analysed by the contour diagram. Hence, ANSYS used to analyse distribution of the stress field. The intensity of the shock wave after detonation was greater than that of the rock strength. Subsequently, the sub-layer shock wave supplemented the shock wave energy that was not enough to break the rock and induced further cracking. This was able to be visualized by the degree of color change post-processing. It was concluded that with the attenuation of the detonation wave energy, the stress exhibited a decreasing trend in this process. According to distribution of the peak effective stress, it was found that the peak value first increased to 10-12 MPa and then showed a downward trend. Overall, the results were validated against the finite element simulation and mathematical analysis.

The Rb/Sr ratio of lake sediments has been widely adopted as an indicator of weathering intensity in studies of past climate change, but the geochemical significance of this ratio varies with timescale. Here we present Rb/Sr data for the last 300 years for sediments collected from Chaonaqiu Lake in the Liupan Mountains of the western Chinese Loess Plateau as a decadal-scale record of weathering intensity. To validate the application of this weathering proxy, we correlated the record with those of other major elements, rock-forming minerals, and paleoclimatic proxies. We found that Rb/Sr ratios are influenced mainly by Sr activity within the lake catchment (where Sr is likely sourced from albite). In addition, higher (lower) Rb/Sr ratios of bulk sediments from Chaonaqiu Lake are correlated with lower (higher) fractions of terrigenous detritus (SiO2, Ti, K2O, Al2O3, and Na2O). Lake bulk sediment Rb/Sr ratios are dominated by the input of terrigenous detritus over decadal timescales. Our data show that physical and chemical weathering in Chaonaqiu Lake catchment have opposing influences on bulk sediment Rb/Sr ratios, competing to dominate these ratios of lake sediments over different timescales, with ratios reflecting the relative importance of the two types of weathering. We infer that the Rb/Sr ratios of Chaonaqiu Lake sediments reflect climatic signals of precipitation and temperature over long timescales and of precipitation over short timescales.

During disaster response, clouds or darkness can prevent the use of optical images for detecting consequences of natural disasters, including landslides. In these situations, radar images can be used to detect changes more rapidly. However, Synthetic Aperture Radar (SAR) backscatter intensity images are underutilized for landslide detection. Unfortunately, there remains a lack of understanding about how to interpret landslide signatures in SAR imagery. In this study, we investigate how the morphometric features and material properties of landslides, and preexisting land cover, control their expression in SAR backscatter intensity change images. Trends in the spatial and temporal signatures of over 1000 landslides in 30 diverse case studies are investigated, using multi-temporal composites and dense time-series of Sentinel-1 C-band SAR backscatter intensity data. The results show that the orientation of landslide surfaces relative to the sensor, pre-existing land cover, and the roughness of the landslide surface, determine whether landslides will produce an increase or decrease in backscatter intensity values. In certain cases, we can identify morphometric features of landslides (e.g. scarps, transit zone, deposits, ponding) and material properties. Generally, we see that landslides appear most clearly with a strong increase in intensity when they occur in herbaceous vegetation or non-vegetated ground surfaces, due to an increase in surface roughness. While in forested or densely vegetated areas, landslides produce a more complex signature with both decreases due to radar shadow and vegetation removal, and an adjacent edge of increased intensity due to double bounce and direct return from vertical tree trunks and convex edges. In most cases, rough deposits produce an increase in intensity, while smooth deposits (e.g. from mudslides) exhibit specular reflection, and thus show decreased values. Landslides are less visible in cases with pre-event very rough ground, or mixed vegetation conditions. The conceptual model developed can aid interpretation of landslides in SAR imagery, and provide domain knowledge needed to train models for automatic landslide detection.

Since the early 1990s, the European (ESA) and Italian (ASI) space agencies have managed and distributed a huge amount of satellite-recorded SAR data to the research community and private industries. Moreover, the availability of advanced cloud computing services implementing different multi-temporal SAR interferometry techniques allows the generation of deformation time series from massive SAR images. We exploit the information provided by a large PS dataset to determine the temporal trend of ground deformation and the relative deformation rate with millimetric accuracy to analyse the spatial and temporal distribution of land subsidence induced by water pumping from a deep confined aquifer in the Northern Valle Umbra Basin (Central Italy), exploiting 24 years of Permanent Scatterers – Interferometric SAR data archives. The SAR images were acquired between 1992 and 2016 by satellite ERS1/2, ENVISAT and Sentinel 1 ESA missions, and the COSMO-SkyMed ASI mission). We observe ground velocities and deformation geometries between 1992 and 2016, with displacements of more than 70 cm and velocities of up to 55 mm/yr and the results suggest that the shape and position of the surface ground displacement are controlled by the fault activity hidden under the valley deposits.

We developed and experimented at sea an innovative fiber optic hydrophone (FOH) capable of operating at very low frequency of the seismic spectrum and for detecting small magnitude earthquakes. The FOH exploits an optical fiber coil wrapped around a sensitive mandrel in a Michelson interferometric configuration. The FOH operated for about seven days at a water depth of 40 m, in the Campi Flegrei volcanic area (Southern Italy), a few meters from a well calibrated PZT hydrophone used as a reference. Thirty-three local earthquakes occurred during the simultaneous operation of the two hydrophones, allowing a straightforward comparison of the recordings. The local earthquakes occurred at an epicentral distance shorter than 2.5 km from the site of recording, and were estimated to be in the range of magnitude from -0.8 to 2.7. The analysis of the recorded earthquake waveforms in the frequency and time domain allowed retrieving the response function of the FOH in the frequency range from 5 to 70 Hz. The FOH responsivity in terms of acoustic pressure reaches about 230 nm/Pa and is flat in the studied frequency range. The high quality of the FOH recordings make this equipment suitable for applications addressing submarine volcanic activity and the background seismicity of active faults in the ocean.

The Roman villa of Pisões (Beja, Portugal), was part of the Lusitanian colony of Pax Iulia. This place stands out for the predominance of the water element in several structures of the villa, highlighting the balneum and the large natatio, one of the largest known in Roman Hispania. The records of the initial excavations that took place since 1967 do not allow the establishment of clear functionalities of the villa. The University of Évora, owner of the site, conceived an action plan for the requalification and enhancement of the archaeological site. One of the tasks aims to investigate using Applied Geophysics. This work analyses the landscape directly related to the villa, given that it is in the flooded area of a river, with a Roman containment dam. It is uncertain whether the water supply comes from this structure or other nearby springs. The use of ground-penetrating radar, combined with unnamed aerial vehicles, all integrated in a geographic information system, allows us to know the location of underground water connections and create a topographic model with high resolution. Considering all the information, we propose a model for the water transport inside the villa and estimate the location of the water supply.

The 1884 Andalusia Earthquake, with an estimated Magnitude between 6.2 and 6.7, is one of the most destructive events that shook the Iberian Peninsula, causing around 1200 casualties. Ac-cording to paleoseismology studies and intensity maps, the earthquake source relates to the nor-mal Ventas de Zafarraya Fault (Granada, Spain). Diverse studies registered and later analyzed hydrological effects, such as landslides, rockfalls, soil liquefaction, all-around surge and loss of springs, alterations in the phreatic level, discharge in springs and brooks, and well levels, along with changes in water properties. Further insight into these phenomena found an interplay be-tween hydro-geomechanical processes and crust surface deformations, conditions, and properties. This study focuses on analyzing and simulating the features involved in the major 1884 event and aims at elucidating the mechanisms concerning the mentioned effects. This ex-post analysis builds on the qualitative effects and visible alterations registered by historical studies. It encompasses conceptual geological and kinematic models, and a 2D finite element simulation to account for the processes undergone by the Zafarraya Fault. The study focuses on the variability of hy-dro-geomechanical features and the time evolution of the ground pore-pressure distribution in both the preseismic and coseismic stages, matching some of the shreds of evidence found by field studies. This methodology can be applied to other events registered in the National Catalogues of Earth-quakes to achieve a deeper insight, further knowledge, and a better understanding of past earth-quakes.

In January and February 2001, El Salvador was hit by two strong earthquakes that triggered thousands of landslides, causing 1,259 fatalities and extensive damage. The analysis of aerial and SPOT-4 satellite images taken a few days after the events allowed us to map 6,491 coesismic landslides, which occurred in 14 study areas extending for about 400 km2. Four different Multivariate Adaptive Regression Splines (MARS) models were produced by using different covariate sets and landslide inventories, the latter containing the slope failures triggered by an extreme rainfall event of November 2009 and those induced by the earthquakes of 2001. Moreover, two validation scenarios were employed, including the 25% and 95% of the mapped landslides, respectively. The results of our experiment revealed that: (i) MARS algorithm provides reliable predictions of coesismic landslides; (ii) models calibrated with rainfall-induced landslides predict with acceptable accuracy landslides caused by deep earthquakes and distributed over vast areas; (iii) the best accuracy is achieved by models trained with both preparatory and trigger variables; (iv) a small portion of the landslides produced by an earthquake can be used to calibrate MARS predictive models that help to identify slopes where yet unreported landslides may have occurred.

InSAR and associated analytic methods enable relative surface deformation measurements from low Earth orbit with a potential accuracy of centimeters to millimeters. However, assessing the actual accuracy can be quite difficult. The analytic methods are complicated enough that naive analytic error propagation is infeasible, and, in many settings, InSAR practitioners lack sufficient ground truth to assess results. Phase noise due to partial decorrelation from changes in the scattering properties of the ground is a prominent source of accuracy loss. In this paper we present a method to assess the loss of precision due to this component of phase noise. The proposed method consists of generating synthetic data stacks whose statistical properties match those of the actual input SAR data stacks, and then using the synthetic data for an ensemble calculation. The spread of the results of the ensemble calculation indicates the loss of precision. We show examples of the ensemble analysis at a mining operation in South Africa, and demonstrate the ability to estimate the precision of two InSAR deformation retrieval methods on a point-by-point and epoch-by-epoch basis.

The Lokpaukwu Uru Quarry was examined geologically, geophysically, and core-wise. The location is between 5056.149'N and 5056.193'N and 7028.312'E and 7028.356'E. The study location may include the Asu River Group and the Eze-Aku Formation. This area has five rock units. In the eastern research region, siltstone forms a "CAP" on the shale. Shale underlies half of the study area. The west has calcareous sandstone. The eastern part of the area is dolerite, the main rock that spans siltstone and shale. The region's geological matter contains iron. Two geological sections were analysed and interpreted to identify the five rock units and their outcrops in the study area. electroresistivity in geophysical research Schlumberger found that the western, northwesterly, and central sections of the research region had a thick sedimentary sequence, whereas the eastern half has an igneous body, the project's main component. Sandstone, siltstone, and shale follow the high-resistivity rock in this location. The rock unit in the region was found in eleven core samples from the east half of the study area. Nine rock-unit core samples were found near Obichioke. The Lokpaukwu area's core data shows the rocks' positions, kinds, minerals, and strengths. Geologic mapping shows that a major fault separates the viable Uru end from the unviable Obichioke lot. Recrystallization dominates the fault track (alcitic matter). Thus, prior to quarrying igneous (basic) units, comparable investigations are advised.

Landslide risk mitigation is limited by data scarcity. This could be improved using continuous landslide detection systems. In order to investigate which image types and machine learning (ML) models are most useful for landslide detection in a Norwegian setting, we compared the performance of five different ML models, for the Jølster case study (30-July-2019), in Western Norway. These included three globally pre-trained models; i) the Continuous Change Detection and Classification (CCDC) algorithm, ii) a combined k-means clustering and Random Forest classification model, and iii) a convolutional neural network (CNN), and two locally-trained models, including; iv) Classification and Regression Trees and v) a U-net CNN model. Images used included Sentinel-1, Sentinel-2, digital elevation model (DEM) and slope. The globally-trained models performed poorly in shadowed areas, and were all outperformed by the locally-trained models. A maximum Matthew’s correlation coefficient (MCC) score of 89% was achieved with model v, using combined Sentinel-1 and -2 images as input. This is one of the first attempts to apply deep-learning to detect landslides with both Sentinel-1 and -2 images. Using Sentinel-1 images only, the locally-trained deep-learning model significantly outperformed the conventional ML model. These findings contribute towards developing a national continuous monitoring system for landslides.

The access road to the powerhouse's hydraulic system's facilities in Guarumales (Azuay, Ecuador) presents a medium-sized landslide. Geophysical tests were conducted in the initial research stage, combining electrical and seismic methods. A SEV and HVSR surveys campaign have been taken as a reference for the analysis of the landslide area. The distribution of these test points has been at three different levels along the landslide where the access road crosses it and trying to cover the area's most extensive possible length and width. In the area, we find the geology dominated by the presence of schists, altered to different degrees and presenting blocks of material with a lower degree of alteration within colluvium formed by a clayey matrix and coarse material of the exact nature. There is also evidence of surface and groundwater in the area. The results obtained through SEV tests have allowed for defining the separation zone of the mobilized or sliding materials compared to the fixed or immobile ones (potentially, the sliding surface was marked). Using the HVSR technique, the natural vibration frequencies of the ground associated with the sliding mass (separation of seismic impedances between a two-layer model: mobile and fixed) have been determined. Previous authors proposed an empirical relationship establishing the exponential relationship, already proposed by previous authors, between sediment thickness and natural frequencies. It has been possible to determine the depth of the position of the loosely compacted sediment zone (and probably moving or mobilized) compared to that of compact materials (immobile) and thereby define the potential rupture surface.

This communication investigates possible anomalies in the lithosphere atmosphere and ionosphere on the occasion of the ML=3.3 earthquake that occurred on 1st January 2023 close to Guidonia Montecelio (Rome, Italy). This earthquake followed another event on 23 December 2022 of magnitude ML=3.1 with a very close epicentre (distance less than 1km). Seismological investigations clearly show an acceleration of seismicity in the last six months in the area. Two solutions of fitting time to failure power law on the Cumulative Benioff strain curve are the more likelihood: the ML3.3 of 1 January is the mainshock of seismic sequence or incoming earthquakes of a magnitude of about 4.1 provides a slightly better fit of the seismic data. Further investigation are necessary to assess if the accumulated stress have been totally released or not. No atmospheric anomalies related to this seismic activity have been identified even if some SO2 emissions could come from tectonic or volcanic sources in the South-Tyrrhenian Sea. Swarm satellites' magnetic data shows an anomalous track on 16 December 2022, which is temporally compatible with the seismic acceleration but other sources for the anomalous signal are also possible.

Drainage divide is a dynamic feature that migrates in response to tectonic activity. The asymmetric uplift between two adjacent basins causes the divide migration from a slower to faster uplift area. Sado Island, Japan, has been affected by southeastward tilting uplift since ca. 300k years. Despite the faster uplift on the northwest, the main divides have existed on the southeast side of the geometric center of the island, with no other feature suggesting tectonic inversion of the tilting direction. In this study, we conducted a DEM-based investigation that focused on divide migration. A spectrum from very inactive to active divide migration in the northwest. Regardless of their position, actively migrating divides are comprehensible, but inactive divides located in a relatively slow uplift area remain unclear. We concluded that some divides slowed down owing to the local balance of erosion rates across the divides, not implying the balance between uplift and river erosion at the basin scale, reflecting disequilibrium in river longitudinal profiles. The main divides of Sado have presumably continued to slowly migrate toward the faster uplift area; however, they are most likely to have never overcome the moving geometric center owing to land expansion at the seacoast due to asymmetric uplift.

As an effective active remote sensing (ARS) technology for shallow underground targets, ground penetrating radar (GPR) is a detection method to obtain the characteristic information of underground targets by transmitting electromagnetic wave from the antenna and analyzing the propagation law of the electromagnetic wave underground. Due to the high frequency(1MHz-3GHz) of GPR, the depth of geological exploration is shallow(0.1m-30m). In order to remote sensing the deep earth, it is necessary to increase the size of the radiation source in order to reduce the radiation frequency. At the same time, for most separated GPRs, a single dipole antenna (SD) is used as the radiation source and another antenna device placed along the electromagnetic wave propagation direction in the far field as a remote sensing sensor (RSS), both of which are horizontally linearly polarized (LP) antennas. In some cases, such a design is apt to cause problems such as multipath effect (ME) and polarization mismatch (PM). When GPR in ARS of deep earth is performed, it often results in increased errors, signal attenuation during data reception and processing. In contrast, at the radiation source, with the use of large aperture multiple-dipole antennas (MD) and multi-channel sequential rotational excitation, the electromagnetic wave can radiate outward in the form of circular polarization at a low frequency. At the RSS, the trouble caused by ME and PM can be reduced even if the LP antennas are used. A novel sequential rotationally excited (SRE) circularly polarized (CP) MD array for separated GPR in ARS of deep earth is proposed in this paper, which uses a large aperture CP MD array instead of a small size LP SD. The analysis and simulation results demonstrate that under the premise of the same transmitting power, comparing circular polarization and linear polarization, by using SRE CP MD antennas array radiation source, a significant enhancement (about 7dB) of the Signal to Noise Ratio (SNR) will occur by collecting the radiant energy at the RSS. More importantly, by reducing the exploration frequency to 10KHz, the exploration depth will also be greatly increased by about 10 times.

The increased accessibility of drone technology and the wide use of Structure from Motion 3D scene reconstruction have transformed the approach for mapping inaccessible slopes undergoing active rockfalls. The Poggio Baldi landslide offers the possibility for many of these techniques to be deployed and integrated with the aim of defining a suitable workflow for the analysis of hazards in mountainous regions. The generation of multitemporal digital slope twins (2016 – 2019), informed a rockfall trajectory analysis that was carried out with a physical-based GIS model. We tested the rockfall scenario reconstructed and calibrated on the analysis of the rock mass characteristics and the geometrical and physical constraints given by the multi-temporal analysis of the SfM point clouds. This time-independent rockfall hazard analysis is a critical component to any subsequent holistic risk analysis on this case study, and any potential similar mountainous setting.

A novel phase-field method based numerical approach to modeling the compaction process of sedi- ments is presented. Water-logged rock or soil sediments are deposited on the water basins over time and the increasing volume of the sediment compacts under its own weight and external pressure. Coupled evolution of mass conservation, Darcy flow, and the viscoelastic constitutive response, in conjunction with the evolution of porosity and permeability, make this problem highly non-linear and involve moving boundaries. We adopt a phase-field approach to represent the moving sediment interface, and the underlying multiphasic model permits for studying compaction in dynamically evolving sediments. This approach does not necessitate a change of domain sizes as is the case of existing traditional models of sedimentation but rather treats sedimentation growth as a problem of interface motion. We first model a classical compaction problem in 1D to compare with existing results in the literature, and then extend the framework to 2D to model the compaction process taking place under the influence of gravity. The model is then extensively applied to understand the effect of the sediment initial state and the sediment material properties on the compaction process and its spatiotemporal evolution. Lastly, a strong validation of our phase-field treatment results against predictions from traditional methods of solving open boundary sediment compaction problems, and a good agreement between the conventional understanding and the numerical predictions of porosity dependence on the fluid pore pressure for various cases of geological interest, are used to demonstrate the applicability and scalability of this novel numerical framework to model more advanced sediment compaction problems and geometries.

The Marchenko method is a data-driven way which makes it possible to calculate Green's functions from virtual points in the subsurface by the reflection data at the surface, only requiring a macro velocity model. This method requires collocated sources and receivers. However, in practice, subsampling of sources or receivers will cause gaps and distortions in the obtained focusing functions and Green's functions. To solve this problem, this paper proposes to integrate sparse inversion into the iterative Marchenko scheme. Specifically, we add sparsity constraints to the Marchenko equations and apply sparse inversion during the iterative process. Our work not only reduces the strict requirements on acquisition geometries, but also avoids the complexity and instability of direct inversion for Marchenko equations. This new method is applied to a two-dimensional numerical example with irregular sampled data. The result shows that it can effectively fill gaps of the obtained focusing functions and Green's functions in the Marchenko method.

Recently on the November 4, 2022 at 03:00 PM, a M_ww 5.6 earthquake struck Cianjur district of West Java province, Indonesia. To our best knowledge, an earthquake will cause a massive damage to the manmade infrastructure and this responsible for a massive casualty. A building damage is affected by various determinant factors. One of important factors is the geological bedrock variations. Here this research aiming to elaborate the associations of geological bedrocks in Cianjur district with residential building damages caused by recent earthquakes. The results confirm that majority of building damages account for 54.67% of the total building was associated with quaternary volcanic bedrocks. The second massive building damages were observed for Andesitic Dacitic volcanic rocks accounted for 32.07%. To conclude, elaborations of geological bedrock-building damage will contribute to the better planning in the future that can reduce the damages in the urban settings.

Las Loras UNESCO Global Geopark (UGGp) is geologically diverse, particularly in relation to water-derived features: springs, karst springs, travertine buildings, waterfalls, caves. In this work, the interactions between geology, geomorphology, structures and hydrogeology are analyzed. These four components are the fundamentals of geodiversity and their interactions provide a conceptual model of hydrogeological functioning at Las Loras UGGp. The most plausible hypothesis is that the system is formed by two superimposed aquifer systems, separated by an aquitard formed by Lower Cretaceous material. The deep lower aquifer formed by the Jurassic limestones only outcrops on the northern and southern edges of the Geopark and in a small arched band to the south of Aguilar de Campoo. It forms a basement subject to intense deformation. The upper aquifer system, formed by outcropping materials from the Upper Cretaceous, is a free aquifer. It is formed by a multilayered aquifer system that is highly compartmentalized by the very disturbed geomorphology of the landscape, constituting each moorland and each lora, an individualized recharge-discharge system. This model explains the base level of the rivers, the abundant number of existing springs and the permanent nature of some rivers, providing the keys to understand the bases of the geoconservation of a rich geological heritage linked to the active processes of the water cycle.

Abstract: The empirical soil loss model, RUSLE, was used with a conjunction of remotely sensed data and geographic information system technology to delineate the soil erosion and watershed priorities in terms of conservation practices at seven boundary sub-watersheds (labeled as SW-00, SW-01, …, SW-06) between Iraq and Iran at Ali Al-Gharbi area, southern Iraq. The six factors of the RUSLE model, ie. the rainfall erosivity, the soil erodibility, the slope steepness length, the crop management, and management practice, were calculated or estimated using information from different data sources such as remotely sensed data and previous studies. The finding showed that the annual soil erosion loss ranges from 0 - 1890 (tons h-1 y-1) with an average of 0.66 (tons h-1 y-1). Values of soil erosion were classified into five classes: very low, low, moderate, high, and very high. The potential soil loss in the high and very high classes ranges from 14.84 to 1890 (tons h-1 y-1), and these classes occupy only 27 km2 of the study area, indicating that the soil loss is very low in the area being examined. In terms of the spatial distribution of soil loss, the northern and northeastern parts (mountains and hills) of the sub-watersheds where the slope is steeper are more likely to erode than the plain area in the southern and southeastern portions, indicating that slope, in addition to rainfall erosivity, has a dominant effect on the soil erosion rate. The study of soil erosion in the watersheds under consideration reveals that only the northern portions of the SW-00, SW-02, and SW-04 watersheds require high priority conservation plans; however, these portions are primarily located in mountain regions, making conservation plans implementation in these areas impractical. Due to low soil loss, other sub-watersheds, particularly SW-01, SW-03, SW-05, and SW-06, are given low priority.

Identifying ambient noise-based (ANb) signatures of streams can help in the estimation of their erosive potential (EP) that promotes reverie landslides and soil losses in the fluvial valleys. This is particularly imperative on flooding or rainy days, leading to stronger erosion-prone conditions (colluvium and boulders) of the valley beds inferred from georadar attribute analysis. Developing such research direction can benefit the local communities, as is the case with the Cerrado region of Brazil, where these phenomena have high destructive potential with social, economic, and climatic implications. For the present study, a seasonal stream in the Federal District of Brazil was investigated by ANb monitoring supported by Ground Penetration Radar (GPR) for site characterization. The ANb monitoring was conducted (at a safe distance) with a seismometer over several durations of dry and rainy conditions. The power spectral density (PSDs) was computed as a function of several variables, including weather conditions (rainfall, wind speed, and pressure), time-frequency spectrograms, and ambient noise displacement root mean square (RMS). This analysis also considered the single station horizontal-to-vertical spectral ratio (HVSR), where rain, wind, pressure, river flow and anthropogenic signatures were evident (at selective frequency ranges). Multi-peaks that emerged on the HVSR curve were further analyzed to identify amplitude and frequency changes, and the three peaks shift on average to a lower position during the rainy period. The GPR amplitude and waveform variation features were attributed to the stratigraphy (i.e., the boundary between valid and invalid regions and coherence value) of the floodplain and regions susceptible to erosion (erosion-prone lithological spots). This approach provides the basis for non-destructive monitoring tools enabling the detection of 'seismic signatures' and weak spots of the fluvial channels for improving environmental management.

This study presents the developments regarding the time-domain-induced polarisation method, as a supporting tool for resistivity soundings during investigations of coastal detrital aquifers that are salinised by marine intrusion. The interpretation of resistivity measurements in such aquifers, which have variable hydrochemistry and lithology, involves uncertainties owing to the presence of low-resistivity lithologies, such as clays. To reduce these uncertainties, the use of other geophysical parameters is necessary; hence, this study focuses on induced polarisation, since it can be measured simultaneously with resistivity. In detail, we propose the determination of induced polarisation using 1D techniques, while developing a different algorithm for processing the induced polarisation data. The aim is to extend the results of this phenomenon, using instead of chargeability, the concepts of polarisability and decay time, which are extracted from the decay curve, given that they represent more intrinsic properties of the various analysed subsurface media. We present results obtained by applying this methodology to a Quaternary aquifer of the Costa del Sol in the SE Iberian Peninsula (in the province of Almería) during two different campaigns, one before and one after winter (i.e., in October and February, respectively). The results reveal the position of the saline front during each campaign, while reflecting the seasonal movement of the marine intrusion.

The Danxia landform area of Jingbian Wave Valley is located in the central part of Ordos Basin. The near-red Danxia landform consisting of sandstone in this area is a new discovery in geological circles at home and abroad, and its depositional age and genesis remain a hot topic of discussion at present. As the material basis for its development, red sandstone is of great importance to the in-depth study on its formation date and origin. The paper explores the origin, tectonic significance and paleogeographic pattern of the red sandstone through field geological investigation and zircon U-Pb dating analysis of the red sandstone of the Luohe Formation. The results show that the original materials of the red sandstone was formed in three main age intervals, 252~456 Ma, 1657~2084 Ma, and 2129~2538 Ma, and regional tectonothermal events in the origin area during the three periods were recorded, among which the events in Mesoproterozoic and Paleoproterozoic periods were the most active. The comparative analysis of zircon U-Pb age spectra shows that for the sediment of red sandstone, the early origin is near-origin and late origin is distant-origin. The main origin areas in the study area are the Xingmeng orogenic belt, the North China massif and the Alashan massif; the vertical section reflects that the origin area has the strongest erosion and transport in the upper deposition period, the smallest erosion and transport in the middle deposition period, and the moderate erosion and transport in the lower deposition period. With origins system dominated by the northern origins and the paleogeographic pattern of multi-period and cyclonic, the area mainly experienced four tectonic movement cycles, such as Wutai, Lvliang, Caledon and Haixi.

The coastal Peru, one of the driest deserts in the world, has recently experienced dramatic landscape changes and asset destruction during precipitation events due to El Niño. Nevertheless, catastrophic explanations for landscape variations and human responses to the so-called "Super" or "Mega" El Niños are recurrent in pre-Hispanic time, even if alternative hypotheses were provided in literature. A deeper understanding of the geological and archaeological record can improve the knowledge on the relationship between such a coupled atmosphere-ocean phenomenon and landscape processes. The bibliographic sources required for this purpose are scattered throughout various disciplines ranging from physical to human sciences, thus comprehensive databases were used to identify and screen the relevant studies. The results provide knowledge synthesis in order to identify critical gaps and suggest specific research goals. Inferred episodes of landscape change due to severe precipitations in late pre-Hispanic time are discussed and consistencies and inconsistencies exposed. Examples of variation in landscape response due to extensive human intervention are reported.

Predicting the tidal wave is essential not only to better understand hydrological cycle at the boundary between land and ocean but also to improve energy production in the coastal area. As affected by various factors such as astronomical, meteorological, and hydrological effects, predicting the tidal wave at the estuary remains uncertain. In this study, we present a novel method to improve short-term tidal wave predictions using a fixed-lag smoother based on sequential data assimilation (DA). The proposed method is implemented for the tidal wave predictions at the estuary of Nakdong River. As a result, the prediction accuracy was improved by 63.9% with DA and the calibration using the regression. Although the accuracy of DA diminished for increasing forecast lead times, the 1-hr lead forecast by DA had still 44.4% improvement over the open loop without DA. Plus, the optimal conditions for the fixed-lag smoother were analyzed in terms of the order of a smoothing function and the length of assimilation window and forecast lead time. It was suggested that the optimal DA configuration could be obtained with the 8th order polynomial as a smoothing function assimilated under 6-hr or longer past and future DA windows.

One of the possible mechanisms causing significant emissions of methane into the atmosphere within the Arctic shelf may be the decomposition of gas hydrates. Their accumulations within the Arctic shelf formed during Ice Age almost simultaneously with the formation of permafrost, which contributed to the emergence of a zone of stable existence of gas hydrates. The subsequent flooding of the Arctic shelf led to the degradation of the permafrost and the violation of the conditions for the existence of hydrates. To assess the state of the stability zone, methods of mathematical numerical modeling are used. Standard seismic methods are widely used to localize gas hydrates, but monitoring their physical state requires the development of fundamentally new approaches based on solving multiparameter inverse seismic problems. In particular, the degree of attenuation of seismic energy is one of the objective parameters for assessing the consolidation of gas hydrates: the closer they are to the beginning of decomposition, the higher the attenuation, and hence the lower the quality factor. Thus, the methods of seismic monitoring of the state of gas hydrates in order to predict the possibility of developing dangerous scenarios should be based on solving a multi-parameter inverse seismic problem. This publication is devoted to the presentation of this approach.

In order to understand the impact of lithologic heterogeneity of continental mixed fine-grained sedimentary rocks on reservoir brittleness in western Qaidam basin, the mechanical properties of the rocks and their correlation with mineral composition and petrographic characteristics were studied by measns of mineralogy, petrography and triaxial stress test. The results show that the reservoir rocks can be divided into 5 different types according to the mechanical properties of the reservoir (characteristics of stress-strain curves), among them Type I and III belong to similar elastoplastic failure model, type II shows a special pulse failure mode for plastic material, type IV shows a failure mode of mixed characteristics, and type V exhibits a typical plastic failure model. The correlation between minerals and mechanical properties indicates that quartz and feldspar, which are often considered brittle minerals, do not contribute much to the brittleness of continental fine-grained sedimentary rocks. The main minerals affecting reservoir brittleness are dolomite and clay minerals, and their contributions to reservoir brittleness are positive and negative, respectively. The petrographic analysis results prove that the abnormal correlation between rock mechanical properties and quartz and feldspar is caused by the different rock fabrics. When dolomite forms a rock skeleton, it typically exhibits greater strength, brittleness and physical properties than other minerals. Based on the results, a brittleness evaluation standard for continental fine-grained sedimentary rock reservoir is proposed, and the validity of the standard is verified by the spatial correlation between the lithology probability model and the micro-seismic monitoring data, indicating that the spatial heterogeneity of dolomite-rich rock is the main controlling factor for the development of ‘sweet spot’ in the Cenozoic continental fine-grained sedimentary rock reservoir in the Western Qaidam Basin.

Sediment thickness increases can cause floodplains and the water level increases. This has the potential to generate a flood. Using electromagnetic waves, Time Domain Electromagnetic (TDEM) detected resistivity or conductivity contrast of lithology in the subsurface. It is measured in the time domain. TDEM method has been developing for decades. Here we tried to develop a 1-D forward modelling program for central loop configuration in the water environment using the Adaptive Born Forward Mapping (ABFM) method. We simulated this program in several water environment conditions (such as freshwater, brackish water and saline water) to know its response. Preventing natural hazards, especially flood hazards which are caused by the floodplain increases is our motivation in this research. Our simulation shows that Central-Loop Configuration Time-Domain Electromagnetic Method is able for imaging the sediment thickness clearly. The response of this method is extremely sensitive in saline water to depth changing than in other water environments.

The PB Sudirman hydroelectric power plant, built in 1988, has fewer than 50 years of operational life. In contrast, the planned operating life for this dam is 50 years. This dam is experiencing intensive sedimentation, which results in siltation. Dredging efforts and catchment area management are conducted to extend the operational life of the dam. Therefore, technical and economic analysis of dredging activities and the potential benefits of dredged sediment is needed to determine the solution's effectiveness. We use the site survey method to conduct an environmental impact analysis and interviews with dam managers. The identification of the extent and distribution of erosion-sedimentation was carried out as part of the technical assessment. In the economic study, we calculate the dredging costs and profits from selling dredged sediments, mainly sand. The analysis results show that the dredging of deposits in the foreset slope area is not profitable. There are 17 areas of bank spoil that have economic benefits. This study concludes that the need for dredging costs on foreset slopes is higher than the profit from selling dredging sediments.

Soil texture is key information in agriculture for improving soil knowledge and crop performance, so the accurate mapping of this crucial feature is imperative for rationally planning cultivations and for targeting interventions. We studied the relationship between radioelements and soil texture in the Mezzano Lowland (Italy), a 189 km² agricultural plain investigated through a dedicated airborne gamma-ray spectroscopy survey. The K and Th abundances were used to retrieve the clay and sand content by means of a multi-approach method. Linear (simple and multiple) and non-linear (machine learning algorithms with deep neural networks) predictive models were trained and tested adopting a 1:50,000 scale soil texture map. The comparison of these approaches highlighted that the non-linear model introduces significant improvements in the prediction of soil texture fractions. The predicted maps of the clay and of the sand content were compared with the regional soil maps. Although the macro-structures were equally present, the airborne gamma-ray data permits us shedding light on finer features. Map areas with higher clay content were coincident with paleo-channels crossing the Mezzano Lowland in Etruscan and Roman periods, confirmed by the hydrographic setting of historical maps and by the geo-morphological features of the study area.

In Qinglong ore field, paleo-oil reservoir is found to be associated with antimony deposits, and they have close genetic relationship. In this study, aromatics geochemistry of paleo-oil reservoir bitumen was studied to further discuss the thermochemical sulfate reduction (TSR) reaction and the mechanism of antimony mineralization. A total of 124 aromatic compounds were identified by Gas chromatography-mass spectrometry (GC–MS) analysis in bitumen samples, including abundant phenanthrene series, dibenzothiophene series, fluoranthene series, chrysene series, and a small number of fluorene series, naphthalene series, dibenzofuran series, biphenyl series, triaromatic steroid series. Aromatic parameters such as trimethylnaphthalene index (TMNr), methylphenanthrene index (MPI), methylphenanthrene distribution fraction (MPDF, F1 and F2), methyldibenzothiophene parameter (MDR), C28TAS-20S/(20R+20S), and benzofluoranthene/benzo[e]pyrene indicate Qinglong paleo-oil reservoir is in over maturity level. The abundance of phenanthrene and chrysene aromatic compounds and a small amount of naphthalene series, benzofluoranthene, fluoranthene, pyrene, anthracene, retene, perylene and biphenyl suggest that the parent material of the paleo-oil reservoir was mainly low aquatic organisms, mixed with a small amount of higher plant. The detected a certain number of compounds, such as retene, triaromatic steroid series and perylene, the ternary diagram of DBF–DBT–F and binary plot of Pr/Ph–DBT/P, DBT/(F+DBT)–DBF/(F+DBF) and Pr/Ph–DBT/DBF reveal that the source rock of the paleo-oil reservoir was formed in the marine environment of weak oxidation and weak reduction. The comprehensive analysis shows that the Qinglong paleo-oil reservoir originated from Devonian source rocks, just like other paleo-oil reservoirs and natural gas reservoirs in the Nanpanjiang basin. Abundant dibenzothiophene series were detected, indicating that the paleo-oil reservoir underwent a certain degree of TSR reaction. We believe that the paleo-gas reservoir formed by the evolution of paleo-oil reservoir participated in antimony mineralization, that is, hydrocarbon organic matter acted as reducing agent and transformed SO42− in oilfield brine into H2S through TSR, providing reduced sulfur and creating environmental conditions for mineralization.

The paper concerns the distribution of apparent magnetic susceptibility in soil profiles located in the areas of topsoil magnetic susceptibility anomalies in Krakow. The type of land use, possible sources of magnetic carriers, and the type of soil were taken into account. Additionally, at each soil profile, a comparison between soil magnetic susceptibility and the results of geochemical analyzes of soil samples was made. The study shows very characteristic changes in magnetic susceptibility with depth, reflecting the interdependencies between natural and anthropogenic factors. A visible magnetic susceptibility maximum at the depth of 10-30 cm is observed at each soil profile. The maximum is associated mainly with the deposition of atmospheric dust and its vertical range depends on the level of anthropopression and natural conditions of soils. At the depth above 40 cm in the eastern part of Krakow, a correlation between the magnetic susceptibility and the soil type (chernozems de-veloped on loess) was found. All indicates that the thicknesses of contaminated upper horizons are not accidental and they depend on human interactions with the environment and the type of soil. An attempt at template establishment with the sources of magnetic particle carriers for different places in the city was made. As the result, in high urbanized sites, the extreme values of magnet-ic susceptibility rapidly change in short vertical distances can identify the richness of anthropo-genic layers with various types of anthropogenic ferrous material and/or additionally Fe-carrying objects buried in soils. In industrialized sites, anthropogenic input plays the most important role in the creation of soil magnetic characteristics. What is more, industrial pollution hides the natural magnetic properties of chernozems. In opposite, the studies at the sites under low anthropopression (mainly located in forests) allow for better insight into magnetic proper-ties arising during pedogenic processes, indirectly giving information about soil conditions. In the forest areas, the lowest values of soil magnetic susceptibility were measured. Additionally, the influence of pedogenic and lithogenic factors on forest soils is manifested in the results. Among the sites concerned, particular attention should be paid to the vicinity of the steel plant because of the agricultural land in the surroundings.

Due to the use of structural and seismic techniques an attempt has been made in this study to determine the hydrocarbon potential and future exploration in Pakistan's Lower Indus basin area and offshore region. In addition to the possibility of hydrocarbon accumulation, the major goals of this research are to analyze structural patterns in the subsurface of the study region, identify the horizons of various formations using surface seismic, and interpret seismic lines. Seismic data has been used to find and integrate the exposed Cretaceous/Tertiary basalt at Ranikot in the lower Indus basin to the well-known Deccan basalts of the Indo-Pak plate, which is the secondary goal of the current investigation. According to the interpretation of seismic data, the Indo-Pakistan plate has experienced rifting because of tectonic activity. Graben and horst structures have been identified in the study area horizon, which indicates that the area has undergone major structural and depositional changes. While faults provide a pathway for hydrocarbons to migrate from their source to reservoir rock, grabens are principally responsible for the accumulation of hydrocarbons. We further deduce that the prior exploration failures in the area were due to a lack of knowledge of the subsurface formations and structural trends, which are essential for characterizing the hydrocarbon play and trap features. Our results can be productive for the local hydrocarbon drilling projects as well as worldwide tectonic stratigraphy studies on passive continental edges.

INT-FFT algorithm presented in this work uses an interpolation methodology to densify 3D-GPR datasets to sharpen images obtained in GPR surveys obtained in an archaeological context. It allows the reconstruction of missing data from the combined use of mathematical transforms (e.g., the Fourier and Curvelet transform) and predictive filters. This technique makes it possible to calculate the missing signal simply by meeting two requirements: the data in the frequency domain must be limited in a range of values and must be able to be represented by a distribution of Fourier coefficients (verified conditions). The INT-FFT algorithm uses an open-access routine (Suinterp, Seismic Unix) to interpolate the GPR profiles based on seismic trace interpolation. This process uses automatic event identification routines by calculating spatial derivatives to identify discontinuities in space by detecting very subtle changes in the signal, thus allowing for more efficient interpolation without artifacts or signal deterioration. We successfully tested the approach using GPR datasets from the Roman Villa of Horta da Torre (Fronteira, Portugal). The results show an increase in the geometric sharpness of the GPR reflectors and have not produced any numerical artifacts. The tests performed to apply the methodology to GPR-3D data allowed for assessing the interpolation efficiency, the level of recovery of missing data, and the level of information lost when one chooses to increase the distance between profiles in the acquisition stage of the data.

InSAR and associated analytic methods can measure surface deformation from low earth orbit with a claimed accuracy of centimeters to millimeters. The realized accuracy depends on the area being measured and on the choice of analytic method, suggesting one choose a method in response to the area being measured. Here we consider a specific fixed analytic method and compare the results it produces to measurements gathered from other means in a variety of settings. In particular we compare Sentinel-1 InSAR with GPS at the Kilauea volcano around the 2018 eruption, with GPS in the city of Arica, Chile, and with public survey data at a decommissioned tailings mine. In addition, we compare two independent Sentinel-1 InSAR analyses for a railway station in Oslo, Norway. Our goal is estimate the accuracy of a fully automated Sentinel-1 InSAR pipeline in various settings. Our conclusions are that centimeter level accuracy is a reasonable claim in many, but not all settings, and that accuracy is typically not lost by using an automated pipeline, instead of hand-selecting and tuning parameters.

In this review, the potential of an emerging field of interdisciplinary climate research, that is the Coastal Boulder Deposits (CBDs) as natural archives for intense storms, is explored with particular reference to the Mediterranean region. First, the identification of the pertinent scientific articles was performed by using Web of Science (WoS) engine. Thus, the selected studies have been analysed to feature CBDs produced and/or activated during the last half century. Then, the meteorological events responsible to the literature reported cases were analysed in some details using the web archives of the Globo-Bolam-Moloch model cascade. The study of synoptical and local characteristics of the storms involved in the documented cases of boulder production/activation proved useful to assess the suitability of selected sites as geomorphological storm proxies. It is argued that a close and fruitful collaboration involving several scientific disciplines is required to develop this climate research field.

A single fracture is the basic unit of fracture medium, and the roughness of fracture wall surface is an important factor influencing hydraulic characteristics of the flow in bedrock fracture. However, effects of the shape and density of roughness elements (various bulges/pits on rough fracture wall surfaces) on water flow in a single rough fracture have not been thoroughly discovered. Thus the water flow in single fracture with different shapes and densities of roughness elements was simulated by using Fluent software in this study. The results show that in wider fractures the flow rate mainly depends on fracture aperture, while in narrow and close fracture medium the surface roughness of fracture wall is the main factor of head loss of seepage; there is a negative power exponential relation between the hydraulic gradient index and the average fracture aperture, i.e. with the increase of fracture aperture, the relative roughness of fracture and the influence weight of hydraulic gradient both decrease; and in symmetrical-uncoupled fractures there is a super-cubic relation between the discharge per unit width and average aperture. Above results would help to deepen the understanding of rough fracture seepage.

The easternmost Mediterranean is a distinct transition zone from the ocean to the continent located at the junction between the largest Earth's lithospheric segments: Eurasian and African. The methodology of paleomagnetic mapping of such transition zones is based on integrating the mapping techniques for both continental and oceanic platforms: paleomagnetic reconstructions, results of radiometric dating of magnetized rocks, tectonic-structural reconstructions, biogeography, and utilization of the results of various geophysical surveys. The geodynamic-paleomagnetic mapping makes it possible to reveal the multilevel structural heterogeneity and display complex elements of the geodynamics of different ages inherent in this transition zone. Northern Israel is obviously the most complex area in the easternmost Mediterranean. For the combined paleomagnetic mapping, well-studied paleomagnetically and radiometrically areas were selected: (1) the Carmel area, (2) the Atlit area (internal part of the Carmel area), (3) the Sea of Galilee with the adjoining zones (primarily, the Kinnarot Valley), and (4) the area of the Hula Basin with adjacent areas of the Golan Plateau, Hermon Mt., and Galilea uplift. The constructed paleomagnetic profiles for the Carmel area (on the top of the accumulative surface of the Lower Cretaceous traps), and the Kinnarot Valley – Sea of Galilee – Hula Basin, evidently indicate the complex history of the paleogeodynamic evolution of the region. These studies demonstrate the effectiveness of paleomagnetic mapping integrated with paleomagnetic profiling, crossing these geologically complex areas.

Seismic monitoring in areas where induced earthquakes could occur is a challenging topic for seismologists due to generally very low signal to noise ratio. Therefore, the seismological com-munity is devoting several efforts to the development of high-quality networks around the areas where fluid injection and storage and geothermal activities take place, also following the national induced seismicity monitoring guidelines. The use of advanced data-mining strategy, such as template matching filters, auto-similarity search and deep-learning approaches is recently further fostering such a monitoring enhancing the seismic catalogues and lowering the magnitude of completeness of these areas. In this framework, we carried out an experiment where a small-aperture seismic array was installed around the gas reservoir of Collalto, in North Italy. The continuous velocimetric data, acquired for 25 days, were analysed through the application of the optimized auto-similarity search technique FAST. The array was conceived as a cost-effective network, aimed at integrating, right above the gas storage site, the permanent high-resolution Collalto Seismic Network. The analysis allowed to detect micro-events down to magnitude Ml=-0.4 within a distance of ~15km from the array. Our results confirmed that the system based on the array installation and the FAST data-analysis might contribute to lower the magnitude of completeness around the site of about 0.7.

The 26 June 1917 tsunamigenic earthquake in Samoa is considered the largest historical event on record to have impacted this region in terms of earthquake magnitude and intensity. Yet, very little is known about the scale and distribution of tsunami impacts compared with the recent 2009 event which originated about 150 km east along the subduction zone bend of the Northern Tonga Trench (NTT). In this study we set out to: 1) reconstruct the 1917 tsunami from source to inundation to understand its hazard risk characteristics in the Samoan islands of Savai’i and Upolu; and 2) assess the hazard implications of tsunamis sourced from different locations along the subduction zone bend of the NTT on present-day exposure of coastal assets relative to the 2009 tsunami benchmark. We use the BG-Flood numerical modelling suite to produce model outputs representing inundation extent and hazard depth intensities at spatially flexible grid resolution (10 m, 20 m and 40 m). These are validated using available tide gauge records in Apia harbour and limited observations of runup that were derived from historical records. We then combine the inundation model with available digital distributions of buildings in the RiskScape multi-hazard risk analysis software, to produce exposure metrics for understanding the likely impacts on present-day coastal asset and population distributions if a similar tsunami were to occur. Results of the tsunami modelling indicate variable modelled-to-observed consistency using available source models, wave and runup validation data. Discrepancies in recorded vs modelled wave arrival time at Apia of between 30—40 mins are observed, with modelled runup underestimated in southeast Upolu Island compared with the rest of the country where runup observations are available (e.g., Savai’i Island). These differences likely reflect complexities in the tsunami source mechanism which might not currently be represented in our modelling. Nevertheless, our results suggest that a larger proportion of people would be exposed in Savai’i island (71% of exposure total), compared with Upolu island if a characteristic 1917-type event were to occur. While this study provides the first detailed inundation model of the 1917 tsunami in the Samoan region, the observed discrepancies suggest that further investigation is required to constrain potential tsunami source complexities which might not be accounted for in this study. Notwithstanding these limitations, our findings help to reinforce an appreciation of the risk to the greater Samoan region faced by local tsunamis sourced at different locations along the subduction zone bend of the NTT.

Archie model is the basis of calculating oil saturation, but there are some limitations when using this model to calculate oil saturation in water-flooded layer. The main reason is that the main parameters, such as rock resistivity and formation water resistivity, are constantly changing dynamic parameters in the underground with the influence of injected water. Considering that water-flooded layers changes with injection waters, influence factors of rock resistivity and formation water resistivity of primary parameters are analyzed. Considering the dynamic data of water cut is the most reactive underground fluid characteristics of real information, combined with dynamic and static, and the “double ratio” model of later development by the same sedimentary layer is established, which realizes the inversion of rock resistivity and formation water resistivity, then modified Archie model relating to water-flooded layers of the same sedimentary explaining formation. The explanation of actual data indicates that the “double ratio” model well considers the dynamic variation of production data, which makes the inversion of the flooded rock resistivity relatively accurate, besides, the modified Archie model can accurately calculate the oil saturation of water-flooded layers with a reasonable result, which offers scientific basis for the predicting of remaining oil distribution rules.

Undisturbed ground temperature (UGT), thermal conductivity (TC) and heat capacity (HC) are essential parameters for the design of borehole heat exchanger (BHE) and borehole thermal energy storage systems. However, field methods to assess the thermal state and properties of the sub-surface are costly and time consuming. Moreover, HC is often not evaluated but arbitrarily selected from literature considering the geological materials intercepted by boreholes. Therefore, this work aims at proposing a field heat tracing method to infer the thermal diffusivity (TD) and HC with assumption of natural transient heat conduction in the subsurface. Empirical equations were developed to reproduce a UGT profile measured along a BHE. Experimental coefficients are found with a non-linear least square solver optimization and used to calculate the damping depth and TD. Subsequently, the TD is used to evaluate HC considering TC obtained from a thermal response test (TRT). Results from this proposed heat tracing method were verified and validated against a set of TRT results and oscillatory TRT analysis using a field dual probe concept to infer HC. The example here described highlights the advantages and novelty of this fast and simple field method relying only on a single UGT profile measured before a TRT.

We examined PT parameters, geochemical peculiarities, and fluid sources of the Ulug-Sair ore occurrence attributable to class of intrusion-related gold deposits and according to ore mineral assemblages corresponding to Au-Bi type with wide Bi minerals (AgBiTe, Bi2Te2Se, Cu3,07BiS3, Bi), tellurides (Au and Ag), Se-tellurides (Ag and Bi), and selenides (Au, Ag, and Hg). We identified that ‘pre-gold’ quartz-tourmaline veins were deposited using an aqueous Mg-Na-K-chloride fluid with a salinity of 8–10 wt % NaCl eq. At 325–370 °C; host breasts were formed due to a CO2-water fluid containing CH4 and N2, with a salinity of 0,18–6,1 wt % NaCl eq. at least 200–400 °C. Gold-bearing mineral assemblages were formed at P ~ 0,75–1,0 kbar (~ 2,3–3 km) due to CO2-water chloride (Na-K±Fe, Mg) fluid with CH4, Na2SO4, and Na2B2O5, and salinities 1,7–12,5 wt % NaCl eq. during the decreasing temperatures from 360 up to 115 °C (gold-sulfide-quartz veins – 360–130 °С, and gold-telluride-sulfide-quartz veins – 330–115 °C) and variations fO2, fS2, fSe2, and fTe2. The isotopic composition of δ34SH2S fluid (-0,7…+2,5 ‰) indicates the juvenile or magmatic origin of fluid and ore elements. The δ34OH2O fluid indicates that, at an early substage, the formation of ore occurrence involved a fluid of magmatic or metamorphic origin (+7,3…+11,4 ‰), and, in the later substage, it mixed with meteoric waters (-2,3…+9,1 ‰).

Surface velocity is traditionally measured with in situ techniques such as velocity probes (in shallow rivers) or Acoustic Doppler Current Profilers (in deeper water). In the last years, researchers have developed remote sensing techniques, both optical (e.g., image-based velocimetry techniques) and microwave (e.g., Doppler radar). These techniques can be deployed from Unmanned Aerial Systems (UAS), which ensure fast and low-cost surveys also in remotely-accessible locations. We compare the results obtained with a UAS-borne Doppler radar and UAS-borne Particle Image Velocimetry (PIV) in different rivers, which presented different hydraulic–morphological conditions (width, slope, surface roughness and sediment material). The Doppler radar was a commercial 24 GHz instrument, developed for static deployment, adapted for UAS integration. PIV was applied with natural seeding (e.g., foam, debris) when possible, or with artificial seeding (woodchips) in the stream where the density of natural particles was insufficient. PIV reconstructed the velocity profile with high accuracy typically in the order of a few cm s⁻¹ and a coefficient of determination (R²) typically larger than 0.7 (in half of the cases larger than 0.85), when compared with acoustic Doppler current profiler (ADCP) or velocity probe, in all investigated rivers. However, UAS-borne Doppler radar measurements show low reliability because of UAS vibrations, large instrument sampling footprint, large required sampling time and difficult-to-interpret quality indicators suggesting that additional research is needed to measure surface velocity from UAS-borne Doppler radar.

The forming of large rivers are the integral consequences of the deep earth process and the surface. In contrast to the hot topics for rivers related to orogenic domains, rift-related large rivers are largely ignored especially in deep time studies. The Cenozoic East Asia margin provides very good opportunity to observe this kind of rivers. It has been believed that basin-and-swell physiography dominated the East Asia margin and impeded the forming of large rives in the early Cenozoic. In this paper, we combined provenance analysis of East China Sea Basin, where is a crucial place to trace the river evolution in East Asia margin, and regional geologic constraints to reveal drainage reorganizations. Detrital zircon U-Pb ages from the Early Eocene sediments of the East China Sea Basin are firstly reported. Our results together with literature data demonstrate that regional provenance changes occurred at the middle Eocene from one singe age peak at ~110 Ma of proximal sources to multiple age spectrum derived from far inland. Source to sink analysis indicated that the North China Block and Korea Peninsular provided the most detritus. Sedimentation and tectonic features of rift basins in the potential source areas indicated that rivers flowed into Bohai Basin and Jianghan Basin cannot provide terrigenous clasts for the lower reaches in the Eocene. Contrastingly, the dominantly fluvial sediments across the Subei-South Yellow Sea Basin suggested external river system and a bypassing region since the middle Eocene, coinciding with provenance change in ECSB. All these demonstrated that a large river (East Asia River) established in east Asia margin in the middle Eocene and flowed southwestward approximately 1500km to the sea in southern ECSB. This river might last to the middle Miocene. The deep earth processes driven by Izanagi-Pacific ridge subduction resulted in the overfilled stage of Subei-South Yellow Sea Basin and the post-rift subsidence in west depression of ECSB, and thus facilitated the initiation of the EAR. Our finding shed new light on the evolving landscape in East Asia and showed how subduction of deep earth process controlled the initiation of rift-related large rivers.

An increase in air temperature leads to a significant transformation of the relief and landscapes of the Arctic. The rate of permafrost degradation, posing a profound change in the Arctic landscape, depends on air temperature, vegetation cover, type of soils, surface and ground waters. The existing international circumpolar programs dedicated to monitoring the temperature state of permafrost TSP (Thermal State Permafrost) and active layer thickness CALM (Circumpolar Active Layer Monitoring) are not sufficient for a comprehensive characterization of geocryological conditions. Yet, no standardized protocol exists for permafrost monitoring and related processes. Here, we propose a novel multi-parameter monitoring protocol and implement for two sites in the European part of the Russian Arctic: the Yary site along the coast of the Baydaratskaya Bay in the Kara Sea (68.9°N) within the continuous permafrost area and the Hanovey site in the Komi Re-public (67.3°N) within the discontinuous permafrost one. The protocol includes drilling boreholes, determining the composition and properties (vegetation cover and soils), snow cover measure-ment, geophysical imaging, active layer estimation and continuous ground temperature meas-urements. Ground temperature measured in 2014 – 2020 revealed that amplitudes of surface tem-perature fluctuation had no significant difference between the Yary and Hanovey sites, while that the mean annual temperatures between the areas had a considerable difference of greater than 3.0 °C. The period of the presence of the active layer changed with the year (e.g., ranging between 135 and 174 days in the Yary site), showing longer when the air temperatures in summer and the preceding winter were higher. Electrical resistivity tomography (ERT) allowed determining the permafrost distribution and active layer thicknesses. Thermometry results were consistent with our geophysical data. Analyzing the composition and properties of frozen soils helped better interpret the data of geophysical and temperature measurements. Our work by integrating the study of the soil properties, ground temperatures, and ERT allowed fully characterize these sites, suggesting that it helps better understand the thermal state in any other research sites in the European north of Russia. Our suggested monitoring protocol enables calibrating and verifying the numerical and analytical models of the heat transfer through the earth’s surface.

A very strong earthquake of magnitude Mw8.2 struck the coastal zone of Alaska (USA), on July 29, 2021. This earthquake was felt around the Gulf of Alaska, on a wide offshore area belonging to USA and Canada. In order to identify an anomalous geomagnetic signal before the onset of this earthquake, we retrospectively analyzed the data collected on the interval June 17 - July 31, 2021, via internet (www.intermagnet.org), at the two geomagnetic observatories, College (CMO) - Alaska and Newport (NEW)-USA, by using the polarization parameter (BPOL) and the strain effect–related to geomagnetic signal identification. Thus, for the both observation sites (CMO and NEW), the daily mean distribution of the BPOL and its standard deviation (STDEV) are carried out using an FFT band-pass filtering in the ULF range (0.001-0.0083Hz). Further on, a statistical analysis based on a standardized random variable equation was applied to emphasize the following: a) the anomalous signature related to Mw8.2 earthquake on the both time series BPOL*(CMO) and BPOL*(NEW); b) the differentiation of the transient local anomalies associated with Mw8.2 earthquake from the internal and external parts of the geomagnetic field, taking the NEW observatory as reference. Finally, on the BPOL*(NEW-CMO) time series, carried out on the interval 07-31 July, 2021, a very clear anomaly of maximum, greater than 1.2 STDEV, was detected on July 22, with 7 days before the onset of Mw8.2 earthquake

Chilean geography is highly variable, not only from a climatic and hydrological point of view, but also a morphological one, showing unpredictable natural patterns with marked contrasts throughout the country, for which sometimes it is considered as a "crazy" geography. In this paper we have investigated this apparent disorganized character by exploring the fractal properties of fluvial networks extracted from basins distributed across the continental territory. Analytical and semi-empirical methods were applied, finding striking patterns of organization in the distributions of Horton parameters and the fractal dimension of the drainage networks. Fractal dimension reveals to be quite dependent on the drainage area of each unit, showing clear groupings by tectonic and climatological factors. Such dimension reveals to be an important geomorphic parameter, if not the only one able to capture the real morphology of a fluvial network. From our results and despite the diversity of landforms, hydrological, climatic and tectonic conditions, Chilean’s geography is perhaps not as crazy and disorganized as believed.

The distributed acoustic sensing (DAS) has great potential for monitoring natural-resource reservoirs and borehole conditions. However, the large volume of data and complicated wavefield add challenges to processing and interpretation. In this study, we demonstrate that seismic interferometry based on deconvolution is a convenient tool for analyzing this complicated wavefield. We extract coherent wave from the observation of a borehole DAS system at the Brady geothermal field in Nevada. Then, we analyze the coherent reverberating waves, which are used for monitoring temporal changes of the system. These reverberations are tirelessly observed in the vertical borehole DAS data due to cable or casing ringing. The deconvolution method allows us to examine the wavefield at different boundary conditions. We interpret the deconvolved wavefields using a simple 1D string model. The velocity of this wave varies with depth, observation time, temperature, and pressure. We find the velocity is sensitive to disturbances in the borehole related to increasing operation intensity. The velocity decreases with rising temperature, which potentially suggests that the DAS cable or the casing are subjected to high temperature. This reverberation can be decomposed into distinct vibration modes in the spectrum. We find that the wave is dispersive, and the the fundamental mode propagate with a large velocity. The method can be useful for monitoring borehole conditions or reservoir property changes. For the later, we need better coupling than through only friction in the vertical borehole to obtain coherent energy from the formation.

Sinkhole flooding is an essential hydrological process to recharge karst aquifer in arid to dry sub-humid regions. On the other hand, the increase of rain extremes is one of the major consequences of the global warming, together with the expansion of drylands. Thus, appropriate runoff regulation in endorheic karst basins in order to reduce the risk of flooding and improve the quantity and quality of the water drained by sinkholes will be more and more crucial. With these premises, a systematic review was performed by using WoS engine to infer the best practices for the karst water management in regions actually or potentially affected by water scarcity. Hydrological models are essential to manage the consequences of climate change on karst water resource, however the review shows that providing the tools necessary for reliable modeling is still challenging. Finally, due to the intrinsic vulnerability of the karst aquifers, pollution reduction and wastewater recycling policy will play key role in the next decades.

(1) The article is devoted to the development of a theoretical and algorithmic basis for numerical modeling of the spontaneous potential method (SP) as applied to the study of sandy-argillaceous reservoirs. (2) In terms of coupled flows, we consider a physical-mathematical model of SP signals from an electrochemical source, with regard to the case of fluid-saturated shaly sandstone. (3) An algorithm for 2D finite-element modeling of SP signals was developed and implemented in software, along with its internal and external testing with analytical solutions. The numerical SP modeling was carried out, with determining the dependences on the reservoir thickness and porosity, the amount of argillaceous material and the type of minerals. We performed a comparative analysis of the simulated and field SP data, using the results of laboratory core examinations taken from wells in a number of fields in the Latitudinal Ob Region of Western Siberia. (4) The results of the study may be used either for the development of the existing SP techniques, by providing them with a consistent computational model, or for the design of new experimental approaches.

We present here seven new zircon (U-Th)/He (ZHe) ages and three new zircon fission track ages (ZFT) analyzed from an age-elevation profile (Machu Picchu, Peru). ZFT data present older ages in comparison with the other thermochronological data, whereas the ZHe data interestingly present similar ages than the ones obtained with apatite (U-Th)/He (AHe). It has been proposed that He retention in zircon is linked to the damage dose, with an evolution of the closure-temperature from low values associated to low α-dose (<1016 α/g), subsequently increasing before decreasing again at very high α-dose (>1018 α /g). Studies have been focused on the He diffusion behavior at high α-dose, but little is known at low dose. We propose that the ZHe closure temperature at α-dose ranging from 0.6×1015 to 4×1016 α/g is in the range of ~60-80°C. This value is lower than the one proposed in the current damage model ZRDAAM and demonstrates that the ZHe and AHe methods could have similar closure temperatures at low α-dose (i.e. similar ages). These new data strengthen our previous geological conclusions and even highlight an about twice more important cooling rate than the one deduced from AHe and apatite fission-track data alone registered at Machu Picchu.

The contamination of areas around solid urban waste dumps is a global challenge for the maintenance of environmental quality in large urban centres in developing countries. This study applied geophysical methods (electrical resistivity) to identify leachate contamina-tion plumes in the subsoil and groundwater, as well as to describe their temporal (2020 and 2021) dynamics in the lithology and groundwater around the Hulene - B waste dump, Maputo, Mozambique. Geophysical methods (electrical resistivity) were applied to identify possible groundwater contamination plumes, their dynamics, mechanisms of their enrichment and dispersion. Eight 400 m electrical resistivity profiles were performed, four profiles in January 2020 and four profiles in May 2021, overlapped, and the data were inverted with RES2D software. The electrical resistivity models indicate an E - W move-ment of large contamination plumes that dilute superficially into the natural surface wa-ter receiving basin and groundwater, creating zones of resistive anomalies. The thickness of the plumes in the subsurface environment was shown to be extensive in summer for profiles 1a and 2b and we associate it with the higher leachate production and migration mechanisms, which are intense in the hot and rainy season. Profile 4b showed the prop-agation of anomalous surface and subsurface areas, which was associated with higher leachate production and migration process in the new deposition zone (west). The spatial distribution of contamination plumes at both stations reduced significantly as we moved further away from the waste deposit, revealing the attenuating effect of groundwater and lithological substrate (Profile 3 a, b, and fig.7).

The origin of the man at the Earth is directly associated with the determination of directions of the flow distribution of the ancient man migration to adjacent territories. In such studies traditionally mainly landscape and climatological changes are considered. We suggest that along with the above factors, regional tectonic-geodynamic factors played a dominant role in the character of migration. The considered African-Levantine-Caucasian region is one of the most complex regions of the world, where collisional and spreading processes of geodynamics converge. First is determined an essential influence of the Akchagylian hydrospheric maximum (about 200 m above the mean sea level) limiting the ancient men migration from Africa to Eurasia. We propose that the Levantine Corridor emerged after the end of the Akchagylian transgression and landscape forming in the Eastern Mediterranean. This corridor location was formed by the movements between the Dead Sea Transform and the boundary of the carbonate platform of the Mesozoic Terrane Belt. Further landscape evolution was largely determined by the geodynamic behavior of the deep mantle rotating structure occurring below the central part of the region under study. All the mentioned events around and in the Levantine Corridor have been studied in detail on the basis of the combined geodynamic, paleogeographic, and paleomagnetic analyzes performed in northern Israel (Carmel uplift and Galilee plateau). Careful studies of the Evron quarry geological section indicate that it is a unique one for the dating of the marine and continental archaeological sequences and sheds light on the movement of the ancient man along the Levantine Corridor.

The article is devoted to the topical problem of estimating water content in water-oil mixtures and porous media they saturate, according to low-field NMR relaxometry and dielectric spectroscopy. The aim of the research is to theoretically substantiate and experimentally validate the capability of joint interpretation of data from these methods to acquire information on the filtration-volumetric properties of drill cuttings, relaxation characteristics of oil-containing fluids, water/oil ratio in water-oil mixtures and saturated with them drill cuttings in order to control the composition of liquids produced from boreholes. The studies were carried out on samples of cuttings and oils taken from fields in the northern and Arctic regions of the West Siberian oil-and-gas province. Based on the experimental data obtained, we evaluated the water content in the water-oil mixtures, determined the main NMR parameters of the mixtures in terms of properties of the constituent oils, and specified the parameters and shapes of NMR and complex dielectric permittivity spectra. The NMR method was found to be effective in examining high-viscosity and medium-viscosity oils, while the dielectric spectroscopy method – in the study of light oils; their integration allows obtaining reliable data for all the samples under study. We also showed how the shapes of NMR and complex dielectric permittivity spectra depend on the rheological properties of oil belonging to the mixture.

Detailed mapping and structural observations on the Cycladic Blueschist Unit (CBU) of Iraklia island integrated by detrital zircon (DZ) U-Pb ages elucidate the Mesozoic pre-subduction evolution and the Cenozoic orogenic events. Field data reveal that the Iraklia tectonostratigraphy includes a heterogeneous Lower Schist Unit juxtaposed against a Variegated Marble Unit and an overlying Upper Schist Unit. The contact is an extensional ductile-to-brittle-ductile, top-to-N shear zone, associated with the Oligo-Miocene exhumation. The DZ spectrum of the Lower Schist Unit characterized by Gondwanan/peri-Gondwanan provenance signatures points to Late Triassic maximum depositional ages (MDAs). A quartz-rich schist layer yielded Precambrian DZ ages exclusively, considered part of the pre-Variscan metasedimentary Cycladic Basement, equivalent to those observed on Ios island. A significant change occurred during the deposition of the Upper Schist Unit, revealing Late Cretaceous MDAs and a high amount of Late Paleozoic DZ ages, attesting to more internal Pelagonian source areas. The imprint of Paleotethyan vs. Neotethyan geodynamic events is revealed in the DZ U-Pb ages record. The Triassic DZ input demonstrates eroded volcanic material related to the final Paleotethys closure and the Pindos/CBU rift basin opening. Late Cretaceous metamorphic/magmatic zircons and ~48-56 Ma zircon rims constrain the onset of subduction and high-pressure metamorphism.

Tephra plumes can cause a significant hazard for surrounding towns, infrastructure, and air traffic. The current work presents the use of a small and compact X-band Multi-Parameter (X-MP) radar for the remote tephra detection and tracking of two eruptive events at Merapi Volcano, Indonesia, in May and June 2018. Tephra detection was done by analysing the multiple parameters of radar: copolar correlation and reflectivity intensity. These parameters were used to cancel unwanted clutter and retrieve tephra properties, which are grain size and concentration. Real-time spatial and temporal forecasting of tephra dispersal was performed by applying an advection scheme (nowcasting) in the manner of Ensemble Prediction System (EPS). Cross-validation was done using field-survey data, radar observations, and Himawari-8 imagery. The nowcasting model computed both the displacement and growth and decaying rate of the plume based on the temporal changes in two-dimensional movement and tephra concentration, respectively. Our results with ground-based data, where the radar-based estimated grain size distribution fell within the range of in-situ data. The uncertainty of real-time forecasted tephra plume depends on the initial condition, which affects the growth-and decaying rate estimation. The EPS improves the predictability rate by reducing the number of missed and false forecasted events. Our findings and the method presented here are suitable for early warning of tephra fall hazard at the local scale.

In response to the COVID-19 pandemic, and resultant cancellation of geoscience fieldwork, as well as outstanding accessibility issues inherent in conducting fieldwork, we developed a virtual geological fieldtrip (VFT) to the Huronian age deposits in the Whitefish Falls area, Ontario, Canada. This region is a geologically significant site in which many Ontario universities conduct undergraduate teaching due to the high-quality exposures. In this contribution we describe and comment on the development of this openly available resource, the motivations in doing so, the challenges faced, its pedagogical impact and relevance, as well as provide suggestions to others in the development of such resources. Our multimedia VFT combines 360° imagery, georeferenced data on integrated maps, and multi-scale imagery (aerial/drone, outcrop, and thin section images). The VFT was built using the Esri Storymaps platform, and thus offers us the opportunity to review the effectiveness of building such resources using this medium, as well as our approach to doing so. We conclude that the Esri Storymaps platform provides a sound medium for the dissemination of multimedia VFTs, but that some aspects of in-person fieldwork remain hard to replicate. Most notably, ‘hands on experience’ and specific activities such as geological mapping. In addition, while VFTs alleviate some accessibility barriers to geoscience fieldwork, substantial barriers remain that should remain the focus of both pedagogical and geoscience work.

We simulate the effects of diagenesis, cementation and compaction on the elastic properties of shales and sandstones with four different petro-elastical theories and a basin-evolution model, based on constant heating and sedimentation rates. We consider shales composed of clay minerals, mainly smectite and illite, depending on the burial depth, and the pore space is assumed to be saturated with water at hydrostatic conditions. Diagenesis in shale (smectite/illite transformation here) as a function of depth is described by a 5th-order kinetic equation, based on an Arrhenius reaction rate. On the other hand, quartz cementation in sandstones is based on a model that estimates the volume of precipitated quartz cement and the resulting porosity loss from the temperature history, using an equation relating the precipitation rate to temperature. Effective pressure effects (additional compaction) are accounted for by using Athy equation and the Hertz-Mindlin model. The petro-elastic models yield similar seismic velocities, despite the different level of complexity and physics approaches, with increasing density and seismic velocities as a function of depth. The methodology provides a simple procedure to obtain the velocity of shales and sandstones versus temperature and pressure due to the diagenesis-cementation-compaction process.

Landslides are serious geological hazards that become a disaster worldwide, causing a large number of casualties and economic losses every year. There are many factors affecting landslide susceptibility, such as rainfall, soil and slope. Each of them has an important role in the process of slope losing stability. In this paper, the effects of rainfall intensity, rainfall pattern, slope gradient and soil type on landslide susceptibility are studied. In the process of rainfall-induced landslide, the relevant physical quantities of soil changes continuously. Their values and processes are closely related to the time of landslide occurrence. Hence, the variation of soil volumetric water content, matrix suction, pore water pressure and total stress throughout the rainfall are measured. As the results, soil type, slope gradient and rainfall intensity have a large influence on landslide susceptibility. The occurrence of landslides has a prerequisite that the slope is greater than or equal to 15°. The rainfall intensity needs to be not less than 80 mm/h. The difference of rainfall pattern also affects the landslide susceptibility. The rainfall pattern with rainfall intensity peak at the later stage is more likely to induce landslide. Coarser soils with gravels are prone to landslides when other conditions are the same. Steeper slopes, stronger rainfall, and coarser soils can all increase the amount of sediment yield.

The term deviation frequency denotes the boundary between the variable part of the amplitude and phase scintillation spectrum and the part of uninformative noises. We suggested the concept of the “characteristic deviation frequency” during the observation period which is defined as the most probable value of the deviation frequency under current local conditions. This work is a case study of the characteristic deviation frequency (fd) registered for GPS, GLONASS, Galileo and SBAS signals under quiet and weakly disturbed geomagnetic conditions in April 2021 at the mid-latitude GNSS station. Our results demonstrated that the fd value for all signal components of GPS, GLONASS and GALILEO signals varied within 15-22 Hz. The characteristic deviation frequency was 20 Hz for the mentioned GNSS signals. In difference, the deviation frequency was limited within 13-20 Hz for SBAS with the lower characteristic deviation frequency at 18 Hz. We assume that the concept of the characteristic deviation frequency can be used to determine the optimal sampling rate of the GNSS carrier phase data for the ionospheric studies. The characteristic deviation frequency can also characterize the state of the regular trans-ionospheric radio channel.

The Sanshui Basin (SSB) is located at the northern continental margin of the South China Sea and characterized by a continental rift basin. The bimodal volcanic rocks in SSB record the early Cenozoic magmatic activity in the South China Block, on the magmatic evolution process of bimodal volcanic rocks are poorly understood. Clinopyroxenes in bimodal volcanic rocks in the SSB provide an opportunity to investigate the magma process during magma ascent. We classified nine types of clinopyroxene phenocrysts according to the different compositions and textures types in cogenetic basalt-trachyandesite-comenditic trachyte, the composition of unzoned clinopyroxene have an evolution sequence of diopside- hedenbergite- aegirine with the decrease of Mg#, and the trace element contents of unzoned clinopyroxenes also increase systematically during magma evolution, indicating that the genesis of clinopyroxene dominated by fractional crystallization in a closed magma system; however, the clinopyroxenes with reverse zoning and multiple zoning record the process of magma mixing and recharge indicating an open magma system. Whilst fractional crystallization is the dominated process, magma mixing, recharge, and crystal settling complicate magma evolution. Thermobarometric calculations show that clinopyroxene phenocrysts in bimodal volcanic rocks of SSB are distributed in the whole crust during magma ascent. We have established a magma plumbing system, which provides a new constrain for the complex magmatic evolution history in the SSB by detailed mineral-scale analysis.

Seismic travel time tomography using surface waves is an effective tool for three-dimensional crustal imaging. Historically, these surface waves are the result of active seismic sources or earthquakes. More recently, however, also surface waves retrieved through the application of seismic interferometry are exploited. Conventionally, two-step inversion algorithms are employed to solve the tomographic inverse problem. That is, a first inversion results in frequency-dependent, two-dimensional maps of phase velocity, which then serve as input for a series of independent, one-dimensional frequency-to-depth inversions. As such, a two-dimensional grid of localized depth-dependent velocity profiles are obtained. Stitching these separate profiles together subsequently yields a three-dimensional velocity model. Relatively recently, a one-step three-dimensional non-linear tomographic algorithm has been proposed. The algorithm is rooted in a Bayesian framework using Markov chains with reversible jumps, and is referred to as transdimensional tomography. Specifically, the three-dimensional velocity field is parameterized by means of a polyhedral Voronoi tessellation. In this study, we investigate the potential of this algorithm for the purpose of recovering the three-dimensional surface-wave-velocity structure from ambient noise recorded on and around the Reykjanes Peninsula, southwest Iceland. To that end, we design a number of synthetic tests that take into account the station configuration of the Reykjanes seismic network. We find that the algorithm is able to recover the 3D velocity structure at various scales in areas where station density is high. In addition, we find that the standard deviation on the recovered velocities is low in those regions. At the same time, the velocity structure is less well recovered in parts of the peninsula sampled by fewer stations. This implies that the algorithm successfully adapts model resolution to the density of rays. Also, it adapts model resolution to the amount of noise on the travel times. Because the algorithm is computationally demanding, we modify the algorithm such that computational costs are reduced while sufficiently preserving non-linearity. We conclude that the algorithm can now be applied adequately to travel times extracted from (time-averaged) station-station cross correlations by the Reykjanes seismic network.

The topic of this paper is the airborne evaluation of ICESat-2 Advanced Topographic Laser Altimeter System (ATLAS) measurement capabilities and surface-height-determination over crevassed glacial terrain, with a focus on the geodetical accuracy of geophysical data collected from a helicopter. To obtain surface heights over crevassed and otherwise complex ice surface, ICESat-2 data are analyzed using the density-dimension algorithm for ice surfaces (DDA-ice), which yields surface heights at the nominal 0.7~m along-track spacing of ATLAS data. As the result of an ongoing surge, Negribreen, Svalbard, provided an ideal situation for the validation objectives in 2018 and 2019, because many different crevasse types and morphologically complex ice surfaces existed in close proximity. Airborne geophysical data, including laser altimeter data (profilometer data at 905~nm frequency), differential Global Positioning System (GPS), Inertial Measurement Unit (IMU) data, on-board-time-lapse imagery and photographs, were collected during two campaigns in summers of 2018 and 2019. Airborne experiment setup, geodetical correction and data processing steps are described here. To date, there is relatively little knowledge of the geodetical accuracy that can be obtained from kinematic data collection from a helicopter. Our study finds that (1)~Kinematic GPS data collection with correction in post-processing yields higher accuracies than Real-Time-Kinematic (RTK) data collection. (2)~Processing of only the rover data using the Natural Resources Canada Spatial Reference System Precise Point Positioning (CSRS-PPP) software is sufficiently accurate for the sub-satellite validation purpose. (3)~Distances between ICESat-2 ground tracks and airborne ground tracks were generally better than 25~m, while distance between predicted and actual ICESat-2 ground track was on the order of 9~m, which allows direct comparison of ice-surface heights and spatial statistical characteristics of crevasses from the satellite and airborne measurements. (4)~The Lasertech Universal Laser System (ULS), operated at up to 300~m above ground level, yields full return frequency (400~Hz) and 0.06-0.08~m on-ice along-track spacing of height measurements. (5)~Cross-over differences of airborne laser altimeter data are 0.1918 $\pm$ 2.385~m along straight paths over generally crevassed terrain, which implies a precision of approximately 2.4~m for ICESat-2 validation experiments. (6)~In summary, the comparatively light-weight experiment setup of a suite of small survey equipment mounted on a Eurocopter (Helicopter AS-350) and kinematic GPS data analyzed in post-processing using CSRS-PPP leads to high accuracy repeats of the ICESat-2 tracks. The technical results (1)-(6) indicate that direct comparison of ice-surface heights and crevasse depths from the ICESat-2 and airborne laser altimeter data is warranted. The final result of the validation is that ICESat-2 ATLAS data, analyzed with the DDA-ice, facilitate surface-height determination over crevassed terrain, in good agreement with airborne data, including spatial characteristics, such as surface roughness, crevasse spacing and depth, which are key informants on the deformation and dynamics of a glacier during surge.

The spectra of the plasma emitted from the studied samples consist of several dozens of narrow bands superimposed on each other. Tables of spectral lines were used to interpret the spectra. It turned out that the largest number of bands corresponds to the radiation of positively charged ions and atoms of elements that make up the crystal lattices of minerals that make up the studied rocks. Thus, the spectra of the plasma emitted from quartz corresponded to the radiation of atoms and positively charged silicon ions, the charge of which varied from 1 to 4, as well as atoms and positively charged oxygen ions, the charge of which varied from 1 to 3. Positively charged ions and atoms of Si, O, K, Ca, Al and Na, which are part of the crystal lattices of quartz and feldspar, flew out of granites. Positively charged ions and Ca, C and O atoms flew out of the calcite.

Focusing geophysics to improve recovery factor of hydrocarbon reservoirs adds value and contributes toward ZERO carbon footprint by increasing the recovery factor by of 30-40 % and thus reducing the cost/carbon emission per produced barrel. Thus, the Enhanced Oil Recovery (EOR) market is expected to grow more than 3.5% annually. This will be even more fueled by the Green-House-Gas (GHG) reduction and subsequent CO2 injection into the reservoirs as they are being produced. Presently, geophysics only ac-counts for a small percentage of this market, thus its growth is inevitable since more deterministic observation lead higher operating efficiency. Imaging the fluids (hydro-carbon, water, and CO2) is a key component to optimized production and injection. We designed a novel electromagnetic (EM) acquisition system that combines mul-ti-physics fluid imaging and acquires surface and borehole data with high fidelity. Borehole calibration is needed to upscale reservoir data and parameters to measurement scale. Multiple electromagnetic methods are used as well as microseismics in one layout for Exploration and Production (E & P) problems. Multi-components in electromagnetics allows resolving oil and water-bearing zones equally well while achieving the best ac-curacy suitable for repeat measurements. Because sedimentary basins are intrinsically anisotropic, considering 3-dimensional anisotropy is essential from measurement and 3D modeling viewpoint. Thus, the results have the better subsurface images. Here, we combine hardware design, methodology, 3D modeling, processing, and interpretations into an integrated technology and demonstrate the success with verifiable case histories.

In the aftermath of pyroclastic-flow –dominated eruptions, lahars are the main geomorphic agent, but at the decadal scale, different sets of processes take place in the volcanic sediment cascade. At Unzen Volcano, in the Gokurakudani Gully we investigated the geomorphologic evolution and how the topographic change and the sediment change over time is controlling this transition. For this purpose, a combination of LiDAR data, aerial photography and photogrammetry, Ground Penetrating Radar and sediment grain-size analysis was done. The results show chocking zones and zones of enlargement of the gully, partly controlled by pre-eruption topography, but also by the overlapping patterns of the pyroclastic flow deposits of 1990 – 1995. The Ground Penetrating Radar revealed that on top of the typical lahar structure at the bottom of the gully, side-wall collapses were trapping finer sandy sediments formed in relatively low-energy deposition environment. This shows that secondary processes are taking place in the sediment transport process, on top of lahar activity, but also that these temporary dams may be a source of sudden sediment and water release, leading to lahars. Finally, the sediments from the gully walls are being preferentially oozed out of the pyroclastic-flow deposit, meaning that over longer period of time, there may be a lack of fines, increasing permeability and reducing internal pore-pressure needed for lahar triggering. It also poses the important question of how much of a past-event one can understand from outcrops in coarse heterometric material, as the deposit structure can remain, even after loosing part of its fine material.

MzSTools is a plugin for QGIS developed by the Institute of Environmental Geology and Geoengineering of the National Research Council (CNR-IGAG). The plugin has been designed as a set of practical and easy-to-use tools to carry out seismic microzonation (SM) studies, by producing standard compliant geographic database and maps, thus making them accurate, homogeneous and uniform for all municipalities in Italy. A geodatabase based on SQLite/SpatiaLite Relational Database Management System (RDBMS) has been designed to collect and store data related to elements such as: geognostic surveys; bedrocks and cover terrains; superficial and buried geomorphological elements; tectonic-structural elements; elements of geological instability such as landslide zones, liquefaction zones and zones affected by active and capable faults; homogeneous microzones in seismic perspective, microzones characterized by a seismic amplification factor. MzSTools assembles in a single software environment a set of useful tools in a configurable QGIS project template, comprising layers, symbol libraries, cartographic styles and print layouts for the SM maps. The plugin is open source and hosted on the GitHub platform, and available via the official QGIS plugins repository (https://plugins.qgis.org/plugins/MzSTools/).

Monitoring of groundwater quality in today's scenario is very much important. Due to urbanization and population pressure regular monitoring of groundwater for drinking as well as irrigation purposes need a major concern. With this aim, a study has been carried out consisting 26 groundwater samples in May 2017, to access the physiochemical characteristic, water quality index (WQI) of groundwater by using GIS software and to find out the groundwater suitableness for drinking as well as for irrigation purpose. The pH is slightly alkaline and the TDS is much more than prescribed limits of BIS. The trend of cations in groundwater are Ca2+>Na+>Mg2+>K+ while anions trend is HCO3->SO42->Cl->NO3->CO32->F-. The Ca-Mg-HCO3 and Na-K-Cl-SO4 types of groundwater facies were dominant. Generally, the chemical changes in groundwater are administered by the evaporation process with ion exchange, and mixing of particles is the significant source of the solute acquisition process. WQI of the study area suggested that the 15% sample is unsuitable, 69% is poor and remaining is good for drinking uses. The potential salinity of the groundwater sample is nearly high although the majority of the sample is suited for irrigation activities.

Temperature-index modeling is used to determine the magnitude of temperature depression on the Blanca Massif, Colorado, required to maintain steady-state mass balances of nine reconstructed glaciers at their extent during the Last Glacial Maximum (LGM). The mean temperature depression thus determined is ~8.6 +0.7/–0.9 °C where the uncertainties account for those inherent in the glacier reconstructions, in model parameters (e.g. melt factors), and possible modest changes in LGM precipitation. Associated equilibrium-line altitudes (ELAs) exhibit a statistically significant directional dependency being lower toward the north and east. Under the assumption that regional temperature change was uniform, required changes in precipitation vary systematically – also exhibiting a directional dependency coinciding with that in ELAs – and indicate increases (over modern) occurred on the eastern side of the massif while decreases occurred on the western side. This disparity represents a strengthening of a precipitation asymmetry, particularly winter precipitation, that exists today. The modern precipitation asymmetry may be a consequence of snow being blown over to the eastern side of the massif (advective transport) by southwesterly flow. Intensification of this flow during the LGM would have enhanced advection, and augmented snow accumulation on glaciers, thus explaining the lower ELAs and increased precipitation on that side of the massif.

With global warming and accelerated degradation of permafrost, the engineering problems caused by the formation of weak zones between the shallow and permafrost layers of soil–rock mixture (S-RM) slopes in permafrost regions have become increasingly prominent. To explore the influence of rock content on the shear strength of the S-RM freezing–thawing interface, the variation in the shear strength for different rock content is studied herein using direct shear tests. In addition, a 3D laser scanner is used for obtaining the topography of the shear failure surface. Combined with the analysis results of the shear band-particle calculation model, the influence of the rock content on the shear strength of the interface is explored. It was found that the impact threshold of the rock content on the interface strength and failure mode is approximately 30%, when the rock content (R) is &gt; 30% and that the shear strength increases rapidly with increasing rock content. When R ≤ 30%, the actual shear plane is similar to waves; when R &gt; 30%, the shear plane appears as gnawing failure. The shear strength of S-RM freezing–thawing interface mainly comes from the bite force and friction between particles. The main reason for the increase in shear strength with increasing rock content is the increase in bite force between particles, which makes the ratio of bite force to friction force approximately 1:1.

We confirm that energy dissipation weighting provides the most accurate approach to determining the effective hydraulic conductivity (Keff) of a binary K grid. A deep learning algorithm (UNET) can infer Keff with extremely high accuracy (R2 > 0.99). The UNET architecture could be trained to infer the energy dissipation weighting pattern from an image of the K distribution with high fidelity, although it was less accurate for cases with highly localized structures that controlled flow. Furthermore, the UNET architecture learned to infer the energy dissipation weighting even if it was not trained on this information directly. However, the weights were represented within the UNET in a way that was not immediately interpretable by a human user. This reiterates the idea that even if ML/DL algorithms are trained to make some hydrologic predictions accurately, they must be designed and trained to provide each user-required output if their results are to be used to improve our understanding of hydrologic systems most effectively.

The paper proposes a novel methodology for the stabilization of shallow foundations, with a simplified model combined with 3D Electrical resistivity tomography (ERT-3D) and conso- lidation injections. To determine its usefulness, the method has been applied in a case located in Estepona (southern Spain). The chosen tomography model is the dipole-dipole configuration, with an optimized distance between electrodes of 0.80 m for a better visualization of the foundation subsoil; with this parameterization, a total of 72 electrodes were installed in the analyzed case. In this work, the depth of the anomaly in the building's supporting subsoil was detected ranging from 2.00 m to 3.90 m deep. The study also delineates areas of high resistivity variations (50-1,000 Ω m) in the middle and eastern end of the field. These data have been validated and corroborated with a field campaign. The results of the ERT-3D monitoring are presented, once the investment data has been processed with the RES3DINV software, from the beginning to the end of the stabilization intervention. The novelty occurs with the interaction between the tomography and the foundation consolidation injections, until the final stabilization; very useful methodology in case of emergency consolidation, where there is a need to minimize damage to the building. Thus, people using this combined system; will be able to practically solve the initial anomalies of the subsoil that caused the damages, in a non-invasive way, considerably lowering the value of the resistivities.

The Ptarmigan and Tom mesothermal gold deposits are located 10 km to the northeast of the city of Yellowknife, Northwest Territories in northern Canada. Both gold deposits comprise a series of en echelon veins that are hosted within upper greenschist to lower amphibolite facies ~2630 Ma (peak) rocks. Supracrustal units across the craton are intruded by the ca. 2610–2605-Ma granodiorite, tonalite, monzodiorite, quartz diorite, and affiliated rocks of the Concession Suite. Hydrothermal apatite is a common accessory mineral in both mineralized and non-mineralized quartz veins in the metasedimentary host rocks that constitute the Ptarmigan and Tom deposits. This study characterizes and compares turbidite-hosted hydrothermal apatite from the Ptarmigan and Tom deposits, non-mineralized veins adjacent to the ore body, and magmatic apatite from proximal LCT-pegmatites. Using electron probe microanalyses (EPMA), laser ablation inductively coupled plasma-mass spectrometry (LA-ICP-MS), micro-XRF, and cathodoluminescence (CL), the major, minor, and trace element abundances have been quantified and mapped. In addition to utilizing this data to determine if the chemistry of apatite can be used to constrain the source of hydrothermal fluids, the apparent age of the apatite is also evaluated utilizing in situ U-Pb dating. The distribution and abundance of major, minor, and trace elements from in situ recovered apatite were studied to characterize the nature of mineralizing fluids. Most apatite from mineralized and non-mineralized veins show different Mn, Sr, and Pb contents, as well as chondrite-normalized rare-earth element (REE) and Y abundance patterns. REEs display five unique chondrite-normalized patterns: (1) negative sloped pattern with slight negative Eu anomaly, (2) a flat pattern with a positive Eu anomaly, (3) a positive slope with a negative Eu anomaly, (iv) light rare earth element (LREE) depleted pattern with positive Eu anomaly, and (v) bell-shaped pattern with a negative Eu anomaly. The REE patterns reflect both the source of the auriferous hydrothermal fluids and, perhaps, co-precipitating mineral phases. Apatite from the Ptarmigan vein occurs with both: (1) a flat pattern with a positive Eu anomaly and (2) bell-shaped pattern with a negative Eu anomaly. The bell-shaped and flat patterns typify orogenic gold deposits. Vein-hosted apatite commonly displays compositional zoning with a characteristic yellow cathodoluminescence (CL) emission spectra with darker cores and brighter rims. The cores have lower REE, whereas the rims are notably higher in REE. It is thought that the darker cores in CL images reflect a transition from an early low REE hydrothermal fluid to one enriched in REE. The hydrothermal apatite age of 2585 ± 15 Ma is consistent with the intrusions of the 2605 and 2590 Ma two-mica granites of the Prosperous Suite and associated LCT pegmatites.

Mediterranean Ridge accretionary complex (MAC) is one of the most critical subduction zones in the world. It is known that the region exhibits a continuous mass change (horizontal/vertical movements). This process is associated with the devastating and tragic earthquakes shaking the MAC for centuries. Here, we investigate the ocean bottom pressure (OBP) anomalies in the MAC derived from the Gravity Recovery and Climate Experiment (GRACE) and GRACE Follow On (GRACE-FO) satellite missions. The OBP time series for the MAC comprises a decreasing trend in addition to 1-, 1.53-, 2.36-, 3.67-, and 9.17-year periodic components partially explained by the atmosphere, oceans, and hydrosphere (AOH) processes, and Earth's pole movement. We noticed that the OBP anomalies appear to link to a rising trend and periods in earthquakes' power time series. This finding sheds new light on the mechanisms controlling the most destructive natural hazard.

An earthquake of Mw6.4 hit the coastal zone of Albania on 26 November 2019, at 02:54:11 UTC. It was intensively felt at about 34km far away, in Tirana City, where damages and lives lost occurred. To emphasize a geomagnetic signature before the onset of this earthquake, the data collected on the interval 15 October–30 November 2019, at the Panagjurishte (PAG)-Bulgaria and Surlari (SUA)-Romania observatories are analyzed by using both the polarization parameter (BPOL)-time invariant in non-seismic conditions, becoming unstable before this seismic event, and the strain effect for geomagnetic signal identification. Consequently, BPOL time series and its standard deviations are performed for the both sites using ULF-FFT band-pass filtering. A statistical analysis, based on a standardized random variable equation, was applied to emphasize on the BPOL*(PAG) and ABS BPOL*(PAG) time series the anomalous signal’s singularity and, to differentiate the transient local anomalies due to Mw6.4earthquake, from the internal and external parts of the geomagnetic field, taken PAG observatory as reference. Finally, the ABS BPOL*(PAG-SUA) time series are obtained on the interval 1-30 November, 2019, where a geomagnetic signature greater than 2.0, was detected on 23 November and the lead time was 3 days before the onset of Mw6.4earthquake.