A Monte Carlo Template based analysis for Air-Cherenkov Arrays

We present a high-performance event reconstruction algorithm: an Image Pixel-wise fit for Atmospheric Cherenkov Telescopes (ImPACT). The reconstruction algorithm is based around the likelihood fitting of camera pixel amplitudes to an expected image template. A maximum likelihood fit is performed to find the best-fit shower parameters. A related reconstruction algorithm has already been shown to provide significant improvements over traditional reconstruction for both the CAT and H.E.S.S. experiments. We demonstrate a significant improvement to the template generation step of the procedure, by the use of a full Monte Carlo air shower simulation in combination with a ray-tracing optics simulation to more accurately model the expected camera images. This reconstruction step is combined with an MVA-based background rejection. Examples are shown of the performance of the ImPACT analysis on both simulated and measured (from a strong VHE source) gamma-ray data from the H.E.S.S. array, demonstrating an improvement in sensitivity of more than a factor two in observation time over traditional image moments-fitting methods, with comparable performance to previous likelihood fitting analyses. ImPACT is a particularly promising approach for future large arrays such as the Cherenkov Telescope Array (CTA) due to its improved high-energy performance and suitability for arrays of mixed telescope types.Comment: 13 pages, 10 figure

A Generic Algorithm for IACT Optical Efficiency Calibration using Muons

Muons produced in Extensive Air Showers (EAS) generate ring-like images in Imaging Atmospheric Cherenkov Telescopes when travelling near parallel to the optical axis. From geometrical parameters of these images, the absolute amount of light emitted may be calculated analytically. Comparing the amount of light recorded in these images to expectation is a well established technique for telescope optical efficiency calibration. However, this calculation is usually performed under the assumption of an approximately circular telescope mirror. The H.E.S.S. experiment entered its second phase in 2012, with the addition of a fifth telescope with a non-circular 600m$^2$ mirror. Due to the differing mirror shape of this telescope to the original four H.E.S.S. telescopes, adaptations to the standard muon calibration were required. We present a generalised muon calibration procedure, adaptable to telescopes of differing shapes and sizes, and demonstrate its performance on the H.E.S.S. II array.Comment: In Proceedings of the 34th International Cosmic Ray Conference (ICRC2015), The Hague, The Netherland

Potential for measuring the longitudinal and lateral profile of muons in TeV air showers with IACTs

Muons are copiously produced within hadronic extensive air showers (EAS) occurring in the Earth's atmosphere, and are used by particle air shower detectors as a means of identifying the primary cosmic ray which initiated the EAS. Imaging Atmospheric Cherenkov Telescopes (IACTs), designed for the detection of gamma-ray initiated EAS for the purposes of Very High Energy (VHE) gamma-ray astronomy, are subject to a considerable background signal due to hadronic EAS. Although hadronic EAS are typically rejected for gamma-ray analysis purposes, single muons produced within such showers generate clearly identifiable signals in IACTs and muon images are routinely retained and used for calibration purposes. For IACT arrays operating with a stereoscopic trigger, when a muon triggers one telescope, other telescopes in IACT arrays usually detect the associated hadronic EAS. We demonstrate for the first time the potential of IACT arrays for competitive measurements of the muon content of air showers, their lateral distribution and longitudinal profile of production slant heights in the TeV energy range. Such information can provide useful input to hadronic interaction models.Comment: 15 pages, 11 figures, 2 tables, accepted for publication in Astroparticle Physic

Fluvio-deltaic avulsions during relative sea-level fall.

Understanding river response to changes in relative sea level (RSL) is essential for predicting fluvial stratigraphy and source-to-sink dynamics. Recent theoretical work has suggested that rivers can remain aggradational during RSL fall, but field data are needed to verify this response and investigate sediment deposition processes. We show with field work and modeling that fluvio-deltaic systems can remain aggradational or at grade during RSL fall, leading to superelevation and continuation of delta lobe avulsions. The field site is the Goose River, Newfoundland-Labrador, Canada, which has experienced steady RSL fall of around 3–4 mm yr⁻¹ in the past 5 k.y. from post-glacial isostatic rebound. Elevation analysis and optically stimulated luminescence dating suggest that the Goose River avulsed and deposited three delta lobes during RSL fall. Simulation results from Delft3D software show that if the characteristic fluvial response time is longer than the duration of RSL fall, then fluvial systems remain aggradational or at grade, and continue to avulse during RSL fall due to superelevation. Intriguingly, we find that avulsions become more frequent at faster rates of RSL fall, provided the system response time remains longer than the duration of RSL fall. This work suggests that RSL fall rate may influence the architecture of falling-stage or forced regression deposits by controlling the number of deposited delta lobes

Optimizing passive acoustic sampling of bats in forests

Passive acoustic methods are increasingly used in biodiversity research and monitoring programs because they are cost-effective and permit the collection of large datasets. However, the accuracy of the results depends on the bioacoustic characteristics of the focal taxa and their habitat use. In particular, this applies to bats which exhibit distinct activity patterns in three-dimensionally structured habitats such as forests. We assessed the performance of 21 acoustic sampling schemes with three temporal sampling patterns and seven sampling designs. Acoustic sampling was performed in 32 forest plots, each containing three microhabitats: forest ground, canopy, and forest gap. We compared bat activity, species richness, and sampling effort using species accumulation curves fitted with the clench equation. In addition, we estimated the sampling costs to undertake the best sampling schemes. We recorded a total of 145,433 echolocation call sequences of 16 bat species. Our results indicated that to generate the best outcome, it was necessary to sample all three microhabitats of a given forest location simultaneously throughout the entire night. Sampling only the forest gaps and the forest ground simultaneously was the second best choice and proved to be a viable alternative when the number of available detectors is limited. When assessing bat species richness at the 1-km(2) scale, the implementation of these sampling schemes at three to four forest locations yielded highest labor cost-benefit ratios but increasing equipment costs. Our study illustrates that multiple passive acoustic sampling schemes require testing based on the target taxa and habitat complexity and should be performed with reference to cost-benefit ratios. Choosing a standardized and replicated sampling scheme is particularly important to optimize the level of precision in inventories, especially when rare or elusive species are expected