Search

Seafloor methane emissions can affect Earth’s climate and ocean chemistry. Vast quantities of methane formed by microbial decomposition of organic matter are locked within gas hydrate and free gas on continental slopes, particularly in large areas with high sediment accumulations such as deep-sea fans. The release of methane in slope environments has frequently been associated with dissociation of gas hydrates near the edge of the gas hydrate stability zone on the upper slope, with discharges in greater water depths less understood. Here we show, using data from the Rio Grande Cone (western South Atlantic), that the intrinsic, gravity-induced downslope collapse of thick slope sediment accumulations creates structures that serve as pathways for gas migration, unlocking methane and causing seafloor emissions via giant gas flares in the water column. The observed emissions in the study region (up to 310 Mg year−1) are three times greater than estimates for the entire US North Atlantic margin and reveal the importance of collapsing sediment accumulations for ocean carbon cycling. Similar outgassing systems on the Amazon and Niger fans suggest that gravity tectonics on passive margins is a common yet overlooked mechanism driving massive seafloor methane emissions in sediment-laden continental slopes.

Several studies have monitored crustal seismic velocity changes and attempted to relate them to the stress state and physical properties in volume embedding fault systems. The aim is to provide constraints on fault system dynamics and earthquake triggering mechanisms. Here, we reconstruct the spatiotemporal (4D) seismic velocity images of volume embedding the Irpinia fault system (IFS, South Italy), which originated the 1980 Ms 6.9 multi-segmented ruptures. By inverting data from more than ten years of continuous seismicity monitoring, we retrieved time-constant velocity anomalies, whose shapes correlate well with crustal lithology, while time-changing (up to 20%) velocity anomalies are mapped in the central region. Here, the Vp-to-Vs changes at depths of 1–5 km and 8–12 km correlate well with groundwater recharge and geodetic displacement during the same time interval. This correlation provides evidence for the existence of pulsating, pore pressure changes induced by groundwater recharge processes in a deep volume (8–12 km of depth), fractured and saturated with a predominant gas phase (likely CO2). We suggest that tomographic measurements of the Vp-to-Vs spatiotemporal changes are a suitable proxy to track the pore pressure evolution at depth in highly sensitive regions of fault systems.

Slope stability is a prominent problem for the efficient application and promotion of highwall mining technology, especially when mining residual coal under high and steep end-slope conditions. This study proposes the concept of target time pillar strength based on the required coal pillar service time. Creep tests were performed to measure the time-varying properties of coal shear strength parameters under different loads, and a time-varying function was established by regression. The highwall mining length is divided into three categories based on discontinuous structural plane theory, including goaf, yielding, and elastic zones, all of which are considered to have resistances against shear stress. The basal coal seam is prone to weakening owing to the weight of overlying strata, which may shift the slope failure mode from circular to sliding along the weak layer. Numerical modeling was used to study the influence of the bearing stress and target time strength on the development of the yielding zone at the coal pillar ribs. The coefficients of the three zones were determined, and the temporal and spatial evolution patterns of the shear strength parameters of the weak layer were acquired. A slope stability calculation method is proposed based on rigid body-limit equilibrium theory that can quantify the influence of highwall mining operations on slope stability, which is significant for popularizing highwall mining technology.

