Constraints on Characteristics and Distribution of Gas Hydrate and Free Gas Using Broad-Band Processing of Three-Dimensional Seismic Data

Gas hydrate drilling expeditions in the Pearl River Mouth Basin,South China Sea,have identified concentrated gas hydrates with variable thickness.Moreover,free gas and the coexistence of gas hydrate and free gas have been confirmed by logging,coring,and production tests in the foraminifera-rich silty sediments with complex bottom-simulating reflectors(BSRs).The broad-band processing is conducted on conventional three-dimensional(3D)seismic data to improve the image and detection accuracy of gas hydratebearing layers and delineate the saturation and thickness of gas hydrate-and free gas-bearing sediments.Several geophysical attributes extracted along the base of the gas hydrate stability zone are used to demonstrate the variable distribution and the controlling factors for the differential enrichment of gas hydrate.The inverted gas hydrate saturation at the production zone is over 40% with a thickness of 90 m,showing the interbedded distribution with different boundaries between gas hydrate-and free gas-bearing layers.However,the gas hydrate saturation value at the adjacent canyon is 70%,with 30-m-thick patches and linear features.The lithological and fault controls on gas hydrate and free gas distributions are demonstrated by tracing each gas hydrate-bearing layer.Moreover,the BSR depths based on broad-band reprocessed 3D seismic data not only exhibit variations due to small-scale topographic changes caused by seafloor sedimentation and erosion but also show the upward shift of BSR and the blocky distribution of the coexistence of gas hydrate and free gas in the Pearl River Mouth Basin.

Geochronology and Geochemistry of the Xingxingxia Triassic A-type Granites in Central Tianshan,NW China:Petrogenesis and Tectonic Implications

The Triassic granitoids in Central Tianshan play a key role in determining the petrogenesis and tectonic evolution on the southern margin of the Central Asian orogenic belt.In this study,we present SHRIMP zircon U-Pb ages,Hf isotopic and geochemical data on the Xingxingxia biotite granite,amazonite granite and granitic pegmatite in Central Tianshan,NW China.Zircon U-Pb dating yielded formation ages of 242 Ma for the biotite granite and 240 Ma for the amazonite granite.These granitoid rocks have high K_(2)O with low MgO and CaO contents.They are enriched in Nb,Ta,Hf and Y,while being depleted in Ba and Sr,showing flat HREE patterns and negative Eu anomalies.They have typical A-type granite geochemical signatures with high Ga/A_(1)(8–13)and TFeO/(TFeO+MgO)ratios,showing an A_(2) affinity for biotite granite and an A_(1) affinity for amazonite granite and granitic pegmatite.Zircon ε_(Hf)(t)values of the granitoids are 0.45–2.66,with Hf model ages of 0.99–1.17 Ga.This suggests that these A-type granites originated from partial melting of the lower crust.We propose that Xingxingxia Triassic A-type granites formed under lithospheric extension from post-orogenic to anorogenic intraplate settings and NE-trending regional strike-slip fault-controlled magma emplacement in the upper crust.

Elemental Compositions of Clinopyroxenes from Southern Okinawa Trough Basalt:Implications for Magmatism and Structural Environment in an Initial Back-Arc Basin

