During a winter thunderstorm on 2017 November 24, a strong burst of gamma rays with energies up to $\sim$10~MeV was detected coincident with a lightning discharge, by scintillation detectors installed at Kashiwazaki-Kariwa Nuclear Power Station at sea level in Japan. The burst had a sub-second duration, which is suggestive of photoneutron productions. The leading part of the burst was resolved into four intense gamma-ray bunches, each coincident with a low-frequency radio pulse. These bunches were separated by 0.7--1.5~ms, with a duration of $\ll$1~ms each. Thus, the present burst may be considered as a ``downward" terrestrial gamma-ray flash (TGF), which is analogous to up-going TGFs observed from space. Although the scintillation detectors were heavily saturated by these bunches, the total dose associated with them was successfully measured by ionization chambers, employed by nine monitoring posts surrounding the power plant. From this information and Monte Carlo simulations, the present downward TGF is suggested to have taken place at an altitude of 2500 $\pm$ 500~m, involving $8^{+8}_{-4} \times 10^{18}$ avalanche electrons with energies above 1~MeV. This number is comparable to those in up-going TGFs.

We report on the highest energy photons from the Crab Nebula observed by the Tibet air shower array with the underground water-Cherenkov-type muon detector array. Based on the criterion of muon number measured in an air shower, we successfully suppress 99.92% of the cosmic-ray background events with energies $E>100$ TeV. As a result, we observed 24 photon-like events with $E>100$ TeV against 5.5 background events, which corresponds to 5.6$\sigma$ statistical significance. This is the first detection of photons with $E>100$ TeV from an astrophysical source.

We examine the possible influence of Earth-directed coronal mass ejections (ECMEs) on the Sun's shadow in the 3~TeV cosmic-ray intensity observed by the Tibet-III air shower (AS) array. We confirm a clear solar-cycle variation of the intensity deficit in the Sun's shadow during ten years between 2000 and 2009. This solar-cycle variation is overall reproduced by our Monte Carlo (MC) simulations of the Sun's shadow based on the potential field model of the solar magnetic field averaged over each solar rotation period. We find, however, that the magnitude of the observed intensity deficit in the Sun's shadow is significantly less than that predicted by MC simulations, particularly during the period around solar maximum when a significant number of ECMEs is recorded. The $\chi^2$ tests of the agreement between the observations and the MC simulations show that the difference is larger during the periods when the ECMEs occur, and the difference is reduced if the periods of ECMEs are excluded from the analysis. This suggests the first experimental evidence of the ECMEs affecting the Sun's shadow observed in the 3~TeV cosmic-ray intensity.

We analyze the Sun's shadow observed with the Tibet-III air shower array and find that the shadow's center deviates northward (southward) from the optical solar disc center in the "Away" ("Toward") IMF sector. By comparing with numerical simulations based on the solar magnetic field model, we find that the average IMF strength in the "Away" ("Toward") sector is $1.54 \pm 0.21_{\rm stat} \pm 0.20_{\rm syst}$ ($1.62 \pm 0.15_{\rm stat} \pm 0.22_{\rm syst}$) times larger than the model prediction. These demonstrate that the observed Sun's shadow is a useful tool for the quantitative evaluation of the average solar magnetic field.