The sandy land leaks water and fertilizer, and is seriously degraded, while the soft rock has a special depression structure, which plays a role in retaining water and fertilizer. The application of soft rock new material to sand reclamation can improve the ecological environment and ensure the quality of basic cultivated land. The soft rock and sand were mixed in different volume ratios (1:0, 11:1, 5:1, 4:1, 3:1, 2:1, 7:5, 1:1, 5:7, 1:2, 1:3, 1:4, 1:5, 1:11, 0:1) to prepare the composite soil, and its Raman spectrum characteristics, microstructure, texture composition and potato yield were studied. The results show that there are more silt and clay particles in the soft rock and more coarse particles in the sand. The peak position of the sand is 464.5 cm−1. With the increase of the content of the soft rock, the peak position decreases gradually. When the content of the soft rock accounts for more than 50%, the soil structure collapses and also becomes compact, at the same time the compressive stress is generated between the soil particles. When the ratio of soft rock to sand is 1:1, the soil texture is loam. The potato yield of the soil with the ratio of 1:5 of soft rock and sand cover increases significantly by 4.89–37.31% and 4.08–35.95% compared with that of 1:1 and 1:2 compound soils. Under the condition of 1:1 and 1:2 compounded ratio of soft rock and sand, there are more cementitious materials between the soil particles generated. The compounded ratio 1:5 is most suitable ratio for potato growth of local economic crop. The results confirmed that the Raman spectroscopy characteristics of SiO2 molecules can be used to study the cementation force between composite soil particles. When the compound ratio is 1:5, the soil improvement of Mu Us sandy land and the high yield of potatoes can be achieved, which could also provides a theoretical basis for sandy land remediation.

Scale dependence of surface roughness of natural rock joints has long been an outstanding issue in rock mechanics. Controversial results were reported by various studies; however, the nature of scale dependency and the underlying mechanism for the conflicting observations remain unclear. Rock joints at different scales characterise two-order asperities, namely, waviness and unevenness; thus understanding how the individual roughness of waviness and unevenness vary as the joint size increases from the laboratory-scale to the large-scale is crucial for revealing the scale effect mystery. Here we digitise three natural granite joint surfaces with the same dimension of 1000 mm × 1000 mm through a high-resolution, three-dimensional scanner. Waviness and unevenness of each rock joint surface are quantitatively separated by selecting an appropriate sampling interval. The respective fractal dimensions of waviness and unevenness of joint surfaces sized from 100 mm × 100 mm to 1000 mm × 1000 mm are estimated through an improved roughness-length method. We find that the fractal dimensions of two-order roughness are scale-dependent but without generalised trends. The stationarity threshold beyond which the scale-dependency of roughness vanishes is absent for all the three joint samples, suggesting that the roughness of natural rock joints be assessed at the specific scale of the rock mass in-situ. We reveal that previous controversial results regarding scale effect are likely due to the composition of the roughness scaling of waviness and unevenness. Thus, accurate stability analysis of rock-engineering projects should consider separate contributions of multi-order asperities across scales to the strength and deformation of jointed rock masses.

Acoustic emission (AE) characterization is an effective technique to indirectly capture the failure process of quasi brittle rock. In previous studies, both experiments and numerical simulations were adopted to investigate the AE characteristics of rocks. However, as the most popular numerical model, the moment tensor model (MTM) cannot be constrained by the experimental result because there is a gap between MTM and experiments in principle, signal processing and energy analysis. In this paper, we developed a particle-velocity-based model (PVBM) that enabled direct monitoring and analysis of the particle velocity in the numerical model and had good robustness. The PVBM imitated the actual experiment and could fill in gaps between the experiment and MTM. AE experiments of marine shale under uniaxial compression were carried out, and the results were simulated by MTM. In general, the variation trend of the experimental result could be presented by MTM. Nevertheless, the magnitudes of AE parameters by MTM presented notable differences of more than several orders of magnitude compared with those by the experiment. We sequentially used PVBM as a proxy to analyse these discrepancies and systematically evaluate the AE characterization of rocks from the experiment to numerical simulation, considering the influence of wave reflection, energy geometrical diffusion, viscous attenuation, particle size and progressive deterioration of rock material. The combination of MTM and PVBM could reasonably and accurately acquire AE characteristics of the actual AE experiment of rocks by making full use of their respective advantages.

