Guangzhou Marine Geological Survey

Clayey silt reservoirs bearing natural gas hydrates (NGH) are considered to be the hydrate-bearing reservoirs that boast the highest reserves but tend to be the most difficult to exploit. They are proved to be exploitable by the first NGH production test conducted in the South China Sea in 2017. Based on the understanding of the first production test, the China Geological Survey determined the optimal target NGH reservoirs for production test and conducted a detailed assessment, numerical and experimental simulation, and onshore testing of the reservoirs. After that, it conducted the second offshore NGH production test in 1225 m deep Shenhu Area, South China Sea (also referred to as the second production test) from October 2019 to April 2020. During the second production test, a series of technical challenges of drilling horizontal wells in shallow soft strata in deep sea were met, including wellhead stability, directional drilling of a horizontal well, reservoir stimulation and sand control, and accurate depressurization. As a result, 30 days of continuous gas production was achieved, with a cumulative gas production of 86.14 ×10⁴ m³. Thus, the average daily gas production is 2.87 ×10⁴ m³, which is 5.57 times as much as that obtained in the first production test. Therefore, both the cumulative gas production and the daily gas production were highly improved compared to the first production test. As indicated by the monitoring results of the second production test, there was no anomaly in methane content in the seafloor, seawater, and atmosphere throughout the whole production test. This successful production test further indicates that safe and effective NGH exploitation is feasible in clayey silt NGH reservoirs. The industrialization of hydrates consists of five stages in general, namely theoretical research and simulation experiments, exploratory production test, experimental production test, productive production test, and commercial production. The second production test serves as an important step from the exploratory production test to experimental production test.

Abstract Combined analyses of deep tow magnetic anomalies and International Ocean Discovery Program Expedition 349 cores show that initial seafloor spreading started around 33 Ma in the northeastern South China Sea (SCS), but varied slightly by 1–2 Myr along the northern continent‐ocean boundary (COB). A southward ridge jump of ∼20 km occurred around 23.6 Ma in the East Subbasin; this timing also slightly varied along the ridge and was coeval to the onset of seafloor spreading in the Southwest Subbasin, which propagated for about 400 km southwestward from ∼23.6 to ∼21.5 Ma. The terminal age of seafloor spreading is ∼15 Ma in the East Subbasin and ∼16 Ma in the Southwest Subbasin. The full spreading rate in the East Subbasin varied largely from ∼20 to ∼80 km/Myr, but mostly decreased with time except for the period between ∼26.0 Ma and the ridge jump (∼23.6 Ma), within which the rate was the fastest at ∼70 km/Myr on average. The spreading rates are not correlated, in most cases, to magnetic anomaly amplitudes that reflect basement magnetization contrasts. Shipboard magnetic measurements reveal at least one magnetic reversal in the top 100 m of basaltic layers, in addition to large vertical intensity variations. These complexities are caused by late‐stage lava flows that are magnetized in a different polarity from the primary basaltic layer emplaced during the main phase of crustal accretion. Deep tow magnetic modeling also reveals this smearing in basement magnetizations by incorporating a contamination coefficient of 0.5, which partly alleviates the problem of assuming a magnetic blocking model of constant thickness and uniform magnetization. The primary contribution to magnetic anomalies of the SCS is not in the top 100 m of the igneous basement.

