Characteristics and triggering mechanisms of early negative Indian Ocean Dipole

Negative Indian Ocean Dipole(nIOD)can exert great impacts on global climate and can also strongly influence the climate in China.Early nIOD is a major type of nIOD,which can induce more pronounced climate anomalies in summer than La Niña-related nIOD.However,the characteristics and triggering mechanisms of early nIOD are unclear.Our results based on reanalysis datasets indicate that the early nIOD and La Niña-related nIOD are the two major types of nIOD,and the former accounts for over one third of all the nIOD events in the past six decades.These two types of nIODs are similar in their intensities,but are different in their spatial patterns and seasonal cycles.The early nIOD,which develops in spring and peaks in summer,is one season earlier than the La Niña-related nIOD.The spatial pattern of the wind anomaly associated with early nIOD exhibits a winter monsoon-like pattern,with strong westerly anomalies in the equatorial Indian Ocean and eastly anomalies in the northern Indian Ocean.Opposite to the triggering mechanism of early positve IOD,the early nIOD is induced by delayed Indian summer monsoon onset.The results of this study are helpful for improving the prediction skill of IOD and its climate impacts.

Long-term ocean temperature trend and marine heatwaves

Marine heatwaves(MHWs)can cause irreversible damage to marine ecosystems and livelihoods.Appropriate MHW characterization remains difficult,because the choice of a sea surface temperature(SST)temporal baseline strongly influences MHW identification.Following a recent work suggesting that there should be a communicating baseline for long-term ocean temperature trends(LTT)and MHWs,we provided an effective and quantitative solution to calculate LTT and MHWs simultaneously by using the ensemble empirical mode decomposition(EEMD)method.The long-term nonlinear trend of SST obtained by EEMD shows superiority over the traditional linear trend in that the data extension does not alter prior results.The MHWs identified from the detrended SST data exhibited low sensitivity to the baseline choice,demonstrating the robustness of our method.We also derived the total heat exposure(THE)by combining LTT and MHWs.The THE was sensitive to the fixed-period baseline choice,with a response to increasing SST that depended on the onset time of a perpetual MHW state(identified MHW days equal to the year length).Subtropical areas,the Indian Ocean,and part of the Southern Ocean were most sensitive to the long-term global warming trend.

The numerical investigation of seasonal variation of the cold water mass in the Beibu Gulf and its mechanisms

A wave-tide-circulation coupled model based on the Princeton Ocean Model is established to explore the seasonal variation of the cold water mass in the Beibu Gulf and its mechanisms. The results show that the cold water mass starts forming in March, reaches the maximum strength during June and July, and fades away since October. Strong mixing in winter transports the cold water from sea surface to bottom. The cold water mass remains in the bottom layer as the thermocline strengthens during spring, except for the shallow water where the themocline is broken by strong tidal mixing, which gradually separate the cold water mass from its surrounding warm water. Further analysis on the ocean current and stream function confirms that the cold water mass in the Beibu Gulf is locally developed, with an anticlockwise circulation caused by a strong temperature gradient. Sensitivity experiments reveal that the cold water mass is controlled by the sea surface heat flux, while the terrain and tidal mixing also play important roles.

Simulations and Projections of Winter Sea Ice in the Barents Sea by CMIP6 Climate Models

Dramatic changes in the sea ice characteristics in the Barents Sea have potential consequences for the weather and climate systems of mid-latitude continents,Arctic ecosystems,and fisheries,as well as Arctic maritime navigation.Simulations and projections of winter sea ice in the Barents Sea based on the latest 41 climate models from the Coupled Model Intercomparison Project Phase 6(CMIP6)are investigated in this study.Results show that most CMIP6 models overestimate winter sea ice in the Barents Sea and underestimate its decreasing trend.The discrepancy is mainly attributed to the simulation bias towards an overly weak ocean heat transport through the Barents Sea Opening and the underestimation of its increasing trend.The methods of observation-based model selection and emergent constraint were used to project future winter sea ice changes in the Barents Sea.Projections indicate that sea ice in the Barents Sea will continue to decline in a warming climate and that a winter ice-free Barents Sea will occur for the first time during 2042-2089 under the Shared Socioeconomic Pathway 585(SSP5-8.5).Even in the observation-based selected models,the sensitivity of winter sea ice in the Barents Sea to global warming is weaker than observed,indicating that a winter ice-free Barents Sea might occur earlier than projected by the CMIP6 simulations.

