The goal of the present study is to investigate, delineate, and evaluate the shallow Miocene groundwater aquifer with vertical and lateral facies variations in the Gara oasis, western desert of Egypt. This oasis represents a typical arid region in North Africa. Through grid texture analysis, lineament detection, edge detection, thresholding, and identifying areas of structural complexity from the filtered residual reduced-to-pole magnetic anomaly map, it is possible to outline the fracture zones that principally control the groundwater aquifers and water flow in the area. The groundwater quality and quantity are examined hydro-geochemically through nine groundwater samples that were gathered from wells and springs distributed throughout the area around Gara Lake. Measurements of physio-chemical parameters (TDS, pH, and EC) are carried out and its spatial distribution is critically studied. The results reveal that the ion-exchange process caused by water-rock interaction is the dominant process. Furthermore, the main ions in the groundwater in the study areas were Na and Cl. This might be due to evaporation or halite dissolution with the upward flow of waters through the fractures from the deep aquifer of the Nubian sandstone to the shallow aquifer of the fractured Miocene limestone.

In this paper, the optimization algorithm based on the population as improved global particle swarm optimization is described and used for inverse modelling of two-dimensional magnetic field data. This algorithm is able to estimate the parameters of depth, shape factor, amplitude coefficient, magnetic inclination angle and origin point coordinates. To evaluate the efficiency of this method, the magnetic field of an artificial model was analysed, with and without added random noise. The results suggest that the proposed algorithm is capable of model parameter estimation with high accuracy. Accordingly, the improved global particle swarm optimization algorithm was used to analyse the magnetic field of the study area in the Ileh region in Iran located in Taybad city. The study area is very rich in terms of iron resources. The estimate for the study area is that the depth of the buried mass centre is about 114.9 m and its approximate shape is similar to a horizontal cylinder based on the calculated shape factor value which is 1.76. The calculated depth is an acceptable match with the average depth of drillings.

This contribution presents a new velocity field describing the crustal motion of Nigeria from more than 5years continuous GNSS data at 8 of 14 permanent stations distributed across the country. GAMIT/GLOBK GNSS processing software was used for processing. The horizontal velocity field of NigNET stations which showed a North-East trend and the selected IGS stations were obtained from cleaned position time series of daily GNSS solutions. The velocity describe the horizontal and vertical motion the selected 8 GNSS stations assuming Flicker +White noise model, which optimally describes the geophysical error source of the adopted GNSS stations. Bland and Altman statistical method shows that residual velocity solutions of our study and others (in ITRF2008) and also MORVEL plate motion model are in agreement.

Groundwater is a significant resource that supports almost one-fifth population globally, but has been is diminishing at an alarming rate in recent years. To delve into this objective more thoroughly, we calculated interannual (2002–2020) GWS (per grid) distribution using GRACE & GRACE-FO (CSR-M, JPL-M and SH) Level 3 RL06 datasets in seven Indian river basins and found comparatively higher negative trends (−20.10 ± 1.81 to −8.60 ± 1.52 mm/yr) in Basin 1–4 than in Basin 5–7 (−7.11 ± 0.64 to −0.76 ± 0.47 mm/yr). After comparing the Groundwater Storage (GWS) results with the CHIRPS (Climate Hazards Group Infrared Precipitation with Stations) derived SPI (Standardized Precipitation Index) drought index, we found that GWS exhausts analogously in the same period (2005–2020) when SPI values show improvement (∼1.89–2), indicating towards wet condition. Subsequently, the GWSA time series is decomposed using the STL (Seasonal Trend Decomposition) (LOESS Regression) approach to monitor long-term groundwater fluctuation. The long term GWS rate (mm/yr) derived from three GRACE & GRACE-FO solutions vary from −20.3 ± 5.52 to −13.19 ± 3.28 and the GWS mass rate (km3/yr) lie in range of −15.17 ± 4.18 to −1.67 ± 0.49 for basins 1–3. Simultaneously, in basin 4–7 the GWS rate observed is −8.56 ± 8.03 to −0.58 ± 7.04 mm/yr, and the GWS mass rate differs by −1.71 ± 0.64$ to −0.26 ± 3.19 km3/yr. The deseasonalized GWS estimation (2002–2020) states that Himalayan River basins 1,2,3 exhibit high GWS mass loss (−260 to −35.12 km3), with Basin 2 being the highest (−260 km3). Whereas the Peninsular River basin 4,6,7 gives moderate mass loss value from −26.72 to −23.58 km3. And in River basin 5, the GWS mass loss observed is the lowest, with a value of −8 km3. 

