Context . Binary asteroids provide key access to fundamental parameters of Solar System remnants and planetary formations. However, the current knowledge of binary asteroids remains strongly biased by observational limitations, and main belt binary systems are still poorly characterised since current techniques preferentially detect either widely separated binaries through direct imaging or close and bright systems via photometry and radar for near-Earth asteroids. In this context, the high-precision astrometry of the Gaia mission has revealed a new population of candidate binaries exhibiting dynamical signatures consistent with unresolved companions. Stellar occultations have therefore emerged as one of the most effective methods to confirm the binary nature of a candidate and improve the current census of intermediate-size systems. Aims . This work is part of the GaiaMoons program, and our aim with it was to characterise a sample of 357 potential binary asteroid targets and confirm or refute their binary nature. The properties of these candidates were derived from the high-precision photometric and astrometric observations provided by the Gaia satellite. Methods . We adopted stellar occultation as the observational method to study these targets. Between October 2023 and February 2026, we successfully carried out 165 observations for 101 targets. We subsequently analysed these events in the context of the available literature and previously reported observations. Results . Out of the 165 observations, 76 led at least to one positive observation. Among these, 33 had at least two positives for 24 objects that have undergone unprecedented occultation observation campaigns, with four objects showing indications of binary or contact binary features, namely (1127) Mimi, (35420) 1998 AG 6 , (206) Hersilia, and (36882) 2000 SW 155 . For the vast majority of these objects, the resulting dataset from all reduced observations provides unique physical and astrometric constraints, as they had never been observed through stellar occultations before. In addition, 89 observations with only negatives allowed the near environment of the targets to be probed. Conclusions . GaiaMoons illustrates how stellar occultation campaigns associated with Gaia observations generate a self-improving cycle to find new binary, thereby probing size and shape to constrain future observations. By standardising this approach, we deliver critical data in unexplored parameter spaces, resolving long-standing observational ambiguities.

Context. The tenuous nitrogen (N2) atmosphere on Pluto undergoes strong seasonal effects due to high obliquity and orbital eccentricity, and has recently (July 2015) been observed by the New Horizons spacecraft. Aims. The main goals of this study are (i) to construct a well calibrated record of the seasonal evolution of surface pressure on Pluto and (ii) to constrain the structure of the lower atmosphere using a central flash observed in 2015. Methods. Eleven stellar occultations by Pluto observed between 2002 and 2016 are used to retrieve atmospheric profiles (density, pressure, temperature) between altitude levels of ~5 and ~380 km (i.e. pressures from ~ 10 μbar to 10 nbar). 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 ~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 ~20%; and/or (d) higher terrains block part of the flash in the Charon facing hemisphere.

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, triaxial ellipsoid, and Jacobi ellipsoid), and we use a Bayesian approach to estimate the corresponding parameters. The spherical model has a radius R=129±3 km. The Maclaurin model has equatorial and polar radii = = + a b 143 km 6 3 and = + c 96 km 4 14 , respectively, with density + 970 kg m 180 300 3. The ellipsoidal model has semiaxes = + a 148 km 4 6 , = + b 132 km 5 6 , and = + c 102 km 8 10 . Finally, the Jacobi model has semiaxes a=157±4km, b=139±4km, and c=86±1km, and density + 796 kg m 4 2 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 mass of Chariklo of 6–8×10 kg places the rings close to 3:1 resonance between the ring mean motion and Chariklo’s rotation period.

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 +/- 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 similar to 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 similar to 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.

ABSTRACT We present results from a multi-chord Pluto stellar occultation observed on 2015 June 29 from New Zealand and Australia. This occurred only two weeks before the NASA New Horizons flyby of the Pluto system and serves as a useful comparison between ground-based and space results. We find that Pluto's atmosphere is still expanding, with a significant pressure increase of 5 ± 2% since 2013 and a factor of almost three since 1988. This trend rules out, as of today, an atmospheric collapse associated with Pluto's recession from the Sun. A central flash, a rare occurrence, was observed from several sites in New Zealand. The flash shape and amplitude are compatible with a spherical and transparent atmospheric layer of roughly 3 km in thickness whose base lies at about 4 km above Pluto's surface, and where an average thermal gradient of about 5 K km−1 prevails. We discuss the possibility that small departures between the observed and modeled flash are caused by local topographic features (mountains) along Pluto's limb that block the stellar light. Finally, using two possible temperature profiles, and extrapolating our pressure profile from our deepest accessible level down to the surface, we obtain a possible range of 11.9–13.7 μbar for the surface pressure.

