We report an estimation of the injected mass composition of ultrahigh energy cosmic rays (UHECRs) at energies higher than 10 EeV. The composition is inferred from an energy-dependent sky distribution of UHECR events observed by the Telescope Array surface detector by comparing it to the Large Scale Structure of the local Universe. In the case of negligible extragalactic magnetic fields (EGMFs), the results are consistent with a relatively heavy injected composition at E similar to 10 EeV that becomes lighter up to E similar to 100 EeV, while the composition at E > 100 EeV is very heavy. The latter is true even in the presence of highest experimentally allowed extragalactic magnetic fields, while the composition at lower energies can be light if a strong EGMF is present. The effect of the uncertainty in the galactic magnetic field on these results is subdominant.

We report an estimation of the injected mass composition of ultrahigh energy cosmic rays (UHECRs) at energies higher than 10 EeV. The composition is inferred from an energy-dependent sky distribution of UHECR events observed by the Telescope Array surface detector by comparing it to the Large Scale Structure of the local Universe. In the case of negligible extragalactic magnetic fields (EGMFs), the results are consistent with a relatively heavy injected composition at E∼10  EeV that becomes lighter up to E∼100  EeV, while the composition at E>100  EeV is very heavy. The latter is true even in the presence of highest experimentally allowed extragalactic magnetic fields, while the composition at lower energies can be light if a strong EGMF is present. The effect of the uncertainty in the galactic magnetic field on these results is subdominant.

We use a new method to estimate the injected mass composition of ultrahigh cosmic rays (UHECRs) at energies higher than 10 EeV. The method is based on comparison of the energy-dependent distribution of cosmic ray arrival directions as measured by the Telescope Array (TA) experiment with that calculated in a given putative model of UHECR under the assumption that sources trace the large-scale structure (LSS) of the Universe. As we report in the companion Letter, the TA data show large deflections with respect to the LSS which can be explained, assuming small extragalactic magnetic fields (EGMF), by an intermediate composition changing to a heavy one (iron) in the highest energy bin. Here we show that these results are robust to uncertainties in UHECR injection spectra, the energy scale of the experiment and galactic magnetic fields. The assumption of weak EGMF, however, strongly affects this interpretation at all but the highest energies E > 100 EeV, where the remarkable isotropy of the data implies a heavy injected composition even in the case of strong EGMF. This result also holds if UHECR sources are as rare as 2 x 10(-5) Mpc(-3), that is the conservative lower limit for the source number density.

On 11 September 2021, two small thunderstorms developed over the Telescope Array Surface Detector (TASD) that produced an unprecedented number of six downward terrestrial gamma ray flashes (TGFs) within one-hour timeframe. The TGFs occurred during the initial stage of negative cloud-to-ground flashes whose return strokes had increasingly large peak currents up to -223 kA, 147 GeV energy deposit in up to 25 1.2 km-spaced surface detectors, and intermittent bursts of gamma-rays with total durations up to 717 mu s. The analyses are based on observations recorded by the TASD network, complemented by data from a 3D lightning mapping array, broadband VHF interferometer, fast electric field change sensor, high-speed video camera, and the National Lightning Detection Network. The TGFs of the final two flashes had gamma fluences of similar or equal to 2 x 10(14) and 8x10(14), logarithmically bridging the gap between previous TASD and satellite-based detections. The observations further emphasize the similarity between upward and downward TGF varieties, suggesting a common mechanism for their production. Plain Language Summary Terrestrial Gamma-Ray Flashes (TGFs) are high-energy lightning-produced radiation events that have primarily been detected by satellites, but are increasingly detected by ground-based instrumentation. While the upward satellite-detected events exhibit extremely high numbers of gamma photons, the downward ground-based versions observed by the Telescope Array Surface Detector have, until now, displayed much weaker numbers by orders of magnitude. In this paper, we report observations of downward TGFs obtained at the large-area cosmic ray detector Telescope Array facility in west-central Utah, showing strengths midway between previously detected downward TGFs and upward satellite-detected events. The results indicate that TGFs span a wide range of fluences, a finding further supported by recent airborne observations over storm tops.

Optical emissions associated with Terrestrial Gamma-ray Flashes (TGFs) have recently become important subjects in space-based and ground-based observations since these emissions play crucial roles in understanding the generation of TGFs during thunderstorms. In this paper, we present the first time-resolved leader spectra of the optical component associated with a downward TGF. The TGF was observed at the Telescope Array Surface Detector (TASD) simultaneously with other lightning detectors, including a Lightning Mapping Array (LMA), an INTerFerometer (INTF), a Fast Antenna (FA), and a spectroscopic system. The spectroscopic system recorded leader spectra at 29,900 frames per second (33.44 $\mu$s time resolution), covering a spectral range from 400 nm to 900 nm, with 2.1 nm per pixel. The recordings of the leader spectra began 11.7 ms before the $-18$ kA return stroke and at a height of 2.37 km above the ground. These spectra reveal that optical emissions of singly-ionized nitrogen and oxygen occur between 167 $\mu$s before and 267 $\mu$s after the TGF detection, while optical emissions of neutrals (H I, 656 nm; N I, 744 nm, and O I, 777 nm) occur right at the moment of the detection. The time-dependent spectra reveal differences in the optical emissions of lightning leaders with and without downward TGFs.

