Unraveling Hydrocarbon Microseepages in Onshore Basins Using Spectral–Spatial Processing of Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER) Data

This study focuses on the spectral characterization of superficial effects caused by natural gaseous hydrocarbon (HC) escapes (microseepages) on onshore basins and means to detect the phenomena using Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER) imagery. The Northern Tucano basin (NTB) (Brazil) displays archetypal onshore microseepages, and it is used here as a case study area. Regional HC geochemical data yielded from NTB soil samples helped to outline the strategies for imagery processing and assessment of results. We firstly devised a generic detection model for onshore basins in which microseepages indicated by HC geochemistry should spatially match a number of surface expressions with key remote sensing signatures, such as bleached soil and rocks (i.e., reduction in Fe³⁺ to Fe²⁺), concentration of specific clays (kaolinite), and carbonates (siderite) and geobotanical markers. These evidences guided the use of remote sensing data and information extraction techniques to find more sites favorable to host HC microseepages in the NTB. ASTER data processed through spectral mixing analysis and supervised and unsupervised neural network systems highlighted sites where features predicted in the detection model concurred seamlessly with geochemical anomalies. The approach also revealed other sites with similar characteristics, but for which no geochemical data exist. These sites are taken as new potential targets for the presence of microseepages and possibly HC reservoirs. The research demonstrated the potential of remote sensing data and spectral–spatial methodologies for economical, onshore detection of HCs for exploration purposes. Potential ramifications of the study include indirect detection of gaseous HC reservoirs in terrestrial planets (fostering the search for signs of life) and depiction of non-negligible, replenishable continental sources of HC emissions on Earth that may contribute to global warming.     

Shallow Crustal Configuration of the Narmada–Son Lineament Transition Zone Near the Sahdol–Katni Area of Central India Using Simultaneous Gravity and Magnetic Observations

Gravity and magnetic studies have been carried out to map the different depth formations of alluvium, Gondwana, Vindhyan, Mahakoshal, and the crystalline basement in the Narmada–Son lineament (NSL) near the Sahdol–Katni area, India. Higher elevations in the northern part of the study area have lower gravity; the southern part of the study area, however, is moderately elevated and also has a higher gravity anomaly, which justifies the isostatic adjustment. This indicates the presence of high-density material, for example Mahakoshal rocks, in the upper crust which causes the higher anomaly value. The Mahakoshal rocks are widely exposed at Mau, Chanaura, Nadawar, Khamaria, Umria, and near the Tala–Barhi area and also extend from Sidhi to Agoni village further east of the area. The Mahakoshal rocks are thinner between the river Son and Tikwa village which causes a different gravity anomaly pattern changing from the NNE–SSW direction to the N–S direction. However, the trend of magnetic anomaly follows the same pattern toward the NNE–SSW direction, mainly because of the effect of the crystalline basement. Two magnetic highs are prominent in the Tikwa and Amarpur regions, 800 and 400 nTesla, respectively, because of the presence of the crystalline basement. The depth of the crystalline basement obtained by 2½ D gravity–magnetic modelling varies from 2.7 to 2.9 km. From spectral analysis the average depth of the crystalline basement varies from 2.83 to 3.05 km. The different crustal depth sections obtained by 2½ D simultaneous gravity–magnetic modelling correlate well with other constrained data.     

Spatial and Spectral Analysis of Refined Gravity Data for Modelling the Crust–Mantle Interface and Mantle-Lithosphere Structure

We analyse spatial and spectral characteristics of various refined gravity data used for modelling and gravimetric interpretation of the crust–mantle interface and the mantle-lithosphere structure. Depending on the purpose of the study, refined gravity data have either a strong or weak correlation with the Moho depths (Moho geometry). The compilation of the refined gravity data is purely based on available information on the crustal density structure obtained from seismic surveys without adopting any isostatic hypothesis. We demonstrate that the crust-stripped relative-to-mantle gravity data have a weak correlation with the CRUST2.0 Moho depths of about 0.02. Since gravitational signals due to the crustal density structure and the Moho geometry are subtracted from gravity field, these refined gravity data comprise mainly the information on the mantle lithosphere and sub-lithospheric mantle. On the other hand, the consolidated crust-stripped gravity data, obtained from the gravity field after applying the crust density contrast stripping corrections, comprise mainly the gravitational signal of the Moho geometry, although they also contain the gravitational signal due to anomalous mass density structures within the mantle. In the absence of global models of the mantle structure, the best possible option of computing refined gravity data, suitable for the recovery/refinement of the Moho interface, is to subtract the complete crust-corrected gravity data from the consolidated crust-stripped gravity data. These refined gravity data, that is, the homogenous crust gravity data, have a strong absolute correlation of about 0.99 with the CRUST2.0 Moho depths due to removing a gravitational signal of inhomogeneous density structures within the crust and mantle. Results of the spectral signal decomposition and the subsequent correlation analysis reveal that the correlation of the homogenous crust gravity data with the Moho depths is larger than 0.9 over the investigated harmonic spectrum up to harmonic degree 90. The crust-stripped relative-to-mantle gravity data correlate substantially with the Moho depths above harmonic degree 50 where the correlation exceeds 0.5.     