Pluto and its main satellite, Charon, occulted the same star on 2008 June 22. This event was observed from Australia and La Reunion Island, providing the east and north Charon Plutocentric offset in the sky plane (J2000): X = + 12,070.5 +/- 4 km (+ 546.2 +/- 0.2 mas), Y = + 4,576.3 +/- 24 km (+ 207.1 +/- 1.1 mas) at 19:20:33.82 UT on Earth, corresponding to JD 2454640.129964 at Pluto. This yields Charon's true longitude L = 153.483 +/- 0 degrees.071 in the satellite orbital plane (counted from the ascending node on J2000 mean equator) and orbital radius r = 19,564 +/- 14 km at that time. We compare this position to that predicted by (1) the orbital solution of Tholen & Buie (the "TB97" solution), (2) the PLU017 Charon ephemeris, and (3) the solution of Tholen et al. (the "T08" solution). We conclude that (1) our result rules out solution TB97, (2) our position agrees with PLU017, with differences of Delta L = + 0.073 +/- 0 degrees.071 in longitude, and Delta r = + 0.6 +/- 14 km in radius, and (3) while the difference with the T08 ephemeris amounts to only Delta L = 0.073 +/- 0 degrees.071 in longitude, it exhibits a significant radial discrepancy of Delta r = 61.3 +/- 14 km. We discuss this difference in terms of a possible image scale relative error of 3.35 x 10(-3) in the 2002-2003 Hubble Space Telescope images upon which the T08 solution is mostly based. Rescaling the T08 Charon semi-major axis, a = 19, 570.45 km, to the TB97 value, a=19636 km, all other orbital elements remaining the same ("T08/TB97" solution), we reconcile our position with the re-scaled solution by better than 12 km (or 0.55 mas) for Charon's position in its orbital plane, thus making T08/TB97 our preferred solution.

Pluto occultations are historically rare events, having been observed in 1988, 2002, 2006, and, as Pluto moves into the crowded Galactic plane, on several occasions in 2007. Here we present six results from our observations of the 2006 June 12 event from several sites in Australia and New Zealand. First, we show that Pluto's 2006 bulk atmospheric column abundance, as in 2002, is over twice the value measured in 1988, implying that nitrogen frost on Pluto's surface is 1.2–1.7 K warmer in 2006 than 1988 despite a 9% drop in incident solar flux. We measure a half-light shadow radius of 1216 ± 8.6 km in 2006, nominally larger than published values of 1213 ± 16 km measured in 2002. Given the current error bars, this latest half-light radius cannot discriminate between continued atmospheric growth or shrinkage, but it rules out several of the volatile transport scenarios modeled by Hansen & Paige. Second, we resolve spikes in the occultation light curve that are similar to those seen in 2002 and model the vertical temperature fluctuations that cause them. Third, we show that Pluto's upper atmosphere appears to hold a steady temperature of ∼100 K, as predicted from the methane thermostat model, even at latitudes where the methane thermostat is inoperative. This implies that energy transport rates are faster than radiational cooling rates. Fourth, this occultation has provided the first significant detection of a non-isothermal temperature gradient in Pluto's upper atmosphere also reported by Elliot et al., possibly the result of CO gas in Pluto's upper atmosphere. Fifth, we show that a haze-only explanation for Pluto's light curve is extremely unlikely; a thermal inversion is necessary to explain the observed light curve. And sixth, we derive an upper limit for the haze optical depth of 0.0023 in the zenith direction at average CCD wavelengths.

The 2006 June 12 occultation of the star P384.2 (2UCAC 26039859) by Pluto was observed from five sites in southeastern Australia with high-speed imaging photometers that produced time-series CCD images. Light curves were constructed from the image time series and fit by least-squares methods with model light curves. A new modeling procedure is presented that allows a simultaneous fit of the atmospheric parameters for Pluto and the astrometric parameters for the occultation to all of the light curves. Under the assumption of a clear atmosphere and using this modeling procedure to establish the upper atmosphere boundary condition, immersion and emersion temperature profiles were derived by inversion of the Siding Spring light curve, which had our best signal-to-noise ratio. Above ∼1230 km radius, atmospheric temperatures are ∼100 K and decrease slightly with altitude—the same as observed in 1988 and 2002. Below 1210 km, the temperature abruptly decreases with altitude (gradients ∼2.2 K km-1), which would reach the expected N2 surface-ice temperature of ∼40 K in the 1158-1184 km radius range. This structure is similar to that observed in 2002, but a much stronger thermal gradient (or stronger extinction) is implied by the 1988 light curve (which shows a "kink" or "knee" at 1210 km). The temperature profiles derived from inversion of the present data show good agreement with a physical model for Pluto's atmosphere selected from those presented by Strobel et al. (1996). Constraints derived from the temperature profiles (and considering the possibility of a deep troposphere) yield a value of 1152 ± 32 km for Pluto's surface radius. This value is compared with surface-radius values derived from the series of mutual occultations and eclipses that occurred in 1985-1989, and the limitations of both types of measurements for determining Pluto's surface radius are discussed. The radius of Pluto's atmospheric shadow at the half-intensity point is 1207.9 ± 8.5 km, the same as obtained in 2002 within measurement error. Values of the shadow radius cast by Pluto's atmosphere in 1988, 2002, and 2006 favor frost migration models in which Pluto's surface has low thermal inertia. Those models imply a substantial atmosphere when New Horizons flies by Pluto in 2015. Comparison of the shape of the stellar occultation light curves in 1988, 2002, and 2006 suggests that atmospheric extinction, which was strong in 1988 (15 months before perihelion), has been dissipating.