Downstream development during ridging South Atlantic Ocean anticyclones

Ridging South Atlantic high pressure systems (ridging highs) are often accompanied by cut-off low (COL) pressure systems aloft, but may also occur without them, in which case a linear baroclinic wave would be observed propagating across the South African domain in the upper troposphere. Using 41 years of ERA-5 reanalysis data, this study documents differences between the characteristics of the prevailing dynamical processes and associated local eddy kinetic energy generation, its downstream transfer and dissipation during these two scenarios. The study shows that when COLs are present then baroclinic conversion is strong and it is confined east of the Greenwich Meridian, whereas it is located downstream of South Africa and it is much weaker, when ridging occurs without COLs. The differences in strength and locations of the baroclinic conversion are associated with the differing jet streak configurations between the two scenarios; which lead to Rossby wave breaking and the absence thereof when there are COLs and when ridging occurs without COLs, respectively. The presence of breaking during COLs leads to trans-ridge downstream development that facilitates energy transfer from the midlatitudes into the South African domain. When there are no COLs present, the trans-trough downstream development is stronger than it is across the upstream ridge. Barotropic conversion from eddy kinetic energy to mean kinetic energy occurs in the South African domain during COLs, but occurs much further downstream when there are no COLs during ridging highs. The difference in the characteristics identified in this study can be traced back to the differences in the potential vorticity anomaly structures, which are largely due to whether the waves break or do not during the evolution of ridging events in the South African domain.


