Extensive air showers, originating from ultra-high energy cosmic rays, have been successfully measured through the use of arrays of water-Cherenkov detectors (WCDs). Sophisticated analyses exploiting WCD data have made it possible to demonstrate that shower simulations, based on different hadronic-interaction models, cannot reproduce the observed number of muons at the ground. The accurate knowledge of the WCD response to muons is paramount in establishing the exact level of this discrepancy. In this work, we report on a study of the response of a WCD of the Pierre Auger Observatory to atmospheric muons performed with a hodoscope made of resistive plate chambers (RPCs), enabling us to select and reconstruct nearly 600 thousand single muon trajectories with zenith angles ranging from 0$^\circ$ to 55$^\circ$. Comparison of distributions of key observables between the hodoscope data and the predictions of dedicated simulations allows us to demonstrate the accuracy of the latter at a level of 2%. As the WCD calibration is based on its response to atmospheric muons, the hodoscope data are also exploited to show the long-term stability of the procedure.

We search for signals of magnetically-induced effects in the arrival directions of ultra-high-energy cosmic rays detected at the Pierre Auger Observatory. We apply two different methods. One is a search for sets of events that show a correlation between their arrival direction and the inverse of their energy, which would be expected if they come from the same point-like source, they have the same electric charge and their deflection is relatively small and coherent. We refer to these sets of events as "multiplets". The second method, called "thrust", is a principal axis analysis aimed to detect the elongated patterns in a region of interest. We study the sensitivity of both methods using a benchmark simulation and we apply them to data in two different searches. The first search is done assuming as source candidates a list of nearby active galactic nuclei and starburst galaxies. The second is an all-sky blind search. We report the results and we find no statistically significant features. We discuss the compatibility of these results with the indications on the mass composition inferred from data of the Pierre Auger Observatory.

We present measurements of the large-scale cosmic-ray anisotropies in right ascension, using data collected by the surface detector array of the Pierre Auger Observatory over more than 14 years. We determine the equatorial dipole component, $\vec{d}_\perp$, through a Fourier analysis in right ascension that includes weights for each event so as to account for the main detector-induced systematic effects. For the energies at which the trigger efficiency of the array is small, the ``East-West'' method is employed. Besides using the data from the array with detectors separated by 1500 m, we also include data from the smaller but denser sub-array of detectors with 750 m separation, which allows us to extend the analysis down to $\sim 0.03$ EeV. The most significant equatorial dipole amplitude obtained is that in the cumulative bin above 8~EeV, $d_\perp=6.0^{+1.0}_{-0.9}$%, which is inconsistent with isotropy at the 6$\sigma$ level. In the bins below 8 EeV, we obtain 99% CL upper-bounds on $d_\perp$ at the level of 1 to 3 percent. At energies below 1 EeV, even though the amplitudes are not significant, the phases determined in most of the bins are not far from the right ascension of the Galactic center, at $\alpha_{\rm GC}=-94^\circ$, suggesting a predominantly Galactic origin for anisotropies at these energies. The reconstructed dipole phases in the energy bins above 4 EeV point instead to right ascensions that are almost opposite to the Galactic center one, indicative of an extragalactic cosmic ray origin.

With the Surface Detector array (SD) of the Pierre Auger Observatory we can detect neutrinos with energy between $10^{17}\,$eV and $10^{20}\,$eV from point-like sources across the sky, from close to the Southern Celestial Pole up to $60^\circ$ in declination, with peak sensitivities at declinations around $\sim -53^\circ$ and $\sim+55^\circ$, and an unmatched sensitivity for arrival directions in the Northern hemisphere. A search has been performed for highly-inclined air showers induced by neutrinos of all flavours with no candidate events found in data taken between 1 Jan 2004 and 31 Aug 2018. Upper limits on the neutrino flux from point-like steady sources have been derived as a function of source declination. An unrivaled sensitivity is achieved in searches for transient sources with emission lasting over an hour or less, if they occur within the field of view corresponding to the zenith angle range between $60^\circ$ and $~95^\circ$ where the SD of the Pierre Auger Observatory is most sensitive to neutrinos.

The determination of the primary energy of extensive air showers using the fluorescence detection technique requires an estimation of the energy carried away by particles that do not deposit all their energy in the atmosphere. This estimation is typically made using Monte Carlo simulations and thus depends on the assumed primary particle mass and on model predictions for neutrino and muon production. In this work we present a new method to obtain the invisible energy from events detected by the Pierre Auger Observatory. The method uses measurements of the muon number at ground level, and it allows us to reduce significantly the systematic uncertainties related to the mass composition and the high energy hadronic interaction models, and consequently to improve the estimation of the energy scale of the Observatory.

