Lateral Variations Research Articles

Machine learning in prenatal MRI predicts postnatal ventricular abnormalities in fetuses with isolated ventriculomegaly.

To evaluate the intracranial structures and brain parenchyma radiomics surrounding the occipital horn of the lateral ventricle in normal fetuses (NFs) and fetuses with ventriculomegaly (FVs), as well as to predict postnatally enlarged lateral ventricle alterations in FVs. Between January 2014 and August 2023, 141 NFs and 101 FVs underwent 1.5 T balanced steady-state free precession (BSSFP), including 68 FVs with resolved lateral ventricles (FVM-resolved) and 33 FVs with stable lateral ventricles (FVM-stable). Demographic data and intracranial structures were analyzed. To predict the enlarged ventricle alterations of FVs postnatally, logistic regression models with 5-fold cross-validation were developed based on lateral ventricle morphology, blended-cortical or/and subcortical radiomics characteristics. Validation of the models' performance was conducted using the receiver operating characteristic curve (ROC), calibration curve, and decision curve analysis (DCA). Significant alterations in cerebral structures were observed between NFs and FVs (p < 0.05), excluding the maximum frontal horn diameter (FD). However, there was no notable distinction between the FVM-resolved and FVM-stable groups (all p > 0.05). Based on subcortical-radiomics on the aberrant sides of FVs, this approach exhibited high efficacy in distinguishing NFs from FVs in the training/validation set, yielding an impressive AUC of 1/0.992. With an AUC value of 0.822/0.743 in the training/validation set, the Subcortical-radiomics model demonstrated its ability to predict lateral ventricle alterations in FVs, which had the greatest predictive advantages indicated by DCA. Microstructural alterations in subcortical parenchyma associated with ventriculomegaly can serve as predictive indicators for postnatal lateral ventricle variations in FVs. It is critical to gain pertinent information from a solitary fetal MRI to anticipate postnatal lateral ventricle alterations in fetuses with ventriculomegaly. This approach holds the potential to diminish the necessity for recurrent prenatal ultrasound or MRI examinations. Fetal ventriculomegaly is a dynamic condition that affects postnatal neurodevelopment. Machine learning and subcortical-radiomics can predict postnatal alterations in the lateral ventricle. Machine learning, applied to single-fetal MRI, might reduce required antenatal testing.

Seismic-aseismic impact on marine depositional dynamics and sedimentary architecture: Study on multilayer-multiscale SSDS in Proterozoic Rohtas Limestone, India

The present study explores a diverse suite of soft-sediment deformation structures (SSDSs) preserved at multiple levels within the Mesoproterozoic carbonate succession of Rohtas Limestone in the Vindhyan Supergroup of India. The carbonate succession shows lithological heterogeneity, characterized by a sharp change of facies association and depositional environment. Three layers of SSDS, traced over 57 km, demarcate the lower lagoonal from upper open shelf deposits. Contorted laminae, load casts and flame structures, convolute laminae, ball and pillow structures, accordion folds and pillar structures within the limestone beds infer liquefaction-fluidization and plastic deformation. The lateral continuity, vertical repetition and confinement of SSDS beds at the interface of lagoon and open shelf deposits indicate the role of seismicity on facies succession. The SSDS beds show lateral variation in thickness and structures, depending on the distance from epicenter. The seismic trigger was accompanied by sudden deepening of the depositional environment and related marine transgression. The coarsening up packages within the lagoonal succession comprise load cast and flame structures, syn-sedimentary faults and intraclastic conglomerates. The open shelf succession shows several fining up cycles, with multiple load and flame structures with intraclastic conglomerate at the top. Origin of SSDS in association of the coarsening up sequence, and fining up cycle infers that intrabasinal tectonics and storms were the triggering agents for lagoon and open shelf successions respectively. The seismite in the studied succession not only reveals the seismicity during the Proterozoic, but also changed the sedimentary facies association, which may use as a diagnostic criterion for the identification of seismite. This study also highlights that seismic and aseismic impacts on sedimentary architecture and paleoceanography evolution.

A Spectral Method to Compute the Tides of Laterally Heterogeneous Bodies