Rapid source inversions of the 2023 SE Türkiye earthquakes with teleseismic and strong-motion data

We conducted rapid inversions of rupture process for the 2023 earthquake doublet occurred in SE Türkiye, the first with a magnitude of M_W7.8 and the second with a magnitude of M_W7.6, using teleseismic and strong-motion data. The teleseismic rupture models of the both events were obtained approximately 88 and 55 minutes after their occurrences, respectively. The rupture models indicated that the first event was an asymmetric bilateral event with ruptures mainly propagating t...     

Exploration and Study of Deep Crustal Structure in the Quanzhou Basin and Its Adjacent Area

INTRODUCTIONQuanzhou Citylocates on the southeast coast of Fujian Province .Tectonically,it locates in thevolcanic fault depression zone in East Fujian betweenthe Wuyi-Daiyunfolded doming-upregion of theSouth China block and the depression zone of the Taiwan Straits . Being on the north segment of thesoutheast seismic belt of China , many destructive earthquakes inthe history affected the area and theregional seismicityis closelyrelated withfaulting.Inthe Quanzhou basin,large-scalelo…     

The Role of Gravity, Magnetic and Electric Prospecting Techniques in the Geological Survey of the Qiangtang Basin

The Qiangtang basin is one of the large-scale basins in China, which is located in the high-altitude and cold area. In 1994 the preliminary gravity, magnetic and electric prospecting work was undertaken in this area. Over three years the gravity, magnetic and electric reconnaissance work had been completed. The integrated interpretation of complex gravity, magnetic and electric data had identified the boundary of the basin. The relationship between the basin and mountains is a south-north ramp structure. The basin structure consisted of one uplift and two sags, five depressions and two arches. The complex evaluation showed that the regions of Baitanhu and Youyiquan in the Qiangtang basin were considered to be the most perspective ones in sense of oil-gas bearing.     

Analysis of faulting destruction and water supply pipeline damage from the first mainshock of the February 6, 2023 Türkiye earthquake doublet

In 2023, two consecutive earthquakes exceeding a magnitude of 7 occurred in Türkiye, causing severe casualties and economic losses. The damage to critical urban infrastructure and building structures, including highways, railroads, and water supply pipelines, was particularly severe in areas where these structures intersected the seismogenic fault. Critical infrastructure projects that traverse active faults are susceptible to the influence of fault movement, pulse velocity, and ground motions. ...     

Edge Detection and Depth Estimation Using a Tilt Angle Map from Gravity Gradient Data of the Kozaklı-Central Anatolia Region,Turkey

In this paper the application of an edge detection technique to gravity data is described. The technique is based on the tilt angle map (TAM) obtained from the first vertical gradient of a gravity anomaly. The zero contours of the tilt angle correspond to the boundaries of geologic discontinuities and are used to detect the linear features in gravity data. I also present that the distance between zero and ±p\mathord

Velocity Structure of the West-Bengal Sedimentary Basin, India along the Palashi-Kandi Profile Using a Travel-time Inversion of Wide-angle Seismic data and Gravity Modeling – An Update

