Mass loss caused by glacier calving is one of the direct contributors to global sea level rise. Reliable calving laws are required for accurate modelling of ice sheet mass balance. Both continuous and discontinuous methods have been used for glacial calving simulations. In this study, the discrete element method (DEM) based on dilated polyhedral elements is introduced to simulate the calving process of a tidewater glacier. Dilated polyhedrons can be obtained from the Minkowski sum of a sphere and a core polyhedron. These elements can be utilized to generate a continuum ice material, where the interaction force between adjacent elements is modeled by constructing bonds at the joints of the common faces. A hybrid fracture model considering fracture energy is introduced. The viscous creep behavior of glaciers on long-term scales is not considered. By applying buoyancy and gravity to the modelled glacier, DEM results show that the calving process is caused by cracks which are initialized at the top of the glacier and spread to the bottom. The results demonstrate the feasibility of using the dilated polyhedral DEM method in glacier simulations, additionally allowing the fragment size of the breaking fragments to be counted. The relationship between crack propagation and internal stress in the glacier is analyzed during calving process. Through the analysis of the Mises stress and the normal stress between the elements, it is found that geometric changes caused by the glacier calving lead to the redistribution of the stress. The tensile stress between the elements is the main influencing factor of glacier ice failure. In addition, the element shape, glacier base friction and buoyancy are studied, the results show that the glacier model based on the dilated polyhedral DEM is sensitive to the above conditions.

The sensitive area of targeted observations for short-term (7 d) prediction of vertical thermal structure (VTS) in summer in the Yellow Sea was investigated. We applied the Conditional Nonlinear Optimal Perturbation (CNOP) method and an adjoint-free algorithm with the Regional Ocean Modeling System (ROMS). We used vertical integration of CNOP-type temperature errors to locate the sensitive areas, where reduction of initial errors is expected to yield the greatest improvement in VTS prediction for the selected verification area. The identified sensitive areas were northeast−southwest orientated northeast to the verification area, which were possibly related to the southwestward background currents. Then, we performed a series of sensitivity experiments to evaluate the effectiveness of the identified sensitive areas. Results show that initial errors in the identified sensitive areas had the greatest negative effect on VTS prediction in the verification area compared to errors in other areas (e.g., the verification area and areas to its east and northeast). Moreover, removal of initial errors through deploying simulated observations in the identified sensitive areas led to more refined prediction than correction of initial conditions in the verification area itself. Our results suggest that implementation of targeted observation in the CNOP-based sensitive areas is an effective method to improve short-term prediction of VTS in summer in the Yellow Sea.

Using the hydrographic data obtained during two nearly simultaneous surveys in June 2015, we carried out semi-diagnostic calculations with the help of a finite element model and a modified inverse method, to study the circulation in the northern South China Sea (NSCS) during the early summer of 2015. A number of new circulation features were found. (1) In most of the observation region, a large, basin-scale anticyclonic gyre appeared south of the 50-m isobath, which contained anticyclonic eddies. One anticyclonic eddy existed from the sea surface to 50-m depth, whose center showed no tilt, while the center of another eddy tilted eastward from the sea surface to 500-m depth. In the eastern part of the observation region, which is west of the Dongsha Islands, there was a sub-basin-scale cyclonic gyre containing a cyclonic eddy whose center tilted southward from the sea surface to 200-m depth. (2) There was a cross-continental slope current (CCSC) in the area southwest of the Dongsha Islands. Its volume transport was about 2.0×10⁶ m³/s. (3) From the estimated order of magnitude of the stream function equation, the joint effect term of the baroclinity and relief (JEBAR) and β-effect term are two important dynamic mechanisms affecting the variation of the circulation in the NSCS. (4) The JEBAR, as a transport-generating term, resulted in the dynamic mechanism determining the pattern of the depth-averaged flow across the contours of potential vorticity fH^–1. Furthermore, we show that the negative values of the JEBAR were the most dominant dynamic mechanism, causing the CCSC southwest of the Dongsha Islands to deflect from the isobaths and veer toward the deep water. The CCSC around the Dongsha Islands was located further southwest during the early summer of 2015 than during the fall of 2005 (revealed by a published study), which suggests that the location of the CCSC around the Dongsha Islands may vary with season.

