The description of the long-term dynamics of highly elliptic orbits under third-body perturbations may require an expansion of the disturbing function in series of the semi-major axes ratio up to higher orders. To avoid dealing with long series in trigonometric functions, we refer the motion to the apsidal frame and efficiently remove the short-period effects of this expansion in vectorial form up to an arbitrary order. We then provide the variation equations of the two fundamental vectors of the Keplerian motion by analogous vectorial recurrences, which are free from singularities and take a compact form useful for the numerical propagation of the flow in mean elements.

The southern hemisphere of comet 67P/Churyumov-Gerasimenko (67P) became observable by the Rosetta mission in March 2015, a few months before cometary southern vernal equinox. The Anhur region in the southern part of the comet's larger lobe was found to be highly eroded, enriched in volatiles, and highly active. We analyze high-resolution images of the Anhur region pre- and post-perihelion acquired by the OSIRIS imaging system on board the Rosetta mission. The Narrow Angle Camera is particularly useful for studying the evolution in Anhur in terms of morphological changes and color variations.}{Radiance factor images processed by the OSIRIS pipeline were coregistered, reprojected onto the 3D shape model of the comet, and corrected for the illumination conditions. We find a number of morphological changes in the Anhur region that are related to formation of new scarps; removal of dust coatings; localized resurfacing in some areas, including boulders displacements; and vanishing structures, which implies localized mass loss that we estimate to be higher than 50 million kg. The strongest changes took place in and nearby the Anhur canyon-like structure, where significant dust cover was removed, an entire structure vanished, and many boulders were rearranged. All such changes are potentially associated with one of the most intense outbursts registered by Rosetta during its observations, which occurred one day before perihelion passage. Moreover, in the niche at the foot of a new observed scarp, we also see evidence of water ice exposure that persisted for at least six months. The abundance of water ice, evaluated from a linear mixing model, is relatively high (> 20%). Our results confirm that the Anhur region is volatile-rich and probably is the area on 67P with the most pristine exposures near perihelion.

The Rosetta space probe accompanied comet 67P/Churyumov-Gerasimenko for more than two years, obtaining an unprecedented amount of unique data of the comet nucleus and inner coma. This work focuses identifying the source regions of faint jets and outbursts and on studying the spectrophotometric properties of some outbursts. We use observations acquired with the OSIRIS/NAC camera during July-October 2015, that is, close to perihelion. More than 200 jets of different intensities were identified directly on the nucleus. Some of the more intense outbursts appear spectrally bluer than the comet dark terrain in the vivible-to-near-infrared region. We attribute this spectral behavior to icy grains mixed with the ejected dust. Some of the jets have an extremely short lifetime. They appear on the cometary surface during the color sequence observations, and vanish in less than some few minutes after reaching their peak. We also report a resolved dust plume observed in May 2016 at a resolution of 55 cm/pixel, which allowed us to estimate an optical depth of $\sim$0.65 and an ejected mass of $\sim$ 2200 kg. We present the results on the location, duration, and colors of active sources on the nucleus of 67P from the medium-resolution (i.e., 6-10 m/pixel) images acquired close to perihelion passage. The observed jets are mainly located close to boundaries between different morphological regions. Jets depart not only from cliffs, but also from smooth and dust-covered areas, from fractures, pits, or cavities that cast shadows and favor the recondensation of volatiles. This study shows that faint jets or outbursts continuously contribute to the cometary activity close to perihelion passage, and that these events are triggered by illumination conditions. Faint jets or outbursts are not associated with a particular terrain type or morphology.

The algorithms used in the construction of a semi-analytical propagator for the long-term propagation of Highly Elliptical Orbits (HEO) are described. The software propagates mean elements and include the main gravitational and non-gravitational effects that may affect common HEO orbits, as, for instance, geostationary transfer orbits or Molniya orbits. Comparisons with numerical integration show that it provides good results even in extreme orbital configurations, as the case of SymbolX.