Lightning and thundercloud are the most dramatic natural particle accelerators on the Earth. Relativistic electrons accelerated by electric fields therein emit bremsstrahlung gamma rays, which have been detected at ground observations, by airborne detectors, and as terrestrial gamma-ray flashes (TGFs) from space. The energy of the gamma rays is sufficiently high to potentially invoke atmospheric photonuclear reactions 14N(gamma, n)13N, which would produce neutrons and eventually positrons via beta-plus decay of generated unstable radioactive isotopes, especially 13N. However, no clear observational evidence for the reaction has been reported to date. Here we report the first detection of neutron and positron signals from lightning with a ground observation. During a thunderstorm on 6 February 2017 in Japan, a TGF-like intense flash (within 1 ms) was detected at our monitoring sites 0.5-1.7 km away from the lightning. The subsequent initial burst quickly subsided with an exponential decay constant of 40-60 ms, followed by a prolonged line emission at about 0.511 megaelectronvolt (MeV), lasting for a minute. The observed decay timescale and spectral cutoff at about 10 MeV of the initial emission are well explained with de-excitation gamma rays from the nuclei excited by neutron capture. The centre energy of the prolonged line emission corresponds to the electron-positron annihilation, and hence is the conclusive indication of positrons produced after the lightning. Our detection of neutrons and positrons is unequivocal evidence that natural lightning triggers photonuclear reactions. No other natural event on the Earth is known to trigger photonuclear reactions. This discovery places lightning as only the second known natural channel on the Earth after the atmospheric cosmic-ray interaction, in which isotopes, such as 13C, 14C, and 15N, are produced.

At the 33rd ICRC, we reported the possible detection of solar gamma rays by a ground level detector and later re-examined this event. On March 7, 2011, the solar neutron telescope (SNT) located at Mt. Sierra Negra, Mexico (4,600 m) observed enhancements of the counting rate from 19:49 to 20:02 UT and from 20:50 to 21:01 UT. The statistical significance was 9.7sigma and 8.5sigma, respectively. This paper discusses the possibility of using this mountain detector to detect solar gamma rays. In association with this event, the solar neutron detector SEDA-FIB onboard the International Space Station has also detected solar neutrons with a statistical significance of 7.5sigma. The FERMI-LAT detector also observed high-energy gamma rays from this flare with a statistical significance of 6.7sigma. We thus attempted to make a unified model to explain this data. In this paper, we report on another candidate for solar gamma rays detected on September 25th, 2011 by the SNT located in Tibet (4,300 m) from 04:37 to 04:47 UT with a statistical significance of 8.0sigma (by the Li-Ma method).

Observations of strong flux of low-energy neutrons were made by $^{3}\mathrm{He}$ counters during thunderstorms [Gurevich et al (Phys. Rev. Lett. 108, 125001, 2012)]. How the unprecedented enhancements were produced remains elusive. To better elucidate the mechanism, a simulation study of surrounding material impacts on measurement by $^{3}\mathrm{He}$ counters was performed with GEANT4. It was found that unlike previously thought, a $^3\mathrm{He}$ counter had a small sensitivity to high-energy gamma rays because of inelastic interaction with its cathode-tube materials (Al or stainless steel). A $^{3}\mathrm{He}$ counter with the intrinsic small sensitivity, if surrounded by thick materials, would largely detect thunderstorm-related gamma rays rather than those neutrons produced via photonuclear reaction in the atmosphere. On the other hand, the counter, if surrounded by thin materials and located away from a gamma-ray source, would observe neutron signals with little gamma-ray contamination. Compared with the Gurevich measurement, the present work allows us to deduce that the enhancements are attributable to gamma rays, if their observatory was very close to or inside a gamma-ray emitting region in thunderclouds.

We report on a clear solar-cycle variation of the Sun's shadow in the 10 TeV cosmic-ray flux observed by the Tibet air shower array during a full solar cycle from 1996 to 2009. In order to clarify the physical implications of the observed solar cycle variation, we develop numerical simulations of the Sun's shadow, using the Potential Field Source Surface (PFSS) model and the Current Sheet Source Surface (CSSS) model for the coronal magnetic field. We find that the intensity deficit in the simulated Sun's shadow is very sensitive to the coronal magnetic field structure, and the observed variation of the Sun's shadow is better reproduced by the CSSS model. This is the first successful attempt to evaluate the coronal magnetic field models by using the Sun's shadow observed in the TeV cosmic-ray flux.