A 2D-numerical modeling of the generation and propagation of internal solitary waves in the Luzon Strait

The wide presence of internal solitary waves （ISWs） in the northern South China Sea （SCS） has been confirmed by both Synthetic Aperture Radar （SAR） images and in situ observations.These ISWs are believed being generated over the varying topography in the Luzon Strait.They typically propagate westwards into the SCS with a diurnal or semidiurnal period.Their generation sites are,however,not yet solidly identified.To obtain a clear picture of the ISWs,we designed numerical experiments to analyze the generation and propagation of the ISWs in the Luzon Strait using a 2-dimensional non-hydrostatic model.The model current is forced by barotropic or baroclinic currents imposed at open boundaries.The experiments show that the tidal current serves as a kind of triggering force for the ISWs over the submarine ridges in the strait.Under the forcing of tidal currents,depressions are formed near the ridges.The ISWs then split from the depressions through a process different from lee-wave generation mechanism.The appearance of the ISWs is influenced by the strength and period of the forcing current： the ISWs are more likely to be generated by a stronger tidal current.That is why the ISWs in the Luzon Strait are frequently observed during spring tide.Compared with diurnal tidal current,the ISWs generated by semidiurnal tidal current with the same amplitude is much more energetic.It is partly because that the wave beams in diurnal frequency have a larger angle with the vertical direction,thus are more likely to be reflected by the topography slope.The impact of the Kuroshio to the ISWs is also analyzed by adding a vertical uniform or shear current at boundaries.A vertically uniform current may generate ISWs directly.On the other hand,a vertically shear current,which is more realistic to represent the Kuroshio branch,seems to have little influence on the generation process and radiating direction of the ISWs in the Luzon Strait.

Performance of CMIP6 models in simulating the dynamic sea level:Mean and interannual variance

本研究采用卫星测高数据与第六次国际耦合模式比较计划(CMIP6)海平面动力进行对比,重点针对40S-40N地区的动力海平面(DSL),评估了模式对其平均态与年际变率的综合模拟能力,结果表明,对于DSL平均态的模拟,模式与观测结果非常吻合,模式之间的差异较小.其中,副热带北大西洋是模拟偏差和模式间差异较为显著的区域,对于DSL年际变率的模拟,模式之间保持较高的一致性,但是,模式与观测结果存在明显差异,模式普遍低估了DSL的年际方差;其中,误差大值区域出现在副热带西边界流附近,模式分辨率会影响CMIP6对中小尺度海洋过程的重现能力,这可能是导致CMIP6历史模拟出现误差的原因之一.

Evaluation of Nonbreaking Wave-Induced Mixing Parameterization Schemes Based on a One-Dimensional Ocean Model

Surface waves have a considerable effect on vertical mixing in the upper ocean.In the past two decades,the vertical mixing induced through nonbreaking surface waves has been used in ocean and climate models to improve the simulation of the upper ocean.Thus far,several nonbreaking wave-induced mixing parameterization schemes have been proposed;however,no quantitative comparison has been performed among them.In this paper,a one-dimensional ocean model was used to compare the performances of five schemes,including those of Qiao et al.(Q),Hu and Wang(HW),Huang and Qiao(HQ),Pleskachevsky et al.(P),and Ghantous and Babanin(GB).Similar to previous studies,all of these schemes can decrease the simulated sea surface temperature(SST),increase the subsurface temperature,and deepen the mixed layer,thereby alleviating the common thermal deviation problem of the ocean model for upper ocean simulation.Among these schemes,the HQ scheme exhibited the weakest wave-induced mixing effect,and the HW scheme exhibited the strongest effect;the other three schemes exhibited roughly the same effect.In particular,the Q and P schemes exhibited nearly the same effect.In the simulation based on observations from the Ocean Weather Station Papa,the HQ scheme exhibited the best performance,followed by the Q scheme.In the experiment with the HQ scheme,the root-mean-square deviation of the simulated SST from the observations was 0.43℃,and the mixed layer depth(MLD)was 2.0 m.As a contrast,the deviations of the SST and MLD reached 1.25℃ and 8.4 m,respectively,in the experiment without wave-induced mixing.

