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.

We report the analysis of the $10-1000$ TeV large-scale sidereal anisotropy of Galactic cosmic rays (GCRs) with the data collected by the Tibet Air Shower Array from October, 1995 to February, 2010. In this analysis, we improve the energy estimate and extend the declination range down to $-30^{\circ}$. We find that the anisotropy maps above 100 TeV are distinct from that at multi-TeV band. The so-called "tail-in" and "loss-cone" features identified at low energies get less significant and a new component appears at $\sim100$ TeV. The spatial distribution of the GCR intensity with an excess (7.2$\sigma$ pre-trial, 5.2$\sigma$ post-trial) and a deficit ($-5.8\sigma$ pre-trial) are observed in the 300 TeV anisotropy map, in a good agreement with IceCube's results at 400 TeV. Combining the Tibet results in the northern sky with IceCube's results in the southern sky, we establish a full-sky picture of the anisotropy in hundreds of TeV band. We further find that the amplitude of the first order anisotropy increases sharply above $\sim100$ TeV, indicating a new component of the anisotropy. All these results may shed new light on understanding the origin and propagation of GCRs.

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.

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.

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.

The Tibet-III air shower array, consisting of 533 scintillation detectors, has been operating successfully at Yangbajing in Tibet, China since 1999. Using the dataset collected by this array from 1999 November through 2005 November, we obtained the energy spectrum of $\gamma$-rays from the Crab Nebula, expressed by a power law as $(dJ/dE) = (2.09\pm0.32)\times10^{-12} (E/{\rm 3 TeV})^{-2.96\pm0.14} {\rm cm}^{-2} {\rm s}^{-1} {\rm TeV}^{-1}$ in the energy range of 1.7 to 40 TeV. This result is consistent with other independent $\gamma$-ray observations by imaging air Cherenkov telescopes. In this paper, we carefully checked and tuned the performance of the Tibet-III array using data on the moon's shadow in comparison with a detailed Monte Carlo simulation. The shadow is shifted to the west of the moon's apparent position as an effect of the geomagnetic field, although the extent of this displacement depends on the primary energy positively charged cosmic rays. This finding enables us to estimate the systematic error in determining the primary energy from its shower size. This error is estimated to be less than $\pm$12% in our experiment. This energy scale estimation is the first attempt among cosmic-ray experiments at ground level. The systematic pointing error is also estimated to be smaller than $0\fdg011$. The deficit rate and position of the moon's shadow are shown to be very stable within a statistical error of $\pm$6% year by year. This guarantees the long-term stability of point-like source observation with the Tibet-III array. These systematic errors are adequately taken into account in our study of the Crab Nebula.

We have already reported the first result on the all-particle spectrum around the knee region based on data from 2000 November to 2001 October observed by the Tibet-III air-shower array. In this paper, we present an updated result using data set collected in the period from 2000 November through 2004 October in a wide range over 3 decades between $10^{14}$ eV and $10^{17}$ eV, in which the position of the knee is clearly seen at around 4 PeV. The spectral index is -2.68 $\pm$ 0.02(stat.) below 1PeV, while it is -3.12 $\pm$ 0.01(stat.) above 4 PeV in the case of QGSJET+HD model, and various systematic errors are under study now.

We present an updated all-particle energy spectrum of primary cosmic rays in a wide range from 10^14 eV to 10^17 eV using 5.5 times 10^7 events collected in the period from 2000 November through 2004 October by the Tibet-III air-shower array located at 4300 m above sea level (atmospheric depth of 606 g/cm^2). The size spectrum exhibits a sharp knee at a corresponding primary energy around 4 PeV. This work uses increased statistics and new simulation calculations for the analysis. We performed extensive Monte Carlo calculations and discuss the model dependences involved in the final result assuming interaction models of QGSJET01c and SIBYLL2.1 and primary composition models of heavy dominant (HD) and proton dominant (PD) ones. Pure proton and pure iron primary models are also examined as extreme cases. The detector simulation was also made to improve the accuracy of determining the size of the air showers and the energy of the primary particle. We confirmed that the all-particle energy spectra obtained under various plausible model parameters are not significantly different from each other as expected from the characteristics of the experiment at the high altitude, where the air showers of the primary energy around the knee reaches near maximum development and their features are dominated by electromagnetic components leading to the weak dependence on the interaction model or the primary mass. This is the highest-statistical and the best systematics-controlled measurement covering the widest energy range around the knee energy region.

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.

