The effect of Oceanic South Atlantic Convergence Zone episodes on regional SST anomalies: the roles of heat fluxes and upper-ocean dynamics

The South Atlantic Convergence Zone (SACZ) is an atmospheric system occurring in austral summer on the South America continent and sometimes extending over the adjacent South Atlantic. It is characterized by a persistent and very large, northwest-southeast-oriented, cloud band. Its presence over the ocean causes sea surface cooling that some past studies indicated as being produced by a decrease of incoming solar heat flux induced by the extensive cloud cover. Here we investigate ocean–atmosphere interaction processes in the Southwestern Atlantic Ocean (SWA) during SACZ oceanic episodes, as well as the resulting modulations occurring in the oceanic mixed layer and their possible feedbacks on the marine atmospheric boundary layer. Our main interests and novel results are on verifying how the oceanic SACZ acts on dynamic and thermodynamic mechanisms and contributes to the sea surface thermal balance in that region. In our oceanic SACZ episodes simulations we confirm an ocean surface cooling. Model results indicate that surface atmospheric circulation and the presence of an extensive cloud cover band over the SWA promote sea surface cooling via a combined effect of dynamic and thermodynamic mechanisms, which are of the same order of magnitude. The sea surface temperature (SST) decreases in regions underneath oceanic SACZ positions, near Southeast Brazilian coast, in the South Brazil Bight (SBB) and offshore. This cooling is the result of a complex combination of factors caused by the decrease of solar shortwave radiation reaching the sea surface and the reduction of horizontal heat advection in the Brazil Current (BC) region. The weakened southward BC and adjacent offshore region heat advection seems to be associated with the surface atmospheric circulation caused by oceanic SACZ episodes, which rotate the surface wind and strengthen cyclonic oceanic mesoscale eddy. Another singular feature found in this study is the presence of an atmospheric cyclonic vortex Southwest of the SACZ (CVSS), both at the surface and aloft at 850 hPa near 24°S and 45°W. The CVSS induces an SST decrease southwestward from the SACZ position by inducing divergent Ekman transport and consequent offshore upwelling. This shows that the dynamical effects of atmospheric surface circulation associated with the oceanic SACZ are not restricted only to the region underneath the cloud band, but that they extend southwestward where the CVSS presence supports the oceanic SACZ convective activity and concomitantly modifies the ocean dynamics. Therefore, the changes produced in the oceanic dynamics by these SACZ events may be important to many areas of scientific and applied climate research. For example, episodes of oceanic SACZ may influence the pathways of pollutants as well as fish larvae dispersion in the region.


