We report six stellar occultations by (Saturn IX) Phoebe, an irregular satellite of Saturn, obtained between mid-2017 and mid-2019. The 2017 July 06 event is the first stellar occultation by an irregular satellite ever observed. The occultation chords were compared to a 3D shape model of the satellite obtained from Cassini observations. The rotation period available in the literature led to a sub-observer point at the moment of the observed occultations where the chords could not fit the 3D model. A procedure was developed to identify the correct sub-observer longitude. It allowed us to obtain the rotation period with improved precision over currently known value from literature. We show that the difference between the observed and the predicted sub-observer longitude suggests two possible solutions for the rotation period. By comparing these values with recently observed rotational light curves and single-chord stellar occultations, we can identify the best solution for Phoebe's rotational period as $9.27365 \pm 0.00002$ h. From the stellar occultations, we also obtained 6 geocentric astrometric positions in the ICRS as realised by the Gaia-DR2 with uncertainties at the 1-mas level.

Progress in astrometry and orbital modelling of planetary moons in the last decade enabled better determinations of their orbits. These studies need accurate positions spread over extended periods. We present the results of the 2014-2015 Brazilian campaign for 40 mutual events from 47 observed light curves by the Galilean satellites plus one eclipse of Amalthea by Ganymede. We also reanalysed and updated results for 25 mutual events observed in the 2009 campaign. All telescopes were equipped with narrow-band filters centred at 889 nm with a width of 15 nm to eliminate the scattered light from Jupiter. The albedos' ratio was determined using images before and after each event. We simulated images of moons, umbra, and penumbra in the sky plane, and integrated their fluxes to compute albedos, simulate light curves and fit them to the observed ones using a chi-square fitting procedure. For that, we used the complete version of the Oren-Nayer reflectance model. The relative satellite positions mean uncertainty was 11.2 mas ($\sim$35 km) and 10.1 mas ($\sim$31 km) for the 2014-2015 and 2009 campaigns respectively. The simulated and observed \textsc{ascii} light curve files are freely available in electronic form at the \textit{Natural Satellites DataBase} (NSDB). The 40/25 mutual events from our 2014-2015/2009 campaigns represent a significant contribution of 17%/15% in comparison with the PHEMU campaigns lead by the IMCCE. Besides that, our result for the eclipse of Amalthea is only the 4$^{th}$ such measurement ever published after the three ones observed by the 2014-2015 international PHEMU campaign. Our results are suitable for new orbital/ephemeris determinations for the Galilean moons and Amalthea.

Context. Bright stellar positions are now known with an uncertainty below 1 mas thanks to Gaia DR2. Between 2019-2020, the Galactic plane will be the background of Jupiter. The dense stellar background will lead to an increase in the number of occultations, while the Gaia DR2 catalogue will reduce the prediction uncertainties for the shadow path. Aims. We observed a stellar occultation by the Galilean moon Europa (J2) and propose a campaign for observing stellar occultations for all Galilean moons. Methods. During a predicted period of time, we measured the light flux of the occulted star and the object to determine the time when the flux dropped with respect to one or more reference stars, and the time that it rose again for each observational station. The chords obtained from these observations allowed us to determine apparent sizes, oblatness, and positions with kilometre accuracy. Results. We present results obtained from the first stellar occultation by the Galilean moon Europa observed on 2017 March 31. The apparent fitted ellipse presents an equivalent radius of 1561.2 $\pm$ 3.6 km and oblatenesses 0.0010 $\pm$ 0.0028. A very precise Europa position was determined with an uncertainty of 0.8 mas. We also present prospects for a campaign to observe the future events that will occur between 2019 and 2021 for all Galilean moons. Conclusions. Stellar occultation is a suitable technique for obtaining physical parameters and highly accurate positions of bright satellites close to their primary. A number of successful events can render the 3D shapes of the Galilean moons with high accuracy. We encourage the observational community (amateurs included) to observe the future predicted events.

Context. Pluto's tenuous nitrogen (N2) atmosphere undergoes strong seasonal effects due to high obliquity and orbital eccentricity, and has been recently (July 2015) observed by the New Horizons spacecraft. Goals are (i) construct a well calibrated record of the seasonal evolution of surface pressure on Pluto and (ii) constrain the structure of the lower atmosphere using a central flash observed in 2015. Method: eleven stellar occultations by Pluto observed between 2002 and 2016 are used to retrieve atmospheric profiles (density, pressure, temperature) between $\sim$5 km and $\sim$380 km altitude levels (i.e. pressures from about 10 microbar to 10 nanobar). Results: (i) Pressure has suffered a monotonic increase from 1988 to 2016, that is compared to a seasonal volatile transport model, from which tight constraints on a combination of albedo and emissivity of N2 ice are derived; (ii) A central flash observed on 2015 June 29 is consistent with New Horizons REX profiles, provided that (a) large diurnal temperature variations (not expected by current models) occur over Sputnik Planitia and/or (b) hazes with tangential optical depth of about 0.3 are present at 4-7 km altitude levels and/or (c) the nominal REX density values are overestimated by an implausibly large factor of about 20% and/or (d) higher terrains block part of the flash in the Charon facing hemisphere.

