We present here the Calar Alto Schmidt-Lemaitre Telescope (CASTLE) concept, a technology demonstrator for curved detectors, that will be installed at the Calar Alto Observatory (Spain). This telescope has a wide field of view (2.36x1.56 deg^2) and a design, optimised to generate a Point Spread Function with very low level wings and reduced ghost features, which makes it considerably less susceptible to several systematic effects usually affecting similar systems. These characteristics are particularly suited to study the low surface brightness Universe. CASTLE will be able to reach surface brightness orders of magnitude fainter than the sky background level and observe the extremely extended and faint features around galaxies such as tidal features, stellar halos, intra-cluster light, etc. CASTLE will also be used to search and detect astrophysical transients such as gamma ray bursts (GRB), gravitational wave optical counterparts, neutrino counterparts, etc. This will increase the number of precisely localized GRBs from 20% to 60% (in the case of Fermi/GMB GRBs).

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.

On 28th January 2018, the large Trans-Neptunian Object (TNO) 2002TC302 occulted a m$_v= $15.3 star with ID 130957813463146112 in the Gaia DR2 stellar catalog. 12 positive occultation chords were obtained from Italy, France, Slovenia and Switzerland. Also, 4 negative detections were obtained near the north and south limbs. This represents the best observed stellar occultation by a TNO other than Pluto, in terms of the number of chords published thus far. From the 12 chords, an accurate elliptical fit to the instantaneous projection of the body, compatible with the near misses, can be obtained. The resulting ellipse has major and minor axes of 543 $\pm$ 18 km and 460 $\pm$ 11 km, respectively, with a position angle of 3 $\pm$ 1 degrees for the minor axis. This information, combined with rotational light curves obtained with the 1.5m telescope at Sierra Nevada Observatory and the 1.23m telescope at Calar Alto observatory, allows us to derive possible 3D shapes and density estimations for the body, based on hydrostatic equilibrium assumptions. The effective area equivalent diameter is $\sim$ 84 km smaller than the radiometrically derived diameter using thermal data from Herschel and Spitzer Space Telescopes. This might indicate the existence of an unresolved satellite of up to $\sim$ 300 km in diameter, to account for all the thermal flux, although the occultation and thermal diameters are compatible within their error bars given the considerable uncertainty of the thermal results. The existence of a potential satellite also appears to be consistent with other ground-based data presented here. From the effective occultation diameter combined with H$_V$ measurements we derive a geometric albedo of 0.147 $\pm$ 0.005, which would be somewhat smaller if 2002TC302 has a satellite. The best occultation light curves do not show any signs of ring features or any signatures of a global atmosphere.

Non-resolved thermal infrared observations enable studies of thermal and physical properties of asteroid surfaces provided the shape and rotational properties of the target are well determined via thermo-physical models. We used calibration-programme Herschel PACS data (70, 100, 160 $\mu$m) and state-of-the-art shape models derived from adaptive-optics observations and/or optical light curves to constrain for the first time the thermal inertia of twelve large main-belt asteroids. We also modelled previously well-characterised targets such as (1) Ceres or (4) Vesta as they constitute important benchmarks. Using the scale as a free parameter, most targets required a re-scaling $\sim$5\% consistent with what would be expected given the absolute calibration error bars. This constitutes a good cross-validation of the scaled shape models, although some targets required larger re-scaling to reproduce the IR data. We obtained low thermal inertias typical of large main belt asteroids studied before, which continues to give support to the notion that these surfaces are covered by fine-grained insulating regolith. Although the wavelengths at which PACS observed are longwards of the emission peak for main-belt asteroids, they proved to be extremely valuable to constrain size and thermal inertia and not too sensitive to surface roughness. Finally, we also propose a graphical approach to help examine how different values of the exponent used for scaling the thermal inertia as a function of heliocentric distance (i.e. temperature) affect our interpretation of the results.

From CCD observations carried out with different telescopes, we present short-term photometric measurements of the large trans-Neptunian object Varuna in 10 epochs, spanning around 19 years. We observe that the amplitude of the rotational light-curve has changed considerably during this period of time from 0.41 to 0.55 mag. In order to explain this variation, we constructed a model in which Varuna has a simple triaxial shape, assuming that the main effect comes from the change of the aspect angle as seen from Earth, due to Varuna's orbital motion in the 19-year time span. The best fits to the data correspond to a family of solutions with axial ratios b/a between 0.56 and 0.60. This constrains the pole orientation in two different ranges of solutions presented here as maps. Apart from the remarkable variation of the amplitude, we have detected changes in the overall shape of the rotational light-curve over shorter time scales. After the analysis of the periodogram of the residuals to a 6.343572 h double-peaked rotational light-curve fit, we find a clear additional periodicity. We propose that these changes in the rotational light-curve shape are due to a large and close-in satellite whose rotation induces the additional periodicity. The peak-to-valley amplitude of this oscillation is in the order of 0.04 mag. We estimate that the satellite orbits Varuna with a period of 11.9819 h (or 23.9638 h), assuming that the satellite is tidally locked, at a distance of ~ 1300 km (or ~ 2000 km) from Varuna, outside the Roche limit.

We present results from three world-wide campaigns that resulted in the detections of two single-chord and one multi-chord stellar occultations by the Plutino object (84922) 2003~VS$_2$. From the single-chord occultations in 2013 and 2014 we obtained accurate astrometric positions for the object, while from the multi-chord occultation on November 7th, 2014, we obtained the parameters of the best-fitting ellipse to the limb of the body at the time of occultation. We also obtained short-term photometry data for the body in order to derive its rotational phase during the occultation. The rotational light curve present a peak-to-peak amplitude of 0.141 $\pm$ 0.009 mag. This allows us to reconstruct the three-dimensional shape of the body, with principal semi-axes $a = 313.8 \pm 7.1$ km, $b = 265.5^{+8.8}_{-9.8}$ km, and $c = 247.3^{+26.6}_{-43.6}$ km, which is not consistent with a Jacobi triaxial equilibrium figure. The derived spherical volume equivalent diameter of $548.3 ^{+29.5}_{-44.6}$ km is about 5\% larger than the radiometric diameter of 2003~VS$_2$ derived from Herschel data of $523 \pm 35$ km, but still compatible with it within error bars. From those results we can also derive the geometric albedo ($0.123 ^{+0.015}_{-0.014}$) and, under the assumption that the object is a Maclaurin spheroid, the density $\rho = 1400^{+1000}_{-300}$ for the plutino. The disappearances and reappearances of the star during the occultations do not show any compelling evidence for a global atmosphere considering a pressure upper limit of about 1 microbar for a pure nitrogen atmosphere, nor secondary features (e.g. rings or satellite) around the main body.

The search for Earth-like planets around late-type stars using ultra-stable spectrographs requires a very precise characterization of the stellar activity and the magnetic cycle of the star, since these phenomena induce radial velocity (RV) signals that can be misinterpreted as planetary signals. Among the nearby stars, we have selected Barnard's Star (Gl 699) to carry out a characterization of these phenomena using a set of spectroscopic data that covers about 14.5 years and comes from seven different spectrographs: HARPS, HARPS-N, CARMENES, HIRES, UVES, APF, and PFS; and a set of photometric data that covers about 15.1 years and comes from four different photometric sources: ASAS, FCAPT-RCT, AAVSO, and SNO. We have measured different chromospheric activity indicators (H$\alpha$, Ca~{\sc II}~HK and Na I D), as well as the FWHM of the cross-correlation function computed for a sub-set of the spectroscopic data. The analysis of Generalized Lomb-Scargle periodograms of the time series of different activity indicators reveals that the rotation period of the star is 145 $\pm$ 15 days, consistent with the expected rotation period according to the low activity level of the star and previous claims. The upper limit of the predicted activity-induced RV signal corresponding to this rotation period is about 1 m/s. We also find evidence of a long-term cycle of 10 $\pm$ 2 years that is consistent with previous estimates of magnetic cycles from photometric time series in other M stars of similar activity levels. The available photometric data of the star also support the detection of both the long-term and the rotation signals.

The study of the visible colours of the trans-Neptunian objects opened a discussion almost 20 years ago which, in spite of the increase in the amount of available data, seems far from subside. Visible colours impose constraints to the current theories of the early dynamical evolution of the Solar System such as the environment of formation, initial surface composition, and how (if) they were scattered to regions closer to the inner planets. In this paper we present an updated version of our database of absolute colours and relative phase coefficients for 117 objects. We define the absolute colours as the difference of the absolute magnitudes $H_V-H_R$, and the relative phase coefficient as the difference of the slopes of the phase curves $\Delta\beta$. These were obtained joining our own observations plus data from the literature. The methodology has been introduced in previous works and here we expand in some interesting results, in particular the strong anti-correlation found between $H_V-H_R$ and $\Delta\beta$, which means that redder objects have steeper phase curves in the R filter, while bluer objects have steeper phase curves in the V filter. We analyse a series of results published in the literature in view of our database, which is free of phase effects, and show that their statistical meaning is not very strong. We point out that phase-colouring and observational errors play an important role in the understanding of these proposed relationships.

We discuss here a lunar impact flash recorded during the total lunar eclipse that occurred on 2019 January 21, at 4h 41m 38.09 +- 0.01 s UT. This is the first time ever that an impact flash is unambiguously recorded during a lunar eclipse and discussed in the scientific literature, and the first time that lunar impact flash observations in more than two wavelengths are reported. The impact event was observed by different instruments in the framework of the MIDAS survey. It was also spotted by casual observers that were taking images of the eclipse. The flash lasted 0.28 seconds and its peak luminosity in visible band was equivalent to the brightness of a mag. 4.2 star. The projectile hit the Moon at the coordinates 29.2 +- 0.3 $^{\circ}$S, 67.5 +- 0.4 $^{\circ}$W. In this work we have investigated the most likely source of the projectile, and the diameter of the new crater generated by the collision has been calculated. In addition, the temperature of the lunar impact flash is derived from the multiwavelength observations. These indicate that the blackbody temperature of this flash was of about 5700 K.

In terms of scientific output, the best way to study solar system bodies is sending spacecraft to make in-situ measurements or to observe at close distance. Probably, the second best means to learn about important physical properties of solar system objects is through stellar occultations. By combining occultation observations from several sites, size and shape can be derived with kilometric accuracy. Also, atmospheric properties can be derived if the body has an atmosphere. Furthermore, the technique can detect rings and even satellites (although rarely) around the main body. Except for the very special cases of Pluto and Charon, stellar occultations by Transneptunian Objects (TNOs) had never been observed until October 2009. This was because the ephemeris of the TNOs have much larger uncertainties than their angular diameters (typically of the order of ~10 milliarcsecond) and also because stellar catalogs were not accurate to the milliarcsecond level. Despite the difficulties, at the time of this writing, 43 occultations by 22 different Trans-Neptunian Objects, and 17 occultations by 5 Centaurs have been detected thanks to the efforts of several teams. Due to the complications of accurately predicting and observing these events, most of the successes have been achieved through wide international collaboration, which is a key issue to succeed in observing stellar occultations by TNOs. Multichord occultations are typically detected at a rate of ~3 per year on average, whereas the majority of the observed occultations are single-chord detections, which means that only one site detects the occultation. In these cases, no tight constraints on size and shape can be derived from those observations alone. Here we review most of the aspects involved in the complex process to successfully observe occultations, and present some of the lessons learned.

A group of trans-Neptunian objects (TNO) are dynamically related to the dwarf planet 136108 Haumea. Ten of them show strong indications of water ice on their surfaces, are assumed to have resulted from a collision, and are accepted as the only known TNO collisional family. Nineteen other dynamically similar objects lack water ice absorptions and are hypothesized to be dynamical interlopers. We have made observations to determine sizes and geometric albedos of six of the accepted Haumea family members and one dynamical interloper. Ten other dynamical interlopers have been measured by previous works. We compare the individual and statistical properties of the family members and interlopers, examining the size and albedo distributions of both groups. We also examine implications for the total mass of the family and their ejection velocities. We use far-infrared space-based telescopes to observe the target TNOs near their thermal peak and combine these data with optical magnitudes to derive sizes and albedos using radiometric techniques. We determine the power-law slope of ejection velocity as a function of effective diameter. The detected Haumea family members have a diversity of geometric albedos $\sim$ 0.3-0.8, which are higher than geometric albedos of dynamically similar objects without water ice. The median geometric albedo for accepted family members is $p_V=0.48_{-0.18}^{+0.28}$, compared to 0.08$_{-0.05}^{+0.07}$ for the dynamical interlopers. In the size range $D=175-300$ km, the slope of the cumulative size distribution is $q$=3.2$_{-0.4}^{+0.7}$ for accepted family members, steeper than the $q$=2.0$\pm$0.6 slope for the dynamical interlopers with D$< $500 km. The total mass of Haumea's moons and family members is 2.4% of Haumea's mass. The ejection velocities required to emplace them on their current orbits show a dependence on diameter, with a power-law slope of 0.21-0.50.