Introduction
The South Atlantic Convergence Zone (SACZ) is an atmospheric system occurring in the warm season of South America from November to March (Casarin and Kousky 1986;Kodama 1992Kodama , 1993Kodama et al. 1997;Gan et al. 2009;Grimm 2011;Quadro et al. 2012). It is characterized by marked low-level winds and a moisture convergence region that persists for at least four days (Carvalho et al. 2004;Rosa et al. 2020), and by an extensive northwestsoutheast-oriented cloud band (NW-SE) (Kodama 1992(Kodama , 1993. It extends from the southern center of the Amazon region, crossing the Central-West and Southeast regions of Brazil. Its northern portion can also reach the central-south Bahia state (near 13°S and 42°W), and its southern portion as far south as the states of Paraná (near 25°S and 50°W) and Santa Catarina (near 27°S and 50°W). The SACZ can also extend towards the Southwestern Atlantic Ocean (SWA) as shown in several previous studies (Carvalho et al. 2002;Ferreira et al. 2004;Chaves and Nobre 2004;Brasiliense et al. 2017). Figure 1a, b show composites of oceanic SACZ cases. Those are composites of oceanic SACZ made with outgoing longwave radiation and moisture flux divergence, based on four cases, which are studied and further discussed here.
In the regions affected by the SACZ, one of the main consequences is the occurrence of high rainfall, particularly in late austral spring and summer months (Grimm 2011;Quadro et al. 2012). One of the main mechanisms contributing to the SACZ configurations is the so-called South American Monsoon System (SAMS), which is responsible for moisture transport from northern Brazil and the Amazon region to central and southeast South America (Casarin and Kousky 1986;Kodama 1992Kodama , 1993Kodama et al. 1997;Gan et al. 2009;Grimm 2011;Quadro et al. 2012). Jorgetti et al. (2014), however, suggest that the SACZ can occur in both active and inactive phases of SAMS, resulting in different intensities and positions of the cloudiness band. The SACZ is therefore an important large-scale feature during summertime in tropical South America region (Kodama 1992;Nougés-Peagle and Mo 1997;Satyamurty et al. 1998;Liebmann et al. 1999).
As a result, SAMS and SACZ modulate the seasonal cycle of precipitation over tropical South America in a region between the equator and 25°S (Silva 2009). Gan et al. (2009) showed that 50% of the annual rainfall over tropical and subtropical South America occurs in the austral summer months (Dec/Jan/Feb) and about 90% of it in the period from October to April. Marengo (2005), analyzing the temporal and spatial variability of the moisture balance in the Amazon basin and surrounding areas, showed that the spring and summer periods present a strong convergence of moisture throughout the SACZ climatological positioning. However, atmospheric blocking events occurring in subtropical South America can prevent SACZ formation (Rodrigues and Woollings 2017; Rodrigues et al. 2019), leading to a deficient rainy season in the Central-West and Southeast regions of Brazil.
Previous studies suggest that reduced cloudiness in late spring (end of dry season) in Central-West and Southeast Brazil would eventually increase net surface solar radiation over the Southeast Brazilian coastal region and offshore Southwest Atlantic Ocean. This condition, in turn, would favor a pressure drop, increasing convergence at low levels and an anomalous atmospheric cyclonic circulation in southeastern Brazil. These characteristics, which are associated with increased convection, tend to intensify precipitation in the Central-West region of Brazil and develop atmospheric configurations that end up favoring the establishment of the SACZ (Grimm et al. 2007).
One of the first articles that mentioned the role of the "Oceanic SACZ'' and its implications for precipitation was developed by Carvalho et al. (2002). They argue that ~ 65% of all extreme rainfall events occur when convective activity in the SACZ is spatially extensive and intense. In ~ 30% of cases where intense precipitation occurred north of São Paulo State, they were associated with an intense SACZ with deep convective activity extending towards the Atlantic Ocean, suggesting that the SACZ probably played an important role in the increase of convection in southeastern Brazil. Regarding large-scale forcings, there is an intensification of convection over the Southwest Atlantic Ocean (Carvalho et al. 2002(Carvalho et al. , 2004Ferreira et al. 2004) during El Niño years and greater convection over the continent in La Niña years. Thus, some have proposed that the oceanic SACZ is sensitive to the SST in the Southwest Atlantic and consequently to the precipitation in South America (Tascheto and Wainer 2008). One of the patterns of variability in SST is the South Atlantic Dipole, which is characterized by the surface thermal gradient between the Tropical and South Atlantic. This influence has been verified during neutral ENSO conditions (Bombardi et al. 2014a, b).
The SACZ can also occur in association with other atmospheric and oceanic phenomena, being influenced by local or remote factors (Kodama 1992;Kodama et al. 1997;Grimm and Silva Dias 1995;Grimm et al. 2007;Nogués-Peagle and Mo 1997;Jones and Horel 1990;Marton 2000;Chaves and Nobre 2004;Hirata and Grimm 2015;Pezzi et al. 2016;Pezzi et al. 2021). These issues have been addressed in numerical modeling studies that simulate the SACZ in its atmospheric and oceanic components (Chaves and Satyamurty 2006;Chaves and Nobre 2004), as well as in reanalysis-based studies dedicated to understanding the spatial and temporal variability of the SACZ (Carvalho et al. 2004;Ferreira et al. 2004;Grimm and Zilli 2009). For example, atmospheric frontal systems in the region of the SACZ may interact with cyclonic sub-synoptic scale high-level vortices (Nobre 1988). Oscillations in the period band of 30-60 days may create atmospheric disturbances that trigger the convection associated with the SACZ (Casarin and Kousky 1986) and explosive convection over central and southern Amazonia, responsible for the generation of the convergence zone at low levels (Figueroa and Nobre 1990). Robertson and Mechoso (2000) have presented observational evidence that positive (negative) SST anomalies in the Southwest Atlantic are associated with a weakening (strengthening) of the SACZ. On the other hand, experiments with atmospheric global circulation models show that precipitation in oceanic regions tends to intensify over warm waters (Barreiro et al. 2002).
It should be noted, however, that two possible cooling mechanisms of the ocean surface can be present during SACZ events; one is here referred to as thermodynamic and the other as dynamic. Both will now be discussed and are summarized in Fig. 1c. Chaves and Nobre (2004) and Almeida et al. (2007) have suggested that the occurrence of a Fig. 1 a Synoptic view for all four cases of oceanic SACZ using outgoing longwave radiation (OLR), units in W m −2 . b Moisture flux divergence (colors in × 10 7 s −1 ) and streamlines in 850 hPa. Both panels are made with CFSv2 reanalysis data set. c Schematic representation of the physical processes involved in the modulation of SST in the oceanic layer when the oceanic SACZ is active. The thermodynamic process is at the top and the dynamic process at the bottom of the figure. The up arrows indicate increase, and the down arrows decrease for each individual physical process negative SST anomaly under the SACZ position is a result of negative feedback between the SACZ, cloud cover, and SST, with the atmosphere forcing the ocean. Atmospheric convection over the SWA increases the amount of clouds. These, in turn, reduce solar radiation at the ocean surface, and promote the cooling of it. Convection over cold waters tends to diminish, dissipating the oceanic portion of the SACZ. The tendency of atmospheric general circulation models (AGCM) to underestimate precipitation over cold waters in the SACZ region, as suggested by Nobre et al. (2012), is attributed to the non-active thermodynamic coupling in this class of models, which differs from ocean-atmosphere coupled models. In this way the authors proposed that the oceanic part of the SACZ is a thermally indirect circulation cell, with intense convection occurring over waters with lower SST values. This is the thermodynamic mechanism responsible for the cooling of the surface of the ocean, Fig. 1c top panel.
In addition to the aforementioned thermodynamic process, Kalnay et al. (1986Kalnay et al. ( , 2004 suggested that in one SACZ case the atmosphere forced negative SST anomalies in the Southwestern Atlantic Ocean via low level cyclonic vorticity, which induces upward flow in the underlying ocean (Ekman pumping), bringing colder subsurface waters to the surface layer (upwelling). At the same time, the SACZcloud-SST feedback intensifies the negative SST anomaly. The cyclonic circulation tends to cease over colder waters due to the Marine Atmospheric Boundary Layer (MABL) vertical mixing adjustment to SST. Over colder waters, air buoyancy and turbulence are reduced, increasing the vertical wind shear (Wallace et al. 1989;Pezzi et al. 2005). The same dynamic feedback of SACZ-Ekman pumping-SST is mentioned in Chaves and Nobre (2004), analyzing December, January, and February averaged SST fields. They found, however, that this dynamic mechanism was one order of magnitude weaker than the thermodynamic one. Up to now, from what is known, the thermodynamic mechanism has been taken as the major mechanism responsible for the oceanic surface cooling (Fig. 1b, lower panel). Our study will challenge this idea by reanalyzing in more detail the dynamic mechanism as further discussed in our results and conclusions section.
An important sector of the SE Brazilian coast is the South Brazil Bight (SBB), a highly productive semi-enclosed marine ecosystem extending from Cabo Frio (23° S, Rio de Janeiro State, RJ) to Cabo de Santa Marta (28° 40′ S, Santa Catarina State, SC), where intense fishing activity is concentrated. Its location can be seen in the studies of Castro (2014) and Soares et al. (2011), and in Fig. 1. This is a region encompassing four Brazilian states and a population of about 82 million people. Despite its importance, the SACZ atmospheric effects on the cooling of this Southwest Atlantic region, and associated impact on regional climate and marine ecosystems, has not yet been given due attention. The SBB has a strong seasonal cycle in its physical oceanography-SST, vertical stratification, thermal and haline fronts, upwelling plumes, and currents (Castro 2014;Cerda and Castro 2014), and its biological and fishery productivity (Soares et al. 2011;Dias et al. 2014;D'Agostini et al. 2015;Endo et al. 2019). This oceanic seasonality seems to be highly associated with the atmospheric seasonality (Castro and Miranda 1998). The most important atmospheric system affecting SBB seasonality is the South Atlantic Subtropical High (SASH), the large-scale semi-permanent pressure center that influences the wind patterns on the Brazilian coast and causes the predominantly northeast wind direction in the SBB in summer (Pezzi and Souza 2009b). Some studies indicate a relationship between the sardine biological cycle and fisheries of SBB to the oceanic and atmospheric physical variability (Matsuura 1998;Sunyé and Servain 1998;Cergole et al. 2002;Gigliotti et al. 2010;Soares et al. 2011;Dias et al. 2014, D'Agostini et al. 2015Endo et al. 2019). These results show the important connection between ocean-atmosphere interactions and the SACZ with the ecosystem dynamics of the ocean in the region.
With the above context in mind, our main objective in this study is to investigate the coupling between the atmosphere and the SWA surface layer waters in terms of air-sea interaction processes. We achieve this by using a coupled ocean-atmosphere model to simulate oceanic SACZ episodes (i.e., with strong convective activity over its oceanic portion) selected from the study made by Rosa et al. (2020). Our main interests and novel results are on verifying how the SACZ influences both dynamic and thermodynamic mechanisms in the oceanic mixed layer (OML) that contribute to the sea surface thermal balance in that region, through changes in the OML velocity field and surface net heat fluxes, respectively.
We organize the article as follows: Sect. 2 introduces the regional numerical coupled system used to simulate four intense oceanic SACZ episodes. In Sect. 3 we present a broad overview of the simulated cases and compare three different model setups: (i) atmospheric model with prescribed SST, (ii) oceanic and atmospheric coupled and (iii) oceanic, atmospheric and wave models all actively coupled. At the end of Sect. 3, we present the dynamic mechanism for the whole set of cases, while Sect. 4 presents the mixed layer heat budget analysis for the SWA. The article finishes with conclusions and final remarks in Sect. 5 2 Methods and data

The regional coupled system
In this section we present the Coupled Ocean-Atmosphere-Wave-Sediment Transport (COAWST) Modeling System v 3.2 (Warner et al. 2010). This modeling system combines three well known numerical models: (i) the atmospheric model is the Weather Research and Forecasting Model v 3.7.1 (WRF) (Skamarock et al. 2005); (ii) the Regional Ocean Modeling System svn 797 (ROMS) McWilliams 2003, 2005;Haidvogel et al. 2008) is the hydrodynamic ocean model; (iii) the wave model is the Simulating Waves Nearshore (SWAN) v 41.01AB (Booij et al. 1999); plus a sediment transport model named Community Sediment Transport Modeling Project (CSTM) . All these models can be actively coupled, but for the purposes of this study CSTM is not used. The coupling is performed using the Model Coupling Toolkit v 2.6.0 (MCT, Larson et al. 2005;Jacob et al. 2005). COAWST is an open-source system that is based on individual open-source models (Warner et al. 2010). Pullen et al. (2017) present an interesting review about coupled modeling frameworks where some of the ocean-atmosphere interaction cases were modeled to investigate and forecast the coupling between tropical cyclones and the ocean, such as the hybrid and first registered hurricane in the South Atlantic, the Catarina event (McTaggart-Cowan et al. 2006).
COAWST is a modular system that allows an active coupling between the ocean and the atmosphere, where the fluxes are actively traded in both directions, from the atmosphere to the ocean and vice versa. It has allowed advances in the representation of coastal dynamics due to the coupling of oceanic and atmospheric circulation models and waves and sediment transport models (Lesser et al. 2004;Warner et al. 2008). This system allows the regional simulation of thermodynamic and dynamic large-scale systems for both atmospheric Pullen et al. 2017) and oceanic phenomena (Mendonça et al. 2017). When regionalized, this system can produce simulations in greater detail, allowing the assessment of their effects on regional scales (Mendonça et al. 2017). Each of the system components has been developed and applied to idealized and realistic scenarios. Some of the ocean-atmosphere applications were made to investigate and forecast the interactions between tropical cyclones and the ocean (Bender and Ginis 2000;Bender et al. 2007;Chen et al. 2007;Warner et al. 2010;Mendonça et al. 2017;Miller et al. 2017;Pullen et al. 2017).
The run of each model begins from its own initial and boundary conditions. After a number of user selected time steps, the models reach a synchronization point. From this point on, the flux exchanges between ocean and atmosphere occur via a Model Coupling Toolkit (MCT) scheme that sends and receives data between the coupled models. Since each model may have a different grid, there is a pre-processing in which the interpolation weights are calculated for each model grid through the remapping Spherical Coordinate Interpolation Package-SCRIP (Jones 1998). Sutil and Pezzi (2020) provide a very detailed description of each model that makes up the COAWST and bring a detailed explanation on how to generate the numerical grids of the models as well as the generation of initial conditions, boundary conditions of each model component. In addition, they show how to install and run COAWST on a personal computer.

Models and experimental design
The ROMS is a 3-D oceanic regional model that has a free surface with terrain following vertical sigma coordinates. For the present study, it was configured with a horizontal grid domain ranging from 55° S to 5° S and 70° W to 20° W, to represent the SWA region. The horizontal resolution is 1/12° (approximately ranging from 9.2 to 6 km) and a vertical discretization of 30 sigma levels, with 40 s and 90 s for barotropic and baroclinic time-steps, respectively. The surface and bottom stretching parameters are 5 and 0.6, respectively, with a 50 m critical depth controlling the stretching and default vertical coordinate transformation. It solves the Navier-Stokes equations by using Reynolds averaging, Boussinesq approximation and hydrostatic vertical momentum balance (Shchepetkin and McWilliams 2005). The model has been widely used for many applications such as ocean circulation forecast (Haidvogel et al. 2008), biophysical modeling D'Agostini et al. 2015), biogeochemical modeling (Arruda et al. 2015), air-sea interaction studies (Putrasahan et al. 2013a, b;Seo et al. 2007;Seo 2017;Miller et al. 2017;Pullen et al. 2017) and South Atlantic coastal studies (Mendonça et al. 2017;Combes and Matano 2014;Palma et al. 2008). In the following, the main physical parameterization options used for ROMS are listed. For momentum, horizontal harmonic viscosity is used (Wajsowicz 1993); for tracers equations, a 3rdorder upstream horizontal advection is invoked (Shchepetkin and McWilliams 2005); the pressure gradient algorithm uses a spline Jacobian density (Shchepetkin and McWilliams 2003); the wave roughness formulation uses a Taylor and Yelland relation (Taylor and Yelland 2001); horizontal mixing of momentum is calculated using constant sigma surfaces (Shchepetkin and McWilliams 2003); horizontal mixing of tracers is determined on geopotential (constant depth) surfaces (Shchepetkin and McWilliams 2003) and for vertical mixing of momentum and tracers a generic lengthscale formulation is used (Warner et al. 2005).
The WRF model provides all surface atmospheric fields to the ocean model, such as long and shortwave radiation surface fluxes, precipitation, atmospheric pressure, relative humidity, air surface temperature, and wind velocity at 10 m. It was configured to match the same oceanic domain with horizontal resolution slightly finer than the oceanic, 6 km, approximately. The model also uses terrain-following vertical-sigma coordinates with 45 levels. The physical parameterization schemes chosen for this study were also used in previous studies such as in Shi et al. (2010), Tao et al. (2011) and Nicholls and Decker (2015). The microphysics scheme is Goddard (Lang et al. 2007); the cumulus convection is Kain-Fritsch (Kain 2004); the longwave and shortwave radiation schemes are the new Goddard scheme (Chou and Suarez 1999); the surface layer is the Eta similarity (Monin and Obukhov 1954); the atmospheric boundary layer is Mellor-Yamada-Janjić (Mellor and Yamada 1982;Janjić 2002) and the land surface scheme is Noah (Chen and Dudhia 2001).
The oceanic open boundaries (north, east, and south) were forced using temperature, salinity, current velocities, and sea surface height 5-day means obtained from Simple Ocean Data Assimilation v 3.3.1 (SODA) (Carton and Giese 2008). The National Centers for Environmental Prediction (NCEP) Final Operational Global Analysis (FNL) provided the atmospheric lateral boundaries. The prescribed SST used on WRF uncoupled experiments is from the Real Time Global (RTG) analyses developed at NCEP of the National Oceanic and Atmospheric Administration (NOAA, 2000) for use in weather prediction and modeling, particularly at high spatial resolution and short time range. The version used here has 0.5 degree of horizontal resolution and daily temporal resolution.
The selected cases for the study of the oceanic SACZ and the periods of numerical integration are based on the work of Rosa (2017) and Rosa et al. (2020), in which an automatic classification algorithm was developed and used. It is based on some parameters, such as cloud cover band defined by OLR values lower than 230 W m −2 , and the life cycle and positioning over the SACZ region. From this information the method classifies the SACZ as continental or oceanic and of the oceanic cases we selected four of them based on their well-defined oceanic configurations. The numerical experiments were carried out for two periods. The first goes from December 1st, 2002 to January 31st, 2003. The second was from December 1st to 31st, 2013. The four oceanic SACZ analyzed episodes were extracted from those two integration periods, as shown in Table 1. The atmospheric and oceanic models were initialized at the beginning of those periods from SODA and NCEP FNL, respectively. The investigation started by verifying if the model setups were able to reproduce the observed oceanic SACZ episodes. This was accomplished by integrating WRF in uncoupled mode for both periods described above, and this experiment was named WRF, as shown in Table 1. Experiments named COA, are for the two-way coupled model configuration, using WRF and ROMS; experiments named COA2 uses WRF, ROMS and SWAM, also shown in Table 1. The idea of simulating four cases of oceanic SACZ using three model (multi-model) configurations was to increase the number and variety of experiments that could reduce uncertainties in the oceanic SACZ simulations. That is, a more robust ensemble of results is obtained with simulations performed with different COAWST setups. It should also be noted here that it is not in the scope of this work to analyze and compare the pros and cons of each COAWST setup. Instead, it was to obtain atmospheric and oceanic simulations that would allow a dynamic and thermodynamic analysis, as shown below.

Auxiliary data
The auxiliary satellite data were used mostly to verify if the model was satisfactorily simulating the oceanic SACZ episodes. The objective of this analysis was to visually compare is a set of new high-resolution precipitation estimates based on blending passive microwave (PMW) sensors and infrared (IR) radiometers data to produce estimates of rainfall over the globe at a spatial resolution of 0.1° latitude/longitude and hourly temporal resolution (Ushio et al. 2009). As remarked by many investigators, satellite data supply crucial information where in situ data are very sparse or nonexistent, which is the case over the Southwest Atlantic (the major focus of this work) and over some parts of the continent.

Main features from Oceanic SACZ simulations
In this section, the main results obtained by the set of oceanic SACZ simulations are presented. First, we analyze and show that our set of simulations can reproduce the expected oceanic SACZ effect on SST. Next, their spatial features are compared against satellite data. Finally, an analysis is performed to evaluate the simulated oceanic SACZ dynamics at SWA.

SST time series
As discussed above, during active oceanic SACZ episodes, SST values are expected to decrease in the SWA region affected by the atmospheric convergence. This result was observed in all simulated oceanic SACZ episodes, as displayed in Fig. 2. Figure Table 1 for the experiment names and more detail). The time series are the SST area averaged over the SWA that is under the SACZ's influence, here defined as the region from 20° to 33° S to and from 35° to 45°W. WRF curves are presented only as a reference since they were not simulated but used as prescribed SST bottom boundary conditions for the atmospheric model. As indicated in the figure, SST values indeed tend to decrease during the SACZ episodes as shown in the area averaged four selected cases (Fig. 2a, b). This is clearly observed for the fully active coupled model runs, COA and COA2. For both observed and simulated SSTs, we see that atmospheric conditions associated with the SACZ, particularly the surface winds and cloud cover, act to modify the oceanic surface thermal condition. An opposite situation that has attracted a great deal of attention is seen in case 5 (Table 1 and Fig. 2a) where the SST increases. This is, however, a case where the SACZ is concentrated over the continent (not shown), thus not strongly affecting either the solar radiation reaching the ocean's surface, nor the atmospheric dynamics over the ocean. This is a case of continental SACZ, which is out of the scope of this study.

SACZ spatial features
In order to analyze the spatial characteristics of the simulated set of oceanic SACZ episodes, at each model grid point a time mean of the four cases was calculated based on the dates of the cases shown in Table 1 (Figs. 3,4). To show that distinct model configurations can simulate oceanic SACZ episodes in a reasonable way it is necessary to compare simulated precipitation and outgoing longwave radiation (OLR) with observed ones. This comparison assures us that the model is simulating cases of oceanic SACZ. However, the purpose here is qualitative similarity, not to quantitatively evaluate model bias or deeply investigate the impacts of each specific model setup on the SACZ simulation. All spatial observed and simulated results presented hereafter were calculated as the mean of the four oceanic SACZ simulated cases, based on the dates of the cases shown in Table 1. Figure 3 shows clearly for the three models the northwest/southeast oriented higher values of precipitation band. In some regions, modelled precipitation reaches 55 to 60 mm day −1 , Fig. 3a-c. Note that the WRF simulation shows a larger area with higher precipitation over the oceanic region; this might be associated with the prescribed bottom boundary conditions over the ocean SST values, which are higher for this simulation (Fig. 5b) compared to the other two simulations made with the fully coupled models, COA and COA2 (Fig. 5e, h). A maximum precipitation of that magnitude is not, however, seen in the observed satellite data (Fig. 3d).
The SACZ cloud band position and intensity can be inferred from lower values of OLR (Fig. 4). OLR lower values (a proxy of cloud cover presence) are observed over the oceanic region (as well as over the adjacent continental region). The different model runs were able to simulate reasonably well the SACZ episodes, especially when compared with satellite derived OLR and cloudiness estimated from infrared satellite images (not shown here). Like what was observed for precipitation, the WRF experiment simulates a less intense OLR belt compared to the other two experiments. However, satellite-derived minimum values of OLR ( Fig. 4d) show a larger extension, indicating a larger area of stronger convection as compared to simulation.
Considering simulated precipitation and OLR (generated under different model setups) and observation, it is fair to say that numerical simulations reproduce reasonably well the oceanic SACZ episodes features. Therefore, we can consider the model results a good proxy data set for the study of oceanic SACZ episodes, even considering the bias seen on precipitation and OLR. We also highlight that our main interest is in detecting the mechanisms through which the  Table 1 for the experiment's names and more details. Note that the WRF SST is prescribed and not simulated OML is modified by oceanic SACZ, not necessarily the absolute magnitude of each dynamic or thermodynamic term of the heat balance.

Dynamic aspects and tendency analysis
Here, we will first analyze the sea level pressure (SLP) and associated surface wind field at 10 m height during oceanic SACZ (Fig. 5a, d, g). The western branch of the South Atlantic Subtropical High (SASH) is indicated in these figures by the E and ENE winds. The SASH is the semi-permanent high-pressure system that is predominant over this region, which migrates northwards (southwards) during austral winter (summer) . Note that because of the SASH, the surface winds are upwelling favorable, flowing southwards and almost parallel to the coast in the region from 21° S to 15° S. However, the most remarkable pattern observed in these results is a surface low pressure cyclonic vortex centered at 24° S and 43° W (approximately) on the Southwestern Atlantic Ocean which appears at the surface (Fig. 5a, d, g) and higher levels, such as on 850 hPa (Fig. 5c, f, i). It encompasses all the Southern Brazilian Continental Shelf (SBCS) and offshore oceanic region, being centered in the southern edge of the oceanic SACZ. This striking feature is consistently seen in the three model experiments. A very similar feature is presented by Brasiliense et al. (2017) in  Rosa et al. (2020) for the composite of all oceanic SACZ cases found by the authors. Following Rosa et al. (2020) we designate this feature as the cyclonic vortex Southwest of oceanic SACZ (CVSS). Both effects, the SASH surface high pressure atmospheric circulation and resulting Cabo Frio upwelling, together with the offshore SBCS upwelling caused by CVSS, are discussed in the following section.
The SST time tendency ( SST tend ) was estimated by subtracting from the average SST, calculated for the days of active oceanic SACZ ( SST episodes ) , the average temperature of the day before the beginning of each individual SACZ episode ( SST d−1 ).
SST tendency superimposed on OLR is presented in Fig. 5b, e, h. These OLR fields are slightly different from those in Fig. 4; here we display only OLR values lower than 250 W m −2 , a threshold chosen to isolate the convection. With this threshold, it is possible to better identify areas where convection associated with the SACZ is well organized, active, and possibly affecting the surface waters. The surface cooling produced by SACZ episodes is well characterized in all experiments (blue colors). Note first that the oceanic model with a higher spatial resolution of about 6 km (COA in Fig. 5e and COA 2 in Fig. 5h) produces an SST field with much more details as compared to the prescribed lower resolution SST field used in the WRF simulation. When compared to the observed prescribed SST for the WRF experiment, the COA and COA2 model show a stronger cooling that is not restricted to the region right underneath the SACZ cloud band. Interestingly, this cooling extends further south from the southern edge of the SACZ, reaching the SBCS region and extending very far offshore.
Moisture flux divergence and the associated vector wind field at 850 hPa during the SACZ episodes are displayed in Fig. 5c, f, i. As indicated by the negative values (blue colors), a convergence of moisture flux is present throughout the oceanic SACZ (850 hPa) region. The shape of the moisture flux divergence field over the ocean (simulated by all experiments) follows a typical oceanic SACZ pattern, where it occurs southwest of the SASH. The humidity transport is the main mechanism responsible for the SACZ precipitation over the ocean and further west over the continent. The CVSS presence during oceanic SACZ has been previously identified (Rosa 2017;Brasiliense et al. 2017;Rosa et al. 2020). In cases of continental SACZ this system is seen as a trough, while in oceanic episodes it appears as a closed cyclonic vortex Rosa et al. 2020). Right underneath the cyclonic vortex circulation a large area of sea surface cooling as well as a decreased humidity is observed, as indicated by the positive moisture flux divergence (red colors).