Seasonal variation of mesoscale eddy intensity in the global ocean

Mesoscale eddies are a prominent oceanic phenomenon that plays an important role in oceanic mass transport and energy conversion.Characterizing by rotational speed,the eddy intensity is one of the most fundamental properties of an eddy.However,the seasonal spatiotemporal variation in eddy intensity has not been examined from a global ocean perspective.In this study,we unveil the seasonal spatiotemporal characteristics of eddy intensity in the global ocean by using the latest satellite-altimetry-derived eddy trajectory data set.The results suggest that the eddy intensity has a distinct seasonal variation,reaching a peak in spring while attaining a minimum in autumn in the Northern Hemisphere and the opposite in the Southern Hemisphere.The seasonal variation of eddy intensity is more intense in the tropical-subtropical transition zones within latitudinal bands between 15°and 30°in the western Pacific Ocean,the northwestern Atlantic Ocean,and the eastern Indian Ocean because baroclinic instability in these areas changes sharply.Further analysis found that the seasonal variation of baroclinic instability precedes the eddy intensity by a phase of 2–3 months due to the initial perturbations needing time to grow into mesoscale eddies.

Current observation and analysis based on mooring systems in the Andaman Sea

The basic structure and intraseasonal evolution of currents in the southeastern Andaman Sea was analyzed based on data collected in 2017 from two subsurface moorings(C1 and C5).Periodic variation in the upper ocean currents of the Andaman Sea was investigated by combining observational and satellite data.Mooring observations show that rapid changes of current speed and direction occurred in May and June,with a significant increase in current velocity at the C1 mooring.In the second half of the year,southward flow dominated at the C1 mooring,and alternating northward and southward flows were evident at the C5 mooring during the same period but the northward flow prevailed in boreal winter.In addition,analysis of the power spectra of the upper currents revealed that the tidal period at both moorings is primarily semidiurnal with weaker energy than that of the low-frequency currents.The upper ocean currents at the C1 and C5 moorings exhibited intraseasonal variation of 30-60 d and 120 d,while the zonal current at the C1 mooring exhibited a notable period of approximately 180 d.Further analysis indicated that the variability of currents in the Andaman Sea is influenced primarily by equatorial Kelvin waves and Rossby wave packets.Moreover,our results suggest that equatorial Kelvin waves from the eastern Indian Ocean entered the Andaman Sea in the form of Wyrtki Jets and propagated primarily along two distinct pathways during the observation period.In addition to coastal boundary Kelvin waves,it was found that a branch of the Wyrtki Jet that directly enters the Andaman Sea and flows northward along the slope of the continental shelf,and reflected Rossby wave packets by topography.

Observation of Low-Level Jets in the Eastern Tropical Indian Ocean Based on Shipborne Coherent Doppler Lidar

In contrast to the Pacific and Atlantic Oceans,the Indian Ocean has lacked in-situ observations of wind profiles over open sea areas for decades.In 2021,a shipborne coherent Doppler lidar(CDL)was used to observe in-situ wind profiles in the eastern tropical Indian Ocean.This equipment successfully captured low-level jets(LLJs)in the region,and their characteristics were thoroughly analyzed.Results reveal that the observed wind speed of LLJs in the eastern Indian Ocean ranges from 6 m s^(-1) to 10 m s^(-1) during the boreal winter and spring seasons,showing a height range of 0.6 to 1 km and two peak times at 0800 and 2000 UTC.This wind shear is weaker than that in land or offshore areas,ranging from 0 s^(-1) to 0.006 s^(-1).Moreover,the accuracy of the CDL data is compared to that of ERA5 data in the study area.The results indicate that the zonal wind from ERA5 data significantly deviated from the CDL measurement data,and the overall ERA5 data are substantially weaker than the in-situ observations.Notably,ERA5 underestimates northwestward LLJs.

Observation of Arctic surface currents using data from a surface drifting buoy