The technique of mutual approximations accurately gives the central instant at the maximum apparent approximation of two moving natural satellites in the sky plane. This can be used in ephemeris fitting to infer the relative positions between satellites with high precision. Only the mutual phenomena -- occultations and eclipses -- may achieve better results. However, mutual phenomena only occur every six years in the case of Jupiter. Mutual approximations do not have this restriction and can be observed at any time along the year as long as the satellites are visible. In this work, we present 104 central instants determined from the observations of 66 mutual approximations between the Galilean moons carried out at different sites in Brazil and France during the period 2016--2018. For 28 events we have at least two independent observations. All telescopes were equipped with a narrow-band filter centred at 889 nm with a width of 15 nm to eliminate the scattered light from Jupiter. The telescope apertures ranged between 25--120 cm. For comparison, the precision of the positions obtained with classical CCD astrometry is about 100 mas, for mutual phenomena it can achieve 10 mas or less and the average internal precision obtained with mutual approximations was 11.3 mas. This new kind of simple, yet accurate observations can significantly improve the orbits and ephemeris of Galilean satellites and thus be very useful for the planning of future space missions aiming at the Jovian system.

Two narrow and dense rings (called C1R and C2R) were discovered around the Centaur object (10199) Chariklo during a stellar occultation observed on 2013 June 3. Following this discovery, we planned observations of several occultations by Chariklo's system in order to better characterize the physical properties of the ring and main body. Here, we use 12 successful occulations by Chariklo observed between 2014 and 2016. They provide ring profiles (physical width, opacity, edge structure) and constraints on the radii and pole position. Our new observations are currently consistent with the circular ring solution and pole position, to within the $\pm 3.3$ km formal uncertainty for the ring radii derived by Braga-Ribas et al. The six resolved C1R profiles reveal significant width variations from $\sim 5$ to 7.5 km. The width of the fainter ring C2R is less constrained, and may vary between 0.1 and 1 km. The inner and outer edges of C1R are consistent with infinitely sharp boundaries, with typical upper limits of one kilometer for the transition zone between the ring and empty space. No constraint on the sharpness of C2R's edges is available. A 1$\sigma$ upper limit of $\sim 20$ m is derived for the equivalent width of narrow (physical width <4 km) rings up to distances of 12,000 km, counted in the ring plane.

We use data from five stellar occultations observed between 2013 and 2016 to constrain Chariklo's size and shape, and the ring reflectivity. We consider four possible models for Chariklo (sphere, Maclaurin spheroid, tri-axial ellipsoid and Jacobi ellipsoid) and we use a Bayesian approach to estimate the corresponding parameters. The spherical model has a radius $R=129\pm3$ km. The Maclaurin model has equatorial and polar radii $a=b=143^{+3}_{-6}$ km and $c=96^{+14}_{-4}$ km, respectively, with density $970^{+300}_{-180}$ kg m$^{-3}$. The ellipsoidal model has semiaxes $a=148^{+6}_{-4}$ km, $b=132^{+6}_{-5}$ km and $c=102^{+10}_{-8}$ km. Finally, the Jacobi model has semiaxes $a$=157$\pm$4 km, $b$=139$\pm$ 4 km and $c$=86$\pm$1 km, and density $796^{+2}_{-4}$ kg m$^{-3}$ . Depending on the model, we obtain topographic features of 6-11 km, typical of Saturn icy satellites with similar size and density. We constrain Chariklo's geometric albedo between 3.1\% (sphere) and 4.9\% (ellipsoid), while the ring $I/F$ reflectivity is less constrained between 0.6\% (Jacobi) and 8.9\% (sphere). The ellipsoid model explains both the optical light curve and the long-term photometry variation of the system, giving a plausible value for the geometric albedo of the ring particles of $10-15\%$. The derived Chariklo's mass of 6-8$\times10^{18}$ kg places the rings close to the 3:1 resonance between the ring mean motion and Chariklo's rotation period.