Oceanic SACZ impacts in OML's dynamics and tendency analysis
In this section, we investigate the coastal and offshore upwelling processes during oceanic SACZ episodes and the associated impact on the SST. Several studies have analyzed SST tendencies associated with oceanic vertical velocity and Ekman pumping dynamics, such as Small et al. (2015) for the Benguela upwelling system, and Seo et al. (2012) for the California Current system. The so-called Cabo Frio upwelling system, located at SBB, is a wind forced coastal upwelling, which has also been extensively studied. Rodrigues and Lorenzzetti (2001) showed via numerical modelling the impact of the continental shelf bottom topography and coastline geometry in controlling the position and intensity of upwelling plumes in the region. Castelão and Barth (2006) showed that upwelling favorable negative wind stress curl prevails throughout the year, but is intensified during the summer, when the wind blows for several consecutive days from the NE direction, affecting the SST seasonality.
Here we analyze, via the linear Ekman pumping vertical velocity ( w E ), the importance of the wind stress curl to the upwelling in the SBCS and adjacent offshore region under oceanic SACZ atmospheric dynamic influence.
where ⃗ is the wind stress, ⃗ k is the unit vector in the vertical direction, is the seawater density and f is the Coriolis parameter. The surface wind stress and its spatial and temporal variability is the main atmospheric forcing of the ocean circulation (the wind-driven component) in the near surface oceanic frictional (turbulent) layer but affecting the ocean interior via mass convergences and divergences in that layer (Talley et al. 2011). These are the main dynamical mechanisms connecting the surface turbulent atmospheric layer to the ocean interior. For the southern hemisphere ( f < 0), the negative wind stress curl is the physical mechanism responsible for the positive surface oceanic divergence that in turn induces upward motion in the Ekman layer, bringing deeper and colder waters to the surface layer. On the other hand, positive wind stress curl causes surface convergence, inducing downward water movement.
Coastal upwelling in the SH is normally associated with an along shore wind field blowing predominantly from a direction having the coastline to its right. This induces an offshore surface Ekman flow, a positive divergence on the coast, and a positive vertical flow (upwelling). In this case, the cooling effect on SST will only be effective if cold bottom waters are present near the coast, which is the case in the region during the spring and summer seasons. Note that the coastal upwelling can occur even for a vorticity free wind flow, but a combination of coastal flow divergence with a negative vorticity large scale wind can result in a strengthened regional upwelling. Related to the effect of Ekman pumping on SST variability in the region, we divided the analysis in two parts, the SBCS continental shelf upwelling system, and the offshore deep ocean upwelling. For the sake of brevity, here we show only the results of COA2 simulations, since WRF and COA (mainly) did not significantly diverge on their results from those found with COA2. Additionally, the COA2 experiment shows the lowest bias for precipitation and OLR as compared with satellite observations (Figs. 3, 4).
Large areas of negative wind stress curl of up to − 1.0 × 10 -6 Nm −3 are seen over the coastal SBCS from 17° S to 23° S, including the Cabo Frio coastal upwelling system region (Fig. 6a). Over this coastal region a relatively strong, along shore, southward blowing, and upwelling favorable, wind stress is observed (Fig. 6a). As indicated in the figure, this alongshore wind field has a strong negative vorticity, associated with an offshore intensification of this negative wind stress. The vorticity of the wind produces an Ekman pumping and an average upwelling flow ( w E ), given by Eq. 2, of up to 1.5 m day −1 (Fig. 6b) in that area. The total oceanic vertical velocity at 40 m depth (typical mixed layer depth) as given by the ROMS simulation ( w ROMS ) is shown in Fig. 6c, with horizontal currents superimposed on it. w ROMS has a similar spatial structure to the Ekman pumping field (Fig. 6b) but shows stronger upward vertical velocities of up to 2.5 m day −1 . These higher vertical velocities are likely being produced by a joint effect of wind vorticity, coastal upwelling, and bottom topography-forced vertical flow (mostly in the shelf/slope region). The similarity of w E and w ROMS fields shows that the atmospheric component alone (the curl of the wind stress) as seen in the w E results is highly important, but from the oceanic side the vertical velocities, as modeled by w ROMS , contain additional dynamic factors that regionally modulate and strengthen the w E flow.
On the southern portion of the SBCS, between 23° S and 28° S (Fig. 6), there is a large offshore area of negative wind stress curl, oriented in the northwest/southeast direction. This is well aligned with the oceanic SACZ, but shifted southward (Fig. 6a), as can be seen comparing it with the SACZ convection belt, inferred from OLR in Fig. 5h. This large negative wind stress curl area is associated with a large area of positive vertical Ekman velocities (Fig. 6b). They are not as intense as those seen near the coast, although they still produce considerable upwelling, with vertical velocities ranging from 0.7 to 1.0 m day −1 . This region is under atmospheric influence of the surface low-pressure system, the CVSS, as derived from the set of simulations (Fig. 5a,  d, g), which is producing a cyclonic wind stress circulation and upward water movement seen in w E and in the vertical velocity simulated by the ocean model Fig. 6c.
The presence of the Brazil Current (BC), the western boundary current of the South Atlantic, can be clearly seen in the surface current field generated by COA2 (Fig. 6c). This southward flow is present throughout the year in the region, normally being present offshore of the outer shelf and shelf break regions. It transports warm and saline tropical waters (Legeckis and Gordon 1982;Lima et al. 1996;Assireu et al. 2003;Souza and Robinson 2004;Goes et al. 2019;Pita et al. 2020), and flows to the south up to about 40°S where it meets the colder and fresher northward flowing Malvinas Current (MC). The region where these two currents collide is called the Brazil Malvinas Confluence (BMC). Strong modulations in MABL stability have been documented in the frontal zone separating warm BC from cold MC waters (Pezzi et al. 2005) and by warm core eddy at the Atlantic Southwest (Pezzi et al. 2021). Another interesting feature seen in the simulated oceanic near-surface currents is a train of cyclonic eddies extending in the same direction and region of the negative wind stress curl areas.
Here we highlight three of them with their centers located at 26° S and 42° W, 27.5° S and 33° W, 28° S and 27° W, approximately (marked as A, B, and C in Fig. 6c). This is a deep oceanic region and therefore without the influence of complex topography, but a number of oceanic eddies are frequently seen in the region. Similar eddies in this region have already been reported (e.g. Calil et al. 2021). The largest cyclonic eddy (26° S and 42° W) is right underneath the CVSS and it is possible that its intensity can be affected by this specific atmospheric phenomenon. An upper ocean vertical section at 26° S (not shown) shows a shoaling thermocline near the shelf slope and decreased temperature between 44° W to 42° W, approximately the region of the core of the cyclonic eddy "A" is, as shown in Fig. 6c. Cyclonic eddies also generate upward movement in their cores; values ranging from 0.6 to 2.5 m day −1 are seen associated with these features. Therefore, we can say that under the presence of the oceanic SACZ, Ekman pumping (negative wind curl) and cyclonic eddies may play a joint (near coast), or separate (offshore), significant role in lowering the heat content of surface waters. That is, both Ekman pumping and eddy dynamics might be in effect in the cooling of surface waters of the region.
In order to complement this analysis, we superpose the horizontal net current effects and wind stress on SST (each as tendencies following Eq. 1) as shown in Fig. 6d, e.   Fig. 5 First column, a, d and g shows WRF, COA, COA2, winds at 10 m height superimposed on SLP, in m s −1 and hPa, respectively. Second column, b, e and h shows WRF, COA, COA2, respectively, SST tendency with OLR superimposed, in °C day −1 and W m −2 . Third column, c, f, i shows for WRF, COA, COA2, respectively the moisture flux divergence with winds superimposed at 850 hPa, in 1 × 10 7 s −1 and m.s −1 . Colors indicate SLP in a, c, e, SST tendency in b, d, f and moisture flux divergence in c, f and i. Note that WRF SST in panel b is prescribed and not simulated ◂ Fig. 6 a Wind stress superimposed on wind stress curl (in N m −2 and × 10 −6 N m −2 ). b Wind stress superimposed on Ekman pumping as given by Eq. 2 (in N m −2 and m day −1 ). c Oceanic horizontal surface currents (m s −1 ) superimposed on oceanic vertical velocity (m day −1 ) where positive (negative) value indicates upward (downward) flow. The letter a, b and c denote three cyclonic eddies on the SACZ region discussed in the text. Difference between composites of oceanic SACZ episodes and 3 days average before for each episode, d SST (°C) and surface wind stress (N m −2 ) tendencies and e SST (°C) and surface oceanic currents (m s −1 ) tendencies. All results in this figure are from COA2, averaged for all four cases of oceanic SACZ However, here we used the 3-day averages before the SACZ episodes on the last term on the r.h.s. of Eq. 1. Figure 6e (currents during SACZ minus previous three days) reveals that during SACZ episodes the difference currents are from the south, meaning that all southward flow, including the Brazil Current, is weakened. It is possible to see that within the entire negative area of SST there is a predominance of anomalously weaker circulation from the north, as indicated by the northwards-oriented current vectors (Fig. 6e).
This modulation analysis of the oceanic currents is complemented by the volume transport calculations over the area shown in Fig. S1, supplementary material. We performed this calculation over three distinct latitudes (22-26°S) within our interest area and for three distinct periods. However here we focus our analysis at 26°S right over the oceanic cyclonic eddy seen also in Fig. 6c. The results obtained from the area in Fig. S1, show that the volume transport is − 11.15 Sv (1 Sv = 10 6 m 3 s −1 ) for all cases of oceanic SACZ ( Figure  S1a), meaning that the southward transport is weakened when the oceanic SACZ is established. This is verified by comparing the volume transport of − 18.37 Sv in Figure S1b, where we used the 3-day averages meridional current prior to the SACZ episodes. An intermediate value of volume transport − 16.98 Sv is found for the three months of simulations used in this study (Fig. S1c). Our transport results agree with those shown in Fig. 1 of Schmid and Majumder (2018), which are based in several scientific articles.
This analysis revealed a strong indication that the horizontal advection of temperature is one of the candidates involved in the oceanic dynamics of surface ocean cooling, and this is analyzed in more detail in the heat budget analysis presented in Sect. 4.
It is important to remark that the colder surface waters are not only underneath the oceanic SACZ position, as given by the OLR plots. They are also present southwards and southwestwards of the oceanic SACZ position, showing that the dynamical effect of the atmospheric surface circulation associated with the oceanic SACZ, is not restricted to the region underneath its cloud band. Also, it should be noted that CVSS presence is one of the causes of the intensification of the oceanic portion of the SACZ as shown in Rosa et al. (2020), as well as the oceanic surface cooling southwest of it.