We report on an observation of the difference between northern and southern skies of the ultrahigh energy cosmic ray energy spectrum with a significance of ∼8σ. We use measurements from the two largest experimentsx2014the Telescope Array observing the northern hemisphere and the Pierre Auger Observatory viewing the southern hemisphere. Since the comparison of two measurements from different observatories introduces the issue of possible systematic differences between detectors and analyses, we validate the methodology of the comparison by examining the region of the sky where the apertures of the two observatories overlap. Although the spectra differ in this region, we find that there is only a 1.8σ difference between the spectrum measurements when anisotropic regions are removed and a fiducial cut in the aperture is applied.

Various hints for anisotropies in the distribution of arrival directions of ultra-high-energy cosmic rays (UHECRs) have been reported. Still, our poor knowledge about intergalactic and Galactic magnetic fields and about the UHECR mass composition makes it non-trivial to interpret such results in terms of possible models of UHECR sources. In this work, we apply the same analyses that have been performed on the Pierre Auger Observatory and the Telescope Array UHECR data to a variety of Monte Carlo simulations generated according to many different combinations of hypotheses about the sources, composition and magnetic deflections of UHECRs. We find that only some of these models can yield results similar to those obtained with the real data.

The present white paper is submitted as part of the "Snowmass" process to help inform the long-term plans of the United States Department of Energy and the National Science Foundation for high-energy physics. It summarizes the science questions driving the Ultra-High-Energy Cosmic-Ray (UHECR) community and provides recommendations on the strategy to answer them in the next two decades.