Neutrinos with energies above $10^{17}$ eV are detectable with the Surface Detector Array of the Pierre Auger Observatory. The identification is efficiently performed for neutrinos of all flavors interacting in the atmosphere at large zenith angles, as well as for Earth-skimming $\tau$ neutrinos with nearly tangential trajectories relative to the earth. No neutrino candidates were found in $\sim\,14.7$ years of data taken up to 31 August 2018. This leads to restrictive upper bounds on their flux. The $90\%$ C.L. single-flavor limit to the diffuse flux of ultra-high-energy neutrinos with an $E_\nu^{-2}$ spectrum in the energy range $1.0 \times 10^{17}~{\rm eV} - 2.5 \times 10^{19}~{\rm eV}$ is $E^2 {\rm d}N_\nu/{\rm d}E_\nu < 4.4 \times 10^{-9}~{\rm GeV~cm^{-2}~s^{-1}~sr^{-1}}$, placing strong constraints on several models of neutrino production at EeV energies and on the properties of the sources of ultra-high-energy cosmic rays.

Contributions of the Pierre Auger Collaboration to the 36th International Cosmic Ray Conference (ICRC 2019), 24 July - 1 August 2019, Madison, Wisconsin, USA.

The profile of the longitudinal development of showers produced by ultra-high energy cosmic rays carries information related to the interaction properties of the primary particles with atmospheric nuclei. In this work, we present the first measurement of the average shower profile in traversed atmospheric depth at the Pierre Auger Observatory. The shapes of profiles are well reproduced by the Gaisser-Hillas parametrization within the range studied, for E > 10^{17.8} eV. A detailed analysis of the systematic uncertainties is performed using 10 years of data and a full detector simulation. The average shape is quantified using two variables related to the width and asymmetry of the profile, and the results are compared with predictions of hadronic interaction models for different primary particles.

An overview of the multi-messenger capabilities of the Pierre Auger Observatory is presented. The techniques and performance of searching for Ultra-High Energy neutrinos, photons and neutrons are described. Some of the most relevant results are reviewed, such as stringent upper bounds that were placed to a flux of diffuse cosmogenic neutrinos and photons, bounds placed on neutrinos emitted from compact binary mergers that were detected by LIGO and Virgo during their first and second observing runs, as well as searches for high energy photons and neutrons from the Galactic center that constrain the properties of the putative Galactic PeVatron, observed by the H.E.S.S.\ collaboration. The observation of directional correlations between ultra-high energy cosmic rays and either high energy astrophysical neutrinos or specific source populations, weighted by their electromagnetic radiation, are also discussed. They constitute additional multi-messenger approaches aimed at identifying the sources of high energy cosmic rays.

We report a comprehensive analysis of the global spectrophotometric properties of Ceres using Dawn Framing Camera images collected from April to June 2015 during the RC3 and Survey mission phases. The single-scattering albedo of Ceres at 555 nm is 0.14$\pm$0.04, the geometric albedo is 0.096$\pm$0.006, and the Bond albedo is 0.037$\pm$0.002. The asymmetry factors calculated from the best-fit two-term Henyey-Greenstein (HG) single-particle phase function (SPPF) show a wavelength dependence, suggesting that the phase reddening of Ceres is dominated by single-particle scattering rather than multiple scattering or small-scale surface roughness. The Hapke roughness parameter of Ceres is derived to be 20$^\circ\pm$6$^\circ$ with no wavelength dependence. The phase function of Ceres shows appreciably strong scattering around 90$^\circ$ phase angle that cannot be fitted with a single-term HG SPPF, suggesting possible stronger forward scattering than other asteroids previously analyzed with spacecraft data. We speculate that such a scattering characteristic of Ceres might be related to its unique surface composition. We grouped the reflectance data into a 1$^\circ$ latitude-longitude grid and fitted each grid independently to study the spatial variations of photometric properties. The albedo and color maps are consistent with previous studies. The SPPF over the surface of Ceres shows stronger backscattering associated with lower albedo and vice versa, consistent with the general trend among asteroids. The Hapke roughness parameter does not vary much across the surface of Ceres, except for the ancient Vendimia Planitia region that has a slightly higher roughness. Based on the wavelength dependence of the SPPF of Ceres, we hypothesize that its regolith grains either contain a considerable fraction of $\lessapprox\mu$m-sized particles, or are strongly affected by internal scatterers of this size.