We studied the rotational properties of the dwarf planet Makemake. The photometric observations were carried out at different telescopes between 2006 and 2017. Most of the measurements were acquired in BVRI broad-band filters of a standard Johnson-Cousins photometric system. We found that Makemake rotates more slowly than was previously reported. A possible lightcurve asymmetry suggests a double-peaked period of P = 22.8266$\pm$0.0001~h. A small peak-to-peak lightcurve amplitude in R-filter A = 0.032$\pm$0.005 mag implies an almost spherical shape or near pole-on orientation. We also measured BVRI colours and the R-filter phase-angle slope and revised the absolute magnitudes. The absolute magnitude of Makemake has remained unchanged since its discovery in 2005. No direct evidence of a newly discovered satellite was found in our photometric data; however, we discuss the possible existence of another larger satellite.

From 1988 to 2016, several stellar occultations have been observed to characterize Pluto's atmosphere and its evolution (Meza et al, 2019). From each stellar occultation, an accurate astrometric position of Pluto at the observation epoch is derived. These positions mainly depend on the position of the occulted star and the precision of the timing. We present Pluto's astrometric positions derived from 19 occultations from 1988 to 2016 (11 from Meza et al. (2019) and 8 from other publications). Using Gaia DR2 for the positions of the occulted stars, the accuracy of these positions is estimated to 2-10~milliarcsec depending on the observation circumstances. From these astrometric positions, we derive an updated ephemeris of Pluto's system barycentre using the NIMA code (Desmars et al., 2015). The astrometric positions are derived by fitting the occultation's light curves by a model of Pluto's atmosphere. The fits provide the observed position of the body's centre for a reference star position. Other publications usually provide circumstances of the occultation such as the coordinates of the stations, the timing, and the impact parameter (i.e. the closest distance between the station and the centre of the shadow). From these parameters, we use a procedure based on the Bessel method to derive an astrometric position. We derive accurate Pluto's astrometric positions from 1988 to 2016. These positions are used to refine the orbit of Pluto'system barycentre providing an ephemeris, accurate to the milliarcsec level, over the period 2000-2020, allowing better predictions for future stellar occultations.

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 studied 93 asteroid pairs. We estimated times elapsed since separation of pair members that are between 7*10^3 and a few 10^6 yr. We derived the rotation periods for all the primaries and a sample of secondaries. We derived the absolute magnitude differences of the asteroid pairs that provide their mass ratios. We refined their WISE geometric albedos and estimated their taxonomic classifications. For 17 pairs, we determined their pole positions. In 2 pairs where we obtained the spin poles for both components, we saw the same sense of rotation for both components and constrained the angles between their original spin vectors at the time of their separation. We found that the primaries of 13 pairs are actually binary or triple systems, i.e., they have one or two bound secondaries (satellites). As by-product, we found 3 new young asteroid clusters (each of them consisting of three known asteroids on highly similar orbits). We compared the obtained asteroid pair data with theoretical predictions and discussed their implications. We found that 86 of the 93 studied pairs follow the trend of primary rotation period vs mass ratio that was found by Pravec et al. (2010). Of the 7 outliers, 3 appear insignificant (may be due to our uncertain or incomplete knowledge), but 4 are high mass ratio pairs that were unpredicted by the theory of asteroid pair formation by rotational fission. We discuss a (remotely) possible way that they could be created by rotational fission of flattened parent bodies followed by re-shaping of the formed components. The 13 pairs with binary primaries are particularly interesting systems that place important constraints on formation and evolution of asteroid pairs. We present two hypotheses for their formation: The pairs having both bound and unbound secondaries could be `failed asteroid clusters', or they could be formed by a cascade primary spin fission process.

A recent occultation measurement of the dwarf planet Haumea (Ortiz et al. 2017) revealed an elongated shape with the longest axis comparable to Pluto's mean diameter. The chords also indicate a ring around Haumea's equatorial plane, where its largest moon, Hi'iaka, is also located. The Haumea occultation size estimate (equivalent diameter 1595 km) is larger than previous radiometric solutions (in the range between 1150 and 1350 km), which lowers the object's density to about 1.8 gcm-3, a value closer to the densities of other large TNOs. We present unpublished and also reprocessed Herschel and Spitzer MIR and FIR measurements. We compare 100 and 160um thermal lightcurve amplitudes - originating from Haumea itself - with models of the total measured system fluxes (ring, satellite, Haumea) from 24-350um. The combination with results derived from the occultation measurements show that Haumea's crystalline water ice surface must have a thermal inertia of ~5 SIu. We also have indications that the satellites (at least Hi'iaka) must have high geometric albedos >=0.5, otherwise the derived thermal amplitude would be inconsistent with the total measured system fluxes at 24, 70, 100, 160, 250, and 350um. The high albedos imply sizes of about 300 and 150 km for Hi'iaka and Namaka, respectively, indicating unexpectedly high densities >1.0 gcm-3 for TNOs this small, and the assumed collisional formation from Haumea's icy crust. We also estimated the thermal emission of the ring for the time period 1980-2030, showing that the contribution during the Spitzer and Herschel epochs was small, but not negligible. Due to the progressive opening of the ring plane, the ring emission will be increasing in the next decade when JWST is operational. In the MIRI 25.5um band it will also be possible to obtain a very high-quality thermal lightcurve to test the derived Haumea properties.

The trans-Neptunian objects (TNOs) are small Solar System bodies at large distances from the Sun. As such, their physical properties are difficult to measure. Accurate determination of their physical parameters is essential to model and theorize the actual composition and distribution of the population, and to improve our understanding of the formation and evolution of the Solar System. The objective of this work is to construct phase curves in two filters, $V$ and $R$, of a large TNO sample obtaining absolute magnitudes ($H$) and phase coefficients ($\beta$), and study possible relations between them and other physical parameters (orbital elements, sizes, and albedos). We used our own data, together with data from the literature, to create the phase curves assuming an overall linear trend. We obtained new magnitudes for 35 TNOs, 27 in the V filter and 35 in the R filter. These magnitudes, together with data from the literature, allowed us to obtain absolutes magnitudes, 114 in the V filter and 113 in the R filter, of which 106 have both. From the search for correlations we found a strong anticorrelation between $H_V-H_R$ and $\Delta \beta=\beta_V-\beta_R$, which is probably more related to surface structure than to composition or size of the objects.

We analyze the outburst experienced by the September $\varepsilon$-Perseid meteor shower on 9 September 2013. As a result of our monitoring the atmospheric trajectory of 60 multi-station events observed over Spain was obtained and accurate orbital data were derived from them. On the basis of these orbits, we have tried to determine the likely parent body of this meteoroid stream by employing orbital dissimilarity criteria. In addition, the emission spectra produced by two events belonging to this meteor shower were also recorded. The analysis of these spectra has provided information about the chemical nature of their progenitor meteoroids. We also present an estimation of the tensile strength for these particles.

We report the first analysis of a flash produced by the impact of a meteoroid on the lunar surface and recorded both in the near-infrared and in the visible. Despite the fact that similar data have been recently published by other team during the refereeing process of our manuscript (Bonanos et al. 2018), our result still forms the first measurement of the temperature of a telescopic lunar impact flash (Madiedo and Ortiz 2016, 2018). The flash exhibited a peak magnitude of 5.1 $\pm$ 0.3 in the near-infrared I band and 7.3 $\pm$ 0.2 in the visible, and the total duration of the event in these bands was 0.20 s and 0.18 s, respectively. The origin of the meteoroid was investigated, and we inferred that the most likely scenario is that the impactor that belonged to the sporadic background. The analysis of this event has provided for the first time an estimation of the emission efficiency in the near-infrared {\eta}I for sporadic meteoroids impacting the Moon. We have determined that this efficiency is around 56% higher than in the visible band and we have found a maximum impact plume temperature of ~4000 K at the initial phase followed by temperatures of around 3200 K after the peak brightness. The size of the crater produced as a consequence of this impact is also calculated.

The centaur 95P/(2060) Chiron is showing comet-like activity since its discovery, but the mass-loss mechanisms triggering its activity remained unexplained. Although the collision rates in the centaur region are expected to be very low, and impacts are thought not to be responsible for the mass-loss, since the recent indications that Chiron might possess a ring similar to Chariklo's, and assuming that there is debris orbiting around, the impact triggered mass-loss mechanism should not be excluded as a possible cause of its activity. From time series observations collected on Calar Alto Observatory in Spain between 2014 and 2016, we found that the photometric scatter in Chiron's data is larger than a control star's scatter, indicating a possible microactivity, possibly caused by debris falling back to Chiron's surface and lifting small clouds of material. We also present rotational light curves, and measurements of Chiron's absolute magnitudes, that are consistent with the models supporting the presumption that Chiron possesses rings. By co-adding the images acquired in 2015, we have detected a $ \sim $5 arcsec long tail, showing a surface brightness of 25.3 mag(V)/arcsec$^{2}$.

Proxima Centauri, the star closest to our Sun, is known to host at least one terrestrial planet candidate in a temperate orbit. Here we report the ALMA detection of the star at 1.3 mm wavelength and the discovery of a belt of dust orbiting around it at distances ranging between 1 and 4 au, approximately. Given the low luminosity of the Proxima Centauri star, we estimate a characteristic temperature of about 40 K for this dust, which might constitute the dust component of a small-scale analog to our solar system Kuiper belt. The estimated total mass, including dust and bodies up to 50 km in size, is of the order of 0.01 Earth masses, which is similar to that of the solar Kuiper belt. Our data also show a hint of warmer dust closer to the star. We also find signs of two additional features that might be associated with the Proxima Centauri system, which, however, still require further observations to be confirmed: an outer extremely cold (about 10 K) belt around the star at about 30 au, whose orbital plane is tilted about 45 degrees with respect to the plane of the sky; and additionally, we marginally detect a compact 1.3 mm emission source at a projected distance of about 1.2 arcsec from the star, whose nature is still unknown.

The combination of visible and thermal data from the ground and astrophysics space missions is key to improving the scientific understanding of near-Earth, main-belt, trojans, centaurs, and trans-Neptunian objects. To get full information on a small sample of selected bodies we combine different methods and techniques: lightcurve inversion, stellar occultations, thermophysical modeling, radiometric methods, radar ranging and adaptive optics imaging. The SBNAF project will derive size, spin and shape, thermal inertia, surface roughness, and in some cases bulk densities and even internal structure and composition, for objects out to the most distant regions in the Solar System. The applications to objects with ground-truth information allows us to advance the techniques beyond the current state-of-the-art and to assess the limitations of each method. We present results from our project's first phase: the analysis of combined Herschel-KeplerK2 data and Herschel-occultation data for TNOs; synergy studies on large MBAs from combined high-quality visual and thermal data; establishment of well-known asteroids as celestial calibrators for far-infrared, sub-millimetre, and millimetre projects; first results on near-Earth asteroids properties from combined lightcurve, radar and thermal measurements, as well as the Hayabusa-2 mission target characterisation. We also introduce public web-services and tools for studies of small bodies in general.

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.

Context. The so-called Barbarian asteroids share peculiar, but common polarimetric properties, probably related to both their shape and composition. They are named after (234) Barbara, the first on which such properties were identified. As has been suggested, large scale topographic features could play a role in the polarimetric response, if the shapes of Barbarians are particularly irregular and present a variety of scattering/incidence angles. This idea is supported by the shape of (234) Barbara, that appears to be deeply excavated by wide concave areas revealed by photometry and stellar occultations. Aims. With these motivations, we started an observation campaign to characterise the shape and rotation properties of Small Main- Belt Asteroid Spectroscopic Survey (SMASS) type L and Ld asteroids. As many of them show long rotation periods, we activated a worldwide network of observers to obtain a dense temporal coverage. Methods. We used light-curve inversion technique in order to determine the sidereal rotation periods of 15 asteroids and the con- vergence to a stable shape and pole coordinates for 8 of them. By using available data from occultations, we are able to scale some shapes to an absolute size. We also study the rotation periods of our sample looking for confirmation of the suspected abundance of asteroids with long rotation periods. Results. Our results show that the shape models of our sample do not seem to have peculiar properties with respect to asteroids with similar size, while an excess of slow rotators is most probably confirmed.

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).

We present a summary of the campaign of remote observations that supported the European Space Agency's Rosetta mission. Telescopes across the globe (and in space) followed comet 67P/Churyumov-Gerasimenko from before Rosetta's arrival until nearly the end of mission in September 2016. These provided essential data for mission planning, large-scale context information for the coma and tails beyond the spacecraft, and a way to directly compare 67P with other comets. The observations revealed 67P to be a relatively `well behaved' comet, typical of Jupiter family comets and with activity patterns that repeat from orbit-to-orbit. Comparison between this large collection of telescopic observations and the in situ results from Rosetta will allow us to better understand comet coma chemistry and structure. This work is just beginning as the mission ends -- in this paper we present a summary of the ground-based observations and early results, and point to many questions that will be addressed in future studies.

