Lateral Heterogeneity Research Articles

Lateral facies heterogeneity and carbonate productivity changes during the Lower Jurassic of the Traras Mountains (NW Algeria): local versus regional controls

Lateral facies heterogeneity and carbonate productivity changes during the Lower Jurassic of the Traras Mountains (NW Algeria): local versus regional controls

Substantial Global Radial Variations of Basalt Content Near the 660‐km Discontinuity

AbstractMid‐ocean ridges generate basalt and harzburgite, which are introduced into the mantle through subduction as a mechanical mixture, contributing to both lateral and radial compositional heterogeneity. The possible accumulation of basalt in the mantle transition zone has been examined, but details of the mantle composition below the 660‐km discontinuity (hereafter d660) remain poorly constrained. In this study, we utilize the subtle waveform details of S660S, the underside shear‐wave reflection off the d660, to interpret the seismic velocity, density, and compositional structure near, and particularly below, the d660. We identify a significant difference in S660S waveform shape in subduction zones compared to other regions. The inversion results reveal globally enriched basalt at the d660, with a notably higher content in subduction zones, consistent with the smaller impedance jump and S660S peak amplitude. The basalt fraction decreases significantly to less than 10% near 800‐km depth, forming a global harzburgite‐enriched layer and resulting in a steep seismic velocity gradient just below the d660, in agreement with 1D global reference models. The striking compositional radial variations near the d660 verify geodynamic predictions and challenge the applicability of homogeneous radial compositional models in the mantle. These variations may also affect the viscosity profile and, consequently, the dynamics at the boundary between the upper and lower mantle.

Hydro-thermal processes and deformation of highway embankment in the active layer in a high-latitude permafrost region of Inner Mongolia in Northeast China

Construction of embankments in the permafrost region significantly changes the heat exchange conditions and hydrothermal transport processes between permafrost and the external environment, causing changes in the state of permafrost under the embankment, which in turn affects the long-term stability of embankment impacts. Considering more complex forest environment and higher technical standard for expressway than ordinary highway, the hydrothermal and deformation characteristics of the embankment are investigated through a full-scale field experimental embankment of the Genhe-Labdalin highway. Further, the study delves into the influence of changes in the active layer thickness, hydrothermal processes, and water above the frozen layer on embankment stability. The main conclusions are as follows: 1) The permafrost table, temperature, moisture and deformation of the embankment showed lateral heterogeneity, with the three-former showing a “concave shape” and the left side (sunny slope) being lower than the right side (shady slope). 2) The permafrost table appears to be unconnected (thawing interlayer), creating preferential flow, thaw zones and even through-thaw zones. 3) Accompanied by the freezing and thawing process of the embankment, the deformation of the pavement is less delayed. These findings will be helpful for better understanding the hydrothermal characteristics of embankments in different frozen ground regions, and for providing important technical guidance to ensure the safe operation of engineering projects.

Wudalianchi volcanism and mantle dynamics in Northeast China: New insight from Pn and Sn tomography

We determine new tomographic models of Pn anisotropic velocity and Sn isotropic velocity in and around the Wudalianchi volcanic area by inverting high-quality Pn and Sn arrival times manually picked from waveforms recorded at the newly deployed high-dense WAVESArray portable seismic stations. Our high-resolution Pn and Sn velocity models reveal strong lateral heterogeneities in the uppermost mantle under the study region. The average Pn and Sn velocities in the uppermost mantle are 8.2 and 4.5 km/s, respectively. Both Pn and Sn velocity models exhibit obvious low-velocity (low-V) anomalies under the Wudalianchi and Keluo volcanoes, whereas under the Songliao basin distinct high-velocity (high-V) zones are revealed. In particular, our Pn model reveals two separate low-V anomalies under the Nuomihe and Halaha volcanic groups, suggesting that they have different deep origins. A large-scale L-shaped low-V zone exists under the Keluo, Wudalianchi, Erkeshan, and Xunke volcanoes, characterized by Pn-wave fast propagation directions (FPDs) parallel with the low-V zone, suggesting that these volcanoes may have the same deep origin. Furthermore, southeast-opened U-shaped Pn FPDs exist around the Wudalianchi volcano, whereas NE-SW FPDs appear under the Great Xing'an range, which are generally consistent with SKS splitting measurements. This feature may reflect lithosphere-asthenosphere coupling beneath the Wudalianchi volcano associated with horizontal flows in the big mantle wedge and compressional tectonics under the Great Xing'an range. These results shed new lights on the Wudalianchi volcanism and mantle dynamics beneath Northeast China.

