Uncertainty Representation of Ocean Fronts Based on Fuzzy-Rough Set Theory

Analysis of ocean fronts＇ uncertainties indicates that they result from indiscemibility of their spatial position and fuzziness of their intensity. In view of this, a flow hierarchy for uncertainty representation of ocean fronts is proposed on the basis of fuzzy-rough set theory. Firstly, raster scanning and blurring are carried out on an ocean front, and the upper and lower approximate sets, the indiscernible relation in fuzzy-rough theories and related operators in fuzzy set theories are adopted to represent its uncertainties, then they are classified into three sets： with members one hundred percent belonging to the ocean front, belonging to the ocean front＇s edge and definitely not belonging to the ocean front. Finally, the approximate precision and roughness degree are utilized to evaluate the ocean front＇s degree of uncertainties and the precision of the representation. It has been proven that the method is not only capable of representing ocean fronts＇ uncertainties, but also provides a new theory and method for uncertainty representation of other oceanic phenomena.

Deep learning to estimate ocean subsurface salinity structure in the Indian Ocean using satellite observations

Accurately estimating the ocean subsurface salinity structure(OSSS)is crucial for understanding ocean dynamics and predicting climate variations.We present a convolutional neural network(CNN)model to estimate the OSSS in the Indian Ocean using satellite data and Argo observations.We evaluated the performance of the CNN model in terms of its vertical and spatial distribution,as well as seasonal variation of OSSS estimation.Results demonstrate that the CNN model accurately estimates the most significant salinity features in the Indian Ocean using sea surface data with no significant differences from Argo-derived OSSS.However,the estimation accuracy of the CNN model varies with depth,with the most challenging depth being approximately 70 m,corresponding to the halocline layer.Validations of the CNN model’s accuracy in estimating OSSS in the Indian Ocean are also conducted by comparing Argo observations and CNN model estimations along two selected sections and four selected boxes.The results show that the CNN model effectively captures the seasonal variability of salinity,demonstrating its high performance in salinity estimation using sea surface data.Our analysis reveals that sea surface salinity has the strongest correlation with OSSS in shallow layers,while sea surface height anomaly plays a more significant role in deeper layers.These preliminary results provide valuable insights into the feasibility of estimating OSSS using satellite observations and have implications for studying upper ocean dynamics using machine learning techniques.

New Record Ocean Temperatures and Related Climate Indicators in 2023

The global physical and biogeochemical environment has been substantially altered in response to increased atmospheric greenhouse gases from human activities.In 2023,the sea surface temperature(SST)and upper 2000 m ocean heat content(OHC)reached record highs.The 0–2000 m OHC in 2023 exceeded that of 2022 by 15±10 ZJ(1 Zetta Joules=1021 Joules)(updated IAP/CAS data);9±5 ZJ(NCEI/NOAA data).The Tropical Atlantic Ocean,the Mediterranean Sea,and southern oceans recorded their highest OHC observed since the 1950s.Associated with the onset of a strong El Niño,the global SST reached its record high in 2023 with an annual mean of~0.23℃ higher than 2022 and an astounding>0.3℃ above 2022 values for the second half of 2023.The density stratification and spatial temperature inhomogeneity indexes reached their highest values in 2023.

Alleviated photoinhibition on nitrification in the Indian Sector of the Southern Ocean

Nitrification,a central process in the marine nitrogen cycle,produces regenerated nitrate in the euphotic zone and emits N_(2)O,a potent greenhouse gas as a by-product.The regulatory mechanisms of nitrification in the Southern Ocean,which is a critical region for CO_(2)sequestration and radiative benefits,remain poorly understood.Here,we investigated the in situ and dark nitrification rates in the upper 500 m and conducted substrate kinetics experiments across the Indian Sector in the Cosmonaut and Cooperation seas in the late austral summer.Our findings indicate that light inhibition of nitrification decreases exponentially with depth,exhibiting a light threshold of 0.53%photosynthetically active radiation.A positive relationship between dark nitrification and apparent oxygen utilization suggests a dependence on substrate availability from primary production.Importantly,an increased NH_(4)^(+) supply can act as a buffer against photo-inhibitory damage.Globally,substrate affinity(α)increases with depth and transitions from light to dark,decreases with increasing ambient NH_(4)^(+)and exhibits a latitudinal distribution,reflecting substrate utilization strategies.We also reveal that upwelling in Circumpolar Deep Water(CDW)stimulates nitrification through the introduction of potentially higher iron and deep diverse nitrifying microorganisms with higherα.We conclude that although light is the primary limiting factor for nitrification in summer,coupling between substrate availability and CDW upwelling can overcome this limitation,thereby alleviating photoinhibition by up to 45%±5.3%.

