The 4.2 ka event that occurred during the period from 4 500–3 900 a BP was characterized by cold and dry climates and resulted in the collapse of civilizations around the world. The cause of this climatic event, however, has been under debate. We collected four corals (Porites lutea) from Yongxing Island, Xisha Islands, South China Sea, dated them with the U-series method, and measured the annual coral growth rates using X-ray technology. The dating results showed that the coral growth ages were from 4 500–3 900 a BP, which coincide well with the period of the 4.2 ka event. We then reconstructed annual sea surface temperature anomaly (SSTA) variations based on the coral growth rates. The growth rate-based SSTA results showed that the interdecadal SSTA from 4 500–3 900 a BP was lower than that during modern times (1961–2008 AD). A spectral analysis showed that the SSTA variations from 4 500–3 900 a BP were under the influence of El Niño-Southern Oscillation (ENSO) activities. From 4 500–4 100 a BP, the climate exhibited La Niña-like conditions with weak ENSO intensity and relatively stable and lower SSTA amplitudes. From 4 100–3 900 a BP, the climate underwent a complicated period of ENSO variability and showed alternating El Niño- or La Niña-like conditions at interdecadal time scales and large SSTA amplitudes. We speculate that during the early and middle stages of the 4.2 ka event, the cold climate caused by weak ENSO activities largely weakened social productivity. Then, during the end stages of the 4.2 ka event, the repeated fluctuations in the ENSO intensity caused frequent extreme weather events, resulting in the collapse of civilizations worldwide. Thus, the new evidence obtained from our coral records suggests that the 4.2 ka event as well as the related collapse of civilizations were very likely driven by ENSO variability.

The Maritime Silk Road is not only a passageway for business and trade, but also the road of friendship between eastern and western civilizations. The Indian Ocean tsunami in 2004 caused major damage to several coastal countries. Tsunami occurrence regularity and hazard analysis are needed to ensure economic and cultural exchange on the Maritime Silk Road. To explore and identify tsunami hazard on the Maritime Silk Road, the spatial and temporal characteristics of historical tsunami events were given out. Some useful information hidden in historical tsunamis was searched from source parameters, such as seismic magnitude, focal depth and water depth. The tsunami possibility in the case of earthquake occurrence was also studied, exploring the probability of tsunami caused by different magnitudes. The analysis result shows that tsunamis on the Maritime Silk Road mainly occurred in 8 major tectonic faults, each of which has different tsunami occurrence regularity. On the basis of statistical analysis, a numerical model was used to simulate the potential tsunamis and show the tsunami hazard levels along the coast of Maritime Silk Road. The research results of this paper can help the tsunami early warning and ensure the safety of economic and cultural exchanges on the Maritime Silk Road.

In this study, the impact of oceanic processes on the sensitivity of transient climate change is investigated using two sets of coupled experiments with and without tidal forcing, which are termed Exp_Tide and Exp_Control, respectively. After introducing tidal forcing, the transient climate response (TCR) decreases from 2.32 K to 1.90 K, and the surface air temperature warming at high latitudes decreases by 29%. Large ocean heat uptake efficiency and heat storage can explain the low TCR in Exp_Tide. Approximately 21% more heat is stored in the ocean in Exp_Tide (1.10×10²⁴J) than in Exp_Control (0.91×10²⁴J). Most of the large ocean warming occurs in the upper 1 000 m between 60°S and 60°N, primarily in the Atlantic and Southern Oceans. This ocean warming is closely related to the Atlantic Meridional Overturning Circulation (AMOC). The initial transport at mid- and high latitudes and the decline in the AMOC observed in Exp_Tide are both larger than those observed in Exp_Control. The spatial structures of AMOC are also different with and without tidal forcing in present experiments. The AMOC in Exp_Tide has a large northward extension. We also investigated the relationship between AMOC and TCR suggested by previous studies using the present experiments.

Wave-current interaction and its effects on the hydrodynamic environment in the Beibu Gulf (BG) have been investigated via employing the Coupled Ocean–Atmosphere–Wave–Sediment Transport (COAWST) modeling system. The model could simulate reasonable hydrodynamics in the BG when validated by various observations. Vigorous tidal currents refract the waves efficiently and make the seas off the west coast of Hainan Island be the hot spot where currents modulate the significant wave height dramatically. During summer, wave-enhanced bottom stress could weaken the near-shore component of the gulf-scale cyclonic-circulation in the BG remarkably, inducing two major corresponding adjustments: Model results reveal that the deep-layer cold water from the southern BG makes critical contribution to maintaining the cold-water mass in the northern BG Basin. However, the weakened background circulation leads to less cold water transported from the southern gulf to the northern gulf, which finally triggers a 0.2°C warming in the cold-water mass area; In the top areas of the BG, the suppressed background circulation reduces the transport of the diluted water to the central gulf. Therefore, more freshwater could be trapped locally, which then triggers lower sea surface salinity (SSS) in the near-field and higher SSS in the far-field.

We analyzed seafloor morphology and geophysical anomalies of the Southeast Indian Ridge (SEIR) to reveal the remarkable changes in magma supply along this intermediate fast-spreading ridge. We found systematic differences of the Australian-Antarctic Discordance (AAD) from adjacent ridge segments with the residual mantle Bouguer gravity anomaly (RMBA) being more positive, seafloor being deeper, morphology being more chaotic, M factors being smaller at the AAD. These systematic anomalies, as well as the observed Na_8.0 being greater and Fe_8.0 being smaller at AAD, suggest relatively starved magma supply and relatively thin crust within the AAD. Comparing to the adjacent ridges segments, the calculated average map-view M factors are relatively small for the AAD, where several Oceanic Core Complexes (OCCs) develop. Close to 30 OCCs were found to be distributed asymmetrically along the SEIR with 60% of OCCs at the northern flank. The OCCs are concentrated mainly in Segments B3 and B4 within the AAD at ~124°–126°E, as well as at the eastern end of Zone C at ~115°E. The relatively small map-view M factors within the AAD indicate stronger tectonism than the adjacent SEIR segments. The interaction between the westward migrating Pacific mantle and the relatively cold mantle beneath the AAD may have caused a reduction in magma supply, leading to the development of abundant OCCs.

