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.

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.

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.

Cosmic rays are energetic charged particles from extraterrestrial sources, with the highest-energy events thought to come from extragalactic sources. Their arrival is infrequent, so detection requires instruments with large collecting areas. In this work, we report the detection of an extremely energetic particle recorded by the surface detector array of the Telescope Array experiment. We calculate the particle's energy as [Formula: see text] (~40 joules). Its arrival direction points back to a void in the large-scale structure of the Universe. Possible explanations include a large deflection by the foreground magnetic field, an unidentified source in the local extragalactic neighborhood, or an incomplete knowledge of particle physics.

We present the first observation by the Telescope Array Surface Detector (TASD) of the effect of thunderstorms on the development of cosmic ray single count rate intensity over a 700 km$^{2}$ area. Observations of variations in the secondary low-energy cosmic ray counting rate, using the TASD, allow us to study the electric field inside thunderstorms, on a large scale, as it progresses on top of the 700 km$^{2}$ detector, without dealing with the limitation of narrow exposure in time and space using balloons and aircraft detectors. In this work, variations in the cosmic ray intensity (single count rate) using the TASD, were studied and found to be on average at the $\sim(0.5-1)\%$ and up to 2\% level. These observations were found to be both in excess and in deficit. They were also found to be correlated with lightning in addition to thunderstorms. These variations lasted for tens of minutes; their footprint on the ground ranged from 6 to 24 km in diameter and moved in the same direction as the thunderstorm. With the use of simple electric field models inside the cloud and between cloud to ground, the observed variations in the cosmic ray single count rate were recreated using CORSIKA simulations. Depending on the electric field model used and the direction of the electric field in that model, the electric field magnitude that reproduces the observed low-energy cosmic ray single count rate variations was found to be approximately between 0.2-0.4 GV. This in turn allows us to get a reasonable insight on the electric field and its effect on cosmic ray air showers inside thunderstorms.

For several decades, the origin of ultra-high-energy cosmic rays (UHECRs) has been an unsolved question of high-energy astrophysics. One approach for solving this puzzle is to correlate UHECRs with high-energy neutrinos, since neutrinos are a direct probe of hadronic interactions of cosmic rays and are not deflected by magnetic fields. In this paper, we present three different approaches for correlating the arrival directions of neutrinos with the arrival directions of UHECRs. The neutrino data is provided by the IceCube Neutrino Observatory and ANTARES, while the UHECR data with energies above $\sim$50 EeV is provided by the Pierre Auger Observatory and the Telescope Array. All experiments provide increased statistics and improved reconstructions with respect to our previous results reported in 2015. The first analysis uses a high-statistics neutrino sample optimized for point-source searches to search for excesses of neutrinos clustering in the vicinity of UHECR directions. The second analysis searches for an excess of UHECRs in the direction of the highest-energy neutrinos. The third analysis searches for an excess of pairs of UHECRs and highest-energy neutrinos on different angular scales. None of the analyses has found a significant excess, and previously reported over-fluctuations are reduced in significance. Based on these results, we further constrain the neutrino flux spatially correlated with UHECRs.

The Telescope Array Collaboration has observed an excess of events with E ≥ 10^19.4 eV in the data which is centered at (RA, dec) = (19^∘, 35^∘). This is near the center of the Perseus-Pisces supercluster (PPSC). The PPSC is about 70 Mpc distant and is the closest supercluster in the Northern Hemisphere (other than the Virgo supercluster of which we are a part). A Li-Ma oversampling analysis with 20^∘-radius circles indicates an excess in the arrival direction of events with a local significance of about 4 standard deviations. The probability of having such excess close to the PPSC by chance is estimated to be 3.5 standard deviations. This result indicates that a cosmic ray source likely exists in that supercluster.

Telescope Array (TA) is the largest experiment in the Northern Hemisphere studying ultrahigh energy cosmic rays. TA measurements of the cosmic ray spectrum using the surface detector have the best statistical power in the experiment, and observe the ankle of the spectrum and the high energy cutoff. When the data are divided into two declination bands, above and below 24.8 degrees, the cutoff appears at 10^19.64 ± 0.04 (10^19.84 ± 0.02) eV in the lower (higher) band, an energy difference of 58%. The global significance of the difference is 4.3 standard deviations. The lack of an instrumental cause of this difference implies it is astrophysical in nature.

The energy of the ultrahigh energy spectral cutoff was measured, integrating over the northern hemisphere sky, by the Telescope Array (TA) collaboration, to be $10^{19.78 \pm 0.06}$ eV, in agreement with the High Resolution Fly's Eye (HiRes) experiment, whereas the Pierre Auger experiment, integrating over the southern hemisphere sky, measured the cutoff to be at 10$^{19.62 \pm 0.02}$ eV. An 11% energy scale difference between the TA and Auger does not account for this difference. However, in comparing the spectra of the Telescope Array and Pierre Auger experiments in the band of declination common to both experiments ($-15.7^{\circ} < \delta < 24.8^{\circ}$) we have found agreement in the energy of the spectral cutoff. While the Auger result is essentially unchanged, the TA cutoff energy has changed to $10^{19.59 \pm 0.06}$ eV. In this paper we argue that this is an astrophysical effect.