Hybrid-driven MRF seismic inversion for gas sand identification: A case study in the Yinggehai Basin

The Yinggehai Basin displays anomalies characterized by heightened levels of temperature and pressure. It represents a depositional model of a submarine fan with gravity-driven flow, demonstrating significant lateral reservoir heterogeneity and intricate spatial distribution of the gas reservoir. Identifying gas formations using elastic parameters as indicators relies on stable seismic inversion results. This requires regularization to alleviate ill-posedness in the inverse problem and to construct model features of the subsurface medium. Markov Random Field (MRF)is an effective soft-constrained regularization method. It enhances the marginal features of inversion results by penalizing the objective function with the gradient of neighborhood points. However, the standard MRF method relies only on the parameter model-driven and has poor applicability in areas with strong lateral inhomogeneity or complex depositional processes. In this research paper, we propose a novel approach to seismic inversion that integrates MRF neighborhood drive and incorporates both seismic data and a parametric model. Multi-order MRF neighborhoods are constructed by using seismic data in the horizontal direction (including horizontal diagonal) and parametric model data in the vertical direction (including vertical diagonal). The model-driven results are also utilized to couple seismic data and improve the stability of the hybrid-driven MRF inversion. In addition, we select the P-impedance as the parameter for inversion due to its heightened sensitivity towards gas formation within the study region. Consequently, we utilize the inversion results to delineate the presence of sandstone in the target layer and discern any indications of gas formation. The implementation of this method in the field has demonstrated its capability to enhance the stability of inversion outcomes, effectively integrating the lateral consistency of seismic data with the vertical precision of parametric model data. This approach significantly improves reservoir heterogeneity characterization and enhances accuracy in identifying sandstone and gas.

Vertical high-velocity structures and seismic activity in western Shandong Rise, China: Case study inspired by double-difference seismic tomography

Abstract In this article, we collected the seismic phase arrival data (115°–120°E, 34°–39°N) of 2,833 local natural earthquakes above magnitude 1.0 recorded by 128 seismograph stations provided by the China Earthquake Networks Center covering the period from January 2008 to October 2023. We extracted the first arrival P- and S-wave arrival time data and obtained 26,351 P- and 26,349 S-wave absolute arrival times and 99,627 P-differential and 99,625 S-differential arrival times. Then, we determined 3-D P and S wave velocity structures under the Luxi region by using the double-difference tomography method. The results show lateral heterogeneities under the Luxi region, and the distribution of surface basins and uplift zones is associated with seismic velocities within the crust. The vertical section shows the widespread presence of vertical high-velocity bodies within the crust of the Luxi region, indicating seismological evidence of mantle upwelling in the region. Meanwhile, our imaging results also clearly indicate the presence of a significant low-velocity anomaly at depths of 13–17 km, which corresponds to the presence of a high-conductivity layer at that depth obtained from electromagnetic bathymetry. The transition zone between high and low velocities exists at a depth of 7–13 km, which plays an important role in the transfer of lithospheric stresses from the deep to the shallow part of the lithosphere. Combined with the results of the previous deep seismic wave velocity structure, petrology, and geological investigations, the deep dynamical background of the formation of slip tectonics and seismic mechanism in the Luxi area can be attributed to the go-slip activity of the Tan-Lu fault Zone, the subduction of the Pacific Plate, and the magmatism generated by mantle upwelling.

A method for structural framework variation of a multilayer field associated with complicated geology

Background. During the probabilistic resource base estimation, the largest uncertainties are associated with the wide range of the gross rock volume. In the course of the formation area work, the resulting maps of gross thickness are constructed incorrectly. Objective. To take into account facies zoning when varying the structural framework if working with the Achimov strata of a multilayer field. Materials and methods. The applied approach implies accounting of the errors in structural constructions in proportion to the power depending on facies zones. The authors proposed an additional stage of mapping using proportion maps to avoid obtaining geologically unrealistic local power anomalies. Results. A set of modified maps of the total thickness of clinocyclites was obtained, reflecting the key features of the geological structure of the modeling object, particularly facies differentiation and related lateral heterogeneity. Conclusions. The chosen approach to constructing maps of gross thickness allowed to obtain geologically realistic maps of effective thickness, which, in turn, formed the basis for mapping of productive volume.