Diamonds grown by high pressure high temperature process (HPHT) are usually characterized by yellow color and high contents of nitrogen. Introduction of Ti decreases nitrogen content in diamond. Understanding the formation of nitrogen-poor diamond is very important not for the progress of HPHT process only, but because these diamond varieties represent the rare natural stones, although their crystallization conditions have not been clarified yet. Here we studied the composition of fluid phase in synthetic diamonds. The experiments were performed using a high-pressure apparatus BARS at pressures 5.5–6.0 GPa and temperatures 1350–1400 °C. It was found that introduction of metallic Ti leads to concentration of nitrogen mainly as nitrogenated hydrocarbons. The hypothesis that elucidates the formation of low-nitrogen diamond in Fe–Ni is proposed: the presence of Ti leads to an increase of hydrogen fugacity in the metal melt which drastically reduces the nitrogen solubility. As a result, nitrogen concentrates in the form of complex hydrocarbon compounds, while diamond grows colorless and characterized by very low nitrogen content. It is suggested that the proposed mechanism acts the same way in the presence of other metals which are strong reducing agents.

Temperature is used to trace ocean density variations, and reveals internal waves and turbulent motions in the deep ocean, called ‘internal motions.’ Ambient temperature detected by geophysical differential pressure gauges (DPGs) may provide year-long, complementary observations. Here, we use data from four DPGs fixed on the ocean bottom and a high-resolution temperature sensor (T-sensor) 13 m above the seafloor as a square-kilometer array deployed offshore ~ 50 km east of Taiwan facing the open Pacific Ocean to examine the impact of temperature on DPG signals related to internal motions. The DPG signals correlate with T-sensor temperature variations between 0.002 and 0.1 mHz, but have time shifts partially caused by slow thermal conduction from the ambient seafloor to the DPG chamber and partially by internal motion propagation time across the array. Applying beamforming-frequency-wavenumber analysis and linear regression to the arrayed T-sensor and DPG data, we estimate the propagating slowness of the internal motions to be between 0.5 and 7.4 s m−1 from the northwest and northeast quadrants of the array. The thermal relaxation time of the DPGs is within 103–104 s. This work shows that a systematic scan of DPG data at frequencies

Landslides are considered as one of the most devastating natural hazards in Iran, causing extensive damage and loss of life. Landslide susceptibility maps for landslide prone areas can be used to plan for and mitigate the consequences of catastrophic landsliding events. Here, we developed a deep convolutional neural network (CNN–DNN) for mapping landslide susceptibility, and evaluated it on the Isfahan province, Iran, which has not previously been assessed on such a scale. The proposed model was trained and validated using training (80%) and testing (20%) datasets, each containing relevant data on historical landslides, field records and remote sensing images, and a range of geomorphological, geological, environmental and human activity factors as covariates. The CNN–DNN model prediction accuracy was tested using a wide range of statistics from the confusion matrix and error indices from the receiver operating characteristic (ROC) curve. The CNN–DNN model was evaluated comprehensively by comparing it to several state-of-the-art benchmark machine learning techniques including the support vector machine (SVM), logistic regression (LR), Gaussian naïve Bayes (GNB), multilayer perceptron (MLP), Bernoulli Naïve Bayes (BNB) and decision tree (DT) classifiers. The CNN–DNN model for landslide susceptibility mapping was found to predict more accurately than the benchmark algorithms, with an AUC = 90.9%, IRs = 84.8%, MSE = 0.17, RMSE = 0.40, and MAPE = 0.42. The map provided by the CNN–DNN clearly revealed a high-susceptibility area in the west and southwest, related to the main Zagros trend in the province. These findings can be of great utility for landslide risk management and land use planning in the Isfahan province.

The Central-South Chile margin is an excellent site to address the changes in the gas hydrate system since the last deglaciation associated with tectonic uplift and great earthquakes. However, the dynamic of the gas hydrate/free gas system along south central Chile is currently not well understood. From geophysical data and modeling analyses, we evaluate gas hydrate/free gas concentrations along a seismic line, derive geothermal gradients, and model past positions of the Bottom Simulating Reflector (BSR; until 13,000 years BP). The results reveal high hydrate/free gas concentrations and local geothermal gradient anomalies related to fluid migration through faults linked to seafloor mud volcanoes. The BSR-derived geothermal gradient, the base of free gas layers, BSR distribution and models of the paleo-BSR form a basis to evaluate the origin of the gas. If paleo-BSR coincides with the base of the free gas, the gas presence can be related to the gas hydrate dissociation due to climate change and geological evolution. Only if the base of free gas reflector is deeper than the paleo-BSR, a deeper gas supply can be invoked.