Time series observations of the dwarf planet Haumea and the Plutinos 2003VS2 and 2003AZ84 with Herschel/PACS are presented in this work. Thermal emission of these trans-Neptunian objects were acquired as part of the TNOs are Cool Herschel Space Observatory key programme. We search for the thermal light curves at 100 and 160um of Haumea and 2003AZ84, and at 70 and 160um for 2003VS2 by means of photometric analysis of the PACS data. The goal of this work is to use these thermal light curves to obtain physical and thermophysical properties of these icy Solar System bodies. Haumea's thermal light curve is clearly detected at 100 and 160um. The effect of the reported dark spot is apparent at 100um. Different thermophysical models were applied to these light curves, varying the thermophysical properties of the surface within and outside the spot. Although no model gives a perfect fit to the thermal observations, results imply an extremely low thermal inertia (< 0.5 MKS) and a high phase integral (> 0.73) for Haumea's surface. We note that the dark spot region appears to be only weakly different from the rest of the object, with modest changes in thermal inertia and/or phase integral. The thermal light curve of 2003VS2 is not firmly detected at 70 and at 160um but a thermal inertia of 2+/-0.5 MKS can be derived from these data. The thermal light curve of 2003AZ84 is not firmly detected at 100um. We apply a thermophysical model to the mean thermal fluxes and to all the Herschel/PACS and Spitzer/MIPS thermal data of 2003AZ84, obtaining a close to pole-on orientation as the most likely for this TNO. For the three TNOs, the thermal inertias derived from light curve analyses or from the thermophysical analysis of the mean thermal fluxes confirm the generally small or very small surface thermal inertias of the TNO population, which is consistent with a statistical mean value of 2.5+/-0.5 MKS.

Context. Precise determination of stellar masses is necessary to test the validity of pre-main-sequence (PMS) stellar evolutionary models, whose predictions are in disagreement with measurements for masses below 1.2 Msun. To improve such a test, and based on our previous studies, we selected the AB Doradus moving group (AB Dor-MG) as the best-suited association on which to apply radio-based high-precision astrometric techniques to study binary systems. Aims. We seek to determine precise estimates of the masses of a set of stars belonging to the AB Dor-MG using radio and infrared observations. Methods. We observed in phase-reference mode with the Very Large Array (VLA) at 5 GHz and with the European VLBI Network (EVN) at 8.4 GHz the stars HD 160934, EK Dra, PW And, and LO Peg. We also observed some of these stars with the near-infrared CCD AstraLux camera at the Calar Alto observatory to complement the radio observations. Results. We determine model-independent dynamical masses of both components of the star HD 160934, A and c, which are 0.70+/-0.07 Msun and 0.45+/-0.04 Msun , respectively. We revised the orbital parameters of EK Dra and we determine a sum of the masses of the system of 1.38+/-0.08 Msun. We also explored the binarity of the stars LO Peg and PW And. Conclusions. We found observational evidence that PMS evolutionary models underpredict the mass of PMS stars by 10%-40%, as previously reported by other authors. We also inferred that the origin of the radio emission must be similar in all observed stars, that is, extreme magnetic activity of the stellar corona that triggers gyrosynchrotron emission from non-thermal, accelerated electrons.

The discovery that the centaur (10199) Chariklo possesses a ring system opens questions about their origin. We here asses the plausibility of different scenarios for the origin of the observed ring system. We first consider the possibility that the material of the ring originated in the disruption of a satellite that had reached a critical distance from the centaur. We discuss the conditions for the putative satellite to approach the centaur as a consequence of tidal interaction. A three-body encounter is also considered as a transport mechanism. In addition, we study the case in which the ring is formed by the ejecta of a cratering collision on the centaur and we constrain the collision parameters and the size of the resulting crater of the event. Finally, we consider that the ring material originates from a catastrophic collision between a background object and a satellite located at a distance corresponding to the the current location of the ring. We compute the typical timescales for these scenarios. We estimate that in order to be tidally disrupted a satellite would have had to be larger than approximately 6.5 km at the location of the rings. However the tidal interaction is rather weak for objects of the size of outer solar system bodies at the ring location, therefore we considered other more effective mechanisms by which a satellite migt have approached the centaur. Collisonal scenarios are both physically plausible for the formation, but semianalytical estimations indicate that the probability of the corresponding collisions is low under current conditions

2007OR10 is currently the third largest known dwarf planet in the transneptunian region, with an effective radiometric diameter of ~1535 km. It has a slow rotation period of ~45 h that was suspected to be caused by tidal interactions with a satellite undetected at that time. Here we report on the discovery of a likely moon of 2007OR10, identified on archival Hubble Space Telescope WFC3/UVIS system images. Although the satellite is detected at two epochs, this does not allow an unambiguous determination of the orbit and the orbital period. A feasible 1.5-5.8x10^21 kg estimate for the system mass leads to a likely 35 to 100 d orbital period. The moon is about 4.2m fainter than 2007OR10 in HST images that corresponds to a diameter of 237 km assuming equal albedos with the primary. Due to the relatively small size of the moon the previous size and albedo estimates for the primary remains unchanged. With this discovery all trans-Neptunian objects larger than 1000 km are now known to harbour satellites, an important constraint for moon formation theories in the young Solar system.

We present time series photometry of Bienor in four observation campaigns from 2013 to 2016 and compare them with previous observations in the literature dating back to 2000. The results show a remarkable decline in the amplitude of the rotational light curve and in the absolute magnitude. This suggests that the angle between the rotation axis and the line of sight has changed noticeably during the last 16 years as Bienor orbits the Sun. From the light curve amplitude data we are able to determine the orientation of the rotation axis of Bienor ($\beta_p=50\pm3^{\circ}$, $\lambda_p=35\pm8^{\circ}$). We are also able to constrain the $b/a$ axial ratio of a triaxial Jacobi ellipsoidal body (with semi-axis $a>b>c$). The best fit is for $b/a=0.45\pm0.05$, which corresponds to a density value of $594^{+47}_{-35}$ kg$\,$m$^{-3}$ under the usual assumption of hydrostatic equilibrium and given that Bienor's rotational period is 9.17 h. However, the absolute magnitude of Bienor at several epochs is not well reproduced. We tested several explanations such as relaxing the hydrostatic equilibrium constraint, a large North-South asymmetry in the surface albedo of Bienor or even a ring system. When a ring system of similar characteristics to those of Chariklo and Chiron is included, we can fit both the light curve amplitude and absolute magnitude. In this case the derived axial ratio is modified to $b/a=0.37\pm0.10$. The implied density is $678^{+209}_{-100}$ kg$\,$m$^{-3}$. Also the existence of a ring is consistent with the spectroscopic detection of water ice on Bienor. Nevertheless the other explanations cannot be discarded.

A stellar occultation by a trans-Neptunian object (TNO) provides an opportunity to probe its size and shape. Very few occultations by TNOs have been sampled simultaneously from multiple locations, while a robust estimation of shadow size has been possible for only two objects. We present the first observation of an occultation by the TNO 2007 UK126 on 15 November 2014, measured by three observers, one nearly on and two almost symmetrical to the shadow's centerline. This is the first multi-chord dataset obtained for a so-called detached object, a TNO subgroup with perihelion distances so large that the giant planets have likely not perturbed their orbits. We revisit Herschel/PACS far-infrared data, applying a new reduction method to improve the accuracy of the measured fluxes. Combining both datasets allows us to comprehensively characterize 2007 UK126. We use error-in-variable regression to solve the non-linear problem of propagating timing errors into uncertainties of the ellipse parameters. Based on the shadow's size and a previously reported rotation period, we expect a shape of a Maclaurin spheroid and derive a geometrically plausible size range. To refine our size estimate of 2007 UK126, we model its thermal emission using a thermophysical model code. We conduct a parametric study to predict far-infrared fluxes and compare them to the Herschel/PACS measurements. The favorable geometry of our occultation chords, combined with minimal dead-time imaging, and precise GPS time measurements, allow for an accurate estimation of the shadow size (best-fitting ellipse with axes 645.80 $\pm$ 5.68 km $\times$ 597.81 $\pm$ 12.74 km) and the visual geometric albedo (15.0 $\pm$ 1.6 %). By combining our analyses of the occultation and the far-infrared data, we can constrain the effective diameter of 2007 UK126 to 599 - 629 km. We conclude that subsolar surface temperatures are $\approx$ 50 - 55 K.

At a distance of 1.295 parsecs, the red-dwarf Proxima Centauri ($\alpha$ Centauri C, GL 551, HIP 70890, or simply Proxima) is the Sun's closest stellar neighbor and one of the best studied low-mass stars. It has an effective temperature of only $\sim$ 3050 K, a luminosity of $\sim$0.1 per cent solar, a measured radius of 0.14 R$_\odot$ and a mass of about 12 per cent the mass of the Sun. Although Proxima is considered a moderately active star, its rotation period is $\sim$ 83 days, and its quiescent activity levels and X-ray luminosity are comparable to the Sun's. New observations reveal the presence of a small planet orbiting Proxima with a minimum mass of 1.3~Earth masses and an orbital period of $\sim$11.2 days. Its orbital semi-major axis is $\sim0.05$ AU, with an equilibrium temperature in the range where water could be liquid on its surface.

We present results derived from the first multi-chord stellar occultation by the trans-Neptunian object (229762) 2007 UK$_{126}$, observed on 2014 November 15. The event was observed by the Research and Education Collaborative Occultation Network (RECON) project and International Occultation Timing Association (IOTA) collaborators throughout the United States. Use of two different data analysis methods obtain a satisfactory fit to seven chords, yelding an elliptical fit to the chords with an equatorial radius of $R=338_{-10} ^{+15}$ km and equivalent radius of $R_{eq}=319_{-7} ^{+14}$ km. A circular fit also gives a radius of $R=324_{-23} ^{+30}$ km. Assuming that the object is a Maclaurin spheroid with indeterminate aspect angle, and using two published absolute magnitudes for the body, we derive possible ranges for geometric albedo between $p_{V}=0.159_{-0.013} ^{+0.007}$ and $p_{R}=0.189_{-0.015}^{+0.009}$, and for the body oblateness between $\epsilon=0.105_{-0.040} ^{+0.050}$ and $\epsilon=0.118_{-0.048} ^{+0.055}$. For a nominal rotational period of 11.05 h, an upper limit for density of $\rho=1740$ kg~m$^{-3}$ is estimated for the body.

Haumea is one of the most interesting and intriguing transneptunian objects (TNOs). It is a large, bright, fast rotator, and its spectrum indicates nearly pure water ice on the surface. It has at least two satellites and a dynamically related family of more than ten TNOs with very similar proper orbital parameters and similar surface properties. The Haumean family is the only one currently known in the transneptunian belt. Various models have been proposed but the formation of the family remains poorly understood. In this work, we have investigated the rotational properties of the family members and unconfirmed family candidates with short-term variability studies, and report the most complete review to date. We present results based on five years of observations and report the short-term variability of five family members, and seven candidates. The mean rotational periods, from Maxwellian fits to the frequency distributions, are 6.27+/-1.19 h for the confirmed family members, 6.44+/-1.16 h for the candidates, and 7.65+/-0.54 h for other TNOs (without relation to the family). According to our study, there is a suggestion that Haumea family members rotate faster than other TNOs, however, the sample of family member is still too limited for a secure conclusion. We also highlight the fast rotation of 2002 GH32. This object has a 0.36+/-0.02 mag amplitude lightcurve and a rotational period of about 3.98 h. Assuming 2002 GH32 is a triaxial object in hydrostatic equilibrium, we derive a lower limit to the density of 2.56 g cm^-3. This density is similar to Haumea's and much more dense than other small TNO densities.

We present results from a multi-chord Pluto stellar occultation observed on 29 June 2015 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$\pm$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~$\mu$bar for the surface pressure.

