Predictability of the upper ocean heat content in a Community Earth System Model ensemble prediction system

Upper ocean heat content(OHC)has been widely recognized as a crucial precursor to high-impact climate variability,especially for that being indispensable to the long-term memory of the ocean.Assessing the predictability of OHC using state-of-the-art climate models is invaluable for improving and advancing climate forecasts.Recently developed retrospective forecast experiments,based on a Community Earth System Model ensemble prediction system,offer a great opportunity to comprehensively explore OHC predictability.Our results indicate that the skill of actual OHC predictions varies across different oceans and diminishes as the lead time of prediction extends.The spatial distribution of the actual prediction skill closely resembles the corresponding persistence skill,indicating that the persistence of OHC serves as the primary predictive signal for its predictability.The decline in actual prediction skill is more pronounced in the Indian and Atlantic oceans than in the Pacific Ocean,particularly within tropical regions.Additionally,notable seasonal variations in the actual prediction skills across different oceans align well with the phase-locking features of OHC variability.The potential predictability of OHC generally surpasses the actual prediction skill at all lead times,highlighting significant room for improvement in current OHC predictions,especially for the North Indian Ocean and the Atlantic Ocean.Achieving such improvements necessitates a collaborative effort to enhance the quality of ocean observations,develop effective data assimilation methods,and reduce model bias.

Seismic monitoring of sub-seafloor fluid processes in the Haima cold seep area using an Ocean Bottom Seismometer (OBS)

The use of ocean bottom seismometers provides an effective means of studying the process and the dynamic of cold seeps by continuously recording micro-events produced by sub-seafloor fluid migration.We deployed a four-component Ocean Bottom Seismometer(OBS)at an active site of the Haima cold seep from 6 November to 19 November in 2021.Here,we present the results of this short-term OBS monitoring.We first examine the OBS record manually to distinguish(by their distinctive seismographic signatures)four types of events:shipping noises,vibrations from our remotely operated vehicle(ROV)operations,local earthquakes,and short duration events(SDEs).Only the SDEs are further discussed in this work.Such SDEs are similar to those observed in other sea areas and are interpreted to be correlated with sub-seafloor fluid migration.In the OBS data collected during the 14-day monitoring period.We identify five SDEs.Compared to the SDE occurrence rate observed in other cold seep regions,five events is rather low,from which it could be inferred that fluid migration,and subsequent gas seepage,is not very active at the Haima site.This conclusion agrees with multi-beam and chemical observations at that site.Our observations thus provide further constraint on the seepage activity in this location.This is the first time that cold seep-related SDEs have been identified in the South China Sea,expanding the list of sea areas where SDEs are now linked to cold seep fluid migration.

A combination of tides and nontidal variations in ocean bottom pressure may generate interannual slip fluctuations in the transition zone along a subduction plate interface

The tidal triggering of earthquakes has been studied for many years.The discovery of slow earthquakes in the early 2000s,including slow slip,has urged scientists to investigate the tidal responses of these earthquakes due to their sensitivity to weak stress perturbations.Previous studies have shown that slow earthquakes correlate with diurnal and semidiurnal tides and seasonal variations in surface loads more clearly than ordinary earthquakes.However,little is known about long-term responses to external stresses.In this paper,based on a widely accepted frictional law for faults,a mechanism is proposed by which nontidal variations in ocean bottom pressure,when combined with tides,pro mote the occurrence of slow earthquakes.Because slow earthquakes accompany a slip on the plate interface,this mechanism allows one to estimate slip modulations.A one-degree-of-freedom slip model is constructed and applied to Ise Bay in the Tonankai region of southwestern Japan,where large-scale ocean mass redistributions have occurred.The model calculated with parameters determined from the observation of tectonic tremors is quantitatively consistent with the slip during 1997-2013 inferred from GNSS data,suggesting that the decrease in the sea-level change in approximately 2006 could cause the accele ration of a slip observed after that.This result implies that the decreases in sea level in approximately 1996 and 2014 could also cause subsequent slip accelerations.These three slip acceleration periods temporally coincide with the increases in background seismicity in a shallower portion of the plate interface.These changes in seismicity are common to shallow earthquakes in the Tokai area,and a similar model can reproduce them.Further studies are expected to reveal causality between shallow earthquakes and long-term slip fluctuations based on modeling that considers changes in the frictional property along the plate interface.

Response of the North Pacific Storm Track Activity in the Cold Season to Multi-scale Oceanic Variations of Kuroshio Extension System: A Statistical Assessment