Nitrate deposits are rare worldwide, especially potassium nitrate deposits; furthermore, their genesis remains disputed. There is a rare salt-lake type potassium nitrate deposit in the Dawadi area of Lop Nor at the eastern margin of the Tarim Basin, and the ore bodies show coexisting solid and liquid phases. Additionally, there are large sulphate-type potash deposits in the adjoining Luobei Depression, south of the Dawadi area. To determine why there are two different types of potash deposits in adjacent depressions with similar climates, field geological surveys were conducted and samples collected. It was found that the Tertiary clastic layer at the periphery of the Dawadi deposit was rich in high-salinity brine, with nitrate contents of up to 495–16,719 mg/L, much higher than those in the Luobei Depression, 1–35 mg/L. Additionally, a type of deep hydrothermal (Ca–Cl) brine was found in the fault zones, with nitrate contents of up to 8044 mg/L, dozens of times greater than that of ordinary groundwater. Using comprehensive analysis and research, we concluded that the Dawadi and Luobei depressions belong to different hydrological systems with no connection between them; thus, the two deposits belong to different metallogenic systems. Furthermore, groundwater played an important role in the mineralization of the potassium nitrate deposit, and a deep source may have been an important source of the ore-forming materials. The fault system widely developed in Lop Nor provides favorable channels for deep hydrothermal recharge, and the groundwater and deep hydrothermal brine could provide the source for the nitrate mineralization in the Dawadi Depression through water–rock reactions.

A large amount of high-quality teleseismic data is used for common conversion point (CCP) stacking of receiver functions in the Longmenshan area. The results show that a large-scale high-velocity anomaly or lithospheric delamination can completely destroy upper mantle discontinuities or erase the phase boundary of olivine, which is a very important finding and can be used to assess stagnant slabs in the mantle transition zone globally. The deepening region of the 660 km discontinuity beneath the Songpan-Ganzi terrane might indicate that the large-scale high-velocity anomaly in the mantle transition zone is a cold domain and can affect the topography of upper mantle discontinuities.

The development of predictive maps for geothermal resources is fundamental for its exploration across Nigeria. In this study, spatial exploration data consisting of geology, geophysics and remote sensing was initially analysed using the Shannon entropy method to ascertain a correlation to known geothermal manifestation. The application of statistical index, frequency ratio and weight of evidence modelling was then used for integrating every predictive data for the generation of geothermal favourability maps. The receiver operating/area under curve (ROC/AUC) analysis was then employed to ascertain the prediction accuracy for all models. Basically, all spatial data displayed a significant statistical correlation with geothermal occurrence. The integration of these data suggests a high probability for geothermal manifestation within the central part of the study location. Accuracy assessment for all models using the ROC/AUC analysis suggests a high prediction capability (above 75%) for all models. Highest prediction accuracy was obtained from the frequency ratio (83.3%) followed by the statistical index model (81.3%) then the weight of evidence model (79.6%). Evidence from spatial and predictive analysis suggests geological data integration is highly efficient for geothermal exploration across the middle Benue trough.

During the blasting excavation of deep underground caverns, the effects of the structural surface on crack propagation are usually considered in addition to the clamping effects of high in situ stress. Based on the notched borehole and timing sequence control (TSC) fracture blasting method, this paper studies the effects of different borehole shapes on the degree of damage of the surrounding rock and profile flatness of the rock anchor beams and the effects of different filled joint characteristics on the blasting crack propagation rules. The results show that the damage depth of the surrounding rocks by round hole smooth blasting is approximately twice that by notched hole smooth blasting, by which the profile formed is flatter. The notched primary borehole (PBH) remains a strong guidance for crack propagation in a rock mass with filled joints, while the stress concentration effects of the round target borehole (TBH) cannot fully guide the cracks until they fall within a certain distance between the PBH and TBH. It is favourable for cracks to propagate along the lines between boreholes with larger filled joint strengths and larger angles between boreholes.