Abstract Coring/logging data and physical property measurements from International Ocean Discovery Program Expedition 349 are integrated with, and correlated to, reflection seismic data to map seismic sequence boundaries and facies of the central basin and neighboring regions of the South China Sea. First‐order sequence boundaries are interpreted, which are Oligocene/Miocene, middle Miocene/late Miocene, Miocene/Pliocene, and Pliocene/Pleistocene boundaries. A characteristic early Pleistocene strong reflector is also identified, which marks the top of extensive carbonate‐rich deposition in the southern East and Southwest Subbasins. The fossil spreading ridge and the boundary between the East and Southwest Subbasins acted as major sedimentary barriers, across which seismic facies changes sharply and cannot be easily correlated. The sharp seismic facies change along the Miocene‐Pliocene boundary indicates that a dramatic regional tectonostratigraphic event occurred at about 5 Ma, coeval with the onsets of uplift of Taiwan and accelerated subsidence and transgression in the northern margin. The depocenter or the area of the highest sedimentation rate switched from the northern East Subbasin during the Miocene to the Southwest Subbasin and the area close to the fossil ridge in the southern East Subbasin in the Pleistocene. The most active faulting and vertical uplifting now occur in the southern East Subbasin, caused most likely by the active and fastest subduction/obduction in the southern segment of the Manila Trench and the collision between the northeast Palawan and the Luzon arc. Timing of magmatic intrusions and seamounts constrained by seismic stratigraphy in the central basin varies and does not show temporal pulsing in their activities.

Deep-sea sediments contain high concentrations of rare earth element (REE) which have been regarded as a huge potential resource. Understanding the marine REE cycle is important to reveal the mechanism of REE enrichment. In order to determine the geochemistry characteristics and migration processes of REE, seawater, porewater and sediment samples were systematically collected from the western Pacific for REE analysis. The results show a relatively flat REE pattern and the HREE (Heavy REE) enrichment in surface and deep seawater respectively. The HREE enrichment distribution patterns, low concentrations of Mn and Fe and negative Ce anomaly occur in the porewater, and high Mn/Al ratios and low U concentrations were observed in sediment, indicating oxic condition. LREE (Light REE) and MREE (Middle REE) enrichment in upper layer and depletion of MREE in deeper layer were shown in porewater profile. This study suggests that porewater flux in the western Pacific basin is a minor source of REEs to seawater, and abundant REEs are enriched in sediments, which is mainly caused by the extensive oxic condition, low sedimentation rate and strong adsorption capacity of sediments. Hence, the removal of REEs of porewater may result in widespread REE-rich sediments in the western Pacific basin.

[1] Gas hydrate saturations were estimated using five different methods in silt and silty clay foraminiferous sediments from drill hole SH2 in the South China Sea. Gas hydrate saturations derived from observed pore water chloride values in core samples range from 10 to 45% of the pore space at 190–221 m below seafloor (mbsf). Gas hydrate saturations estimated from resistivity (Rt) using wireline logging results are similar and range from 10 to 40.5% in the pore space. Gas hydrate saturations were also estimated by P wave velocity obtained during wireline logging by using a simplified three-phase equation (STPE) and effective medium theory (EMT) models. Gas hydrate saturations obtained from the STPE velocity model (41.0% maximum) are slightly higher than those calculated with the EMT velocity model (38.5% maximum). Methane analysis from a 69 cm long depressurized core from the hydrate-bearing sediment zone indicates that gas hydrate saturation is about 27.08% of the pore space at 197.5 mbsf. Results from the five methods show similar values and nearly identical trends in gas hydrate saturations above the base of the gas hydrate stability zone at depths of 190 to 221 mbsf. Gas hydrate occurs within units of clayey slit and silt containing abundant calcareous nannofossils and foraminifer, which increase the porosities of the fine-grained sediments and provide space for enhanced gas hydrate formation. In addition, gas chimneys, faults, and fractures identified from three-dimensional (3-D) and high-resolution two-dimensional (2-D) seismic data provide pathways for fluids migrating into the gas hydrate stability zone which transport methane for the formation of gas hydrate. Sedimentation and local canyon migration may contribute to higher gas hydrate saturations near the base of the stability zone.