Body tides reveal information about planetary interiors and affect their evolution. Most models to compute body tides rely on the assumption of a spherically symmetric interior. However, several processes can lead to lateral variations of interior properties. We present a new spectral method to compute the tidal response of laterally heterogeneous bodies. Compared to previous spectral methods, our approach is not limited to small-amplitude lateral variations; compared to finite element codes, this approach is more computationally efficient. While the tidal response of a spherically symmetric body has the same wavelength as the tidal force; lateral heterogeneities produce an additional tidal response with a spectra that depends on the spatial pattern of such variations. For Mercury, the Moon, and Io, the amplitude of this signal is as high as 1%–10% of the main tidal response for long-wavelength shear modulus variations higher than ∼10% of the mean shear modulus. For Europa, Ganymede, and Enceladus, shell-thickness variations of 50% of the mean shell thickness can cause an additional signal of ∼1% and ∼10% for the Jovian moons and Encelaudus, respectively. Future missions, such as BepiColombo and JUICE, might measure these signals. Lateral variations of viscosity affect the distribution of tidal heating. This can drive the thermal evolution of tidally active bodies and affect the distribution of active regions.

Tomographic imaging of South Rewa basin, Central India: Implications of Gondwana rifting and Late Cretaceous volcanism

We have imaged thick (3.0–5.0 km) sub-trappean Gondwana sediments deposited in a typical rift-environment of the south Rewa basin in Central India masked by the Deccan basalts. The basalts intrude through the basin forming several dykes, sills, and other intrusions along with large-scale undulations of the basement showing complex geology and tectonic settings due to Late Cretaceous volcanism (∼65 Ma). To delineate complex subsurface geological structures and understand the tectonic settings, we use long-offset seismic reflection data along the N-S trending 127 km Kuvari-Shahdol profile of the basin. The conventional stack and pre-stack time migration techniques fail to image the complex subsurface structures dominated by steeply dipping faults. The steeply dipping faults with small-scale structures are influenced by strong lateral and vertical velocity variations that cause problems in imaging the basin. Hence, a robust tomographic inversion with pre-stack depth migration technique is adopted to image fine-scale subtle subsurface geological structures with smooth velocity variations. The tomographic velocity model and pre-stack depth migration seismic image could decipher basaltic trap thickness, numerous dyke intrusions, alternate horsts and grabens with deposition of thick hydrocarbon-bearing sub-trappean Gondwana rocks along with the basement configuration. The deep basinal faults with highly distorted basement structures represent the intense tectonic activity of this Gondwana rift-basin. The deep-seated hydrocarbon reservoir is discovered with an excellent trapping mechanism forming both pre-rift and syn-rift settings of the rift-basin. The seismic image is further constrained and corroborated by the inversion of residual Bouguer gravity anomaly data to obtain corresponding density model derived from tomographic velocity model. The large basement faults, horsts and grabens, several intrusives like dyke, sill, laccolith, and thick sub-trappean Gondwana formations with Deccan trap cover indicate complex geology and tectonic settings of south Rewa basin forming the Late Cretaceous Volcanic Province of India.

A Common Diffusional Mechanism for Creep and Grain Growth in Polymineralic Rocks: Application to Lower Mantle Viscosity Estimates