We searched for pulsed gamma-ray emissions from the Crab pulsar using data of the Tibet-III air shower array from November 1999 through November 2005. No evidence for the pulsed emissions was found in our analysis. Upper limits at different energies were calculated for a $3 \sigma$ confidence level in the energy range of multi-TeV to several hundred TeV.

We have observed the shadowing of galactic cosmic ray flux in the direction of the moon, the so-called moon shadow, using the Tibet-III air shower array operating at Yangbajing (4300 m a.s.l.) in Tibet since 1999. Almost all cosmic rays are positively charged; for that reason, they are bent by the geomagnetic field, thereby shifting the moon shadow westward. The cosmic rays will also produce an additional shadow in the eastward direction of the moon if cosmic rays contain negatively charged particles, such as antiprotons, with some fraction. We selected 1.5 x10^{10} air shower events with energy beyond about 3 TeV from the dataset observed by the Tibet-III air shower array and detected the moon shadow at $\sim 40 \sigma$ level. The center of the moon was detected in the direction away from the apparent center of the moon by 0.23$^\circ$ to the west. Based on these data and a full Monte Carlo simulation, we searched for the existence of the shadow produced by antiprotons at the multi-TeV energy region. No evidence of the existence of antiprotons was found in this energy region. We obtained the 90% confidence level upper limit of the flux ratio of antiprotons to protons as 7% at multi-TeV energies.

The flux upper limits of the diffuse gamma rays, from the inner and outer Galactic planes, are revised by factors of 4.0$\sim$3.7 for mode energies 3$\sim$10 TeV, respectively, by using the simulation results of the effective area ratios for gamma-ray induced showers and cosmic-ray induced ones in the Tibet air shower array. In our previous work, (Amenomori et al., ApJ, 580, 887, 2002) the flux upper limits were deduced only from the flux ratio of air showers generated by gamma rays versus cosmic rays. The details of the simulation are given in the paper (Amenomori et al., Advances in Space Research, 37, 1932, 2006). The present result using the same data as in ApJ suggests that the spectral index of source electrons is steeper than 2.2 and 2.1 for the inner and outer Galactic planes, respectively.

We propose to build a large water-Cherenkov-type muon-detector array (Tibet MD array) around the 37,000 m$^{2}$ Tibet air shower array (Tibet AS array) already constructed at 4,300 m above sea level in Tibet, China. Each muon detector is a waterproof concrete pool, 6 m wide $\times$ 6 m long $\times$ 1.5 m deep in size, equipped with a 20 inch-in-diameter PMT. The Tibet MD array consists of 240 muon detectors set up 2.5 m underground. Its total effective area will be 8,640 m$^{2}$ for muon detection. The Tibet MD array will significantly improve gamma-ray sensitivity of the Tibet AS array in the 100 TeV region (10-1000 TeV) by means of gamma/hadron separation based on counting the number of muons accompanying an air shower. The Tibet AS+MD array will have the sensitivity to gamma rays in the 100 TeV region by an order of magnitude better than any other previous existing detectors in the world.

The intensity of Galactic cosmic rays is nearly isotropic because of the influence of magnetic fields in the Milky Way. Here, we present two-dimensional high-precision anisotropy measurement for energies from a few to several hundred teraelectronvolts (TeV), using the large data sample of the Tibet Air Shower Arrays. Besides revealing finer details of the known anisotropies, a new component of Galactic cosmic ray anisotropy in sidereal time is uncovered around the Cygnus region direction. For cosmic-ray energies up to a few hundred TeV, all components of anisotropies fade away, showing a corotation of Galactic cosmic rays with the local Galactic magnetic environment. These results have broad implications for a comprehensive understanding of cosmic rays, supernovae, magnetic fields, and heliospheric and Galactic dynamic environments.

We searched for steady PeV gamma-ray emission from the Monogem ring region with the Tibet air shower array from 1997 February to 2004 October. No evidence for statistically significant gamma-ray signals was found in a region 111$\degr$ $\leq$ R.A. $<$ 114$\degr$, 12$\fdg$5 $\leq$ decl. $<$ 15$\fdg$5 in the Monogem ring where the MAKET-ANI experiment recently claimed a positive detection of PeV high-energy cosmic radiation, although our flux sensitivity is approximately 10 times better than MAKET-ANI's. We set the most stringent integral flux upper limit at a 99% confidence level of 4.0 $\times$ 10$^{-12}$ cm$^{-2}$ s$^{-1}$ sr$^{-1}$ above 1 PeV on diffuse gamma rays extended in the 3$^{\circ}$ $\times$ 3$^{\circ}$ region.