The Influence of the Closure of the East Paleo-Tethys Ocean on Southern South China:Evidences from Kinematics and^(40)Ar/^(39)Ar Geochronology of the Rongxian Ductile Shear Zone in Southeastern Guangxi

The Triassic was a crucial period in the tectonic evolution of the South China Block.Research on tectonic deformation during this period provides information on intracontinental orogenic mechanisms in South China.In this study,alongside thermochronological analyses,we examine the macroscopic and microscopic structural features of the Rongxian ductile shear zone,located south of the Darongshan granite in the southeastern part of Guangxi Province,on the southern margin of South China.Sinistral shear is indicated by the characteristics of rotatedσ-type feldspar porphyroclasts,stretching lineations defined by elongated quartz grains and the orientations of quartz c-axes.LA-ICP-MS U-Pb dating of zircons from two samples of granitic mylonite and one of granite yielded ages of ca.256 Ma.Furthermore,two samples of granitic mylonite yield muscovite^(40)Ar/^(39)Ar plateau ages of 249-246 Ma.These results indicate that the Rongxian ductile shear zone resulted from Early Triassic deformation of the late Permian Darongshan granite.This deformation was likely related to the closure of the eastern Paleo-Tethys Ocean and the subsequent collision of the South China and Indochina blocks,during the early stage of the Indosinian orogeny.

Networking Observation and Applications of Chinese Ocean Satellites

This paper presents the networking observation capabilities of Chinese ocean satellites and their diverse applications in ocean disaster prevention,ecological monitoring,and resource development.Since the inaugural launch in 2002,China has achieved substantial advancements in ocean satellite technology,forming an observation system composed of the HY-1,HY-2,and HY-3 series satellites.These satellites are integral to global ocean environmental monitoring due to their high resolution,extensive coverage,and frequent observations.Looking forward,China aims to further enhance and expand its ocean satellite capabilities through ongoing projects to support global environmental protection and sustainable development.

Microearthquake reveals the lithospheric structure at midocean ridges and oceanic transform faults

Mid-ocean ridge and oceanic transforms are among the most prominent features on the seafloor surface and are crucial for understanding seafloor spreading and plate tectonic dynamics,but the deep structure of the oceanic lithosphere remains poorly understood.The large number of microearthquakes occurring along ridges and transforms provide valuable information for gaining an indepth view of the underlying detailed seismic structures,contributing to understanding geodynamic processes within the oceanic lithosphere.Previous studies have indicated that the maximum depth of microseismicity is controlled by the 600-℃isotherm.However,this perspective is being challenged due to increasing observations of deep earthquakes that far exceed this suggested isotherm along mid-ocean ridges and oceanic transform faults.Several mechanisms have been proposed to explain these deep events,and we suggest that local geodynamic processes(e.g.,magma supply,mylonite shear zone,longlived faults,hydrothermal vents,etc.)likely play a more important role than previously thought.

Impacts of Ice-Ocean Stress on the Subpolar Southern Ocean:Role of the Ocean Surface Current

The mechanical influences involved in the interaction between the Antarctic sea ice and ocean surface current(OSC)on the subpolar Southern Ocean have been systematically investigated for the first time by conducting two simulations that include and exclude the OSC in the calculation of the ice-ocean stress(IOS), using an eddy-permitting coupled ocean-sea ice global model. By comparing the results of these two experiments, significant increases of 5%, 27%, and 24%, were found in the subpolar Southern Ocean when excluding the OSC in the IOS calculation for the ocean surface stress,upwelling, and downwelling, respectively. Excluding the OSC in the IOS calculation also visibly strengthens the total mechanical energy input to the OSC by about 16%, and increases the eddy kinetic energy and mean kinetic energy by about38% and 12%, respectively. Moreover, the response of the meridional overturning circulation in the Southern Ocean yields respective increases of about 16% and 15% for the upper and lower branches;and the subpolar gyres are also found to considerably intensify, by about 12%, 11%, and 11% in the Weddell Gyre, the Ross Gyre, and the Australian-Antarctic Gyre, respectively. The strengthened ocean circulations and Ekman pumping result in a warmer sea surface temperature(SST), and hence an incremental surface heat loss. The increased sea ice drift and warm SST lead to an expansion of the sea ice area and a reduction of sea ice volume. These results emphasize the importance of OSCs in the air-sea-ice interactions on the global ocean circulations and the mass balance of Antarctic ice shelves, and this component may become more significant as the rapid change of Antarctic sea ice.

