Uranus and Neptune, and their diverse satellite and ring systems, represent the least explored environments of our Solar System, and yet may provide the archetype for the most common outcome of planetary formation throughout our galaxy. Ice Giants will be the last remaining class of Solar System planet to have a dedicated orbital explorer, and international efforts are under way to realise such an ambitious mission in the coming decades. In 2019, the European Space Agency released a call for scientific themes for its strategic science planning process for the 2030s and 2040s, known as Voyage 2050. We used this opportunity to review our present-day knowledge of the Uranus and Neptune systems, producing a revised and updated set of scientific questions and motivations for their exploration. This review article describes how such a mission could explore their origins, ice-rich interiors, dynamic atmospheres, unique magnetospheres, and myriad icy satellites, to address questions at the heart of modern planetary science. These two worlds are superb examples of how planets with shared origins can exhibit remarkably different evolutionary paths: Neptune as the archetype for Ice Giants, whereas Uranus may be atypical. Exploring Uranus' natural satellites and Neptune's captured moon Triton could reveal how Ocean Worlds form and remain active, redefining the extent of the habitable zone in our Solar System. For these reasons and more, we advocate that an Ice Giant System explorer should become a strategic cornerstone mission within ESA's Voyage 2050 programme, in partnership with international collaborators, and targeting launch opportunities in the early 2030s.

Among the four known transneptunian dwarf planets, Haumea is an exotic, very elongated, and fast rotating body. In contrast to the other dwarf planets, its size, shape, albedo, and density are not well constrained. Here we report results of a multi-chord stellar occultation, observed on 2017 January 21. Secondary events observed around the main body are consistent with the presence of a ring of opacity 0.5, width 70 km, and radius 2,287$_{-45}^{+75}$ km. The Centaur Chariklo was the first body other than a giant planet to show a ring system and the Centaur Chiron was later found to possess something similar to Chariklo's rings. Haumea is the first body outside the Centaur population with a ring. The ring is coplanar with both Haumea's equator and the orbit of its satellite Hi'iaka. Its radius places close to the 3:1 mean motion resonance with Haumea's spin period. The occultation by the main body provides an instantaneous elliptical limb with axes 1,704 $\pm$ 4 km x 1,138 $\pm$ 26 km. Combined with rotational light-curves, it constrains Haumea's 3D orientation and its triaxial shape, which is inconsistent with a homogeneous body in hydrostatic equilibrium. Haumea's largest axis is at least 2,322 $\pm$ 60 km, larger than thought before. This implies an upper limit of 1,885 $\pm$ 80 kg m$^{-3}$ for Haumea's density, smaller and less puzzling than previous estimations, and a geometric albedo of 0.51 $\pm$ 0.02, also smaller than previous estimations. No global N$_2$ or CH$_4$ atmosphere with pressures larger than 15 and 50 nbar (3-$\sigma$ limits), respectively, is detected.

Planetary-scale waves are thought to play a role in powering the yet-unexplained atmospheric superrotation of Venus. Puzzlingly, while Kelvin, Rossby and stationary waves manifest at the upper clouds (65--70 km), no planetary-scale waves or stationary patterns have been reported in the intervening level of the lower clouds (48--55 km), although the latter are probably Lee waves. Using observations by the Akatsuki orbiter and ground-based telescopes, we show that the lower clouds follow a regular cycle punctuated between 30$^{\circ}$N--40$^{\circ}$S by a sharp discontinuity or disruption with potential implications to Venus's general circulation and thermal structure. This disruption exhibits a westward rotation period of $\sim$4.9 days faster than winds at this level ($\sim$6-day period), alters clouds' properties and aerosols, and remains coherent during weeks. Past observations reveal its recurrent nature since at least 1983, and numerical simulations show that a nonlinear Kelvin wave reproduces many of its properties.

On 7th August 2019, an impact flash lasting $\sim1$s was observed on Jupiter. The video of this event was analysed to obtain the lightcurve and determine the energy release and initial mass. We find that the impactor released a total energy of $96-151$ kilotons of TNT, corresponding to an initial mass between $190-260$ metric tonnes with a diameter between $4-10$m. We developed a fragmentation model to simulate the atmospheric breakup of the object and reproduce the lightcurve. We model three different materials: cometary, stony and metallic at speeds of $60$, $65 $ and $70$ km/s to determine the material makeup of the impacting object. The slower cases are best fit by a strong, metallic object while the faster cases require a weaker material.