We report the first observation of 3$-$30 MeV prolonged gamma-ray emission that was abruptly terminated by lightning. The gamma-ray detection was made during winter thunderstorms on December 30, 2010 by the Gamma-Ray Observation of Winter THunderclouds (GROWTH) experiment carried out in a coastal area along the Sea of Japan. The gamma-ray flux lasted for less than 3 min, continuously hardening closer to the lightning occurrence. The hardening at energies of 3$-$10 MeV energies was most prominent. The gamma-ray flux abruptly ceased less than 800 ms before the lightning flash that occurred over 5 km away from the experimental site. In addition, we observed a clear difference in the duration of the 3$-$10 MeV gamma rays and those $>$10 MeV, suggesting that the area of $>$10 MeV gamma-ray emission is considerably smaller than that of the lower-energy gamma rays. This work may give a manifestation that a local region emitting prolonged gamma rays connects with a distant region to initiate lightning.

A new hybrid experiment has been started by AS{\gamma} experiment at Tibet, China, since August 2011, which consists of a low threshold burst-detector-grid (YAC-II, Yangbajing Air shower Core array), the Tibet air-shower array (Tibet-III) and a large underground water Cherenkov muon detector (MD). In this paper, the capability of the measurement of the chemical components (proton, helium and iron) with use of the (Tibet-III+YAC-II) is investigated by means of an extensive Monte Carlo simulation in which the secondary particles are propagated through the (Tibet-III+YAC-II) array and an artificial neural network (ANN) method is applied for the primary mass separation. Our simulation shows that the new installation is powerful to study the chemical compositions, in particular, to obtain the primary energy spectrum of the major component at the knee.

The Gamma-Ray Observation of Winter THunderclouds (GROWTH) experiment, consisting of two radiation-detection subsystems, has been operating since 2006 on the premises of Kashiwazaki-Kariwa nuclear power plant located at the coastal area of Japan Sea. By 2010 February, GROWTH detected 7 long-duration $\gamma$-rays emissions associated with winter thunderstorms. Of them, two events, obtained on 2007 December 13 and 2008 December 25, are reported.On both occasions, all inorganic scintillators (NaI, CsI, and BGO) of the two subsystems detected significant gamma-ray signals lasting for >1 minute. Neither of these two events were associated with any lightning. In both cases, the gamma-ray energy spectra extend to 10 MeV, suggesting that the detected gamma-rays are produced by relativistic electrons via bremsstrahlung. Assuming that the initial photon spectrum at the source is expressed by a power-law function,the observed photons can be interpreted as being radiated from a source located at a distance of 290-560 m for the 2007 event and 110-690 m for the 2008 one, both at 90% confidence level.Employing these photon spectra, the number of relativistic electrons is estimated as 10^9 - 10^{11}. The estimation generally agrees with those calculated based on the relativistic runaway electron avalanche model. A GROWTH photon spectrum, summed over 3 individual events including the present two events and another reported previously, has similar features including a cut-off energy, to an averaged spectrum of terrestrial gamma-ray flashes.

We analyze the large-scale two-dimensional sidereal anisotropy of multi-TeV cosmic rays by Tibet Air Shower Array, with the data taken from 1999 November to 2008 December. To explore temporal variations of the anisotropy, the data set is divided into nine intervals, each in a time span of about one year. The sidereal anisotropy of magnitude about 0.1% appears fairly stable from year to year over the entire observation period of nine years. This indicates that the anisotropy of TeV Galactic cosmic rays remains insensitive to solar activities since the observation period covers more than a half of the 23rd solar cycle.