The crystallization of the magma ocean resulted in the present layered structure of the Earth’s mantle. An open question is the electronic spin state of iron in bridgmanite (the most abundant mineral on Earth) crystallized from a deep magma ocean, which has been neglected in the crystallization history of the entire magma ocean. Here, we performed energy-domain synchrotron Mössbauer spectroscopy measurements on two bridgmanite samples synthesized at different pressures using the same starting material (Mg0.78Fe0.13Al0.11Si0.94O3). The obtained Mössbauer spectra showed no evidence of low-spin ferric iron (Fe3+) from the bridgmanite sample synthesized at relatively low pressure of 25 gigapascals, while that directly synthesized at a higher pressure of 80 gigapascals contained a relatively large amount. This difference ought to derive from the large kinetic barrier of Fe3+ rearranging from pseudo-dodecahedral to octahedral sites with the high-spin to low-spin transition in experiments. Our results indicate a certain amount of low-spin Fe3+ in the lower mantle bridgmanite crystallized from an ancient magma ocean. We therefore conclude that primordial bridgmanite with low-spin Fe3+ dominated the deeper part of an ancient lower mantle, which would contribute to lower mantle heterogeneity preservation and call for modification of the terrestrial mantle thermal evolution scenarios.

The evolution of the Pearl River Estuary (PRE), China in recent decades has been dominated by human activities. Historical admiralty charts and remote sensing images indicated that from 1936 to 2017, the tidal flat area and water area decreased by 23.6 × 107 m2 and 60.7 × 107 m2, respectively. The average advancing rate of the coastline of the PRE to the sea from 1972 to 2017 reached approximately 64.8 m/year, which is several times or even dozens of times that since the mid-Holocene. Land reclamation was the main reason for the dramatic changes in the water area and coastline. Although the water volume of the PRE showed a decreasing trend from 1936 to 2017, the water volume reduction rates for 1996–2005 and 2005–2017 were only 29% (1.27 × 107 m3/year) and 12% (0.53 × 107 m3/year), respectively, of that for 1936–1972. The combined influences of channel dredging, sand mining, and sediment load reduction caused by dam construction have contributed to this change. From the perspective of the filling up of the estuary, channel dredging, sand mining, and dam construction in the river basin are beneficial for prolonging the life of the estuary.

Landslide dam outburst floods have a significant impact on landform evolution in high mountainous areas. Historic landslide dams on the Yigong River, southeastern Tibet, generated two outburst superfloods > 105 m3/s in 1902 and 2000 AD. One of the slackwater deposits, which was newly found immediately downstream of the historic dams, has been dated to 7 ka BP. The one-dimensional backwater stepwise method gives an estimate of 225,000 m3/s for the peak flow related to the paleo-stage indicator of 7 ka BP. The recurrence of at least three large landslide dam impoundments and super-outburst floods at the exit of Yigong Lake during the Holocene greatly changed the morphology of the Yigong River. More than 0.26 billion m3 of sediment has been aggraded in the dammed lake while the landslide sediment doubles the channel slope behind the dam. Repeated landslide damming may be a persistent source of outburst floods and impede the upstream migration of river knickpoints in the southeastern margin of Tibet.

Geo-materials may present varying mechanical properties under different stress paths, especially for tunnel excavation, which is typically characterized by the decreased radial stress and increased axial stress during the complex loading and unloading process. This study carried out a comparative analysis between the loading and unloading model testing, which was then combined with PFC2D simulation, aiming to reveal the fracture propagation pattern, microscopic stress and force chain distribution of the rock mass surrounding the tunnel. Comparisons of extents and development of tensile strain between loading and unloading testing results were made. The overall stability, the integrity of rock mass, and the failure pattern transition under loading and unloading processes were systematically examined. In addition, for the two unloading cases with different vertical stresses imposed, the failure patterns were both identified as the collapse of the V − shaped extruded sidewall, due to the coupling of the shear failure and the vertical tensile failure in the sidewall wedge.