The 2-D shallow velocity structure along the north-south Palashi-Kandi profile in the West Bengal sedimentary basin has been updated by travel-time inversion of seismic refraction, wide-angle reflection and gravity data. A six-layer shallow model up to a depth of about 7 km has been derived. The first layer, which has an average velocity of 2.0 kms^?1, represents the alluvium deposit, which rests over the shale formation with average velocity of 3.0 kms^?1. The thin (200 m) Sylhet limestone, observed at a nearby Palashi well, remains hidden in the present data set. Hence a 200-m thin layer with a velocity of 3.7 kms^?1, corresponding to the Sylhet limestone, has been assumed to be present throughout the profile. The fourth layer with a velocity of 4.5–4.7 kms^?1 at a depth of 1.7–2.4 km represents the Rajmahal traps. The ‘skip’ phenomenon and rapid amplitude decay of first arrivals indicate a low-velocity zone (LVZ) in the study area. Using the ‘skip’ phenomena and wide-angle reflection data, identified on seismograms, the LVZ with a velocity of 4.0 kms^?1, indicating the Gondwana sediments, has been delineated below the Rajmahal traps. The next layer with a velocity 5.4–5.6 kms^?1 overlying the crystalline basement (5.8–6.25 kms^?1) may be associated with the Singhbhum group of meta volcanic rock that has been exposed in the western part of the basin. The basement lies at a variable depth of 4.9 to 6.8 km. The overall uncertainties of various velocity and boundary nodes are ± 0.15 kms^?1 and ± 0.5 km, respectively. The elevated basement feature in the north might have acted as a structural barrier for the deposition of Sylhet limestone during the Eocene epoch. The seismically derived shallow structure correctly explains the observed Bouguer gravity anomaly along the profile. The addition of reflections in the present analysis provides a stronger control on the depths and velocities of basement and overlying sedimentary formations, compared to the earlier model derived mainly by the first arrival seismic data.     

Imaging the Seismic Crustal Structure of the Western Mexican Margin between 19°N and 21°N

Three thousand kilometres of multichannel (MCS) and wide-angle seismic profiles, gravity and magnetic, multibeam bathymetry and backscatter data were recorded in the offshore area of the west coast of Mexico and the Gulf of California during the spring 1996 (CORTES survey). The seismic images obtained off Puerto Vallarta, Mexico, in the Jalisco subduction zone extend from the oceanic domain up to the continental shelf, and significantly improve the knowledge of the internal crustal structure of the subduction zone between the Rivera and North American (NA) Plates. Analyzing the crustal images, we differentiate: (1) An oceanic domain with an important variation in sediment thickness ranging from 2.5 to 1 km southwards; (2) an accretionary prism comprised of highly deformed sediments, extending for a maximum width of 15 km; (3) a deformed forearc basin domain which is 25 km wide in the northern section, and is not seen towards the south where the continental slope connects directly with the accretionary prism and trench, thus suggesting a different deformational process; and (4) a continental domain consisting of a continental slope and a mid slope terrace, with a bottom simulating reflector (BSR) identified in the first second of the MCS profiles. The existence of a developed accretionary prism suggests a subduction–accretion type tectonic regime. Detailed analysis of the seismic reflection data in the oceanic domain reveals high amplitude reflections at around 6 s [two way travel time (twtt)] that clearly define the subduction plane. At 2 s (twtt) depth we identify a strong reflection which we interpret as the Moho discontinuity. We have measured a mean dip angle of 7° ± 1° at the subduction zone where the Rivera Plate begins to subduct, with the dip angle gently increasing towards the south. The oceanic crust has a mean crustal thickness of 6.0–6.5 km. We also find evidence indicating that the Rivera Plate possibly subducts at very low angles beneath the Tres Marias Islands.     

Examination of Historical Data of the 1125 Lanzhou M7.0 Earthquake and Analysis of Seismogenic Structure

Detailed examination of historical data of earthquakes and field investigations of loess landslide caused by the earthquake and tracing of active faults in Lanzhou area indicate that the Yijitanpu town,one of six towns of Jincheng city,was devastated by the 1125 Lanzhou earthquake.The citly is now located in the Vinylon Factory south of Hekou(River Mouth)in the Xigu distict of Lanzhou city.We delermined that the six old towns mentioned in historical records lie in an area stretching from the south of Xigu district to Hekou in Lanzhou.This is consistent with the distibution of loess landslides caused by the earthquake,the extension of Holocene active faults,and the distribution of traces of the seismic rupture zone.A comprehensive analysis shows that the seismogenic structure for the 1125 Lanzhou M7.0 earthquake should be the Xianshuigou fault segment at the western termination of the north-border active fault zone of the Maxianshan Mountains which are located in south of Lanzhou city with the distance of only 4km.     

Structure and Significance of S-wave Velocity and Poisson’s Ratio in the Crust beneath the Eastern Side of the Qinghai-Tibet Plateau