The Europlanet-2020 programme, which ended on Aug 31st, 2019, included an activity called VESPA (Virtual European Solar and Planetary Access), which focused on adapting Virtual Observatory (VO) techniques to handle Planetary Science data. This paper describes some aspects of VESPA at the end of this 4-years development phase and at the onset of the newly selected Europlanet-2024 programme starting in 2020. The main objectives of VESPA are to facilitate searches both in big archives and in small databases, to enable data analysis by providing simple data access and online visualization functions, and to allow research teams to publish derived data in an interoperable environment as easily as possible. VESPA encompasses a wide scope, including surfaces, atmospheres, magnetospheres and planetary plasmas, small bodies, helio-physics, exoplanets, and spectroscopy in solid phase. This system relies in particular on standards and tools developed for the Astronomy VO (IVOA) and extends them where required to handle specificities of Solar System studies. It also aims at making the VO compatible with tools and protocols developed in different contexts, for instance GIS for planetary surfaces, or time series tools for plasma-related measurements. An essential part of the activity is to publish a significant amount of high-quality data in this system, with a focus on derived products resulting from data analysis or simulations.

Moist convective storms in Jupiter develop frequently and can trigger atmospheric activity of different scales, from localized storms to planetary-scale disturbances including convective activity confined inside a larger meteorological system. In February 2018 a series of convective storms erupted in Jupiter's South Temperate Belt (STB) (planetocentric latitudes from -23$^{\circ}$ to -29.5$^{\circ}$). This occurred inside an elongated cyclonic region known popularly as the STB Ghost, close to the large anticyclone Oval BA, resulting in the clouds from the storms being confined to the cyclone. The initial storms lasted only a few days but they generated abundant enduring turbulence. They also produced dark features, possibly partially devoid of clouds, that circulated around the cyclone over the first week. The subsequent activity developed over months and resulted in two main structures, one of them closely interacting with Oval BA and the other one being expelled to the west. Here we present a study of this meteorological activity based on daily observations provided by the amateur community, complemented by observations obtained from PlanetCam UPV/EHU at Calar Alto Observatory, the Hubble Space Telescope and by JunoCam on the Juno spacecraft. We also perform numerical simulations with the EPIC General Circulation Model to reproduce the phenomenology observed. The successful simulations require a complex interplay between the Ghost, the convective eruptions and Oval BA, and they demonstrate that water moist convection was the source of the initial storms. A simple scale comparison with other moist convective storms that can be observed in the planet in visible and methane absorption band images strongly suggests that most of these storms are powered by water condensation instead of ammonia.

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.

Remote sensing observations suffer significant limitations when used to study the bulk atmospheric composition of the giant planets of our solar system. This impacts our knowledge of the formation of these planets and the physics of their atmospheres. A remarkable example of the superiority of in situ probe measurements was illustrated by the exploration of Jupiter, where key measurements such as the determination of the noble gases' abundances and the precise measurement of the helium mixing ratio were only made available through in situ measurements by the Galileo probe. Here we describe the main scientific goals to be addressed by the future in situ exploration of Saturn, Uranus, and Neptune, placing the Galileo probe exploration of Jupiter in a broader context. An atmospheric entry probe targeting the 10-bar level would yield insight into two broad themes: i) the formation history of the giant planets and that of the Solar System, and ii) the processes at play in planetary atmospheres. The probe would descend under parachute to measure composition, structure, and dynamics, with data returned to Earth using a Carrier Relay Spacecraft as a relay station. An atmospheric probe could represent a significant ESA contribution to a future NASA New Frontiers or flagship mission to be launched toward Saturn, Uranus, and/or Neptune.

The Cassini spacecraft made a gravity assist flyby of Jupiter in December 2000. The Imaging Science Subsystem (ISS) acquired images of the planet that covered the visual range with filters sensitive to the distribution of clouds and hazes, their altitudes and color. We use a selection of these images to build high-resolution cylindrical maps of the planet in 9 wavelengths. We explore the spatial distribution of the planet reflectivity examining the distribution of color and altitudes of hazes as well as their relation. A variety of analyses is presented: a) Principal Component Analysis (PCA); b) color-altitude indices; and c) chromaticity diagrams (for a quantitative characterization of Jupiter "true" colors as they would be perceived by a human observer). PCA of the full dataset indicates that six components are required to explain the data. These components are likely related to the distribution of cloud opacity at the main cloud, the distribution of two types of hazes, two chromophores or coloring processes and the distribution of convective storms.

The Venusian atmosphere is covered by clouds with super-rotating winds whose accelerating mechanism is still not well understood. The fastest winds, occurring at the cloud tops ($\sim$70 km height), have been studied for decades thanks to their visual contrast in dayside ultraviolet images. The middle clouds ($\sim$50-55 km) can be observed at near-infrared wavelengths (800-950 nm), although with very low contrast. Here we present the first extensive analysis of their morphology and motions at lower latitudes along 2016 with 900-nm images from the IR1 camera onboard Akatsuki. The middle clouds exhibit hemispherical asymmetries every 4-5 days, sharp discontinuities in elongated "hook-like" stripes, and large contrasts (3-21%) probably associated with large changes in the optical thickness. Zonal winds obtained with IR1 images and with ground-based observations reveal mean zonal winds peaking at the equator, while their combination with Venus Express unveils long-term variations of 20 m s$^{-1}$ along 10 years.