In this paper,a statistical method called Generalized Equilibrium Feedback Analysis(GEFA)is used to investigate the responses of the North Pacific Storm Track(NPST)in the cold season to the multi-scale oceanic variations of the Kuroshio Extension(KE)system,including its large-scale variation,oceanic front meridional shift,and mesoscale eddy activity.Results show that in the cold season from the lower to the upper troposphere,the KE large-scale variation significantly weakens the storm track activity over the central North Pacific south of 30°N.The northward shift of the KE front significantly strengthens the storm track activity over the western and central North Pacific south of 40°N,resulting in a southward shift of the NPST.In contrast,the NPST response to KE mesoscale eddy activity is not so significant and relatively shallow,which only shows some significant positive signals near the dateline in the lower and middle troposphere.Furthermore,it is found that baroclinicity and baroclinic energy conversion play an important role in the formation of the NPST response to the KE multi-scale oceanic variations.

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.

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.

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.

Sedimentary architecture of submarine channel-lobe systems under different seafloor topography:Insights from the Rovuma Basin offshore East Africa

Seafloor topography plays an important role in the evolution of submarine lobes.However,it is still not so clear how the shape of slope affects the three-dimensional(3-D)architecture of submarine lobes.In this study,we analyze the effect of topography factors on different hierarchical lobe architectures that formed during Pliocene to Quaternary in the Rovuma Basin offshore East Africa.We characterize the shape,size and growth pattern of different hierarchical lobe architectures using 3-D seismic data.We find that the relief of the topographic slope determines the location of preferential deposition of lobe complexes and single lobes.When the topography is irregular and presents topographic lows,lobe complexes first infill these depressions.Single lobes are deposited preferentially at positions with higher longitudinal(i.e.across-slope)slope gradients.As the longitudinal slope becomes higher,the aspect ratio of the single lobes increases.Lateral(i.e.along-slope)topography does not seem to have a strong influence on the shape of single lobe,but it seems to affect the overlap of single lobes.When the lateral slope gradient is relatively high,the single lobes tend to have a larger overlap surface.Furthermore,as the average of lateral slope and longitudinal slope gets greater,the width/thickness ratio of the single lobe is smaller,i.e.sediments tend to accumulate vertically.The results demonstrate that the shape of slopes more comprehensively influences the 3-D architecture of lobes in natural deep-sea systems than previously other lobe deposits and analogue experiments,which helps us better understand the development and evolution of the distal parts of turbidite systems.

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.

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.

An Overview on the Composition and Age of Upper Crust of Proto-Tethyan Lajishan Intra-oceanic Arc,NE Tibet Plateau

Identification and anatomy of oceanic arcs within ancient orogenic belt are significant for better understanding the tectonic framework and closure process of paleo-ocean basin.This article summarizes the geological,geochemical,and geochronological characteristics of upper crust of Proto-Tethyan Lajishan intra-oceanic arc and provides new data to constrain the subduction evolution of the South Qilian Ocean.The intra-oceanic arc volcanic rocks,including intermediate-mafic lava,breccia,tuff,and minor felsic rocks,are distributed along southern part of the Lajishan ophiolite belt.Geochemical and isotopic compositions indicate that the intermediate-mafic lava were originated from depleted mantle contaminated by sediment melts or hydrous fluids,whereas the felsic rocks were likely generated by partial melting of juvenile mafic crust in intra-oceanic arc setting.Zircons from felsic rocks yield consistent and concordant ages ranging from 506 to 523 Ma,suggesting these volcanic rocks represent the relicts of upper crust of the Cambrian intra-oceanic arc.Combined with the Cambrian forearc ophiolite and accretionary complex,we suggest that the Cambrian intra-oceanic arc in the Lajishan ophiolite belt is belonging to the intra-oceanic arc system which was generated by south-directed subduction in the South Qilian Ocean at a relatively short interval between approximately 530 and 480 Ma.

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.

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.

Unveiling three-dimensional sea surface signatures caused by internal solitary waves:insights from the surface water ocean topography mission

Internal solitary waves(ISW),characterized by large amplitude and long propagation distance,are widespread in global oceans.While remote sensing images have played an essential role in studying ISWs,they mainly exploit two-dimensional image information.However,with the launch of the surface water ocean topography(SWOT)satellite on December 16,2022,a unique opportunity has emerged to capture wide-swath three-dimensional ISW-induced sea surface information.In this study,we examine ISWs in the Andaman Sea using data from the Ka-band Radar Interferometer(KaRIN),a crucial sensor onboard SWOT.KaRIN not only provides backscattering satellite images but also employs synthetic aperture interferometry techniques to retrieve wide-swath two-dimensional sea surface height measurements.Our observations in the Andaman Sea revealed the presence of ISWs characterized by dark-bright strips and surface elevation solitons.The surface soliton has an amplitude of 0.32 m,resulting in an estimation of ISW amplitude of approximately 60 m.In contrast to traditional two-dimensional satellite images or nadir-looking altimetry data,the SWOT mission’s capability to capture threedimensional sea surface information represents a significant advancement.This breakthrough holds substantial promise for ISW studies,particularly in the context of ISW amplitude inversion.

