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.

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.

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.

The orthogonal supersegment of the ultraslow-spreading Southwest Indian Ridge at 16°–25°E is characterized by significant along-axis variations of mantle potential temperature. A detailed analysis of multibeam bathymetry, gravity, and magnetic data were performed to investigate its variations in magma supply and crustal accretion process. The results revealed distinct across-axis variations of magma supply. Specifically, the regionally averaged crustal thickness reduced systematically from around 7 Ma to the present, indicating a regionally decreasing magma supply. The crustal structure is asymmetric in regional scale between the conjugate ridge flanks, with the faster-spreading southern flank showing thinner crust and greater degree of tectonic extension. Geodynamic models of mantle melting suggested that the observed variations in axial crustal thickness and major element geochemistry can be adequately explained by an eastward decrease in mantle potential temperature of about 40°C beneath the ridge axis. In this work, a synthesized model was proposed to explain the axial variations of magma supply and ridge segmentation stabilities. The existence of large ridge-axis offsets may play important roles in controlling melt supply. Several large ridge-axis offsets in the eastern section (21°–25°E) caused sustained along-axis focusing of magma supply at the centers of eastern ridge segments, enabling quasi-stable segmentation. In contrast, the western section (16°–21°E), which lacks large ridge-axis offsets, is associated with unstable segmentation patterns.

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.

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.

The methane bubble plume attracts interest because it offers direct evidence of seafloor gas leakage and plays an indirect role in the exploration and identification of natural gas hydrate. In this study, based on established plume models and their migration sections, three amplitude-class attributes were extracted from three formations for the migration sections of five plumes, and the correlation between the gas content and seismic attribute was obtained. As the gas content increases, the amplitude attribute correspondingly increases, and the linear correlation is relatively good. Moreover, correlation coefficients between gas content and amplitude attributes are close to 1.0. By using linear fitting, the relation model between the gas content of the plume and the seismic attribute was obtained. The relation model was subsequently used to invert the gas content from a real data-bearing plume. Comparison of the gas content section of the plume with the attribute section and real seismic section reveals common distribution characteristics, namely, the color of the section in the lower right corner is dark. If the amplitude value is large in the seismic section of the real plume, the amplitude attribute value is also large in the corresponding attribute section, and the inverted value of the gas content is also large (because gas content and amplitude are linearly correlated), which indicates that the plume bubbles of the section in the lower right corner is intensively distributed. Finally, the obtained gas content section of the plume can reflect the distribution of the plume bubble content more simply and intuitively, from which the distribution law of seafloor bubbles can be deduced, and this lays a foundation for the further estimation of the gas content of the plume and hydrate reserves.

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.

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 statistical characteristics and vertical thermohaline properties of mesoscale eddies in the Bay of Bengal are studied from the view of satellite altimetry data and Argo profiles. Eddy propagation preferences in different lifetimes, eddy evolution process, and geographical distribution of eddy kinetic properties are analyzed in this area. Eddies exist principally in the western Bay of Bengal, and most of them propagate westward. There is a clear southward (equatorward) preference for eddies with long lifetimes, especially for cyclones. Moreover, the eddies in different areas of the bay show different north-southward preferences. Evolution of eddy kinetic properties with lifetime shows that eddies have the significant three-stage feature: the growth period in the former one-fifth lifetime, the stable period in the middle two-fifth to four-fifth lifetime, and the dying period in the last one-fifth lifetime. Large-amplitude and high-intensity eddies occur only in the relatively confined regions of highly unstable currents, such as the East Indian Coastal Current and eastern Sri Lanka. Based on Argo profile data and climatology data, the eddy synthesis method was used to construct three-dimensional temperature and salt structures of eddies in this area. The mean temperature anomaly is negative/positive to the cyclonic/anticyclonic eddies in the upper 300×10⁴ Pa, and below this depth, the anomaly becomes weak. The salinity structures of positive anomalies inside cyclonic eddies and negative anomalies inside anticyclonic eddies in the Bay of Bengal are not consistent with other regions. Due to the special characteristics of the water mass in the bay, especially under the control of the low-salinity Bay of Bengal water at the surface and the Indian equatorial water in the deep ocean, the salinity of seawater shows a monotonic increase with depth. For regional varieties of temperature and salinity structures, as the eddies move westward, the temperature anomaly induced by the eddies increases, the effecting depth of the eddies deepens, and the salinity structures are more affected by inflows. In the north-south direction, the salinity structures of the eddies are associated with the local water masses, which comprise low-salinity water in the northern bay due to the inflow of freshwater from rivers and salty water in the southern bay due to the invasion of Arabian Sea high-salinity water from the north Indian Ocean.