We present measurements of the wind speeds at the nightside lower clouds of Venus from observations by JAXA's mission Akatsuki during 2016, complemented with new wind measurements from ground-based observations acquired with TNG/NICS in 2012 and IRTF/SpeX in 2015 and 2017. Zonal and meridional components of the winds were measured from cloud tracking on a total of 466 Akatsuki images of Venus acquired by the camera IR2 using the 2.26-$\mathrm{\mu m}$ filter, with spatial resolutions ranging 10--80 km per pixel and covering from 2016 March 22 to October 31. More than 149,000 wind vectors were obtained with an automatic technique of template matching, and 2,947 wind vectors were inferred with the manual procedure. The meridional profiles for both components of the winds are found to be consistent with results from the Venus Express mission during 2006--2008, although stronger wind variability is found for the zonal component at equatorial latitudes where Akatsuki observations have better viewing geometry than Venus Express. The zonal winds at low latitudes also suggest a zonal variability that could be associated with solar tides or vertically propagating orographic waves. Finally, the combination of our wind measurements from TNG/NICS, IRTF/SpeX and Akatsuki images with previously published and based in data from 1978 to 2017 suggests variations of up to 30 m s$^{-1}$ in the winds at the lower clouds of the Venus nightside.

We report the discovery of a large ($\sim$8500 km diameter) infrared-bright storm at Neptune's equator in June 2017. We tracked the storm over a period of 7 months with high-cadence infrared snapshot imaging, carried out on 14 nights at the 10 meter Keck II telescope and 17 nights at the Shane 120 inch reflector at Lick Observatory. The cloud feature was larger and more persistent than any equatorial clouds seen before on Neptune, remaining intermittently active from at least 10 June to 31 December 2017. Our Keck and Lick observations were augmented by very high-cadence images from the amateur community, which permitted the determination of accurate drift rates for the cloud feature. Its zonal drift speed was variable from 10 June to at least 25 July, but remained a constant $237.4 \pm 0.2$ m s$^{-1}$ from 30 September until at least 15 November. The pressure of the cloud top was determined from radiative transfer calculations to be 0.3-0.6 bar; this value remained constant over the course of the observations. Multiple cloud break-up events, in which a bright cloud band wrapped around Neptune's equator, were observed over the course of our observations. No "dark spot" vortices were seen near the equator in HST imaging on 6 and 7 October. The size and pressure of the storm are consistent with moist convection or a planetary-scale wave as the energy source of convective upwelling, but more modeling is required to determine the driver of this equatorial disturbance as well as the triggers for and dynamics of the observed cloud break-up events.

Small-scale waves were observed along the boundary between Jupiter's North Equatorial Belt and North Tropical Zone, ~16.5{\deg} N planetographic latitude in Hubble Space Telescope data in 2012 and throughout 2015 to 2018, observable at all wavelengths from the UV to the near IR. At peak visibility, the waves have sufficient contrast (~10%) to be observed from ground-based telescopes. They have a typical wavelength of about 1.2{\deg} (1400 km), variable-length wave trains, and westward phase speeds of a few m/s or less. New analysis of Voyager 2 data shows similar wave trains over at least 300 hours. Some waves appear curved when over cyclones and anticyclones, but most are straight, but tilted, shifting in latitude as they pass vortices. Based on their wavelengths, phase speeds, and faint appearance at high-altitude sensitive passbands, the observed NEB waves are consistent with inertia-gravity waves at the 500-mbar pressure level, though formation altitude is not well constrained. Preliminary General Circulation Model simulations generate inertia-gravity waves from vortices interacting with the environment and can reproduce the observed wavelengths and orientations. Several mechanisms can generate these waves, and all may contribute: geostrophic adjustment of cyclones; cyclone/anticyclone interactions; wind interactions with obstructions or heat pulses from convection; or changing vertical wind shear. However, observations also show that the presence of vortices and/or regions of convection are not sufficient by themselves for wave formation, implying that a change in vertical structure may affect their stability, or that changes in haze properties may affect their visibility.