Quantifying the Role of the Eddy Transfer Coefficient in Simulating the Response of the Southern Ocean Meridional Overturning Circulation to Enhanced Westerlies in a Coarse-resolution Model

This study assesses the capability of a coarse-resolution ocean model to replicate the response of the Southern Ocean Meridional Overturning Circulation(MOC) to intensified westerlies,focusing on the role of the eddy transfer coefficient(κ).κ is a parameter commonly used to represent the velocities induced by unresolved eddies.Our findings reveal that a stratification-dependent κ,incorporating spatiotemporal variability,leads to the most robust eddy-induced MOC response,capturing 82% of the reference eddy-resolving simulation.Decomposing the eddy-induced velocity into its vertical variation(VV) and spatial structure(SS) components unveils that the enhanced eddy compensation response primarily stems from an augmented SS term,while the introduced VV term weakens the response.Furthermore,the temporal variability of the stratification-dependent κ emerges as a key factor in enhancing the eddy compensation response to intensified westerlies.The experiment with stratification-dependent κ exhibits a more potent eddy compensation response compared to the constant κ,attributed to the structure of κ and the vertical variation of the density slope.These results underscore the critical role of accurately representing κ in capturing the response of the Southern Ocean MOC and emphasize the significance of the isopycnal slope in modulating the eddy compensation mechanism.

Hydrographic Characteristics and Oceanic Heat Flux in the Upper Arctic Ocean over the Alpha Ridge Observed by the DTOP Platform in 2018 and 2021

In 2018 and 2021,the Drift-Towing Ocean Profilers(DTOP)provided extensive temperature and salinity data on the upper 120m ocean through their drifts over the Alpha Ridge north of the Canada Basin.The thickness and temperature maximum of Alaska Coastal Water(ACW)ranged from 20m to 40m and-1.5℃to-0.8℃,respectively,and the salinity generally maintained from 30.2 to 32.5.Comparison with World Ocean Atlas 2018’s climatology manifested a 40m-thick and warm ACW roughly ex-ceeding the temperature maximum by 0.4–0.5℃in June–August 2021.This anomalously warm ACW was highly related to the ex-pansion of the Beaufort Gyre in the negative Arctic Oscillation phase.During summer,the under-ice oceanic heat flux F_(w)^(OHF)was elevated,with a maximum value of above 25Wm^(-2).F_(w)^(OHF)was typically low in the freezing season,with an average value of 1.2Wm^(-2).The estimates of upward heat flux contributed by ACW to the sea ice bottom F_(w)^(OHF)were in the range of 3–4Wm^(-2)in June–August 2021,when ACW contained a heat content of more than 80MJm^(-2).The heat loss over this period was driven by a weak stratification upon the ACW layer associated with a surface mixed layer(SML)approaching the ACW core.After autumn,F_(w)^(OHF)was reduced(<2 Wm^(-2))except during rare events when it elevated F_(w)^(OHF)slightly.In addition,the intensive and widespread Ekman suction,which created a violent upwelling north of the Canada Basin,was largely responsible for the substantial cooling and thinning of the ACW layer in the summer of 2021.

Contrasting physical mechanisms of yellowfin tuna fluctuations between the western and eastern Indian Ocean

