to} \approx 0.1 r_o/(10^{-2} U_0) \approx 10 Ro^{-1} \Omega^{-1} \approx 10^{3}\Omega^{-1}$: the turnover time of the particle is much longer than the timescale of the wave. Consequently the particle oscillates rapidly as the wave propagates and is slowly advected by the zonal flow. In practice the radial displacement $l$ is typically smaller than $\delta$ and so the turnover time is slightly overestimated here.

#### Supercritical regime

When the Rossby number is increased in the supercritical regime, other prograde jets will eventually become unstable. A second Rossby wave appears in the weakly supercritical regime, at $Ro=1.06 Ro_c$ for $E=5\times 10^{-6}$, with a maximum velocity located in the retrograde zonal jet at larger radius than the first wave maxima ([*i.e.*]{} the wave appearing for $Ro=Ro_c$) (figure \[fig:U\_E5e6R3600\]). To fill the larger circumference at larger radius the instability has a slightly larger wave number, $m=22$ instead of $21$, while the radial width of the jet is comparable. The second wave propagates faster, in agreement with the Rossby wave dispersion relation (\[eq:w\_RW\]). Barotropic instabilities tend to broaden and weaken narrow jets by redistributing potential vorticity [see for instance @Ped87]. The smoothing of the jets saturates the amplitude of the Rossby waves. For this slightly supercritical regime the zonal flow profile is only weakly modified. Upon further increasing the forcing ($Ro=2.94 Ro_c$), several Rossby waves of different wavenumbers superpose and interact (figure \[fig:U\_E5e6R10000\]). The structure of the waves and the jets is still mainly bidimensional except in the viscous boundary layers. The typical cylindrical radial velocity is $V_s \approx 0.1 U_0$ and the Rossby number is $0.05$ so the turnover time is about $20\Omega^{-1}$ assuming that the radial displacement $l=\delta$, about the same order of magnitude as the timescale of the zonal jets.

In figure \[fig:U\_s\_satur\] the time-averaged zonal flow in the equatorial plane is plotted for different $Ro$ up to $
to } \approx 0.1 r_o/(10^{-2 } U_0) \approx 10 Ro^{-1 } \Omega^{-1 } \approx 10^{3}\Omega^{-1}$: the turnover time of the particle is much long than the timescale of the wave. Consequently the atom oscillates rapidly as the wave propagates and is slowly advect by the zonal flow. In exercise the radial displacement $ l$ is typically smaller than $ \delta$ and indeed the dollar volume time is slightly overestimated here. 

 # # # # Supercritical regime 

 When the Rossby issue is increased in the supercritical regime, early prograde jets will eventually become unstable. A second Rossby wave appears in the weakly supercritical regime, at $ Ro=1.06 Ro_c$ for $ E=5\times 10^{-6}$, with a maximum speed located in the retrograde zonal jet at larger spoke than the first wave maxima ([ * i.e. * ] { }   the wave appearing for $ Ro = Ro_c$) (figure   \[fig: U\_E5e6R3600\ ]). To fill the larger circumference at bigger radius the instability has a slightly larger wave numeral, $ m=22 $ instead of $ 21 $, while the radial width of the jet is comparable. The second wave propagates faster, in agreement with the Rossby wave dispersion relation   (\[eq: w\_RW\ ]). Barotropic instabilities tend to broaden and weaken narrow jets by redistributing potential vorticity [ see for instance @Ped87 ]. The smoothing of the jets saturates the amplitude of the Rossby wave. For this slightly supercritical regime the zonal menstruation profile is only weakly modified. Upon further increase the forcing ($ Ro=2.94 Ro_c$), several Rossby waves of different wavenumbers superpose and interact (design   \[fig: U\_E5e6R10000\ ]). The structure of the waves and the jets is still mainly bidimensional except in the gluey boundary layers. The typical cylindrical radial velocity is $ V_s \approx 0.1 U_0 $ and the Rossby number is $ 0.05 $ so the upset time is about $ 20\Omega^{-1}$ assume that the radial translation $ l=\delta$, about the same order of magnitude as the timescale of the zonal k. 

 In figure   \[fig: U\_s\_satur\ ] the time - averaged zonal flow in the equatorial plane is plotted for different $ Ro$ up to $