We characterise the origin and evolution of a mesoscale wave pattern in Jupiter's North Equatorial Belt (NEB), detected for the first time at 5 $\mu$m using a 2016-17 campaign of `lucky imaging' from the VISIR instrument on the Very Large Telescope and the NIRI instrument on the Gemini observatory, coupled with M-band imaging from Juno's JIRAM instrument during the first seven Juno orbits. The wave is compact, with a $1.1-1.4^\circ$ longitude wavelength (wavelength 1,300-1,600 km, wavenumber 260-330) that is stable over time, with wave crests aligned largely north-south between $14$ and $17^\circ$N (planetographic). The waves were initially identified in small ($10^\circ$ longitude) packets immediately west of cyclones in the NEB at $16^\circ$N, but extended to span wider longitude ranges over time. The waves exhibit a 7-10 K brightness temperature amplitude on top of a $\sim210$-K background at 5 $\mu$m. The thermal structure of the NEB allows for both inertio-gravity waves and gravity waves. Despite detection at 5 $\mu$m, this does not necessarily imply a deep location for the waves, and an upper tropospheric aerosol layer near 400-800 mbar could feature a gravity wave pattern modulating the visible-light reflectivity and attenuating the 5-$\mu$m radiance originating from deeper levels. Strong rifting activity appears to obliterate the pattern, which can change on timescales of weeks. The NEB underwent a new expansion and contraction episode in 2016-17 with associated cyclone-anticyclone formation, which could explain why the mesoscale wave pattern was more vivid in 2017 than ever before.

Video observations of Jupiter obtained by amateur astronomers over the past eight years have shown five flashes of light of 1-2 s. The first three of these events occurred on 3 June 2010, 20 August 2010, and 10 September 2012. Previous analyses showed that they were caused by the impact of objects of 5-20 m in diameter, depending on their density, with a released energy comparable to superbolides on Earth of the class of the Chelyabinsk airburst. The most recent two flashes on Jupiter were detected on 17 March 2016 and 26 May 2017 and are analyzed here. We characterize the energy involved together with the masses and sizes of the objects that produced these flashes. The rate of similar impacts on Jupiter provides improved constraints on the total flux of impacts on the planet, which can be compared to the amount of exogenic species detected in the upper atmosphere of Jupiter. We extracted light curves of the flashes and calculated the masses and sizes of the impacting objects. An examination of the number of amateur observations of Jupiter as a function of time allows us to interpret the statistics of these detections. The cumulative flux of small objects (5-20 m or larger) that impact Jupiter is predicted to be low (10-65 impacts per year), and only a fraction of them are potentially observable from Earth (4-25 per year in a perfect survey). More impacts will be found in the next years, with Jupiter opposition displaced toward summer in the northern hemisphere. Objects of this size contribute negligibly to the exogenous species and dust in the stratosphere of Jupiter when compared with the continuous flux from interplanetary dust punctuated by giant impacts. Flashes of a high enough could produce an observable debris field on the planet. We estimate that a continuous search for these impacts might find these events once every 0.4 to 2.6 years.

At the cloud top level of Venus (65-70 km altitude) the atmosphere rotates 60 times faster than the underlying surface, a phenomenon known as superrotation. Whereas on Venus's dayside the cloud top motions are well determined and Venus general circulation models predict a mean zonal flow at the upper clouds similar on both day and nightside, the nightside circulation remains poorly studied except for the polar region. Here we report global measurements of the nightside circulation at the upper cloud level. We tracked individual features in thermal emission images at 3.8 and 5.0 $\mathrm{\mu m}$ obtained between 2006 and 2008 by the Visible and Infrared Thermal Imaging Spectrometer (VIRTIS-M) onboard Venus Express and in 2015 by ground-based measurements with the Medium-Resolution 0.8-5.5 Micron Spectrograph and Imager (SpeX) at the National Aeronautics and Space Administration Infrared Telescope Facility (NASA/IRTF). The zonal motions range from -110 to -60 m s$^{-1}$, consistent with those found for the dayside but with larger dispersion. Slow motions (-50 to -20 m s$^{-1}$) were also found and remain unexplained. In addition, abundant stationary wave patterns with zonal speeds from -10 to +10 m s$^{-1}$ dominate the night upper clouds and concentrate over the regions of higher surface elevation.

We provide measurements of the absolute reflectivity of Jupiter and Saturn along their central meridians in filters covering a wide range of visible and near-infrared wavelengths (from 0.38 to 1.7 $\mu$m) that are not often presented in the literature. We also give measurements of the geometric albedo of both planets and discuss the limb-darkening behavior and temporal variability of their reflectivity values for a period of four years (2012-2016). This work is based on observations with the PlanetCam-UPV/EHU instrument at the 1.23 m and 2.2 m telescopes in Calar Alto Observatory (Spain). The instrument simultaneously observes in two channels: visible (VIS; 0.38-1.0 $\mu$m) and short-wave infrared (SWIR; 1.0--1.7 $\mu$m). We obtained high-resolution observations via the lucky-imaging method. We show that our calibration is consistent with previous independent determinations of reflectivity values of these planets and, for future reference, provide new data extended in the wavelength range and in the time. Our results have an uncertainty in absolute calibration of 10--20\%. We show that under the hypothesis of constant geometric albedo, we are able to detect absolute reflectivity changes related to planetary temporal evolution of about 5-10\%.