We present a detailed study of the large-scale anisotropies of cosmic rays with energies above 4 EeV measured using the Pierre Auger Observatory. For the energy bins [4,8] EeV and $E\geq 8$ EeV, the most significant signal is a dipolar modulation in right ascension at energies above 8 EeV, as previously reported. In this paper we further scrutinize the highest-energy bin by splitting it into three energy ranges. We find that the amplitude of the dipole increases with energy above 4 EeV. The growth can be fitted with a power law with index $\beta=0.79\pm 0.19$. The directions of the dipoles are consistent with an extragalactic origin of these anisotropies at all the energies considered. Additionally we have estimated the quadrupolar components of the anisotropy: they are not statistically significant. We discuss the results in the context of the predictions from different models for the distribution of ultrahigh-energy sources and cosmic magnetic fields.

The surface reflectance of planetary regoliths may increase dramatically towards zero phase angle, a phenomenon known as the opposition effect (OE). Two physical processes that are thought to be the dominant contributors to the brightness surge are shadow hiding (SH) and coherent backscatter (CB). The occurrence of shadow hiding in planetary regoliths is self-evident, but it has proved difficult to unambiguously demonstrate CB from remote sensing observations. One prediction of CB theory is the wavelength dependence of the OE angular width. The Dawn spacecraft observed the OE on the surface of dwarf planet Ceres. We characterize the OE over the resolved surface, including the bright Cerealia Facula, and to find evidence for SH and/or CB. We analyze images of the Dawn framing camera by means of photometric modeling of the phase curve. We find that the OE of most of the investigated surface has very similar characteristics, with an enhancement factor of 1.4 and a FWHM of 3{\deg} (broad OE). A notable exception are the fresh ejecta of the Azacca crater, which display a very narrow brightness enhancement that is restricted to phase angles $< 0.5${\deg} (narrow OE); suggestively, this is in the range in which CB is thought to dominate. We do not find a wavelength dependence for the width of the broad OE, and lack the data to investigate the dependence for the narrow OE. The prediction of a wavelength-dependent CB width is rather ambiguous. The zero-phase observations allow us to determine Ceres' visible geometric albedo as $p_V = 0.094 \pm 0.005$. A comparison with other asteroids suggests that Ceres' broad OE is typical for an asteroid of its spectral type, with characteristics that are primarily linked to surface albedo. Our analysis suggests that CB may occur on the dark surface of Ceres in a highly localized fashion.

With the Auger Engineering Radio Array (AERA) of the Pierre Auger Observatory, we have observed the radio emission from 561 extensive air showers with zenith angles between 60$^\circ$ and 84$^\circ$. In contrast to air showers with more vertical incidence, these inclined air showers illuminate large ground areas of several km$^2$ with radio signals detectable in the 30 to 80\,MHz band. A comparison of the measured radio-signal amplitudes with Monte Carlo simulations of a subset of 50 events for which we reconstruct the energy using the Auger surface detector shows agreement within the uncertainties of the current analysis. As expected for forward-beamed radio emission undergoing no significant absorption or scattering in the atmosphere, the area illuminated by radio signals grows with the zenith angle of the air shower. Inclined air showers with EeV energies are thus measurable with sparse radio-antenna arrays with grid sizes of a km or more. This is particularly attractive as radio detection provides direct access to the energy in the electromagnetic cascade of an air shower, which in case of inclined air showers is not accessible by arrays of particle detectors on the ground.

Asteroid (162173) Ryugu is the target object of Hayabusa2, an asteroid exploration and sample return mission led by Japan Aerospace Exploration Agency (JAXA). Ground-based observations indicate that Ryugu is a C-type near-Earth asteroid with a diameter of less than 1 km, but the knowledge of its detailed properties is still very limited. This paper summarizes our best understanding of the physical and dynamical properties of Ryugu based on remote sensing and theoretical modeling. This information is used to construct a design reference model of the asteroid that is used for formulation of mission operations plans in advance of asteroid arrival. Particular attention is given to the surface properties of Ryugu that are relevant to sample acquisition. This reference model helps readers to appropriately interpret the data that will be directly obtained by Hayabusa2 and promotes scientific studies not only for Ryugu itself and other small bodies but also for the Solar System evolution that small bodies shed light on.