Crustal structures and anisotropy in southeastern China: Continental collision and intraplate earthquakes

This study investigates the role of crustal structures and anisotropy in southeastern China regarding intraplate earthquakes within the context of concentrated and perturbed tectonic stresses. The region is characterized by a complex multiphase tectonic history,including the ongoing collision between the Philippine Plate and the South China Sea margin. Utilizing the teleseismic receiver function method, we offer a comprehensive depiction of the crustal structures and layered anisotropy in Guangdong and Fujian, two earthquake-prone areas. Our results reveal that the crustal anisotropy is stratified into two layers, with a slow axis perpendicular to the orogenic belt in the upper-middle crust, whereas the fast axis aligns perpendicularly at lower crust levels. The lateral heterogeneity and anisotropy of the crustal structure emphasize both the boundary delineation and differential movement between continental blocks. The observed increase in density and intensity of intraplate earthquakes across this boundary supports our conclusion regarding a connection between decoupled crustal deformation and seismicity. Specifically, the foliation of metamorphic rocks during the Mesozoic collision generated orogen-parallel fabrics in the upper-middle crust, while ductile shearing resulted in orogen-perpendicular fast directions at lower crust levels. In contrast, mismatch evolution between continental blocks reactivated by present collision provides stress conditions conducive to earthquake occurrence in southeastern China. Our findings provide valuable insights into mechanisms behind intraplate earthquakes as well as tectonic evolution within continental collision zones.Plain Language SummaryThe mechanism of intraplate earthquakes is a topic of great interest in the study of continental collision and regional tectonic evolution. The controlling factors behind these earthquakes are highly complex and uncertain. In this study, we aim to investigate the role of crustal structures and anisotropy in southeastern China within the context of intraplate seismic activity. This region, which is currently experiencing a collision between the Philippine Plate and the South China Sea Margin, has a multiphase tectonic history that includes remnants of paleo-Pacific subduction. By utilizing the receiver function method, we provide a detailed picture of crustal structures and layered anisotropy in Guangdong and Fujian, two earthquake-prone regions in southeastern China. Our findings reveal that anisotropy in the crust can be stratified into two layers: an upper-middle crust and a lower crust with a slow axis and a fast axis, both perpendicular to the orogeny, respectively. The lateral heterogeneity and anisotropy of the crustal structure highlight the boundary and differential movement between the continental blocks. The enhanced density and intensity of the intraplate earthquakes across this boundary provide evidence for a connection between decoupled crustal deformation and seismicity. Specifically, our observations suggest that (1) the Mesozoic continental collision induced metamorphic foliation in the upper-middle crust, (2) ductile shearing is characteristic of the lower crust, and (3) ongoing collision has led to mismatch evolution between continental blocks, which creates stress conditions favorable for earthquake occurrence.

Towards a new generation of fluvial facies models for the interpretation of ancient deposits, based on inter‐annual peak discharge variance of modern rivers

ABSTRACTFluvial facies models based on river planform are flawed and not fully fit for purpose. An alternative approach to classifying fluvial deposits is based on river discharge characteristics. A recent paper showed that the co‐efficient of variance of annual peak discharge (CVQp) correlates well with the characteristics of preserved alluvium for a suite of modern rivers. Here, the database of rivers is expanded and augmented, and the utility of this statistic confirmed. For modern rivers, plotting the running‐value of CVQp for a set interval (for example, 30 years) over a longer discharge record shows that the statistic varies both spatially and temporally. The CVQp can rise abruptly at a given location in response to a single flow event if that event is significantly larger than is common in that river. Such abrupt changes in CVQp are superimposed on longer‐term patterns of more gradual change stimulated by environmental factors or anthropogenic effects such as dam construction. If the period of calculation for CVQp is relatively short, ‘anomalous’ events appear more frequently. Sequential studies of modern rivers indicate that such anomalous events have a high probability of being preserved in the alluvial record of the river, and suggest that, over the timeframes of channel filling and abandonment, they may constitute the dominant (characteristic) portion of the record in high CVQp systems. Since flow events with strongly peaked hydrographs (flashy discharges) are more common in high CVQp systems and leave a record dominated by upper flow regime stratification styles, it follows that the alluvial record of high CVQp rivers will be dominated by upper flow regime stratification. In contrast, in lower CVQp systems the occurrence of an event that changes the running CVQp value significantly may form a distinct event deposit, and the temporarily changed CVQp is unrepresentative for the duration of the channel's life. The dominance of upper flow regime stratification, together with decreased preservation of identifiable macroform structure, increased lateral lithological heterogeneity, and in many cases increased preservation of in situ tree fossils in channel deposits of higher CVQp rivers, forms the basis for a series of new facies models for low, moderate, high and very high/ultra‐high CVQp alluvial systems and deposits. The new models can be applied to ancient successions either in combination with traditional, planform‐based models, or not. These models provide information about ancient palaeoenvironmental conditions, and changes in such conditions over time at various scales via stratigraphic investigations.