Since 2013, observations of Neptune with small telescopes have resulted in several detections of long-lived bright atmospheric features that have also been observed by large telescopes such as Keck II or Hubble. The combination of both types of images allows the study of the long term evolution of major cloud systems in the planet. In 2013 and 2014 two bright features were present on the planet at southern mid latitudes. These may have merged in late 2014, possibly leading to the formation of a single bright feature observed during 2015 at the same latitude. This cloud system was first observed in January 2015 and nearly continuously from July to December 2015 in observations with telescopes in the 2 to 10 meter class and in images from amateur astronomers. These images show the bright spot as a compact feature at 40.1 deg South planetographic latitude well resolved from a nearby bright zonal band that extended from 42 deg South to 20 deg South. Tracking its motion from July to November 2015 suggests a longitudinal oscillation of 16 deg in amplitude with a 90 day period, typical of dark spots on Neptune and similar to the Great Red Spot oscillation in Jupiter. The limited time covered by high-resolution observations only covers one full oscillation and other interpretations of the changing motions could be possible. HST images in September 2015 show the presence of a dark spot at short wavelengths in the southern flank of the bright cloud observed throughout 2015.

The ice giants Uranus and Neptune are the least understood class of planets in our solar system but the most frequently observed type of exoplanets. Presumed to have a small rocky core, a deep interior comprising ~70% heavy elements surrounded by a more dilute outer envelope of H2 and He, Uranus and Neptune are fundamentally different from the better-explored gas giants Jupiter and Saturn. Because of the lack of dedicated exploration missions, our knowledge of the composition and atmospheric processes of these distant worlds is primarily derived from remote sensing from Earth-based observatories and space telescopes. As a result, Uranus's and Neptune's physical and atmospheric properties remain poorly constrained and their roles in the evolution of the Solar System not well understood. Exploration of an ice giant system is therefore a high-priority science objective as these systems (including the magnetosphere, satellites, rings, atmosphere, and interior) challenge our understanding of planetary formation and evolution. Here we describe the main scientific goals to be addressed by a future in situ exploration of an ice giant. An atmospheric entry probe targeting the 10-bar level, about 5 scale heights beneath the tropopause, would yield insight into two broad themes: i) the formation history of the ice giants and, in a broader extent, that of the Solar System, and ii) the processes at play in planetary atmospheres. The probe would descend under parachute to measure composition, structure, and dynamics, with data returned to Earth using a Carrier Relay Spacecraft as a relay station. In addition, possible mission concepts and partnerships are presented, and a strawman ice-giant probe payload is described. An ice-giant atmospheric probe could represent a significant ESA contribution to a future NASA ice-giant flagship mission.

The VESPA data access system focuses on applying Virtual Observatory (VO) standards and tools to Planetary Science. Building on a previous EC-funded Europlanet program, it has reached maturity during the first year of a new Europlanet 2020 program (started in 2015 for 4 years). The infrastructure has been upgraded to handle many fields of Solar System studies, with a focus both on users and data providers. This paper describes the broad lines of the current VESPA infrastructure as seen by a potential user, and provides examples of real use cases in several thematic areas. These use cases are also intended to identify hints for future developments and adaptations of VO tools to Planetary Science.

Since 2003 the Planetary Virtual Observatory and Laboratory (PVOL) has been storing and serving publicly through its web site a large database of amateur observations of the Giant Planets (Hueso et al., 2010a). These images are used for scientific research of the atmospheric dynamics and cloud structure on these planets and constitute a powerful resource to address time changing phenomena in their atmospheres. Advances over the last decade in observation techniques, and a wider conscience by professional astronomers of the quality of amateur observations, have resulted in the necessity to upgrade this database. We here present major advances in the PVOL database that has evolved into a full virtual planetary observatory encompassing also observations of Mercury, Venus, Mars, the Moon and the Galilean satellites. Besides the new objects, the images can be tagged and the database allows simple and complex searches over the data. The new web service: PVOL2 is available online in http://pvol2.ehu.eus/ , contains a fully functional search engine and constitutes one of the many services included in VESPA (Virtual Europan Solar and Planetary Access). Data from PVOL2 can be served from the VESPA portal using the EPN-TAP protocol. PVOL2 also provides long-term storage to amateur observations containing about 30,000 amateur observations starting in the year 2000. Current and past observations from the amateur community provide a global view of the Solar System planets over the years with several possibilities for scientific analysis and amateur astronomers involvement in planetary science.

