Distance transform: a tool for the study of animal colour patterns

Summary The information in animal colour patterns plays a key role in many ecological interactions; quantification would help us to study them, but this is problematic. Comparing patterns using human judgement is subjective and inconsistent. Traditional shape analysis is unsuitable as patterns do not usually contain conserved landmarks. Alternative statistical approaches also have weaknesses, particularly as they are generally based on summary measures that discard most or all of the spatial information in a pattern. We present a method for quantifying the similarity of a pair of patterns based on the distance transform of a binary image. The method compares the whole pattern, pixel by pixel, while being robust to small spatial variations among images. We demonstrate the utility of the distance transform method using three ecological examples. We generate a measure of mimetic accuracy between hoverflies (Diptera: Syrphidae) and wasps (Hymenoptera) based on abdominal pattern and show that this correlates strongly with the perception of a model predator (humans). We calculate similarity values within a group of mimetic butterflies and compare this with proposed pairings of Müllerian comimics. Finally, we characterise variation in clypeal badges of a paper wasp (Polistes dominula) and compare this with previous measures of variation. While our results generally support the findings of existing studies that have used simpler ad hoc methods for measuring differences between patterns, our method is able to detect more subtle variation and hence reveal previously overlooked trends

Azimuthal asymmetry in the risetime of the surface detector signals of the Pierre Auger Observatory

The azimuthal asymmetry in the risetime of signals in Auger surface detector stations is a source of information on shower development. The azimuthal asymmetry is due to a combination of the longitudinal evolution of the shower and geometrical effects related to the angles of incidence of the particles into the detectors. The magnitude of the effect depends upon the zenith angle and state of development of the shower and thus provides a novel observable, $(\sec \theta)_\mathrm{max}$, sensitive to the mass composition of cosmic rays above $3 \times 10^{18}$ eV. By comparing measurements with predictions from shower simulations, we find for both of our adopted models of hadronic physics (QGSJETII-04 and EPOS-LHC) an indication that the mean cosmic-ray mass increases slowly with energy, as has been inferred from other studies. However, the mass estimates are dependent on the shower model and on the range of distance from the shower core selected. Thus the method has uncovered further deficiencies in our understanding of shower modelling that must be resolved before the mass composition can be inferred from $(\sec \theta)_\mathrm{max}$.Comment: Replaced with published version. Added journal reference and DO

Ultrahigh-energy neutrino follow-up of Gravitational Wave events GW150914 and GW151226 with the Pierre Auger Observatory

On September 14, 2015 the Advanced LIGO detectors observed their first gravitational-wave (GW) transient GW150914. This was followed by a second GW event observed on December 26, 2015. Both events were inferred to have arisen from the merger of black holes in binary systems. Such a system may emit neutrinos if there are magnetic fields and disk debris remaining from the formation of the two black holes. With the surface detector array of the Pierre Auger Observatory we can search for neutrinos with energy above 100 PeV from point-like sources across the sky with equatorial declination from about -65 deg. to +60 deg., and in particular from a fraction of the 90% confidence-level (CL) inferred positions in the sky of GW150914 and GW151226. A targeted search for highly-inclined extensive air showers, produced either by interactions of downward-going neutrinos of all flavors in the atmosphere or by the decays of tau leptons originating from tau-neutrino interactions in the Earth's crust (Earth-skimming neutrinos), yielded no candidates in the Auger data collected within $\pm 500$ s around or 1 day after the coordinated universal time (UTC) of GW150914 and GW151226, as well as in the same search periods relative to the UTC time of the GW candidate event LVT151012. From the non-observation we constrain the amount of energy radiated in ultrahigh-energy neutrinos from such remarkable events.Comment: Published version. Added journal reference and DOI. Added Report Numbe

Multi-resolution anisotropy studies of ultrahigh-energy cosmic rays detected at the Pierre Auger Observatory

We report a multi-resolution search for anisotropies in the arrival directions of cosmic rays detected at the Pierre Auger Observatory with local zenith angles up to $80^\circ$ and energies in excess of 4 EeV ($4 \times 10^{18}$ eV). This search is conducted by measuring the angular power spectrum and performing a needlet wavelet analysis in two independent energy ranges. Both analyses are complementary since the angular power spectrum achieves a better performance in identifying large-scale patterns while the needlet wavelet analysis, considering the parameters used in this work, presents a higher efficiency in detecting smaller-scale anisotropies, potentially providing directional information on any observed anisotropies. No deviation from isotropy is observed on any angular scale in the energy range between 4 and 8 EeV. Above 8 EeV, an indication for a dipole moment is captured; while no other deviation from isotropy is observed for moments beyond the dipole one. The corresponding $p$-values obtained after accounting for searches blindly performed at several angular scales, are $1.3 \times 10^{-5}$ in the case of the angular power spectrum, and $2.5 \times 10^{-3}$ in the case of the needlet analysis. While these results are consistent with previous reports making use of the same data set, they provide extensions of the previous works through the thorough scans of the angular scales.Comment: Published version. Added journal reference and DOI. Added Report Numbe