Using the Tibet-III air shower array, we search for TeV gamma-rays from 27 potential Galactic sources in the early list of bright sources obtained by the Fermi Large Area Telescope at energies above 100 MeV. Among them, we observe 7 sources instead of the expected 0.61 sources at a significance of 2 sigma or more excess. The chance probability from Poisson statistics would be estimated to be 3.8 x 10^-6. If the excess distribution observed by the Tibet-III array has a density gradient toward the Galactic plane, the expected number of sources may be enhanced in chance association. Then, the chance probability rises slightly, to 1.2 x 10^-5, based on a simple Monte Carlo simulation. These low chance probabilities clearly show that the Fermi bright Galactic sources have statistically significant correlations with TeV gamma-ray excesses. We also find that all 7 sources are associated with pulsars, and 6 of them are coincident with sources detected by the Milagro experiment at a significance of 3 sigma or more at the representative energy of 35 TeV. The significance maps observed by the Tibet-III air shower array around the Fermi sources, which are coincident with the Milagro >=3sigma sources, are consistent with the Milagro observations. This is the first result of the northern sky survey of the Fermi bright Galactic sources in the TeV region.

We develop a model anisotropy best-fitting to the two-dimensional sky-map of multi-TeV galactic cosmic ray (GCR) intensity observed with the Tibet III air shower (AS) array. By incorporating a pair of intensity excesses in the hydrogen deflection plane (HDP) suggested by Gurnett et al., together with the uni-directional and bi-directional flows for reproducing the observed global feature, this model successfully reproduces the observed sky-map including the "skewed" feature of the excess intensity from the heliotail direction, whose physical origin has long remained unknown. These additional excesses are modeled by a pair of the northern and southern Gaussian distributions, each placed ~50 degree away from the heliotail direction. The amplitude of the southern excess is as large as ~0.2 %, more than twice the amplitude of the northern excess. This implies that the Tibet AS experiment discovered for the first time a clear evidence of the significant modulation of GCR intensity in the heliotail and the asymmetric heliosphere.

Aiming to observe cosmic gamma rays in the 10 - 1000 TeV energy region, we propose a 10000 m^2 underground water-Cherenkov muon-detector (MD) array that operates in conjunction with the Tibet air-shower (AS) array. Significant improvement is expected in the sensitivity of the Tibet AS array towards celestial gamma-ray signals above 10 TeV by utilizing the fact that gamma-ray-induced air showers contain far fewer muons compared with cosmic-ray-induced ones. We carried out detailed Monte Carlo simulations to assess the attainable sensitivity of the Tibet AS+MD array towards celestial TeV gamma-ray signals. Based on the simulation results, the Tibet AS+MD array will be able to reject 99.99% of background events at 100 TeV, with 83% of gamma-ray events remaining. The sensitivity of the Tibet AS+MD array will be ~20 times better than that of the present Tibet AS array around 20 - 100 TeV. The Tibet AS+MD array will measure the directions of the celestial TeV gamma-ray sources and the cutoffs of their energy spectra. Furthermore, the Tibet AS+MD array, along with imaging atmospheric Cherenkov telescopes as well as the Fermi Gamma-ray Space Telescope and X-ray satellites such as Suzaku and MAXI, will make multiwavelength observations and conduct morphological studies on sources in the quest for evidence of the hadronic nature of the cosmic-ray acceleration mechanism.

During thunderstorms on 2008 September 20, a simultaneous detection of gamma rays and electrons was made at a mountain observatory in Japan located 2770 m above sea level. Both emissions, lasting 90 seconds, were associated with thunderclouds rather than lightning. The photon spectrum, extending to 10 MeV, can be interpreted as consisting of bremsstrahlung gamma rays arriving from a source which is 60 - 130 m in distance at 90% confidence level. The observed electrons are likely to be dominated by a primary population escaping from an acceleration region in the clouds.

The amplitude of the Compton-Getting (CG) anisotropy contains the power-law index of the cosmic-ray energy spectrum. Based on this relation and using the Tibet air-shower array data, we measure the cosmic-ray spectral index to be $-3.03 \pm 0.55_{stat} \pm < 0.62_{syst}$ between 6 TeV and 40 TeV, consistent with $-$2.7 from direct energy spectrum measurements. Potentially, this CG anisotropy analysis can be utilized to confirm the astrophysical origin of the ``knee'' against models for non-standard hadronic interactions in the atmosphere.