As an economically critical pelagic migratory species,yellowfin tuna(Thunnus albacores,YFT)is very sensible to physical and environmental conditions,such as sea surface temperature(SST),ocean heat content(OHC),and the mixed layer depth(MLD).We investigated the impact of SST,OHC,and MLD on fluctuations of YFT catch in the western/eastern Indian Ocean using the long time series of 63-year environmental and YFT datasets.We found that the impact of SST on YFT was heavily overestimated in the past,and MLD plays a more critical role in the YFT catch fluctuation.When the MLD deepens(>34.8 m),SST was more influential in predicting the catches of YFT than OHC in the western Indian Ocean,and OHC was more critical to YFT than SST in the eastern Indian Ocean.However,when the MLD shallows(<34.8 m),MLD was more vital to predict the catch per unit effort(CPUE)of YFT than SST/OHC in the western.After 2000,there was an asynchronous pattern of YFT CPUE induced by higher frequency variations and ocean hiatus of SST/OHC signals in the western and eastern Indian Oceans basins.The impact of the subsurface hiatus may induce the decrease of YFT in the eastern Indian Ocean.The above findings clarified a non-stationary relationship between the environmental factors and catches of YFT and provided new insights into variations in YFT abundance.

The sudden ocean warming and its potential influences on earlyfrozen landfast ice in the Prydz Bay, East Antarctica

The ocean conditions beneath the ice cover play a key role in understanding the sea ice mass balance in the polar regions.An integrated high-frequency ice-ocean observation system,including Acoustic Doppler Velocimeter,Conductivity-Temperature-Depth Sensor,and Sea Ice Mass Balance Array(SIMBA),was deployed in the landfast ice region close to the Chinese Zhongshan Station in Antarctica.A sudden ocean warming of 0.14℃(p<0.01)was observed beneath early-frozen landfast ice,from(−1.60±0.03)℃during April 16-19 to(−1.46±0.07)℃during April 20-23,2021,which is the only significant warming event in the nearly 8-month records.The sudden ocean warming brought a double rise in oceanic heat flux,from(21.7±11.1)W/m^(2) during April 16-19 to(44.8±21.3)W/m^(2) during April 20-23,2021,which shifted the original growth phase at the ice bottom,leading to a 2 cm melting,as shown from SIMBA and borehole observations.Simultaneously,the slowdown of ice bottom freezing decreased salt rejection,and the daily trend of observed ocean salinity changed from+0.02 d^(-1) during April 16-19,2021 to+0.003 d^(-1) during April 20-23,2021.The potential reasons are increased air temperature due to the transit cyclones and the weakened vertical ocean mixing due to the tide phase transformation from semi-diurnal to diurnal.The high-frequency observations within the ice-ocean boundary layer enhance the comprehensive investigation of the ocean’s influence on ice evolution at a daily scale.

Distribution and influencing factors of microeukaryote in different water layers of the southwestern Indian Ocean Ridge

Microeukaryotes play a vital role in shaping marine ecosystems,especially in marine productivity,the microbial food web,and carbon cycle.The Indian Ocean is one of the largest oligotrophic areas in the world,but little is known about the biodiversity of microeukaryotes in the area.The community composition and geographical distribution of microeukaryotes collected from the surface(SUR)and deep chlorophyll maximum(DCM)layers in the southwestern Indian Ocean were studied using high-throughput sequencing of the 18S rRNA gene.The metagenomic data helped quantify the impact of environmental factors on microeukaryotic communities.The relative abundance of different taxa groups exhibited distinct patterns between SUR and DCM layers,except for the most dominant Dinoflagellata that accounted for more than 40.6%abundance in each sample.Radiolaria was much more abundant in the nutrient-rich DCM layer than the SUR layer.The community similarity of microeukaryotes decreased with increasing of geographic distance,whereas the temperature and inorganic nitrogen were the most important environmental parameters to community structure.Abundant communities were more influenced by dispersal limitations and rare communities were more responsive to environmental factors.Correlation network analyses revealed strong biotic interactions indicative of parasitism,predation and competition,and their contribution to microeukaryotic population in diverse environments.Overall,this study provided insights into the biodiversity of microeukaryotes by characterizing the differences between water layers and identifying the driving factors in the ocean.

Geochronology and Petrogenesis of Late Carboniferous to Early Permian Basalts in the Central Lhasa Subterrane, Southern Tibet: Implications for the Evolution of the Sumdo Paleo-Tethys Ocean