The receiver functions of body waves of distant earthquakes obtained for the regions beneath 41 digital stations (Lhasa and GANZ in Tibet, Mandalay and Rangoon in Myanmar, SHIO in India, CHTO in Thailand, and station network in Sichuan and Yunnan) were used to invert for S-wave structure in the crust and upper mantle in Sichuan, Yunnan, and their surrounding areas. Meanwhile the distribution characteristics of the Poisson’s ration and the crustal thickness in Sichuan and Yunnan areas were also obtained. Results indicate that the depth of Moho beneath the eastern side of Qinghai-Tibetan plateau varies strikingly. It is obvious that the greatest changes in crustal thickness occur in a north-south direction. The crustal thickness decreases from north to south, being as thick as 70 km in eastern Tibet, the northern portion of our area of interest, and less than 30 km in Chaing Mai and Rangoon, the southern portion of our area. There are, however, exceptions regarding the trend. The thickness exhibits an east-west variation trend in the area from Ma’erkong-Kongding in Sichuan to Lijiang in Yunnan. In general the Jinpingshan-Longmenshan fault and Anninghe fault can be taken as the boundaries of this exception area. The thickness in Kongding in the west is 68 km, while it is only 39 km in Yongchuan in the east. Moreover the Poisson’s ratio values in the blocks of central Sichuan and Sichuan-Yunnan Diamond are high, and a low velocity layer in the crust of this area can be obviously detected. The distribution characteristics of the high Poisson’s ratio and the low velocity of the crust in this block correspond to the tectonic structure, being in contrast with the surrounding areas. Combining with the distribution features of the modern tectonic stress field, it is deduced that the Sichuan-Yunnan area is probably the channel through which the materials of the lithosphere flow eastward.     

Determination of Q β(f) in Different Parts of Kumaon Himalaya from the Inversion of Spectral Acceleration Data

This paper presents the results of a modified two-step inversion algorithm approach to find S wave quality factor Q _β(f) given by Joshi (Bull Seis Soc Am 96:2165–2180, 2006). Seismic moment is calculated from the source displacement spectra of the S wave using both horizontal components. Average value of seismic moment computed from two horizontal components recorded at several stations is used as an input to the first part of inversion together with the spectra of S phase in the acceleration record. Several values of the corner frequency have been selected iteratively and are used as inputs to the inversion algorithm. Solution corresponding to minimum root mean square error (RMSE) is used for obtaining the final estimate of Q _β(f) relation. The estimates of seismic moment, corner frequency and Q _β(f) from the first part of inversion are further used for obtaining the residual of theoretical and observed source spectra which are treated as site amplification terms. The acceleration record corrected for the site amplification term is used for determination of seismic moment from source spectra by using Q _β(f) obtained from first part of inversion. Corrected acceleration record and new estimate of seismic moment are used as inputs to the second part of the inversion scheme which is similar to the first part except for use of input data. The final outcome from this part of inversion is a new Q _β(f) relation together with known values of seismic moment and corner frequency of each input. The process of two-step inversion is repeated for this new estimate of seismic moment and goes on until minimum RMSE is obtained which gives final estimate of Q _β(f) at each station and corner frequency of input events. The Pithoragarh district in the state of Uttarakhand in India lies in the border region of India and Nepal and is part of the seismically active Kumaon Himalaya zone. A network of eight strong motion recorders has been installed in this region since March, 2006. In this study we have analyzed data from 18 local events recorded between March, 2006 and October, 2010 at various stations. These events have been located using HYPO71 and data has been used to obtain frequency-dependent shear-wave attenuation. The Q _β(f) at each station is calculated by using both the north-south (NS) and east-west (EW) components of acceleration records as inputs to the developed inversion algorithm. The average Q _β(f) values obtained from Q _β(f) values at different stations from both NS and EW components have been used to compute a regional average relationship for the Pithoragarh region of Kumaon Himalaya of form Q _β(f)?=?(29?±?1.2)f ^{(1.1 ± 0.06)}.     

Crustal Movements and Deformations in Eastern Turkey in Connection with the Van Earthquake (October 23, 2011, Мw = 7.2): Study from GPS Data

Izvestiya, Physics of the Solid Earth - Abstract—A devastating earthquake with magnitude Mw = 7.2 occurred in the Van region, Eastern Turkey, on October 23, 2011. This study analyzes the...     

The Rotational Structure along the Eastern Segment of the Altyn Tagh Fault and Its Implications in Dynamics