The drilling recovered high-concentration methane hydrates (maximum 26–48%) in a disseminated form in silty clay sediments in Shenhu area of Pearl River Mouth Basin, South China Sea. Combining the geochemical data, the gas hydrate-bearing sediments are 10 m to 43 m in thickness and located just above the base of the gas hydrate stability zone. The methane content is 96.10–99.91% with small amount of ethane and propane. The baseline chlorinity of pore waters shows 10% lower than that of shallow sediments below and inside the gas hydrate zone. The methane/ethane ratios are higher than 1000 above the gas hydrate zone and less than 1000 at the interval of gas hydrate zone. The depth of sulphate methane interface varies from site to site as 17 to 27 mbsf. These results show that the methane of gas hydrate was mainly originated from microbial activity and the upward methane flux is minor. This is evidenced by the <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M1"><mml:mi>δ</mml:mi><mml:mtext>C</mml:mtext><mml:mn>13</mml:mn><mml:mrow><mml:mtext>CH4</mml:mtext></mml:mrow></mml:math> values of headspace gases from the gravity piston cores and released gases from pressure cores, which range from −74.3‰ PDB to −46.2‰ PDB, with the majority less than −55%‰ PDB. The hydrate deposit is a distributed gas hydrate system in Shenhu area.

Based on the comprehensive interpretation of cores, loggings and 2D/3D seismic data of Shenhu GMGS3 drilling area in the northern South China Sea, the distribution characteristics, differential accumulation mechanism and reservoir forming mechanism of diffusion type natural gas hydrate with high saturation discovered from clayey silt reservoirs were investigated. The following findings are reached through the research: (1) Gas hydrate with high saturation often displays high resistivity, low interval transit time, and strong bottom-simulating reflectors (BSRs), and accompanies with fluid seepage phenomena beneath BSRs, such as mud diapiric structure and gas chimney. (2) The gas hydrate reservoirs are dominated by fine grained clayey silt sediments, and the reservoirs have higher porosity and permeability in local parts. (3) The gas hydrate is largely type I, whereas type II gas hydrate may exist below the type I gas hydrate. (4) The gas sources are mixed microbial and thermogenic gases, and the thermogenic gas originated from the deep formation in the center of Baiyun Sag migrated into shallow strata through faults, mud diapirs and gas chimneys, then was mixed with microbial gas in situ and continued to migrate until they accumulated in the temperature and pressure stability zone and formed diffusion type gas hydrate with high saturation finally. (5) The fluid migration system influenced and controlled the differential distribution of gas hydrate with high saturation.

In traditional soil heavy metal (HM) pollution assessment, spatial interpolation analysis is often carried out on the limited sampling points in the study area to get the overall status of heavy metal pollution. Unfortunately, in many machine learning spatial information enhancement algorithms, the additional spatial information introduced fails to reflect the hierarchical heterogeneity of the study area. Therefore, we designed hierarchical regionalization labels based on three interpolation techniques (inverse distance weight, ordinary kriging, and trend surface interpolation) as new spatial covariates for a machine learning (ML) model. It was demonstrated that regional spatial information improved the prediction performance of the model ( R 2 > 0.7). On the basis of the prediction results, the status of HM pollution in the Pearl River Delta (PRD) region was evaluated: cadmium (Cd) and copper (Cu) were the most serious pollutants in the PRD (the point overstandard rates are 18.77% and 12.95%, respectively). The analysis of index importance and bivariate local indicators of spatial association (LISA) shows that the key factors affecting the spatial distribution of heavy metals are geographical and climatic conditions [namely, altitude, humidity index, and normalized vegetation difference index (NDVI)] and some industrial activities (such as metal processing, printing and dyeing, and electronics industry). This study develops a novel approach to improve existing spatial interpolation techniques, which will enable more precise and scientific soil environmental management.

© The Author(s), 2019. This article is distributed under the terms of the Creative Commons Attribution License. The definitive version was published in Sun, Z., Lin, J., Qiu, N., Jian, Z., Wang, P., Pang, X., Zheng, J., &amp; Zhu, B. The role of magmatism in the thinning and breakup of the South China Sea continental margin: Special Topic: the South China Sea Ocean Drilling. National Science Review, 6(5), (2019): 871-876, doi:10.1093/nsr/nwz116.