Characteristics and hydrocarbon accumulation model of Paleogene whole petroleum system in western depression of Qaidam Basin, NW China

Based on the oil and gas exploration in western depression of the Qaidam Basin, NW China, combined with the geochemical, seismic, logging and drilling data, the basic geological conditions, oil and gas distribution characteristics, reservoir-forming dynamics, and hydrocarbon accumulation model of the Paleogene whole petroleum system (WPS) in the western depression of the Qaidam Basin are systematically studied. A globally unique ultra-thick mountain-style WPS is found in the western depression of the Qaidam Basin. Around the source rocks of the upper member of the Paleogene Lower Ganchaigou Formation, the structural reservoir, lithological reservoir, shale oil and shale gas are laterally distributed in an orderly manner and vertically overlapped from the edge to the central part of the lake basin. The Paleogene WPS in the western depression of the Qaidam Basin is believed unique in three aspects. First, the source rocks with low organic matter abundance are characterized by low carbon and rich hydrogen, showing a strong hydrocarbon generating capacity per unit mass of organic carbon. Second, the saline lake basinal deposits are ultra-thick, with mixed deposits dominating the center of the depression, and strong vertical and lateral heterogeneity of lithofacies and storage spaces. Third, the strong transformation induced by strike-slip compression during the Himalayan resulted in the heterogeneous enrichment of oil and gas in the mountain-style WPS. As a result of the coordinated evolution of source-reservoir-caprock assemblage and conducting system, the Paleogene WPS has the characteristics of “whole process” hydrocarbon generation of source rocks which are low-carbon and hydrogen-rich, “whole depression” ultra-thick reservoir sedimentation, “all direction” hydrocarbon adjustment by strike-slip compressional fault, and “whole succession” distribution of conventional and unconventional oil and gas. Due to the severe Himalayan tectonic movement, the western depression of the Qaidam Basin evolved from depression to uplift. Shale oil is widely distributed in the central lacustrine basin. In the sedimentary system deeper than 2 000 m, oil and gas are continuous in the laminated limy-dolomites within the source rocks and the alga limestones neighboring the source kitchen, with intercrystalline pores, lamina fractures in dolomites and fault-dissolution bodies serving as the effective storage space. All these findings are helpful to supplement and expand the WPS theory in the continental lake basins in China, and provide theoretical guidance and technical support for oil and gas exploration in the Qaidam Basin.

Combined application of FTAN and cross-spectra analysis to ambient noise recorded by a microseismic monitoring network

A case study of seismic interferometry applied to a small microseismic monitoring network is here presented. The main objectives of this study are (i) to quantify the lateral variability of shear-wave velocities in the studied area, and (ii) to investigate the bias produced by noise directionality and non-stationarity in the velocity estimate. Despite the limited number of stations and the short-period character of the seismic sensors, the empirical Green’s functions were retrieved for all station pairs using two years of passive data. Both group and phase velocities were derived, the former using the widespread frequency-time analysis, the latter through the analysis of the real part of the cross-spectra. The main advantage of combining these two methods is a more accurate identification of higher modes, resulting in a reduction of ambiguity during picking and data interpretation. Surface wave tomography was run to obtain the spatial distribution of group and phase velocities for the same wavelengths. The low standard deviation of the results suggests that the sparse character of the network does not limit the applicability of the method, for this specific case. The obtained maps highlight the presence of a lower velocity area that extends from the centre of the network towards southeast. Group and phase velocity dispersion curves have been jointly inverted to retrieve as many shear-wave velocity profiles as selected station pairs. While the average model can be used for a more accurate location of the local natural seismicity, the associated standard deviations give us an indication of the lateral heterogeneity of seismic velocities as a function of depth. Finally, the same velocity analysis was repeated for different time windows in order to quantify the error associated to variations in the noise field. Errors as large as 4% have been found, related to the unfavorable orientation of the receiver pairs with respect to strongly directional noise sources, and to the very short time widows. It was shown that using a one-year time window these errors are reduced to 0.3%.

