The Telescope Array (TA) observatory utilizes fluorescence detectors and surface detectors (SDs) to observe air showers produced by ultra high energy cosmic rays in Earth’s atmosphere. Cosmic-ray events observed in this way are termed hybrid data. The depth of air shower maximum is related to the mass of the primary particle that generates the shower. This paper reports on shower maxima data collected over 8.5 yr using the Black Rock Mesa and Long Ridge fluorescence detectors in conjunction with the array of SDs. We compare the means and standard deviations of the observed X max distributions with Monte Carlo X max distributions of unmixed protons, helium, nitrogen, and iron, all generated using the QGSJet II-04 hadronic model. We also perform an unbinned maximum likelihood test of the observed data, which is subjected to variable systematic shifting of the data X max distributions to allow us to test the full distributions, and compare them to the Monte Carlo to see which elements are not compatible with the observed data. For all energy bins, QGSJet II-04 protons are found to be compatible with TA hybrid data at the 95% confidence level after some systematic X max shifting of the data. Three other QGSJet II-04 elements are found to be compatible using the same test procedure in an energy range limited to the highest energies where data statistics are sparse.

In this paper we report the first close, high‐resolution observations of downward‐directed terrestrial gamma‐ray flashes (TGFs) detected by the large‐area Telescope Array cosmic ray observatory, obtained in conjunction with broadband VHF interferometer and fast electric field change measurements of the parent discharge. The results show that the TGFs occur during strong initial breakdown pulses (IBPs) in the first few milliseconds of negative cloud‐to‐ground and low‐altitude intracloud flashes and that the IBPs are produced by a newly identified streamer‐based discharge process called fast negative breakdown. The observations indicate the relativistic runaway electron avalanches (RREAs) responsible for producing the TGFs are initiated by embedded spark‐like transient conducting events (TCEs) within the fast streamer system and potentially also by individual fast streamers themselves. The TCEs are inferred to be the cause of impulsive sub‐pulses that are characteristic features of classic IBP sferics. Additional development of the avalanches would be facilitated by the enhanced electric field ahead of the advancing front of the fast negative breakdown. In addition to showing the nature of IBPs and their enigmatic sub‐pulses, the observations also provide a possible explanation for the unsolved question of how the streamer to leader transition occurs during the initial negative breakdown, namely, as a result of strong currents flowing in the final stage of successive IBPs, extending backward through both the IBP itself and the negative streamer breakdown preceding the IBP.

We report on a measurement of the cosmic ray composition by the Telescope Array Low-Energy Extension (TALE) air fluorescence detector (FD). By making use of the Cherenkov light signal in addition to air fluorescence light from cosmic ray (CR) induced extensive air showers, the TALE FD can measure the properties of the cosmic rays with energies as low as ∼ 2 PeV and exceeding 1 EeV. In this paper, we present results on the measurement of X max distributions of showers observed over this energy range. Data collected over a period of ∼ 4 years was analyzed for this study. The resulting X max distributions are compared to the Monte Carlo (MC) simulated data distributions for primary cosmic rays with varying composition and a 4-component fit is performed. The comparison and fit are performed for energy bins, of width 0.1 or 0.2 in log 10 (E/eV), spanning the full range of the measured energies. We also examine the mean X max value as a function of energy for cosmic rays with energies greater than 10 15.8 eV. Below 10 17.3 eV, the slope of the mean X max as a function of energy (the elongation rate) for the data is significantly smaller than that of all elements in the models, indicating that the composition is becoming heavier with energy in this energy range. This is consistent with a rigidity-dependent cutoff of events from galactic sources. Finally, an increase in the X max elongation rate is observed at energies just above 10 17 eV indicating another change in the cosmic rays composition.

The Pierre Auger Collaboration (Auger) recently reported a correlation between the arrival directions of cosmic rays with energies above 39 EeV and the flux pattern of 23 nearby starburst galaxies (SBGs). In this Letter, we tested the same hypothesis using cosmic rays detected by the Telescope Array experiment (TA) in the 9 yr period from 2008 May to 2017 May. Unlike the Auger analysis, we did not optimize the parameter values but kept them fixed to the best-fit values found by Auger, namely 9.7% for the anisotropic fraction of cosmic rays assumed to originate from the SBGs in the list and 12.°9 for the angular scale of the correlations. The energy threshold that we adopted is 43 EeV, corresponding to 39 EeV in Auger when taking into account the energy-scale difference between two experiments. We find that the TA data is compatible with isotropy to within 1.1σ and with the Auger result to within 1.4σ, meaning that it is not capable to discriminate between these two hypotheses.

One of the uncertainties in interpretation of ultra-high energy cosmic ray (UHECR) data comes from the hadronic interaction models used for air shower Monte Carlo (MC) simulations. The number of muons observed at the ground from UHECR-induced air showers is expected to depend upon the hadronic interaction model. One may therefore test the hadronic interaction models by comparing the measured number of muons with the MC prediction. In this paper, we present the results of studies of muon densities in UHE extensive air showers obtained by analyzing the signal of surface detector stations which should have high muon purity. The muon purity of a station will depend on both the inclination of the shower and the relative position of the station. In 7 years' data from the Telescope Array experiment, we find that the number of particles observed for signals with an expected muon purity of ∼65% at a lateral distance of 2000 m from the shower core is 1.72 ± 0.10(stat.) ± 0.37(syst.) times larger than the MC prediction value using the QGSJET II-03 model for proton-induced showers. A similar effect is also seen in comparisons with other hadronic models such as QGSJET II-04, which shows a 1.67 ± 0.10 ± 0.36 excess. We also studied the dependence of these excesses on lateral distances and found a slower decrease of the lateral distribution of muons in the data as compared to the MC, causing larger discrepancy at larger lateral distances.