During the 10th Chinese Arctic scientific expedition carried out in the summer of 2019,the surface current in the high-latitude areas of the Arctic Ocean was observed using a self-developed surface drifting buoy,which was initially deployed in the Chukchi Sea.The buoy traversed the Chukchi Sea,Chukchi Abyssal Plain,Mendeleev Ridge,Makarov Basin,and Canada Basin over a period of 632 d.After returning to the Mendeleev Ridge,it continued to drift toward the pole.Overall,the track of the buoy reflected the characteristics of the transpolar drift and Chukchi Slope Current,as well as the inertial flow,cross-ridge surface flow,and even the surface disorganized flow for some time intervals.The results showed that:(1)the transpolar drift mainly occurs in the Chukchi Abyssal Plain,Mendeleev Ridge,and western Canada Basin to the east of the ridge where sea ice concentration is high,and the average northward flow velocity in the region between 79.41°N and 86.32°N was 5.1 cm/s;(2)the average surface velocity of the Chukchi Slope Current was 13.5 cm/s,and while this current moves westward along the continental slope,it also extends northwestward across the continental slope and flows to the deep sea;and(3)when sea ice concentration was less than 50%,the inertial flow was more significant(the maximum observed inertial flow was 26 cm/s,and the radius of the inertia circle was 3.6 km).

Simulation and analysis on seasonal variability of average salinity in the Yellow Sea

The CTD (conductivity, temperature and depth) data collected by six China-Korea joint cruises during 1996-1998 and the climatological data suggest that the seasonal variability of average salinity in the Yellow Sea (Sa) presents a general sinusoid pattern. To study the mechanism of the variability, annual cycles of Sa were simulated and a theoretical analysis based on the governing equations was reported.Three main factors are responsible for the variability: the Yellow Sea Warm Current (YSWC), the Changji-ang (Yangtze) River diluted water (YRDW) and the evaporation minus precipitation (E-P). From December to the next May, the variability of Sa is mainly controlled by the salt transportation of the YSWC. But in early July, the YSWC is overtaken and replaced by the YRDW which then becomes the most important controller in summer. From late September to November, the E-P gradually took the lead. The mass exchange north of the 37癗 line is not significant.

Retrieval of inherent optical properties of the Yellow Sea and East China Sea using a quasi-analytical algorithm

We tested and modified the quasi-analytical algorithm （QAA） using 57 groups of field data collected in the spring of 2003 in the Yellow Sea and East China Sea. The QAA performs well in deriving total absorption coefficients of typical coastal waters. The average percentage difference （APD） is in a range of 13.9%-38.5% for the total absorption coefficient （13.9% at 440 nm）, and differences in particle backscattering coefficient bbp（2） are less than 50% （in the case of the updated QAA）. To obtain improved results, we modified the QAA by adjusting the empirical relationships. The modified algorithm is then applied to the field data to test its performance. The APDs were 44.7%-46.6% for bbp（λ） and 9.9%-32.8% （9.9% at 555 nm） for the total absorption coefficient. This indicates that the modified QAA derives better results. We also used the modified model to derive phytoplankton pigment absorption （aph） and detritus and CDOM absorption （aug） coefficients. The APDs for aph and a dg at 440 nm are 37.1% and 19.8%. In this paper, we discuss error sources using the measured dataset. More independent field data can improve this algorithm and derive better results.

Development of a fine-resolution atmosphere-wave-ocean coupled forecasting model for the South China Sea and its adjacent seas

A 72-h fine-resolution atmosphere-wave-ocean coupled forecasting system was developed for the South China Sea and its adjacent seas. The forecasting model domain covers from from 15°S to 45°N in latitude and 99°E to135°E in longitude including the Bohai Sea, the Yellow Sea, the East China Sea, the South China Sea and the Indonesian seas. To get precise initial conditions for the coupled forecasting model, the forecasting system conducts a 24-h hindcast simulation with data assimilation before forecasting. The Ensemble Adjustment Kalman Filter(EAKF) data assimilation method was adopted for the wave model MASNUM with assimilating Jason-2 significant wave height(SWH) data. The EAKF data assimilation method was also introduced to the ROMS model with assimilating sea surface temperature(SST), mean absolute dynamic topography(MADT) and Argo profiles data. To improve simulation of the structure of temperature and salinity, the vertical mixing scheme of the ocean model was improved by considering the surface wave induced vertical mixing and internal wave induced vertical mixing. The wave and current models were integrated from January 2014 to October 2015 driven by the ECMWF reanalysis 6 hourly mean dataset with data assimilation. Then the coupled atmosphere-wave-ocean forecasting system was carried out 14 months operational running since November 2015. The forecasting outputs include atmospheric forecast products, wave forecast products and ocean forecast products. A series of observation data are used to evaluate the coupled forecasting results, including the wind, SHW, ocean temperature and velocity.The forecasting results are in good agreement with observation data. The prediction practice for more than one year indicates that the coupled forecasting system performs stably and predict relatively accurate, which can support the shipping safety, the fisheries and the oil exploitation.

