Sea surface temperature (SST) measurements from 26 coastal hydrological stations of China during 1960–2015 were homogenized and analyzed in this study. The homogenous surface air temperature (SAT) series from meteorological stations which were highly correlated to SST series was used to construct the reference series. Monthly mean SST series were then derived and subjected to a statistical homogeneity test, called penalized maximal t test. Homogenized monthly mean SST series were obtained by adjusting all significant change points which were supported by historic metadata information. Results show that the majority of break points are caused by instrument change and station relocation, which accounts for about 61.3% and 24.2% of the total break points, respectively. The regionally averaged annual homogeneous SST series from the 26 stations shows a warming trend (0.19°C per decade). This result is consistent with that based on the homogenized annual mean SAT at the same region (0.22°C per decade), while the regionally averaged mean original SST series from the same stations shows a much weaker warming of 0.09°C per decade for 1960–2015. This finding suggests that the effects of artificial change points on the result of trend analysis are remarkable, and the warming rate from original SST observations since 1960 may be underestimated. Thus a high quality homogenized observation is crucial for robust detection and assessment of regional climate change. Furthermore, the trends of the seasonal mean homogenized SST were also analyzed. This work confirmed that there was an asymmetric seasonal temperature trends in the Chinese coastal water in the past decades, with the largest warming rate occurring in winter. At last, the significant warming in winter and its relationships to the variability of three large-scale atmospheric modes were investigated.

The newly discovered East Longjing-2 hydrothermal field (ELHF-2) is located on the Dragon Horn oceanic core complex of the ultraslow-spreading Southwest Indian Ridge, approximately 12 km from the ridge axis. This study measured the chemical compositions of pyrite from ELHF-2 using a laser ablation inductively coupled plasma mass spectrometry (LA-ICP-MS) to investigate the genesis of the field. Three generations of pyrite were classified, and found that: Py1 and Py2, rich in V, Mn, U, and Se, occur in altered basalt debris and the silica alteration matrix, respectively. Py3 was mainly intergrown with chalcopyrite in quartz veins and had higher Cu, In, Ag, Sb, and Au contents than Py1 and Py2. Some elements, such as Au, Se, and Pb, are likely presented as direct substitution with Fe²⁺ in pyrite, while Cu, Zn, Co, Ni, and Ag probably occur both as direct substitution with Fe and as distributed micro- to nanoparticle-sized sulfides. Meanwhile, the occurrence of V, Mn, and U is likely presented as oxide inclusions. Trace element geochemistry suggested that the pyrite was formed under high-temperature conditions, and the ore forming elements were likely derived from ultramafic rocks. In addition, Py1 and Py2 were formed under higher water/rock ratio and higher temperature conditions, with more seawater involvement compared with Py3. The formation of ELHF-2 was probably driven by exothermic serpentinization reactions with an additional magmatic heat. This study shows that high-temperature hydrothermal circulation driven by magmatic activity can be developed on distal rift flank areas of magma-starved ultraslow-spreading ridges.

Satellite altimeter needs to be calibrated to evaluate the accuracy of sea surface height data. The dedicated altimeter calibration field needs to establish a special calibration strategy and needs to evaluate its calibration ability. This paper describes absolute calibration of HY-2B altimeter SSH using the GPS calibration method at the newly Wanshan calibration site, located in the Wanshan Islands, China. There are two HY-2B altimeter passes across the Wanshan calibration site. Pass No. 362 is descending and the ground track passes the east of Dan’gan Island. Pass No. 375 is ascending and crosses the Zhiwan Island. The GPS data processing strategy of Wanshan calibration site was established and the accuracy of GPS calibration method of Wanshan calibration site was evaluated. Meanwhile, the processing strategies of the HY-2B altimeter for the Wanshan calibration site were established, and a dedicated geoid model data were used to benefit the calibration accuracy. The time-averaged HY-2B altimeter bias was approximately 2.12 cm with a standard deviation of 2.08 cm. The performance of the HY-2B correction microwave radiometer was also evaluated in terms of the wet troposphere path delay and showed a mean difference −0.2 cm with a 1.4 cm standard deviation with respect to the in situ GPS radiosonde.