During the Holocene there has been a gradual increase in the influence of humans on Earth systems. High-resolution sedimentary records can help us to assess how erosion and weathering have evolved in response to recent climatic and anthropogenic disturbances. Here we present data from a highresolution (75 cm/k.y.) sedimentary archive from the South China Sea. Provenance data indicate that the sediment was derived from the Red River, and can be used to reconstruct the erosion and/or weathering history in this river basin. Accelerator mass spectrometry 14C dating provides direct age control and reveals coherent variations in clay mineralogy, geochemistry, and terrigenous flux, indicative of strong chemical weathering and physical erosion during the mid-Holocene warm period (6400-4000 cal [calibrated] yr B.P.), followed by weakening from ca. 4000-1800 cal yr B.P., and renewed intensification since 1800 cal yr B.P.. Comparison with climatic records from China indicates that precipitation and temperature controlled both physical erosion and chemical weathering intensity before 1800 cal yr B.P.. However, weathering proxies in the offshore sediment indicate recent increased soil erosion. We suggest that enhanced human activity (deforestation, cultivation, and mining) since the end of the Chinese Han Dynasty (220 CE) has overwhelmed the natural climatic controls on erosion in the Red River.

Shenhu Area is located in the Baiyun Sag of Pearl River Mouth Basin, which is on the northern continental slope of the South China Sea. Gas hydrates in this area have been intensively investigated, achieving a wide coverage of the three-dimensional seismic survey, a large number of boreholes, and detailed data of the seismic survey, logging, and core analysis. In the beginning of 2020, China has successfully conducted the second offshore production test of gas hydrates in this area. In this paper, studies were made on the structure of the hydrate system for the production test, based on detailed logging data and core analysis of this area. As to the results of nuclear magnetic resonance (NMR) logging and sonic logging of Well GMGS6-SH02 drilled during the GMGS6 Expedition, the hydrate system on which the production well located can be divided into three layers: (1) 207.8–253.4 mbsf, 45.6 m thick, gas hydrate layer, with gas hydrate saturation of 0–54.5% (31% av.); (2) 253.4–278 mbsf, 24.6 m thick, mixing layer consisting of gas hydrates, free gas, and water, with gas hydrate saturation of 0–22% (10% av.) and free gas saturation of 0–32% (13% av.); (3) 278–297 mbsf, 19 m thick, with free gas saturation of less than 7%. Moreover, the pore water freshening identified in the sediment cores, taken from the depth below the theoretically calculated base of methane hydrate stability zone, indicates the occurrence of gas hydrate. All these data reveal that gas hydrates, free gas, and water coexist in the mixing layer from different aspects.

Previous research studies have demonstrated that the relationship between remote sensing-derived parameters and aboveground biomass (AGB) could vary across different species types. However, there are few studies that calibrate reliable statistical models for mangrove AGB. This study quantifies the differences of accuracy in AGB estimation between the results obtained with and without the consideration of species types using Worldview-2 images and field surveys. A Back Propagation Artificial Neural Network (BP ANN) based model is developed for the accurate estimation of uneven-aged and dense mangrove forest biomass. The contributions of the input variables are further quantified using a “Weights” method based on BP ANN model. Two types of mangrove species, Sonneratia apetala (S. apetala) and Kandelia candel (K. candel), are examined in this study. Results show that the species type information is the most important variable for AGB estimation, and the red edge band and the associated vegetation indices from WorldView-2 images are more sensitive to mangrove AGB than other bands and vegetation indices. The RMSE of biomass estimation at the incorporation of species as a dummy variable is 19.17% lower than that of the mixed species level. The results demonstrate that species type information obtained from the WorldView-2 images can significantly improve of the accuracy of the biomass estimation.

