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.

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.

Variations in incoming shortwave radiation influence the net surface heat flux, contributing to the formation of a temperature inversion. The effects of shortwave radiation on the temperature inversions in the Bay of Bengal and eastern equatorial Indian Ocean have never been investigated. Thus, a high-resolution (horizontal resolution of 0.07°×0.07° with 50 vertical layers) Regional Ocean Modeling System (ROMS) model is utilized to quantify the contributions of shortwave radiation to the temperature inversions in the study domain. Analyses of the mixed layer heat and salt budgets are performed, and different model simulations are compared. The model results suggest that a 30% change in shortwave radiation can change approximately 3% of the temperature inversion area in the Bay of Bengal. Low shortwave radiation reduces the net surface heat flux and cools the mixed layer substantially; it also reduces the evaporation rate, causing less evaporative water vapor losses from the ocean than the typical situation, and ultimately enhances haline stratification. Thus, the rudimentary outcome of this research is that a decrease in shortwave radiation produces more temperature inversion in the study region, which is primarily driven by the net surface cooling and supported by the intensive haline stratification. Moreover, low shortwave radiation eventually intensifies the temperature inversion layer by thickening the barrier layer. This study could be an important reference for predicting how the Indian Ocean climate will respond to future changes in shortwave radiation.

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.

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.

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.

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.