The Europlanet 2020 Research Infrastructure will include new Planetary Space Weather Services (PSWS) that will extend the concepts of space weather and space situational awareness to other planets in our Solar System and in particular to spacecraft that voyage through it. PSWS will make five entirely new toolkits accessible to the research community and to industrial partners planning for space missions: a general planetary space weather toolkit, as well as three toolkits dedicated to the following key planetary environments: Mars, comets, and outer planets. This will give the European planetary science community new methods, interfaces, functionalities and/or plugins dedicated to planetary space weather in the tools and models available within the partner institutes. It will also create a novel event-diary toolkit aiming at predicting and detecting planetary events like meteor showers and impacts. A variety of tools are available for tracing propagation of planetary and/or solar events through the Solar System and modelling the response of the planetary environment (surfaces, atmospheres, ionospheres, and magnetospheres) to those events. But these tools were not originally designed for planetary event prediction and space weather applications. PSWS will provide the additional research and tailoring required to apply them for these purposes. PSWS will be to review, test, improve and adapt methods and tools available within the partner institutes in order to make prototype planetary event and space weather services operational in Europe at the end of the programme. To achieve its objectives PSWS will use a few tools and standards developed for the Astronomy Virtual Observatory (VO). This paper gives an overview of the project together with a few illustrations of prototype services based on VO standards and protocols.

We report Venus image observations around the two maximum elongations of the planet at June and October 2015. From these images we describe the global atmospheric dynamics and cloud morphology in the planet before the arrival of JAXA Akatsuki mission on December the 7th. The majority of the images were acquired at ultraviolet wavelengths (380-410 nm) using small telescopes. The Venus dayside was also observed with narrow band filters at other wavelengths (890 nm, 725-950 nm, 1.435 {\mu}m CO2 band) using the instrument PlanetCam-UPV/EHU at the 2.2m telescope in Calar Alto Observatory. In all cases, the lucky imaging methodology was used to improve the spatial resolution of the images over the atmospheric seeing. During the April-June period, the morphology of the upper cloud showed an irregular and chaotic texture with a well developed equatorial dark belt (afternoon hemisphere), whereas during October-December the dynamical regime was dominated by planetary-scale waves (Yhorizontal, C-reversed and {\psi}-horizontal features) formed by long streaks, and banding suggesting more stable conditions. Measurements of the zonal wind velocity with cloud tracking in the latitude range from 50$^{\circ}$N to 50$^{\circ}$S shows agreement with retrievals from previous works.

The Hera Saturn entry probe mission is proposed as an M--class mission led by ESA with a contribution from NASA. It consists of one atmospheric probe to be sent into the atmosphere of Saturn, and a Carrier-Relay spacecraft. In this concept, the Hera probe is composed of ESA and NASA elements, and the Carrier-Relay Spacecraft is delivered by ESA. The probe is powered by batteries, and the Carrier-Relay Spacecraft is powered by solar panels and batteries. We anticipate two major subsystems to be supplied by the United States, either by direct procurement by ESA or by contribution from NASA: the solar electric power system (including solar arrays and the power management and distribution system), and the probe entry system (including the thermal protection shield and aeroshell). Hera is designed to perform in situ measurements of the chemical and isotopic compositions as well as the dynamics of Saturn's atmosphere using a single probe, with the goal of improving our understanding of the origin, formation, and evolution of Saturn, the giant planets and their satellite systems, with extrapolation to extrasolar planets. Hera's aim is to probe well into the cloud-forming region of the troposphere, below the region accessible to remote sensing, to the locations where certain cosmogenically abundant species are expected to be well mixed. By leading to an improved understanding of the processes by which giant planets formed, including the composition and properties of the local solar nebula at the time and location of giant planet formation, Hera will extend the legacy of the Galileo and Cassini missions by further addressing the creation, formation, and chemical, dynamical, and thermal evolution of the giant planets, the entire solar system including Earth and the other terrestrial planets, and formation of other planetary systems.

This white paper examines the benefit of the upcoming James Webb Space Telescope for studies of the Solar System's four giant planets: Jupiter, Saturn, Uranus, and Neptune. JWST's superior sensitivity, combined with high spatial and spectral resolution, will enable near- and mid-infrared imaging and spectroscopy of these objects with unprecedented quality. In this paper we discuss some of the myriad scientific investigations possible with JWST regarding the giant planets. This discussion is preceded by the specifics of JWST instrumentation most relevant to giant planet observations. We conclude with identification of desired pre-launch testing and operational aspects of JWST that would greatly benefit future studies of the giant planets.