Causes for different spatial distributions of minimum Arctic sea-ice extent in 2007 and 2012

Satellite records show the minimum Arctic sea ice extents （SIEs） were observed in the Septembers of 2007 and 2012, but the spatial distributions of sea ice concentration reduction in these two years were quite different. Atmospheric circulation pattern and the upper-ocean state in summer were investigated to explain the difference. By employing the ice-temperature and ice-specific humidity （SH） positive feedbacks in the Arctic Ocean, this paper shows that in 2007 and 2012 the higher surface air temperature （SAT） and sea level pressure （SLP） accompanied by more surface SH and higher sea surface temperature （SST）, as a consequence, the strengthened poleward wind was favorable for melting summer Arctic sea ice in different regions in these two years. SAT was the dominant factor influencing the distribution of Arctic sea ice melting. The correlation coefficient is -0.84 between SAT anomalies in summer and the Arctic SIE anomalies in autumn. The increase SAT in different regions in the summers of 2007 and 2012 corresponded to a quicker melting of sea ice in the Arctic. The SLP and related wind were promoting factors connected with SAT. Strengthening poleward winds brought warm moist air to the Arctic and accelerated the melting of sea ice in different regions in the summers of 2007 and 2012. Associated with the rising air temperature, the higher surface SH and SST also played a positive role in reducing summer Arctic sea ice in different regions in these two years, which form two positive feedbacks mechanism.

The improvement of the one-dimensional Mellor-Yamada and K-profile parameterization turbulence schemes with the non-breaking surface wave-induced vertical mixing

Both the level 2.5 Mellor-Yamada turbulence closure scheme （MY） and K-profile parameterization （KPP） are popularly used by the ocean modeling community. The MY and the KPP are improved through including the non-breaking surface wave-induced vertical mixing （Bv）, and the improved schemes were tested by using continuous data at the Papa ocean weather station （OWS） during 1961-1965. The numerical results showed that the Bv can make the temperature simulations fit much better with the continuous data from Papa Sta- tion. The two improved schemes overcame the shortcomings of predicting too shallow upper mixed layer depth and consequently overheated sea surface temperature during summertime, which are in fact com- mon problems for all turbulence closure models. Statistical analysis showed that the Bv effectively reduced the mean absolute error and root mean square error of the upper layer temperature and increased the corre- lation coefficient between simulation and the observation. Furthermore, the performance of vertical mixing induced by shear instability and the Bv is also compared. Both the temperature structure and its seasonal cycle significantly improved by including the Bv, regardless of whether shear instability was included or not, especially for the KPP mixing scheme, which suggested that Bv played a dominant role in the upper ocean where the mean current was relatively weak, such as at Papa Station. These results may provide a clue to improve ocean circulation models.

Historical simulation and twenty-first century prediction of oceanic CO_2 sink and pH change

A global ocean carbon cycle model based on the ocean general circulation model POP and the improved biogeochemical model OCMIP-2 is employed to simulate carbon cycle processes under the historically observed atmospheric CO 2 concentration and different future scenarios （called Rep- resentative Concentration Pathways, or RCPs）. The RCPs in this paper follow the design of Inter- governmental Panel on Climate Change （IPCC） for the Fifth Assessment Report （AR5）. The model results show that the ocean absorbs CO 2 from atmosphere and the absorbability will continue in the 21st century under the four RCPs. The net air-sea CO 2 flux increased during the historical time and reached 1.87 Pg/a （calculated by carbon） in 2005; however, it would reach peak and then decrease in the 21st century. The ocean absorbs CO 2 mainly in the mid latitude, and releases CO 2 in the equator area. However, in the Antarctic Circumpolar Current （ACC） area the ocean would change from source to sink under the rising CO 2 concentration, including RCP4.5, RCP6.0, and RCP8.5. In 2100, the anthropogenic carbon would be transported to the 40 S in the Atlantic Ocean by the North Atlantic Deep Water （NADW）, and also be transported to the north by the Antarctic Bottom Water （AABW） along the Antarctic continent in the Atlantic and Pacific oceans. The ocean pH value is also simulated by the model. The pH decreased by 0.1 after the industrial revolution, and would continue to decrease in the 21st century. For the highest concentration sce- nario of RCP8.5, the global averaged pH would decrease by 0.43 to reach 7.73 due to the absorption of CO 2 from atmosphere.