A series of sub-parallel linear glacial scours are identified on the crest of the Baoshi Seamount in the Northwind Abyssal Plain by compiling new multibeam data acquired during the 9th Chinese Arctic Research Expedition (CHINARE-Arc9) in 2018 and previously published data. The new data reveal scours that developed at water depths of 850–1 030 m with an orientation of about 75°/255°. The maximum water depth occurs in the southernmost scour and is deeper than that from previous investigations, which showed a maximum scouring depth of about 900 m on the seamount. The topographic and geomorphological characteristics suggest that these scours resulted from erosion by the ice shelf extending from the Chukchi margin and/or Laurentide Ice Sheet that grounded on the crest of the seamount and moved in a NE–SW direction. Other possibilities of their genesis include armadas of large icebergs/multi-keel icebergs calved from the Chukchi Shelf or the Laurentide Ice Sheet. The new data provide new constraints for assessing the extent and volume of the ice sheet in the Chukchi area during glacial maxima.

Long term in situ atmospheric observation of the landfast ice nearby Zhongshan Station in the Prydz Bay was performed from April to November 2016. The in situ observation, including the conventional meteorological elements and turbulent flux, enabled this study to evaluate the sea ice surface energy budget process. Using in situ observations, three different reanalysis datasets from the European Centre for Medium-Range Weather Forecasts Interim Re-analysis (ERA-Interim), National Centers for Environmental Prediction Reanalysis2 (NCEP R2), and Japanese 55-year Reanalysis (JRA55), and the Los Alamos sea ice model, CICE, output for surface fluxes were evaluated. The observed sensible heat flux (SH) and net longwave radiation showed seasonal variation with increasing temperature. Air temperature rose from the middle of October as the solar elevation angle increased. The ice surface lost more energy by outgoing longwave radiation as temperature increased, while the shortwave radiation showed obvious increases from the middle of October. The oceanic heat flux demonstrated seasonal variation and decreased with time, where the average values were 21 W/m² and 11 W/m², before and after August, respectively. The comparisons with in situ observations show that, SH and LE (latent heat flux) of JRA55 dataset had the smallest bias and mean absolute error (MAE), and those of NCEP R2 data show the largest differences. The ERA-Interim dataset had the highest spatial resolution, but performance was modest with bias and MAE between JRA55 and NCEP R2 compare with in situ observation. The CICE results (SH and LE) were consistent with the observed data but did not demonstrate the amplitude of inner seasonal variation. The comparison revealed better shortwave and longwave radiation stimulation based on the ERA-Interim forcing in CICE than the radiation of ERA-Interim. The average sea ice temperature decreased in June and July and increased after September, which was similar to the temperature measured by buoys, with a bias and MAE of 0.9°C and 1.0°C, respectively.

An ensemble optimal interpolation (EnOI) data assimilation method is applied in the BCC_CSM1.1 to investigate the impact of ocean data assimilations on seasonal forecasts in an idealized twin experiment framework. Pseudo-observations of sea surface temperature (SST), sea surface height (SSH), sea surface salinity (SSS), temperature and salinity (T/S) profiles were first generated in a free model run. Then, a series of sensitivity tests initialized with predefined bias were conducted for a one-year period; this involved a free run (CTR) and seven assimilation runs. These tests allowed us to check the analysis field accuracy against the “truth”. As expected, data assimilation improved all investigated quantities; the joint assimilation of all variables gave more improved results than assimilating them separately. One-year predictions initialized from the seven runs and CTR were then conducted and compared. The forecasts initialized from joint assimilation of surface data produced comparable SST root mean square errors to that from assimilation of T/S profiles, but the assimilation of T/S profiles is crucial to reduce subsurface deficiencies. The ocean surface currents in the tropics were better predicted when initial conditions produced by assimilating T/S profiles, while surface data assimilation became more important at higher latitudes, particularly near the western boundary currents. The predictions of ocean heat content and mixed layer depth are significantly improved initialized from the joint assimilation of all the variables. Finally, a central Pacific El Niño was well predicted from the joint assimilation of surface data, indicating the importance of joint assimilation of SST, SSH, and SSS for ENSO predictions.

