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.

Anthropogenic Influence on Decadal Changes in Concurrent Hot and Dry Events over China around the Mid-1990s

The frequency and duration of observed concurrent hot and dry events(HDEs) over China during the growing season(April–September) exhibit significant decadal changes across the mid-1990s. These changes are characterized by increases in HDE frequency and duration over most of China, with relatively large increases over southeastern China(SEC), northern China(NC), and northeastern China(NEC). The frequency of HDEs averaged over China in the present day(PD,1994–2011) is double that in the early period(EP, 1964–81);the duration of HDEs increases by 60%. Climate experiments with the Met Office Unified Model(MetUM-GOML2) are used to estimate the contributions of anthropogenic forcing to HDE decadal changes over China. Anthropogenic forcing changes can explain 60%–70% of the observed decadal changes,suggesting an important anthropogenic influence on HDE changes over China across the mid-1990s. Single-forcing experiments indicate that the increase in greenhouse gas(GHG) concentrations dominates the simulated decadal changes,increasing the frequency and duration of HDEs throughout China. The change in anthropogenic aerosol(AA) emissions significantly decreases the frequency and duration of HDEs over SEC and NC, but the magnitude of the decrease is much smaller than the increase induced by GHGs. The changes in HDEs in response to anthropogenic forcing are mainly due to the response of climatological mean surface air temperatures. The contributions from changes in variability and changes in climatological mean soil moisture and evapotranspiration are relatively small. The physical processes associated with the response of HDEs to GHG and AA changes are also revealed.

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.

The Role of Underlying Boundary Forcing in Shaping the Recent Decadal Change of Persistent Anomalous Activity over the Ural Mountains

Observational analyses demonstrate that the Ural persistent positive height anomaly event(PAE) experienced a decadal increase around the year 2000, exhibiting a southward displacement afterward. These decadal variations are related to a large-scale circulation shift over the Eurasian Continent. The effects of underlying sea ice and sea surface temperature(SST) anomalies on the Ural PAE and the related atmospheric circulation were explored by Atmospheric Model Intercomparison Project(AMIP) experiments from the Coupled Model Intercomparison Project Phase 6 and by sensitivity experiments using the Atmospheric General Circulation Model(AGCM). The AMIP experiment results suggest that the underlying sea ice and SST anomalies play important roles. The individual contributions of sea ice loss in the Barents-Kara Seas and the SST anomalies linked to the phase transition of the Pacific Decadal Oscillation(PDO) and Atlantic Multidecadal Oscillation(AMO) are further investigated by AGCM sensitivity experiments isolating the respective forcings.The sea ice decline in Barents-Kara Seas triggers an atmospheric wave train over the Eurasian mid-to-high latitudes with positive anomalies over the Urals, favoring the occurrence of Ural PAEs. The shift in the PDO to its negative phase triggers a wave train propagating downstream from the North Pacific. One positive anomaly lobe of the wave train is located over the Ural Mountains and increases the PAE there. The negative-to-positive transition of the AMO phase since the late-1990s causes positive 500-h Pa height anomalies south of the Ural Mountains, which promote a southward shift of Ural PAE.

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%.

Strengthened Regulation of the Onset of the South China Sea Summer Monsoon by the Northwest Indian Ocean Warming in the Past Decade

Traditionally,a delayed(early)onset of the South China Sea summer monsoon(SCSSM)has been observed to follow a warm(cold)El Niño-Southern Oscillation(ENSO)event in winter,supporting high seasonal predictability of SCSSM onset.However,the empirical seasonal forecasting skill of the SCSSM onset,solely based on ENSO,has deteriorated since 2010.Meanwhile,unexpected delayed onsets of the SCSSM have also occurred in the past decade.We attribute these changes to the Northwest Indian Ocean(NWIO)warming of the sea surface.The NWIO warming has teleconnections related to(1)suppressing the seasonal convection over the South China Sea,which weakens the impacts of ENSO on SCSSM onset and delays the start of SCSSM,and(2)favoring more high-frequency,propagating moist convective activities,which enhances the uncertainty of the seasonal prediction of SCSSM onset date.Our results yield insight into the predictability of the SCSSM onset under the context of uneven ocean warming operating within the larger-scale background state of global climate change.

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.

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.

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.

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.