An analysis of the water level and current data taken in Qiongzhou Strait in the South China Sea (SCS) over the last 37 years (1963 to 1999) was made to examine the characteristics of tidal waves and residual flow through the strait and their roles in the seasonal variation of the SCS circulation. The observations reveal that Qiongzhou Strait is an area where opposing tidal waves interact and a source of water transport to the Gulf of Beibu (Gulf of Tonkin), SCS. A year-round westward mean flow with a maximum speed of 10-40 cm s 1 is found in Qiongzhou Strait. This accounts for water transport of 0.2-0.4 Sv and 0.1-0.2 Sv into the Gulf of Beibu in winter-spring and summer-autumn, respectively. The outflow from Qiongzhou Strait may cause up to 44% of the gulf water to be refreshed each season, suggesting that it has a significant impact on the seasonal circulation in the Gulf of Beibu. This finding is in contrast to our current understanding that the seasonal circulation patterns in the South China Sea are primarily driven by seasonal winds. Several numerical experiments were conducted to examine the physical mechanisms responsible for the formation of the westward mean flow in Qiongzhou Strait. The model provides a reasonable simulation of semidiurnal and diurnal tidal waves in the strait and the predicted residual flow generally agrees with the observed mean flow. An analysis of the momentum equations indicates that the strong westward flow is driven mainly by tidal rectification over variable bottom topography. Both observations and modeling suggest that the coastal physical processes associated with tidal rectification and buoyancy input must be taken into account when the mass balance of the SCS circulation is investigated, especially for the regional circulation in the Gulf of Beibu.

To accurately estimate leaf area index (LAI) in mangrove areas, the selection of appropriate models and predictor variables is critical. However, there is a major challenge in quantifying and mapping LAI using multi-spectral sensors due to the saturation effects of traditional vegetation indices (VIs) for mangrove forests. WorldView-2 (WV2) imagery has proven to be effective to estimate LAI of grasslands and forests, but the sensitivity of its vegetation indices (VIs) has been uncertain for mangrove forests. Furthermore, the single model may exhibit certain randomness and instability in model calibration and estimation accuracy. Therefore, this study aims to explore the sensitivity of WV2 VIs for estimating mangrove LAI by comparing artificial neural network regression (ANNR), support vector regression (SVR) and random forest regression (RFR). The results suggest that the RFR algorithm yields the best results (RMSE = 0.45, 14.55% of the average LAI), followed by ANNR (RMSE = 0.49, 16.04% of the average LAI), and then SVR (RMSE = 0.51, 16.56% of the average LAI) algorithms using 5-fold cross validation (CV) using all VIs. Quantification of the variable importance shows that the VIs derived from the red-edge band consistently remain the most important contributor to LAI estimation. When the red-edge band-derived VIs are removed from the models, estimation accuracies measured in relative RMSE (RMSEr) decrease by 3.79%, 2.70% and 4.47% for ANNR, SVR and RFR models respectively. VIs derived from red-edge band also yield better accuracy compared with other traditional bands of WV2, such as near-infrared-1 and near-infrared-2 band. Furthermore, the estimated LAI values vary significantly across different mangrove species. The study demonstrates the utility of VIs of WV2 imagery and the selected machine-learning algorithms in developing LAI models in mangrove forests. The results indicate that the red-edge band of WV2 imagery can help alleviate the saturation problem and improve the accuracy of LAI estimation in a mangrove area.

Large-scale step-like features within the South Taiwan Shoal and West Penghu submarine canyons on the northeastern continental slope of the South China Sea are investigated by integrating high-resolution multibeam bathymetric data and multichannel seismic profiles. These step-like features, ranging from 1.2 to 10.0 km in wavelength and 5.4–80.9 m in wave height, are mostly interpreted as cyclic steps formed by turbidity currents flowing through the canyons, based on their characteristic step-like morphology, in-train alignment, large wavelengths and aspect ratios (ratio of wavelength to wave height), and typical upstream-sloping backset bedding, among others. A train of 19 continuous steps delineated along the thalweg of the South Taiwan Shoal canyon measures up to 100 km and may be the longest ever reported. Nine short trains of scours identified on a terrace of the South Taiwan Shoal canyon are oriented parallel to the distributaries draining over the terrace and roughly perpendicular to the main canyon thalweg, indicating a complicated flow pattern within the canyon valley. Two trains of scours separated by an intracanyon high in the steeper middle reach of the West Penghu canyon are interpreted as transitional bed forms between antidunes and cyclic steps, which develop downstream into a train of five net-depositional cyclic steps with typical backset bedding in the gentler-sloping lower reach of the canyon. Average slope gradients for the canyon reaches with cyclic steps range from 0.26° to 1.24°. Along each thalweg step train, a slope break is identified to separate the net-erosional cyclic steps in the steeper upstream segment from the net-depositional ones in the gentler downstream segment. Rough estimations indicate that the paleoflows are 100 to 300 m thick with maximum velocities of up to 10 m s–1. The estimated flow depths match well with those inferred from geomorphologic analysis. Estimated paleodischarges of ?7–23 × 105 m3 s–1 are equivalent to ten times the discharge of the modern Amazon River.