Context: Accurate measurements of diameters of trans-Neptunian objects are extremely complicated to obtain. Thermal modeling can provide good results, but accurate absolute magnitudes are needed to constrain the thermal models and derive diameters and geometric albedos. The absolute magnitude, Hv, is defined as the magnitude of the object reduced to unit helio- and geocentric distances and a zero solar phase angle and is determined using phase curves. Phase coefficients can also be obtained from phase curves. These are related to surface properties, yet not many are known. Aims: Our objective is to measure accurate V band absolute magnitudes and phase coefficients for a sample of trans-Neptunian objects, many of which have been observed, and modeled, within the 'TNOs are cool' program, one of Herschel Space Observatory key projects. Methods: We observed 56 objects using the V and R filters. These data, along with those available in the literature, were used to obtain phase curves and measure V band absolute magnitudes and phase coefficients by assuming a linear trend of the phase curves and considering magnitude variability due to rotational light-curve. Results: We obtained 237 new magnitudes for the 56 objects, six of them with no reported previous measurements. Including the data from the literature we report a total of 110 absolute magnitudes with their respective phase coefficients. The average value of Hv is 6.39, bracketed by a minimum of 14.60 and a maximum of -1.12. In the case of the phase coefficients we report 0.10 mag per degree as the median value and a very large dispersion, ranging from -0.88 up tp 1.35 mag per degree.

In this paper we investigate the opportunities provided by the James Webb Space Telescope (JWST) for significant scientific advances in the study of solar system bodies and rings using stellar occultations. The strengths and weaknesses of the stellar occultation technique are evaluated in light of JWST's unique capabilities. We identify several possible JWST occultation events by minor bodies and rings, and evaluate their potential scientific value. These predictions depend critically on accurate a priori knowledge of the orbit of JWST near the Sun-Earth Lagrange-point 2 (L2). We also explore the possibility of serendipitous stellar occultations by very small minor bodies as a by-product of other JWST observing programs. Finally, to optimize the potential scientific return of stellar occultation observations, we identify several characteristics of JWST's orbit and instrumentation that should be taken into account during JWST's development.

We analyze lunar impact flashes recorded by our team during runs in December 2007, 2011, 2013 and 2014. In total, 12 impact flashes with magnitudes ranging between 7.1 and 9.3 in V band were identified. From these, 9 events could be linked to the Geminid stream. Using these observations the ratio of luminous energy emitted in the flashes with respect to the kinetic energy of the impactors for meteoroids of the Geminid stream is estimated. By making use of the known Geminids meteoroid flux on Earth we found this ratio to be 2.1x10$^{-3}$ on average. We compare this luminous efficiency with other estimations derived in the past for other meteoroid streams and also compare it with other estimations that we present here for the first time by making use of crater diameter measurements. We think that the luminous efficiency has to be revised downward, not upward, at least for sporadic impacts. This implies an increase in the influx of kilogram-sized and larger bodies on Earth that has been derived thus far through the lunar impact flash monitoring technique.

From two observing runs during the 2014 summer at Calar Alto Observatory in Almer\'ia (Spain) and at Sierra Nevada Observatory in Granada (Spain), we were able to derive CCD photometry of the Trans-Neptunian Object 2008 OG$_{19}$. We analyzed the time series and obtained a double-peaked light curve with a peak to valley amplitude of (0.437 $\pm$ 0.011) mag and a rotational period of (8.727$\pm$ 0.003) h. This implies that this object is very elongated, closely resembling Varuna's case. The photometry also allowed us to obtain an absolute magnitude in R-band of (4.39 $\pm$ 0.07) mag. From this result we estimated an equivalent diameter of 2008 OG$_{19}$ which is 619$^{+56}_{-113}$ km using an average albedo for Scattered Disk Objects. Finally we interpreted the results under the assumption of hydrostatic equilibrium and found a lower limit for the density of 544$^{+42}_{-4}$ kg$\,$m$^{-3}$. However, a more likely density is (609 $\pm$ 4) kg$\,$m$^{-3}$ using an aspect angle of 60$^\circ$, which corresponds to the most likely configuration for the spin axis with respect to the observer assuming random orientations.

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.

Since 2009 we are running a project to identify flashes produced by the impact of meteoroids on the surface of the Moon. For this purpose we are employing small telescopes and high-sensitivity CCD video cameras. To automatically identify these events a software package called MIDAS was developed and tested. This package can also perform the photometric analysis of these flashes and estimate the value of the luminous efficiency. Besides, we have implemented in MIDAS a new method to establish which is the likely source of the meteoroids (known meteoroid stream or sporadic background). The main features of this computer program are analyzed here, and some examples of lunar impact events are presented.

We present the results of our Moon impact flashes detection campaigns performed around the maximum activity period of the Perseid meteor shower in 2012 and 2013. Just one flash produced by a Perseid meteoroid was detected in 2012 because of very unfavourable geometric conditions, but 12 of these were confirmed in 2013. The visual magnitude of the flashes ranged between 6.6 and 9.3. A luminous efficiency of 1.8 $\times$ 10$^{-3}$ has been estimated for meteoroids from this stream. According to this value, impactor masses would range between 1.9 and 190 g. In addition, we propose a criterion to establish, from a statistical point of view, the likely origin of impact flashes recorded on the lunar surface.

We propose that several short duration events observed in past stellar occultations by Chiron were produced by rings material. From a reanalysis of the stellar occultation data in the literature we determined two possible orientations of the pole of Chiron's rings, with ecliptic coordinates l=(352+/-10) deg, b=(37+/-10) deg or l=(144+/-10) deg, b=(24+/-10) deg . The mean radius of the rings is (324 +/- 10) km. One can use the rotational lightcurve amplitude of Chiron at different epochs to distinguish between the two solutions for the pole. Both imply lower lightcurve amplitude in 2013 than in 1988, when the rotational lightcurve was first determined. We derived Chiron's rotational lightcurve in 2013 from observations at the 1.23-m CAHA telescope and indeed its amplitude is smaller than in 1988. We also present a rotational lightcurve in 2000 from images taken at CASLEO 2.15-m telescope that is consistent with our predictions. Out of the two poles the l=(144+/-10) deg, b=(24+/-10) deg solution provides a better match to a compilation of rotational lightcurve amplitudes from the literature and those presented here. We also show that using this preferred pole, Chiron's long term brightness variations are compatible with a simple model that incorporates the changing brightness of the rings as the tilt angle with respect to the Earth changes with time. Also, the variability of the water ice band in Chiron's spectra in the literature can be explained to a large degree by an icy ring system whose tilt angle changes with time and whose composition includes water ice, analogously to the case of Chariklo. We present several possible formation scenarios for the rings from qualitative points of view and speculate on the reasons why rings might be common in centaurs. We speculate on whether the known bimodal color distribution of centaurs could be due to presence of rings and lack of them.

We observed comet C/2012 S1 (ISON) during six nights in February 2013 when it was at 4.8 AU from the sun. At this distance and time the comet was not very active and it was theoretically possible to detect photometric variations likely due to the rotation of the cometary nucleus. The goal of this work is to obtain differential photometry of the comet inner coma using different aperture radii in order to derive a possible rotational period. Large field of view images were obtained with a 4k x 4k CCD at the f/3 0.77m telescope of La Hita Observatory in Spain. Aperture photometry was performed in order to get relative magnitude variation versus time. Using calibrated star fields we also obtained ISON's R-magnitudes versus time. We applied a Lomb-Scargle periodogram analysis to get possible periodicities for the observed brightness variations, directly related with the rotation of the cometary nucleus. The comet light curve obtained is very shallow, with a peak-to-peak amplitude of 0.03 $\pm$ 0.02 mag. A tentative synodic rotational period (single-peaked) of 14.4 $\pm$ 1.2 hours for ISON's nucleus is obtained from our analysis, but there are other possibilities. We studied the possible effect of the seeing variations in the obtained periodicities during the same night, and from night to night. These seeing variations had no effect on the derived periodicity. We discuss and interpret all possible solutions for the rotational period of ISON's nucleus.

Until now, rings have been detected in the Solar System exclusively around the four giant planets. Here we report the discovery of the first minor-body ring system around the Centaur object (10199) Chariklo, a body with equivalent radius 124$\pm$9 km. A multi-chord stellar occultation revealed the presence of two dense rings around Chariklo, with widths of about 7 km and 3 km, optical depths 0.4 and 0.06, and orbital radii 391 and 405 km, respectively. The present orientation of the ring is consistent with an edge-on geometry in 2008, thus providing a simple explanation for the dimming of Chariklo's system between 1997 and 2008, and for the gradual disappearance of ice and other absorption features in its spectrum over the same period. This implies that the rings are partially composed of water ice. These rings may be the remnants of a debris disk, which were possibly confined by embedded kilometre-sized satellites.

