We examine Saturn's atmosphere with observations from ground-based telescopes and Hubble Space Telescope (HST). We present a detailed analysis of observations acquired during 2018. A system of polar storms that appeared in the planet in March 2018 and remained active with a complex phenomenology at least until Sept. is analyzed elsewhere (Sanchez-Lavega et al., in press , 2019). Many of the cloud features in 2018 are long-lived and can be identified in images in 2017, and in some cases, for up to a decade using also Cassini ISS images. Without considering the polar storms, the most interesting long-lived cloud systems are: i) A bright spot in the EZ that can be tracked continuously since 2014 with a zonal velocity of 444 m/s in 2014 and 452 m/s in 2018. This velocity is different from the zonal winds at the cloud level at its latitude during the Cassini mission, and is closer to zonal winds obtained at the time of the Voyager flybys and zonal winds from Cassini VIMS infrared images of the lower atmosphere. ii) A large Anticyclone Vortex, here AV, that formed after the GWS of 2010-2011. This vortex has changed significantly in visual contrast, drift rate and latitude with minor changes in size over the last years. iii) A system of subpolar vortices present at least since 2011. These vortices follow drift rates consistent with zonal winds obtained by Cassini. We also present the positions of the vertices of the North polar hexagon from 2015 to 2018 compared with previous analyses during Cassini (2007-2014), observations with HST, and Voyager data in 1980-1981 to explore the long-term hexagon's drift rate. Variations in the drift rate cannot be fit by seasonal changes. Instead, the different drift rates reinforce the role of the North Polar Spot that was present in the Voyager epoch to cause a faster drift rate of the hexagon at that time compared with the current one.

Jupiter's banded structure undergoes strong temporal variations, changing the visible and infrared appearance of the belts and zones in a complex and turbulent way due to physical processes that are not yet understood. In this study we use ground-based 5-$\mu$m infrared data captured between 1984 and 2018 by 8 different instruments mounted on the Infrared Telescope Facility in Hawai'i and on the Very Large Telescope in Chile to analyze and characterize the long-term variability of Jupiter's cloud-forming region at the 1-4 bar pressure level. The data show a large temporal variability mainly at the equatorial and tropical latitudes, with a smaller temporal variability at mid-latitudes. We also compare the 5-$\mu$m-bright and -dark regions with the locations of the visible zones and belts and we find that these regions are not always co-located, specially in the southern hemisphere. We also present Lomb-Scargle and Wavelet Transform analyzes in order to look for possible periodicities of the brightness changes that could help us understand their origin and predict future events. We see that some of these variations occur periodically in time intervals of 4-8 years. The reasons of these time intervals are not understood and we explore potential connections to both convective processes in the deeper weather layer and dynamical processes in the upper troposphere and stratosphere. Finally we perform a Principal Component analysis to reveal a clear anticorrelation on the 5-$\mu$m brightness changes between the North Equatorial Belt and the South Equatorial Belt, suggesting a possible connection between the changes in these belts.

We analyse the potential vorticity of Saturn's polar regions, as it is a fundamental dynamical tracer that enables us to improve our understanding of the dynamics of these regions and their seasonal variability. In particular, we present zonally averaged quasi-geostrophic potential vorticity maps between 68{\deg} planetographic latitude and the poles at altitudes between 500 mbar and 1mbar for three different epochs: (i) June 2013 (early northern summer) for the north polar region, (ii) December 2008 (late northern winter) for both polar regions and (iii) October 2006 (southern summer) for the south, computed using temperature profiles retrieved from Cassini Composite Infrared Spectrometer (CIRS) data and wind profiles obtained from Cassini's Imaging Science Subsystem (ISS). The results show that quasi-geostrophic potential vorticity maps are very similar at all the studied epochs, showing positive vorticities at the north and negative at the south, indicative of the dominance of the Coriolis parameter 2${\Omega}sin{\phi}$ at all latitudes, except near the pole. The meridional gradients of the quasi-geostrophic potential vorticity show that dynamical instabilities, mainly due to the barotropic term, could develop at the flanks of the Hexagon at 78{\deg}N, the jet at 73.9{\deg}S and on the equatorward flank of both polar jets. There are no differences in potential vorticity gradients between the two hemispheres that could explain why a hexagon forms in the north and not in the south. No seasonal variability of the potential vorticity and its meridional gradient has been found, despite significant changes in the atmospheric temperatures over time.