Endosymbiosis with chemosynthetic bacteria has enabled many deep-sea invertebrates to thrive at hydrothermal vents and cold seeps, but most previous studies on this mutualism have focused on the bacteria only. Vesicomyid clams dominate global deep-sea chemosynthesis-based ecosystems. They differ from most deep-sea symbiotic animals in passing their symbionts from parent to offspring, enabling intricate coevolution between the host and the symbiont. Here, we sequenced the genomes of the clam Archivesica marissinica (Bivalvia: Vesicomyidae) and its bacterial symbiont to understand the genomic/metabolic integration behind this symbiosis. At 1.52 Gb, the clam genome encodes 28 genes horizontally transferred from bacteria, a large number of pseudogenes and transposable elements whose massive expansion corresponded to the timing of the rise and subsequent divergence of symbiont-bearing vesicomyids. The genome exhibits gene family expansion in cellular processes that likely facilitate chemoautotrophy, including gas delivery to support energy and carbon production, metabolite exchange with the symbiont, and regulation of the bacteriocyte population. Contraction in cellulase genes is likely adaptive to the shift from phytoplankton-derived to bacteria-based food. It also shows contraction in bacterial recognition gene families, indicative of suppressed immune response to the endosymbiont. The gammaproteobacterium endosymbiont has a reduced genome of 1.03 Mb but retains complete pathways for sulfur oxidation, carbon fixation, and biosynthesis of 20 common amino acids, indicating the host's high dependence on the symbiont for nutrition. Overall, the host-symbiont genomes show not only tight metabolic complementarity but also distinct signatures of coevolution allowing the vesicomyids to thrive in chemosynthesis-based ecosystems.

In this article, we present an adaptive reinforcement learning optimal tracking control (RLOTC) algorithm for an underactuated surface vessel subject to modeling uncertainties and time-varying external disturbances. By integrating backstepping technique with the optimized control design, we show that the desired optimal tracking performance of vessel control is guaranteed due to the fact that the virtual and actual control inputs are designed as optimized solutions of every subsystem. To enhance the robustness of vessel control systems, we employ neural network (NN) approximators to approximate uncertain vessel dynamics and present adaptive control technique to estimate the upper boundedness of external disturbances. Under the reinforcement learning framework, we construct actor-critic networks to solve the Hamilton-Jacobi-Bellman equations corresponding to subsystems of surface vessel to achieve the optimized control. The optimized control algorithm can synchronously train the adaptive parameters not only for actor-critic networks but also for NN approximators and adaptive control. By Lyapunov stability theorem, we show that the RLOTC algorithm can ensure the semiglobal uniform ultimate boundedness of the closed-loop systems. Compared with the existing reinforcement learning control results, the presented RLOTC algorithm can compensate for uncertain vessel dynamics and unknown disturbances, and obtain the optimized control performance by considering optimization in every backstepping design. Simulation studies on an underactuated surface vessel are given to illustrate the effectiveness of the RLOTC algorithm.