Shallow crustal structure of the northern Longmen Shan fault zone revealed by a dense seismic array with ambient noise analysis

The Longmen Shan (LMS) fault zone located in the southwest of China is not only the boundary tectonic zone of the Tibetan Plateau and South China block, but also an important part of the north–south seismic belt of China. To investigate the detailed crustal structure of the LMS fault zone, we deployed a dense seismic array of 151 short-period temporary stations in its northern section from March 16 to April 6, 2023. Continuous vertical component ambient noise waveforms recorded by these stations were cross-correlated for all station pairs, enabling the extraction of Rayleigh wave phase velocity dispersion curves across periods ranging from 0.5 s to 6 s. We then inverted these data to derive the crustal shear wave velocity using direct surface wave tomography. Meanwhile, a linear sub-array consisting of 54 stations and crossing the surface rupture of the 2008 Wenchuan earthquake was used to image the fault structure by applying the Transmitted Surface Wave Reverse Time Migration method (TSW-RTM). The ambient noise tomography results show that strong lateral heterogeneity exists in the shallow crust in the study area, with a high shear velocity in the northwest and a low shear velocity in the southeast. The low shear velocity is generally distributed between the Yingxiu-Beichuan fault (YBF) and the Guanxian-Jiangyou fault (GJF). The results of TSW-RTM indicate that the Wenchuan-Maoxian fault (WMF) and the YBF are both characterized with a steep dip to the northwest. Our results could serve as the important basis for the study of segmentation characteristics of the LMS fault and seismic hazard preparation and mitigation.

Using a 1D Radially Symmetric Coda Envelope Model for Robust Moment Magnitude (Mw) Estimation in Iraq’s Tectonically Diverse Zones

ABSTRACT Robust estimation of moment magnitude (Mw) can be challenging for Iraq due to the strong lateral heterogeneity across diverse tectonic zones. We aim to improve moment magnitude estimation by investigating the reliability of using a 1D coda envelope model in diverse tectonic zones of different lateral effects and offer a way forward for reliable estimates of Mw for small events that are difficult to waveform model. Iraq comprises two main tectonic zones: (1) the Outer platform, consisting of the northwestern Zagros fold-thrust belt and the Mesopotamian foredeep, and (2) the Inner Arabian platform which is overlain by the Iraqi desert. A simple 1D coda envelope model was used because coda waves have a low sensitivity to the source and path heterogeneity. Three separate coda calibrations were conducted to investigate the robustness of a single 1D calibration to fit the country: Whole-region calibration, Zagros calibration, and Mesopotamia calibration. In the whole-region calibration, we used stations from both the Zagros and Mesopotamia zones. In the two other calibration models, we used only stations that were in those particular zones. Ground-truth reference spectra derived from the coda spectral ratio method were used to constrain high-frequency site terms. There was no drastic difference when comparing the moment magnitudes calculated from the waveform modeling and the three calibration models. The results show that the 1D coda envelope model is a reliable method even for a region with diverse tectonic zones. Hence, we recommend using the whole-region calibration model for moment magnitude estimation that provides more complete path coverage and avoids biases introduced by path correction failures. The proposed calibration is a fundamental step in updating the comprehensive earthquake catalog and probabilistic seismic hazard assessments for Iraq.

Stress heterogeneity in the eastern Tibetan Plateau and implications for the present-day plateau expansion