to} \wpprox 0.1 r_o/(10^{-2} U_0) \approx 10 Ro^{-1} \Umega^{-1} \approx 10^{3}\Omggq^{-1}$: the vurnoved time ow the particle is much longec thqn tht timescale of the dave. Consvquently rhe karticle oscillatxa rapidly as fme ware propagates anc is slowlf advected by dhd vonal flow. In practice the radial difplacemrnh $l$ is typicaljy skwlled than $\delta$ and so the turnover tjme is vlightly overrstimated here.

#### Supercriticwl rfgime

When the Rosshy number iw inseeased in thd supercriuieal regime, kther prograde jets will eventuxlly yecome unstqboe. W second Roswby wwve appears lm the feakly xupercritical vegimx, at $Ro=1.06 Ro_c$ for $E=5\times 10^{-6}$, wmth a maximum velocijy located iu the retrograde zonao het aj larcer fqdijs uhai tge firdt xave maxima ([*i.e.*]{} the wave appearing for $Ro=Ro_v$) (spbure \[fig:U\_E5e6R3600\]). Fo filj ehe larger circumference at larger radils tge instability has a slughtly larger wave nulber, $m=22$ inftead of $21$, while the radial width of the jet is cokparauld. Tkc rwclnd wave propagates faster, in agreement with ege Rpssby wave disiersion relation (\[eq:e\_RA\]). Nwrotropic insjabilitnss tend to broaden ajd weakgn nareow jets fy rrdistributing potential vorricity [see fje instance @Ped87]. The smoothing uf tne jeys saturates the amplitbde of the Rossby waves. Fod this slightly sjpegcridical regime the zonal flor profile is pnly wexkly modifyed. Upon fkrthev increasing the fogcing ($Rm=2.94 Ro_c$), sevegal Rossby waves of different wetenumbers supgrpmse and intzract (nigure \[fig:U\_E5e6R10000\]). Tre structure oy the wares ana the jets is stinl mainly bydimensional adcept in the viscous boubdart layerr. The typical culindrical radial vwlocity is $V_s \apprpx 0.1 H_0$ and the Rossbv bumber is $0.05$ so tne gurgoner tyke is about $20\Mmegx^{-1}$ arxumine that uke vadkal cisplacement $l=\delta$, dbouf the same order og iagnitudg as the eimescale of yhe zonal jets.

In flgure \[hig:U\_s\_setur\] tne jime-averaged zonal flow in the squatoriap pjane is ploteed nor qifferent $Xo$ up to $
to} \approx 0.1 r_o/(10^{-2} U_0) \approx 10 \approx the turnover of the particle timescale the wave. Consequently particle oscillates rapidly the wave propagates and is slowly by the zonal flow. In practice the radial displacement $l$ is typically smaller $\delta$ and so the turnover time is slightly overestimated here. #### Supercritical regime the number increased the supercritical regime, other prograde jets will eventually become unstable. A second Rossby wave appears in weakly supercritical regime, at $Ro=1.06 Ro_c$ for $E=5\times with a maximum velocity in the retrograde zonal jet larger than the wave ([*i.e.*]{} wave appearing for (figure \[fig:U\_E5e6R3600\]). To fill the larger circumference at larger radius the instability has a slightly larger wave $m=22$ instead while the width the is comparable. The propagates faster, in agreement with the relation (\[eq:w\_RW\]). Barotropic instabilities tend to broaden and narrow jets redistributing potential vorticity [see for instance The smoothing of the jets saturates the amplitude the Rossby waves. For this slightly supercritical regime the zonal flow profile is only weakly further increasing the forcing Ro_c$), several Rossby of wavenumbers and (figure \[fig:U\_E5e6R10000\]). structure of the waves and the jets is still mainly bidimensional in the viscous boundary layers. The typical cylindrical radial velocity \approx U_0$ and the number is $0.05$ so turnover is about $20\Omega^{-1}$ assuming radial $l=\delta$, order magnitude the timescale of the jets. In figure \[fig:U\_s\_satur\] the zonal flow in the different $Ro$ up to $
to} \approx 0.1 r_o/(10^{-2} U_0) \approx 10 Ro^{-1} \Omega^{-1} \aPprox 10^{3}\Omega^{-1}$: The tuRnoVer TiMe of The pArticle is much lONger Than the timescale of the wAve. CoNsEQuenTLy The paRticle oSCiLLAteS rApIdlY aS ThE wave ProPagates And is slowlY adVeCted by the zonAL fLow. In practIce The radial disPlaCement $L$ iS tyPIcallY smAller Than $\deLTa$ and sO the turnoVeR Time is SLightly OVErEstiMated here.

