1. Reconstructing air shower parameters with LOFAR using event specific GDAS atmosphere
- Author
-
Brian Hare, Anna Nelles, Stijn Buitink, Katie Mulrey, Heino Falcke, S. ter Veen, T. N. G. Trinh, Laura Rossetto, H. Pandya, Jörg P. Rachen, G. K. Krampah, Antonio Bonardi, Tim Huege, Jörg R. Hörandel, Pragati Mitra, Arthur Corstanje, Olaf Scholten, Tobias Winchen, Faculty of Sciences and Bioengineering Sciences, Physics, Elementary Particle Physics, Astronomy, and Research unit Astroparticle Physics
- Subjects
Meteorology ,Effects of humidity ,Astronomy ,Astrophysics::High Energy Astrophysical Phenomena ,Monte Carlo method ,FOS: Physical sciences ,Cosmic ray ,Atmospheric model ,01 natural sciences ,Cosmic Ray ,Atmosphere ,Data assimilation ,0103 physical sciences ,Instrumentation and Methods for Astrophysics (astro-ph.IM) ,010303 astronomy & astrophysics ,Physics::Atmospheric and Oceanic Physics ,astro-ph.HE ,High Energy Astrophysical Phenomena (astro-ph.HE) ,Physics ,Atmospheric models ,Radio detection technique ,010308 nuclear & particles physics ,Astrophysics::Instrumentation and Methods for Astrophysics ,Index of refraction ,X-max reconstruction ,Astronomy and Astrophysics ,LOFAR ,EAS ,GDAS ,COSMIC-RAYS ,Air shower ,13. Climate action ,RADIO-EMISSION ,Astrophysics - High Energy Astrophysical Phenomena ,Astrophysics - Instrumentation and Methods for Astrophysics ,astro-ph.IM - Abstract
The limited knowledge of atmospheric parameters like humidity, pressure, temperature, and the index of refraction has been one of the important systematic uncertainties in reconstructing the depth of the shower maximum from the radio emission of air showers. Current air shower Monte Carlo simulation codes like CORSIKA and the radio plug-in CoREAS use various averaged parameterized atmospheres. However, time-dependent and location-specific atmospheric models are needed for the cosmic ray analysis method used for LOFAR data. There, dedicated simulation sets are used for each detected cosmic ray, to take into account the actual atmospheric conditions at the time of the measurement. Using the Global Data Assimilation System (GDAS), a global atmospheric model, we have implemented time-dependent, realistic atmospheric profiles in CORSIKA and CoREAS. We have produced realistic event-specific atmospheres for all air showers measured with LOFAR, an event set spanning several years and many different weather conditions. A complete re-analysis of our data set shows that for the majority of data, our previous correction factor performed rather well; we found only a small systematic shift of 2 g/cm$^2$ in the reconstructed $X_{\rm max}$. However, under extreme weather conditions, for example, very low air pressure, the shift can be up to 15 g/cm$^2$. We provide a correction formula to determine the shift in $X_{\rm max}$ resulting from a comparison of simulations done using the US-Std atmosphere and the GDAS-based atmosphere., Comment: Accepted for publication in Astroparticle Physics. arXiv admin note: text overlap with arXiv:1911.02859
- Published
- 2020