Ocean mixed layer heat balance
In addition to the dynamical effects discussed above, airsea radiative (short and long wave radiation) and turbulent (latent and sensible) heat fluxes at the surface, or their net result ( Q net ) are key elements contributing thermodynamically to the temperature of surface layer waters. In an effort to understand and evaluate the SACZ impacts on the ocean surface waters, the ocean mixed layer heat balance for this area was calculated following similar approaches used in McPhaden and Hayes (1991), Wilkin (2008), Marchesiello et al. (2010), Foltz et al. (2013) and Tamsitt et al. (2016). To our knowledge, this is the first ocean mixed layer heat budget done for this region, for a case of an oceanic SACZ. The vertically integrated mixed layer heat storage equation used is given as (Dijkstra 2008;Moisin and Niiler 1998): where Q net is the surface net heat flux (positive is heat flux into the ocean), T m , u m , v m andw m are, respectively, the vertically averaged temperature, zonal, meridional and vertical velocity components in the mixed layer. c p h T m t is the heat storage (W m −2 ) or temperature tendency T m t ( o C s −1 ), is the water density (kg m −3 ), c p is the water specific heat at constant pressure (J kg −1 o C −1 ). The first two terms on the right-hand side of the heat storage equation are lumped into a single horizontal heat advection term. The horizontal diffusive mixing (horizontal heat diffusion) is the fourth term, and the vertical diffusive mixing (vertical heat diffusion) is the fifth term, and Q net is applied as the surface boundary condition to this term. Specific heat of water at constant pressure ( C p ) and seawater density ( ) are assumed constant and as 3.94 × 10 3 Jkg −1 °C −1 and 10 3 kg m −3 , respectively. The mixed layer depth h was calculated as the depth corresponding to a potential density ( ) change of Δ = 0.125 kg m −3 relative to the most superficial level of the ocean model. This is the same methodology used by Monterey and Levitus (1997) and Cirano et al. (2006) in their global and Southwestern Atlantic Ocean analysis, respectively. A mean ocean mixed layer heat balance for these cases was estimated by time averaging the terms in Eq. 3 over the number of days of active oceanic SACZ episodes.
The terms of Eq. 3 and their components are presented in Figs. 7 and 8, except for the horizontal diffusion term, which was found to be very small (three orders of magnitude smaller than the other terms) and consequently was left out of this analysis. Examining the heat balance results from the individual SACZ cases (figures not shown), it was noticed that cases 1-3 were very similar in their patterns, leading us to conclude that they describe canonical features of oceanic SACZ cases. Based on that, it was decided to average the period of cases 1-3 (as shown in Table 1) considering this as a robust representation of the oceanic SACZ cases (Figs. 7, 8).
(3) The major oceanic heat loss and cooling effect (represented by negative values, upward flux) is due to the evaporative latent heat flux (Fig. 7a), which is followed by long wave radiation (Fig. 7c), and by a smaller contribution from the sensible heat flux (Fig. 7b). In normal conditions (no oceanic ZCAS) these three Q net components are in large part compensated by the incoming short wave solar radiation (Fig. 7d). However, in the area strongly affected by the oceanic SACZ there is a strong attenuation of short-wave  Table 1. a Heat storage. b Vertical mixing (vertical heat diffusion + vertical heat advection). c Horizon-tal heat advection. All terms are integrated from surface down to the mixed layer depth (-h) and multiplied by C p . The panels in this figure were obtained from COA2 experiment. All units are in W m −2 heat flux at surface of the order of 100-130 W m −2 due to the high cloud cover, as clearly seen in Fig. 7d.
The sum of the individual surface heat flux terms (Q net ), is presented in Fig. 7e. Under the oceanic SACZ's offshore position, the surface net heat flux is either positive but near zero, or it is negative (white/light green shades in Fig. 7e). This seems to be a combination of a reduction of incident shortwave radiation caused by the cloud cover associated with oceanic SACZ and mostly an increase of outgoing latent and sensible heat fluxes caused by a joint effect of lower humidity (negative moisture flux divergence; see Fig. 5i) and stronger offshore winds. Additionally, the slightly negative values around 26°S and 42°W in Fig. 7e match with the cyclonic atmospheric vortex position (see Fig. 5d) and are also in the southern edge of the cloud cover (see Fig. 5e), indicating that the heat balance over this region is a result of a complex interplay among Eq. 3 terms, not simply a consequence of the oceanic SACZ cloud cover blocking the shortwave to reach the ocean surface (thermodynamic mechanism). Some regions of higher negative values of surface sensible and latent heat fluxes can be observed south and close to the cloud band of the SACZ. It is likely that they are caused by an enhanced air-sea temperature difference induced by sea surface cooling produced by Ekman pumping and mesoscale eddies, also by regions of lower surface humidity modulated by moisture flux divergences and the presence of stronger surface winds (see e.g., wind plot in Figures and 5i) factors which enhance both, the sensible and the latent heat fluxes. To some extent, this is like what is shown in Fig. 10 of Tirabassi et al. (2015) where near the coast there is a marked loss of heat through the air-sea fluxes. In addition, in their study it also extends farther south of the SACZ position. Our results also show that the net surface heat fluxes tend to cool the ocean underneath the cloud band (see Fig. 5b, e, h) but the heat storage loss (Fig. 8a) southwest of SACZ is due to a combination of net surface heat fluxes from other terms in the heat balance, related to the dynamic mechanism that will be further analyzed.
The sea surface cooling in Figs. 5e, 6d, e underneath the SACZ is unquestionably due to the reduction of shortwave incidence at sea surface (shown in Fig. 7d), as a consequence of increased cloud cover. This cooling mechanism is the same as suggested by Chaves and Nobre (2004) and Nobre et al. (2012) as the negative thermodynamic mechanism, as illustrated in Fig. 1. This can be seen in Fig. 7, where the negative values of net surface heat fluxes ( Q net ) are evident. This is the case when thermodynamics feedback is playing a role on the sea surface cooling. Chaves and Nobre (2004) have argued that the thermodynamic negative feedback is the main cooling process of the sea surface due to the SACZ's cloud band presence. In their study they argue that this thermodynamic effect is one order of magnitude (~ 0.4 °C/ month) larger than the oceanic dynamic effect (~ 0.02 °C/ month), caused by Ekman pumping. However, Chaves and Nobre (2004) calculated the effect of the dynamic heating/cooling of the SST considering only the vertical heat advection term caused by Ekman Pumping, whereas here we consider advection (horizontal and vertical) and a changing mixed layer depth. In addition, their work was made over 120 consecutive days, which includes days of active SACZ and days without SACZ, while here, through case selection, only active SACZ days are used. This may explain the difference from our results, where the heat budget clearly shows that vertical mixing and horizontal advection in Fig. 8b, c are as important as the air-sea heat fluxes in controlling the SST variability, as well as the negative feedbacks associated with it.
Those two terms (Fig. 8b, c) have the same order of magnitude as the surface heat fluxes (Fig. 7e). The stronger green shaded negative values shown in Fig. 8b, c are approximately located where Ekman pumping is active (yellow/red values in Fig. 6b), showing the importance of this dynamic process (i.e., vertical mixing) in controlling the surface water temperatures. We can see that both vertical mixing and horizontal advection are contributing to intensify the surface temperature decreases seen in Fig. 8a via the oceanic mixed layer heat budget. This behavior occurs for all the three sets of simulations (not shown).
Negative heat advection effects seen in Fig. 8c are present both in offshore open ocean regions and along the coast, except for some parts of São Paulo and Rio de Janeiro states. Considering its importance, it is worthwhile to discuss the possible ways in which the horizontal advection of temperature can be changed. In contrast to the other terms of heat budget equation, which are mostly local (they can be estimated at each grid point regardless of the neighboring grid points), the horizontal advection depends not only on the local vector current, but also on the temperature gradient field, which is determined by the neighboring grid points. Even in a situation in which the flow field is stationary, changes in the temperature field (caused inhomogeneous surface heat fluxes or dynamical ocean processes) will change the temperature gradient field and consequently the scalar product of the flow field with it and the horizontal advection. Another important way of modulating the horizontal advection is by the changes of horizontal flow field induced by the presence of mesoscale eddies. Regions where the mesoscale flow is up the temperature gradient are regions of negative horizontal advection and will contribute to cooling. Therefore, the cooling seen underneath the SACZ position (and south of it) has a contribution from the horizontal advection (Fig. 8c), which is produced by a combination of the mesoscale eddy field flow changes, and temperature modulation produced by eddy induced upwelling at their cores and vertical flow produced by Ekman pumping dynamics associated with the curl of surface wind.
A decrease of surface water velocity of the Brazil Current and its southward heat transport, as discussed in Sect. 3, will impact the heat balance, but the cause of this velocity decrease and its precise magnitude is still a matter for further studies. Near the shelf break and open oceanic regions negative advection terms of about − 200 W m −2 are seen in a southwestward shifted region from the SACZ's position (Fig. 8c), where Q net values are positive (Fig. 7e) and coincident with the absence of the cloud cover band region (Figs. 3, 4), around 42° W to 39° W and 30° S to 26° S. However, it must be noted that at this region, vertical mixing caused by CVSS is also playing a role on sea surface cooling, as seen in Fig. 8b.
Our analysis of the ocean's thermal response based on the oceanic mixed layer heat balance shows that ocean-atmosphere interaction under the oceanic SACZ is very complex and cannot be simplified to being viewed as a merely a consequence of the reduced incident shortwave due to the cloud cover presence. Instead, ocean cooling processes are produced by both dynamics and thermodynamics. Our findings suggest that, contrary to what has been previously proposed, the thermodynamic processes through surface ocean-atmosphere fluxes, while important underneath the SACZ, are of the same magnitude as the ocean dynamic processes in controlling the SST of that region. In this region the horizontal advection plays a significant role. Moreover, to the south of it, both horizontal advection and Ekman pumping cause cooling, especially where the CVSS is located, as shown on the sequence of our results.