The eastward expansion of the Tibetan Plateau has resulted in different earthquake types in the eastern Tibetan Plateau, but the mechanism remains unclear. Here, we construct a three-dimensional visco-elastoplastic finite element model considering the topography to investigate the influence of fault geometry and rheological heterogeneity on stress fields. In our best-fitting model, the minimum principal stress is nearly vertical around the southern Huya fault zone, which is adjacent to the Longmen Shan fault zone, due to the significant mid-lower crust lateral rheological heterogeneity, and the thrust stress regime accounts for the reverse fault and thrust-dominated earthquakes. In this scenario, the eastward horizontal motion of the mid-lower crust is obstructed and facilitates thrust faulting, suggesting the limited eastward expansion of the Tibetan Plateau. In contrast, the northern Huya fault zone, one of the terminal branches of the East Kunlun fault, accommodates the continuous eastward extrusion of the East Kunlun fault, where the stress regime under a more homogenized crust favors the strike-slip faulting process, along with the dominant strike-slip earthquakes. Moreover, the best-fitting of stress regime explains the thrust-dominated 2008 Ms. 8.0 Wenchuan and 2013 Ms. 7.0 Lushan earthquakes on the Longmen Shan fault zone. Combining geophysical and geodetic observations and model analyses, we propose that the hybrid deformation mode in the eastern Tibetan Plateau is accommodated by upper crustal shear and thrusting deformation and mid-lower crustal thickening driven by the gravitational potential energy gradient. Our results elucidate the mechanism for differences in strong historical earthquakes and, more importantly, isolate the effect of fault geometry from those of heterogeneous viscosity on crustal deformation and stress heterogeneity in the eastern Tibetan Plateau.

FEATURES OF THE LITHOLOGICAL STRUCTURE AND DEVELOPMENT OF RESERVOIR ROCKS AND HYDROCARBON RESERVOIRS IN THE MIDDLE DEVONIAN SEDIMENTS OF THE SOUTHERN PART OF THE TUZLIV DEPRESSION (DOBRUDJA FOREDEEP)

The aim of the work was to clarify the features of the lithological structure of the Middle Devonian sediments of the Zhovtoyarska-Tuzlivska and Biloliska-Zarichnenska prospective areas of the Dobrudja Foredeep and to determine its influence on the formation of reservoir rocks, traps and hydrocarbon reservoirs. The research was based on the results of geophysical research of wells (radioactive methods) in combination with fragmentary geological (lithological) data. As a result, 9 lithocycles of regressive nature were identified in the section for the first time (e-1 – e-5; g-1 – g-4). Each lithocycle is characterized by a two-membered structure and represents a separate productive horizon (PH). The lower parts of the latter are composed of porous and fractured reservoir rocks (limestones, dolomites, siltstones, sandstones), and the upper parts are composed of fluid-resistant packs (marls, anhydrites). It has been shown that cyclites are grouped into larger units (mesocyclites) of re- or progressive nature, which in different wells (in separate intervals) exhibit a certain similarity or difference in structure. Such spatial and age variability is caused, on the one hand, by the placement of wells in different facies zones of the sulfate-carbonate shelf, and on the other, by the processes of interaction of synsedimentary tectonic movements and paleoceanographic (sea level changes, etc.) factors. The consequence of this is the lateral lithological heterogeneity, which in the petrogeological aspect is manifested in the morphological features of the reservoirs and the character of the development of various types of reservoirs and ultimately determines, taking into account secondary rock changes, the prospects of specific productive horizons. Thus, for the Zhovtoyarska-Tuzlivska area, the vaulted trap is localized near borehole Zhovtoyarska-1 (e-1, e-2 productive horizons), according to other horizons it is fixed near borehole Zhovtoyarska-2. The development of lithological traps is predicted in e-1, e-3 and e-4 productive horizons, which is associated with the wedging of pore reservoirs in the direction from borehole Tuzlivska-2 to borehole Zhovtoyarska -2 and -1. For most productive horizons, the quality of reservoir rocks is expected to deteriorate in the same direction, which suggests the possibility of catagenetic screening of fluids. For the Biloliska-Zarichnenska area, along all productive horizons, trap vaults are localized near borehole Zarichnenska-1 or further south. The most capacious are the traps of horizons e-1, e-4 and e-5, composed mainly of pore-type collectors, the quality of which is low due to anhydritization. The best reservoir rocks are predicted in the Givetian sediments (g-2 – g-4 productive horizons) of borehole Zarichnenska -1.

Raft-like lipid mixtures in the highly coarse-grained Cooke membrane model.