Cold seeps are widespread chemosynthetic ecosystems in the deep-sea environment, and cold seep microbial communities of the South China Sea are poorly constrained. Here we report on the archaeal communities, particularly those involved in methane metabolization, in sediments of a newly discovered cold seep (named 'Haima') on the northwest slope of the South China Sea. Archaeal diversity, abundance and distribution were investigated in two piston cores collected from a seep area (QDN-14B) and a non-seep control site (QDN-31B). Geochemical investigation of the QDN-14B core identified an estimated sulfate-methane transition zone (Estimated SMTZ) at 300-400 cm below sea floor (cmbsf), where a high abundance of anaerobic methane-oxidizing archaea (ANME) occurred, as revealed by analysis of the 16S rRNA gene and the gene (mcrA) encoding the α-subunit of the key enzyme methyl-coenzyme M reductase. ANME-2a/b was predominant in the upper and middle layers of the estimated SMTZ, whereas ANME-1b outcompeted ANME-2 in the sulfate-depleted bottom layers of the estimated SMTZ and the methanogenic zone. Fine-scale phylogenetic analysis further divided the ANME-1b group into three subgroups with different distribution patterns: ANME-1bI, ANME-1bII and ANME-1bIII. Multivariate analyses indicated that dissolved inorganic carbon and sulfate may be important factors controlling the composition of the methane-metabolizing community. Our study on ANME niche separation and interactions with other archaeal groups improves our understanding of the metabolic diversity and flexibility of ANME, and the findings further suggest that ANME subgroups may have evolved diversified/specified metabolic capabilities other than syntrophic anaerobic oxidation of methane coupled with sulfate reduction in marine sediments.

Gas hydrates distribution in the Shenhu area, northern South China Sea: comparisons between the eight drilling sites with gashydrate petroleum systemThe results of the first marine gas hydrate drilling expedition of Guangzhou Marine Geological Survey (GMGS-1) in northern continental slope of the South China Sea revealed a variable distribution of gas hydrates in the Shenhu area.In this study, comparisons between the eight sites with gas-hydrate petroleum system were used to analyze and re-examine hydrate potential.In the Shenhu gas hydrate drilling area, all the sites were located in a suitable low-temperature, high-pressure environment.Biogenic and thermogenic gases contributed to the formation of hydrates.Gas chimneys and some small-scale faults (or micro-scale fractures) compose the migration pathways for gas-bearing fluids.Between these sites, there are three key differences: the seafloor temperatures and pressures; geothermal gradient and sedimentary conditions.Variations of seafloor temperatures and pressures related to water depths and geothermal gradient would lead to changes in the thickness of gas hydrate stability zones.Although the lithology and grain size of the sediments were similar, two distinct sedimentary units were identified for the first time through seismic interpretation, analysis of deep-water sedimentary processes, and the Cm pattern (plotted one-percentile and median values from grain-size analyses), implying the heterogeneous sedimentary conditions above Bottom Simulating Reflectors (BSRs).Based on the analyses of forming mechanisms and sedimentary processes, these two fine-grained sedimentary units have different physical properties.Fine-grained turbidites (Unit I) with thin-bedded chaotic reflectors at the bottom acted as the host rocks for hydrates; whereas, finegrained sediments related to soft-sediment deformation (Unit II) characterized by thick continuous reflectors at the top would serve as regional homogeneous caprocks.Low-flux methane that migrated upwards along chimneys could be enriched preferentially in fine-grained turbidites, resulting in the formation of hydrates within Unit I.

The rare earth elements and yttrium (REY) in bioapatite from deep-sea sediments are potential proxies for reconstructing paleoenvironmental conditions. However, the REY enrichment mechanism and the reliability of this tracer remain elusive because of the lack of key information from ambient pore water. Here, we report high-resolution geochemical data for pore water, bottom water, and bioapatite from deep-sea sites in the western Pacific. Our results reveal that the benthic flux of REY from the deep sea is less substantial than from the shallow marine realm, resulting in REY-rich sediment. The depth distribution of REY in pore water is opposite to that of bioapatite, and REY patterns and neodymium isotopic compositions are not uniformly distributed within bioapatite. These results indicate alteration of REY and neodymium isotopic compositions during early diagenesis. Therefore, we infer that REY from bioapatite are not robust recorders of the deep marine environment through Earth's history.