We analyse spatiotemporal gravity changes observed on the Ischia island (Italy) accompanying the destructive earthquake of 21 August 2017. The 29 May 2016 to 22 September 2017 time-lapse gravity changes observed at 18 benchmarks of the Ischia gravimetric network are first corrected for the gravitational effect of the surface deformation using the deformation-induced topographic effect (DITE) correction. The co-seismic DITE is computed by Newtonian volumetric integration using the Toposk software, a high-resolution LiDAR DEM and the co-seismic vertical displacement field derived from Sentinel-1 InSAR data. We compare numerically the DITE field with its commonly used Bouguer approximation over the island of Ischia with the outcome that the Bouguer approximation of DITE is adequate and accurate in this case. The residual gravity changes are then computed at gravity benchmarks by correcting the observed gravity changes for the planar Bouguer effect of the elevation changes at benchmarks over the same period. The residual gravity changes are then inverted using an inversion approach based on model exploration and growing source bodies, making use of the Growth-dg inversion tool. The found inversion model, given as subsurface time-lapse density changes, is then interpreted as mainly due to a co-seismic or post-seismic disturbance of the hydrothermal system of the island. Pros and weak points of such interpretation are discussed.

To understand the variation of stress levels in the region 80°E – 89°E and 26°N – 31°N, the statistical analysis of earthquake frequency-magnitude distribution and spatio-temporal variation of fractal correlation dimension of earthquake epicenter distribution are estimated. The analysis is carried out on declusterised catalogue containing 1185 events of 56 years from February 1964 to November 2020. The study area is divided into three regions the western Nepal and vicinity (Region A), central Nepal and vicinity (Region B) and eastern Nepal and vicinity (Region C), respectively. The magnitude of completeness (Mc) varies from 3.6 to 4.0 for the study period. The spatial fractal dimension (Dc) and b-value are calculated as 1.89 ± 0.02 and 0.68 ± 0.03 for the western Nepal, 1.76 ± 0.01 and 0.60 ± 0.05 for the central Nepal, whereas they are estimated as 1.85 ± 0.02 and 0.63 ± 0.03 for the eastern part of the Nepal. The b-values obtained for all three regions are very low comparing to global average value of 1. The time clustering of the events in the respective regions are 0.26 ± 0.003, 0.31 ± 0.004 and 0.26 ± 0.02 as indicated by temporal fractal dimension (Dt). The higher Dc, lower b and Dt values associated with the regions indicate high stress concentration and stronger epicenter clustering in these regions. The strongly increasing trend of fractal dimension and strongly decreasing trend of b-value show the high probabilities of occurring the large earthquake in both central Nepal (82.5°E – 85.5°E and 27.5°N – 30°N) and eastern Nepal (85.5°E – 88.2°E and 26.45°N – 28.6°N) as compared to western Nepal (80°E – 82.5°E and 28°N – 30.5°N). This statistical analysis of spatial and temporal characteristics of the earthquake activity may give significant signs of the future seismic hazard along central Himalaya region.