A new analysis of the dataset from the Pierre Auger Observatory provides evidence for anisotropy in the arrival directions of ultra-high-energy cosmic rays on an intermediate angular scale, which is indicative of excess arrivals from strong, nearby sources. The data consist of 5514 events above 20 EeV with zenith angles up to 80 deg recorded before 2017 April 30. Sky models have been created for two distinct populations of extragalactic gamma-ray emitters: active galactic nuclei from the second catalog of hard Fermi-LAT sources (2FHL) and starburst galaxies from a sample that was examined with Fermi-LAT. Flux-limited samples, which include all types of galaxies from the Swift-BAT and 2MASS surveys, have been investigated for comparison. The sky model of cosmic-ray density constructed using each catalog has two free parameters, the fraction of events correlating with astrophysical objects and an angular scale characterizing the clustering of cosmic rays around extragalactic sources. A maximum-likelihood ratio test is used to evaluate the best values of these parameters and to quantify the strength of each model by contrast with isotropy. It is found that the starburst model fits the data better than the hypothesis of isotropy with a statistical significance of 4.0 sigma, the highest value of the test statistic being for energies above 39 EeV. The three alternative models are favored against isotropy with 2.7-3.2 sigma significance. The origin of the indicated deviation from isotropy is examined and prospects for more sensitive future studies are discussed.

We present a new method for probing the hadronic interaction models at ultra-high energy and extracting details about mass composition. This is done using the time profiles of the signals recorded with the water-Cherenkov detectors of the Pierre Auger Observatory. The profiles arise from a mix of the muon and electromagnetic components of air-showers. Using the risetimes of the recorded signals we define a new parameter, which we use to compare our observations with predictions from simulations. We find, firstly, inconsistencies between our data and predictions over a greater energy range and with substantially more events than in previous studies. Secondly, by calibrating the new parameter with fluorescence measurements from observations made at the Auger Observatory, we can infer the depth of shower maximum for a sample of over 81,000 events extending from 0.3 EeV to over 100 EeV. Above 30 EeV, the sample is nearly fourteen times larger than currently available from fluorescence measurements and extending the covered energy range by half a decade. The energy dependence of the average depth of shower maximum is compared to simulations and interpreted in terms of the mean of the logarithmic mass. We find good agreement with previous work and extend the measurement of the mean depth of shower maximum to greater energies than before, reducing significantly the statistical uncertainty associated with the inferences about mass composition.

Contributions of the Pierre Auger Collaboration to the 35th International Cosmic Ray Conference (ICRC 2017), 12-20 July 2017, Bexco, Busan, Korea.

We present a novel method to measure precisely the relative spectral response of the fluorescence telescopes of the Pierre Auger Observatory. We used a portable light source based on a xenon flasher and a monochromator to measure the relative spectral efficiencies of eight telescopes in steps of 5 nm from 280 nm to 440 nm. Each point in a scan had approximately 2 nm FWHM out of the monochromator. Different sets of telescopes in the observatory have different optical components, and the eight telescopes measured represent two each of the four combinations of components represented in the observatory. We made an end-to-end measurement of the response from different combinations of optical components, and the monochromator setup allowed for more precise and complete measurements than our previous multi-wavelength calibrations. We find an overall uncertainty in the calibration of the spectral response of most of the telescopes of 1.5% for all wavelengths; the six oldest telescopes have larger overall uncertainties of about 2.2%. We also report changes in physics measureables due to the change in calibration, which are generally small.

We collect observational evidence that supports the scheme of mass transfer on the nucleus of comet 67P/Churyumov-Gerasimenko. The obliquity of the rotation axis of 67P causes strong seasonal variations. During perihelion the southern hemisphere is four times more active than the north. Northern territories are widely covered by granular material that indicates back fall originating from the active south. Decimetre sized chunks contain water ice and their trajectories are influenced by an anti-solar force instigated by sublimation. OSIRIS observations suggest that up to 20 % of the particles directly return to the nucleus surface taking several hours of travel time. The back fall covered northern areas are active if illuminated but produce mainly water vapour. The decimetre chunks from the nucleus surface are too small to contain more volatile compounds such as CO 2 or CO. This causes a north-south dichotomy of the composition measurements in the coma. Active particles are trapped in the gravitational minimum of Hapi during northern winter. They are "shock frozen" and only reactivated when the comet approaches the sun after its aphelion passage. The insolation of the big cavity is enhanced by self-heating, i. e. reflection and IR radiation from the walls. This, together with the pristinity of the active back fall, explains the early observed activity of the Hapi region. Sobek may be a role model for the consolidated bottom of Hapi. Mass transfer in the case of 67P strongly influences the evolution of the nucleus and the interpretation of coma measurements.