Gomes-J\'unior et al. (2015) published 3613 positions for the 8 largest irregular satellites of Jupiter and 1787 positions for the largest irregular satellite of Saturn, Phoebe. These observations were made between 1995 and 2014 and have an estimated error of about 60 to 80 mas. Based on this set of positions, we derived new orbits for the eight largest irregular satellites of Jupiter: Himalia, Elara, Pasiphae, Carme, Lysithea, Sinope, Ananke and Leda. For Phoebe we updated the ephemeris from Desmars et al. (2013) using 75% more positions than the previous one. Due to their orbital characteristics, it is common belief that the irregular satellites were captured by the giant planets in the early Solar System, but there is no consensus for a single model explaining where they were formed. Size, shape, albedo and composition would help to trace back their true origin, but these physical parameters are yet poorly known for irregular satellites. The observation of stellar occultations would allow for the determination of such parameters. Indeed Jupiter will cross the galactic plane in 2019-2020 and Saturn in 2018, improving a lot the chances of observing such events in the near future. Using the derived ephemerides and the UCAC4 catalogue we managed to identify 5442 candidate stellar occultations between January 2016 and December 2020 for the 9 satellites studied here. We discussed how the successful observation of a stellar occultation by these objects is possible and present some potential occultations.

Typically we can deliver astrometric positions of natural satellites with errors in the 50-150 mas range. Apparent distances from mutual phenomena, have much smaller errors, less than 10 mas. However, this method can only be applied during the equinox of the planets. We developed a method that can provide accurate astrometric data for natural satellites -- the mutual approximations. The method can be applied when any two satellites pass close by each other in the apparent sky plane. The fundamental parameter is the central instant $t_0$ of the passage when the distances reach a minimum. We applied the method for the Galilean moons. All observations were made with a 0.6 m telescope with a narrow-band filter centred at 889 nm with width of 15 nm which attenuated Jupiter's scattered light. We obtained central instants for 14 mutual approximations observed in 2014-2015. We determined $t_0$ with an average precision of 3.42 mas (10.43 km). For comparison, we also applied the method for 5 occultations in the 2009 mutual phenomena campaign and for 22 occultations in the 2014-2015 campaign. The comparisons of $t_0$ determined by our method with the results from mutual phenomena show an agreement by less than 1-sigma error in $t_0$, typically less than 10 mas. This new method is particularly suitable for observations by small telescopes.

The prediction of stellar occultations by Transneptunian objects and Centaurs is a difficult challenge that requires accuracy both in the occulted star position as for the object ephemeris. Until now, the most used method of prediction involving tens of TNOs/Centaurs was to consider a constant offset for the right ascension and for the declination with respect to a reference ephemeris. This offset is determined as the difference between the most recent observations of the TNO and the reference ephemeris. This method can be successfully applied when the offset remains constant with time. This paper presents an alternative method of prediction based on a new accurate orbit determination procedure, which uses all the available positions of the TNO from the Minor Planet Center database plus sets of new astrometric positions from unpublished observations. The orbit determination is performed through a numerical integration procedure (NIMA), in which we develop a specific weighting scheme. The NIMA method was applied for 51 selected TNOs/Centaurs. For this purpose, we have performed about 2900 new observations during 2007-2014. Using NIMA, we succeed in predicting the stellar occultations of 10 TNOs and 3 Centaurs between 2013 and 2015. By comparing the NIMA and JPL ephemerides, we highlighted the variation of the offset between them with time. Giving examples, we show that the constant offset method could not accurately predict 6 out of the 13 observed positive occultations successfully predicted by NIMA. The results indicate that NIMA is capable of efficiently refine the orbits of these bodies. Finally, we show that the astrometric positions given by positive occultations can help to further refine the orbit of the TNO and consequently the future predictions. We also provide the unpublished observations of the 51 selected TNOs and their ephemeris in a usable format by the SPICE library.