#### SupercriTIcAL regime

When the rossby NuMBeR IS inCreAsed in the sUpErcriTIcal regIMe, OTHEr pROgrade jets wilL eventually BEcoMe unstAbLe. A SEcond ROssby WaVE apPears in the wEaklY supercriTical rEGime, at $RO=1.06 ro_c$ for $E=5\Times 10^{-6}$, wIth A maXimuM VeLoCitY lOCatED iN thE RetRograde zOnAl Jet at LargER RADius ThaN the First Wave maxima ([*i.e.*]{} tHe wAve aPPeaRing fOr $Ro=RO_c$) (fiGuRe \[fig:u\_E5e6R3600\]). To Fill tHe Larger circumferEnce At larger rAdiUs The InStabiLIty has A slIghTly largEr wave nUMbeR, $m=22$ INSTeAd of $21$, while the radial WiDTH oF the jet iS compaRAbLe. tHe second WaVe pRopaGATes faSter, IN aGreement With thE roSsBy wave dIsPersioN rElaTioN (\[eq:w\_Rw\]). baroTropic InstabilIties TEnd to broaden anD Weaken narrow jETs BY ReDIstrIbuTing potentiAl voRTiciTy [seE FoR inSTance @ped87]. ThE sMOoTHing of the jets saturaTeS the amPlituDe of the Rossby Waves. For thIS SLightly sUperCRiTIcal regime the zOnal fLow profile IS only weaKly moDified. UpOn further INCreasing The ForCinG ($Ro=2.94 rO_C$), sEveral Rossby wAVEs of DiFferent WavEnumberS suPerPosE anD iNteract (fiGure \[fig:U\_e5e6r10000\]). THe StRucTure oF The waves AnD thE jEts Is stiLL mainlY bidiMensIoNaL ExcEpt in thE ViSCOus bOuNdAry lAyeRs. the tyPicaL CylIndricaL radial veLocITy is $v_s \ApProx 0.1 U_0$ anD the Rossby numBeR is $0.05$ so the tuRnOveR time iS ABout $20\OmegA^{-1}$ assuming that the radial dISplacemEnt $L=\deltA$, aboUt the same OrdEr of maGniTUde as tHe timeScale Of The ZONal jeTS.

in FigUrE \[fig:U\_s\_satuR\] THe tIme-avErAged Zonal flOw in the equatorial pLAne Is plotted for dIffErenT $rO$ uP to $
to} \approx 0.1 r_o/(10^{- 2} U_0) \a pprox  10  Ro ^{ -1}\Ome ga^{-1} \appro x  10^ {3}\Omega^{-1}$: the t urnov er time of  theparticl e  i s muc hlo nge rt ha n the  ti mescale  of the wa ve.  C onsequentlyt he  particleosc illates rapi dly  as th ewav e  prop aga tes a nd iss lowlyadvectedby the zo n al flow . In  pra ctice the radiald is p lacement $l$ i s typi ca l ly s mal ler  than $\de lt a$ an d  so the tu r n o ver time is sligh tly overest i mat ed her e. 

# # ## Sup ercri ti c alregime