The discovery of almost 2000 exoplanets has revealed an unexpectedly diverse planet population. Observations to date have shown that our Solar System is certainly not representative of the general population of planets in our Milky Way. The key science questions that urgently need addressing are therefore: What are exoplanets made of? Why are planets as they are? What causes the exceptional diversity observed as compared to the Solar System? EChO (Exoplanet Characterisation Observatory) has been designed as a dedicated survey mission for transit and eclipse spectroscopy capable of observing a large and diverse planet sample within its four-year mission lifetime. EChO can target the atmospheres of super-Earths, Neptune-like, and Jupiter-like planets, in the very hot to temperate zones (planet temperatures of 300K-3000K) of F to M-type host stars. Over the next ten years, several new ground- and space-based transit surveys will come on-line (e.g. NGTS, CHEOPS, TESS, PLATO), which will specifically focus on finding bright, nearby systems. The current rapid rate of discovery would allow the target list to be further optimised in the years prior to EChO's launch and enable the atmospheric characterisation of hundreds of planets. Placing the satellite at L2 provides a cold and stable thermal environment, as well as a large field of regard to allow efficient time-critical observation of targets randomly distributed over the sky. A 1m class telescope is sufficiently large to achieve the necessary spectro-photometric precision. The spectral coverage (0.5-11 micron, goal 16 micron) and SNR to be achieved by EChO, thanks to its high stability and dedicated design, would enable a very accurate measurement of the atmospheric composition and structure of hundreds of exoplanets.

Remote sensing observations meet some limitations when used to study the bulk atmospheric composition of the giant planets of our solar system. A remarkable example of the superiority of in situ probe measurements is illustrated by the exploration of Jupiter, where key measurements such as the determination of the noble gases abundances and the precise measurement of the helium mixing ratio have only been made available through in situ measurements by the Galileo probe. This paper describes the main scientific goals to be addressed by the future in situ exploration of Saturn placing the Galileo probe exploration of Jupiter in a broader context and before the future probe exploration of the more remote ice giants. In situ exploration of Saturn's atmosphere addresses two broad themes that are discussed throughout this paper: first, the formation history of our solar system and second, the processes at play in planetary atmospheres. In this context, we detail the reasons why measurements of Saturn's bulk elemental and isotopic composition would place important constraints on the volatile reservoirs in the protosolar nebula. We also show that the in situ measurement of CO (or any other disequilibrium species that is depleted by reaction with water) in Saturn's upper troposphere would constrain its bulk O/H ratio. We highlight the key measurements required to distinguish competing theories to shed light on giant planet formation as a common process in planetary systems with potential applications to most extrasolar systems. In situ measurements of Saturn's stratospheric and tropospheric dynamics, chemistry and cloud-forming processes will provide access to phenomena unreachable to remote sensing studies. Different mission architectures are envisaged, which would benefit from strong international collaborations.

Amateur contributions to professional publications have increased exponentially over the last decades in the field of Planetary Astronomy. Here we review the different domains of the field in which collaborations between professional and amateur astronomers are effective and regularly lead to scientific publications. We discuss the instruments, detectors, softwares and methodologies typically used by amateur astronomers to collect the scientific data in the different domains of interest. Amateur contributions to the monitoring of planets and interplanetary matter, characterization of asteroids and comets, as well as the determination of the physical properties of Kuiper Belt Objects and exoplanets are discussed.

We investigated the Venus O2 visible nightglow with imagery from the Venus Monitoring Camera on Venus Express. Drawing from data collected between April 2007 and January 2011, we study the global distribution of this emission, discovered in the late 70s by the Venera 9 and 10 missions. The inferred limb-viewing intensities are on the order of 150 kiloRayleighs at the lower latitudes and seem to drop somewhat towards the poles. The emission is generally stable, although there are episodes when the intensities rise up to 500 kR. We compare a set of Venus Monitoring Camera observations with coincident measurements of the O2 nightglow at 1.27 {\mu}m made with the Visible and Infrared Thermal Imaging Spectrometer, also on Venus Express. From the evidence gathered in this and past works, we suggest a direct correlation between the instantaneous emissions from the two O2 nightglow systems. Possible implications regarding the uncertain origin of the atomic oxygen green line at 557.7 nm are noted.

A dedicated mission to investigate exoplanetary atmospheres represents a major milestone in our quest to understand our place in the universe by placing our Solar System in context and by addressing the suitability of planets for the presence of life. EChO -the Exoplanet Characterisation Observatory- is a mission concept specifically geared for this purpose. EChO will provide simultaneous, multi-wavelength spectroscopic observations on a stable platform that will allow very long exposures. EChO will build on observations by Hubble, Spitzer and groundbased telescopes, which discovered the first molecules and atoms in exoplanetary atmospheres. EChO will simultaneously observe a broad enough spectral region -from the visible to the mid-IR- to constrain from one single spectrum the temperature structure of the atmosphere and the abundances of the major molecular species. The spectral range and resolution are tailored to separate bands belonging to up to 30 molecules to retrieve the composition and temperature structure of planetary atmospheres. The target list for EChO includes planets ranging from Jupiter-sized with equilibrium temperatures Teq up to 2000 K, to those of a few Earth masses, with Teq ~300 K. We have baselined a dispersive spectrograph design covering continuously the 0.4-16 micron spectral range in 6 channels (1 in the VIS, 5 in the IR), which allows the spectral resolution to be adapted from several tens to several hundreds, depending on the target brightness. The instrument will be mounted behind a 1.5 m class telescope, passively cooled to 50 K, with the instrument structure and optics passively cooled to ~45 K. EChO will be placed in a grand halo orbit around L2. We have also undertaken a first-order cost and development plan analysis and find that EChO is easily compatible with the ESA M-class mission framework.