We analyze two multi-chord stellar occultations by Pluto observed on July 18th, 2012 and May 4th, 2013, and monitored respectively from five and six sites. They provide a total of fifteen light-curves, twelve of them being used for a simultaneous fit that uses a unique temperature profile, assuming a clear (no-haze) and pure N_2 atmosphere, but allowing for a possible pressure variation between the two dates. We find a solution that fits satisfactorily (i.e. within the noise level) all the twelve light-curves, providing atmospheric constraints between ~1,190 km (pressure ~ 11 \mubar) and ~ 1,450 km (pressure ~0.1 \mubar) from Pluto's center. Our main results are: (1) the best-fitting temperature profile shows a stratosphere with strong positive gradient between 1,190 km (at 36 K, 11 \mubar) and r = 1,215 km (6.0 \mubar), where a temperature maximum of 110 K is reached; above it is a mesosphere with negative thermal gradient of -0.2 K/km up to ~ 1,390 km (0.25 \mubar), where, the mesosphere connects itself to a more isothermal upper branch around 81 K; (2) the pressure shows a small (6 %) but significant increase (6-\sigma level) between the two dates; (3) without troposphere, Pluto's radius is found to be R_P = 1,190 +/- 5km. Allowing for a troposphere, R_P is constrained to lie between 1,168 and 1,195 km; (4) the currently measured CO abundance is too small to explain the mesospheric negative thermal gradient. Cooling by HCN is possible, but only if this species is largely saturated; Alternative explanations like zonal winds or vertical compositional variations of the atmosphere are unable to explain the observed mesospheric trend.

We determine accurate positions of the main satellites of Uranus: Miranda, Ariel, Umbriel, Titania, and Oberon. Positions of Uranus, as derived from those of these satellites, are also determined. The observational period spans from 1992 to 2011. All runs were made at the Pico dos Dias Observatory, Brazil. We used the software called Platform for Reduction of Astronomical Images Automatically (PRAIA) to minimise (digital coronography) the influence of the scattered light of Uranus on the astrometric measurements and to determine accurate positions of the main satellites. The positions of Uranus were then indirectly determined by computing the mean differences between the observed and ephemeris positions of these satellites. A series of numerical filters was applied to filter out spurious data. These filters are mostly based on the comparison between the positions of Oberon with those of the other satellites and on the offsets as given by the differences between the observed and ephemeris positions of all satellites. We have, for the overall offsets of the five satellites, -29 (+/-63) mas in right ascension and -27 (+/-46) mas in declination. For the overall difference between the offsets of Oberon and those of the other satellites, we have +3 (+/-30) mas in right ascension and -2 (+/-28) mas in declination. Ephemeris positions for the satellites were determined from DE432+ura111. Comparisons using other modern ephemerides for the solar system -INPOP13c- and for the motion of the satellites -NOE-7-2013- were also made. They confirm that the largest contribution to the offsets we find comes from the motion of the barycenter of the Uranus system around the barycenter of the solar system, as given by the planetary ephemerides. Catalogues with the observed positions are provided.

The irregular satellites of the giant planets are believed to have been captured during the evolution of the solar system. Knowing their physical parameters, such as size, density, and albedo is important for constraining where they came from and how they were captured. The best way to obtain these parameters are observations in situ by spacecrafts or from stellar occultations by the objects. Both techniques demand that the orbits are well known. We aimed to obtain good astrometric positions of irregular satellites to improve their orbits and ephemeris. We identified and reduced observations of several irregular satellites from three databases containing more than 8000 images obtained between 1992 and 2014 at three sites (Observat\'orio do Pico dos Dias, Observatoire de Haute-Provence, and European Southern Observatory - La Silla). We used the software PRAIA (Platform for Reduction of Astronomical Images Automatically) to make the astrometric reduction of the CCD frames. The UCAC4 catalog represented the International Celestial Reference System in the reductions. Identification of the satellites in the frames was done through their ephemerides as determined from the SPICE/NAIF kernels. Some procedures were followed to overcome missing or incomplete information (coordinates, date), mostly for the older images. We managed to obtain more than 6000 positions for 18 irregular satellites: 12 of Jupiter, 4 of Saturn, 1 of Uranus (Sycorax), and 1 of Neptune (Nereid). For some satellites the number of obtained positions is more than 50\% of what was used in earlier orbital numerical integrations. Comparison of our positions with recent JPL ephemeris suggests there are systematic errors in the orbits for some of the irregular satellites. The most evident case was an error in the inclination of Carme.

We report the results of two multi-chord stellar occultations by the dwarf planet (1) Ceres that were observed from Brazil on 2010 August 17, and from the USA on 2013 October 25. Four positive detections were obtained for the 2010 occultation, and nine for the 2013 occultation. Elliptical models were adjusted to the observed chords to obtain Ceres' size and shape. Two limb fitting solutions were studied for each event. The first one is a nominal solution with an indeterminate polar aspect angle. The second one was constrained by the pole coordinates as given by Drummond et al. Assuming a Maclaurin spheroid, we determine an equatorial diameter of 972 $\pm$ 6 km and an apparent oblateness of 0.08 $\pm$ 0.03 as our best solution. These results are compared to all available size and shape determinations for Ceres made so far, and shall be confirmed by the NASA's Dawn space mission.