AbstractIn a previous study (Okamoto &amp; Hiraga, 2022, https://doi.org/10.1029/2022jb024638), we concluded that diffusion creep and grain growth in polymineralic rocks proceed by a common diffusional mechanism. Here, we built on that finding and estimated lower mantle grain size and viscosity during a single mantle convection cycle dominated by diffusion creep. We approximated the lower mantle as a two‐phase material consisting of bridgmanite + ferropericlase and post‐perovskite + ferropericlase, depending on depth. We used previously reported self‐diffusivities for bridgmanite and post‐perovskite. We predict a bridgmanite grain‐size of tens to hundreds of microns shortly after the phase transition at ∼660 km depth. This size remains relatively constant until the mantle material enters the post‐perovskite zone, which is marked by significant grain growth up to ∼9 mm just prior to upwelling. This size is sufficient to prevent further grain growth until the mantle material reaches the top of the lower mantle. These grain sizes combined with the diffusivities yield viscosities that vary laterally and with depth. At a lateral temperature difference of up to 800 K in the lower mantle, fine‐grained cold downwelling mantle is almost as viscous as, or more likely to be softer than coarse‐grained hot upwelling mantle. The lateral viscosity variations cannot be more than 2 orders of magnitude, and we estimated viscosities of 1018–1020 Pa · s in the upper lower mantle, 5 × 1020–5 × 1022 Pa · s in the lower bridgmanite zone, and 1017–1019 Pa · s in the post‐perovskite zone, which compare well with the values estimated in previous geophysical modeling studies.

Integrated analysis of wireline logs analysis, seismic interpretation, and machine learning for reservoir characterisation: Insights from the late Eocene McKee Formation, onshore Taranaki Basin, New Zealand

This study aims to integrate seismic and well data for re-assessment and delineating the geological and petrophysical characteristics of the Late Eocene McKee Formation. Despite the McKee Formation being a hydrocarbon-producing field, there have been no previous integrated studies conducted to characterize its deposits in the Taranaki Basin. To execute these tasks, the process involves the interpretation of seismic data, wireline analysis, and 3D porosity volume prediction using multi-layered feed-forward neural network models. The seismic data was interpreted and found to contain low displacement faults due to intricate geometric complexities, the contour map that was generated shows lateral extension and thickness variation, and the region of direct hydrocarbon indicators was detected using the root-mean-square. Additionally, the seismic inversion provides a comprehensive understanding of the acoustic impedance variations from low to high value, which signifies the presence of mudstone clasts content. Furthermore, in the wireline log analysis, the reservoir was found to be a Silici-clastic and carbonate field owing to the presence of dolomite, calcite, and quartz. After applying cutoffs, petrophysical parameters were evaluated. The estimated values for effective porosity, reservoir thickness, shale volume, and water and hydrocarbon saturation were 15.04 % to 17.29 %, 85.44 m to 122.41 m, 8.99 to 14.87 % and 12.97–27.10 %, and 72.9 % to 87.03 %, respectively. The outcome of this extensive research aids in accurately characterizing this region for the future prospectively of the field. This integrated study can be applied and adapted to other basins in New Zealand and the southwestern Pacific. The proposed workflow is also suitable for analogous scenarios in various basins worldwide.

Reconstruction of the Subsurface of Al-Hassa Oasis Using Gravity Geophysical Data

Al-Hassa city, located in Eastern Saudi Arabia, boasts the world’s largest oasis and the most expansive naturally irrigated lands. Historically, a total of 280 natural springs facilitated significant groundwater discharge and irrigation of agricultural land. Furthermore, the water in certain springs formerly had a high temperature. The spatial variability of the water quality was evident. At the same time, Al-Hassa Oasis is situated on the northeastern side of the Ghawar field, which is the largest conventional onshore oil field in the world in terms of both reserves and daily output (approximately 3.8 mmb/d). The aforementioned traits suggest an intricate subsurface that has not yet been publicly and thoroughly characterized. Due to the presence of significant cultural noise caused by agricultural and nearby industrial activities, a robust, easy-to-use, and accurate geophysical method (gravity) was used to cover an area of 350 km2, producing the 3D subsurface model of the study area. A total of 571 gravity stations were collected, covering the whole Al-Hassa Oasis and parts of the nearby semi-urban areas. The gravity data were corrected and processed, and a 3D inversion was applied. The resulting 3D geophysical subsurface modeling unveiled an intricate subterranean configuration and revealed lateral variations in density, indicating the presence of a potential salt dome structure, as well as fracture zones that serve as conduits or obstacles to the flow of the subsurface fluids. This comprehensive modeling approach offers valuable insights into the subsurface dynamics of the broader study area, enhancing our understanding of its qualitative tectonic and hydraulic features and their impacts on the area’s natural resources, such as groundwater and hydrocarbons.

Subsurface stratigraphy suggested by the layered ejecta craters in the Martian northern planitiae

The formation of layered ejecta craters (LECs) is generally thought to be related to the subsurface water-ice. The distribution of LECs in the northern regions of Mars (0 to 70° N) features both local and latitudinal variations, e.g., there are generally few LECs within volcanic areas, and the observed minimum and median diameters of the LECs decreases with increasing latitude. The observed minimum diameters of the LECs in the six planitiae within the northern lowlands suggest the roof depth of the LECs-forming layer has clear lateral variability, ranging between ∼90 m in the Arcadia Planitia and ∼280 m in the Isidis Planitia. The observed LECs are suggested to mainly form during the Amazonian, implying that the LEC-forming water-ice-bearing strata might be still present. The variation might reflect the comprehensive influence on the subsurface water-ice by the climate conditions and geologic events during the Amazonian period. The proposed subsurface stratigraphies in the six planitiae suggest that the total thickness of the LECs-related strata in the Chryse Planitia is the largest (∼4.6 km) and in the Arcadia Planitia is the smallest (∼1.6 km). There could be as much as ∼3.5–16 million km3 subsurface water stored in these strata, resulting in a ∼24–110 m global equivalent layer (GEL) of water. Given that water is a fundamental ingredient to life, its migration from surface to subsurface shifts the focus of its exploration. Both Arcadia and Utopia Planitiae represent ideal candidates as locations where both shallow and deep subsurface water-ice might be abundant, they represent high-priority targets for future missions at mid-latitudes for biosignature searching and in-situ resource utilization.

Sagittal Full-Spine vs. Sectional Cervical Lateral Radiographs: Are the Measurements of Cervical Alignment Interchangeable?

(1) Background: This study assessed the relationship between cervical spine parameters taken on standing full-spine lateral radiographic images compared to sectional lateral cervical radiographs. (2) Methods: Full-spine (FS) and sectional lateral cervical (LC) radiographs from four spine treatment facilities across the USA retrospectively provided data collected on 220 persons to assess the comparison of three sagittal cervical radiographic measurements between the two views. The measures included cervical lordosis using the absolute rotation angle from C2-C7, sagittal cervical translation of C2-C7, and atlas plane angle to horizontal. Linear correlation and R2 models were used for statistical comparison of the measures for the two views. (3) Results: The mean values of the three measurements were statistically different from each other: C2-C7 translation (FS = 19.84 ± 11.98 vs. LC = 21.18 ± 11.8), C2-C7 lordosis (FS = -15.3 ± 14.63 vs. LC = -18.32 ± 13.16), and atlas plane (FS = -19.99 ± 8.88 vs. LC = -22.56 ± 8.93), where all values were p < 0.001. Weak-to-moderate-to-strong correlations existed between the full-spine and sectional lateral cervical radiographic variables. The R2 values varied based on the measurement were R2 = 0.768 (p < 0.001) for sagittal cervical translation of C2-C7 (strong), R2 = 0.613 (p < 0.001) for the absolute rotation angle C2-C7 (moderate), and R2 = 0.406 (p < 0.001) for the atlas plane line (weak). Though a linear correlation was identified, there were consistent intra-person differences between the measurements on the full spine versus sectional lateral cervical radiographic views, where the full-spine view consistently underestimated the magnitude of the variables. (4) Conclusion: Key sagittal cervical radiographic measurements on the full spine versus sectional lateral cervical radiographic views show striking intra-person differences. The findings of this study confirm that full spine versus sectional lateral cervical radiographic views provide different biomechanical magnitudes of cervical sagittal alignment, and caution should be exercised by health care providers as these are not interchangeable. We recommend the LC view for measurement of cervical sagittal alignment variables.

An efficient ultrasonic full-matrix imaging method for industrial curved-surface components defect detection

Curved-surface components are ubiquitous in various industries. However, defects inside the components badly degenerate the mechanical properties and lead to premature failure. Ultrasonic full-matrix imaging is a prominent method in nondestructive testing, but there is still a lack of efficient imaging methods for cases with curved interfaces. To amend this problem, this article proposes an efficient ultrasonic full-matrix imaging method to detect the defects in the industrial curved-surface components. The wavefield extrapolation for full matrix capture (FMC) dataset is taken as the combined effect of the downward migration of the source and receiver wavefields in the wavenumber-frequency domain. Compensation terms corresponding to the source and receiver wavefields act on the wavefield of FMC in the spatial-frequency domain to make up for lateral variations of the sound speed. After the wavefield reconstruction, an efficient imaging condition is implemented to transform the wavefield into the defect image. The best compensation term is determined and the influence factors is investigated in the simulation experiments. Laboratory experiments furtherly validate the superiority of the proposed method over the existing methods, especially for the imaging efficiency which is twice of the phase shift migration (PSM) based method.

New fossiliferous sites from the mid-Cretaceous Tendrara dome (High Plateaus, Morocco): biostratigraphical, paleoenvironmental and paleogeographical implications

Sedimentological and stratigraphical studies in the Cretaceous series of the Tendrara dome led to the discovery of fossiliferous levels in the pre-Cenomanian and Cenomanian deposits, providing significant biostratigraphic, paleoenvironmental, and paleogeographic results. In the SW part of the Tendrara dome, the predominantly terrigenous deposits at the base of the pre-Cenomanian (Dekkar 1 Formation) yielded dinosaur eggshell fragments and charophytes, in particular Clavatoraceans, pointing to the Barremian-Aptian. Above this unit, two beds with fish remains were discovered in the Dekkar 2 Formation and at the base of the Dekkar 3 Formation, respectively. To the NE of the Tendrara dome, a fossiliferous site with diversified benthic fauna and abundant fish remains was discovered in a thin marly unit unconformably overlying the Middle Jurassic basement. Elasmobranch micro-remains indicate a Cenomanian age for this unit. The SW-NE correlations indicate a marked thickness reduction and lateral facies variations, implying a strong asymmetry in the Cretaceous paleogeography of the dome. The first continental and lagoonal basins of the Barremian-Aptian and Albian?-Cenomanian are located in the southern part of the Tendrara dome. The Cenomanian transgression, initiated from the south, gradually covered the dome, depositing reduced coastal plain elasmobranch-rich facies in its northern part. The Tendrara dome constituted a paleogeographic barrier, limiting the first transgressions of the Cenomanian sea. This paleostructure is part of an emerged area (Idrissides High) located between the Tethyan Ocean and the Saharan epicontinental sea.

Crustal structure of the central and southern Kyushu-Palau Ridge: Implications for intra-oceanic arc evolution

As an intra-oceanic arc, the Kyushu–Palau Ridge (KPR) is an exceptional place for studying intra-oceanic arc evolution and related magmatism using its interior velocity structure. Two wide-angle seismic profiles along the KPR provide an opportunity to unveil the crustal structure of its central and southern parts. Strong along-ridge variations are found either in P-wave velocity or crustal thickness of the central and southern KPR. The crustal thickness is between 6 and 12 km, with the velocity increasing from 4.0 km/s to 7.0 km/s from top to bottom. The upper crust has a large velocity gradient of ∼1 s−1 and a thickness between 2 and 4 km, while the lower crust has a lower velocity gradient of <0.2 s−1 and is characterized by inhomogeneous velocity anomalies and a substantial lateral variation in thickness. By comparing with the mature arc in the Izu-Bonin-Mariana arc in the east, we suggest the central and southern KPR is a thickened oceanic crust rather than a mature arc crust. Partly thick crust, velocity anomalies in the lower crust and upper mantle imply the crust is in a transition stage from oceanic crust to mature arc one. Compared with the quantitatively calculated seamount volumes, the lower crustal thickness is generally correlated well with the distribution of the seamount volumes along the central and southern KPR. These observations imply that the central and southern parts of the KPR are at the juvenile stage of arc evolution, with crustal growth primarily attributed to the lower crust thickening.

Investigation of leachate infiltration on groundwater using georesistivity and natural electric field method around Ojoou-Olayanju’s dumpsite, Ada, southwestern Nigeria

The investigation of groundwater within the dumpsite environment is highly important in geophysical study. This is because the extent of interaction between the aquiferous medium and the contaminated zone could pose a serious threat to the end users especially humans when consumed. This research aimed at investigating leachate infiltration and its potential influence on groundwater at Ojoou Olayanju's Dumpsite using combining geo-resistivity and natural electric field (NEF) methods. In this study, five dipole-dipole and five NEF measurements were obtained using the Omega resistivity meters and PQWT-150 equipment respectively. The dipole-dipole method was deployed to obtain a 2D near-surface pseudo-section, and the NEF method was used to obtain the frequency curve and profile maps of electric potential difference. The dipoledipole results revealed the lateral variation in the resistivity along the traverses, suggesting that the materials within this near-surface are heterogeneous, and the closely spaced contours' varying gradients indicate fracture, which would facilitate potential leachate filtration. The NEF results revealed curves, and a subsurface image with respect to depth and profile distance. The points of convergence signals on the frequency model correspond to a medium with low resistivity on the profile map. The conductive medium is seen as being saturated with leachate, which suggests that very large portion of the study area around the dumpsite has been contaminated by leachate. Conclusively, it was revealed that leachate filtration is evidence especially at the topsoil and due to the presence of fractured zones, the groundwater quality is at risk of contamination by continuous filtration of leachate

Application of seismic refraction and MASW methods for investigating the Spillway Fault trace along the western side of the Aswan High Dam, Egypt

AbstractAn earthquake of local magnitude ML = 4.6 occurred on November 7, 2010, 4.5 km northwest of the Aswan High Dam on the Spillway Fault. In the Aswan metropolitan region this earthquake was felt intensely. As no surface rupture was found, the focal mechanism and the distribution of seismic activity was one of the tools used for finding fault dimensions. The composite fault-plane solutions for the observed events on the Spillway Fault showed a left lateral strike-slip faulting with normal-fault component striking NNW-SSE. Also, remote sensing techniques were applied for the detection and identification of the geomorphology and geometry of the Spillway Fault. In this research, sub-surface layers and structures are delineated utilizing near-surface seismic techniques. Furthermore, the area’s supposed path and position of the Spillway Fault are also investigated. Two active seismic techniques, Seismic Refraction and Multi-Channel Analysis of Surface Waves (MASW), are utilized for recording near-surface seismic wave data at 9 sites. The seismic refraction profiles are conducted as a 2D cross-section on the trace of the detected Spillway Fault in the study area to evaluate the maximum depth of penetration of the P-wave for fault investigation. The constructed 2D seismic and structural sections from P-wave results show that the obtained average depth of about 30 m. In addition, the estimated P-wave velocities extend from 600 m/s to over 6500 m/s. Some lateral variation in the seismic wave velocities in all layers may represent fault zones. Moreover, the 1D MASW technique is conducted to estimate the velocities of the shear wave for the upper 30 m (Vs30) to provide the site classes and soil characteristics along both sides of the detected Spillway Fault trace in the study area. The calculated Vs30 values emphasized the idea of the existence of a normal dip-slip fault trace which divides the study area into two different lithological parts. The first part is located on the eastern side and characterized by almost class B (hard rock, according to NEHRP classification), while the other part is located to the west, and shows almost class type C (denoted as dense soil and soft rock soil).

Heterogeneous mantle effects on the behaviour of SmKS waves and outermost core imaging

SUMMARY Seismic traveltime anomalies of waves that traverse the uppermost 100–200 km of the outer core have been interpreted as evidence of reduced seismic velocities (relative to radial reference models) just below the core–mantle boundary (CMB). These studies typically investigate differential traveltimes of SmKS waves, which propagate as P waves through the shallowest outer core and reflect from the underside of the CMB m times. The use of SmKS and S(m-1)KS differential traveltimes for core imaging are often assumed to suppress contributions from earthquake location errors and unknown and unmodelled seismic velocity heterogeneity in the mantle. The goal of this study is to understand the extent to which differential SmKS traveltimes are, in fact, affected by anomalous mantle structure, potentially including both velocity heterogeneity and anisotropy. Velocity variations affect not only a wave's traveltime, but also the path of a wave, which can be observed in deviations of the wave's incoming direction. Since radial velocity variations in the outer core will only minimally affect the wave path, in contrast to other potential effects, measuring the incoming direction of SmKS waves provides an additional diagnostic as to the origin of traveltime anomalies. Here we use arrays of seismometers to measure traveltime and direction anomalies of SmKS waves that sample the uppermost outer core. We form subarrays of EarthScope's regional Transportable Array stations, thus measuring local variations in traveltime and direction. We observe systematic lateral variations in both traveltime and incoming wave direction, which cannot be explained by changes to the radial seismic velocity profile of the outer core. Moreover, we find a correlation between incoming wave direction and traveltime anomaly, suggesting that observed traveltime anomalies may be caused, at least in part, by changes to the wave path and not solely by perturbations in outer core velocity. Modelling of 1-D ray and 3-D wave propagation in global 3-D tomographic models of mantle velocity anomalies match the trend of the observed traveltime anomalies. Overall, we demonstrate that observed SmKS traveltime anomalies may have a significant contribution from 3-D mantle structure, and not solely from outer core structure.

Two‐dimensional SH‐wave and acoustic P‐wave full waveform inversion: A Midland Basin case study

AbstractBuilding a reliable S‐wave velocity model remains challenging from P‐wave land seismic data as multi‐parameter elastic full waveform inversion which updates P‐ and S‐wave velocities simultaneously is not yet widely adopted due to its computational cost, highly nonlinear nature and noise contamination from land seismic data. To overcome this challenge, we propose implementing the SH‐wave full waveform inversion to obtain the S‐wave velocity model. The value of this approach has been examined by applying two‐dimensional time‐domain SH‐ and P‐wave acoustic full waveform inversion to a nine‐component three‐dimensional seismic survey acquired in the Midland Basin. The nine‐component seismic survey was acquired by using both vertical and horizontal component vibrators and receivers, which enables us to apply both SH‐ and P‐wave full waveform inversion to the raw and preprocessed shot gathers along a two‐dimensional line. The SH‐ and P‐wave full waveform inversion applied to the raw shot gathers provide more detailed P‐ and S‐wave velocities compared to the velocities from the stacking velocity analysis. Additionally, there are no problematic artefacts in the inverted S‐wave velocity model, which shows the stability of the SH‐wave inversion. Through the SH‐wave full waveform inversion from the preprocessed shot gathers, we reveal the lateral velocity variations in the Grayburg–San Andres interval, which coincides with the depositional environment changes suggested by the previous regional study mainly based on well log and cores. These variations demonstrate that SH‐wave full waveform inversion can identify additional geological features that are not observed in the inverted P‐wave velocity model. The inverted P‐ and SH‐wave velocities are validated by the flattening of the events observed in the common image gathers after Kirchhoff depth migration. The SH‐wave migration image further reveals the irregular geometries of carbonates in the Spraberry formation. Vp/Vs ratios calculated from the independently inverted P‐ and S‐wave velocity models show strong lateral variations in the Wolfcamp interval (key‐producing interval), which could be caused by the organic content variations between the reservoir and non‐reservoir rocks. The inversion results demonstrate that the SH‐wave full waveform inversion can be used to provide an S‐wave velocity model that is comparable to the P‐wave velocity model derived from acoustic full waveform inversion, which is widely used for P‐wave velocity model building from P‐wave land seismic data.

Deciphering the sedimentary imprint of tsunamis and storms in the Lesser Antilles (Saint Martin): A 3500-year record in a coastal lagoon

Our study aims at identifying and characterizing tsunami and storm deposits by combining sedimentological, geochemical, and radiocarbon dating analyses. X-ray computed microtomography (micro-CT) is used to characterize the sedimentary fabric of the deposits. The study was conducted on a transect of oriented short sediment cores from a small lagoon at Saint Martin (Lesser Antilles), where both event types could be observed. Using micro-CT data of different event deposits, we were able to derive paleo-flow orientations and directions and compare the dynamics of the events to their environmental setting. The results showed that seven extreme-wave event (EWE) deposits, originating from hurricanes or tsunamis, had the highest Ca/Fe ratios, distinctive for marine input to the lagoon, and were characterized by coarser grains. Some also presented above the sandy part of the event deposit a siliciclastic geochemical signal rich in fine siliciclastic sediment with organic matter, which might correspond to backwash deposits. The thickness and frequency of these sandy EWE deposits in our cores decreased from the proximal to the distal zone in the center of the lagoon, with high lateral and vertical variability in sediment thickness. The youngest EWE deposit relates to the unprecedentedly powerful Category 5 Hurricane Irma, which barely left a fine layer of sand in the lagoon, demonstrating the lagoon's resilience to tropical storms. Five of the six other EWEs, characterized by micro-CT-derived sedimentary fabric, likely originated from tsunamis and occurred over the last 3500 years BP with a recurrence interval of 300–400 years, based on radiocarbon dating. A well-studied tsunami deposit from Pre-Columbian times (∼400 years BP) showed oriented sedimentary fabric, which could be used for paleo-flow reconstruction and EWE-type distinction in the Lesser Antilles, suggesting the use of an additional method for identifying EWE-type deposits in the geological record rather than other methods alone. To develop effective strategies for mitigating the natural hazards faced by communities in the coastal areas of the Lesser Antilles, it will be crucial to examine EWEs in the Caribbean region beyond historic times successfully.

Introducing a tailored site delineation approach to optimize the design of managed aquifer recharge surface spreading infrastructure

Managed aquifer recharge (MAR) is an emerging solution to effectively replenish overused groundwater resources, but high associated costs often hinder its uptake. There are several parameters to consider when determining the cost of MAR, including recharge types, scale of MAR schemes, land acquisition, operation and maintenance period, and the hydrogeological setting. Hydraulic conductivity and its spatial variability are the most significant hydrogeological parameters for predicting infiltration rates at MAR sites, but limited data availability makes accurate predictions often challenging. Hence, it is essential to increase data availability to optimize MAR efficiency and to better assess its economic performance. Therefore, a novel approach for MAR planning is presented in this article to efficiently enhance the data availability during the conceptual stage of MAR planning by combining site exploration, spatially resolved infiltration rate estimation, and cost analysis. The approach was applied at an actual MAR site in Vincenza, Italy, and incorporated an electromagnetic induction survey and direct push profiling as across-scale investigation methods to enhance site delineation and subsequent parameterization of the identified zones using laboratory analysis of collected samples. This MAR site is particularly suitable for its heterogeneous subsurface structure and potential to counteract groundwater depletion. The performed investigation resulted in an image of the lateral variation of subsurface conditions and allowed the delineation of three zones with different infiltration behavior, with 99.5% of infiltration occurring in just two zones, representing only 22.1% of the site. In this particular scenario, the overall cost per cubic meter of recharged water can be reduced by 59.1% by identifying and eliminating unfavorable zones. This study provides scientists and practitioners with a useful tool for MAR planning that can be applied to a wide range of sites with complex geology, thereby reducing water costs, minimizing the environmental impact of infrastructure construction, and reducing land use conflicts.

Continental Hotspots Tracks From an Analysis of GOCE Gravity Gradients Data

AbstractThis study aims to detect residual thermal anomalies along the path of a hotspot beneath the continental lithosphere using gradiometry data. We calculate perturbations in gravity vectors and gradients by employing a simple thermal anomaly model. The results highlight that the horizontal components of Bouguer gravity gradients, sensitive to source structure and elongation, exhibit substantial magnitudes, exceeding current detection levels. To improve signals from sources of different sizes and orientations, we apply a filtering approach involving wavelet transforms and rotated local frames. We also assess the impact of lateral crustal variations on the gravity field by introducing random density anomalies and a seismic crustal model, testing the method's ability to distinguish between these anomalies and thermal sources. Using Gravity field and steady‐state Ocean Circulation Explorer satellite data, we generated scale‐orientation diagrams aiming to identify signals aligned with the direction of plate movement and with a spatial scale of a few hundred kilometers (scale of the plumes). Maps of filtered Bouguer gravity gradients aligned with hotspot trajectories are generated for all continental hot plumes, with the position and the age of the intraplate volcanism to support the interpretation. While we do not show a clear signal in regions that are too tectonically complex or with lateral variations in the crust or at the lithospheric boundary correlated with the hotspot trace, we found observed signals with expected sign and amplitude over Hoggar, Tibesti, Darfour and Cameroon tracks in Africa and over Iceland and Jan Mayen in Greenland. However, determining corresponding mass source depths remains a complex task.

An Updated Fault Coupling Model Along Major Block‐Bounding Faults on the Eastern and Northeastern Tibetan Plateau From a Stress‐Constrained Inversion of GPS and InSAR Data

AbstractLarge block‐bounding faults on the Tibetan plateau are significant geological structures that accommodate tectonic movements and accumulate stress, leading to large earthquakes. Quantifying the interseismic slip deficit rate helps to better assess the earthquake potential. We combine available InSAR (2015–2020) and interseismic GPS data to determine fault coupling along 14 major block‐bounding faults. Spatially dense InSAR measurements remarkably improve the resolution of the coupling model. Combined with a GPS‐constrained block model, we examine the performance of the inversion approach with the stress constraint and the common Laplacian smoothing based on both synthetic tests and real data. We suggest that, for continental strike‐slip faults, adding the stress constraint can mitigate unphysical coupling distributions due to unreasonable assumptions or modeling artifacts, reducing the model uncertainty and improving the accuracy of the coupling model. This is particularly useful for segments featured by a highly heterogeneous distribution of coupling along the transition zone from locking to creeping region, partially‐coupling segment, and junction zone between main and subsidiary faults. We present a large‐scale fault coupling map along the major block‐bounding faults on the northeastern and eastern Tibetan plateau, highlighting the distinct degrees of fault coupling and lateral variations. The collage of coupling maps along different faults demonstrates the kinematic features over a broad time scale during earthquake cycles ranging from early to late interseismic phases, such as the segments ruptured during the 2001 Kokoxili earthquake and the 1920 Haiyuan earthquake.