Basalts from the Late Carboniferous to Early Permian are extensively developed in the central Lhasa subterrane, southern Tibet. Studying the petrogenesis of these rocks may have implications for the late Paleozoic arc magmatism along the central Lhasa subterrane uncovering more of the evolution of the Sumdo Paleo-Tethys Ocean and its dynamic mechanism. Basalt samples from the Luobadui Formation in the Leqingla area, NW of Linzhou City in the central Lhasa subterrane, southern Tibet exhibit arc-like geochemical signatures in a subduction-zone tectonic setting characterized by high Al_(2)O_(3) and low TiO_(2) contents, fractionated REE patterns with low Nb/La ratios and high LREE concentrations, and negative HFSE anomalies. Based on their higher Th/Ce, Nb/Zr, and lower Ba/Th, Pb/Nd ratios, slightly negative to positive ε_(Nd)(t) values, and the relatively high Sr-Pb isotopic compositions, these samples were probably derived from partial melting of a depleted mantle source of garnet + spinel lherzolite, metasomatized by subducted sediments around 297 Ma. Modeling of the trace elements indicates that these basalts experienced fractional crystallization of olivine, clinopyroxene and minor plagioclase during magma ascent and eruption. It is proposed that these Late Carboniferous–Early Permian basalts are associated with the northward subduction of the Sumdo Paleo-Tethys Ocean seafloor along the southern margin of the central Lhasa subterrane.

Effects of Ocean Acidification on Nitrogen Metabolism of Skeletonema costatum

Ocean acidification(OA),caused by the rising concentration of atmospheric CO_(2),leads to changes in the marine carbonate system.This,in turn,affects the physiological processes of phytoplankton.In response to increased pCO_(2) levels,marine microalgae modulate their physiological responses to meet their energy and metabolic requirements.Nitrogen metabolism is a critical metabolic pathway,directly affecting the growth and reproductive capacity of marine microorganisms.Understanding the molecular mechanisms that regulate nitrogen metabolism in microalgae under OA conditions is therefore crucial.This study aimed to investi-gate how OA affects the expression profiles of key genes in the nitrogen metabolic pathway of the marine diatom Skeletonema costatum.Our findings indicate that OA upregulates key genes involved in the nitrogen metabolic pathway,specifically those related to nitrate assimilation and glutamate metabolism.Moreover,pCO_(2) has been identified as the predominant factor affecting the expression of these genes,with a more significant impact than pH variations in S.costatum.This research not only advances our understanding of the adaptive mechanisms of S.costatum in response to OA but also provides essential data for predicting the ecological consequences of OA on marine diatoms.

Investigating the Impact of Ocean Acidification on Anti-Stress Mechanisms in Sepia esculenta Larvae Based on Transcriptome Profiling

With the rapid development of oil,energy,power and other industries,CO_(2) emissions rise sharply,which will cause a large amount of CO_(2) in the air be absorbed by the ocean and lead to ocean acidification.The growth and development of organisms can be seriously affected by acidified seawater.Sepia esculenta is a mollusk with high nutritional and economic value and is widely cultured in offshore waters of China.Larvae are the early life forms of the organism and are more vulnerable to changes in the external environment.Too low pH will lead to some adverse reactions in larvae,which will affect metabolism,immune response and other life activities.In this study,we sequenced the transcriptome of S.esculenta subjected to acidified seawater stress and identified 1072differentially expressed genes(DEGs).The detected atypical expression of DEGs substantiates cellular malformation and translocation in S.esculenta under low pH stimulation.Simultaneously,this also substantiates the notable impact of ocean acidification on mollusks.These DEGs were used for functional enrichment analysis of GO and KEGG,and the top twenty items of the biological process classification in GO terms and 11 KEGG signaling pathways were significantly enriched.Finally,the constructed proteinprotein interaction network(PPI)was used to analyze protein-protein interactions,and 12 key DEGs and 3 hub genes were identified.The reliability of 12 genes was verified by quantitative RT-PCR.A comprehensive analysis of the KEGG signaling pathway and PPI revealed that ocean acidification leads to abnormalities in lipid metabolism in S.esculenta larvae,which can lead to cancer development and metastasis,accompanied by some degree of inflammation.The results of the study will help to further investigate the physiological processes of S.esculenta when stimulated by ocean acidification,and provide a reference to cope with the captive breeding of S.esculenta affected by acidification.