The non-integrability of the Hill problem makes that its global dynamics must be necessarily approached numerically. However, the analytical approach is feasible in the computation of relevant solutions. In particular, the nonlinear dynamics of the Hill problem close to the origin, and the libration point dynamics have been thoroughly investigated by perturbation methods. Out of the Hill sphere, the analytical approach is also feasible, at least in the case of distant retrograde orbits. Previous analytical investigations of this last case succeeded in the qualitative description of the dynamics, but they commonly failed in providing accurate results. This is a consequence of the essential dependance of the dynamics on elliptic functions, a fact that makes to progress in the perturbation approach beyond the lower orders of the solution really difficult. We propose an alternative perturbation approach that allows us to provide a very simple low order analytical solution in trigonometric functions, on the one hand, and, while still depending on special functions, to compute higher orders of the solution, on the other.

The Southern hemisphere of the 67P/Churyumov-Gerasimenko comet has become visible from Rosetta only since March 2015. It was illuminated during the perihelion passage and therefore it contains the regions that experienced the strongest heating and erosion rate, thus exposing the subsurface most pristine material. In this work we investigate, thanks to the OSIRIS images, the geomorphology, the spectrophotometry and some transient events of two Southern hemisphere regions: Anhur and part of Bes. Bes is dominated by outcropping consolidated terrain covered with fine particle deposits, while Anhur appears strongly eroded with elongated canyon-like structures, scarp retreats, different kinds of deposits, and degraded sequences of strata indicating a pervasive layering. We discovered a new 140 m long and 10 m high scarp formed in the Anhur/Bes boundary during/after the perihelion passage, close to the area where exposed CO$_2$ and H$_2$O ices were previously detected. Several jets have been observed originating from these regions, including the strong perihelion outburst, an active pit, and a faint optically thick dust plume. We identify several areas with a relatively bluer slope (i.e. a lower spectral slope value) than their surroundings, indicating a surface composition enriched with some water ice. These spectrally bluer areas are observed especially in talus and gravitational accumulation deposits where freshly exposed material had fallen from nearby scarps and cliffs. The investigated regions become spectrally redder beyond 2 au outbound when the dust mantle became thicker, masking the underlying ice-rich layers.

We present the results of the photometric and spectrophotometric properties of the 67P/Churyumov-Gerasimenko nucleus derived with the OSIRIS instrument during the closest fly-by over the comet, which took place on 14 th February 2015 at a distance of {\~} 6 km from the surface. Several images covering the 0{\deg}-33{\deg} phase angle range were acquired, and the spatial resolution achieved was 11 cm/pxl. The flown-by region is located on the big lobe of the comet, near the borders of the Ash, Apis and Imhotep regions. Our analysis shows that this region features local heterogeneities at the decimetre scale. We observed difference of reflectance up to 40{\%} between bright spots and sombre regions, and spectral slope variations up to 50{\%}. The spectral reddening effect observed globally on the comet surface by Fornasier et al. (2015) is also observed locally on this region, but with a less steep behaviour. We note that numerous metre-sized boulders, which exhibit a smaller opposition effect, also appear spectrally redder than their surroundings. In this region, we found no evidence linking observed bright spots to exposed water-ice-rich material. We fitted our dataset using the Hapke 2008 photometric model. The region overflown is globally as dark as the whole nucleus (geometric albedo of 6.8{\%}) and it has a high porosity value in the uppermost-layers (86{\%}). These results of the photometric analysis at a decimetre scale indicate that the photometric properties of the flown-by region are similar to those previously found for the whole nucleus.

In 1883, on the 12th and 13th of August, Mexican astronomer Jose A. y Bonilla observed several objects passing in front of the solar disk. In 1886 in the L'Astronomie magazine, he reported his observations without providing a hypothesis explaining the registered phenomena. Our objective in this work is to interpret, with current knowledge, what he observed in Zacatecas. Our working hypothesis is that what Bonilla observed in 1883 was a highly fragmented comet, in an approach almost flush to the Earth's surface. The fragmentation of the comet's nucleus is a phenomenon known since the XIX century. Using the results reported by Bonilla, we can estimate the distance at which the objects approach to the Earth's surface, their size, their mass and total mass of the comet before fragmentation. According to our calculations, the distance at which the objects passed over the Earth's surface, was between 538 km and 8,062 km, the width of the objects was between 46 m and 795 m and its length between 68 m and 1,022 m, the object's mass was between 5.58e8 kg and 2.5e12 kg. Finally, the mass of the original comet, before fragmentation, was between 1.83e12 and 8.19e15 kg, i.e., between 2e-3 and 8.19 times the mass of Halley Comet.