Lipid rafts are nanoscopic assemblies of sphingolipids, cholesterol, and specific membrane proteins. They are believed to underlie the experimentally observed lateral heterogeneity of eukaryotic plasma membranes and implicated in many cellular processes, such as signaling and trafficking. Ternary model membranes consisting of saturated lipids, unsaturated lipids, and cholesterol are common proxies because they exhibit phase coexistence between a liquid-ordered (lo) and liquid-disordered (ld) phase and an associated critical point. However, plasma membranes are also asymmetric in terms of lipid type, lipid abundance, leaflet tension, and corresponding cholesterol distribution, suggesting that rafts cannot be examined separately from questions about elasticity, curvature torques, and internal mechanical stresses. Unfortunately, it is challenging to capture this wide range of physical phenomenology in a single model that can access sufficiently long length- and time scales. Here we extend the highly coarse-grained Cooke model for lipids, which has been extensively characterized on the curvature-elastic front, to also represent raft-like lo/ld mixing thermodynamics. In particular, we capture the shape and tie lines of a coexistence region that narrows upon cholesterol addition, terminates at a critical point, and has coexisting phases that reflect key differences in membrane order and lipid packing. We furthermore examine elasticity and lipid diffusion for both phase separated and pure systems and how they change upon the addition of cholesterol. We anticipate that this model will enable significant insight into lo/ld phase separation and the associated question of lipid rafts for membranes that have compositionally distinct leaflets that are likely under differential stress-like the plasma membrane.

Effects of 2.5-D ultra-low and ultra-high velocity zones on flip-reverse-stacking (FRS) of the ScS wavefield

SUMMARY Within the last decade, thin ultra-low velocity zone (ULVZ) layering, sitting directly on top of the core–mantle boundary (CMB), has begun to be investigated using the flip-reverse-stack (FRS) method. In this method, pre- and post-cursor arrivals that are symmetrical in time about the ScS arrival, but with opposite polarities, are stacked. This same methodology has also been applied to high velocity layering, with indications that ultra-high velocity zones (UHVZs) may also exist. Thus far, studies using the FRS technique have relied on 1-D synthetic predictions to infer material properties of ULVZs. 1-D ULVZ models predominantly show a SdS precursor that reflects off the top of the ULVZ and an ScscS reverberation within the ULVZ that arrives as a post-cursor. 1-D UHVZ models are more complex and have a different number of arrivals depending on a variety of factors including UHVZ thickness, velocity contrast, and lateral extent. 1-D modelling approaches assume that lower mantle heterogeneity is constant and continuous everywhere across the lower mantle. However, lower mantle features display lateral heterogeneity and are either finite in extent or display local thickness variations. We examine the interaction of the ScS wavefield with ULVZs and UHVZs in 2.5-D geometries of finite extent. We show that multiple additional arrivals exist that are not present in 1-D predictions. In particular, multipath ScS arrivals as well as additional post-cursor arrivals are generated. Subsequent processing by the FRS method generates complicated FRS traces with multiple peaks. Furthermore, post-cursor arrivals can be generated even when the ScS ray path does not directly strike the heterogeneity from above. Analysing these predictions for 2.5-D models using 1-D modelling techniques demonstrates that a cautious approach must be adopted in utilization and interpretation of FRS traces to determine if the ScS wavefield is interacting with a ULVZ or UHVZ through a direct strike on the top of the feature. In particular, traveltime delays or advances of the ScS arrival should be documented and symmetrical opposite polarity arrivals should be demonstrated to exist around ScS. The latter can be quantified by calculation of a time domain multiplication trace. Because multiple post-cursor arrivals are generated by finite length heterogeneities, interpretation should be confined to single layer models rather than to interpret the additional peaks as internal layering. Furthermore, strong trade-offs exist between S-wave velocity perturbation and thickness making estimations of ULVZ or UHVZ elastic parameters highly uncertain. We test our analysis methods using data from an event occurring in the Fiji-Tonga region recorded in North America. The ScS bounce points for this event sample the CMB region to the southeast of Hawaii, in a region where ULVZs have been identified in several recent studies. We see additional evidence for a ULVZ in this region centred at 14°N and 153°W with a lateral scale of at least 250 km × 360 km. Assuming a constant S-wave velocity decrease of −10 or −20 per cent with respect to the PREM model implies a ULVZ thickness of up to 16 or 9 km, respectively.

Geochemical characteristics and classification of Oligocene source rocks with different facies in Bozhong Sag, Bohai Bay Basin, East China