Cosmic collisions on planets cause detectable optical flashes that range from terrestrial shooting stars to bright fireballs. On June 3, 2010 a bolide in Jupiter's atmosphere was simultaneously observed from the Earth by two amateur astronomers observing Jupiter in red and blue wavelengths. The bolide appeared as a flash of 2 s duration in video recording data of the planet. The analysis of the light curve of the observations results in an estimated energy of the impact of 0.9-4.0x10^{15} J which corresponds to a colliding body of 8-13 m diameter assuming a mean density of 2 g cm^{-3}. Images acquired a few days later by the Hubble Space Telescope and other large ground-based facilities did not show any signature of aerosol debris, temperature or chemical composition anomaly, confirming that the body was small and destroyed in Jupiter's upper atmosphere. Several collisions of this size may happen on Jupiter on a yearly basis. A systematic study of the impact rate and size of these bolides can enable an empirical determination of the flux of meteoroids in Jupiter with implications for the populations of small bodies in the outer Solar System and may allow a better quantification of the threat of impacting bodies to Earth. The serendipitous recording of this optical flash opens a new window in the observation of Jupiter with small telescopes.

On 2009 July 19, we observed a single, large impact on Jupiter at a planetocentric latitude of 55^{\circ}S. This and the Shoemaker-Levy 9 (SL9) impacts on Jupiter in 1994 are the only planetary-scale impacts ever observed. The 2009 impact had an entry trajectory opposite and with a lower incidence angle than that of SL9. Comparison of the initial aerosol cloud debris properties, spanning 4,800 km east-west and 2,500 km north-south, with those produced by the SL9 fragments, and dynamical calculations of pre-impact orbit, indicate that the impactor was most probably an icy body with a size of 0.5-1 km. The collision rate of events of this magnitude may be five to ten times more frequent than previously thought. The search for unpredicted impacts, such as the current one, could be best performed in 890-nm and K (2.03-2.36 {\mu}m) filters in strong gaseous absorption, where the high-altitude aerosols are more reflective than Jupiter's primary cloud.

It has been proposed that the enrichment in noble gases found by Galileo in Jupiter's atmosphere can be explained by their delivery inside cold planetesimals. We propose instead that this is a sign that the planet formed in a chemically evolved disk and that noble gases were acquired mostly in gaseous form during the planet's envelope capture phase. We show that the combined settling of grains to the disk midplane in the cold outer layers, the condensation of noble gases onto these grains at temperatures below 20-30K, and the evaporation from high disk altitudes effectively lead to a progressive, moderate enrichment of the disk. The fact that noble gases are vaporized from the grains in the hot inner disk regions (e.g. Jupiter formation region) is not a concern because a negative temperature gradient prevents convection from carrying the species into the evaporating region. We show that the ~2 times solar enrichment of Ar, Kr, Xe in Jupiter is hence naturally explained by a continuous growth of the planet governed by viscous diffusion in the protosolar disk in conjunction with an evaporation of the disk and its progressive enrichment on a million years timescale.

We present a one-dimensional model of the formation and viscous evolution of protoplanetary disks. The formation of the early disk is modeled as the result of the gravitational collapse of an isothermal molecular cloud. The disk's viscous evolution is integrated according to two parameterizations of turbulence: The classical $\alpha$ representation and a $\beta$ parameterization, representative of non-linear turbulence driven by the keplerian shear. We apply the model to DM Tau and GM Aur, two classical T-Tauri stars with relatively well-characterized disks, retrieving the evolution of their surface density with time. We perform a systematic Monte-Carlo exploration of the parameter space (i.e. values of the $\alpha$-$\beta$ parameters, and of the temperature and rotation rate in the molecular cloud) to find the values that are compatible with the observed disk surface density distribution, star and disk mass, age and present accretion rate. We find that the observations for DM Tau require $0.001<\alpha<0.1$ or $2\times 10^{-5}<\beta<5\times 10^{-4}$. For GM Aur, we find that the turbulent viscosity is such that $4\times 10^{-4}<\alpha<0.01$ or $2\times 10^{-6}<\beta<8\times 10^{-5}$. These relatively large values show that an efficient turbulent diffusion mechanism is present at distances larger than $\sim 10 $AU. This is to be compared to studies of the variations of accretion rates of T-Tauri stars versus age that mostly probe the inner disks, but also yield values of $\alpha\sim 0.01$. We show that the mechanism responsible for turbulent diffusion at large orbital distances most probably cannot be convection because of its suppression at low optical depths.