CHASM (CHerenkov Air Shower Model) is a python package which leverages the universality of charged particles in an extensive air shower to produce a deterministic prediction of the Cherenkov light signal for a given shower profile and geometry. At sampled points throughout the domain of all shower development stages and altitudes, the angular and yield distributions of Cherenkov light have been calculated at an array of distances from a shower axis. Chasm accesses and interpolates between these distributions at runtime to produce the aggregate signal from the whole shower at user defined telescope locations. This paper gives a detailed description of the methods used to compute the Cherenkov distribution tables from universal charged particle energy and angular distributions. It also describes the workflow of CHASM itself, how a shower and axis is constructed, and how a user can customize the simulation. Finally, this paper describes how to install CHASM both from source code on GitHub and using pip (python's built in package manager) and how to use it.

Complementing the overview contribution about the whitepaper on ultra-high-energy cosmic rays (UHECR) prepared for the Snowmass community survey in the U.S. [Astroparticle Physics 149 (2023) 102819 - arXiv:2205.05845], this contribution focuses on Chapter 6, the ‘Instrumentation Roadmap’ for UHECR physics in the next decades. In addition to an increase in statistics, a higher measurement accuracy of cosmic-ray air showers is needed to answer open questions regarding the astrophysics and particle physics related to UHECR. The needed boost in exposure can be provided by space-borne fluorescence detectors with POEMMA or by huge ground arrays using a single cost-effective technique, such as the giant radio arrays envisioned with GRAND. These observatories maximizing the exposure need to be complemented by ground arrays featuring an event-by-event resolution of the rigidity of the primary particle, which is the essence of GCOS. The required high mass resolution demands the simultaneous measurement of the electromagnetic (energy and $X_\mathrm{max}$) and muonic shower components, possibly by combining layered water-Cherenkov with radio detectors and next-generation fluorescence telescopes, together with novel analysis techniques, such as neural networks. The higher accuracy for air-shower measurements is also important for UHECR particle physics because it will enable stricter tests of hadronic interaction models and will help to identify ultra-high-energy photons or neutrinos. This contribution will give an overview of the instrumentation needed for the future of UHECR physics in the context of the next generation experiments discussed in the whitepaper.

Using the Earth as a neutrino converter, tau neutrino fluxes from astrophysical point sources can be detected by tau-lepton-induced extensive air showers (EASs). Both muon neutrino and tau neutrino induced upward-going EAS signals can be detected by terrestrial, sub-orbital and satellite-based instruments. The sensitivity of these neutrino telescopes can be evaluated with the nuSpaceSim package, which includes the nuPyProp simulation package. The nuPyProp package propagates neutrinos ($\nu_\mu$, $\nu_\tau$) through the Earth to produce the corresponding charged leptons (muons and tau-leptons). We use nuPyProp to quantify the uncertainties from Earth density models, tau depolarization effects and photo-nuclear electromagnetic energy loss models in the charged lepton exit probabilities and their spectra. The largest uncertainties come from electromagnetic energy loss modeling, with as much as a 20-50% difference between the models. We compare nuPyProp results with other simulation package results.

In the study of cosmic rays, the measurement of the energy spectrum of the primaries is one of the main issues and provides fundamental information on the most energetic phenomena in the Universe. At ultrahigh energies, beyond 1018 eV, the cosmic rays are studied by the two largest observatories built so far, the Pierre Auger Observatory and the Telescope Array. Both observatories are based on a hybrid design and reported a measurement of the energy spectrum using the high duty cycle of the surface detector and the calorimetric estimation of the energy scale provided by the fluorescence detector. The differences among the reported spectra are scrutinized by a working group made by members of the Pierre Auger and Telescope Array Collaborations. The two measurements have been found well in agreement below 1019 eV while, at higher energies, they show an energy-dependent difference that is beyond the systematic uncertainties associated to the energy scale. In this contribution we review the status and perspectives of the working group activities including new studies aiming at addressing the impact on the flux measurement at the highest energies of potential biases in the estimation of the shower size.

Ultra-high-energy neutrinos serve as messengers of some of the highest energy astrophysical environments. Given that neutrinos are neutral and only interact via weak interactions, neutrinos can emerge from sources, traverse astronomical distances, and point back to their origins. Their weak interactions require large target volumes for neutrino detection. Using the Earth as a neutrino converter, terrestrial, sub-orbital, and satellite-based instruments are able to detect signals of neutrino-induced extensive air showers. In this paper, we describe the software code nuPyProp that simulates tau neutrino and muon neutrino interactions in the Earth and predicts the spectrum of the τ-leptons and muons that emerge. The nuPyProp outputs are lookup tables of charged lepton exit probabilities and energies that can be used directly or as inputs to the nuSpaceSim code designed to simulate optical and radio signals from extensive air showers induced by the emerging charged leptons. We describe the inputs to the code, demonstrate its flexibility and show selected results for τ -lepton and muon exit probabilities and energy distributions. The nuPyProp code is open source, available on github.

The present white paper is submitted as part of the "Snowmass" process to help inform the long-term plans of the United States Department of Energy and the National Science Foundation for high-energy physics. It summarizes the science questions driving the Ultra-High-Energy Cosmic-Ray (UHECR) community and provides recommendations on the strategy to answer them in the next two decades.

The Telescope Array (TA) Cosmic Ray Observatory is the largest cosmic ray detector in the northern hemisphere. TA was built to study ultra-high-energy cosmic rays (UHECRs), cosmic rays with energies above 1018 eV. TA is a hybrid detector, employing two distinct detection methods: a surface detector array and a set of fluorescence telescopes. We will present a measurement of the cosmic ray energy spectrum for energies above 1017.5 eV using only the fluorescence telescopes. A novel weather classification scheme using machine learning was used to select data parts with good weather to ensure the quality of the fluorescence data. The data from the Black Rock Mesa (BRM) and Long Ridge (LR) fluorescence telescope sites were analyzed separately in monocular mode, with the calculated fluxes combined into a single spectrum. The 10-year monocular combined cosmic ray energy spectrum is observed to be in excellent agreement with previous measurements from the northern hemisphere. We present fits of the combined spectrum to a series of broken power law models. The thrice-broken power law was observed to be the best fit considering the Poisson deviance per degrees of freedom. The three breaks suggest an additional feature of the spectrum between the previously observed Ankle feature at 1018.7 eV and the GZK suppression at 1019.8 eV.

The reconstruction of cosmic-ray-induced extensive air showers with a non-imaging Cherenkov detector array requires knowledge of the Cherenkov yield of any given air shower for a given set of shower parameters. Although air showers develop in a stochastic cascade, certain characteristics of the particles in the shower have been shown to come from universal probability distributions, a property known as shower universality. Both the energy and the angular distributions of charged particles within a shower have been parameterized. One can use these distributions to calculate the Cherenkov photon yield as an angular distribution from the Cherenkov cones of charged particles at various stages of shower development. This Cherenkov photon yield can then be tabulated for use in the reconstruction of air showers. In this work, we develop the calculation of both the Cherenkov angular distribution and Cherenkov yield per shower particle, and show how a look-up table was constructed to capture the relevant features of these distributions for general use. We compare the results of our calculations with the results of full, particle-stack, Monte Carlo simulation of the Cherenkov light produced in extensive air showers using CORSIKA-IACT. We make comparisons of both the lateral distribution of the Cherenkov photon flux amongst several detectors and of the arrival-time distribution of the Cherenkov photons in a single detector.

In this paper, we present the first high-speed video observation of a cloud-to-ground lightning flash and its associated downward-directed Terrestrial Gamma-ray Flash (TGF). The optical emission of the event was observed by a high-speed video camera running at 40,000 frames per second in conjunction with the Telescope Array Surface Detector, Lightning Mapping Array, interferometer, electric-field fast antenna, and the National Lightning Detection Network. The cloud-to-ground flash associated with the observed TGF was formed by a fast downward leader followed by a very intense return stroke peak current of -154 kA. The TGF occurred while the downward leader was below cloud base, and even when it was halfway in its propagation to ground. The suite of gamma-ray and lightning instruments, timing resolution, and source proximity offer us detailed information and therefore a unique look at the TGF phenomena.