We study the spectrophotometric properties of dwarf planet Ceres in the VIS-IR spectral range by means of hyper-spectral images acquired by the VIR imaging spectrometer on board the NASA Dawn mission. Disk-resolved observations with a phase angle within the $7^{\circ}<\alpha<132^{\circ}$ interval were used to characterize Ceres' phase curve in the 0.465-4.05 $\mu$m spectral range. Hapke's model was applied to perform the photometric correction of the dataset, allowing us to produce albedo and color maps of the surface. The $V$-band magnitude phase function of Ceres was fitted with both the classical linear model and H-G formalism. The single-scattering albedo and the asymmetry parameter at 0.55$\mu$m are $w=0.14\pm0.02$ and $\xi=-0.11\pm0.08$, respectively (two-lobe Henyey-Greenstein phase function); the modeled geometric albedo is $0.094\pm0.007$; the roughness parameter is $\bar{\theta}=29^{\circ}\pm6^{\circ}$. Albedo maps indicate small variability on a global scale with an average reflectance of $0.034 \pm 0.003$. Isolated areas such as the Occator bright spots, Haulani, and Oxo show an albedo much higher than average. We measure a significant spectral phase reddening, and the average spectral slope of Ceres' surface after photometric correction is $1.1\%k\AA^{-1}$ and $0.85\%k\AA^{-1}$ at VIS and IR wavelengths, respectively. Broadband color indices are $V-R=0.38\pm0.01$ and $R-I=0.33\pm0.02$. H-G modeling of the $V$-band magnitude phase curve for $\alpha<30^{\circ}$ gives $H=3.14\pm0.04$ and $G=0.10\pm0.04$, while the classical linear model provides $V(1,1,0^{\circ})=3.48\pm0.03$ and $\beta=0.036\pm0.002$. The comparison with spectrophotometric properties of other minor bodies indicates that Ceres has a less back-scattering phase function and a slightly higher albedo than comets and C-type objects. However, the latter represents the closest match in the usual asteroid taxonomy.

