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.

We present results derived from four stellar occultations by the plutino object (208996) 2003~AZ$_{84}$, detected at January 8, 2011 (single-chord event), February 3, 2012 (multi-chord), December 2, 2013 (single-chord) and November 15, 2014 (multi-chord). Our observations rule out an oblate spheroid solution for 2003~AZ$_{84}$'s shape. Instead, assuming hydrostatic equilibrium, we find that a Jacobi triaxial solution with semi axes $(470 \pm 20) \times (383 \pm 10) \times (245 \pm 8)$~km % axis ratios $b/a= 0.82 \pm 0.05$ and $c/a= 0.52 \pm 0.02$, can better account for all our occultation observations. Combining these dimensions with the rotation period of the body (6.75~h) and the amplitude of its rotation light curve, we derive a density $\rho=0.87 \pm 0.01$~g~cm$^{-3}$ a geometric albedo $p_V= 0.097 \pm 0.009$. A grazing chord observed during the 2014 occultation reveals a topographic feature along 2003~AZ$_{84}$'s limb, that can be interpreted as an abrupt chasm of width $\sim 23$~km and depth $> 8$~km or a smooth depression of width $\sim 80$~km and depth $\sim 13$~km (or an intermediate feature between those two extremes).

The orbital displacement of the Moon is known to a precision on the order of centimeters from the lunar laser ranging data, while the lunar profile's confidence level is still only about 200 m. The lunar motion is used to measure the solar diameter during central eclipses with accurate timing of Baily's beads at the umbral path limits. The onset or disappearance of a Baily's bead is due to the alignment between lunar limb's valleys and solar photosphere: this timing is determined outside our atmosphere. Due to scintillation, naked eye evaluations of such phenomena can be systematically different with respect to telescopes. We analyze the case of Venus' occultations as bead-like situation, in order to infer considerations on naked eye reports in historical total eclipse reports.