There have been many studies conducted on the Eocene source rocks (Shahejie Formation Mbrs 1 and 3) (E2s1 and E2s3) of the Bozhong Sag, Bohai Bay Basin in East China. However, Oligocene source rocks (Dongying Formation Mbrs 2 and 3) (E3d2 and E3d3) have also made significant contributions to petroleum accumulations in this region. Their organic matter (OM) enrichment mechanisms, environmental and ecological variations, and laterally variable facies remain subjects of debate, particularly in the lower Dongying Formation Mbrs 2 (E3d2L). Thus, we investigated the geochemical properties, depositional environments, and OM origins of E3d2L and E3d3 source rocks from seven prospecting wells in southwestern Bozhong Sag. Our findings indicate that E3d2L source rocks are of poor-to-medium quality, whereas E3d3 source rocks are classified as good-to-excellent quality. Notable differences in the geochemical characteristics between E3d2L and E3d3 source rocks were observed. Hierarchical cluster analysis identified five genetic organic facies. The E3d2L source rocks comprises facies D2-Ⅰ, D2-Ⅱ, D2-Ⅲ, and D2–D3, whereas E3d3 contains facies D3 and D2–D3, with D2–D3 being common to both source rocks. Additionally, an integrated OM enrichment model was established to evaluate Oligocene source rocks. The five genetic organic facies exhibit distinctive differences and some similarities regarding redox conditions, water salinity, and OM origins. Their geochemical characteristics within the same intervals display pronounced lateral heterogeneity from deep source rocks deposits to more marginal water depths. The established model provides guidance for the prediction and exploration of effective Oligocene source rocks in the Bohai Sea and research on the Paleogene and Neogene petroleum systems. Moreover, the integration of biomarkers and hierarchical cluster analysis methods for organic facies analysis has significant implications for predicting of source-rock heterogeneity in other lacustrine rift basins.

Nanoscale lipid domains determine the dynamic molecular portraits of mixed DOPC/DOPS bilayers in a fluid phase: A computational insight

Lateral heterogeneity, or mosaicity, is a fundamental property inherent to cell membranes that is crucial for their functioning. While microscopic inhomogeneities (e.g. rafts) are easily detected experimentally, lipid domains with nanoscale dimensions (nanoclusters of nanodomains, NDs) resist reliable characterization by instrumental methods. In such a case, important insight can be gained via computer modeling. Here, NDs composed of lipid's head groups in the mixed zwitterionic dioleoylphosphatidylcholine (DOPC) and negatively charged dioleoylphosphatidylserine (DOPS) bilayers were studied by molecular dynamics. A new algorithm has been developed to identify NDs. Unlike most similar methods, it implicitly considers the heterogeneous distribution of lipid head atomic density and does not require subjectively chosen parameters. In DOPS-rich membranes, lipids form more compact and stable NDs due to strong interlipid interactions. In DOPC-rich systems, NDs arise due to the “packing” effect of weakly bound lipid heads. The clustering picture is related to the physical properties of the bilayer surface: DOPS-rich systems show more pronounced surface heterogeneity of hydrophilic/hydrophobic regions compared to DOPC-rich ones. The results obtained are important for the effective quantitative characterization of the “dynamic molecular portrait” of a membrane surface – its “fingerprint” characterizing dynamical distribution of its physicochemical properties.

Machine Learning‐Based Hydrofacies Classification: Effects of Noise and Regularization

AbstractThis study utilizes an unsupervised ML approach, the expectation‐maximization (EM) algorithm using Gaussian Mixture Models (GMM), to integrate near‐surface geophysics measurements for hydrofacies classification. We examined the impact of noise and noise estimation on classification across two synthetic models with varying lateral heterogeneity, simulating and inverting resistivity and seismic data with noise levels ranging from minimal to very high. The algorithm proved robust in accurately reconstructing hydrofacies when noise was correctly estimated or not significantly underestimated, showing minimal misclassification in shallow hydrofacies. However, severe underestimation of noise during inversion led to increased misclassifications and artifact‐laden hydrofacies models, especially in shallow regions. Higher lateral heterogeneity lessened the negative impact of noise, slightly improving algorithm performance when noise was correctly estimated. We also explored the influence of geophysical measurement uncertainties on classification uncertainty through hydrofacies probability maps, noting the greatest impact when noise was underestimated. Additionally, we investigated the effect of the regularization trade‐off parameter on the hydrofacies classification and show how the performance of the algorithm can be evaluated in the absence of ground truth data using the average silhouette score of the classification with data obtained from a basement complex field site. We found that moderate regularization (λ = 200) yielded the best hydrofacies model, as indicated by the highest average silhouette score. Our findings underscore the effectiveness of unsupervised ML for facies classification and emphasizes the critical role of accurate noise characterization in geophysical data processing for enhancing the integration of subsurface heterogeneity information into hydrological models.