We obtained new upper limits on the diffuse gamma rays from the inner Galactic (IG) and outer Galactic (OG) planes in 3.10 TeV region, using the Tibet air shower data and new Monte Carlo simulation results. A difference of the effective area of the air-shower array for observing gamma rays and cosmic rays was carefully taken into account in this analysis, resulting in that the flux upper limits of the diffuse TeV gamma rays were reduced by factors of 4.0.3.7 for 3.10 TeV than those in our previous results (Amenomori, M., Ayabe, S., Cui, S.W., et al. Observation of multi-TeV diffuse gamma rays from the Galactic plane with the Tibet air shower array. Astrophys. J. 580, 887.895, 2002.). This new result suggests that the inverse power index of the energy spectrum of source electrons responsible for generating diffuse TeV gamma rays through inverse Compton effect should be steeper than 2.2 and 2.1 for IG and OG planes, respectively, with 99%C.L.

A hybrid experiment consisting of emulsion chambers, burst detectors and the Tibet II air-shower array was carried out at Yangbajing (4,300 m a.s.l., 606 g/cm$^2$) in Tibet to obtain the energy spectra of primary protons and heliums. From three-year operation, these energy spectra are deduced between $10^{15}$ and $10^{16}$ eV by triggering the air showers associated with a high energy core and using a neural network method in the primary mass separation. The proton spectrum can be expressed by a single power-law function with a differential index of $-3.01 \pm 0.11$ and $-3.05 \pm 0.12$ based on the QGSJET+HD and SIBYLL+HD models, respectively, which are steeper than that extrapolated from the direct observations of $-2.74 \pm 0.01$ in the energy range below $10^{14}$ eV. The absolute fluxes of protons and heliums are derived within 30% systematic errors depending on the hadronic interaction models used in Monte Carlo simulation. The result of our experiment suggests that the main component responsible for the change of the power index of the all-particle spectrum around $3 \times 10^{15}$ eV, so-called ``knee'', is composed of nuclei heavier than helium. This is the first measurement of the differential energy spectra of primary protons and heliums by selecting them event by event at the knee energy region.

Results on steady TeV $\gamma$-ray point source search using data taken from the Tibet HD (Feb. 1997 -- Sep. 1999) and Tibet III (Nov. 1999 -- Oct. 2001) arrays are presented. From $0^{\circ}$ to $60^{\circ}$ in declination, significant excesses from the well-known steady source Crab Nebula and the high state of the flare type source Markarian 421 are observed. Because the levels of significance from other positions are not sufficiently high, 90% confidence level upper limits on the flux are set assuming different power law spectra. To allow crosschecking, two independently developed analyses are used in this work.

We present the large-scale sidereal anisotropy ofgalactic cosmic-ray intensity in the multi-TeV region observed with the Tibet-IIIair shower array during the period from 1999 through 2003. The sidereal daily variation of cosmic rays observed in this experiment shows an excess of relative intensity around $4\sim7 $ hours local sidereal time, as well as a deficit around 12 hours local sidereal time. While the amplitude of the excess is not significant when averaged over all declinations, the excess in individual declinaton bands becomes larger and clearer as the viewing direction moves toward the south. The maximum phase of the excess intensity changes from $\sim$7 at the northern hemisphere to $\sim$4 hours at the equatorial region. We also show that both the amplitude and the phase of the first harmonic vector of the daily variation are remarkably independent of primary energy in the multi-TeV region. This is the first result determining the energy and declination dependences of the full 24-hour profiles of the sidereal daily variation in the multi-TeV region with a single air shower experiment.

We report on the solar diurnal variation of the galactic cosmic-ray intensity observed by the Tibet III air shower array during the period from 1999 to 2003. In the higher-energy event samples (12 TeV and 6.2 TeV), the variations are fairly consistent with the Compton-Getting anisotropy due to the terrestrial orbital motion around the sun, while the variation in the lower-energy event sample (4.0 TeV) is inconsistent with this anisotropy. This suggests an additional anisotropy superposed at the multi-TeV energies, e.g. the solar modulation effect. This is the highest-precision measurement of the Compton-Getting anisotropy ever made.