In this study, the short-term offshore extension of Brahmaputra-Ganges (BG) and Irrawaddy freshwater plumes to the central northern Bay of Bengal (BoB) was investigated based on in situ and satellite observations. In the summer and winter of 2015, two significant freshening events with periods of weeks were observed from a moored buoy at 15°N, 90°E in the BoB. Soil Moisture Active Passive (SMAP) satellite sea surface salinity compares well with the in situ data and shows that these freshening events are directly related to the short-term offshore extension of the BG and Irrawaddy freshwater, respectively. These data combined with the altimeter sea level anomaly data show that the offshore extending plumes result from freshwater modulated by eddies. During summer, the BG freshwater is modulated by a combination of three closely located eddies: a large anticyclonic eddy (ACE) off the northwestern BoB coast and two cyclonic eddies in the northern BoB. Consequently, the freshwater extends offshore from the river mouth and forms a long and narrow tongue-shaped plume extending southwestward to the central BoB. During winter, the Irrawaddy freshwater is modulated by two continuous ACEs evolved from Rossby wave propagating westward from the Irrawaddy Delta off Myanmar, forming a tongue-shaped plume extending to the central BoB. Strong salinity fronts are formed along the boundaries of these tongue-shaped plumes. These findings confirm good capability of the SMAP data to investigate the short-term offshore extension of the BG and Irrawaddy freshwater. This study provides direct evidences of the pathways of the offshore extension of the BG and Irrawaddy freshwater and highlights the role of eddies in the northern BoB freshwater plume variability.

Vegetation in wetlands is a large-scale nature-based resource that can provide multiple benefits to human beings and the environment, such as wave attenuation in coastal zones. Traditionally, there are two main calibration approaches to calculate the attenuation of wave driven by vegetation. The first method is a straightforward one based on the exponential attenuation of wave height in the direction of wave transmission, which, however, overlooks the crucial drag coefficient (C_D). The other method is in accordance with more complicate equations for predicting the damping factor, which is regarded as a function of C_D. In this study, a new relation, combining these above two conventional approaches, is proposed to predict the C_D in an operative approach. Results show that values yielded by the new assessment method perform a strong linear relationship with a collection of historical observations, with a promising R² value of 0.90. Besides, the linear regression derives a new predictive equation for the bulk drag coefficient. Additionally, a calibrated value of 4 for the empirical plant drag coefficient (C_P) is revealed. Overall, this new equation, with the superiority of the convenient exponential regression, is expected to be a rapid assessment method for calculating wave attenuation by vegetation and predicting the drag coefficient.

The response of the mixed layer depth (MLD) and subduction rate in the subtropical Northeast Pacific to global warming is investigated based on 9 CMIP5 models. Compared with the present climate in the 9 models, the response of the MLD in the subtropical Northeast Pacific to the increased radiation forcing is spatially non-uniform, with the maximum shoaling about 50 m in the ensemble mean result. The inter-model differences of MLD change are non-negligible, which depend on the various dominated mechanisms. On the north of the MLD front, MLD shallows largely and is influenced by Ekman pumping, heat flux, and upper-ocean cold advection changes. On the south of the MLD front, MLD changes a little in the warmer climate, which is mainly due to the upper-ocean warm advection change. As a result, the MLD front intensity weakens obviously from 0.24 m/km to 0.15 m/km (about 33.9%) in the ensemble mean, not only due to the maximum of MLD shoaling but also dependent on the MLD non-uniform spatial variability. The spatially non-uniform decrease of the subduction rate is primarily dominated by the lateral induction reduction (about 85% in ensemble mean) due to the significant weakening of the MLD front. This research indicates that the ocean advection change impacts the MLD spatially non-uniform change greatly, and then plays an important role in the response of the MLD front and the subduction process to global warming.

The optimum multiparameter (OMP) method was often used to determine the percentages of water masses based on temperature, salinity and other parameters, like nutrient or dissolved oxygen (DO). There are a number of water masses in the East China Sea (ECS), a marginal sea of the western Pacific Ocean. However, it is difficult to clarify the proportion of water masses using traditional parameters, such as temperature, salinity, nutrient or DO because of the occurring of intensive biogeochemical processes in the near shore and shelf areas. Here, we reported the use of ²³⁴U/²³⁸U activity ratio embedded in the OMP method. The results indicate that seawater in the northern ECS mainly consisted of the estuarine water of Changjiang River (CEW), Kuroshio water (KW), and Yellow Sea Coastal Current (YSCC). In March 2017, the CEW only influenced the offshore waters shallower than 30 m; the KW affected the east edge and the YSCC contributed more than 75% in the northern ECS.