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.

The stress state and rock mechanical properties govern the growth of faults and fractures, which constitute shallow hydrothermal pathways and control the distribution of seafloor massive sulfide (SMS) mounds in the seafloor hydrothermal field. The stress field has an important influence on the formation and persistence of hydrothermal pathways. Based on multibeam bathymetric data from the Trans-Atlantic Geotraverse (TAG) field, we establish two three-dimensional geological models with different scales to simulate the stress field, which investigate the characteristics of hydrothermal pathways and associated SMS mounds. The simulation results show that oblique faults and fissures form in the tensile stress zone and that mounds, including active and inactive hydrothermal mounds form in the compressive stress zone. Fault activity, which is related to the stress field, affects the opening and closing of hydrothermal channels and changes the permeability structure of subseafloor wall rock. Therefore, the stress field controls the development and persistence of shallow hydrothermal pathways. The features of shallow hydrothermal pathways in the stress field can provide geomechanical information that is useful for identifying favorable zone for SMS deposit formation.

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.

Smaller mesoscale eddies (SMEs) have an important effect on the transmission of ocean temperatures, salinity, energy, and marine biochemical processes. However, traditional altimeters, the dominant sensors used to identify and track eddies, have made it challenging to observe SMEs accurately due to resolution limitations. Eddies drive local upwelling or downwelling, leaving signatures on sea surface temperatures (SSTs) and chlorophyll concentrations (Chls). SST can be observed by spaceborne infrared sensors, and Chl can be measured by ocean color remote sensing. Therefore, multisatellite observations provide an opportunity to obtain information to characterize SMEs. In this paper, an eddy detection algorithm based on SST and Chl images is proposed, which identifies eddies by characterizing the spatial and temporal distribution of SST and Chl data. The algorithm is applied to characterize and analyze SMEs in the Kuroshio Extension. Statistical results on their distribution and seasonal variability are shown, and the formation processes are preliminarily discussed. SMEs generation may be contributed by horizontal strain instability, the interaction of topographic obstacles and currents, and wind stress curl.

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.

This study presents an efficient Boussinesq-type wave model accelerated by a single Graphics Processing Unit (GPU). The model uses the hybrid finite volume and finite difference method to solve weakly dispersive and nonlinear Boussinesq equations in the horizontal plane, enabling the model to have the shock-capturing ability to deal with breaking waves and moving shoreline properly. The code is written in CUDA C. To achieve better performance, the model uses cyclic reduction technique to solve massive tridiagonal linear systems and overlapped tiling/shared memory to reduce global memory access and enhance data reuse. Four numerical tests are conducted to validate the GPU implementation. The performance of the GPU model is evaluated by running a series of numerical simulations on two GPU platforms with different hardware configurations. Compared with the CPU version, the maximum speedup ratios for single-precision and double-precision calculations are 55.56 and 32.57, respectively.

Information on the Arctic sea ice climate indicators is crucial to business strategic planning and climate monitoring. Data on the evolvement of the Arctic sea ice and decadal trends of phenology factors during melt season are necessary for climate prediction under global warming. Previous studies on Arctic sea ice phenology did not involve melt ponds that dramatically lower the ice surface albedo and tremendously affect the process of sea ice surface melt. Temporal means and trends of the Arctic sea ice phenology from 1982 to 2017 were examined based on satellite-derived sea ice concentration and albedo measurements. Moreover, the timing of ice ponding and two periods corresponding to it were newly proposed as key stages in the melt season. Therefore, four timings, i.e., date of snow and ice surface melt onset (MO), date of pond onset (PO), date of sea ice opening (DOO), and date of sea ice retreat (DOR); and three durations, i.e., melt pond formation period (MPFP, i.e., MO–PO), melt pond extension period (MPEP, i.e., PO–DOR), and seasonal loss of ice period (SLIP, i.e., DOO–DOR), were used. PO ranged from late April in the peripheral seas to late June in the central Arctic Ocean in Bootstrap results, whereas the pan-Arctic was observed nearly 4 days later in NASA Team results. Significant negative trends were presented in the MPEP in the Hudson Bay, the Baffin Bay, the Greenland Sea, the Kara and Barents seas in both results, indicating that the Arctic sea ice undergoes a quick transition from ice to open water, thereby extending the melt season year to year. The high correlation coefficient between MO and PO, MPFP illustrated that MO predominates the process of pond formation.

In this paper, we propose a hybrid forecasting model to improve the forecasting accuracy for depth-averaged current velocities (DACVs) of underwater gliders. The hybrid model is based on a discrete wavelet transform (DWT), a deep belief network (DBN), and a least squares support vector machine (LSSVM). The original DACV series are first decomposed into several high- and one low-frequency subseries by DWT. Then, DBN is used for high-frequency component forecasting, and the LSSVM model is adopted for low-frequency subseries. The effectiveness of the proposed model is verified by two groups of DACV data from sea trials in the South China Sea. Based on four general error criteria, the forecast performance of the proposed model is demonstrated. The comparison models include some well-recognized single models and some related hybrid models. The performance of the proposed model outperformed those of the other methods indicated above.

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.