The selection of a choice electrode is pertinent to attenuating noise and improving geophysical tomographic inversion results. Besides, the detailed understanding of the geodynamic condition of subsurface formation is crucial to sustainable potable groundwater abstraction. Hence, the subsurface lithostratigraphic units and groundwater potential of two sites (i.e., Site 1 and Site 2) within the Universiti Sains Malaysia were evaluated using borehole-constrained electrical resistivity tomography (ERT) and induced polarisation (IP) tomography. Both methods employed the resolution capacities of stainless-steel and copper electrodes at dual-spacing. The ERT and IP field data and inversion results for copper electrodes were generally robust due to the generated higher positive data points and lower RMS errors, percentage relative differences, and mean absolute percentage errors (MAPE) than the stainless-steel electrodes, especially at Site 1 with a profile length of 200 m and an electrode spacing of 5 m. However, both electrodes tend to produce inversion models with almost the same parameters at Site 2, using half the profile length and electrode spacing of Site 1, i.e., 100 m and 2.5 m, respectively. Thus, the sensitivities and resolution capacities of the tomographic electrodes are heavily influenced by electrode spacing, profile length, amount of injected current, and depth of investigation. The borehole lithostratigraphic units, typically sandy silt, sand, and silty sand, have good correlations with the ERT and IP inversion results. The variability in observed resistivity and chargeability values were due to heterogeneous weathered materials and saturating water fills within the fractured and deeply-weathered granitic bedrock, with <200 Ωm and a chargeability of >1.8 msec. The models' median depth of >40 m mapped for the weathered and/or fractured sections was suggestive of high groundwater-yielding capacity in boreholes to sustain a part of the university community. 

The Zagros suture zone is seismically active region in Iranian plateau. This region is of high importance in terms of seismicity, since it is a vast and populated region and in recent years the earthquakes with high intensities have frequently occurred and have caused extensive destruction and heavy human loss. The study of the focal mechanism is very important in understanding the seismotectonic characteristics. Focal mechanisms of Zagros were collected over a period of 20 years and they were classified by FMC software. Seven groups were considered for the type of faulting and Zagros was divided into three zones. For each zone, the frequency percentage of each group of faults was determined. The most of faulting are of the reverse and compression type with the strike-slip component. Finally, the role of nodal plane selection in determining the type of faulting was discussed and it was found that the selection of each nodal plane in determining the type of faulting has the same result.

Al-Gor area is a part of Southwestern Sinai of Egypt. It is considered as one of the most promising areas for mineralization in Egypt, being rich in many mineral deposits of: manganese, iron, copper, zinc, lead, cobalt, nickel, silver, gibbsite, and uranium. Besides, some industrial ore minerals such as kaolin and glass sand,… etc. are also found in this area. The area was studied by Gama-ray spectrometry to trace the radioactive anomalies, their concentrations and their relationship to the existing rocks, and by magnetic survey to study the relationship of radioactive anomalies and their trends with the trends of geological structures. The gamma-ray spectrometric maps show different levels over the surveyed area, which reflect contrasting radioelement contents for the exposed various rock types. The highest radiospectrometric levels are located in the northwest southeast direction and some scattered parts all over the study area. They are mainly associated with Um-Bogma Formation, bearing gibbsite. The study area possesses radiospectrometric ranging between 0.6 and 110.9 Ur as a total-count, 0.1 to 1.8 % for K, 0.1 to 99 ppm for eU and 0.1 to 23 ppm for eTh. The qualitative analyses of magnetic data show the existences of a number of different magnetic anomalies, with different amplitudes and frequencies as well as trends. From the application of spectral analyses of magnetic data, the regional and residual depths of magnetic anomalies can be computed. The first depth represents the regional (deep-seated) anomalies, at about 75 m and the residual (shallow-seated) anomalies, at about 20 m. The trends of the structures as derived from the spectrometric and ratio maps correspond to those inferred from the residual-component magnetic map, which reflects the effect of structures on the concentration of radioactive elements and, consequently mineralization.

In this contribution we present results from a case-study, which was performed in collaboration between geophysicists and explosive ordnance disposal technicians at the Rohožník military training range in SW Slovakia. The aim of this study was to locate a deep-penetrated unexploded Mk-82 aerial bomb using high-definition digital magnetometry. The location where this bomb had entered the ground was known but its final position needed to be determined so that a safe excavation and disposal could be conducted. However, the detection of this unexploded ordnance object was complicated by the presence of intense magnetic interference from a number of near surface ferrous items including non-explosive test bombs, fragmentation and other iron junk. These items contributed a localised, high amplitude of magnetic clutter masking any deeper source. Our strategy was to approach the problem in three stages. First, we used magnetic data to locate the near surface items. After the detection and before the excavation of the searched objects, two quantitative interpretation methods were used. These involved an optimised modelling of source bodies and the application of a 3D Euler deconvolution. Both methods yielded acceptable results, but the former was found to be more accurate. After the interpretation phase, many of the items were then safely excavated and removed individually. A second magnetic mapping was then performed and from this data which was now significantly less cluttered, we were able to identify but not quantify, two deep source items and to confirm that all remaining near surface items were significantly smaller in size than a Mk-82 bomb. As the remaining near surface sources were interpreted as being contained within the surface one metre of soil and being small they could be assured to be non-explosive, it was considered most practical to mechanically excavate and remove this soil and the remaining objects contained.