As a result of the left-lateral strike-slipping of the Altyn Tagh fault in Neotectonic period, a contra-rotational structure, namely the Zhaobishan vortex structure, has developed at the juncture of the main Altyn Tagh fault and the northern fringe fault of the Qilian Mountains.Preliminary analysis on the deformation and evolution of the Zhaobishan vortex structure. In combination with the previous data, suggests that the tectonic transform between the Altyn Tagh fault and the northern fringe fault of the Qilian Mountains attributes to the deformation of the rotational structure. The existence of a series of rotational structures along the Altyn Tagh fault and on the northeastern edge of the Qinghai-Xizang(Tibet) plateau indicate that as the substance in the northern Qinghal-Xizang (Tibet) plateau moves clockwise around the eastern tectonic knot of the Himalayas, rotational structures become the principal mode on the northern marginal zone of the Plateau of transforming and absorbing tectonic deformation.     

Gravity Data from the Teboursouk Area (“Diapirs Zone”, Northern Tunisia): Characterization of Deep Structures and Updated Tectonic Pattern

Located between eastern segments of the Atlas and Tell-Rif oro-genic belts, the “Dome zone” of northern Tunisia is characterized by the juxtaposition of various structures that mainly controlled the long geody-namic history of this part of the south-Tethyan Margin. To better understand the organization and deep extension of these structures, gravity data from the Teboursouk key area are proposed. These data include the plotting of Bouguer anomaly map and related parameters such as vertical and horizontal gradients, upward continuation and Euler solution. Compared to geological and structural maps available, they allow the identification of new deep structures and greater precision regarding the characteristics and organization of known ones; consequently, an updated structural pattern is proposed.     

Comparison of Factorial Kriging Analysis Method and Upward Continuation Filter to Recognize Subsurface Structures — A Case Study: Gravity Data from a Hydrocarbon Field in the Southeast Sedimentary Basins of the East Vietnam Sea

To interpret geophysical anomaly maps, it is necessary to filter out regional and sometimes noise components. Each measured value in a gravity survey consists of different components. Upward continuation (UC) is one of the most widely used filters. The shortcoming of this filter is not to consider the spatial structure of the data, and also the fact that the trial and error approach and expert’s judgment are needed to adjust it. This study aims to compare the factorial kriging analysis (FKA) and UC filters for separation of local and regional anomalies in the gravity data of a hydrocarbon field in the southeast sedimentary basins of the East Vietnam Sea. As shown in this paper, FKA method permits to filter out all of the identified structures, while the UC filter does not possess this capability. Therefore, beside general and classic filtering methods, the FKA method can be used as a strong method in filtering spatial structures and anomaly component.     

Pulse-like ground motion observed during the 6 February 2023 MW7.8 Pazarcık Earthquake (Kahramanmaraş, SE Türkiye)

In this study, we analyzed 100 three-component strong ground motion records observed within 200 km of the causative fault of the 6 February 2023 M_W7.8 Pazarcık (Kahramanmaraş) Earthquake in SE Türkiye. The wavelet method was utilized to identify and analyze the characteristics of pulse-like ground motions in the near-fault region, while considering the uncertainty of the pulse orientation during the analysis. Our investigation focused on the effects of the focal mechanism and rupture process on the spatial distribution, pulse orientation, and maximum pulse direction of the observed pulse-like ground motion. We also analyzed the amplitude and period of the observed ground pulses and the effect of long-period amplification on the ground motion response spectra. Our results indicated the following: (1) A total of 21 typical ground velocity pulses were observed during this earthquake, exhibiting complex characteristics due to the influence of the strike-slip mechanism and rupture directivity. Most ground pulses (17 out of 21) were recorded within 20 km of the fault, in a wide range of orientations, including normal and parallel to the fault direction. The waveforms exhibited unidirectional features, indicating the effects of left-lateral fault slip. Distinct pulses observed more than 20 km from the fault were mainly oriented normal to the fault. The waveforms were bidirectional with double- or multi-round trips as a result of rupture directivity. (2) The amplitudes of the observed pulses ranged from 30.5 to 220.0 cm/s, with the largest peak velocity of 220.0 cm/s observed at Station 3138. The pulse periods ranged from 2.3 to 14.5 s, with the longest pulse period of 14.5 s observed at Station 3116. The amplitude and period of the pulses observed during this earthquake were comparable to those of similar-magnitude global earthquakes. The amplitude of the pulses decreased significantly with increasing fault distance, whereas the pulse period was not significantly affected by the fault distance. (3) Compared with non-pulse records, the velocity pulse records had a pronounced amplification effect on the acceleration response spectra near the pulse period, with factors ranging from 2.1 to 5.8. The larger velocity pulses also significantly amplified the velocity response spectra, particularly over the long periods. This significant amplification effect of the pulses on the response spectra leads to empirical models underestimating the long-period earthquake ground motion.