Several strong TeV gamma-ray flares were detected from Markarian 421 in the years 2000 and 2001 by the Tibet III air shower array at a level of statistical significance of 5.1 sigma. Markarian 421 was unprecedentedly active at X-ray and TeV gamma-ray energies during this period, and a positive correlation was found between the change of the ASM/RXTE X-ray flux and the Tibet TeV gamma-ray flux. When a power-law energy spectrum for gamma-rays from this source is assumed, the spectral index is calculated to be -3.24+-0.69 at the most active phase in 2001. The spectral index observed by the Tibet air shower array is consistent with those obtained via imaging air Cerenkov telescopes.

Since 1996, a hybrid experiment consisting of the emulsion chamber and burst detector array and the Tibet-II air-shower array has been operated at Yangbajing (4300 m above sea level, 606 g/cm^2) in Tibet. This experiment can detect air-shower cores, called as burst events, accompanied by air showers in excess of about 100 TeV. We observed about 4300 burst events accompanied by air showers during 690 days of operation and selected 820 proton-induced events with its primary energy above 200 TeV using a neural network method. Using this data set, we obtained the energy spectrum of primary protons in the energy range from 200 to 1000 TeV. The differential energy spectrum obtained in this energy region can be fitted by a power law with the index of -2.97 $\pm$ 0.06, which is steeper than that obtained by direct measurements at lower energies. We also obtained the energy spectrum of helium nuclei at particle energies around 1000 TeV.

We have shown that the Sun's shadow by high energy cosmic rays moves year by year and its behavior is correlated with a time variation of the large-scale structure of the solar and interplanetary magnetic fields. The solar activity was near minimum in the period from 1996 through 1997. Using the data obtained with the Tibet air shower array, we examined the shadowing of cosmic rays by the Sun in this quiet phase of solar cycle, and found that the Sun's shadow was just in the apparent direction of the Sun, though it was observed at the position considerably away from the Sun to the south-west in the period between 1990 and 1993. It is known that the magnetic pole of equivalent solar dipole was reversed during the previous active phase, and near solar minimum the dipole was aligned with the rotating axis, preserving its N-pole on the north pole side of the Sun. This causes the solar magnetic field to shift the Sun's shadow to the east. Thus, the observed results suggest that the shift of the Sun's shadow due to the solar magnetic field was pushed back by the effect of the geomagnetic field, since the geomagnetic field always make the shadow shift to the west. We discuss the Sun's shadow observed during the period near solar minimum in 1996-1997 and compare it with the simulation results.

Since 1996, a hybrid experiment consisting of an emulsion chamber and a burst detector array and the Tibet-II air shower array has been operated at Yangbajing (4300 m above sea level) in Tibet. This experiment can detect air-shower cores, called burst events, accompanied by air showers in excess of about 100 TeV. Using the burst event data observed by this experiment, we discuss the primary cosmic ray composition around the knee in comparison with the Monte Carlo simulations. In this paper, we show that all the features of burst events are wholly compatible with the heavy enriched composition in the knee energy region.

In 1997, the BL Lac Object Mrk 501 entered a very active phase and was the brightest source in the sky at TeV energies, showing strong and frequent flaring. Using the data obtained with a high density air shower array that has been operating successfully at Yangbajing in Tibet since 1996, we searched for gamma-ray signals from this source during the period from February through August in 1997. Our observation detected multi-TeV $\gamma$-ray signals at the 3.7-Sigma level during this period. The most rapid increase of the excess counts was observed between April 7 and June 16 and the statistical significance of the excess counts in this period was 4.7-Sigma. Among several observations of flaring TeV gamma-rays from Mrk 501 in 1997, this is the only observation using a conventional air shower array. We present the energy spectrum of gamma-rays which will be worthy to compare with those obtained by imaging atmospheric Cerenkov telescopes.

The Tibet experiment, operating at Yangbajing (4,300 m above sea level), is the lowest energy air shower array and the new high density array constructed in 1996 has sensitivity to $\gamma$-ray air showers at energies as low as 3 TeV. With this new array, the Crab Nebula was observed in multi-TeV $\gamma$-rays and a signal was detected at the 5.5 $\sigma$ level. We also obtained the energy spectrum of $\gamma$-rays in the energy region above 3 TeV which partially overlaps those observed with imaging atmospheric Cherenkov telescopes. This is the first observation of $\gamma$-ray signals from point sources with a conventional air shower array using scintillation detectors.