Satellite-tracked Lagrangian drifters are used to investigate the transport pathways of near-surface water around the Luzon Strait. Particular attention is paid to the intrusion of Pacific water into the South China Sea (SCS). Results from drifter observations suggest that except for the Kuroshio water, other Pacific water that carried by zonal jets, Ekman currents or eddies, can also intrude into the SCS. Motivated by this origin problem of the intrusion water, numerous simulated trajectories are constructed by altimeter-based velocities. Quantitative estimates from simulated trajectories suggest that the contribution of other Pacific water to the total intrusion flux in the Luzon Strait is approximately 13% on average, much smaller than that of Kuroshio water. Even so, over multiple years and many individual intrusion events, the contribution from other Pacific water is quite considerable. The interannual signal in the intrusion flux of these Pacific water might be closely related to variations in a wintertime westward current and eddy activities east of the Luzon Strait. We also found that Ekman drift could significantly contribute to the intrusion of Pacific water and could affect the spreading of intrusion water in the SCS. A case study of an eddy-related intrusion is presented to show the detailed processes of the intrusion of Pacific water and the eddy-Kuroshio interaction.

At the end of May 2008, a massive bloom of macroalgae occurred in the western Yellow Sea off China and lasted for nearly two months, and annual blooms have occurred since then on. During bloom period, the surface-drifting macroalgae have showed an interesting pattern dominated by a banded structure, and the distance between neighboring bands ranged from hundreds of meters to about 6 km with a peak at 1–1.5 km, which is an order of higher than the scale of Langmuir circulation of 50–100 m. In order to explain this new phenomenon, ocean current data obtained from a Doppler current profiler off Qingdao was used to implement stability analysis. By numerically solving the resulting differential Orr-Sommerfeld equation, the secondary circulation induced from the instability of the Emkan current was found to fit well with the observed spatial scale of the surface-drifting macroalgae’s banded structure. As the wind driven Emkan current exist universally in the global ocean, it is reasonable to conclude that the banded structure with kilometers distance between adjoining bands is ubiquitous. We found a new circulation in the upper ocean which is important for exchange of energy, materials and gas between the upper ocean and subsurface layer.

Based on a coupled ocean-sea ice model, this study investigates how changes in the mean state of the atmosphere in different CO₂ emission scenarios (RCP 8.5, 6.0, 4.5 and 2.6) may affect the sea ice in the Bohai Sea, China, especially in the Liaodong Bay, the largest bay in the Bohai Sea. In the RCP 8.5 scenario, an abrupt change of the atmospheric state happens around 2070. Due to the abrupt change, wintertime sea ice of the Liaodong Bay can be divided into 3 periods: a mild decreasing period (2021–2060), in which the sea ice severity weakens at a near-constant rate; a rapid decreasing period (2061–2080), in which the sea ice severity drops dramatically; and a stabilized period (2081–2100). During 2021–2060, the dates of first ice are approximately unchanged, suggesting that the onset of sea ice is probably determined by a cold-air event and is not sensitive to the mean state of the atmosphere. The mean and maximum sea ice thickness in the Liaodong Bay is relatively stable before 2060, and then drops rapidly in the following decade. Different from the RCP 8.5 scenario, atmospheric state changes smoothly in the RCP 6.0, 4.5 and 2.6 scenarios. In the RCP 6.0 scenario, the sea ice severity in the Bohai Sea weakens with time to the end of the twenty-first century. In the RCP 4.5 scenario, the sea ice severity weakens with time until reaching a stable state around the 2070s. In the RCP 2.6 scenario, the sea ice severity weakens until the 2040s, stabilizes from then, and starts intensifying after the 2080s. The sea ice condition in the other bays of the Bohai Sea is also discussed under the four CO₂ emissions scenarios. Among atmospheric factors, air temperature is the leading one for the decline of the sea ice extent. Specific humidity also plays an important role in the four scenarios. The surface downward shortwave/longwave radiation and meridional wind only matter in certain scenarios, while effects from the zonal wind and precipitation are negligible.

