The dispersion characteristics of shallow water can be described by the dispersion curves, which contain substantial ocean parameter information. A fast ocean parameter inversion method based on dispersion curves with a single hydrophone is presented in this paper. The method is achieved through Bayesian theory. Several sets of dispersion curves extracted from measured data are used as the input function. The inversion is performed by matching a replica calculated with a dispersion formula. The bottom characteristics can be described by the bottom reflection phase shift parameter P. The propagation range and the depth can be inverted quickly when the seabed parameters are represented by on parameter P. The inversion results improve the inversion efficiency of the seabed parameters. Consequently, the inversion efficiency and accuracy are improved while the number of inversion parameters is decreased and the computational speed of replica is increased. The inversion results have lower error than the reference values, and the dispersion curves calculated with inversion parameters are also in good agreement with extracted curves from measured data; thus, the effectiveness of the inversion method is demonstrated.

In variational methods, coupled parameter optimization (CPO) often needs a long minimization time window (MTW) to fully incorporate observational information, but the optimal MTW somehow depends on the model nonlinearity. The analytical four-dimensional ensemble-variational (A-4DEnVar) considers model nonlinearity well and avoids adjoint model. It can theoretically be applied to CPO. To verify the feasibility and the ability of the A-4DEnVar in CPO, “twin” experiments based on A-4DEnVar CPO are conducted for the first time with the comparison of four-dimensional variational (4D-Var). Two algorithms use the same background error covariance matrix and optimization algorithm to control variates. The experiments are based on a simple coupled ocean-atmosphere model, in which the atmospheric part is the highly nonlinear Lorenz-63 model, and the oceanic part is a slab ocean model. The results show that both A-4DEnVar and 4D-Var can effectively reduce the error of state variables through CPO. Besides, two methods produce almost the same results in most cases when the MTW is less than 560 time steps. The results are similar when the MTW is larger than 560 time steps and less than 880 time steps. The largest MTW of 4D-Var and A-4DEnVar are 1 200 time steps. Moreover, A-4DEnVar is not sensitive to ensemble size when the MTW is less than 720 time steps. A-4DEnVar obtains satisfactory results in the case of highly nonlinear model and long MTW, suggesting that it has the potential to be widely applied to realistic CPO.

The velocity and direction of internal waves (IWs) are important parameters of the ocean, however, traditional observation methods can only obtain the average parameters of IWs for a single location or large area. Herein, a new method based on optical flow is proposed to derive the phase velocity vectors of IWs from X-band marine radar images. First, the X-band marine radar image sequence is averaged, and ramp correction is used to reduce the attenuation of gray values with increasing radial range. Second, the average propagation direction of the IWs is determined using the two-dimensional Fourier transform of the radar images; two radial profiles along this direction are selected from two adjacent radar images; and then, the average phase velocity of the IWs is estimated from these radial profiles. Third, the averaged radar images are processed via histogram equalization and binarization to reduce the influence of noise on the radar images. Fourth, a weighting factor is determined using the average phase velocity of a reference point; the phase velocities on the wave crest of the IWs are subsequently estimated via the optical flow method. Finally, the proposed method is validated using X-band marine radar image sequences observed on an oil platform in the South China Sea, and the error of the phase velocity is calculated to be 0.000 3–0.073 8 m/s. The application conditions of the proposed method are also discussed using two different types of IW packets.

Laboratory visual detection on the hydrate accumulation process provides an effective and low-cost method to uncover hydrate accumulation mechanisms in nature. However, the spatial hydrate distribution and its dynamic evolutionary behaviors are still not fully understood due to the lack of methods and experimental systems. Toward this goal, we built a two-dimensional electrical resistivity tomography (ERT) apparatus capable of measuring spatial and temporal characteristics of hydrate-bearing porous media. Beach sand (0.05–0.85 mm) was used to form artificial methane hydrate-bearing sediment. The experiments were conducted at 1°C under excess water conditions and the ERT data were acquired and analyzed. This study demonstrates the utility of the ERT method for hydrate mapping in laboratory-scale. The results indicate that the average electrical conductivity decreases nonlinearly with the formation of the hydrate. At some special time-intervals, the average conductivity fluctuates within a certain scope. The plane conductivity fields evolve heterogeneously and the local preferential hydrate-forming positions alternate throughout the experimental duration. We speculate that the combination of hydrate formation itself and salt-removal effect plays a dominant role in the spatial and temporal hydrate distribution, as well as geophysical parameters changing behaviors during hydrate accumulation.