Environmental characteristics of trace metals in seawater from the Ninety East Ridge in the Indian Ocean

The Ninety East Ridge in the Indian Ocean has complex and unique characteristics.The concentrations and distribution characteristics of 10 trace metals(V,Cr,Mn,Fe,Co,Ni,Cu,Cd,Pb,and U)in seawater from the Ninety East Ridge in the Indian Ocean were investigated.Results show that the average concentrations of different trace metals in all the collected seawater samples were 1.134μg/L for V,0.158μg/L for Cr,0.489μg/L for Mn,0.427μg/L for Fe,0.011μg/L for Co,0.395μg/L for Ni,0.403μg/L for Cu,0.097μg/L for Cd,0.139μg/L for Pb,and 3.470μg/L for U.Differences in the horizontal and vertical distributions of all measured trace metals were revealed,and the occurrence of high concentrations was nonuniform.In addition,the significant differences in the concentration distribution of different trace metals in seawater on both sides of the Ninety East Ridge present regional segmentation in the area for various trace metals in deep sea water.This study provided basic data for future investigations on the environmental and ecological impact of trace metals in the Indian Ocean and the potential water mass transport mechanism.

Early Triassic Silicic Volcanics of South China:Petrogenesis and Constraints on the Geodynamic Evolution of the Paleo-Tethys Ocean Region

The only occurrence of Lower Triassic silicic volcanic rocks within the South China Block is in the Qinzhou Bay area of Guangxi Province.LA-ICP-MS zircon U-Pb dating reveals that volcanic rocks of the Beisi and Banba formations formed between 248.8±1.6 and 246.5±1.3 Ma,coeval with peraluminous granites of the Qinzhou Bay Granitic Complex.The studied rhyolites and dacites are characterized by high SiO_(2),K_(2)O,and Al_(2)O_(3),and low MgO,CaO,and P_(2)O_(5) contents and are classified as high-K calc-alkaline S-type rocks,with A/CNK=0.98-1.19.The volcanic rocks are depleted in high field strength elements,e.g.,Nb,Ta,Ti,and P,and enriched in large ion lithophile elements,e.g.,Rb,K,Sr,and Ba.Although the analyzed volcanic rocks have extremely enriched zircon Hf isotopic compositions(ε_(Hf)(t)=-29.1 to-6.9),source discrimination indicators and high calculated Ti-in-zircon temperatures(798-835℃)reveal that magma derived from enriched lithospheric mantle not only provided a heat source for anatectic melting of the metasedimentary protoliths but was also an endmember component of the S-type silicic magma.The studied early Triassic volcanics are inferred to have formed immediately before closure of the Paleo-Tethys Ocean in this region,as the associated subduction would have generated an extensional setting in which the mantle-derived upwelling and volcanic activity occurred.

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.

Assessment of the Southern Ocean Sea Surface Temperature Biases in CMIP5 and CMIP6 Models

This work evaluates the performances of climate models in simulating the Southern Ocean(SO)sea surface temperature(SST)by a large ensemble from phases 5 and 6 of the Coupled Model Intercomparison Project(CMIP5 and CMIP6).By combining models from the same community sharing highly similar SO SST biases and eliminating the effect of global-mean biases on local SST biases,the results reveal that the ensemble-mean SO SST bias at 70°-30°S decreases from 0.38℃ in CMIP5 to 0.28℃ in CMIP6,together with increased intermodel consistency.The dominant mode of the intermodel variations in the zonal-mean SST biases is characterized as a meridional uniform warm bias pattern,explaining 79.1% of the intermodel variance and exhibiting positive principal values for most models.The ocean mixed layer heat budget further demonstrates that the SST biases at 70°-50°S primarily result from the excessive summertime heating effect from surface net heat flux.The biases in surface net heat flux south of 50°S are largely impacted by surface shortwave radiation from cloud and clear sky components at different latitudes.North of 50°S,the underestimated westerlies reduce the northward Ekman transport and hence northward cold advection in models,leading to warm SST biases year-round.In addition,the westerly biases are primarily traced back to the atmosphere-alone model simulations forced by the observed SST and sea ice.These results disclose the thermal origin at the high latitude and dynamical origin at the low latitude of the SO SST biases and underscore the significance of the deficiencies of atmospheric models in producing the SO SST biases.