We present a global spectrophotometric characterization of the Ceres surface using Dawn Framing Camera (FC) images. We identify the photometric model that yields the best results for photometrically correcting images. Corrected FC images acquired on approach to Ceres were assembled into global maps of albedo and color. Generally, albedo and color variations on Ceres are muted. The albedo map is dominated by a large, circular feature in Vendimia Planitia, known from HST images (Li et al., 2006), and dotted by smaller bright features mostly associated with fresh-looking craters. The dominant color variation over the surface is represented by the presence of "blue" material in and around such craters, which has a negative spectral slope over the visible wavelength range when compared to average terrain. We also mapped variations of the phase curve by employing an exponential photometric model, a technique previously applied to asteroid Vesta (Schr\"oder et al., 2013b). The surface of Ceres scatters light differently from Vesta in the sense that the ejecta of several fresh-looking craters may be physically smooth rather than rough. High albedo, blue color, and physical smoothness all appear to be indicators of youth. The blue color may result from the desiccation of ejected material that is similar to the phyllosilicates/water ice mixtures in the experiments of Poch et al. (2016). The physical smoothness of some blue terrains would be consistent with an initially liquid condition, perhaps as a consequence of impact melting of subsurface water ice. We find red terrain (positive spectral slope) near Ernutet crater, where De Sanctis et al. (2017) detected organic material. The spectrophotometric properties of the large Vendimia Planitia feature suggest it is a palimpsest, consistent with the Marchi et al. (2016) impact basin hypothesis. The central bright area in Occator crater, Cerealia...

The Descent Imager/Spectral Radiometer aboard the Huygens probe successfully acquired images and spectra of the surface of Titan. To counter the effects of haze and atmospheric methane absorption it carried a Surface Science Lamp to illuminate the surface just before landing. We reconstruct the reflectance spectrum of the landing site in the 500-1500 nm range from Downward Looking Visual and Infrared Spectrometer data that show evidence of lamp light. Our reconstruction is a followup to the analysis by Tomasko et al. (2005), who scaled their result to the ratio of the up- and down flux measured just before landing. Instead, we use the lamp flux from the calibration experiment, and find a significantly higher overall reflectance. We attribute this to a phase angle dependance, possibly representing the opposition surge commonly encountered on solar system bodies. The reconstruction in the visible wavelength range is greatly improved. Here, the reflectance spectrum features a red slope, consistent with the presence of organic material. We confirm the blue slope in the near-IR, featureless apart from a single shallow absorption feature at 1500 nm. We agree with Tomasko et al. that the evidence for water ice is inconclusive. By modeling of absorption bands we find a methane mixing ratio of 4.5% +/- 0.5% just above the surface. There is no evidence for the presence of liquid methane, but the data do not rule out a wet soil at a depth of several centimeters.

While landing on Titan, several instruments onboard Huygens acquired measurements that indicate the probe did not immediately come to rest. Detailed knowledge of the probe's motion can provide insight into the nature of Titan's surface. Combining accelerometer data from the Huygens Atmospheric Structure Instrument (HASI) and the Surface Science Package (SSP) with photometry data from the Descent Imager/Spectral Radiometer (DISR) we develop a quantitative model to describe motion of the probe, and its interaction with the surface. The most likely scenario is the following. Upon impact, Huygens created a 12 cm deep hole in the surface of Titan. It bounced back, out of the hole onto the flat surface, after which it commenced a 30-40 cm long slide in the southward direction. The slide ended with the probe out of balance, tilted in the direction of DISR by around 10 degrees. The probe then wobbled back and forth five times in the north-south direction, during which it probably encountered a 1-2 cm sized pebble. The SSP provides evidence for movement up to 10 s after impact. This scenario puts the following constraints on the physical properties of the surface. For the slide over the surface we determine a friction coefficient of 0.4. While this value is not necessarily representative for the surface itself due to the presence of protruding structures on the bottom of the probe, the dynamics appear to be consistent with a surface consistency of damp sand. Additionally, we find that spectral changes observed in the first four seconds after landing are consistent with a transient dust cloud, created by the impact of the turbulent wake behind the probe on the surface. The optical properties of the dust particles are consistent with those of Titan aerosols from Tomasko et al. (P&SS 56, 669). We suggest that the surface at the landing site was covered by a dust layer, possibly the 7 mm layer of...

The OSIRIS camera onboard Rosetta successfully acquired images of asteroid 2867 Steins through a variety of color filters during the flyby on 5 September 2008. The best images of this 5 km diameter asteroid have a resolution of 78 meters per pixel. We process the images by deconvolving with the point spread function and enlarging through the Mitchell-Netravali filter. The enhanced set is analyzed by means of various techniques (PCA, band ratios, stereo anaglyphs) to study surface morphology and search for variegation. We identify a landslide, which supports a YORP origin for Steins' unusual diamond shape. In addition, we find that the interior of one of two large craters on the south pole is bluer than the rest of the body.

During its year-long orbital mission, the Dawn spacecraft has mapped the surface of main-belt asteroid Vesta multiple times at different spatial resolutions and illumination and viewing angles. The onboard Framing Camera has acquired thousands of clear filter and narrow band images, which, with the availability of high-resolution global shape models, allows for a photometric characterization of the surface in unprecedented detail. We analyze clear filter images to retrieve the photometric properties of the regolith. In the first part of the paper we evaluate different photometric models for the global average. In the second part we use these results to study variations in albedo and steepness of the phase curve over the surface. Maps of these two photometric parameters show large scale albedo variations, which appear to be associated with compositional differences. They also reveal the location of photometrically extreme terrains, where the phase curve is unusually shallow or steep. We find that shallow phase curves are associated with steep slopes on crater walls and faults, as calculated from a shape model. On the other hand, the phase curve of ejecta associated with young impact craters is steep. We interpret these variations in phase curve slope in terms of physical roughness of the regolith. The lack of rough ejecta around older craters suggests that initially rough ejecta associated with impact craters on Vesta are smoothed over a relatively short time of several tens of Myr. We propose that this process is the result of impact gardening, and as such represents a previously unrecognized aspect of Vesta space weathering (Pieters et al., 2012). If this type of space weathering is common, we may expect to encounter this photometric phenomenon on other main belt asteroids.

The NASA Dawn spacecraft acquired thousands of images of asteroid Vesta during its year-long orbital tour, and is now on its way to asteroid Ceres. A method for calibrating images acquired by the onboard Framing Camera was described by Schr\"oder et al. (2013; Icarus 226, 1304). However, their method is only valid for point sources. In this paper we extend the calibration to images of extended sources like Vesta. For this, we devise a first-order correction for in-field stray light, which is known to plague images taken through the narrow band filters, and revise the flat fields that were acquired in an integrating sphere before launch. We used calibrated images of the Vesta surface to construct simple photometric models for all filters, that allow us to study how the spectrum changes with increasing phase angle (phase reddening). In combination with these models, our calibration method can be used to create near-seamless mosaics that are radiometrically accurate to a few percent. Such mosaics are provided in JVesta, the Vesta version of the JMARS geographic information system.

We present a method for calibrating images acquired by the Dawn Framing Camera (FC) that is based on the results of an in-flight calibration campaign performed during the cruise from Earth to Vesta. We describe this campaign and the data analysis in full. Both the primary camera FC2 and the backup camera FC1 are radiometrically and geometrically calibrated through observations of standard stars, star fields, and solar system objects. The calibration in each spectral filter is accurate to within a few percent for point sources. Geometric distortion, small by design, is characterized with high accuracy. Dark current, monitored on a regular basis, is very low at flight operational temperatures. Out-of-field stray light was characterized using the Sun as a stray light source. In-field stray light is confirmed in narrow-band filter images of Vesta. Its magnitude and distribution are scene-dependent, and expected to contribute significantly to images of extended objects. Description of a method for in-field stray light correction is deferred to a follow-up paper, as is a discussion of the closely related topic of flat-fielding.

The visible reflectance spectrum of many solar system bodies changes with changing viewing geometry for reasons not fully understood. It is often observed to redden (increasing spectral slope) with increasing solar phase angle, an effect known as phase reddening. Only once, in an observation of the Martian surface by the Viking 1 lander, was reddening observed up to a certain phase angle with bluing beyond, making the reflectance ratio as a function of phase angle shaped like an arch. However, in laboratory experiments this arch-shape is frequently encountered. To investigate this, we measured the bidirectional reflectance of particulate samples of several common rock types in the 400-1000 nm wavelength range and performed ray-tracing simulations. We confirm the occurrence of the arch for surfaces that are forward scattering, i.e. are composed of semi-transparent particles and are smooth on the scale of the particles, and for which the reflectance increases from the lower to the higher wavelength in the reflectance ratio. The arch shape is reproduced by the simulations, which assume a smooth surface. However, surface roughness on the scale of the particles, such as the Hapke and van Horn (1963) fairy castles that can spontaneously form when sprinkling a fine powder, leads to monotonic reddening. A further consequence of this form of microscopic roughness (being indistinct without the use of a microscope) is a flattening of the disk function at visible wavelengths, i.e. Lommel-Seeliger-type scattering. The experiments further reveal monotonic reddening for reflectance ratios at near-IR wavelengths. The simulations fail to reproduce this particular reddening, and we suspect that it results from roughness on the surface of the particles. Given that the regolith of atmosphereless solar system bodies is composed of small particles, our results indicate that the prevalence of monotonic (...)

During the flyby in 2010, the OSIRIS camera on-board Rosetta acquired hundreds of high-resolution images of asteroid Lutetia's surface through a range of narrow-band filters. While Lutetia appears very bland in the visible wavelength range, Magrin et al. (2012) tentatively identified UV color variations in the Baetica cluster, a group of relatively young craters close to the north pole. As Lutetia remains a poorly understood asteroid, such color variations may provide clues to the nature of its surface. We take the color analysis one step further. First we orthorectify the images using a shape model and improved camera pointing, then apply a variety of techniques (photometric correction, principal component analysis) to the resulting color cubes. We characterize variegation in the Baetica crater cluster at high spatial resolution, identifying crater rays and small, fresh impact craters. We argue that at least some of the color variation is due to space weathering, which makes Lutetia's regolith redder and brighter.

After coming to rest on the night side of comet 67P/Churyumov-Gerasimenko, the ROLIS camera on-board Rosetta's Philae lander acquired five images of the surface below the lander, four of which were with the aid of LED illumination of different colors. The images confirm that Philae was perched on a sloped surface. A local horizon is visible in one corner of the image, beyond which we can see the coma. Having spent a full day on the surface Philae was commanded to lift and rotate, after which a final, sixth, LED image was acquired. The change in perspective allowed us to construct a shape model of the surface. The distance to the foreground was about 80 cm, much larger than the nominal 30 cm. This caused stray light, rather than directly reflected LED light, to dominate the image signal, complicating the analysis. The images show a lumpy surface with a roughness of apparently fractal nature. Its appearance is completely different from that of the first landing site, which was characterized by centimeter to meter-sized debris (Mottola et al., 2015). We recognize neither particles nor pores at the image resolution of 0.8 mm per pixel and large color variations are absent. The surface has a bi-modal brightness distribution that can be interpreted in terms of the degree of consolidation, a hypothesis that we support with experimental evidence. We propose the surface below the lander to consist of smooth, cracked plates with unconsolidated edges, similar to terrain seen in CIVA images.

We present an extensive study of the BL Lac object Mrk 501 based on a data set collected during the multi-instrument campaign spanning from 2009 March 15 to 2009 August 1 which includes, among other instruments, MAGIC, VERITAS, Whipple 10-m, Fermi-LAT, RXTE, Swift, GASP-WEBT and VLBA. We find an increase in the fractional variability with energy, while no significant interband correlations of flux changes are found in the acquired data set. The higher variability in the very high energy (>100 GeV, VHE) gamma-ray emission and the lack of correlation with the X-ray emission indicate that the highest-energy electrons that are responsible for the VHE gamma-rays do not make a dominant contribution to the ~1 keV emission. Alternatively, there could be a very variable component contributing to the VHE gamma-ray emission in addition to that coming from the synchrotron self-Compton (SSC) scenarios. The space of SSC model parameters is probed following a dedicated grid-scan strategy, allowing for a wide range of models to be tested and offering a study of the degeneracy of model-to-data agreement in the individual model parameters. We find that there is some degeneracy in both the one-zone and the two-zone SSC scenarios that were probed, with several combinations of model parameters yielding a similar model-to-data agreement, and some parameters better constrained than others. The SSC model grid-scan shows that the flaring activity around 2009 May 22 cannot be modeled adequately with a one-zone SSC scenario, while it can be suitably described within a two-independent-zone SSC scenario. The observation of an electric vector polarization angle rotation coincident with the gamma-ray flare from 2009 May 1 resembles those reported previously for low frequency peaked blazars, hence suggesting that there are many similarities in the flaring mechanisms of blazars with different jet properties.

Previous observations suggested that Ceres has active but possibly sporadic water outgassing, and possibly varying spectral characteristics in a time scale of months. We used all available data of Ceres collected in the past three decades from the ground and the Hubble Space Telescope, and the newly acquired images by Dawn Framing Camera to search for spectral and albedo variability on Ceres, in both a global scale and local regions, particularly the bright spots inside Occator crater, over time scales of a few months to decades. Our analysis has placed an upper limit on the possible temporal albedo variation on Ceres. Sporadic water vapor venting, or any possibly ongoing activity on Ceres, is not significant enough to change the albedo or the area of the bright features in Occator crater by >15%, or the global albedo by >3% over various time scales that we searched. Recently reported spectral slope variations can be explained by changing Sun-Ceres-Earth geometry. The active area on Ceres is less than 1 km$^2$, too small to cause global albedo and spectral variations detectable in our data. Impact ejecta due to impacting projectiles of tens of meters in size like those known to cause observable changes to the surface albedo on Asteroid Scheila cannot cause detectable albedo change on Ceres due to its relatively large size and strong gravity. The water vapor activity on Ceres is independent of Ceres' heliocentric distance, rulling out the possibility of comet-like sublimation process as a possible mechanism driving the activity.

Dawn spacecraft orbited Vesta for more than one year and collected a huge volume of multispectral, high-resolution data in the visible wavelengths with the Framing Camera. We present a detailed disk-integrated and disk-resolved photometric analysis using the Framing Camera images with the Minnaert model and the Hapke model, and report our results about the global photometric properties of Vesta. The photometric properties of Vesta show weak or no dependence on wavelengths, except for the albedo. At 554 nm, the global average geometric albedo of Vesta is 0.38+/-0.04, and the Bond albedo range is 0.20+/-0.02. The bolometric Bond albedo is 0.18+/-0.01. The phase function of Vesta is similar to those of S-type asteroids. Vesta's surface shows a single-peaked albedo distribution with a full-width-half-max ~17% relative to the global average. This width is much smaller than the full range of albedos (from ~0.55x to >2x global average) in localized bright and dark areas of a few tens of km in sizes, and is probably a consequence of significant regolith mixing on the global scale. Rheasilvia basin is about 10% brighter than the global average. The phase reddening of Vesta measured from Dawn Framing Camera images is comparable or slightly stronger than that of Eros as measured by the Near Earth Asteroid Rendezvous mission, but weaker than previous measurements based on ground-based observations of Vesta and laboratory measurements of HED meteorites. The photometric behaviors of Vesta are best described by the Hapke model and the Akimov disk- function, when compared with the Minnaert model, Lommel-Seeliger model, and Lommel- Seeliger-Lambertian model. The traditional approach for photometric correction is validated for Vesta for >99% of its surface where reflectance is within +/-30% of global average.