The canyon system, including 17 small slope-confined canyons in the Shenhu area, northern South China Sea, is significantly characterized by mounded or undulating features on the canyon flanks and canyon heads. However, the mechanism underlying the formation of these features has yet to be elucidated. In previous studies, most of them were interpreted as sediment deformation on the exploration seismic profiles. In this paper, we collected high-resolution bathymetric data, chirp profiles and geotechnical test data to investigate their detailed morphology, internal structures, and origin. The bathymetric data indicated that most mounded seismic units have smooth seafloors and are separated by grooves or depressions. The distance between two adjacent mounded units is only hundreds of meters. On chirp profiles, mounded seismic units usually exhibit chaotic reflections and wavy reflections, of which the crests migrate upslope. The slope stability analysis results revealed that the critical angle of the soil layers in the study area tends to be 9°, indicating that most mounded seismic units on the canyon flanks and heads are stable at present. The terrain characteristics and seismic configurations combined with the slope stability analysis results indicated that most mounded seismic units are not sediment deformation but depositional structures or mixed systems composed of deformation and depositional structures.

Time series measurements (2010–2017) from the Research Moored Array for African–Asian–Australian Monsoon Analysis and Prediction (RAMA) moorings at 15°N, 90°E and 12°N, 90°E are used to investigate the effect of the seasonal barrier layer (BL) on the mixed-layer heat budget in the Bay of Bengal (BoB). The mixed-layer temperature tendency (∂T/∂t) is primarily controlled by the net surface heat flux that remains in the mixed layer (${Q}^{{'}}$) from March to October, while both ${Q}^{{'}}$ and the vertical heat flux at the base of the mixed layer ($ {Q}_{h} $), estimated as the residual of the mixed-layer heat budget, dominate during winter (November–February). An inverse relation is observed between the BL thickness and the mixed-layer temperature ($ \mathrm{M}\mathrm{L}\mathrm{T} $). Based on the estimations at the moorings, it is suggested that when the BL thickness is ≥25 m, it exerts a considerable influence on ∂T/∂t through the modulation of $ {Q}_{h} $ (warming) in the BoB. The cooling associated with $ {Q}_{h} $ is strongest when the BL thickness is ≤10 m with the $ \mathrm{M}\mathrm{L}\mathrm{T} $ exceeding 29°C, while the contribution from $ {Q}_{h} $ remains nearly zero when the BL thickness varies between 10 m and 25 m. Temperature inversion is evident in the BoB during winter when the BL thickness remains ≥25 m with an average MLT<28.5°C. Furthermore, $ {Q}_{h} $ follows the seasonal cycle of the BL at these RAMA mooring locations, with r>0.72 at the 95% significance level.

Basal melting is an important factor affecting the stability of the ice shelf. The basal channel is formed from uneven melting, which also has an important impact on the stability of the ice shelf. Therefore, it has important scientific value to study the basal channel changes. This study combined datasets of Mosaics of Antarctica, Reference Elevation Model of Antarctica (REMA) and Operation IceBridge to study the temporal and spatial changes of basal channels at the Getz Ice Shelf in Antarctica. The relationships between the cross-sectional area and width of basal channel and those of its corresponding surface depression were statistically analyzed. Then, the changes of the basal channels of Getz Ice Shelf were derived from the ICESat observations and REMA digital elevation models (DEMs). After a detailed analysis of the factors affecting the basal channel changes, we found that the basal channels of Getz Ice Shelf were mainly concentrated in the eastern of the ice shelf, and most of them belonged to the ocean-sourced basal channel. From 2009 to 2016, the total length of the basal channel has increased by approximately 60 km. Affected by the warm Circumpolar Deep Water (CDW), significant changes in the basal channel occurred in the middle reaches of the Getz Ice Shelf. The change of the basal channels at the edge of the Getz Ice Shelf is significantly weaker than that in its middle and upper reaches. Especially in 2005–2012, the eastward wind on the ocean wind field and the westward wind around the continental shelf caused the invasion and upwelling of CDW. Meanwhile, the continuous warming of deep seawater also caused the deepening of the basal channel. During from 2012 to 2020, the fluctuations of the basal channels seem to be caused by the changes in temperature of CDW.

The Dongfang1-1 gas field (DF1-1) in the Yinggehai Basin is currently the largest offshore self-developed gas field in China and is rich in oil and gas resources. The second member of the Pliocene Yinggehai Formation (YGHF) is the main gas-producing formation and is composed of various sedimentary types; however, a clear understanding of the sedimentary types and development patterns is lacking. Here, typical lithofacies, logging facies and seismic facies types and characteristics of the YGHF are identified based on high-precision 3D seismic data combined with drilling, logging, analysis and testing data. Based on 3D seismic interpretation and attribute analysis, the origin of high-amplitude reflections is clarified, and the main types and evolution characteristics of sedimentary facies are identified. Taking gas formation upper II (IIU) as an example, the plane distribution of the delta front and bottom current channel is determined; finally, a comprehensive sedimentary model of the YGHF second member is established. This second member is a shallowly buried “bright spot” gas reservoir with weak compaction. The velocity of sandstone is slightly lower than that of mudstone, and the reflection has medium amplitude when there is no gas. The velocity of sandstone decreases considerably after gas accumulation, resulting in an increase in the wave impedance difference and high-amplitude (bright spot) reflection between sandstone and mudstone; the range of high amplitudes is consistent with that of gas-bearing traps. The distribution of gas reservoirs is obviously controlled by dome-shaped diapir structural traps, and diapir faults are channels through which natural gas from underlying Miocene source rocks can enter traps. The study area is a delta front deposit developed on a shallow sea shelf. The lithologies of the reservoir are mainly composed of very fine sand and coarse silt, and a variety of sedimentary structural types reflect a shallow sea delta environment; upward thickening funnel type, strong toothed bell type and toothed funnel type logging facies are developed. In total, 4 stages of delta front sand bodies (corresponding to progradational reflection seismic facies) derived from the Red River and Blue River in Vietnam have developed in the second member of the YGHF; these sand bodies are dated to 1.5 Ma and correspond to four gas formations. During sedimentation, many bottom current channels (corresponding to channel fill seismic facies) formed, which interacted with the superposed progradational reflections. When the provenance supply was strong in the northwest, the area was dominated by a large set of delta front deposits. In the period of relative sea level rise, surface bottom currents parallel to the coastline were dominant, and undercutting erosion was obvious, forming multistage superimposed erosion troughs. Three large bottom current channels that developed in the late sedimentary period of gas formation IIU are the most typical.

Coastal wetlands are located in the ecotone of interaction between the land surface and sea, and anthropogenic activities extensively interfere with these wetlands through the reclamation of large tidal wetlands and destruction of the function of the ecosystems. In this study, we investigated the dynamic evolutionary characteristics of the Bohai Rim coastal area over the past 40 years using the Modified Normalized Difference Water Index, the fractal dimension, object-oriented classification, the land-use transfer trajectory, and regression analysis. Additionally, we quantified and monitored the evolution of reclamation and analyzed the correlation between reclamation and coastal wetlands based on 99 Landsat-2, -5, and -8 images (at 60 m and 30 m spatial resolution) over the period 1980–2019. The results are as follows. (1) The coastline of the Bohai Rim increased by 1 631.2 km from 1980 to 2019 with a zigzag variation. The artificial coastline increased by 2 946.1 km, whereas the natural coastline decreased by 90%. (2) The area of man-made wetlands increased by 3 736.9 km², the area of construction land increased by 1 008.4 km², and the natural wetland area decreased by 66%. The decrease of tidal flats is the main contributor to the decrease of natural wetland area (takes account for 91.1%). Coastal areas are affected by intense human disturbance, which was taken place across a large area of tidal flats and caused the landscape to fragment and be more heterogeneous. The coastal zone development activities were primarily concentrated in the southern Laizhou Bay, the Yellow River Delta, the Bohai Bay, the northern Liaodong Bay, and the Pulandian Bay. The solidified shorelines and increase in sea level have resulted in intertidal wetlands decreasing and impaired wetland ecology. (3) There is a good agreement between reclamation and the size of the coastal wetlands. Both land reclamation and the reduction in coastal wetland areas are significantly related to the population size, fishery output value, and urbanization rate. In summary, human activities, such as the construction of aquaculture ponds and salt pans, industrialization, and urbanization, are the primary forces that influence the environmental changes in the coastal region. This study is beneficial for establishing and improving the systems for the rational development and utilization of natural resources, and provides theoretical references for restoring wetland ecology and managing future reclamation activities in other coastal zone-related areas.