Using 110 newly revised Holocene sea level indicators categorized into three types, sediments (67), shelly cheniers (27) and oyster reefs (16), this paper firstly provides a Holocene relative sea level curve, based on multiple approaches of litho- and biostratigraphies and sedimentary faces analysis, for the west coast of Bohai Bay, China. Following considerations, including indicative meaning, the paleo tidal pattern and range and conversion from mean tidal level to mean sea level, an apparent relative mean sea level (RMSL) curve was further reconstructed. After systematical calibration using CALIB, those of the 48 reworked samples were further corrected for the residence-time effect. Similarly, the younger ages for another 35 samples were chosen at the subsample level. These result in a younger-oriented shift for about 0.5 ka. Three local spatial factors, including neotectonic subsidence (average rate about 0.1 mm/a), self-compaction of unconsolidated sediments (between a few decimeters to about 6 m) and subsidence due to groundwater withdrawal (between a few centimeters to about 2.5 m), were quantitatively corrected. Finally, the amended RMSL curve after eliminating all these local temporo-spatial factors is very likely to show non-existence of mid-Holocene highstand and imply potential influences of both ice-volume equivalent sea level and regional glacial isostatic adjustment. Although it is still unable to divide both global and regional factors, the slowdown of sea level rise, in 7.5–6.8 ka with a maximum height less than +1 m, may constrain the model effort in the near future.

The Mozambique Ridge (MOZR) is one of the basement high structures located in the Southwest Indian Ocean, parallel to the Southeast African continental margin. It was formed as a result of the tectono-magmatic evolution of the Gondwana breakup. The origin of the MOZR has been highly debated, with models suggesting either continental or oceanic origin. With new free-air gravity anomaly and multichannel seismic (MCS) reflection data, we present results of 2D density modeling along two seismic profiles acquired by R/V Xiangyanghong 10 at the northern Mozambique Ridge (N-MOZR) between 26°S and 28°S. We observed high free-air gravity anomaly and strong positive magnetic anomaly related to the emplaced seaward dipping reflectors (SDR) and high density lower crustal body (HDLCB), and high Bouguer gravity anomaly associated with the thinning of the continental crust underneath the N-MOZR over a distance of ~82 km. This suggests a thinned and intruded continental crust bound by the Mozambique Fracture Zone (MFZ) that is characterized by gravity low and negative magnetic anomaly. This fracture zone marks the continent-ocean boundary (COB) while the N-MOZR is the transform margin high, i.e., marks the continent-ocean transition (COT) of the Southern Mozambique margin, following the definition of transform margins. We suggest that the N-MOZR was formed by continental extension and subsequent breakup of the MFZ, accompanied by massive volcanism during the southward movement of the Antarctica block. The presence of SDR, HDLCB, and relatively thick oceanic crust indicates the volcanic nature of this transform margin.

A high resolution one-dimensional thermodynamic snow and ice (HIGHTSI) model was used to model the annual cycle of landfast ice mass and heat balance near Zhongshan Station, East Antarctica. The model was forced and initialized by meteorological and sea ice in situ observations from April 2015 to April 2016. HIGHTSI produced a reasonable snow and ice evolution in the validation experiments, with a negligible mean ice thickness bias of (0.003±0.06) m compared to in situ observations. To further examine the impact of different snow conditions on annual evolution of first-year ice (FYI), four sensitivity experiments with different precipitation schemes (0, half, normal, and double) were performed. The results showed that compared to the snow-free case, the insulation effect of snow cover decreased bottom freezing in the winter, leading to 15%–26% reduction of maximum ice thickness. Thick snow cover caused negative freeboard and flooding, and then snow ice formation, which contributed 12%–49% to the maximum ice thickness. In early summer, snow cover delayed the onset of ice melting for about one month, while the melting of snow cover led to the formation of superimposed ice, accounting for 5%–10% of the ice thickness. Internal ice melting was a significant contributor in summer whether snow cover existed or not, accounting for 35%–56% of the total summer ice loss. The multi-year ice (MYI) simulations suggested that when snow-covered ice persisted from FYI to the 10th MYI, winter congelation ice percentage decreased from 80% to 44% (snow ice and superimposed ice increased), while the contribution of internal ice melting in the summer decreased from 45% to 5% (bottom ice melting dominated).

Wave radiation stress is the main driving force of wave-induced near-shore currents. It is directly related to the hydrodynamic characteristics of near-shore current whether the calculation of wave radiation stress is accurate or not. Irregular waves are more capable of reacting wave motion in the ocean compared to regular waves. Therefore, the calculation of the radiation stress under irregular waves will be more able to reflect the wave driving force in the actual near-shore current. Exact solution and approximate solution of the irregular wave radiation stress are derived in this paper and the two kinds of calculation methods are compared. On the basis of this, the experimental results are used to further verify the calculation of wave energy in the approximate calculation method. The results show that the approximate calculation method of irregular wave radiation stress has a good accuracy under the condition of narrow-band spectrum, which can save a lot of computing time, and thus improve the efficiency of calculation. However, the exact calculation method can more accurately reflect the fluctuation of radiation stress at each moment and each location.