The Tibet air shower array, which has an effective area of 37,000 square meters and is located at 4300 m in altitude, has been observing air showers induced by cosmic rays with energies above a few TeV. We are planning to add a large muon detector array to it for the purpose of increasing its sensitivity to cosmic gamma rays in the 100 TeV (10 - 1000 TeV) energy region by discriminating them from cosmic-ray hadrons. We report on the possibility of detection of gamma rays in the 100 TeV energy region in our field of view, based on the improved sensitivity of our air shower array deduced from the full Monte Carlo simulation.

The Tibet air shower array, which has an effective area of 37,000 square meters and is located at 4300 m in altitude, has been observing air showers induced by cosmic rays with energies above a few TeV. We have a plan to add a large muon detector array to it for the purpose of increasing its sensitivity to cosmic gamma rays in the 100 TeV energy region by discriminating them from cosmic-ray hadrons. We have deduced the attainable sensitivity of the muon detector array using our Monte Carlo simulation. We report here on the detailed procedure of our Monte Carlo simulation.

A report is made on a comprehensive observation of a burst-like $\gamma$-ray emission from thunderclouds on the Sea of Japan, during strong thunderstorms on 2007 January 6. The detected emission, lasting for $\sim$40 seconds, preceded cloud-to-ground lightning discharges. The burst spectrum, extending to 10 MeV, can be interpreted as consisting of bremsstrahlung photons originating from relativistic electrons. This ground-based observation provides first clear evidence that strong electric fields in thunderclouds can continuously accelerate electrons beyond 10 MeV prior to lightning discharges.

The purpose of this work is to search the Yangbajing neutron monitor data obtained between 1998 October and 2000 June for solar neutrons associated with solar flares. Using the onset times of 166 BATSE-detected flares with the GOES peak flux (1 -- 8 \AA) higher than $1.0 \times 10^{-5}$ $\mathrm{Wm^{-2}}$, we prepare for each flare a light curve of the Yangbajing neutron monitor, spanning $\pm$ 1.5 hours from the BATSE onset time. Based on the light curves, a systematic search for solar neutrons in energies above 100 MeV from the 166 flares was performed. No statistically significant signals due to solar neutrons were found in the present work. Therefore, we put upper limits on the $>$ 100 MeV solar-neutron flux for 18 events consisting of 2 X and 16 M class flares. The calculation assumed a power-law shaped neutron energy spectrum and three types of neutron emission profiles at the Sun. Compared with the other positive neutron detections associated with X-class flares, typical 95% confidence level upper limits for the two X-class flares are found to be comparable to the lowest and second lowest neutron fluxes at the top of the atmosphere.In addition, the upper limits for M-class flares scatter in the range of $10^{-2}$ to 1 neutrons $\mathrm{cm^{-2}s^{-1}}$. This provides the first upper limits on the solar-neutron flux from M-class solar flares, using space observatories as well as ground-based neutron monitors.

Solar neutrons have been detected using the neutron monitor located at Mt. Chacaltaya, Bolivia, in association with a large solar flare on November 24, 2000. This is the first detection of solar neutrons by the neutron monitor that have been reported so far in solar cycle 23. The statistical significance of the detection is 5.5 sigma. In this flare, the intense emission of hard X-rays and gamma-rays has been observed by the Yohkoh Hard X-ray Telescope (HXT) and Gamma Ray Spectrometer (GRS), respectively. The production time of solar neutrons is better correlated with those of hard X-rays and gamma-rays than with the production time of soft X-rays. The observations of the solar neutrons on the ground have been limited to solar flares with soft X-ray class greater than X8 in former solar cycles. In this cycle, however, neutrons were detected associated with an X2.3 solar flare on November 24, 2000. This is the first report of the detection of solar neutrons on the ground associated with a solar flare with its X-ray class smaller than X8.