The 2019 Ridgecrest, California seismic sequence, including an Mw6.4 foreshock and Mw7.1 mainshock, represent the largest regional seismic events within the past 20 years. To obtain accurate coseismic fault-slip distribution, we used precise positioning data of small earthquakes from January 2019 to October 2020 to determine the dip parameters of the eight fault geometry, and used the Interferometric Synthetic Aperture Radar (InSAR) data processed by Xu et al. (Seismol Res Lett 91(4):1979–1985, 2020) at UCSD to constrain inversion of the fault-slip distribution of both earthquakes. The results showed that all faults were sinistral strike-slips with minor dip-slip components, exception for dextral strike-slip fault F2. Fault-slip mainly occurred at depths of 0–12 km, with a maximum slip of 3.0 m. The F1 fault contained two slip peaks located at 2 km of fault S4 and 6 km of fault S5 depth, the latter being located directly above the Mw7.1hypocenter. Two slip peaks with maximum slip of 1.5 m located 8 and 20 km from the SW endpoint of the F2 fault were also identified, and the latter corresponds to the Mw6.4 earthquake. We also analyzed the influence of different inversion parameters on the fault slip distribution, and found that the slip momentum smoothing condition was more suitable for the inversion of the earthquakes slip distribution than the stress-drop smoothing condition.

Detrital zircon grains within four (4) deposits including two (2) Cretaceous and two (2) Paleogene/Neogene kaolins in Nigeria were analysed using U–Pb LA-SFICP-MS to determine their provenance. The zircon textures were dominated by xenocrystic cores and oscillatory zoning in the Cretaceous and Paleogene/Neogene kaolins, respectively. The Th/U ratios obtained for the detrital zircon grains within the kaolins were predominantly within known values for rocks with magmatic origin. The age populations obtained for the detrital zircon grains were dominated by values from 529 to 978 Ma within the Neoproterozoic, followed by values from 1754 to 2497 Ma of the Paleoproterozoic. Detrital zircon ages obtained between 553.2 ± 6.2 and 583.5 ± 2.0 Ma represent part of the minimum provenance ages for the primary minerals that were kaolinised. The Cretaceous–Paleogene/Neogene kaolins were derived from parent rocks of Eburnean and Pan African ages within the Western and Northern Nigeria Basements.

We image the lateral variations in the Moho depths and average crustal composition across the Kumaon–Garhwal (KG) Himalaya, through the H–K stacking of 1400 radial PRFs from 42 three-component broadband stations. The modelled Moho depth, average crustal Vp/Vs, and Poisson’s ratio estimates vary from 28.3 to 52.9 km, 1.59 to 2.13 and 0.17 to 0.36, respectively, in the KG Himalaya. We map three NS to NNE trending transverse zones of significant thinning of mafic crust, which are interspaced by zones of thickening of felsic crust. These mapped transverse zones bend toward the north to form a NE dipping zone of maximum changes in Moho depths, below the region between Munsiari and Vaikrita thrusts. The 1991 Mw6.6 Uttarakashi and 1999 Mw6.4 Chamoli earthquakes have occurred on the main Himalayan thrust (MHT), lying just above the mapped zone of maximum changes in Moho depths. Modelled large values of average crustal Vp/Vs (> 1.85) could be attributed to the high fluid (metamorphic fluids) pressure associated with the mid-crustal MHT. Additionally, the serpentinization of the lowermost crust resulted from the continent–continent Himalayan collision process could also contribute to the increase of the average crustal Vp/Vs ratio in the region.