Calibration of the Logarithmic-Periodic Dipole Antenna (LPDA) Radio Stations at the Pierre Auger Observatory using an Octocopter

An in-situ calibration of a logarithmic periodic dipole antenna with a frequency coverage of 30 MHz to 80 MHz is performed. Such antennas are part of a radio station system used for detection of cosmic ray induced air showers at the Engineering Radio Array of the Pierre Auger Observatory, the so-called Auger Engineering Radio Array (AERA). The directional and frequency characteristics of the broadband antenna are investigated using a remotely piloted aircraft (RPA) carrying a small transmitting antenna. The antenna sensitivity is described by the vector effective length relating the measured voltage with the electric-field components perpendicular to the incoming signal direction. The horizontal and meridional components are determined with an overall uncertainty of 7.4^{+0.9}_{-0.3} % and 10.3^{+2.8}_{-1.7} % respectively. The measurement is used to correct a simulated response of the frequency and directional response of the antenna. In addition, the influence of the ground conductivity and permittivity on the antenna response is simulated. Both have a negligible influence given the ground conditions measured at the detector site. The overall uncertainties of the vector effective length components result in an uncertainty of 8.8^{+2.1}_{-1.3} % in the square root of the energy fluence for incoming signal directions with zenith angles smaller than 60{\deg}.Comment: Published version. Updated online abstract only. Manuscript is unchanged with respect to v2. 39 pages, 15 figures, 2 table

Evidence for a mixed mass composition at the `ankle' in the cosmic-ray spectrum

We report a first measurement for ultra-high energy cosmic rays of the correlation between the depth of shower maximum and the signal in the water Cherenkov stations of air-showers registered simultaneously by the fluorescence and the surface detectors of the Pierre Auger Observatory. Such a correlation measurement is a unique feature of a hybrid air-shower observatory with sensitivity to both the electromagnetic and muonic components. It allows an accurate determination of the spread of primary masses in the cosmic-ray flux. Up till now, constraints on the spread of primary masses have been dominated by systematic uncertainties. The present correlation measurement is not affected by systematics in the measurement of the depth of shower maximum or the signal in the water Cherenkov stations. The analysis relies on general characteristics of air showers and is thus robust also with respect to uncertainties in hadronic event generators. The observed correlation in the energy range around the `ankle' at $\lg(E/{\rm eV})=18.5-19.0$ differs significantly from expectations for pure primary cosmic-ray compositions. A light composition made up of proton and helium only is equally inconsistent with observations. The data are explained well by a mixed composition including nuclei with mass $A > 4$. Scenarios such as the proton dip model, with almost pure compositions, are thus disfavoured as the sole explanation of the ultrahigh-energy cosmic-ray flux at Earth.Comment: Published version. Added journal reference and DOI. Added Report Numbe

Cosmic-Ray Anisotropies in Right Ascension Measured by the Pierre Auger Observatory

We present measurements of the large-scale cosmic-ray (CR) anisotropies in R.A., using data collected by the surface detector array of the Pierre Auger Observatory over more than 14 yr. We determine the equatorial dipole component, through a Fourier analysis in R.A. 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 subarray of detectors with 750 m separation, which allows us to extend the analysis down to ∼0.03 EeV. The most significant equatorial dipole amplitude obtained is that in the cumulative bin above 8 EeV, %, which is inconsistent with isotropy at the 6σ level. In the bins below 8 EeV, we obtain 99% CL upper bounds on d ⊥ at the level of 1%-3%. 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 R.A. of the Galactic center, at GC =-94°, suggesting a predominantly Galactic origin for anisotropies at these energies. The reconstructed dipole phases in the energy bins above 4 EeV point instead to R.A. that are almost opposite to the Galactic center one, indicative of an extragalactic CR origin

Update on the Combined Analysis of Muon Measurements from Nine Air Shower Experiments

Over the last two decades, various experiments have measured muon densities in extensive air showers over several orders of magnitude in primary energy. While some experiments observed differences in the muon densities between simulated and experimentally measured air showers, others reported no discrepancies. We will present an update of the meta-analysis of muon measurements from nine air shower experiments, covering shower energies between a few PeV and tens of EeV and muon threshold energies from a few 100 MeV to about 10GeV. In order to compare measurements from different experiments, their energy scale was cross-calibrated and the experimental data has been compared using a universal reference scale based on air shower simulations. Above 10 PeV, we find a muon excess with respect to simulations for all hadronic interaction models, which is increasing with shower energy. For EPOS-LHC and QGSJet-II.04 the significance of the slope of the increase is analyzed in detail under different assumptions of the individual experimental uncertainties