Lateral Variations Research Articles

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 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.

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 (SSDS) 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 comprises 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 Common Diffusional Mechanism for Creep and Grain Growth in Polymineralic Rocks: Application to Lower Mantle Viscosity Estimates

AbstractIn a previous study (Okamoto & 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.

Lateral transport controls the tidally averaged gravitationally driven estuarine circulation: Tidal mixing effects

Abstract In classic models of the tidally averaged gravitationally driven estuarine circulation, denser salty oceanic water moves up the estuary near the bottom while less dense riverine water flows towards the ocean near the surface. Traditionally, it is assumed that the associated pressure gradient forces and salt advection are balanced by vertical mixing. This study, however, demonstrates that lateral (across the estuary width) transport processes are essential for maintaining the estuarine circulation. This is because for realistic estuarine bathymetry, the depth-integrated salt transport up the estuary is enhanced in the deeper estuary channel. A closed salt budget then requires the lateral transport of this excess salt in the deeper channel towards the estuarine flanks. To understand how such lateral transport affects the estuarine salt and momentum balances, we devise an idealized model with explicit lateral transport focusing on tidally-averaged lateral mixing effects. Solutions for the along-estuary velocity and salinity are non-dimensionalized to depend only on one single non-dimensional parameter, referred to as Fischer number, which describes the relative importance of lateral to vertical tidal mixing. For relatively strong lateral tidal mixing (greater Fischer number), salinity and velocity variations are predominantly vertical. For relatively weak lateral tidal mixing (smaller Fischer number), salinity and velocity variations are predominantly lateral. Overall, lateral transport greatly affects the estuarine circulation and controls the estuarine salinity intrusion length, which is demonstrated to scale inversely with Fischer number.

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.

Pulse-echo ultrasound attenuation tomography.

Objective.We present the first fully two-dimensional attenuation imaging technique developed for pulse-echo ultrasound systems. Unlike state-of-the-art techniques, which use line-by-line acquisitions, our method uses steered emissions to constrain attenuation values at each location with multiple crossing wave paths, essential to resolve the spatial variations of this tissue property.Approach.At every location, we compute normalized cross-correlations between the beamformed images that are obtained from emissions at different steering angles. We demonstrate that their log-amplitudes provide the changes between attenuation-induced amplitude losses undergone by the different incident waves. This allows us to formulate a linear tomographic problem, which we efficiently solve via a Tikhonov-regularized least-squares approach.Main results.The performance of our tomography technique is first validated in numerical examples and then experimentally demonstrated in custom-made tissue-mimicking phantoms with inclusions of varying size, echogenicity, and attenuation. We show that this technique is particularly good at resolving lateral variations in tissue attenuation and remains accurate in media with varying echogenicity.Significance.Based on a similar principle, this method can be easily combined with computed ultrasound tomography in echo mode (CUTE) for speed-of-sound imaging, paving the way towards a multi-modal ultrasound tomography framework characterizing multiple acoustic tissue properties simultaneously.

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.

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.

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.

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, MAR schemes scale, 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.

2-D Sn wave attenuation tomography beneath the Eastern Himalaya

SUMMARY The Arunachal and Bhutan Himalaya, which are tectonically distinct from other regions of the Himalaya, have a structure that is quite intricate. The eastern Himalayan segment is a component of the region where the Indian and Eurasian plates collided 50 Myr. The Indian Plate goes beneath the Eurasian Plate in the north, and in the eastern part of the region, the Indian Plate subducts under the Burmese Plate. Here, we studied the seismic attenuation of the uppermost mantle by measuring the quality factor of the Sn wave (SnQ) to understand the dynamics of the lithospheric mantle and the cause of the seismic anomalies found in this area. The upper mantle Q structure has significant lateral differences in Arunachal and the Bhutan Himalaya. Arunachal Himalaya’s central region is characterized by a very low $Q\, (\le 150)$. The successive low-high-low SnQ values in eastern Arunachal Himalaya near Siang region have been observed. The western Arunachal region, close to the Bhutan border, exhibits a contrast in Q values. We notice that low Q values (≤200) predominate in the central to eastern Bhutan Himalaya. The western part of Bhutan Himalaya exhibits relatively high $Q\, (\ge 200)$ values, mostly near Paro and Thimpu. Interestingly, a clear boundary between low and high Q has been observed near Kakthang thrust (KT) in the Bhutan Himalaya. We found significant lateral variation of frequency dependent parameter (η) across the study region. They range from 0.25 to 0.75, with low values (≤0.5) found mostly in the central Bhutan Himalaya and in a few isolated areas of the Arunachal Himalaya. Low Q and a relatively higher $\eta \, (\ge 0.5)$ might suggest that the scattering attenuation is the controlling mechanism for Sn wave attenuation in the upper mantle beneath Arunachal Himalaya. On the contrary, dominant low Q values across the central segment of the Bhutan Himalaya, along with a low to moderate body wave velocity and dominating low η values, subsequently corroborate that intrinsic attenuation is the dominant factor in the upper mantle of the central Bhutan Himalaya.

N-S variations of crustal structure beneath the central Tarim Basin from joint inversion of receiver functions, ambient seismic noise surface wave dispersion, and magnetotelluric data

To understand the complex crustal structure beneath the central Tarim Basin, we perform joint inversions of P-wave receiver functions, ambient seismic noise Rayleigh wave dispersion, and magnetotelluric data along a nearly N-S oriented linear array. In this contribution, we present a multistep flowchart to enhance the efficiency and reliability of the joint inversion. Synthetic examples indicate that this method is effective and appropriate for realistic conditions. By using these three complementary datasets, we can construct accurate 1-D models, particularly in the presence of thick overlying low-velocity layers. We apply this approach to data recorded at 45 stations within the Tarim Basin and generate quasi-2-D images that reveal distinct lateral variations along the N-S profile. In the near-surface region, we observe that the conductive and low-velocity sedimentary layer is thinner beneath uplift units and thicker beneath structural depressions, which correlates well with surface geology and previous studies. Additionally, we interpret several major faults, concealed crustal detachment layers, and discuss the geological implications. A conspicuous feature of the integrated images is the identification of a localized and deep-rooted doming structure beneath the Bachu uplift, characterized by high-velocity and high-resistivity anomalies. This structure also exhibits elevated density, intense magnetism and the thinnest overlying sedimentary layer. Consequently, we infer that these anomalies are probably linked to a hybrid basement consisting of plume-related intrusions and mafic dyke swarms.

DAS for 2-D MASW imaging: a case study on the benefits of flexible subarray processing

SUMMARY Distributed acoustic sensing (DAS) is a relatively new technology for recording the propagation of seismic waves, with promising applications in both engineering and geophysics. DAS's ability to simultaneously collect high spatial resolution waveforms over long arrays suggests that it is well-suited for near-surface imaging applications such as 2-D multichannel analysis of surface waves (MASWs), which require, at a minimum, long, linear arrays of single-component receivers. The 2-D MASW method uses a large number of sensor subarrays deployed along a linear alignment to produce 1-D shear-wave velocity (VS) profiles beneath each subarray. The 1-D VS profiles are then combined to form a pseudo-2-D VS image beneath the entire linear alignment that can be used for the purpose of identifying and characterizing lateral variations in subsurface layering. Traditionally, 2-D MASW is conducted using arrays consisting of either 24 or 48 geophones. While additional receivers could easily be incorporated into the testing configuration, it is rare for researchers and practitioners to have access to greater numbers of seismographs and geophones. When a limited number of geophones are available for deployment, there is a need to pre-determine the geophone spacing and subarray length prior to field data acquisition. Studies examining how the choice of subarray geometry impacts the resulting pseudo-2-D VS cross-sections have been largely limited to synthetic data. In response, this study utilizes DAS data to examine the effects of using various subarray lengths by comparing pseudo-2-D VS cross-sections derived from active-source waveforms collected at a well-characterized field site. DAS is particularly useful for 2-D MASW applications because the subarray geometry does not need to be determined prior to field data acquisition. We organize the DAS waveforms into multiple sets of overlapping MASW subarrays of differing lengths, ranging from 11 to 47 m, along the same alignment, allowing for direct comparison of the derived pseudo-2-D VS results at the site. We show that the length of the individual MASW subarrays has a significant effect on the resulting VS cross-sections, including the resolved location of a strong impedance contrasts at our study site, and evaluate the results relative to ground truth from invasive testing. Our results suggest that the choice of subarray length is important and should be carefully chosen to meet project-specific goals. Furthermore, analysts may consider using multiple subarray geometries during the data processing stage, as is made possible by DAS, to properly evaluate the uncertainty of 2-D MASW results. This study demonstrates the potential of using DAS to collect data for 2-D MASW in a manner that is efficient and flexible, and can be easily scaled up for use with very long arrays.

Antarctic Geothermal Heat Flow, Crustal Conductivity and Heat Production Inferred From Seismological Data

AbstractGeothermal heat flow is a key parameter in governing ice dynamics, via its influence on basal melt and sliding, englacial rheology, and erosion. It is expected to exhibit significant lateral variability across Antarctica. Despite this, surface heat flow derived from Earth's interior remains one of the most poorly constrained parameters controlling ice sheet evolution. To obtain a continent‐wide map of Antarctic heat supply at regional‐scale resolution, we estimate upper mantle thermomechanical structure directly from VS. Until now, direct inferences of Antarctic heat supply have assumed constant crustal composition. Here, we explore a range of crustal conductivity and radiogenic heat production values by fitting thermodynamically self‐consistent geotherms to their seismically inferred counterparts. Independent estimates of crustal conductivity derived from VP are integrated to break an observed trade‐off between crustal parameters, allowing us to infer Antarctic geothermal heat flow and its associated uncertainty.