Plate bending-related normal faults (i.e. bend-faults) develop at the outer trench-slope of the oceanic plate incoming into the subduction zone. Numerous geophysical studies and numerical simulations suggest that bend-faults play a key role by providing pathways for seawater to flow into the oceanic crust and the upper mantle, thereby promoting hydration of the oceanic plate. However, deep penetration of seawater along bend-faults remains controversial because fluids that have percolated down into the mantle are difficult to detect. This report presents anomalously high helium isotope (3He/4He) ratios in sediment pore water and seismic reflection data which suggest fluid infiltration into the upper mantle and subsequent outflow through bend-faults across the outer slope of the Japan trench. The 3He/4He and 4He/20Ne ratios at sites near-trench bend-faults, which are close to the isotopic ratios of bottom seawater, are almost constant with depth, supporting local seawater inflow. Our findings provide the first reported evidence for a potentially large-scale active hydrothermal circulation system through bend-faults across the Moho (crust-mantle boundary) in and out of the oceanic lithospheric mantle.

A transcrustal mush system has been recognized beneath Dominica (Lesser Antilles) with different magma ponding zones that generated a series of pumiceous eruptions from Morne Trois Pitons–Micotrin volcano. Here, the latest, large, pumiceous eruption (Grand Fond - 24 kyrs cal BP) and four, smaller, Plinian eruptions (18-9 kyrs cal BP) are investigated. Pre-eruptive magma dynamics within the mush are unraveled through orthopyroxene phenocrysts by combining a Crystal System Analysis approach (on unzoned and zoned orthopyroxenes) and timescale estimates derived by intracrystalline Fe-Mg interdiffusion modeling. Two magmatic environments are recognized in the mush and have mixed, more or less vigorously, before the successive eruptions. Few interactions between the two magmas began 15-34 years prior to the small Plinian eruptions, but the sustained mixing occurred in the last 2 years. This contrasts with longer timescales (2-80 years) obtained for the larger eruption of Grand Fond with magmas stored deeper. These magma mixing timescales have significant implications for volcanic risk mitigation, with a growing reactivation signal that could be registered at the surface few years prior to the eruptions.

This paper resolve the salinity-dependent interactions of polar components of crude oil at calcite-brine interface in atomic resolution. Molecular dynamics simulations carried out on the present study showed that ordered water monolayers develop immediate to a calcite substrate in contact with a saline solution. Carboxylic compounds, herein represented by benzoic acid (BA), penetrate into those hydration layers and directly linking to the calcite surface. Through a mechanism termed screening effect, development of hydrogen bonding between –COOH functional groups of BA and carbonate groups is inhibited by formation of a positively-charged Na+ layer over CaCO3 surface. Contrary to the common perception, a sodium-depleted solution potentially intensifies surface adsorption of polar hydrocarbons onto carbonate substrates; thus, shifting wetting characteristic to hydrophobic condition. In the context of enhanced oil recovery, an ion-engineered waterflooding would be more effective than injecting a solely diluted saltwater.

Light-blue barite from Jebel Ouichane in Morocco forms blade-like tabular crystals (up to ca. 10 cm) with superb transparency and lustre and represents one of the most spectacular gem-quality worldwide. The barite is hosted by iron-ore-bearing skarns, developed within Jurassic-Cretaceous limestones, and occurs in close spatial association with calcite. The crystals have their cores enriched in Sr and contain abundant monophase (liquid) fluid inclusions of primary and pseudosecondary origin. The barite probably precipitated slowly at a relatively low supersaturation and under the control of a surface reaction precipitation mechanism. However, there were some episodes during its formation with a fast growth rate and the coupled dissolution and recrystallization processes. A combination of fluid inclusion data and stable δ18O value for barite (+ 6.71‰ VSMOW) suggests that low-salinity barite-forming solutions resulted from the mixing of strongly-diluted meteoric waters (enriched in light oxygen isotope) with magmatic-hydrothermal fluids under low-temperature conditions (42- derived from magmatic-hydrothermal SO2 reacted with sulphates that originate from the oxidation of H2S under near-surface conditions.