An improved algorithm for retrieving thin sea ice thickness in the Arctic Ocean from SMOS and SMAP L-band radiometer data

The aim of this study was to develop an improved thin sea ice thickness(SIT)retrieval algorithm in the Arctic Ocean from the Soil Moisture Ocean Salinity and Soil Moisture Active Passive L-band radiometer data.This SIT retrieval algorithm was trained using the simulated SIT from the cumulative freezing degree days model during the freeze-up period over five carefully selected regions in the Beaufort,Chukchi,East Siberian,Laptev and Kara seas and utilized the microwave polarization ratio(PR)at incidence angle of 40°.The improvements of the proposed retrieval algorithm include the correction for the sea ice concentration impact,reliable reference SIT data over different representative regions of the Arctic Ocean and the utilization of microwave polarization ratio that is independent of ice temperature.The relationship between the SIT and PR was found to be almost stable across the five selected regions.The SIT retrievals were then compared to other two existing algorithms(i.e.,UH_SIT from the University of Hamburg and UB_SIT from the University of Bremen)and validated against independent SIT data obtained from moored upward looking sonars(ULS)and airborne electromagnetic(EM)induction sensors.The results suggest that the proposed algorithm could achieve comparable accuracies to UH_SIT and UB_SIT with root mean square error(RMSE)being about 0.20 m when validating using ULS SIT data and outperformed the UH_SIT and UB_SIT with RMSE being about 0.21 m when validatng using EM SIT data.The proposed algorithm can be used for thin sea ice thickness(<1.0 m)estimation in the Arctic Ocean and requires less auxiliary data in the SIT retrieval procedure which makes its implementation more practical.

Migration and role of zinc in biogeochemical cycles in the Antarctic Ice Sheet and the Southern Ocean

Zinc(Zn),a widespread metal in the Earth’s crust,serves as a crucial nutrient in the Southern Ocean’s primary production.Studies on Zn in Antarctic snow and ice offer insights into the origins of this metal and its transport routes,as well as its impact on the biogeochemical processes within the Antarctic atmosphere–land–ocean system.This review examines research on the spatial and temporal distribution of Zn in Antarctic snow and ice,as well as in Southern Ocean waters.It includes an overview of advanced methods for sampling and analyzing Zn,along with explanations for the observed variations.The review also discusses various sources of Zn as a nutrient to the Southern Ocean.Finally,it addresses prospective issues related to the use of Zn isotopes in identifying atmospheric sources and their biogeochemical effects on the development of the Southern Ocean ecosystem.

Assessment of the Intertropical Convergence Zone over the Atlantic Ocean through an Algorithm Based on Precipitation

The Intertropical Convergence Zone(ITCZ)is a key atmospheric system on a global scale,primarily driven by trade wind convergence near the equator.The ITCZ plays a crucial role in modulating the climate of the borders of tropical continental areas.For instance,Northeastern Brazil experiences a climate influenced by the ITCZ over the Atlantic Ocean.In some periods,the ITCZ exhibits double bands,known as the double ITCZ.While the features of the ITCZ have been described using various approaches and atmospheric variables,there is still a lack of regional studies focusing on the ITCZ and double ITCZ in the Atlantic Ocean.In this context,the main goals of this study are(1)to describe a simple algorithm based on precipitation to identify the ITCZ and double ITCZ,(2)to present a climatology(1997-2022)of the position,width,and intensity of these two convective bands,and(3)to investigate variabilities in the ITCZ characteristics associated with anomalies of sea surface temperature(SST)in the tropical Pacific and Atlantic oceans.The double ITCZ typically occurs southward of the main cloud band,and between February and April,both bands are more distant(~4.5°).In the western sector of the Atlantic Ocean,the ITCZ and its double band extend to more southerly latitudes in austral autumn.Considering the entire Atlantic basin,the annual mean of the latitudinal position,width,and intensity of the ITCZ is 4.9°N,4.2°,and 11 mm/day,respectively,while for the double ITCZ,it is 0.4°N,2.6°,10.3 mm/day,respectively.While the SST anomalies in the Pacific Ocean(El Niño and La Niña episodes)affect more the ITCZ width,the SST anomalies in the Tropical South Atlantic affect both its position and width.