GPS measurements have proved extremely useful in quantifying strain accumulation rate and assessing seismic hazard in a region. Continuous GPS measurements provide estimates of secular motion used to understand the earthquake and other geodynamic processes. GNSS stations extending from the South of India to the Higher Himalayan region have been used to quantify the strain build-up rate in Central India and the Himalayan region to assess the seismic hazard potential in this realm. Velocity solution has been determined after the application of Markov noise estimated from GPS time series data. The recorded GPS data are processed along with the closest International GNSS stations data for estimation of daily basis precise positioning. The baseline method has been used for the estimation of the linear strain rate between the two stations. Whereas the principal strain axes, maximum shear strain, rotation rate, and crustal shortening rate has been calculated through the site velocity using an independent approach; least-square inversion approach-based triangulation method. The strain rate analysis estimated by the triangulation approach exhibits a mean value of extension rate of 26.08 nano-strain/yr towards N131°, the compression rate of –25.38 nano-strain/yr towards N41°, maximum shear strain rate of 51.47 nano-strain/yr, dilation of –37.57 nano-strain/yr and rotation rate of 0.7°/Ma towards anti-clockwise. The computed strain rate from the Baseline method and the Triangulation method reports an extensive compression rate that gradually increases from the Indo-Gangetic Plain in South to Higher Himalaya in North. The slip deficit rate between India and Eurasia Plate in Kumaun Garhwal Himalaya has been computed as 18±1.5 mm/yr based on elastic dislocation theory. Thus, in this study, present-day surface deformation rate and interseismic strain accumulation rate in the Himalayan region and the Central Indian region have been estimated for seismic hazard analysis using continuous GPS measurements.

We used 2D integrated geophysical modelling approach to calculate the temperature distribution in the lithosphere along profile VII passing through the Eastern Carpathians. With assigned rheological parameters of rocks and obtained temperature field, we derived the rheological model of the lithosphere along the studied profile. We have calculated the strength distribution in the lithosphere, based on the brittle and ductile deformation, for compressional and extensional regimes and the vertically integrated strength along the profile. To illustrate the strength distribution in different tectonic units, we have calculated the yield strength envelopes for chosen lithospheric columns. Ours results show that the dominant regime is compressional and the largest strength occurs on the boundary between the upper and lower crust. Along the studied profile, the strength decreases from its high values in the European platform towards its minimum in the Trans-European Suture Zone (TESZ). In the Eastern Carpathians, the strength increases, reaches two maxima, the first in the Outer Eastern Carpathians, and the second in the Inner Eastern Carpathians, where the highest values of strength can be observed. Another local maximum along the profile can be observed in the Apuseni Mountains, while the minimal strength is observed in the Transylvanian Basin. The diverse rheological behaviour of studied tectonic units seems to be in accordance with their lithospheric structure and tectonics.

Gravitational effect of surface deformation is in 4D microgravimetry treated as the deformation-induced topographic effect (DITE). The DITE field is computed using Newtonian volumetric integration which requires high resolution digital elevation model (DEM) and vertical displacement field in areal form. If only elevation changes on benchmarks of the gravimetric network are available, instead of the vertical displacement field, the DITE on benchmarks can be evaluated only approximately, using a planar Bouguer or a normal free-air-effect (nFAE) approximation. Here we analyze the adequacy and accuracy of these two approximations in a case study for the December 2018 fissure eruption on Etna accompanied by significant surface deformation caused primarily by a relatively shallow dyke. The outcome is that in volcanic areas of similar morphology as that over the Etna summit area, and for surface deformation fields due to relatively shallow dykes, neither the Bouguer nor the nFAE approximation of the DITE is accurate enough. In such situations the residual gravity changes should be computed with both the Bouguer and nFAE corrections and interpreted as two marginal cases. In addition we analyze also a correction for the effect of benchmark elevation change based on the topographically modelled (predicted) vertical gradient of gravity (VGG) meant to approximate the in-situ VGG values at benchmarks. This correction does not appear suitable to approximate the DITE in conditions of our case study or in broader sense.