Numerical wave tanks are widely-acknowledged tools in studying waves and wave-structure interactions. They can generate waves under realistic scales and offers more information on the fluid field. However, most numerical wave tanks suffer from issues known as the numerical dissipation and numerical dispersion. The former causes wave energy to be slowly dissipated and the latter shifts wave frequencies during wave propagation. This paper proposes a simple method of depressing numerical dissipation effects on the basis of solving Euler equations using the finite difference method (FDM). The wave propagation solutions are solved analytically taking into account the influence of the damping terms. The main idea of the method is to append a source term to the momentum equation, whose strength is determined by how strong the numerical damping effect is. The method is verified by successfully depressing numerical effects during the simulation of regular linear waves, Stokes waves and irregular waves. By applying the method, wave energy is able to be close to its initial value after long distance of travel.

We used satellite altimetry data to investigate the Kuroshio Current because of the higher resolution and wider range of observations. In previous studies, satellite absolute geostrophic velocities were used to study the spatio-temporal variability of the sea surface velocity field along the current, and extraction methods were employed to detect the Kuroshio axes and paths. However, sea surface absolute geostrophic velocity estimated from absolute dynamic topography should be regarded as the geostrophic component of the actual surface velocity, which cannot represent a sea surface current accurately. In this study, mathematical verification between the climatic absolute geostrophic and bin-averaged drifting buoy velocity was established and then adopted to correct the satellite absolute geostrophic velocities. There were some differences in the characteristics between satellite geostrophic and drifting buoy velocities. As a result, the corrected satellite absolute geostrophic velocities were used to detect the Kuroshio axis and path based on a principal-component detection scheme. The results showed that the detection of the Kuroshio axes and paths from corrected absolute geostrophic velocities performed better than those from satellite absolute geostrophic velocities and surface current estimations. The corrected satellite absolute geostrophic velocity may therefore contribute to more precise day-to-day detection of the Kuroshio Current axis and path.

This paper overviews research progress in observation, theoretical analysis and numerical modeling of submesoscale dynamic processes in the South China Sea (SCS) particularly during recent five years. The submesoscale processes are defined according to both spatial and dynamic scales, and divided into four subcategories as submesoscale waves, submesoscale vortexes, submesoscale shelf processes, and submesoscale turbulence. The major new findings are as follows. (1) Systematic mooring observations provide new insights into the solitary waves (ISWs) and the typhoon-forced near-inertial waves (NIWs), of which a new type of ISWs with period of 23 h was observed in the northern SCS (NSCS), and the influences of background vorticity, summer monsoon onset, and deep meridional overturning circulation on the NIWs, as well as nonlinear wave-wave interaction between the NIWs and internal tides, are better understood. On the other hand, satellite altimeter sea surface height data are used to reveal the internal tide radiation patterns and provide solid evidence for that the ISWs in the northeastern SCS originate from the Luzon Strait. (2) Submesoscale offshore jets and associated vortex trains off the Vietnam coast in the western boundary of the SCS were observed from satellite chlorophyll concentration images. Spiral trains with the horizontal scale of 15–30 km and the spacing of 50–80 km were identified.(3) 3-D vertical circulation in the upwelling region east of Hainan Island was theoretically analyzed. The results show that distribution patterns of all the dynamic terms are featured by wave-like structures with horizontal wavelength scale of 20–40 km. (4) Numerical models have been used for the research of submesoscale turbulence. Submesoscale vertical pump of an anticyclonic eddy and the spatiotemporal features of submesoscale processes in the northeastern SCS are well modeled.

An extensive study collected in situ data along the Yellow Sea (YS) and East China Sea (ECS) to assess the radiometric properties and the concentration of the water constituents derived from Moderate Resolution Imaging Spectroradiometer (MODIS). Thirteen high quality match-ups were obtained for evaluating the MODIS estimates of R_rs(λ), chlorophyll a (Chl a) and concentrations of suspended particulate sediment matter (SPM). For MODIS R_rs(λ), the mean absolute percentage difference (APD) was in the range of 20%–36%, and the highest uncertainty appeared at 412 nm, whereas the band ratio of R_rs(λ) at 488 nm compared with that at 547 nm was highly consistent, with an APD of 7%. A combination of near-infrared bands and shortwave infrared wavelengths atmosphere correction algorithm (NIR-SWIR algorithm) was applied to the MODIS data, and the estimation accuracy of R_rs were improved at most of the visible spectral bands except 645 nm, 667 nm and 678 nm. Two ocean-colour empirical algorithms for Chl a estimation were applied to the processed data, the results indicated that the accuracy of the derived Chl a values was obviously improved, the four-band algorithms outperformed the other algorithm for measured and simulated datasets, and the minimum APD was 35%. The SPM was also quantified. Two regional and two coastal SPM algorithms were modified according to the in situ data. By comparison, the modified Tassan model had a higher accuracy for the application along the YS and ECS with an APD of 21%. However, given the limited match-up dataset and the potential influence of the aerosol properties on atmosphere correction, further research is required to develop additional algorithms especially for the low Chl a coastal water.