Introduction
The subtropical South Atlantic Ocean is characterised by a quasi-stationary anticyclonic circulation that exhibits pronounced variability in intensity and position at intra-annual and inter-annual (Sun et al. 2017) and multidecadal (Reason 2000) time scales. This circulation may potentially be impacted by future climate change (Engelbrecht et al. 2009;Reboita et al. 2019) and therefore expected to change at decadal and centennial time scales. Under certain stratospheric conditions (Ndarana et al. 2018), the anticyclone extends east and ridges around the South African landmass to induce an onshore flow into the country (Ndarana et al. 2021a). This eastward extension process is referred to as ridging. Ridging South Atlantic Ocean high pressure systems (hereafter referred to as ridging highs) draw moisture from the South Atlantic and Southwest Indian Oceans (Rapolaki et al. 2020) by means of geostrophic and ageostrophic moisture fluxes (Ndarana et al. 2021a) and, by so doing, contribute to the moisture budget of the country. These oceans are, however, not the only sources of moisture for South Africa, because water vapour may originate from the Equatorial South Atlantic Ocean (Vigaud et al. 2007) and Equatorial Indian Ocean (D'Abreton and Tyson 1995) and transported into tropical Africa. The Angola low (Reason and Jagadheesha 2005;Crétat et al. 2019;Howard and Washington 2018) then organises it in a south easterly direction (Hart et al. 2010) and tropical/extra-tropical interaction events then transport the moisture further south (Ratna et al. 2013;Macron et al. 2014) into South Africa.
Ridging highs are consequently important in terms of the dynamics of rain-producing weather systems of South Africa because they are low level processes that co-occur with moving troughs (Piva et al. 2008;Ndarana et al. 2022) that may develop into cut-off low (COL) pressure systems (Taljaard 1985;Singleton and Reason 2007) when the waves break (Ndarana and Waugh 2010;Reyers and Shao 2019;Barnes et al. 2021). Heavy rainfall that sometimes leads to floods which cause extensive damage to property and loss of life often results from ridging highs that occur in tandem with COLs (Trigaardt et al. 1988;Singleton and Reason 2006;Pinto et al. 2022). Moreover, ridging highs may also be associated with tropical-temperate troughs (Hart et al. 2010) whose characteristic outgoing longwave radiation anomaly pattern (Fauchereau et al. 2009) being negative and reaching maximum values as the ridging process matures (Ndarana et al. 2021a). It follows then that understanding the dynamics associated with this process at both the lower and upper levels may play an important role in informing predictability studies of extreme rainfall events in the South Atlantic Ocean/South African sector.
One approach that previous studies have employed in understanding upper level dynamical processes is the local energetics framework (Orlanski and Katzfey 1991) to diagnose the evolution of baroclinic waves (Orlanski and Sheldon 1995), which gives rise to the notion of downstream development Sheldon 1993, 1995). The local energetics framework has since been applied to a broad range of atmospheric dynamical settings, beyond the evolution of baroclinic waves, including intensifying jet streaks (Lackmann et al. 1999), surface cyclolysis (McLay and Martin 2002), life cycles of conservative extratropical cyclones that exhibit explosive deepening behaviour (Decker and Martin 2005), extratropical transitions of tropical cyclones (Harr and Dea 2009), travelling mid-tropospheric troughs (Piva et al. 2010), the evolution of cut-off low (COL) pressure systems Piva 2013, 2016;Ndarana et al. 2021b;Pinheiro et al. 2022), dynamical processes that lead to the formation of frost events (Müller et al. 2017), the processes that maintain jet wave trains (Li 2021) and dust storms (Rivière et al. 2015); to name a few. Some studies (e.g. Lackmann et al. 1999;Danielson et al. 2006) have also linked the energetics framework with wave activity fluxes. It might therefore find some applications to upper level dynamical processes that are associated with ridging pressure high systems as well, as these processes are baroclinic waves. Ndarana et al. (2022) showed that there are at least two types of ridging high evolution that occur in the South African domain. These were categorised according to which latitude they occur, relative to the 40 • S latitude line. Type-S (Type-N) ridging events occur south (north) of this latitude line. They also showed that the different types of ridging highs are associated with baroclinic waves aloft, whose troughs pass over South Africa as the ridging process is initiated. Whilst other studies (e.g. O'Brien and Reeder 2017) have analysed the evolution of these waves in the South Atlantic Ocean/South African sector using their interaction with the jet as a basis for their analysis; little is known about the dynamical characteristics of these waves in the South African domain from a downstream development perspective. Recent studies (e.g. Ndarana et al. 2021b;Pinheiro et al. 2022) showed that COLs are characterised by such downstream development, which is facilitated by ageostrophic geopotential fluxes that draw eddy kinetic energy from the mid-latitudes into the COL's formation region. This in turn induces the closed cyclonic circulation of the COL, together with the effect of the high potential vorticity (PV) anomaly that does the same (Hoskins et al. 1985). This high PV anomaly is normally caused by Rossby wave breaking (RWB), defined as the rapid and irreversible overturning of PV contours which was first identified during boreal winter stratosphere (McIntrye and Palmer 1983). RWB has been shown in the upper troposphere to be associated with COLs (Ndarana and Waugh 2010;Reyers and Shao 2019;Portmann et al. 2021;Barnes et al. 2021) and rainfall (de Vries 2021).
It is well known from forecasting experience in southern Africa that COLs are often accompanied by ridging highs at the surface, but the opposite is not necessarily true. This is expected to be the case as Ivanciu et al. (2022) showed that 44% of ridging highs are associated with RWB events. It follows then that as not all RWB cases are associated with COLs (Ndarana and Waugh 2010;Barnes et al. 2021), it may then be expected that even fewer ridging highs will have COL association. When COLs do occur in the domain, the associated RWB events are observed somewhere in the region near the Greenwich Meridian to about 25 • E and they often have an anticyclonic orientation either with an equatorward (Thorncroft et al. 1993) or a poleward (Peters and Waugh 1996) bias. In this study the proportion of ridging highs that are associated with COLs will be quantified.
Given this, there are many instances when ridging highs occur with ordinary trough/ridge systems aloft that do not exhibit the wave breaking that is observed during COLs. Piva et al. (2010) analysed the energetics of such troughs in the South American domain and found that the presence of the Andes mountains amplifies the troughs as they propagate across the land and affects the energetics. In the South African domain, where the mountains are much lower than the Andes, only the energetics of COLs have been considered (Ndarana et al. 2021b) and therefore the effect of the absence of breaking baroclinic waves that lead to isolated potential vorticity anomalies (or cutoffs) has not been considered yet. As such the questions raised in the study are • How does the absence of COLs (and associated RWB) affect the upper level dynamical processes as the linear baroclinic waves propagate through the South Atlantic Ocean/South African sector during the evolution of ridging highs? • Do the prevailing upper level ridging highs dynamical processes lead to different eddy kinetic energy generation, transfer and dissipation?
The questions raised above will complement recent energetics studies in the South African domain (Ndarana et al. 2021b) and the Southern Hemisphere in general (Pinheiro et al. 2022) by considering the case of linear baroclinic wave energetics that might lead to rainfall over South Africa due to their association with the ridging process. The rest of the paper is structured as follows: in the next section a brief review of the energetics diagnostic framework is presented, followed by Data and Methods in Sect. 3. The results are presented in Sect. 4 and the Concluding remarks are given in Sect. 5.

Review of flow processes and downstream development in baroclinic waves
To characterise downstream development, the variables of interest are first divided into their basic state and perturbation parts, where the former are represented by capital letter symbols with a subscript m and the perturbation fields are the small letters so that the instantaneous horizontal velocity and geopotential fields are separated into their basic state and perturbation portions as V = V m + v and Φ = Φ m + , respectively. The horizontal flow may be decomposed into its zonal and meridional components as V = Ui + Vj and the three dimensional flow is represented by U = V + k , where is the vertical velocity in pressure coordinates, so that u becomes its perturbation. Note that i, j and k are unit vectors in the zonal, meridional and vertical directions, respectively. On this basis, the eddy kinetic energy is defined as K e = 1 2 v ⋅ v and the flux form (Orlanski and Sheldon 1993;McLay and Martin 2002) of the equation that describes its evolution is then given by The first and second terms on the right hand side of Eq.
(1) are the horizontal and vertical eddy kinetic energy flux convergence, respectively. The third term can be interpreted as the advection of the perturbation geopotential field by the horizontal perturbation velocity field. The fourth term is the Reynolds stress term, which is the conversion of K e to mean kinetic energy and the "Residual" is obtained by subtracting the terms just described from t K e . In all these terms the subscript p indicate that the horizontal gradient and divergence operators are evaluated on isobaric surfaces. Note that ∇ p (⋅) and ∇(⋅) are the two and three dimensional gradient operators, respectively, so that ∇(⋅) = ∇ p (⋅) + k p (⋅) . The pressure work term, v ⋅ ∇ p , plays a dominant role in this eddy kinetic energy budget equation and various processes are embedded in it (Orlanski and Katzfey 1991). These processes are most clearly seen when it is decomposed in the following manner, Here = − p and the rest of the symbols have already been defined. Once decomposed, it becomes clear that −v ⋅ ∇ p is a dominant term in Eq. (1) because it embodies all the K e generation processes. The first term on the righthand side, − , is the baroclinic conversion, which when it is positive (i.e. − > 0 ), represents the conversion from eddy available potential to eddy kinetic energy. There clearly needs to be vertical motion so that baroclinic conversion is essentially facilitated by vertical heat fluxes (Li 2021). The second term is the convergence of ageostrophic geopotential flux of K e which is a horizontal process and the fluxes themselves, (v ) a , as represented by the red curved arrows in Fig. 1, serve to radiate energy in the direction of the ageostrophic flow, across a ridge axis, so that a downstream, K e , may then develop at the inflection point immediately east of the ridge axis. Across a trough axis, the downstream directed ageostrophic flux is caused by the fact that the flow is subgeostrophic (with a ageostrophic flow directed upstream) and < 0 . The term p ( ) is the vertical geopotential flux convergence, and it is considered small when the fields are vertically integrated from 1000 hPa to 100 hPa.
Based on Orlanski and Sheldon (1995) for baroclinic waves in general and for the COL case in the South Atlantic Ocean/southern African sector (Ndarana et al. 2021b) and Southern Hemisphere in general (Pinheiro et al. 2022), downstream development is a three stage process. During the first stage a K e centre I is generated by means of baroclinic conversion ( − > 0 ), as represented by the red hatched oval shape in Fig. 1. At the bottom of the wave, across the ridge axis, the supergeostrophic flow, as represented by the downstream directed ageostrophic flow, combined with > 0 , combine to produce the ageostrophic geopotential fluxes that begin to transport K e downstream.
As K e I decreases in strength from stage 1 to stage 2, K e II increases in strength to reach a maximum during the latter. This may largely be due to the increase in positive perturbation geopotential anomaly in the ridge and the perturbation ageostrophic flow, thus increasing transport of energy downstream via the ageostrophic geopotential fluxes, which are particularly strong in COLs. This strength is caused by anticyclonically breaking waves which induce deep geopotential anomalies and a strong supergeostrophic flow ahead of the closed circulation (Ndarana et al. 2021a;Reyers and Shao 2019;Ndarana et al. 2021b). The negative perturbation field across the trough axis (downstream of energy centre II) begins to appear and the upstream ageostrophic flow lead to energy transfer toward the inflection point beyond the trough axis to cause K e III to develop. As was the case with the strengthening ridge, as the trough deepens, the flux increases; thus increasing K e centre III.

Data
We use the Fifth Generation European Centre for Medium-Range Weather Forecasts Reanalysis (ERA5, Hersbach et al. 2020) products from 1979 to 2020. The horizontal grid spacing of each dataset is 2.5 • at 6-hourly time intervals. Even though these products are available at finer mesh grids than this, we deem the chosen grid spacing sufficient because ridging high pressure systems are synoptic scale processes, the horizontal scale of which is ∼ 10 6 m (Holton and Hakim 2014). The variables used are mean sea level pressure (MSLP), the zonal and meridional wind components and geopotential heights to calculate the geostrophic(v g i + v g j ) and ageostrophic ( v a i + v a j ) wind fields at standard pressure levels. In this study we assumed a variable f geostrophic flow (Blackburn 1985) so that it is not non-divergent (Cook 1999). Unlike other reanalysis products, ERA5 comes with potential vorticity (PV), which is a derived field from the flow and potential temperature. In order to linearly interpolate PV from isobaric to isentropic surfaces, the later is calculated directly from the temperature output. Shematic representation of a propagating baroclinic wave (black curves) together with its evolving eddy kinetic energy ( K e ) density centres (yellow and brown shaded oval shapes and labelled in Roman numerals I, II and III-the brown shade represents increasing strength of eddy kinetic energy) and geopotential perturbation ( ) patterns (oval shapes with no shading. The thin black solid (dashed) contours represent area in the wave where > 0 ( < 0 ). The black curved arrows represent the ageostrophic flow ( v a ) and the red curved arrows represent the ageostrophic geopotential fluxes ( v a ). The solid (dashed) straight line marked as A (B) represent the ridge (trough) axis. The increasing size of the arrows ( v a , v a ) and increasing number of contours ( K e , ), represent increasing intensity of these fields. The red (blue) unshaded oval shapes represent areas in the baroclinic wave where the perturbation meridional velocity field is negative (positive). The hacthed region in eddy labelled I represent baroclinic conversion from eddy available potential energy to eddy kinetic energy ( − > 0 ). [The schematic is adapted from Orlanski and Sheldon (1995) but also draws from the findings of Ndarana et al. (2021b) and Pinheiro et al. (2022) for COLs in the South Atlantic Ocean/southern African sector and Southern Hemisphere, respectively.]

Methods
We objectively identify ridging highs using a simple algorithm consisting of three steps. Its details are provided in Ndarana et al. (2018) and only a brief description of it is provided here. In the first step, closed contours in the 6-hourly MSLP fields are identified within a domain bounded by 40 • W and 60 • E. We then group these closed contours so that concentric contours in the South Atlantic Ocean represent the South Atlantic Ocean High (SAOH). This is the second step. In the third and final step, if the outermost contour extends east of the 25 • E longitude line, we record such instances as the ridging process having occurred. If this condition is met on consecutive time steps, without any breaks in between, then this constitutes a single ridging event and is the basis on which the duration of the events is determined.
We also identify COLs at the 250 hPa isobaric level using the method of Barnes et al. (2021). In the first step closed contours at 1 gpm geopoential intervals and their centres are identified and concentric contours are grouped together, in a similar fashion as the ridging high identification algorithm discussed above. This step produces closed circulations, that may be cyclonic or anticyclonic, whose centres are considered to be the mean of the longitudes and latitudes of the concentric contours. To discard the former, each of these centres are then required to be 10 gpm lower than at least 6 of the surrounding 8 grid points on a 2 × 2 grid box, so that the retained closed circulations are low pressure systems. To impose the closed core condition, we employ the 250-500 hPa thickness, which was found in Ndarana and Waugh (2010) to perform better than temperature fields. If the thickness at the centre of the circulation is lower than at least 5 of the surrounding 8 grid points, then the cold core condition is considered to have been met. Following Barnes et al. (2021), gaps are filled. After tracking the COL events, if two of them occur within 12 h and 500 km of one another, these are considered to be the same event. As in Singleton and Reason (2007), COLs that affect South Africa are considered to be those that occur within the grid box bounded by 10 • E and 40 • E, and 20 • S and 40 • S for at least 24 h.
To establish the climatological behaviour of the ridging events and associated fields, we calculate composite means. To understand the energetics associated with these events more clearly, we isolate cases when these occur without COLs and compare them to cases when the former occurs in the South African domain. To do this, for each ridging event we search for a COL that occurs during the evolution of the former. If no COL is found, then we group the ridging events as "ridging without COLs". Composites for this group are then generated using the duration of the ridging events as the only basis for the calculations. Corresponding times during the evolution of the ridging events are brought together and averaged. For instance, all the fields occurring at the time step that ridging occurs (e.g. t = 0 h) occurs are averaged, then those at the second time step are brought together and averaged (t = +6 h) and so on. The COL events are further sub-divided into groups that occur in four quadrands of the South African domain, following Singleton and Reason (2007) and refer to these as south-western, northwestern, north-eastern and south-eastern sub-regions (see cf. Fig. 1 in Singleton and Reason 2007). Following the categorisation, the same procedure of calculating composites is followed, but using the day of COL formation as t = 0 hrs instead of that of the ridging highs. The reason for this is to bring out the characteristics of the COLs more clearly so that a clearer comparison between the two categories of events may then be possible. The statistical significance of the composite anomaly fields is established following Brown and Hall (1999).
The diagnostics described in Sect. 2 are first calculated for each COL and ridging event and then composited in the manner just described. As noted in that section, the variables are decomposed into their basic state and perturbation components. In this study, the basic state variables are defined as the 31-day mean for each event, centred on the first time that COLs and ridging events without COLs are identified. Note that all the energetics diagnostics are integrated vertically from 1000 to 100 hPa, as previous studies have done, unless otherwise stated.

General baroclinic wave characteristics
A total of 2423 ridging events were objectively identified during the study period and 21% (497 events) of these occurred when there are COLs aloft. This proportion is less than that observed for RWB (Ivanciu et al. 2022) for reasons already mentioned in Sect. 1. As noted in Sect. 3.2, the COLs that occur with ridging are subdivided into four categories, according to where they occur in the South African domain and 212, 107, 102 and 76 of them were located in the south-western, north-western, north-eastern and south-eastern quadrants of the domain, respectively. In all the composite figures shown, for brevity only the south-western COLs events are presented in the left panels. This quadrant is also of particular interest since it contains twice as many events as any of the other quadrants and because COLs here play an important role in helping to mitigate severe winter drought (De Kock et al. 2021). The results shown below are broadly similar between events in each quadrant, with only a few small differences in the composite structure and evolution of the MSLP fields. Note that in all the figures, the left hand panels composites are presented with t = 0 hrs corresponding to the day of the formation of the closed circulation of the COL and on the right panels t = 0 hrs corresponds to the first time that ridging occurs.
To illustrate the differences between the circulation patterns when ridging highs occur with COLs and when they occur without COLs aloft, Fig. 2 shows the 340 K PV field represented by the thin black contours for the entire isentropic surface. The PV fields were chosen over the geopotential heights because of the former's better ability to show the undulations more clearly to illustrate wave breaking (as the contours turn back on themselves). As noted in Peters and Waugh (1996), this process is most clearly seen by examining PV on isentropic surfaces. Another added advantage for using PV is that the location of the closed PV anomaly relative to the dynamical tropopause is more readily observed, which is more intuitive than the geopotential heighs, in a dynamical sense. The evolution of the PV on the left and right panels of Fig. 2 shows that the amplitude of the trough (equatorward disturbance east of the thick black curve) grows more rapidly when there are COLs aloft, than when there are no COLs. In former case it actually grows into non-linear regimes, with the PV contours turning back on themselves indicating RWB. This is clearly indicated by the north-western/south-eastern orientation of the trough axis, meaning that the perturbation horizontal flow has a strongly negative zonal component in the midlatitudes, just upstream of the South African mainland (Fig. 2b). This demonstrates that there are significant differences in the circulation patterns aloft between the two scenarios.
As noted in Sect. 1, the overall aim of this study is to characterise downstream development in the South Atlantic Ocean/Southern African domain during ridging high pressure systems, by contrasting upper level dynamics when COLs occur aloft and when ridging events are not associated with COLs. Composites of the COL events were calculated for the sub-regions referred to in Sect. 3.2 but only the evolution of composite fields associated with south western COLs are shown in the left panels of all the figures presented in the paper. South western COLs were chosen because they have the highest frequency of occurrence than in the other three quadrants, as noted above. Furthermore, many high impact COLs (e.g. Thoithi et al. 2023) are first identified here and then propagate eastward to impact the eastern coast. In Fig. 3 the MSLP (shaded), geopotential height anomalies (thin solid and dashed black contours representing positive and negative anomalies, respectively) and PV = − 2.5 PVU (thick blue contours), PV = -3 PVU (thick red contours) and PV = − 3.5 PVU (thick black contours) on the 340 K isentropic surface, are shown from time lags t = -24 h to t = +48 h, in 24 h intervals. Note that we use PV on the 340 K isentropic surface in the case of south-western and south-eastern COLs and on the 360 K isentropic surface for the other two categories. Note also that the anomalies are calculated as the difference between the total variables and the 31-day mean centred on the day of the events. In using the 31-day mean as the basic state flow, this study follows previous studies (e.g. Orlanski and Katzfey 1991). Using daily climatologies of the variables would not change the patterns that have been observed in the study. During the evolution of the COLs (left panels in Fig. 3), the positive geopotential height anomalies are stronger across the ridge axes where there is clear PV overturning, indicating RWB. The negative geopotential height anomalies east of the PV overturning region are weaker but exhibit characteristics of a reflected Rossby wave, which were also observed in O' Brien and Reeder (2017) in the different context of PV based interactions of the waves with the jet stream. These observations apply to all categories of COLs (not shown). The wave breaking is deeper for southern COLs, in particular the ones upstream. This has implications for energy transfer as Gan and Piva (2013), Ndarana et al. (2021b) and Pinheiro et al. (2022) all showed that downstream development is most significant across the breaking ridge, because of the combined strength of the supergeostrophic flow (see Sect. 2 and Fig. 1) and the geopotential height anomalies found in that component of the baroclinic wave.
In contrast, the evolution of the baroclinic waves during ridging high pressure systems that do not have COLs associated with them (right panels of Fig. 3) have a stronger negative geopotental height anomalies across the trough axis (dashed contour in Fig. 3e) than the trailing positive geopotential height anomaly field, and attain absolute maximum values at time lag t = 0 hrs. This signals the inception of the ridging process. The waves that manifest under these circumstances are linear because they neither break nor is there any PV intrusion of significance. This difference in the characteristics of the baroclinic wave structure between instances when COLs occur and when there are no COLs during the ridging process is the first indication that there would be significant differences between downstream development processes associated with the two atmospheric situations. Whilst there are also differences at the surface, with the MSLP anomalies associated with south-western and south-eastern COLs showing much stronger values, all the waves exhibit a baroclinic structure, as they lag westward with height.

Dynamical processes associated with the wave trains
A place to start in explaining the possible dynamical processes that cause the differences described in the previous section are the meridional vertical structures of PV and geopotential height anomalies. These are shown in Fig. 4 as shaded fields and thin black contours, respectively, plotted with PV = -2 PVU (thick dashed blue contour) to indicate where the dynamical tropopause might be located. The PV The geopotential height anomalies are the circulation response to the PV anomalies, which is anticyclonic and cyclonic for positive and negative anomalies, respectively; by the invertibility principle (Hoskins et al. 1985;Davis and Emmanuel 1991;Røsting and Kristjánsson 2012). As in Fig. 3, the left panels show composite evolution of these fields for south-western COLs and the right panels show the evolution of fields associated with ridging events when there are no COLs aloft, plotted for the same time lags as in Fig. 3. Consistently with the deformation of the PV contours in Fig. 3, the south-western COLs exhibit higher positive PV anomalies, which can be seen by comparing the eddies labelled B in Fig. 4a and e. This in turn leads to deeper geopotential height perturbation fields that also more clearly extend to the surface for southern COLs. This is consistent with the stronger surface response that is observed for COLs (left panels of Fig. 3) noted previously. Also consistent with Fig. 3, is the fact that the fields associated with eddy B are stronger than those associated with eddy C, for all COL categories. The opposite is seen when ridging occurs without COLs. The right panels of Fig. 4, and in particular at t = 0 h, show that the negative anomaly downstream at C is stronger than the positive one at B. This is the case right through the evolution of ridging Areas where there no dots are significant at the 95% confidence level highs and most clearly seen at t = 0 h (Fig. 4f). Another difference to note between the left and right panels of Fig. 4 is that the COL PV perturbation fields are broader (larger zonal extent) than that which the upper level waves associated with ridging events with no COLs exhibit.
Therefore a stronger perturbation flow and associated eddy kinetic energy during COLs aloft can be expected.
To explain the likely cause of the breaking waves indicated by the overturning PV contours in Fig. 3b and the differences in the PV anomaly structures, we make use of  Figure 5 shows zonal isotachs (thin black contours) that are larger than 35 m s −1 to highlight jet streaks, the meridional perturbation velocity v = −3 m s −1 (thick red contours) and v = 3 m s −1 (thick blue contours), which, as indicated schematically in Fig. 1, show the structure of the baroclinic waves, with the ridge axes indicated by the thick solid black curve. Note that the baroclinic eddies of interest are labelled with Roman numerals I, II and III and that in the right panels there is a fourth component of the baroclinic wave labelled IV. The shading represents the magnitude of strain calculated using ( √ 2a) −1 | u + u tan | (Akahori and Yoden 1997). Early idealised modelling experiments of RWB processes (e.g. Nakamura and Plumb 1994) showed that increasing the rate of strain causes RWB to occur. As such, the left panels of Fig. 5 show that when COLs occur in the South African domain, the baroclinic waves break anticyclonically as the magnitude of strain increases and becomes largest at t = 0 hrs (Fig. 3b on the anticyclonic side of the jet. The strain rate, in turn is caused by the presence of the jet streaks which strides the ridge axes. Note that the axes have a pronounced north-western/south-eastern tilt, as would be expected in an anticyclonically breaking wave. In clear contrast, the corresponding jet streak in the right panels of Fig. 5 is much weaker, and the associated magnitude of strain is much weaker. It is therefore expected that the likelihood of RWB It is well known from classic jet streak theory (Keyser and Shapiro 1986) that the cyclonic jet entrance region is characterised by a negative vorticity anomaly field. So, the fact that the relatively stronger downstream jet, which, according to these composite calculations, is more pronounced when there are no COLs extends back to the trough axis means that it induces a stronger negative vorticity anomaly there, and its associated PV anomaly. This stems from the fact that the cyclonic entrance of a jet streak is characterised by a negative vorticity anomaly field. A stronger negative PV anomaly results in the part of the wave marked C in the right panels of Fig. 4.
The three dimensional view of the dynamical processes involved here that are presented in Figs. 3, 4 and 5 show that there are significant differences in upper tropospheric dynamical processes during ridging highs in the South African domain when COLs are present and they are not present. We may then expect that there would be differences in the energy generation and the transfer thereof and these will be discussed next.

The eddy kinetic budget
We now discuss the spatial structure of K e and consider the terms that contribute to its rate of change according to Equation (1). The spatial distribution of K e is shown in Fig. 6 as shaded fields, together with statistically significant v (thick grey) vectors, which indicate the statistical significance at the 95% confidence level of the K e fields, and the U = 40 m s −1 (very thick black contours) isotach to highlight the position of the jet streaks. Again, the eddies of interest are marked I, II and III in both the top left and right panels of Fig. 6 and eddy IV is shown in the latter. The dynamical processes discussed in the previous sections clearly have implications for the distribution and evolution of K e in Fig. 6.  Figs. 3, 4, and 5 all show that there are marked differences in the upper level dynamical characteristics during ridging high pressure systems where there are COLs and no COLs present. These differences also manifest in the eddy kinetic energy fields shown in Fig. 6. In particular, the perturbation PV fields induced by the more clearly defined RWB leads to a stronger perturbation flow for COLs located in the southern quadrants, so the eddy kinetic energy is relatively higher for these COLs as well. As it was the case with the other fields, there are marked differences between the left and right panels of Fig. 6. The stronger K e field is located in eddy I during COLs, but it is located in component III for the case of ridging events that have no COLs associated with them. This is caused by the differences in the PV anomaly structure presented in Fig. 4, in which the negative PV anomaly is stronger than in the COL case, thus inducing a stronger perturbation flow in eddy II (left panels of Fig. 6). Note also that K e IV also grows beyond t = 0 hrs. Ndarana et al. (2021b) showed that the kinetic energy in eddy I is converted from eddy available potential energy, as the jet streak (see left panels of Fig. 5) propagates across the South Atlantic Ocean, causing a strong eastward energy flux, thus translating the energy centre in the same direction of the movement of the midlatitude jet streak. Pinheiro et al. (2022) confirmed that this conversion is located upstream of the PV overturning region, and that it is stronger for strong COLs. Pinheiro et al. (2022) also noted that this conversion, which, according to Eq. 1., means that − > 0 ⇒ t K e > 0 , occurs in the thermally direct circulation region of jet streaks. Figure 7a shows that this is the case because according to these calculations, baroclinic conversion is concentrated in eddy I, which Figs. 5 and 6 show to be co-located with the jet entrance. Note that, as in previous figures, Fig. 7a is for south-western COLs. When COLs are not present, there are two jet streaks, but the one upstream is much weaker thus the weaker conversion values in eddies I and II in Fig. 7b. The presence of the relatively strong downstream jet in the right panels of Fig. 5 and also in the context of eddy kinetic energy in Fig. 6 leads to the baroclinic conversion in eddy III in Fig. 7b, which is weaker than that which is observed for COLs because the associated jet is also weaker. Thus the thermally direct transverse circulation would also be weaker. Note that in this case this process is concentrated beyond the inception of the ridging process in eddy II at around 40 • E, whereas in the COL case, it is located upstream of the Greenwich Meridian. Note also that Fig. 7b shows that there is some baroclinic conversion in eddies I and II, even though it is weak, and again, consistent with the weak upstream jet streak.
Baroclinic conversion constitutes one aspect of the pressure work term (see Eq. 2). The other aspect is the ageostrophic geopotential convergence term. As Ndarana et al. (2021b) and Pinheiro et al. (2022) note, the downstream transfer of K e by means of ageostrophic geopotential flux is most clearly seen across the ridge (i.e. thick black curve between eddies I and II, shown in Fig. 5). This curve corresponds to the black solid line which schematically shows the ridge axis in Fig. 1. Note that the orientation of the flux divergence/convergence and the flux vectors themselves is such that energy is transported from the midlatitudes into the South African domain. This energy is transferred from eddy I into II. This, of course, is due to the wave breaking which induces a strong supergeostrophic flow across the ridge axis. When the flow is combined with the strong geopotential anomalies (see Fig. 5 and the associated discussion in the previous sections) produces a strong ageostrophic geopotential flux. Because the waves do not break in the right panels of Fig. 8, the ageostrophic fluxes across the ridge are much weaker, but are stronger across the trough, where the subgeostrophic flow is stronger than the same for COLs. Again, the relatively stronger easterly ageostrophic flow combined with the relatively stronger negative anomalies, lead to a more clearly defined flux field, across the trough axis between II and III in the right panels of Fig. 8. For this reason, −∇ ⋅ (v ) a < 0 is stronger in II when there are no COLs in the South African domain, instead of eddy I. Note that there is also energy transfer across the ridge between III and IV, so that the latter grows by means of downstream development, as suggested by the growth of K e . Finally, Fig. 9 shows the Reynolds stress term, which may be decomposed into the barotropic and baroclinic conversions of K e to mean kinetic energy (Rivière et al. 2015), which signals the dissipation of the wave. This term is dominated by the former and is concentrated in eddy II for COLs but is located in component IV for the case of ridging highs without COLs. This is different from what Wirth et al. (2018) suggest. In that study, they proposed that baroclinic conversion occurs at the rear end of a wave packet.

Discussion of downstream development characteristics
In Sects. 4.1 and 4.2 it was shown that there are very clear differences in the flow across the ridge and trough axes The thick solid red and blue contours represent v = − 3 m s −1 and +3 m s −1 , resepctively, as in Fig. 4. The components of the baroclinic wave that propagate past the South African domain during ridging are labelled with Roman numerals I, II and III in the top panels, as in Fig. 4. In the right panels the eastern most eddy is labelled IV between the two categories of events. The supergeostrophic and subgeostrophic flows are respectively stronger for the cases of ridging with COLs and ridging without COLs. These differences are due to the PV anomaly structures that induce the flow which are associated with the jet streak structures. In turn, the characteristics of ageostrophic geopotential fluxes, and their divergence patterns, shown in Fig. 8, arise. Downstream development is best explained in terms of (1) the behaviour of K e of the different eddies of a baroblinic wave and (2) downstream energy transfer by means of ageostrophic geopotential fluxes ( v a ) and their divergence and convergence patterns.
Beginning with the evolution of K e , it is shown in Fig. 9a that the eddy kinetic energy I (green curve) grows in a parabolic fashion to maximise at about 24 h before the formation of the COLs, as found to be the case by Ndarana et al. (2021b). The energy grows by means of converting eddy available potential energy, as Fig. 7 demonstrates. The blue curve in Fig. 9a shows that the eddy kinetic energy II also grows, but starts gradually at first and then more steeply when the one upstream has saturated. Furthermore it grows to surpass the maximum that eddy kinetic energy attains and then saturates at t = 0 h, when the COLs form. The facts that (1) there is no significant baroclinic conversion associated with the eddy marked II in Fig. 7 and (2) there is no significant barotropic conversion in eddy I in Fig. 10 are together a clear indication that eddy II can only grow at the expense of eddy I, according to Eqs. (1) and (2), by means of transridge ageostrophic geopotential fluxes (Gan and Piva 2013;Ndarana et al. 2021b;Pinheiro et al. 2022), as suggested above. Note that the structure of the green curve (upstream eddy) vs. blue (downstream eddy) graphs in Fig. 9a is similar to Orlanski and Sheldon (1995) (cf. their Fig. 1a). The down turn in eddy kinetic energy II, beyond the formation of COLs, is caused by barotropic conversion (Fig. 10a).
The dynamical issues associated with the above, and the differences of these dynamics between COLs and ridging events that occur without COLs have already been discussed. For one thing, the supergeostrophic and subgeostrophic flow strengths are different because of the differences in the PV anomalies that induce it. The most significant downstream development processes occur downstream between eddies II and III across the trough, as the structure of the blue curve relative to the red one in Fig. 10b indicates. The fact that the growth of eddy I and II grow to maximise at lower values of K e compared to when there COLs aloft is consistent with the weaker baroclinic conversion found there. Eddy III also grows from baroclinic conversion, but dissipates through trans-ridge downstream development in the Indian Ocean, after the commencement of the ridging process. This downstream energy transfer gives rise to the growth of eddy IV, which does not have a counterpart during COLs. Eddy IV can only grow from energy transferred from eddy III because the former has no baroclinic conversion associated with it (Fig. 7b). It dissipates via barotropic processes as clearly indicated in Fig. 10e. It follows from this discussion that downstream development during COLs occurs upstream so that energy in transferred from the midlatitudes into the South African domain, whilst when ridging high occur without COLs aloft, it occurs further downstream as the associated baroclinic wave propagates past the South African domain.

Concluding remarks
Using ERA 5 reanalysis data and building on Gan and Piva (2013), Ndarana et al. (2021b) and Pinheiro et al. 2022, this study shown that the behaviour of energy transfer in the South Atlantic Ocean/South African/Indian Ocean sector differs between days when there are COLs affecting South Africa and when they do not form during ridging high pressure systems. Figure 11 summarises these findings. The use of composite analysis in a geographic context shows that very strong baroclinic conversion occurs west of the Greenwich Meridian during COL pressure systems (indicated by the letters VS-BCC in Fig. 11a), and it is much stronger  Fig. 7 Hovmöller diagrams of the composite meridional perturbation velocity v (thin black contours) and baroclinic conversion from eddy available potential energy to eddy kinetic energy ( − > 0 , shaded) plotted in 1 × 10 4 m 2 s −3 for (a) south-western, (b) ridging highs without COLs. The solid (dashed) black contours represent the positive (negative) v values plotted at 2 m s −1 contour intervals and the thick solid red and blue contours highlight v = -3 m s −1 and +3 m s −1 , resepctively, as in Fig. 4. The components of the baroclinic wave that propagate past the South African domain during ridging are labelled with Roman numerals I, II and III in the top panels, as in Fig. 4. In the bottom panel the eastern most eddy is labelled IV than when there are no COLs and it is found in eddy I (indicated by the letters W-BCC in Fig. 11d and e). By way of contrast, the weaker upstream baroclinic conversion found when no COLs occur during the ridging process but we also observe a slightly stronger conversion field in eddy II downstream of South Africa (indicated by the letters S-BCC in Fig. 11e), which appears to start occurring as the ridging process begins. Baroclinic conversion appears to generally occur at the rear end of the waves, which is consistent with Fig. 10 in Wirth et al. (2018). However, what this study has found is that barotropic conversion occurs downstream (i.e. eddy II in the case of COLs and eddy IV when there are no COLs, as indicated by the letters BTC ub Fig. 11c and f, respectively).
The differences in the energy transfer associated with these contrasting atmospheric conditions are due to the differences in the prevailing upper level dynamical processes. As Pinheiro et al. (2022) noted, baroclinic conversion occurs in the thermally direct region of the jet streaks embedded in the baroclinic waves. This jet streak is much stronger for COLs, hence the stronger baroclinic conversion. The presence of the baroclinic conversion downstream from South Africa in eddy III is informed by the presence of the downstream jet, which does not have a counterpart in that when COLs occur. Relative position composite fields in Ndarana and Waugh (2010), Ndarana et al. (2021b) and Pinheiro et al. (2022) suggest that it would be located over the South African mainland as suggested Fig. 11b. The stronger jet associated with COLs causes an increase the rate of strain, which (a) t = -24 hours (with southwestern COLs)  Fig. 8 Same as in Fig. 2 but for −∇ ⋅ (v ) a (shaded) plotted in 1 × 10 4 m 2 s −3 and the geopotential height anomalies are plotted in 30 gpm contour intervals and the arrows represent (v a ). The black arrows are statistically significant at the 95% confidence level then causes the waves to break anticyclonically (Nakamura and Plumb 1994), in the vicinity of the Greenwich Meridian. The breaking has no counterpart when COLs do not form, because the upstream jet streak is weak. The breaking and lack thereof then leads to stronger ageostrophic geopotential fluxes (as indicated by the double red curved red arrows in Fig. 11b), which are oriented so that they transport eddy kinetic energy in a north-easterly direction towards South Africa. The strong supergeostrophic flow is induced by the stronger geopotential height anomalies, which do not materialise during ridging highs, when COLs do not form. Even though the trans-trough axis (dashed line in Fig. 11e) downstream energy transport is weaker than when there is breaking, it is stronger than its upstream counterpart. This is largely due to the fact that the downstream jet in Fig. 11e marked "DJ" extends westward so that its entrance, and thus the negative relative vorticity anomaly at its entrance increases the negative anomalies of across the trough axis. This causes downstream propagation of the wave because it causes eddy III to grow, which itself dissipates by transferring energy into eddy IV by means of downstream development across the downstream ridge axis. The findings in this paper will have implications for medium range forecasting (MRF) time scale predictability studies of ridging highs and COLs in the South African domain. One way these findings could be implemented in the South African domain is to characterise downstream development in MRF ensemble systems to indicate whether or affirm that a COL is looming or not using the dynamical characteristics identified in this study. Author contributions TN, LEL and TSR conceptualised the research and TN prepared the schematics Figs. 1 and 11 and prepared the first draft of the manuscripts. LEL and TSR prepared the rest of the plots. All the authors contributed to the analysis and reviewed the manuscript.
Funding Open access funding provided by University of Pretoria.

Availability of data and materials
This study has been supported by the Water Research Foundation (Grant number: C2020/2021-00653) of South Africa and it is part of a larger project that aims to understand the underlying dynamical processes of the important synoptic systems in the South Atlantic Ocean/South Africa basin. All the data used in the study was obtained are the Fifth Generation European Centre for Medium-Range Weather Forecasts reanalysis (ERA5) can be obtained from. https:// www. ecmwf. int/ en/ forec asts/ datas ets/ reana lysis-datas ets/ era5 Declaration Conflict of interest The authors neither have conflict of interest nor competing interests.

Fig. 11
Schematic of K e during the evolution of COLs (left panels) and when there are no COLs (right panels) in the South African domain. The yellow and brown oval shapes represent the K e , with the brown shade representing higher values of K e . The dashed oval shapes represent the jet streaks. Areas marked W-BCC, S-BCC and VS-BCC indicate where weak, strong and very strong baroclinic conversion is located, respectively. The curved red arrows represent the ageostrophic geopotential fluxes and BTC represents where barotropic conversion is located. The eddies are labelled according to the schematic in Fig. 1 Ethics approval Not applicable.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http:// creat iveco mmons. org/ licen ses/ by/4. 0/.