Among important parameters in simulation of earthquake data in high frequencies are the high frequency spectral amplitude decay and the Quality factor. Amplitude spectral decay is determined by the Kappa parameter (K) and the Quality factor (Q) which is usually expressed by a power relation of frequency (f) as Q = Q0 f n, where Q0 is Q at 1 Hz. The 2017 Sarpol-e-Zahab earthquake with magnitude Mw = 7.3 in Kermanshah province near the Iran-Iraq border caused extensive destruction and heavy human loss. Thus, the study of different aspects of this event is of high importance. In this paper an attempt is made to partly explain the attenuation properties of this region in Zagros suture zone by determining the Kappa and the Quality factors in this region. In this study, accelerograph records of aftershocks of the above-mentioned earthquake have been analysed. The best linear fit for the Kappa, based on the distance (R) in km, is estimated as: K = 0.0005 R + 0.034 for the horizontal component, which exhibits increase with increasing epicentral distance. The correlation of the Quality factor was also found as Q = 88.6 f 0.8, which is in accordance with an active tectonic region.

Hydrometeors (rain, fog and ice crystals) affect the transmission of electromagnetic signals. Previous research showed that alterations in the signal (amplitude and phase) are affected by the composition of the atmosphere, e.g. the presence of hydrometeors. The majority of hydrometeorological detecting methods are based on the attenuation of electromagnetic signals as they penetrate the atmosphere. Novel methods based on monitoring of parameters of the signal appeared in recent time. This article presents the first results from our investigation of how hydrometeors affect the phase differences in signals transmitted by BTS stations. Cell phone operators transmit electromagnetic signals in the 1 GHz frequency band. This paper describes a novel concept of how phase differences between two signals arriving at two different antennas can be used to detect hydrometeors. Although the described concept is assumed to be independent from the signal strength, the analysed signal must be detectable. The primary advantage of the proposed passive method is that the signal is almost ubiquitous and does not require demodulation. In densely populated areas, the network of BTS stations reaches a spatial density of 1 station per 1 km2 which gives excellent opportunity to use the signal for detection purposes.

Tuticorin, located in the southeastern part of Tamil Nadu, is a coastal aquifer on which Vertical Electrical Sounding (VES) was conducted, thereby analysing the interpreted subsurface resistivity layer parameters covering around 112 km2 area of the study region. VES is an essential tool for investigating hard rock terrains of coastal aquifers and perceive an idea about the groundwater quality. In this study, Dar-Zarrouk (D-Z) parameters like longitudinal conductance (Sc), transverse resistance (Tr) and anisotropy (λ) are analysed as these are well-established parameters in delineating the occurrence and distribution of both fresh and saline water aquifers. These parameters are also very persuasive in investigating complex subsurface parameters (resistivity and conductivity) within saline water intruded coastal region environment. After conducting a thorough survey, the resistivity results reflect that the sediments are enriched with saltwater, clay with moderate freshwater and freshwater-bearing formations. The analysis shows that the D-Z parameters offer a helpful and assured answer in demarcating the saline, moderate fresh, and freshwater aquifers. Therefore, the behaviour and patterns of the D-Z parameters in space established the existence of saline water and freshwater aquifer structures in the coastal aquifers over a vast area.

Several surface geophysical methods were used to study the displays of undermined spaces in their physical fields and to try to distinguish undermined from non-undermined parts of underground mining area. The studied area of the Čáry lignite mine in the Western Slovakia represents an actively subsided place with high risk to the population. Despite very low radioactivity and relatively high gas permeability of building geological formations, the results of soil radon emanometry show the possibility of radon gas accumulation inside the undermined spaces, but their permanent subsidence causes loosing of overlying material and escape of radon gas. The boundaries (edges) of undermined and sunken areas were identified as the only places with increased values of 222Rn activity, probably due to the presence of vertical supporting mine walls allowing radon gas accumulation and upward movement. Thus, the soil radon emanometry clearly indicates the borders between undermined or sunken and non-undermined parts.