The major and trace elemental compositions of clinopyroxene from basalt were used to characterize the nature of the primitive magma and structural environment beneath the southern Okinawa Trough(SOT),which is an initial back-arc basin at a continental margin.The clinopyroxenes in the basalt were augite with variable Mg^(#)contents(73.37-78.22).The regular variations in major oxide contents(i.e.,CaO,FeO,TiO_(2),and Cr)with decreasing Mg#implied that the clinopyroxenes evolved from being enriched in Mg,Ca,and Cr to being enriched in Fe and Ti.The clinopyroxenes had relatively low rare earth element concentrations(7.51×10^(-6) to 12.68×10^(-6))and negative Eu anomalies(δEu=0.67-0.95).The Kd_(cpx) values of clinopyroxenes(0.2-0.26),which were used to examine whether the clinopyroxene was equilibrated with its host basalt,demonstrate that these clinopyroxene phenocrysts were not captured crystals but were instead produced by crystallization differentiation of the magma.The calculated clinopyroxene crystallization temperatures showed a narrow range of 990-1061℃,and their crystallization pressures ranged from 2.0 to 3.2 kbar.The geochemistry features of these clinopyroxenes indicated that the parent magma belonged to the subalkaline tholeiitic magma series and suggested that the magma experienced crystallization differentiation of olivine,plagioclase,and clinopyroxene,where the crystallization of plagioclase occurred earlier than that of clinopyroxene.Combined with geophysical data,this research on primitive magma and its crystallization differentiation from clinopyroxene indicates that the SOT is in the stage of‘seafloor spreading’and that basaltic rocks produced from tholeiitic magma represent the generation of oceanic crust.

Correlation Study on Grain Size Characteristics and Geotechnical Properties of Surface Sediments in Qingdao Offshore Area

The comprehension of sediment grain size parameters and the corresponding sedimentary environment holds paramount importance in elucidating the engineering geological attributes of the subaqueous seabed.This study delineated the sedimentary environment zoning in the northern sea area of Qingdao through cluster analysis of grain size parameters derived from 123 surface sediment samples.The study analyzed the correlation between sediment geotechnical indices and grain size parameters across diverse sedimentary environments.A correlation equation was established for samples exhibiting a strong correlation.The study found four distinct sedimentary environments in the study area:coastal,transitional,shallow sea,and residual.Within the same sedimentary environment,the average grain size and sorting coefficient exhibit significant correlations with geotechnical indices such as water content,density,shear strength,plastic limit,liquid limit,and plastic index.However,notable disparities in the correlation between grain size parameters and geotechnical indices emerge across different sedimentary environments.

Effect of particle composition and consolidation degree on the wave-induced liquefaction of soil beds

The wave-induced liquefaction of seabed is responsible for causing damage to marine structures.Particle composition and consolidation degree are the key factors affecting the pore water pressure response and liquefaction behavior of the seabed under wave action.The present study conducted wave flume experiments on silt and silty fine sand beds with varying particle compositions.Furthermore,a comprehensive analysis of the differences and underlying reasons for liquefaction behavior in two different types of soil was conducted from both macroscopic and microscopic perspectives.The experimental results indicate that the silt bed necessitates a lower wave load intensity to attain the liquefaction state in comparison to the silty fine sand bed.Additionally,the duration and development depth of liquefaction are greater in the silt bed.The dissimilarity in liquefaction behavior between the two types of soil can be attributed to the variation in their permeability and plastic deformation capacity.The permeability coefficient and compression modulus of silt are lower than those of silty fine sand.Consequently,silt is more prone to the accumulation of pore pressure and subsequent liquefaction under external loading.Prior research has demonstrated that silt beds with varying consolidation degrees exhibit distinct initial failure modes.Specifically,a dense bed undergoes shear failure,whereas a loose bed experiences initial liquefaction failure.This study utilized discrete element simulation to examine the microscopic mechanisms that underlie this phenomenon.

Pre-Drilling Prediction Techniques on the High-Temperature High-Pressure Hydrocarbon Reservoirs Offshore Hainan Island,China

Decreasing the risks and geohazards associated with drilling engineering in high-temperature high-pressure(HTHP) geologic settings begins with the implementation of pre-drilling prediction techniques(PPTs). To improve the accuracy of geopressure prediction in HTHP hydrocarbon reservoirs offshore Hainan Island, we made a comprehensive summary of current PPTs to identify existing problems and challenges by analyzing the global distribution of HTHP hydrocarbon reservoirs, the research status of PPTs, and the geologic setting and its HTHP formation mechanism. Our research results indicate that the HTHP formation mechanism in the study area is caused by multiple factors, including rapid loading, diapir intrusions, hydrocarbon generation, and the thermal expansion of pore fluids. Due to this multi-factor interaction, a cloud of HTHP hydrocarbon reservoirs has developed in the Ying-Qiong Basin, but only traditional PPTs have been implemented, based on the assumption of conditions that do not conform to the actual geologic environment, e.g., Bellotti's law and Eaton's law. In this paper, we focus on these issues, identify some challenges and solutions, and call for further PPT research to address the drawbacks of previous works and meet the challenges associated with the deepwater technology gap. In this way, we hope to contribute to the improved accuracy of geopressure prediction prior to drilling and provide support for future HTHP drilling offshore Hainan Island.

Grain-Size Distribution of Surface Sediments in the Bohai Sea and the Northern Yellow Sea: Sediment Supply and Hydrodynamics

The grain-size distribution of surface sediments in the Bohai Sea(BS) and the northern Yellow Sea(NYS), and its relationship with sediment supply and hydrodynamic environment were investigated based on grain-size compositions of surface sediments and modern sedimentation rates. The results showed that the surface sediments in the BS and the NYS were primarily composed of silty sand and clayey silt with a dominant size of silt. In addition, the Yellow River delivered high amount of water and sediments to the BS, and they are dominated in surface sediments(mainly silt) in the Bohai Bay, the Yellow River mouth, the center of the BS, and the north coast of Shandong Peninsula. The coarse-grained sediments were mainly deposited at the river mouth due to the estuarine filtration and physical sorting. Meanwhile, there was a significant relationship among the modern sedimentation rate, the surface sediment grain size distribution and sediment transport pattern. The areas with coarser surface sediments generally corresponded low sedimentation rates because of strong erosion;whereas the sedimentation rate was relatively high at the place that the surface sediments were fine-grained. Furthermore, the grain-size trend analysis showed that the areas with fine-grained surface sediments such as the mud area in the central BS and the upper Liaodong Bay were the convergent centers of surface sediments, except for the Bohai Bay and the subaqueous Yellow River Delta where offshore sediment transport was evident.

Sensitivity Analysis of P-Waves and S-Waves to Gas Hydrate in the Shenhu Area Using OBS

Compared to towed streamers, ocean-bottom seismometers(OBS) obtain both S-wave data and richer wavefield information. In this paper, the induced polarization method is used to conduct wavefield separation on OBS data obtained from the Shenhu area in the South China Sea. A comparison of the changes in P-and S-waves, and a comprehensive analysis of geological factors within the area, enable analysis and description of the occurrence of natural gas hydrate in the study area. Results show an increase in P-wave velocity when natural gas hydrate exists in the formation, whereas the S-wave velocity remains almost constant, as S-waves can only propagate through the rock skeleton. Therefore, the bottom-simulating reflection(BSR) response of the P-wave is better than that of the S-wave in the frequency analysis profile. In a wide-angle section, the refractive wave of the hydrate layer is evident when using P-wave components but identification is difficult with S-wave components. This velocity model illustrates the sensitivity of Pand S-wave components to gas hydrate. The use of this polarization method and results of analysis provide technical and theoretical support for research on hydrate deposits and other geological features in the Shenhu area.

Heavy metals in Changjiang estuarine and offshore sediments:responding to human activities

The Changjiang （Yangtze） estuarine and offshore sediments were analyzed for total heavy metals concentrations and chemical fractions. Distributions of heavy metals show typical banded diffusion pattern, with high concentrations near the river mouth and following a decreasing trend in the offshore direction. According to chemical fractions, Fe/Mn oxide fraction is the major non-residual fraction in the Changjiang estuarine and offshore sediments, Higher percentage of non-residual fraction of Pb implies that, the industrial contaminations transported via the atmosphere and river input, may affect the non-residual fraction of heavy metals. Over past fifteen years, the concentration of Pb normalizing to A1 presents significant increasing trend, corresponding to the effect of human activities. By comparison of heavy metals fractions in 2003 to 2006, it has been realized that increasing water and sediment may cause a higher percentage non-residual fraction of Cu in the southern part of offshore muddy sediments.

Evolution history and trend of the modern Huanghe River Delta

The evolvement history of modern Huanghe River Delta tidal flat, coastline and underwater terrain were studied based on the analysis of remote sensing images and water depth data. Based on the analysis of seafloor terrain evolution on different historical stages, a formula simulating the erosion and deposition evolvement model of subaqueous Huanghe River Delta slope was proposed, and the evolvement trend of the subaqueous delta terrain was predicted. The result shows that the equilibrium transition zone is near the water depth of 12 m with seabed erosion in shallower water and accumulation in deeper water during the first 150 a after the river channel was deserted. In the meantime, the underwater slope became gentler and the coastal erosion rate became slow gradually. Then, the subaqueous delta slope changed to up concave from upper convex, and the shape of subaqueous delta disappeared. The coast type changed to silt-mud coast about 100-150 a after the river course was deserted. The erosion depth in the foot of the seawall is calculated based on the formula.

Temporal-Spatial Variation of Surface Suspended Matter and Controlling Factors in the Inner Shelf of the East China Sea in Winter

Based on the observed surface suspended matter in the East China Sea in February 2007 and June 2015, an empirical model was established using L1 b's band 4 data to retrieve surface suspended matter from the Moderate Resolution Imagine Spectroradiometer Terra imagery. The squared correlation coefficient is 0.8358, and the root mean square error is 0.4285 mg L-1. The model reflects the distribution characteristics of surface suspended matter in the inner shelf of the East China Sea. In this paper, the satellite images of the study area were retrieved in January from 2001 to 2015, and the monthly distribution of surface suspended matter were obtained. The inter-annual distribution of the study area is similar, and the concentration of surface suspended matter is higher near the shore than offshore. A large amount of surface suspended matter is transported southeast under the influence of Zhejiang and Fujian coastal current and Taiwan warm current. Only a small amount of surface suspension can reach the Kuroshio area. The surface suspended matter concentration changes obviously near the estuary because of the effect of differences in the flux of the Yangtze River. Meanwhile, winter monsoon, temperature front, El Ni?o events, and other factors affect the distribution of surface suspended matter in 100 m isobath to coastal water but minimally influence the distribution in 100 m isobath to deep sea.

Effects of Uncertainties in the Position and Orientation of Both the Transmitter and Receivers on Marine Controlled-Source Electromagnetic Data

Simulation and interpretation of marine controlled-source electromagnetic(CSEM) data often approximate the transmitter source as an ideal horizontal electric dipole(HED) and assume that the receivers are located on a flat seabed.Actually,however,the transmitter dipole source will be rotated,tilted and deviated from the survey profile due to ocean currents.And free-fall receivers may be also rotated to some arbitrary horizontal orientation and located on sloping seafloor.In this paper,we investigate the effects of uncertainties in the transmitter tilt,transmitter rotation and transmitter deviation from the survey profile as well as in the receiver's location and orientation on marine CSEM data.The model study shows that the uncertainties of all position and orientation parameters of both the transmitter and receivers can propagate into observed data uncertainties,but to a different extent.In interpreting marine data,field data uncertainties caused by the position and orientation uncertainties of both the transmitter and receivers need to be taken into account.

Geochemical Characteristics of the Trace and Rare Earth Elements in Reef Carbonates from the Xisha Islands(South China Sea): Implications for Sediment Provenance and Paleoenvironment

Based on the concentrations of the trace elements,rare earth elements(REE),and Sr isotopic compositions in reef carbonates from the well‘Xike-1’reef core of the Xisha Islands,the constraints on sediment provenance and paleoenvironment were defined.Variations of the terrigenous input into the paleoseawater were traced in detail and the paleoenvironment and sediment provenance were further investigated.The results show that the HREE/LREE values in the reef carbonates are negatively associated with their Th and Al concentrations;however,their Al and Th concentrations show positive correlation.The lowest 87 Sr/86 Sr values in the reef carbonates generally coincide with the lowest values of Al,Th concentrations and the highest values of HREE/LREE.These data indicate that the HREE/LREE,Al concentrations,and Th concentrations of the reef carbonates are useful indexes for evaluating the influence of the terrigenous inputs on the seawater composition in the study area.From top to bottom,the changing process of the HREE/LREE values and Al,Th concentrations can be divided into 6 intervals;they are H1(0–89.30 m,about 0–0.11 Myr),L1(89.30–198.30 m,about 0.11–2.2 Myr),H2(198.30–374.95 m,about 2.2–5.3 Myr),D(374.95–758.40 m,about 5.3–13.6 Myr),L2(758.40–976.86 m,about 13.6–15.5 Myr),and H3(976.86–1200.00 m,about 15.5–21.5 Myr).Moreover,the changing trend of the HREE/LREE values coincides with that of the seawater δ^13C values recorded by benthonic foraminiferal skeletons from the drill core of ODP site 1148 in the South China Sea(SCS),but not with that of the seawaterδ18O values.The high uplifting rates of the Qinghai-Tibet Plateau coincide with the high Th and Al concentrations and the low HREE/LREE values in the reef carbonates.These data indicate that the main factors controlling the changes of terrigenous flux in the SCS are the tectonic activities associated with Qinghai-Tibet Plateau uplifting and the variations of uplifting rates rather than paleoclimatic changes.

Geochemical and Grain-Sized Implications for Provenance Variations of the Central Yellow Sea Muddy Area Since the Middle Holocene

Based on high-resolution analysis to a 280-cm long sediment core obtained from the muddy area in the central Yellow Sea, we examined the provenance of muddy sediments and discussed the changing marine sedimentary environment since the middle Holocene. The results indicated that fine-grained sediments in the muddy area were mainly derived from the Huanghe(Yellow River) and Changjiang(Yangtze River) with considerable stepwise variations during the past 6.6 kyr. The Yellow Sea Warm Current was initiated at 6 kyr when the sea level was high together with the enhanced East Asian Winter Monsoon. These in combination established the framework of shelf circulation in the Yellow Sea that began to trap the river-derived fine-grained sediments. From 4.9 kyr to 2.8 kyr, both the Kushiro Current and East Asian Monsoon were significantly weakened, reducing the delivery of Changjiang sediments to the muddy area. As a result, the sediments were mainly originated from the Huanghe. From 2.8 kyr to 1.5 kyr the continuously weakened East Asian Winter Monsoon and enhanced Yellow Sea Warm Current entrapped more fine-grain sediments. Whereas the enhanced East Asian Winter Monsoon and the human caused increase in sediment load of the Huanghe since 1.5 kyr, and direct delivery of Huanghe sediments to the Yellow Sea during 1128–1855 AD might dominated the sedimentation in the study area. The stepwise variations of the sediment provenance and composition of the Central Yellow Sea muddy sediments are of importance to understanding the formation of muddy deposit in the central Yellow Sea and the associated variations of marine environment since the middle Holocene.

Multi-beam and seismic investigations of the active Haima cold seeps,northwestern South China Sea

To confirm the seabed fluid flow at the Haima cold seeps,an integrated study of multi-beam and seismic data reveals the morphology and fate of four bubble plumes and investigates the detailed subsurface structure of the active seepage area.The shapes of bubble plumes are not constant and influenced by the northeastward bottom currents,but the water depth where these bubble plumes disappear(630–650 m below the sea level)(mbsl)is very close to the upper limit of the gas hydrate stability zone in the water column(620 m below the sea level),as calculated from the CTD data within the study area,supporting the“hydrate skin”hypothesis.Gas chimneys directly below the bottom simulating reflectors,found at most sites,are speculated as essential pathways for both thermogenic gas and biogenic gas migrating from deep formations to the gas hydrate stability zone.The fracture network on the top of the basement uplift may be heavily gas-charged,which accounts for the chimney with several kilometers in diameter(beneath Plumes B and C).The much smaller gas chimney(beneath Plume D)may stem from gas saturated localized strong permeability zone.High-resolution seismic profiles reveal pipe-like structures,characterized by stacked localized amplitude anomalies,just beneath all the plumes,which act as the fluid conduits conveying gas from the gas hydrate-bearing sediments to the seafloor,feeding the gas plumes.The differences between these pipe-like structures indicate the dynamic process of gas seepage,which may be controlled by the build-up and dissipation of pore pressure.The 3D seismic data show high saturated gas hydrates with high RMS amplitude tend to cluster on the periphery of the gas chimney.Understanding the fluid migration and hydrate accumulation pattern of the Haima cold seeps can aid in the further exploration and study on the dynamic gas hydrate system in the South China Sea.

Active Depths of Main Faults in the Ying-Qiong Basin Investigated by Multi-Scale Wavelet Decomposition of Bouguer Gravity Anomalies and Power Spectral Methods

The Ying-Qiong Basin is located on the northwestern margin of the South China Sea and at the junction of the South China Block and the Indochina Block.It is characterized by complex geological structures.The existing seismic data in the study area is sparse due to the lack of earthquake activities.Because of the limited source energy and poor coverage of seismic data,the knowledge of deep structures in the area,including the spatial distribution of deep faults,is incomplete.Contrarily,satellite gravity data cover the entire study area and can reveal the spatial distribution of faults.Based on the wavelet multi-scale decomposition method,the Bouguer gravity field in the Ying-Qiong Basin was decomposed and reconstructed to obtain the detailed images of the first-to sixth-order gravitational fields.By incorporating the known geological features,the gravitational field responses of the main faults in the Ying-Qiong Basin were identified in the detailed fields,and the power spectrum analysis yielded the depths of 1.4,8,15,26.5,and 39 km for the average burial depths of the bottom surfaces from the first-to fifth-order detailed fields,respectively.The four main faults in the Yinggehai Basin all have a large active depth range:fault A(No.1)is between 5 and 39 km,fault B is between 26.5 and 39 km,and faults C and D are between 15 and 39 km.However,the depth of active faults in the Qiongdongnan Basin is relatively shallow,mainly between 8 and 26.5 km.

Research Progress of Seafloor Pockmarks in Spatio-Temporal Distribution and Classification

Seafloor pockmarks are important indicators of submarine methane seepages and slope instabilities.In order to promote the understanding of submarine pockmarks and their relationship with sediment instabilities and climate changes,here we summarize the research results of pockmarks in the spatio-temporal distributions and shaping factors.Most of pockmarks occur along active or passive continental margins during the last 25 kyr B.P..Circular and ellipse are the most common forms of pockmarks,whereas pockmarks in a special crescent or elongated shape are indicators of slope instabilities,and ring-shape pockmarks are endemic to the gas hydrate zones.Further researches should be focused on the trigger mechanism of climate changes based on the pockmarks in the high latitudes formed during the deglaciation periods,and the role of gas hydrates in the seafloor evolution should be elucidated.In addition,the feature of pockmarks at their early stage(e.g.,developing gas chimneys and gas driving sedimentary doming)and the relations between pockmarks and mass movements,mud diapirs could be further studied to clarify the influences of rapid methane release from submarine sediments on the atmospheric carbon contents.

Mineralogical and Geochemical Constraints on the Mantle Source Characteristics of Basaltic Lavas from the Central Mariana Trough

The composition of mantle-derived basalts reflects the nature of their mantle source regions,which constrain magma generation and composition.Here we present a new whole rock major and trace elements and phenocryst composition of the basaltic lava in the central Mariana Trough.These data provide insights into the mantle source characteristics affected by subduction components.The rocks range from basalts to basaltic andesites,which have high subduction-mobile element contents(e.g.,K,U,Th,LREE)related to N-MORB.The calculated temperature and depth of magma generation are about 1300℃and 30 km,respectively.Although the results above suggest that the addition of hydrous fluid and/or a melt derived from a slab decreases the temperature of mantle partial melting and mildly modifies the composition of a mantle source,the mantle source lithology from which primary magma is generated remains to be peridotite.

Grain size characteristics of the carbonate-free fraction of surface sediments from the Southwest Indian Ridge area and its geological significance

The carbonate-free fraction of 20 surface sediments collected from the ultraslow-spreading Southwest Indian Ridge(SWIR) was studied by grain size analysis and mineralogical analysis with X-ray powder diffraction(XRD),stereo microscopy and scanning electron microscopy(SEM). The characteristics of the carbonate-free fraction of the sediments were obtained, and related influential factors were discussed. The results show that the mean grain size of this fraction is in 1.96Φ–8.19Φ, with poorly sorting and unimodal, bimodal or irregular bimodal distribution patterns. Four grain size end members of the fraction are derived with the End Member Model method. The finest end member EM1 shows a significant contribution of terrigenous materials of the aeolian input and sediment carried by the bottom current. End member EM2 with medium size mainly reflects sediment of a siliceous bioclast origin. EM3 and EM4 are interpreted as representing the coarser volcanic materials related to bedrock weathering or volcanic activities. Multi-provenance is the dominant factor controlling the grain size pattern of the carbonate-free fraction of the sediments in that area. In addition, sediment transport processes such as the bottom current and wind are the minor factors that influence the grain size distribution of the carbonate-free fraction sediments.

Stratigraphic Sequence and Sedimentary Systems in the Middle-Southern Continental Slope of the East China Sea from Seismic Reflection Data: Exploration Prospects of Gas Hydrate

Many evidences for gas hydrate bearing sediments had been found in the continental slope of the East China Sea,such as bottom simulating reflections(BSRs),undersea gas springs,pyrite associated with methane leakage,mud diapirs/mud volcanos,bottom-water methane anomalies and so on.In this study,six key stratigraphic interfaces including T_0(seafloor),T_1(LGM,23 kyr B.P.),T_2(2.58 Myr),T_3(5.33 Myr),T_4(11.02 Myr)and T_5(16.12 Myr)were identified,and then five third-order sequences of SQIII1 to SQIII5 were divided.However,T5 in southern continental slope is not found,which shows that the middle-northern Okinawa Trough had begun to rift in the early Miocene,earlier than the southern segment.Four system tracts including lowstand systems tract(LST),transgressive systems tract(TST),highstand systems tract(HST)and falling stage systems tract(FSST)are further divided.The marine erosion interface of 11.02 Myr and regressive unconformity interface of 23 kyr B.P.indicate two large-scale sea level drop events in the research area.Seven typical seismic facies identified in the continental slope are continental shelf-edge deltas,littoral fluvial-delta plains,incised channels or submarine canyons,slope fans,submarine fans or coastal sandbars,littoral-neritic finegrained sediments,mud volcanos and some other geological bodies respectively.The minimum water depth for hydrate occurrence in the Okinawa Trough is 630 m,and the thickness of gas hydrate stability zone in continental slope is between 0 and 590 m.The calculated bottom boundary of hydrate stability zone is slightly deeper than BSRs on the seismic sections.The re-depositional turbidite sand bodies,such as canyon channels,slope fans and submarine fans developed in Quaternary strata,are the predominant hydrate reservoirs.According to developing process,the dynamic accumulation of hydrate systems can be divided into three evolutionary stages including canyon erosion and hydrate stability zone migration stage,sediments destabilizing and methane leakage stage,and channel filling and hydrate re-occurrence stage.