Decadal Changes in Dry and Wet Heatwaves in Eastern China:Spatial Patterns and Risk Assessment

Under global warming,understanding the long-term variation in different types of heatwaves is vital for China’s preparedness against escalating heat stress.This study investigates dry and wet heatwave shifts in eastern China over recent decades.Spatial trend analysis displays pronounced warming in inland midlatitudes and the Yangtze River Valley,with increased humidity in coastal regions.EOF results indicate intensifying dry heatwaves in northern China,while the Yangtze River Valley sees more frequent dry heatwaves.On the other hand,Indochina and regions north of 25°N also experience intensified wet heatwaves,corresponding to regional humidity increases.Composite analysis is conducted based on different situations:strong,frequent dry or wet heatwaves.Strong dry heatwaves are influenced by anticyclonic circulations over northern China,accompanied by warming SST anomalies around the coastal midlatitudes of the western North Pacific(WNP).Frequent dry heatwaves are related to strong subsidence along with a strengthened subtropical high over the WNP.Strong and frequent wet heatwaves show an intensified Okhotsk high at higher latitudes in the lower troposphere,and a negative circumglobal teleconnection wave train pattern in the upper troposphere.Decaying El Niño SST patterns are observed in two kinds of wet heatwave and frequent dry heatwave years.Risk analysis indicates that El Niño events heighten the likelihood of these heatwaves in regions most at risk.As global warming continues,adapting and implementing mitigation strategies toward extreme heatwaves becomes crucial,especially for the aforementioned regions under significant heat stress.

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.

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 Dynamic Impact of Ocean on Continent

Around 71% of the Earth’s surface is covered by oceans with depths that exceed several kilometers, while continents are geographically enclosed by these vast bodies of water. The principle of fluid mechanics stipulates that water yields pressure everywhere in the container that holds it, and the water pressure against the wall of container generates force. Ocean basins are naturally gigantic containers of water, in which continents form the walls of the containers. In this study, we present that the ocean water pressure against the walls of continents generates enormous force, and determine the distribution of this force around continents and estimate its amplitude to be of the order of 1017 N per kilometer of continent width. Our modelling suggests that the stresses yielded by this force are mostly concentrated on the upper part of the continental crust, and their magnitudes reach up to 2.0 - 6.0 MPa. Our results suggest that the force may have significantly impacted the dynamics of continent (lithospheric plate) and its evolution.

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.

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.

Changes in area fraction of sediment-laden sea ice in the Arctic Ocean during 2000 to 2021

Sediment-laden sea ice plays an important role in Arctic sediment transport and biogeochemical cycles,as well as the shortwave radiation budget and melt onset of ice surface.However,at present,there is a lack of efficient observation approach from both space and in situ for the coverage of Arctic sediment-laden sea ice.Thus,both spatial distribution and long-term changes in area fraction of such ice floes are still unclear.This study proposes a new classification method to extract Arctic sediment-laden sea ice on the basic of the difference in spectral characteristics between sediment-laden sea ice and clean sea ice in the visible band using the MOD09A1 data with the resolution of 500 m,and obtains its area fraction over the pan Arctic Ocean during 2000−2021.Compared with Landsat-8 true color verification images with a resolution of 30 m,the overall accuracy of our classification method is 92.3%,and the Kappa coefficient is 0.84.The impact of clouds on the results of recognition and spatiotemporal changes of sediment-laden sea ice is relatively small from June to July,compared to that in May or August.Spatially,sediment-laden sea ice mostly appears over the marginal seas of the Arctic Ocean,especially the continental shelf of Chukchi Sea and the Siberian seas.Associated with the retreat of Arctic sea ice extent,the total area of sediment-laden sea ice in June-July also shows a significant decreasing trend of 8.99×10^(4) km^(2) per year.The occurrence of sediment-laden sea ice over the Arctic Ocean in June-July leads to the reduce of surface albedo over the ice-covered ocean by 14.1%.This study will help thoroughly understanding of the role of sediment-laden sea ice in the evolution of Arctic climate system and marine ecological environment,as well as the heat budget and mass balance of sea ice itself.

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.