BeiDou Navigation Satellite System (BDS) is composed of satellites in geostationary Earth orbit (GEO), medium Earth orbit (MEO) and inclined geosynchronous orbit (IGSO). However, the orbit determination of geostationary Earth orbits and of geosynchronous orbits (GSO) with small inclination angle and small eccentricity is a challenging task that is addressed in this paper using Extended Kalman Filter (EKF). The satellite positions were predicted in Earth-centred inertial (ECI) reference frame when propagated through Keplerian model and perturbation force model for different values of right ascension of ascending node (RAAN). Root mean square (RMS) errors of 9.61 cm, 6.73 cm and 11.46 cm were observed in ECI X, Y and Z satellite position coordinates of GSO respectively, whereas, the RMS errors for GEO satellite were 8.89 cm, 7.92 cm, and 0.93 cm respectively in ECI X, Y and Z coordinates; for perturbation force model with maximum value of RAAN when compared with dynamic orbit determination model. Kolmogorov-Smirnov test for EKF reported a p-value > 0.05, indicating a good fit of perturbation force model for orbit propagation. Orbit determination using EKF with perturbation force model were compared with that using EKF with Kepler's model. Wilcoxon Rank Sum test was used to compare the residuals from EKF algorithm through Kepler's model and perturbation force model. EKF with Perturbation force model showed improvement in predicting the satellite positions as compared to Kepler's model. EKF with Perturbation force model was further applied to International GNSS Service (IGS) station data and kilometre level accuracy was achieved. RMS errors of 0.75 km, 2.53 km and 1.91 km were observed in ECI X, Y and Z satellite position coordinates of GSO, respectively, whereas, the RMS errors for GEO satellite were 3.89 km, 4.20 km and 6.66 km respectively in ECI X, Y and Z coordinates for perturbation force model.

3D post-stack time migrated seismic data and a suite of composite well log data from six wells drilled within the “AFUN” field Niger delta were used to effect a detailed interpretation of the field. This was with a view to delineating architectural elements that control reservoir quality of a deepwater turbidite reservoir. The data analyses were done using the Petrel software. LAS file of logs were imported into the Petrel software as well as SEG.Y. seismic data. Fault interpretation and horizon mapping were based on the well-seismic tie from the generated seismogram. Time and depth structure maps were created. Thirty faults which include growth faults, reverse faults, collapsed crest structure and as well as faults that are synthetic and antithetic to the growth faults were mapped. The growth faults are believed to act as pathways for the updip movement of hydrocarbon from the Akata Formation to Agbada Formation. The structural interpretation showed that the area has been subjected to compressional deformation which resulted in reverse faulting system in toe thrust zone influenced by shale diapirs. The maps revealed contour closures that belong to an anticlinal structure which is forming traps in the reservoirs. The structures are faulted North-South trending rollover anticlines. It has also been shown that the distribution and type of architectural elements i.e. fractures within the fan system have major impact upon the reservoir distribution, continuity and connectivity of sand/shale bodies. The study concluded that structural style and facies architecture are the two fundamental elements that defined the reservoir heterogeneity of the “AFUN” Field.

Interpretation and inversion of microgravity anomalies belong to important tasks of near-surface geophysics, mostly in cavities detection in engineering, environmental and archaeological applications. One of the mostly used concepts of inversion in applied gravimetry is based on the approximation of the model space by means of 2D or 3D elementary sources with the aim to estimate their densities by means of the solution of a corresponding linear equation system. There were published several approaches trying to obtain correct and realistic results, which describe real parameters of the sources. In the proposed contribution we analyse the properties of two additional functionals, which describe additional properties of the searched solution – namely so-called L2-smoothing and minimum support focusing stabilizers. For the inversion itself, we have used the regularized conjugate gradient method. We have studied properties of these two stabilizers in the case of one synthetic model and one real-world dataset (microgravity data from St. Nicholas church in Trnava). Results have shown that proposed algorithm with the minimum support stabilizer can generate satisfactory model results, from which we can describe real geometry, dimensions and physical properties of interpreted cavities.