We have analyzed the meteor activity associated with meteoroids of fresh dust trails of Comet 209P/LINEAR, which produced an outburst of the Camelopardalid meteor shower (IAU code #451, CAM) in May 2014. With this aim, we have employed an array of high-sensitivity CCD video devices and spectrographs deployed at 10 meteor observing stations in Spain in the framework of the Spanish Meteor Network (SPMN). Additional meteoroid flux data were obtained by means of two forward-scatter radio systems. The observed peak zenithal hourly rate (ZHR) was much lower than expected, of around 20 meteors h-1. Despite of the small meteor flux in the optical range, we have obtained precise atmospheric trajectory, radiant and orbital information for 11 meteor and fireball events associated with this stream. The ablation behaviour and low tensile strength calculated for these particles reveal that Camelopardalid meteoroids are very fragile, mostly pristine aggregates with strength similar to that of the Orionids and the Leonids. The mineral grains seem to be glued together by a volatile phase. We also present and discuss two unique emission spectra produced by two Camelopardalid bright meteors. These suggest a non-chondritic nature for these particles, which exhibit Fe depletion in their composition.

In this work we aim to study if the variability in the absolute magnitude of Chariklo and the temporal variation of the spectral ice feature, even its disappearance in 2007, can be explained by an icy ring system whose aspect angle changes with time. We modeled the light reflected by a system as the one described above to explain the variations on the absolute magnitude of Chariklo and its rings. Using X-Shooter at VLT we obtained a new reflectance spectra, here we compared this new set of data with the ones available in the literature. We showed how the water ice feature is visible in 2013 in accordance with the ring configuration, which had an opening angle of nearly 34$^o$ in 2013. Finally we also used models of the scattering of light to fit the visible and near-infrared spectra showing different characteristic to obtain information on the composition of Chariklo and its rings. {We showed that past absolute photometry of Chariklo from the literature and new photometric data that we obtained in 2013 can be explained by a ring of particles whose opening angle changes as a function of time. We used the two possible pole solutions for the ring system and found that only one of them, $\alpha$=151.30$\pm0.5$, $\delta=41.48\pm0.2$ $^o$ ($\lambda=137.9\pm0.5$, $\beta=27.7\pm0.2$ $^o$) provides the right variation of the aspect angle with time to explain the photometry, whereas the other possible pole solution fails to explain the photometry. From spectral modeling, using the result on the pole solution, we derived the composition of Chariklo surface and of that of the rings. Chariklo surface is composed by nearly 60% of amorphous carbon, 30% of silicates and 10\% of organics, no water ice was found on the surface. Whereas the ring contains 20% of water ice, 40-70% of silicates and 10-30% of tholins and small quantities of amorphous carbon.

We present results for the short-term variability of Binary Trans-Neptunian Objects (BTNOs). We performed CCD photometric observations using the 3.58 m Telescopio Nazionale Galileo, the 1.5 m Sierra Nevada Observatory telescope, and the 1.23 m Centro Astronomico Hispano Aleman telescope at Calar Alto Observatory. We present results based on five years of observations and report the short-term variability of six BTNOs. Our sample contains three classical objects: 2003MW12, or Varda, 2004SB60, or Salacia, and 2002 VT130; one detached disk object: 2007UK126; and two resonant objects: 2007TY430 and 2000EB173, or Huya. For each target, possible rotational periods and/or photometric amplitudes are reported. We also derived some physical properties from their lightcurves, such as density, primary and secondary sizes, and albedo. We compiled and analyzed a vast lightcurve database for Trans-Neptunian Objects (TNOs) including centaurs to determine the lightcurve amplitude and spin frequency distributions for the binary and non-binary populations. The mean rotational periods, from the Maxwellian fits to the frequency distributions, are 8.63+/-0.52 h for the entire sample, 8.37+/-0.58 h for the sample without the binary population, and 10.11+/-1.19 h for the binary population alone. Because the centaurs are collisionally more evolved, their rotational periods might not be so primordial. We computed a mean rotational period, from the Maxwellian fit, of 8.86+/-0.58 h for the sample without the centaur population, and of 8.64+/-0.67 h considering a sample without the binary and the centaur populations. According to this analysis, regular TNOs spin faster than binaries, which is compatible with the tidal interaction of the binaries. Finally, we examined possible formation models for several systems studied in this work and by our team in previous papers.

We analyze here two very bright fireballs produced by the ablation in the atmosphere of two large meteoroids in 2009 and 2010. These slow-moving and deep-penetrating events were observed over Spain in the framework of our Spanish Fireball Network continuous meteor monitoring campaign. The analysis of the emission spectrum imaged for one of these fireballs has provided the first clues about the chemical nature of the progenitor meteoroids. The orbital parameters of these particles suggests a likely association with the recently identified July rho-Herculid meteoroid stream. In addition, considerations about the likely parent body of this stream are also made on the basis of orbital dissimilarity criteria. This orbital analysis reveals that both meteoroids and PHA 2007LQ19 exhibit a similar evolution during a time period of almost 8,000 years, which suggests that either this NEO is the potential parent of these particles or that this NEO and both meteoroids had a common progenitor in the past.

The classical Kuiper belt contains objects both from a low-inclination, presumably primordial, distribution and from a high-inclination dynamically excited population. Based on a sample of classical TNOs with observations at thermal wavelengths we determine radiometric sizes, geometric albedos and thermal beaming factors as well as study sample properties of dynamically hot and cold classicals. Observations near the thermal peak of TNOs using infra-red space telescopes are combined with optical magnitudes using the radiometric technique with near-Earth asteroid thermal model (NEATM). We have determined three-band flux densities from Herschel/PACS observations at 70.0, 100.0 and 160.0 $\mu$m and Spitzer/MIPS at 23.68 and 71.42 $\mu$m when available. We have analysed 18 classical TNOs with previously unpublished data and re-analysed previously published targets with updated data reduction to determine their sizes and geometric albedos as well as beaming factors when data quality allows. We have combined these samples with classical TNOs with radiometric results in the literature for the analysis of sample properties of a total of 44 objects. We find a median geometric albedo for cold classical TNOs of 0.14 and for dynamically hot classical TNOs, excluding the Haumea family and dwarf planets, 0.085. We have determined the bulk densities of Borasisi-Pabu (2.1 g/cm^3), Varda-Ilmare (1.25 g/cm^3) and 2001 QC298 (1.14 g/cm^3) as well as updated previous density estimates of four targets. We have determined the slope parameter of the debiased cumulative size distribution of dynamically hot classical TNOs as q=2.3 +- 0.1 in the diameter range 100<D<500 km. For dynamically cold classical TNOs we determine q=5.1 +- 1.1 in the diameter range 160<D<280 km as the cold classical TNOs have a smaller maximum size.

Context. Ionization plays a central role in the gas-phase chemistry of molecular clouds. Since ions are coupled with magnetic fields, which can in turn counteract gravitational collapse, it is of paramount importance to measure their abundance in star-forming regions. Aims. We use spectral line observations of the high-mass star-forming region NGC 6334 I to derive the abundance of two of the most abundant molecular ions, HCO+ and N2H+, and consequently, the cosmic ray ionization rate. In addition, the line profiles provide information about the kinematics of this region. Methods. We present high-resolution spectral line observations conducted with the HIFI instrument on board the Herschel Space Observatory of the rotational transitions with Jup > 5 of the molecular species C17O, C18O, HCO+, H13CO+, and N2H+. Results. The HCO+ and N2H+ line profiles display a redshifted asymmetry consistent with a region of expanding gas. We identify two emission components in the spectra, each with a different excitation, associated with the envelope of NGC 6334 I. The physical parameters obtained for the envelope are in agreement with previous models of the radial structure of NGC 6334 I based on submillimeter continuum observations. Based on our new Herschel/HIFI observations, combined with the predictions from a chemical model, we derive a cosmic ray ionization rate that is an order of magnitude higher than the canonical value of 10^(-17) s-1. Conclusions. We find evidence of an expansion of the envelope surrounding the hot core of NGC 6334 I, which is mainly driven by thermal pressure from the hot ionized gas in the region. The ionization rate seems to be dominated by cosmic rays originating from outside the source, although X-ray emission from the NGC 6334 I core could contribute to the ionization in the inner part of the envelope.

Centaurs are the transitional population between trans-Neptunian objects (TNOs) and Jupiter-family comets. For this reason it is possible to access the smaller ones, which is more difficult to do with the TNO population. The goal of this work is to characterize a set of 16 Centaurs in terms of their size, albedo, and thermal properties. We study the correlations, for a more extended sample obtained from the literature, of diameter, albedo, orbital parameters, and spectral slopes. We performed three-band photometric observations using Herschel-PACS and used a consistent method for the data reduction and aperture photometry of this sample to obtain monochromatic flux densities at 70, 100, and 160 $\mu$m. Additionally, we used Spitzer-MIPS flux densities at 24 and 70 $\mu$m when available. We also included in our Centaur sample scattered disk objects (SDOs), a dynamical family of TNOS, using results previously published by our team, and some Centaurs observed only with the Spitzer/MIPS instrument. We have determined new radiometric sizes and albedos of 16 Centaurs. The first conclusion is that the albedos of Centaur objects are not correlated with their orbital parameters. Similarly, there is no correlation between diameter and orbital parameters. Most of the objects in our sample are dark (pv < 7%) and most of them are small (D < 120km). However, there is no correlation between albedo and diameter, in particular for the group of the small objects as albedo values are homogeneously distributed between 4 to 16%. The correlation with the color of the objects showed that red objects are all small (mean diameter 65 km), while the gray ones span a wide range of sizes (mean diameter 120 km). Moreover, the gray objects tend to be darker, with a mean albedo of 5.6%, compared with a mean of 8.5% (ranging from 5 to 15%) for the red objects.

On 2013 September 11 at 20h07m28.68 $\pm$ 0.01 s UTC, two telescopes operated in the framework of our lunar impact flashes monitoring project recorded an extraordinary flash produced by the impact on the Moon of a large meteoroid at selenographic coordinates 17.2 $\pm$ 0.2 $^{\circ}$ S, 20.5 $\pm$ 0.2 $^{\circ}$ W. The peak brightness of this flash reached 2.9 $\pm$ 0.2 mag in V and it lasted over 8 seconds. The estimated energy released during the impact of the meteoroid was 15.6 $\pm$ 2.5 tons of TNT under the assumption of a luminous efficiency of 0.002. This event, which is the longest and brightest confirmed impact flash recorded on the Moon thus far, is analyzed here. The likely origin of the impactor is discussed. Considerations in relation to the impact flux on Earth are also made.

We analyze the likely link between the recently discovered Near Earth Object 2012XJ112 and a bright fireball observed over the south of Spain on December 27, 2012. The bolide, with an absolute magnitude of -9 +- 1, was simultaneously imaged during the morning twilight from two meteor stations operated by the SPanish Meteor Network (SPMN). It was also observed by several casual witnesses. The emission spectrum produced during the ablation of the meteoroid in the atmosphere was also recorded. From its analysis the chemical nature of this particle was inferred. Although our orbital association software identified several potential parent bodies for this meteoroid, the analysis of the evolution of the orbital elements performed with the Mercury 6 symplectic integrator supports the idea that NEO 2012XJ112 is the source of this meteoroid. The implications of this potential association are discussed here. In particular, the meteoroid bulk chemistry is consistent with a basaltic achondrite, and this emphasizes the importance to deduce from future Earth approaches the reflectance spectrum and taxonomic nature of 2012XJ112.

It is well known that some Main Belt asteroids show comet-like features. A representative example is the first known Main Belt comet 133P/(7968) Elst-Pizarro. If the mechanisms causing this activity are too weak to develop visually evident comae or tails, the objects stay unnoticed. We are presenting a novel way to search for active asteroids, based on looking for objects with deviations from their expected brightnesses in a database. Just by using the MPCAT-OBS Observation Archive we have found five new candidate objects that possibly show a type of comet-like activity, and the already known Main Belt comet 133P/(7968) Elst-Pizarro. Four of the new candidates, (315) Constantia, (1026) Ingrid, (3646) Aduatiques, and (24684) 1990 EU4, show brightness deviations independent of the object's heliocentric distance, while (35101) 1991 PL16 shows deviations dependent on its heliocentric distance, which could be an indication of a thermal triggered mechanism. The method could be implemented in future sky survey programmes to detect outbursts on Main Belt objects almost simultaneously with their occurrence.

Maser lines from different molecular species, including water, hydroxyl, and methanol, are common observational phenomena associated with massive star forming regions. In particular, the methanol maser appears as an ideal tool to study the early phases of massive star formation. However, it is difficult to establish the exact start of the methanol maser phase, and it would then be interesting to detect and study low-flux density methanol masers (i.e., < 0.1 Jy or even << 0.1 Jy), in order to determine if they can effectively be used to mark a specific evolutionary phase in high-mass star formation. Past surveys have been unable to systematically detect many low-flux density methanol masers, and thus we do not yet know how many such masers exist in the Galaxy and what is their physical nature. Out of a sample of 107 observed Hi-GAL sources we detected a total of 32 methanol masers, with 22 sources being new and weak (median peak flux density 0.07 Jy) detections, in the Galactic longitude range [32.0, 59.8]deg. We also detected 12 6.035-GHz OH maser, with 9 objects being new detections. Our survey covers a similar range of source distances as the "Arecibo Methanol Maser Galactic Plane Survey", but the methanol masers detected by us are clearly shifted towards lower integrated flux densities. The newly detected methanol masers are mostly of low-luminosity and, except for some sources, their weakness is not due to distance effects or positional offsets. No specific correlation is found with the physical parameters of the Hi-GAL clumps, except for sources with both CH3OH and OH masers which tend to have higher mass and luminosity. The intensity of the methanol masers correlates well with the velocity range of the maser emission, which suggests that the low brightness of these masers is related to the number of maser spots in the emitting region and their evolution with time.

A superbolide with an estimated absolute magnitude of -20+-1 was seen on July 13, 2012 over the center and south of Spain. This extraordinary event, which was witnessed by numerous casual observers, was recorded in the framework of the continuous fireball monitoring and meteor spectroscopy campaign performed by the SPanish Meteor Network (SPMN). Thus, because of optimal weather conditions, the bolide was imaged from ten meteor observing stations. Here we present the analysis of this magnificent event, which is the brightest fireball ever recorded by our team. The atmospheric trajectory of the bolide and the orbit in the Solar System of the parent meteoroid were obtained. The emission spectrum produced during the ablation of this particle is also discussed. We found that the meteoroid, which was following a Halley Type Comet orbit, was depleted in Na and had a tensile strength one order of magnitude higher than that corresponding to typical cometary materials. By means of orbital analysis tools we have investigated the likely parent body of this particle and the results suggest that the progenitor is a damocloid. The impact area of the hypothetical remnants of the meteoroid is also given and a search for meteorites was performed, but none were found.

A Geminid fireball with an absolute magnitude of -13 was observed over the south of Spain on Dec. 15, 2009. This extraordinarily bright event (the brightest Geminid ever recorded by our team) was imaged from two meteor observing stations operated by the SPanish Meteor Network (SPMN). The bolide exhibited fast and quasi-periodic variations in brightness, a behaviour typically associated to the rotation of the parent meteoroid. The inferred tensile strength of this particle was found to be significantly higher that the typical values obtained for Geminid meteoroids. The fireball penetrated in the atmosphere till a final height of about 25 km above the ground level and a non-zero terminal mass was calculated at the ending point of the luminous trajectory. In this way, the observational evidence points to the existence of a population of meteoroids in the higher end of the Geminid mass distribution capable of producing meteorites. From the relative chemical abundances inferred from the emission spectrum of this bolide we conclude that the Geminid-forming materials are similar to some primitive carbonaceous chondrite groups. Then, we conclude that in meteorite collections from cold deserts, capable to preserve meteorites of few tens of grams, some rare groups of carbonaceous chondrites could be coming from the Geminid parent body: (3200) Phaeton.

Combining stellar occultation observations probing Pluto's atmosphere from 1988 to 2013 and models of energy balance between Pluto's surface and atmosphere, we conclude that Pluto's atmosphere does not collapse at any point in its 248-year orbit. The occultation results show an increasing atmospheric pressure with time in the current epoch, a trend present only in models with a high thermal inertia and a permanent N2 ice cap at Pluto's north rotational pole.

We present the analysis of the atmospheric trajectory and orbital data of four bright bolides observed over Spain, one of which is a potential meteorite dropping event. Their absolute magnitude ranges from -10 to -11. Two of these are of sporadic origin, although a Geminid and a kappa-Cygnid fireball are also considered. These events were recorded in the framework of the continuous fireball monitoring and spectroscopy campaigns developed by the SPanish Meteor Network (SPMN) between 2010 and 2012. The tensile strength of the parent meteoroids is estimated and the abundances of the main rock-forming elements in these particles are calculated from the emission spectrum obtained for three of these events. This analysis revealed a chondritic nature for these meteoroids.

Near-Earth asteroid 2012 DA14 made its closest approach on February 15, 2013, when it passed at a distance of 27,700 km from the Earth's surface. It was the first time an asteroid of moderate size was predicted to approach that close to the Earth, becoming bright enough to permit a detailed study from ground-based telescopes. Asteroid 2012 DA14 was poorly characterized before its closest approach. We acquired data using several telescopes on four Spanish observatories: the 10.4m Gran Telescopio Canarias (GTC) and the 3.6m Telescopio Nazionale Galileo (TNG), both in the El Roque de los Muchachos Observatory (ORM, La Palma); the 2.2m CAHA telescope, in the Calar Alto Observatory (Almeria); the f/3 0.77m telescope in the La Hita Observatory (Toledo); and the f/8 1.5m telescope in the Sierra Nevada Observatory (OSN, Granada). We obtained visible and near-infrared color photometry, visible spectra and time-series photometry. Visible spectra together with color photometry of 2012 DA14 show that it can be classified as an L-type asteroid, a rare spectral type with a composition similar to that of carbonaceous chondrites. The time-series photometry provides a rotational period of 8.95 +- 0.08 hours after the closest approach, and there are indications that the object suffered a spin-up during this event. The large amplitude of the light curve suggests that the object is very elongated and irregular, with an equivalent diameter of around 18m. We obtain an absolute magnitude of H_R = 24.5 +- 0.2, corresponding to H_V = 25.0 +- 0.2. The GTC photometry also gives H_V = 25.29 +- 0.14. Both values agree with the value listed at the Minor Planet Center shortly after discovery. From the absolute photometry, together with some constraints on size and shape, we compute a geometric albedo of p_V = 0.44 +- 0.20, which is slightly above the range of albedos known for L-type asteroids (0.082 - 0.405).

On October 8th, 2011 the Earth crossed the dust trails left by comet 21P/Giacobini-Zinner during its XIX and XX century perihelion approaches with the comet being close to perihelion. The geometric circumstances of that encounter were thus favorable to produce a meteor storm, but the trails were much older than in the 1933 and 1946 historical encounters. As a consequence the 2011 October Draconid display exhibited several activity peaks with Zenithal Hourly Rates of about 400 meteors per hour. In fact, if the display had been not forecasted, it could have passed almost unnoticed as was strongly attenuated for visual observers due to the Moon. This suggests that most meteor storms of a similar nature could have passed historically unnoticed under unfavorable weather and Moon observing conditions. The possibility of obtaining information on the physical properties of cometary meteoroids penetrating the atmosphere under low-geocentric velocity encounter circumstances motivated us to set up a special observing campaign. Added to the Spanish Fireball Network wide-field all-sky and CCD video monitoring, other high-sensitivity 1/2" black and white CCD video cameras were attached to modified medium-field lenses for obtaining high resolution orbital information. The trajectory, radiant, and orbital data of 16 October Draconid meteors observed at multiple stations are presented. The results show that the meteors appeared from a geocentric radiant located at R.A.=263.0+-0.4 deg. and Dec.=+55.3+-0.3 deg. that is in close agreement with the radiant predicted for the 1873-1894 and the 1900 dust trails. The estimated mass of material from 21P/Giacobini-Zinner delivered to Earth during the six-hours outburst was around 950+-150 kg.

On October 8, 2011 the Earth crossed dust trails ejected from comet 21P/Giacobini-Zinner in the late 19th and early 20th Century. This gave rise to an outburst in the activity of the October Draconid meteor shower, and an international team was organized to analyze this event. The SPanish Meteor Network (SPMN) joined this initiative and recorded the October Draconids by means of low light level CCD cameras. In addition, spectroscopic observations were carried out. Tens of multi-station meteor trails were recorded, including an extraordinarily bright October Draconid fireball (absolute mag. -10.5) that was simultaneously imaged from three SPMN meteor ob-serving stations located in Andalusia. Its spectrum was obtained, showing a clear evolution in the relative intensity of emission lines as the fireball penetrated deeper into the atmosphere. Here we focus on the analysis of this remarkable spectrum, but also discuss the atmospheric trajectory, atmospheric penetration, and orbital data computed for this bolide which was probably released during 21P/Giacobini-Zinner return to perihelion in 1907. The spectrum is discussed together with the tensile strength for the October Draconid meteoroids. The chemical profile evolution of the main rocky elements for this extremely bright bolide is compared with the elemental abundances obtained for 5 October Draconid fireballs also recorded during our spectroscopic campaign but observed only at a single station. Significant chemical heterogeneity between the small meteoroids is found as we should expect for cometary aggregates being formed by diverse dust components.

2012 DR30 is a recently discovered Solar System object on a unique orbit, with a high eccentricity of 0.9867, a perihelion distance of 14.54 AU and a semi-major axis of 1109 AU, in this respect outscoring the vast majority of trans-Neptunian objects. We performed Herschel/PACS and optical photometry to uncover the size and albedo of 2012 DR30, together with its thermal and surface properties. The body is 185 km in diameter and has a relatively low V-band geometric albedo of ~8%. Although the colours of the object indicate that 2012 DR30 is an RI taxonomy class TNO or Centaur, we detected an absorption feature in the Z-band that is uncommon among these bodies. A dynamical analysis of the target's orbit shows that 2012 DR30 moves on a relatively unstable orbit and was most likely only recently placed on its current orbit from the most distant and still highly unexplored regions of the Solar System. If categorised on dynamical grounds 2012 DR30 is the largest Damocloid and/or high inclination Centaur observed so far.

We present our latest results about the short-term variability of trans-Neptunian objects (TNOs). We performed broad-band CCD photometric observations using several telescopes in Spain and Chile. We present results based on three years of observations and report the short-term variability of 10 TNOs. Our sample of studied targets contains classical objects: (275809) 2001 QY297, (307251) 2002 KW14, (55636) 2002 TX300, 2004 NT33, (230965) 2004 XA192, and (202421) 2005 UQ513, a resonant body: (84522) 2002 TC302, a scattered target: (44594) 1999 OX3, and two detached objects: (145480) 2005 TB190, and (40314) 1999 KR16. For each target, light curves as well as possible rotation periods and photometric amplitudes are reported. The majority of the observed objects present a low peak-to-peak amplitude, <0.15mag. Only two objects exhibit light curve amplitudes higher than 0.15mag: (275809) 2001 QY297, and (307251) 2002 KW14. We remark two biases in the literature, previously studied in Thirouin et al. (2010) and confirmed by this new study: a bias towards objects with a small amplitude light curve and a second one against objects with a long rotational period in the data base of published rotational periods. We derived constraints on physical properties of some targets. We also report the solar phase curves of (40314) 1999 KR16, and (44594) 1999 OX3 for solar phase angles from 0 degree to around 2 degrees. Part of our discussion is focused on the study of (275809) 2001 QY297 which turned out to be an asynchronous binary system.

We present estimates of the basic physical properties (size and albedo) of (90377) Sedna, a prominent member of the detached trans-Neptunian object population and the recently discovered scattered disk object 2010 EK139, based on the recent observations acquired with the Herschel Space Observatory, within the "TNOs are Cool!" key programme. Our modeling of the thermal measurements shows that both objects have larger albedos and smaller sizes than the previous expectations, thus their surfaces might be covered by ices in a significantly larger fraction. The derived diameter of Sedna and 2010 EK139 are 995 +/- 80 km and 470 +35/-10 km, while the respective geometric albedos are pV 0.32 +/- 0.06 and 0.25 +0.02/-0.05. These estimates are based on thermophysical model techniques.

Trans-Neptunian objects (TNO) represent the leftovers of the formation of the Solar System. Their physical properties provide constraints to the models of formation and evolution of the various dynamical classes of objects in the outer Solar System. Based on a sample of 19 classical TNOs we determine radiometric sizes, geometric albedos and beaming parameters. Our sample is composed of both dynamically hot and cold classicals. We study the correlations of diameter and albedo of these two subsamples with each other and with orbital parameters, spectral slopes and colors. We have done three-band photometric observations with Herschel/PACS and we use a consistent method for data reduction and aperture photometry of this sample to obtain monochromatic flux densities at 70.0, 100.0 and 160.0 \mu m. Additionally, we use Spitzer/MIPS flux densities at 23.68 and 71.42 \mu m when available, and we present new Spitzer flux densities of eight targets. We derive diameters and albedos with the near-Earth asteroid thermal model (NEATM). As auxiliary data we use reexamined absolute visual magnitudes from the literature and data bases, part of which have been obtained by ground based programs in support of our Herschel key program. We have determined for the first time radiometric sizes and albedos of eight classical TNOs, and refined previous size and albedo estimates or limits of 11 other classicals. The new size estimates of 2002 MS4 and 120347 Salacia indicate that they are among the 10 largest TNOs known. Our new results confirm the recent findings that there are very diverse albedos among the classical TNOs and that cold classicals possess a high average albedo (0.17 +/- 0.04). Diameters of classical TNOs strongly correlate with orbital inclination in our sample. We also determine the bulk densities of six binary TNOs.

We present Herschel PACS photometry of 18 Plutinos and determine sizes and albedos for these objects using thermal modeling. We analyze our results for correlations, draw conclusions on the Plutino size distribution, and compare to earlier results. Flux densities are derived from PACS mini scan-maps using specialized data reduction and photometry methods. In order to improve the quality of our results, we combine our PACS data with existing Spitzer MIPS data where possible, and refine existing absolute magnitudes for the targets. The physical characterization of our sample is done using a thermal model. Uncertainties of the physical parameters are derived using customized Monte Carlo methods. The correlation analysis is performed using a bootstrap Spearman rank analysis. We find the sizes of our Plutinos to range from 150 to 730 km and geometric albedos to vary between 0.04 and 0.28. The average albedo of the sample is 0.08 \pm 0.03, which is comparable to the mean albedo of Centaurs, Jupiter Family comets and other Trans-Neptunian Objects. We were able to calibrate the Plutino size scale for the first time and find the cumulative Plutino size distribution to be best fit using a cumulative power law with q = 2 at sizes ranging from 120-400 km and q = 3 at larger sizes. We revise the bulk density of 1999 TC36 and find a density of 0.64 (+0.15/-0.11) g cm-3. On the basis of a modified Spearman rank analysis technique our Plutino sample appears to be biased with respect to object size but unbiased with respect to albedo. Furthermore, we find biases based on geometrical aspects and color in our sample. There is qualitative evidence that icy Plutinos have higher albedos than the average of the sample.

Physical characterization of Trans-Neptunian objects, a primitive population of the outer solar system, may provide constraints on their formation and evolution. The goal of this work is to characterize a set of 15 scattered disk (SDOs) and detached objects, in terms of their size, albedo, and thermal properties. Thermal flux measurements obtained with the Herschel-PACS instrument at 70, 100 and 160 \mu m, and whenever applicable, with Spitzer-MIPS at 24 and 70 \mu m, are modeled with radiometric techniques, in order to derive the objects' individual size, albedo and when possible beaming factor. Error bars are obtained from a Monte-Carlo approach. We look for correlations between these and other physical and orbital parameters. Diameters obtained for our sample range from 100 to 2400 km, and the geometric albedos (in V band) vary from 3.8 % to 84.5 %. The unweighted mean V geometric albedo for the whole sample is 11.2 % (excluding Eris); 6.9 % for the SDOs, and 17.0 % for the detached objects (excluding Eris). We obtain new bulk densities for three binary systems: Ceto/Phorcys, Typhon/Echidna and Eris/Dysnomia. Apart from correlations clearly due to observational bias, we find significant correlations between albedo and diameter (more reflective objects being bigger), and between albedo, diameter and perihelion distance (brighter and bigger objects having larger perihelia). We discuss possible explanations for these correlations.

We present a 3-mm and 1.3-cm spectral line survey conducted with the Mopra 22-m and Parkes 64-m radio telescopes of a sample of 40 cold dust cores, previously observed with BLAST, including both starless and proto-stellar sources. 20 objects were also mapped using molecular tracers of dense gas. To trace the dense gas we used the molecular species NH3, N2H+, HNC, HCO+, H13CO+, HCN and H13CN, where some of them trace the more quiescent gas, while others are sensitive to more dynamical processes. The selected cores have a wide variety of morphological types and also show physical and chemical variations, which may be associated to different evolutionary phases. We find evidence of systematic motions in both starless and proto-stellar cores and we detect line wings in many of the proto-stellar cores. Our observations probe linear distances in the sources >~0.1pc, and are thus sensitive mainly to molecular gas in the envelope of the cores. In this region we do find that, for example, the radial profile of the N2H+(1-0) emission falls off more quickly than that of C-bearing molecules such as HNC(1-0), HCO+(1-0) and HCN(1-0). We also analyze the correlation between several physical and chemical parameters and the dynamics of the cores. Depending on the assumptions made to estimate the virial mass, we find that many starless cores have masses below the self-gravitating threshold, whereas most of the proto-stellar cores have masses which are near or above the self-gravitating critical value. An analysis of the median properties of the starless and proto-stellar cores suggests that the transition from the pre- to the proto-stellar phase is relatively fast, leaving the core envelopes with almost unchanged physical parameters.

We present several lines of evidence based on different kinds of observations to conclude that rotational fission has likely occurred for a fraction of the known Trans-Neptunian Objects (TNOs). It is also likely that a number of binary systems have formed from that process in the trans-neptunian belt. We show that Haumea is a potential example of an object that has suffered a rotational fission. Its current fast spin would be a slight evolution of a primordial fast spin, rather than the result of a catastrophic collision, because the percentage of objects rotating faster than 4 hours would not be small in a maxwellian distribution of spin rates that fits the current TNO rotation database. On the other hand, the specific total angular momentum of Haumea and its satellites falls close to that of the high size ratio asteroid binaries, which are thought to be the result of rotational fissions or mass shedding. We also present N-body simulations of rotational fissions applied to the case of Haumea, which show that this process is feasible, might have generated satellites, and might have even created a "family" of bodies orbitally associated to Haumea. The orbitally associated bodies may come from the direct ejection of fragments according to our simulations, or through the evolution of a proto-satellite formed during the fission event. Also, the disruption of an escaped fragment after the fission might create the orbitally related bodies. If any of those mechanisms are correct, other rotational fission families may be detectable in the trans-neptunian belt in the future, and perhaps even TNO pairs might be found (pairs of bodies sharing very similar orbital elements, but not bound together).

We report on the follow-up and recovery of 100 program NEAs, PHAs and VIs using the ESO/MPG 2.2m, Swope 1m and INT 2.5m telescopes equipped with large field cameras. The 127 fields observed during 11 nights covered 29 square degrees. Using these data, we present the incidental survey work which includes 558 known MBAs and 628 unknown moving objects mostly consistent with MBAs from which 58 objects became official discoveries. We planned the runs using six criteria and four servers which focus mostly on faint and poorly observed objects in need of confirmation, follow-up and recovery. We followed 62 faint NEAs within one month after discovery and we recovered 10 faint NEAs having big uncertainties at their second or later opposition. Using the INT we eliminated 4 PHA candidates and VIs. We observed in total 1,286 moving objects and we reported more than 10,000 positions. All data were reduced by the members of our network in a team effort, and reported promptly to the MPC. The positions of the program NEAs were published in 27 MPC and MPEC references and used to improve their orbits. The O-C residuals for known MBAs and program NEAs are smallest for the ESO/MPG and Swope and about four times larger for the INT whose field is more distorted. The incidental survey allowed us to study statistics of the MBA and NEA populations observable today with 1--2m facilities. We calculate preliminary orbits for all unknown objects, classifying them as official discoveries, later identifications and unknown outstanding objects. The orbital elements a, e, i calculated by FIND_ORB software for the official discoveries and later identified objects are very similar with the published elements which take into account longer observational arcs; thus preliminary orbits were used in statistics for the whole unknown dataset. (CONTINUED)

We present Cygnus X in a new multi-wavelength perspective based on an unbiased BLAST survey at 250, 350, and 500 micron, combined with rich datasets for this well-studied region. Our primary goal is to investigate the early stages of high mass star formation. We have detected 184 compact sources in various stages of evolution across all three BLAST bands. From their well-constrained spectral energy distributions, we obtain the physical properties mass, surface density, bolometric luminosity, and dust temperature. Some of the bright sources reaching 40 K contain well-known compact H II regions. We relate these to other sources at earlier stages of evolution via the energetics as deduced from their position in the luminosity-mass (L-M) diagram. The BLAST spectral coverage, near the peak of the spectral energy distribution of the dust, reveals fainter sources too cool (~ 10 K) to be seen by earlier shorter-wavelength surveys like IRAS. We detect thermal emission from infrared dark clouds and investigate the phenomenon of cold ``starless cores" more generally. Spitzer images of these cold sources often show stellar nurseries, but these potential sites for massive star formation are ``starless" in the sense that to date there is no massive protostar in a vigorous accretion phase. We discuss evolution in the context of the L-M diagram. Theory raises some interesting possibilities: some cold massive compact sources might never form a cluster containing massive stars; and clusters with massive stars might not have an identifiable compact cold massive precursor.

From CCD observations of a fixed and large star field that contained the binary TNO Orcus, we have been able to derive high-precision relative astrometry and photometry of the Orcus system with respect to background stars. The RA residuals of an orbital fit to the astrometric data revealed a periodicity of 9.7+-0.3 days, which is what one would expect to be induced by the known Orcus companion. The residuals are also correlated with the theoretical positions of the satellite with regard to the primary. We therefore have revealed the presence of Orcus' satellite in our astrometric measurements. The photocenter motion is much larger than the motion of Orcus around the barycenter, and we show here that detecting some binaries through a carefully devised astrometric technique might be feasible with telescopes of moderate size. We also analyzed the system's mid-term photometry to determine whether the rotation could be tidally locked to the satellite's orbital period. We found that a photometric variability of 9.7+-0.3 days is clear in our data, and is nearly coincident with the orbital period of the satellite. We believe this variability might be induced by the satellite's rotation. There is also a slight hint for an additional small variability in the 10 hr range that was already reported in the literature. This short-term variability would indicate that the primary is not tidally locked and therefore the system would not have reached a double synchronous state. Implications for the basic physical properties of the primary and its satellite are discussed. From angular momentum considerations we suspect that the Orcus satellite might have formed from a rotational fission. This requires that the mass of the satellite would be around 0.09 times that of the primary, close to the value that one derives by using an albedo of 0.12 for the satellite and assuming equal densities for both objects.

The goal of the paper is to present new results on light echoes from young stellar objects. Broad band CCD images were obtained over three months at one-to-two week intervals for the field of NGC 6726, using the large field-of-view remotely-operated telescope on top of Cerro Burek. We detected scattered light echoes around two young, low-amplitude, irregular variable stars. Observations revealed not just one, but multiple light echoes from brightness pulses of the T Tauri star S CrA and the Herbig Ae/Be star R CrA. Analysis of S CrA's recurring echoes suggests that the star is located 138 +/- 16 pc from Earth, making these the closest echoes ever detected. The environment that scatters the stellar light from S CrA is compatible with an incomplete dust shell or an inclined torus some 10,000 AU in radius and containing $\sim$ $2 \times 10^{-3}$ $M_{\sun}$ of dust. The cause of such concentration at $\sim$ 10,000AU from the star is unknown. It could be the remnant of the envelope from which the star formed, but the distance of the cloud is remarkably similar to the nominal distance of the Oort cloud to the Sun, leading us to also speculate that the dust (or ice) seen around S CrA might have the same origin as the Solar System Oort cloud.

The dust ejecta of Main-Belt Comet P/2010 A2 (LINEAR) have been observed with several telescopes at the at the Observatorio del Roque de los Muchachos on La Palma, Spain. Application of an inverse dust tail Monte Carlo method to the images of the dust ejecta from the object indicates that a sustained, likely water-ice driven, activity over some eight months is the mechanism responsible for the formation of the observed tail. The total amount of dust released is estimated to be 5E7 kg, which represents about 0.3% of the nucleus mass. While the event could have been triggered by a collision, this cannot be decided from the currently available data.

A prototype of a low cost Adaptive Optics (AO) system has been developed at the Instituto de Astrofisica de Andalucia (CSIC) and tested at the 2.2m telescope of the Calar Alto observatory. We present here the status of the project, which includes the image stabilization system and compensation of high order wavefront aberrations with a membrane deformable mirror. The image stabilization system consists of magnet driven tip-tilt mirror. The higher order compensation system comprises of a Shack-Hartmann sensor, a membrane deformable mirror with 39 actuators and the control computer that allows operations up to 420Hz in closed loop mode. We have successfully closed the high order AO loop on natural guide stars. An improvement of 4 times in terms of FWHM was achieved. The description and the results obtained on the sky are presented in this paper.

Thermal emission from Kuiper Belt object (136108) Haumea was measured with Herschel-PACS at 100 and 160 micrometers for almost a full rotation period. Observations clearly indicate a 100-micrometer thermal lightcurve with an amplitude of a factor of ~ 2, which is positively correlated with the optical lightcurve. This confirms that both are primarily due to shape effects. A 160-micrometer lightcurve is marginally detected. Radiometric fits of the mean Herschel- and Spitzer- fluxes indicate an equivalent diameter D ~ 1300 km and a geometric albedo p_v ~ 0.70-0.75. These values agree with inferences from the optical lightcurve, supporting the hydrostatic equilibrium hypothesis. The large amplitude of the 100-micrometer lightcurve suggests that the object has a high projected a/b axis ratio (~ 1.3) and a low thermal inertia as well as possible variable infrared beaming. This may point to fine regolith on the surface, with a lunar-type photometric behavior. The quality of the thermal data is not sufficient to clearly detect the effects of a surface dark spot.

The goal of the Herschel Open Time Key programme "TNOs are Cool!" is to derive the physical and thermal properties for a large sample of Centaurs and trans-Neptunian objects (TNOs), including resonant, classical, detached and scattered disk objects. We present results for seven targets either observed in PACS point-source, or in mini scan-map mode. Spitzer-MIPS observations were included for three objects. The sizes of these targets range from 100 km to almost 1000 km, five have low geometric albedos below 10%, (145480) 2005 TB190 has a higher albedo above 15%. Classical thermal models driven by an intermediate beaming factor of $\eta$=1.2 or $\eta$-values adjusted to the observed colour temperature fit the multi-band observations well in most cases. More sophisticated thermophysical models give very similar diameter and albedo values for thermal inertias in the range 0-25 Jm-2s-0.5K-1, consistent with very low heat conductivities at temperatures far away from the Sun. The early experience with observing and model strategies will allow us to derive physical and thermal properties for our complete Herschel TNO sample of 140 targets as a benchmark for understanding the solar system debris disk, and extra-solar ones as well.

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.

We present results of 6 years of observations, reduced and analyzed with the same tools in a systematic way. We report completely new data for 15 objects, for 5 objects we present a new analysis of previously published results plus additional data and for 9 objects we present a new analysis of data already published. Lightcurves, possible rotation periods and photometric amplitudes are reported for all of them. The photometric variability is smaller than previously thought: the mean amplitude of our sample is 0.1mag and only around 15% of our sample has a larger variability than 0.15mag. The smaller variability than previously thought seems to be a bias of previous observations. We find a very weak trend of faster spinning objects towards smaller sizes, which appears to be consistent with the fact that the smaller objects are more collisionally evolved, but could also be a specific feature of the Centaurs, the smallest objects in our sample. We also find that the smaller the objects, the larger their amplitude, which is also consistent with the idea that small objects are more collisionally evolved and thus more deformed. Average rotation rates from our work are 7.5h for the whole sample, 7.6h for the TNOs alone and 7.3h for the Centaurs. All of them appear to be somewhat faster than what one can derive from a compilation of the scientific literature and our own results. Maxwellian fits to the rotation rate distribution give mean values of 7.5h (for the whole sample) and 7.3h (for the TNOs only). Assuming hydrostatic equilibrium we can determine densities from our sample under the additional assumption that the lightcurves are dominated by shape effects, which is likely not realistic. The resulting average density is 0.92g/cm^3 which is not far from the density constraint that one can derive from the apparent spin barrier that we observe.

We report on optical observations of Gamma-Ray Bursts (GRBs) followed up by our collaboration with the 1.23m telescope located at the Calar Alto observatory. The 1.23m telescope is an old facility, currently undergoing upgrades to enable fully autonomous response to GRB alerts. We discuss the current status of the control system upgrade of the 1.23m telescope. The upgrade is being done by the ARAE our group, based on members of IAA (Instituto de Astrofiisica de Andalucia). Currently the ARAE group is responsible to develop the BOOTES network of robotic telescopes based on the Remote Telescope System, 2nd Version (RTS2), which controls the available instruments and interacts with the EPICS database of Calar Alto. Currently the telescope can run fully autonomously or under observer supervision using RTS2. The fast reaction response mode for GRB reaction (typically with response times below 3 minutes from the GRB onset) still needs some development and testing. The telescope is usually operated in legacy interactive mode, with periods of supervised autonomous runs under RTS2. We show the preliminary results of several GRBs followed up with observer intervention during the testing phase of the 1.23m control software upgrade.

We analyze a vast light curve database by obtaining mean rotational properties of the entire sample, determining the spin frequency distribution and comparing those data with a simple model based on hydrostatic equilibrium. For the rotation periods, the mean value obtained is 6.95 h for the whole sample, 6.88 h for the Trans-neptunian objects (TNOs) alone and 6.75 h for the Centaurs. From Maxwellian fits to the rotational frequencies distribution the mean rotation rates are 7.35 h for the entire sample, 7.71 h for the TNOs alone and 8.95 h for the Centaurs. These results are obtained by taking into account the criteria of considering a single-peak light curve for objects with amplitudes lower than 0.15 mag and a double-peak light curve for objects with variability >0.15mag. The best Maxwellian fits were obtained with the threshold between 0.10 and 0.15mag. The mean light-curve amplitude for the entire sample is 0.26 mag, 0.25mag for TNOs only, and 0.26mag for the Centaurs. The amplitude versus Hv correlation clearly indicates that the smaller (and collisionally evolved) objects are more elongated than the bigger ones. From the model results, it appears that hydrostatic equilibrium can explain the statistical results of almost the entire sample, which means hydrostatic equilibrium is probably reached by almost all TNOs in the H range [-1,7]. This implies that for plausible albedos of 0.04 to 0.20, objects with diameters from 300km to even 100km would likely be in equilibrium. Thus, the great majority of objects would qualify as being dwarf planets because they would meet the hydrostatic equilibrium condition. The best model density corresponds to 1100 kg/m3.

Across all disciplines that work with image data - from astrophysics to medical research and historic preservation - there is a growing need for efficient ways to browse and inspect large sets of high-resolution images. We present the development of a visualization software for solar physics data based on the JPEG 2000 image compression standard. Our implementation consists of the JHelioviewer client application that enables users to browse petabyte-scale image archives and the JHelioviewer server, which integrates a JPIP server, metadata catalog and an event server. JPEG 2000 offers many useful new features and has the potential to revolutionize the way high-resolution image data are disseminated and analyzed. This is especially relevant for solar physics, a research field in which upcoming space missions will provide more than a terabyte of image data per day. Providing efficient access to such large data volumes at both high spatial and high time resolution is of paramount importance to support scientific discovery.

Photometric surveys of transNeptunian objects (TNOs) and Centaurs have suggested possible correlations between some orbital parameters and surface colors of classical objects, scattered disk objects (SDOs), and Centaurs. However, larger sample sizes are needed in order to corroborate or rule out the possible correlations and find some possible new ones. We use VLT-FORS images through BVRI filters of 32 Kuiper Belt Objects (KBOs) and obtain their colors after proper reduction and calibration. We study the possible correlations merging these new measurements with the VLT published results from the ESO large program and with the latest published results of the Meudon Multicolor Survey via non-parametric statistical tests. We obtain a large dataset of 116 objects (classical, SDOs and Centaurs) and, in addition to confirming most of the correlations and conclusions reached in the literature, some possible new correlations are found. The most interesting ones are some correlations of color vs. orbital parameters for the different dynamical groups. We find that some correlations in the classical group, as well as the (dynamically) cold and hot subgroups depend on the size of the objects. As a by-product of our study, we were able to identify new candidates for light curve studies and found that ~55% of the objects showed variability above 0.15 mags. This is a higher value than what is found in other studies. Since our sample contains smaller objects than samples from other studies, this result might be an indication that the smaller TNOs are more elongated than the larger ones.

The fluorescence technique has been successfully used to detect ultrahigh energy cosmic rays by indirect measurements. The underlying idea is that the number of charged particles in the atmospheric shower, i.e, its longitudinal profile, can be extracted from the amount of emitted nitrogen fluorescence light. However the influence of shower fluctuations and the very possible presence of different nuclear species in the primary cosmic ray spectrum makes the estimate of the shower energy from the fluorescence data analysis a difficult task. We investigate the potential of shower size at maximum depth as estimator of shower energy. The detection of the fluorescence light is simulated in detail and the reconstruction biases are carefully analyzed. We extend our calculations to both HiRes and EUSO experiments. This approach has shown some advantages to the reconstruction of the energy when compared to the standard analysis procedure.

Several methods can be used to perform statistical inference of primary composition of cosmic rays measured with water Cerenkov detectors as those in use at the Pierre Auger Southern Observatory. In the present work we assess the impact of additional information about the number of muons in the air shower, on the problem of statistical primary mass discrimination. Several tools are studied, including neural networks, principal component analysis and traditional methods in current use in the field. For our case study we use hypothetical plastic scintillators as muon counters, buried at the side and outside the shade of the water Cerenkov tanks. The study is extended to protons and Fe nuclei impinging on an array with two different spacings, 750 and 1500 m and, therefore, suitable to the 1-10 EeV energy range. A prototype of such a detector is under construction.

Along the years HiRes and AGASA experiments have explored the fluorescence and the ground array experimental techniques to measure extensive air showers, being both essential to investigate the ultra-high energy cosmic rays. However, such Collaborations have published contradictory energy spectra for energies above the GZK cut-off. In this article, we investigate the acceptance of fluorescence telescopes to different primary particles at the highest energies. Using CORSIKA and CONEX shower simulations without and with the new pre-showering scheme, which allows photons to interact in the Earth magnetic field, we estimate the aperture of the HiRes-I telescope for gammas, iron nuclei and protons primaries as a function of the number of simulated events and primary energy. We also investigate the possibility that systematic differences in shower development for hadrons and gammas could mask or distort vital features of the cosmic ray energy spectrum at energies above the photo-pion production threshold. The impact of these effects on the true acceptance of a fluorescence detector is analyzed in the context of top-down production models.

The HiRes and AGASA Collaborations have published contradictory energy spectra for energies above the GZK cut-off. In this article, we investigate the acceptance of fluorescence telescopes to different primary particles at the highest energies. Using CORSIKA shower simulations without and with the new pre-showering scheme, which allows photons to interact in the Earth magnetic field, we estimate the aperture of the HiRes I telescope for gammas and protons primaries. We calculate the dependence of the telescope sensitivity to primary particle identity. We also investigate the possibility that systematic differences in shower development for hadrons and gammas could mask or distort vital features of the cosmic ray energy spectrum at energies above the photo-pion production threshold. The impact of these effects on the true acceptance of a fluorescence detector is analyzed in the context of top-down production models.

The ultrahigh energy tail of the cosmic ray spectrum has been explored with unprecedented detail. For this reason, new experiments are exerting a severe pressure on extensive air shower modeling. Detailed fast codes are in need in order to extract and understand the richness of information now available. In this sense we explore the potential of SENECA, an efficient hybrid tridimensional simulation code, as a valid practical alternative to full Monte Carlo simulations of extensive air showers generated by ultrahigh energy cosmic rays. We discuss the influence of this approach on the main longitudinal characteristics of proton, iron nucleus and gamma induced air showers for different hadronic interaction models. We also show the comparisons of our predictions with those of CORSIKA code.

The fluorescence technique has been successfully used to detect ultrahigh energy cosmic rays by indirect measurements. The underlying idea is that the number of charged particles in the atmospheric shower, i.e, its longitudinal profile, can be extracted from the amount of emitted nitrogen fluorescence light. However the influence of shower fluctuations and the very possible presence of different nuclear species in the primary cosmic ray spectrum make the estimate of the shower energy from the fluorescence data analysis a difficult task. We investigate the potential of shower size at maximum depth as estimator of shower energy. The detection of the fluorescence light is simulated in detail and the reconstruction biases are carefully analyzed. We extend our calculations to both Auger and EUSO experiments. This kind of approach is of particular interest for showers that are not fully contained inside the field of view of the detector.

The impact of the lateral distribution of light in extensive air showers on the detection and reconstruction of shower profiles is investigated for the Auger fluorescence telescopes. Based on three-dimensional simulations, the capability of the Auger telescopes to measure the lateral distribution of light is evaluated. The ability to infer the actual lateral distribution is confirmed by the comparison of detailed simulations with real data. The contribution of pixels located far from the axis of the shower image is calculated and the accepted signal is rescaled in order to reconstruct a correct shower profile. The analysis presented here shows that: (a) the Auger telescopes are able to observe the lateral distribution of showers and (b) the energy corrections to account for the signal in outlying pixels can exceed 10%, depending on shower geometry.

New experiments, exploring the ultra-high energy tail of the cosmic ray spectrum with unprecedented detail, are exerting a severe pressure on extensive air hower modeling. Detailed fast codes are in need in order to extract and understand the richness of information now available. Some hybrid simulation codes have been proposed recently to this effect (e.g., the combination of the traditional Monte Carlo scheme and system of cascade equations or pre-simulated air showers). In this context, we explore the potential of SENECA, an efficient hybrid tridimensional simulation code, as a valid practical alternative to full Monte Carlo simulations of extensive air showers generated by ultra-high energy cosmic rays. We extensively compare hybrid method with the traditional, but time consuming, full Monte Carlo code CORSIKA which is the de facto standard in the field. The hybrid scheme of the SENECA code is based on the simulation of each particle with the traditional Monte Carlo method at two steps of the shower development: the first step predicts the large fluctuations in the very first particle interactions at high energies while the second step provides a well detailed lateral distribution simulation of the final stages of the air shower. Both Monte Carlo simulation steps are connected by a cascade equation system which reproduces correctly the hadronic and electromagnetic longitudinal profile. We study the influence of this approach on the main longitudinal characteristics of proton-induced air showers and compare the predictions of the well known CORSIKA code using the QGSJET hadronic interaction model.

The air fluorescence technique has long been used to detect extensive air showers and to reconstruct its geometry and energy. The fluorescence photon yield of an electron in air is of main importance in the reconstruction procedure. Historically, the fluorescence yield used in the reconstruction of the showers is approximated at all energies by that of an electron with kinetic energy of 80 MeV, because these are the most abundant in a shower. In this paper, we calculate the fluorescence yield taking into account the energy spectrum of electrons in showers initiated by proton and iron nuclei as a function of height. We compare our results with previous calculations based on mono energetic electrons (80 MeV) and a difference in excess of 8% is found. Finally, the influence of a more realistic fluorescence yield in the shower energy reconstruction is also discussed.

We study the influence of shower fluctuations, and the possible presence of different nuclear species in the primary cosmic ray spectrum, on the experimental determination of both shower energy and the proton air inelastic cross section from studies of the longitudinal development of atmospheric showers in fluorescence experiments. We investigate the potential of track length integral and shower size at maximum as estimators of shower energy. We find that at very high energy (~10^19-10^20 eV) the error of the total energy assignment is dominated by the dependence on the hadronic interaction model, and is of the order of 5%. At lower energy (~10^17-10^18 eV), the uncertainty of the energy determination due to the limited knowledge of the primary cosmic ray composition is more important. The distribution of depth of shower maximum is discussed as a measure of the proton-air cross section. Uncertainties in a possible experimental measurement of this cross section introduced by intrinsic shower fluctuations, the model of hadronic interactions, and the unknown mixture of primary nuclei in the cosmic radiation are numerically evaluated.

We explore the constant intensity cut method that is widely used for the derivation of the cosmic ray energy spectrum, for comparisons of data obtained at different atmospheric depths, for measuring average shower profiles, and for estimates of the proton-air cross section from extensive air shower data. The constant intensity cut method is based on the selection of air showers by charged particle or muon size and therefore is subject to intrinsic shower fluctuations. We demostrate that, depending on the selection method, shower fluctuations can strongly influence the characteristics of the selected showers. Furthermore, a mixture of different primaries in the cosmic ray flux complicates the interpretation of measurements based on the method of constant intensity cuts. As an example we consider data published by the Akeno Collaboration. The interpretation of the Akeno measurements suggests that more than 60-70% of cosmic ray primaries in the energy range 10^{16}-19^{17} eV are heavy nuclei. Our conclusions depend only weakly on the hadronic interaction model chosen to perform the simulations, namely SIBYLL and QGSjet.

We report on coordinated thermal and optical measurements of trans-Neptunian object (20000) Varuna obtained in January-February 2002, respectively from the IRAM 30-m and IAA 1.5 m telescopes. The optical data show a lightcurve with a period of 3.176+/-0.010 hr, a mean V magnitude of 20.37+/-0.08 and a 0.42+/-0.01 magnitude amplitude. They also tentatively indicate that the lightcurve is asymmetric and double-peaked. The thermal observations indicate a 1.12+/-0.41 mJy flux, averaged over the object's rotation. Combining the two datasets, we infer that Varuna has a mean 1060(+180/-220) km diameter and a mean 0.038(+0.022/-0.010) V geometric albedo, in general agreement with an earlier determination using the same technique.

We present a fast one dimensional hybrid method to efficiently simulate extensive air showers up to the highest observed energies. Based on precalculated pion showers and a bootstrap technique, our method predicts the average shower profile, the number of muons at detector level above several energy thresholds as well as the fluctuations of the electromagnetic and hadronic components of the shower. We study the main characteristics of proton-induced air showers up to ultra-high energy, comparing the predictions of three different hadronic interaction models: SIBYLL 1.7, SIBYLL 2.1 and QGSjet98. The influence of the hadronic interaction models on the shower evolution, in particular the elongation rate, is discussed and the applicability of analytical approximations is investigated.