Final remarks and conclusions
Our results show that oceanic SACZ episodes strongly modulate the SST by a complex interplay of the combined effects from significant changes in both oceanic dynamics and thermodynamics in the mixed layer. The schematic diagram presented in Fig. 9 summarizes the main physical ideas of ocean-atmosphere interactions present during ocean SACZ episodes, as well as the findings of our investigation that are discussed below.
The low-level atmospheric circulation and associated airsea flux modification caused by the convergence of humidity flux, together with changes in the Ekman pumping and eddy dynamics, and the presence of an extensive cloud cover band over Southwest Atlantic Ocean result in sea surface water cooling through dynamic and thermodynamic mechanisms (Fig. 1). As shown in the ocean mixed layer heat balance and Ekman pumping calculations, horizontal advection and vertical mixing (as part of dynamic mechanism) have the same order of magnitude as the net surface heat fluxes (as part of the thermodynamic mechanism). The surface atmospheric circulation patterns established during oceanic SACZ cases seems to promote a reduction of the heat transport to the south associated with the Brazil Current and neighboring offshore waters. An enhancement of the mesoscale eddy field seems to be produced with the net effect of reversing the direction of the horizontal heat transport in various regions affected by the oceanic SACZ. This can be seen from the heat budget's results (Sect. 5) combined with the joint analysis of the current and the SST in Fig. 6c, and mainly from their tendencies in Fig. 6e (Sect. 4).
On the region right underneath the oceanic SACZ, contrary to what was previously thought and reported in the literature, upwelling is not a major player on causing surface cooling, even though it is caused by well-established dynamic atmospheric circulation (Fig. 6b) and has shown its presence in dynamic oceanic analysis (Fig. 8b). We suggest that the surface cooling by Ekman upwelling underneath the cloud band does not occur because the water just below the mixed layer bottom in that region is not colder than the first layers of the ocean as our simulations indicate (figures not shown).
A striking feature found in this study is the presence of a cyclonic atmospheric vortex (previously defined as CVSS) both at the surface and aloft at 850 hPa, south of the SACZ. The surface atmospheric circulation established by CVSS produces a surface water cooling southwest of the SACZ position that seems to be produced by a divergent surface layer Ekman transport. This cyclonic atmospheric system is rarely discussed in the literature: right underneath it, the model results show a large area of cooled surface waters, as well as decreased air humidity, indicated by the positive humidity flux divergence. Rosa et al. (2020) found that during continental SACZ episodes there is a trough near the Uruguayan coast but during oceanic SACZ episodes this trough closes into a cyclonic vortex positioned southwest of the SACZ, accompanied by an intensification and northward displacement of the SACZ cloud band. According to our present results, the proposed mechanism of shortwave-SST negative thermodynamic feedback ( Fig. 1b and represented on schematic in Fig. 9) apparently is not acting alone on oceanic surface cooling, as previously argued in the literature (e.g. Chaves and Nobre (2004); Chaves and Satyamurty 2006;Almeida et al. 2007;Nobre et al. 2012, among others). The analysis of intense oceanic SACZ episodes here presented shows that oceanic dynamics, as seen in our ocean mixed layer heat balance analysis through vertical mixing and horizontal advection, significantly contributes to oceanic surface cooling underneath of SACZ position. In addition, this dynamical mechanism is dominating the cooling to the southwest of the SACZ, (Fig. 5h) where the CVSS (Fig. 5g) is producing a cyclonic wind stress circulation (Fig. 6b) that induces horizontal advection and upward water movement seen in w E and in the vertical velocity simulated by the ocean model (Fig. 6c). The oceanic analysis of this phenomenon was complemented by the calculation of the thermal balance terms in the ocean mixed layer, which elucidated the role of horizontal advection (Figs. 6c-e, 7c) as an important player in the surface cooling of the ocean. This cooling caused by horizontal advection in that region is a noteworthy new result of this study. It is our belief that the results obtained here do not totally reveal the full interplay of ocean-atmosphere interactions in oceanic SACZ episodes modulating the oceanic dynamics and thermodynamics. Since BC shows a high intrinsic variability in the intraseasonal scale, a reasonable point could be made on how significant are the transport anomalies which were observed in our model results relative to the mean variability. It is difficult to come to any conclusions on this regard from our free run experiments spanning only 3 months in total. However, as explained before, we consider our results more as a kind of sensitivity and processes-oriented experiments in which we try to highlight the different atmospheric and oceanic processes operating when oceanic SACZ is present or not.
Nevertheless, our results expand upon previous knowledge about the SACZ, and bring out new points for further investigations, both in modeling and observations, and stressing the need for considering the SST changes due to modifications in the ocean dynamics. This subject Fig. 9 Schematic representation of Oceanic SACZ processes, showing the ocean dynamic and thermodynamic mechanisms that decrease the SST, which is represented by the darker blue area. The convergence of surface winds and divergence aloft (curved red arrows around dashed green line) causes the weakening of the Brazil current (BC) (represented by the thicker darker red arrow parallel to the continent) and consequently decreasing the heat transport towards to the south, likely caused by this weakened current. This generates cold water anomalies, as indicated by dark blue thicker northward blue arrows. This CB's representation should be interpreted as the whole current action which extends far offshore but is not represented in its totality here for the sake of schematic's clarity. Upwelling is caused by well-established dynamic atmospheric circulation and oceanic eddy activity (darker blue spiral). The thermodynamic mechanism which is a reduction of the shortwave radiation reaching the sea surface is represented by the thick cloud cover presence (gray clouds), which is a remarkable SACZ feature. The presence of a Cyclonic Vortex Southwest of SACZ (CVSS) is represented by the thicker curved red arrow