Study presented in this paper is focused on comparison and statistical assessment of differences between the selected Level 2 products of the satellite mission Gravity Recovery and Climate Experiment (GRACE). Global monthly gravity field models in terms of spherical harmonic coefficients produced by three institutes of GRACE Science Data System are compared with the partially independent MASCON global gravity field model. Detailed comparison and statistical analysis of differences is performed in 5 selected river basins: Amazon, Congo, Danube, Yenisei and Lena. For each spherical harmonic solution, 8 different filtrations available at International Center for Global Gravity Field Models (ICGEM) are tested over the time span from April 2002 to July 2016. Fischer test at two significance levels 10% and 5% has been performed in order to qualify the statistical significance between the particular solutions.

Santorini is located in the central part of the Hellenic Volcanic Arc (South Aegean Sea) and is well known for the Late-Bronze-Age “Minoan” eruption that may have been responsible for the decline of the great Minoan civilization on the island of Crete. To use gravity to probe the internal structure of the volcano and to determine whether there are temporal variations in gravity due to near surface changes, we construct two gravity maps. Dionysos Satellite Observatory (DSO) of the National Technical University of Athens (NTUA) carried out terrestrial gravity measurements in December 2012 and in September 2014 at selected locations on Thera, Nea Kameni, Palea Kameni, Therasia, Aspronisi and Christiana islands. Absolute gravity values were calculated using raw gravity data at every station for all datasets. The results were compared with gravity measurements performed in July 1976 by DSO/NTUA and absolute gravity values derived from the Hellenic Military Geographical Service (HMGS) and other sources. Marine gravity data that were collected during the PROTEUS project in November and December 2015 fill between the land gravity datasets. An appropriate Digital Elevation Model (DEM) with topographic and bathymetric data was also produced. Finally, based on the two combined datasets (one for 2012–2014 and one for the 1970s), Free air and complete Bouguer gravity anomaly maps were produced following the appropriate data corrections and reductions. The pattern of complete Bouguer gravity anomaly maps was consistent with seismological results within the caldera. Finally from the comparison of the measurements made at the same place, we found that, within the caldera, the inner process of the volcano is ongoing both before, and after, the unrest period of 2011–2012.

A detailed knowledge of the thickness of crust and upper mantle structure is important for understanding a plate tectonics and geodynamics in the region. We use body wave for detecting details of the subsurface structure. The information in this research is collected from a seismic linear profile that extends across the Sanandaj-Sirjan metamorphic zone in seismic states of Central Iran and Zagros. We compute P receiver functions to investigate crustal and upper mantle discontinuities. We use teleseismic events (mb ≥ 5.5, 30° < Δ < 95°) registered between 1996 and 2018 and recorded at 10 short-period stations with 3 components and 17 broadband stations with high signal to noise ratio. The observed depth of Moho in the study area is approximately 50 km and rises to 70 km at the end of the seismic linear profile beneath Sanandaj-Sirjan zone. In Central Iran, depths discontinuities in the transition zone are shown by the reference model of deviation, which can be attributed to the convergence of Arabian plate with the Central Iran plateau. Also, the study area was identified as geothermal susceptibility by SUNA and this observation was confirmed.

A quasi 3D electrical resistivity (ERT) survey was undertaken at a UNESCO World Heritage site, Saqqara, Giza, Egypt, during a joint archaeological-geophysical mission from Cairo University. The main objective is to detect the locations of the subsurface archaeological tombs/or crypts and to allocate any possible archaeological bodies/features buried underneath the study area. In this survey, SYSCAL Pro system with 24 electrodes and a multi-core cable is used for automatic data acquisition of profiling data. The dipole–dipole array was used to enhance resolution, 14 resistivity lines are conducted during this Survey. The processed data were analysed in order to produce resistivity tomography (ERT) for qualitative and quantitative interpretations. Inversion of the ERT data identified variation of resistivity values and the expected locations of the underground galleries and highlight the presence of regular shape structures probably due to features of archaeological interest. Excavations made accordingly in the study area led to an interesting discovery of a tomb of the Great Army General, Iwrhya. The tomb is approximately 2000 years old as it covers the reigns of both Kings Seth I and Ramesses II. Using the 3D resistivity tomography with such a multi-electrode technique proved its efficiency and applicability for non-invasive archaeo-geophysical prospecting.