Water characteristics of abyssal and hadal zones in the southern Yap Trench observed with the submersible Jiaolong

Using observations in an applied cruise of the submersible Jiaolong,water characteristics,geostrophic transport,and turbulent mixing in abyssal and hadal zones of the southern Yap Trench were studied.The spatial structures of deep water show that the abyssal water is cold,saline,and oxygen rich.The hadal water has very small changes in potential temperature and potential density,and a little decrease in salinity and obvious decrease in oxygen.The isotherm,isopycnal,and isohaline are depressed in abyss over the central trench.The turbulent mixing is enhanced in the near-bottom zone and the hadal water on the trench slope,especially at the steep slope,the dissipation rate and diffusivity is strong,which weakens the stratification.The geostrophic flows move southward in the western region of the trench and northward in the eastern region,indicating cyclonic circulation.In the central region of the trench,the water transport is^1.74 Sv southward.In the hadal zone,the northward and southward transports are balanced.Our analysis suggests that the abyssal water in the southern Yap Trench is from Lower Circumpolar Water(LCPW)and the hadal water seems to be of the isolated local water rather than LCPW.

A model for the transportation and distribution of jellyfish Rhopilema esculentum for stock enhancement in the Liaodong Bay, China

A numerical model for jellyfish Rhopilema esculentum stock enhancement is developed for the first time. The model is based on an operational ocean circulation-surface wave coupled forecasting system for the seas off China and adjacent areas and uses a Lagrangian particle-tracking scheme to track the trajectories of released jellyfish. The Jellyfish are modeled as particles with diel vertical migration and are passively drifted by the current and dispersion due to the sub-grid processes. A comparison between the simulation and survey results demonstrate that the model can capture the primary distribution patterns of the released jellyfish. The model results show that the ocean current and indirect wind impact are the main drivers controlling the jellyfish transport. A connectivity matrix between the release sites and fishing grounds indicates the top of the bay is better than the eastern and western coasts for jellyfish fishing. The matrix also shows that only 45% and 27% of the jellyfish released from Wafangdian(WFD) can enter the fishing ground in 2008 and 2010; thus, the site near WFD is not an advisable location for jellyfish release. A Lagrangian probability density function based on a nine-year tracing experiment validates the results and further provides a "climatology" distribution of the released jellyfish.Several experiments are conducted to examine the sensitivity of the model to random walk schemes and to release conditions. The model requires a random walk but is insensitive to the random walk scheme. The experiments with different habitat depths show that if the jellyfish are fixed on the bottom of the water, most of them will be transported to the center, or even out of the bay, by the bottom circulation.

Ocean current observation and spectrum analysis in central Chukchi Sea during the summer of 2008

During the summer of 2008, the third CHINARE Arctic Expedition was carried out on board of Xuelong Icebreaker in the central Chukchi Sea. A submersible mooring system was deployed and recovered at Station CN-01 （71.40.024＇N, 167.58.910＇W） with 33 days of the current profile records, and continuous observation of temperature and salinity data were collected. This mooring station locates in the blank of similar observation area and it is the first time for our Chinese to finish this kind of long-term mooring work in this area. This mooring system finished integrated hydrological observations with long-term continuous record of the whole profile velocity for the first time. Based on time series analysis of temperature, salinity, velocity and flow direction, we get the following main results. （1） During the observation period, the mean surface current velocity is 70.2 cm/s eastward, and velocity reaches its maximum in average at 3 m level with magnitude 90.0 cm/s, direction 206.. （2） In 9-30 m layers, the semidiurnal period variation is the most obvious, the flow direction is quite stable, and the flow is synchronous and consistent vertically. （3） Besides the semidiurnal period variation, the main variation in the upper layer is in 11-d period, with variations in period 5.5, 5.5, and 3.7 d, which reflect the influences of sea surface wind change and maintenance. （4） Near the bottom the temperature change is correlated and synchronized with the conductivity.