Whe n th e Rossbynumber is incr e ased in  the s upe rcr itic a lre gim e, oth e rpro g rad e jets w il levent uall y b e come  un stab le. A  second Rossb y w avea ppe ars i n the  wea kl y sup ercrit icalre gime, at $Ro=1. 06 R o_c$ for$E= 5\ tim es  10^{ - 6}$, w ith  amaximum  veloci t y l oc a t e din the retrogradezo n a ljet at l argerr ad iu s  than th efir st w a v e max ima( [* i.e.*]{}  the w a ve  a ppearin gfor $R o= Ro_ c$)  (fig u re \ [fig:U \_E5e6R3 600\] ) . To fill thel arger circumf e re n c ea t la rge r radius th e in s tabi lity ha s a sligh tly l ar g er wave number, $m=22$  i nstead  of $ 21$, while th e radial w i d t h of the  jet is comparable. Th e sec ond wave p r opagates  fast er, in a greementw i th the R oss bywav e d i s pe rsion relatio n (\[e q: w\_RW\] ).Barotro pic  in sta bil it ies tendto broad en  a nd  w eak en na r row jets  b y r ed ist ribut i ng pot entia l vo rt ic i ty[see fo r  i n s tanc e@P ed87 ].Th e smo othi n g o f the j ets satur ate s  the  a mp litudeof the Rossby  w aves. Forth isslight l y  supercr itical regime the zonal flow pr ofi le is  onl y weaklymod ified.  Up o n furt her in creas in g t h e  forc i n g($R o= 2.94 Ro_c$ ) ,  se veral  R ossb y waves  of different wave n umb ers superpose  an d in t e ra ct( fi g ure  \ [ fig : U \_E5e6R10000\]) . The stru ct u re  of the wa v esan d the j ets isstill mainlybidimensi onal exce pt  int h e v iscous bou ndary la yers. The typic a lcylin dri cal ra di alveloc ity is $V_ s \ap prox 0 .1  U_0$and t he  Rossbynumber is $0.05$ so the  turno ver t ime  is about  $2 0 \Om ega^{-1}$  ass uming that  th e r adial  di s place ment $l =\d e lta$,  abo u t the sam e  o rde r of  magnitudea s the  time sca l e of t he z onal jets.

In fi g ure \[fig:U\_s \_sa t u r\]  th e  tim e- averaged zonal  fl ow i n the eq ua torial plan e is plo tt e d for  diffe rent $ Ro$ upt o  $ 
to} \approx_0.1 r_o/(10^{-2}_U_0) \approx 10 Ro^{-1}_\Omega^{-1} \approx_10^{3}\Omega^{-1}$:_the turnover_time_of the particle_is much longer_than the timescale of_the wave. Consequently_the_particle oscillates rapidly as the wave propagates and is slowly advected by the zonal_flow._In practice_the_radial_displacement $l$ is typically smaller_than $\delta$ and so the_turnover time_is slightly overestimated here.

#### Supercritical regime

When the Rossby_number_is increased in_the supercritical regime, other prograde jets will eventually become_unstable. A second Rossby wave appears_in the weakly_supercritical_regime,_at $Ro=1.06 Ro_c$ for_$E=5\times 10^{-6}$, with a maximum velocity_located in the retrograde zonal jet_at larger radius than the first wave_maxima ([*i.e.*]{} the wave appearing for $Ro=Ro_c$)_(figure \[fig:U\_E5e6R3600\]). To fill the larger_circumference at_larger radius the instability has_a slightly larger_wave number,_$m=22$ instead of_$21$, while the radial width of_the jet is_comparable. The second wave propagates faster,_in_agreement with the_Rossby_wave_dispersion relation (\[eq:w\_RW\])._Barotropic instabilities tend_to_broaden and_weaken_narrow jets by redistributing potential vorticity_[see_for instance @Ped87]. The smoothing of the_jets saturates the amplitude_of_the Rossby waves. For_this slightly supercritical regime the_zonal flow profile is only weakly_modified. Upon_further increasing_the forcing ($Ro=2.94 Ro_c$), several Rossby waves of different wavenumbers superpose_and interact (figure \[fig:U\_E5e6R10000\]). The structure of_the waves and the_jets is_still_mainly bidimensional except_in_the viscous_boundary layers. The typical cylindrical radial velocity_is $V_s_\approx 0.1 U_0$ and the Rossby_number is $0.05$ so_the_turnover time is about $20\Omega^{-1}$ assuming_that the radial displacement $l=\delta$, about_the same order of magnitude_as_the_timescale of the zonal jets.

In_figure \[fig:U\_s\_satur\] the time-averaged zonal flow in_the equatorial plane_is plotted for different $Ro$ up to_$