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.

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.

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.