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 investigated 64 pairs of interacting-CME events identified by the simultaneous observations of SOHO and STEREO spacecraft from 2010 January to 2014 August, to examine the relationship between the large SEP events in the energy of 25-60MeV and the properties of the interacting CMEs. We found that during CME interactions the large SEP events in this study were all generated by CMEs with the presence of enhanced type II radio bursts, which also have wider longitudinal distributions comparing to events with the absence of type II radio burst or its enhancement (almost associated with small SEP events). It seems that the signature of type II radio bursts enhancement is a good discriminator between large SEP and small or none SEP event producers during CME interactions. The type II radio burst enhancement is more likely to be generated by CME interactions, with the main CME having larger speed (v), angular width (WD), mass (m) and kinetic energy (E_k), that taking over the preceding CMEs which also have higher v, WD, m and E_k, than those preceding CMEs in CME pairs missing the type II radio bursts or enhancements. Our analysis also revealed that the intensities of associated SEP events correlate negatively with the intersection height of the two CMEs. Most of type-II-enhanced events and SEP events are coincidentally and almost always made by the fast and wide main CMEs that sweeping fully over a relatively slower and narrower preceding CMEs. We suggest that a fast CME with enough energy completely overtaking a relatively narrower preceding CME, especially in low height, can drive more energetic shock signified by the enhanced type II radio bursts. The shock may accelerate ambient particles and lead to large SEP event more easily.

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.

It has been reported that a 5.7sigma directional muon excess coincident with the 2000 July 14 solar flare was registered by the L3 precision muon spectrometer [Ruiguang Wang, Astroparticle Phys., 31(2009) 149]. Using a same analysis method and similar criteria of event selection, we have analyzed the L3 precision muon spectrometer data during November 2000. The result shows that a 4.7sigma muon excess appeared at a time coincident with the solar flare of 8 of November 2000. This muon excess corresponds to above 40 GeV primary protons which came from a sky cell of solid angle 0.048 sr. The probability of being a background fluctuation is estimated to be about 0.1%. It has been convinced that solar protons could be accelerated to tens of GeV in those Class X solar flares which usually arose solar cosmic ray ground level enhancement (GLE) events. However, whether a Class M solar flare like the non-GLE event of 8 November 2000 may also accelerate solar protons to such high energies? It is interesting and noteworthy.

It has been recently suggested that large solar energetic particle (SEP) events are often caused by twin CMEs. In the twin-CME scenario, the preceding CME is to provide both an enhanced turbulence level and enhanced seed population at the main CME-driven shock. In this work, we study the effect of the preceding CMEs on the seed population. We examine event-integrated abundance of iron to oxygen ratio (Fe/O) at energies above 25 MeV/nuc for large SEP events in solar cycle 23. We find that the Fe/O ratio (normalized to the reference coronal value of $0.134$) $\leq2.0$ for almost all single-CME events and these events tend to have smaller peak intensities. In comparison, the Fe/O ratio of twin-CME events scatters in a larger range, reaching as high as $8$, suggesting the presence of flare material from perhaps preceding flares. For extremely large SEP events with peak intensity above $1000$ pfu, the Fe/O drop below $2$, indicating that in these extreme events the seed particles are dominated by coronal material than flare material. The Fe/O ratios of Ground level enhancement (GLE) events, all being twin-CME events, scatter in a broad range. For a given Fe/O ratio, GLE events tend to have larger peak intensities than non-GLE events. Using velocity dispersion analysis (VDA), we find that GLE events have lower solar particle release (SPR) heights than non-GLE events, \red{agreeing with earlier results by Reames 2009b.

Recently it has been suggested that the "twin-CME" scenario Li.etal2012 may be a very effective mechanism in causing extreme Solar Energetic Particle (SEP) events and in particular Ground Level Enhancement (GLE) events. Ding.etal2013 performed a statistical examination of the twin-CME scenario with a total of $126$ fast and wide western Coronal Mass Ejections (CMEs). They found that CMEs having a preceding CME with a speed $>$ 300 $km/s$ within $9$ hours from the same active region have larger probability of leading to large SEP events than CMEs that do not have preceding CMEs. The choice of $9$ hours being the time lag $\tau$ between the preceding CME and the main CME was based on some crude estimates of the decay time of the turbulence downstream of the shock driven by the preceding CME. In this work, we examine this choice. For the $126$ fast wide CMEs examined in Ding.etal2013, we vary the time lag $\tau$ from $1$ hour to $24$ hours with an increment of $1$ hour. By considering three quantities whose values depend on the choice of this time lag $\tau$, we show that the choice of $13$ hours for $\tau$ is more appropriate. Our study confirms our earlier result that twin CMEs are more likely to lead to large SEP events than single fast CMEs. The results shown here are of great relevance to space weather studies.

On 2012 March 7, the STEREO Ahead and Behind spacecraft, along with the near-earth spacecraft (e.g. SOHO, Wind) situated between the two STEREO spacecraft, observed an extremely large global solar energetic particle (SEP) event in Solar Cycle 24. Two successive coronal mass ejections (CMEs) have been detected close in time. From the multi-point in-situ observations, it can be found that this SEP event was caused by the first CME, and the second one was not involved. Using the velocity dispersion analysis (VDA), we find that for well magnetically connected point, the energetic protons and electrons are released nearly at the same time. The path lengths to STEREO-B(STB) of protons and electrons have distinct difference and deviate remarkably from the nominal Parker spiral path length, which is likely due to the presence of interplanetary magnetic structures situated between the source and the STB. Also the VDA method seems only to obtain reasonable results at well-connected locations and the inferred energetic particles release times in different energy channels are similar. We suggest that good-connection is crucial for obtaining both accurate release time and path length simultaneously, agreeing with the modeling result of Wang & Qin (2015).

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 EAS hybrid experiment has been designed by constructing a YAC (Yangbajing Air shower Core) detector array inside the existing Tibet-III air shower array. The first step of YAC, called "YAC-I", consists of 16 plastic scintillator units (4 rows times 4 columns) each with an area of 40 cm * 50 cm which is used to check hadronic interaction models used in AS simulations. A Monte Carlo study shows that YAC-I can record high energy electromagnetic component in the core region of air showers induced by primary particles of several tens TeV energies where the primary composition is directly measured by space experiments. It may provide a direct check of the hadronic interaction models currently used in the air shower simulations in the corresponding energy region. In present paper, the method of the observation and the sensitivity of the characteristics of the observed events to the different interaction models are discussed.

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.

Using the data taken from Tibet II High Density (HD) Array (1997 February-1999 September) and Tibet-III array (1999 November-2005 November), our previous northern sky survey for TeV $\gamma-$ray point sources has now been updated by a factor of 2.8 improved statistics. From $0.0^{\circ}$ to $60.0^{\circ}$ in declination (Dec) range, no new TeV $\gamma-$ray point sources with sufficiently high significance were identified while the well-known Crab Nebula and Mrk421 remain to be the brightest TeV $\gamma-$ray sources within the field of view of the Tibet air shower array. Based on the currently available data and at the 90% confidence level (C.L.), the flux upper limits for different power law index assumption are re-derived, which are approximately improved by 1.7 times as compared with our previous reported limits.

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 new type of cosmic-ray experiment is under construction in the Pyh\"asalmi mine in the underground laboratory of the University of Oulu, Finland. It aims to study the composition of cosmic rays at and above the knee region. The experiment, called EMMA, will cover approximately 150 square-metres of detector area. The array is capable of measuring the multiplicity and the lateral distribution of underground muons, and the arrival direction of the air shower. The full-size detector is expected to run by the end of 2007.

The DAMA experiments are running deep underground in the Gran Sasso National Laboratory. Several interesting results have been achieved so far. Here a maximum likelihood method to search for the WIMP annual modulation signature is discussed and applied to a set of preliminary test data collected with large mass highly radiopure NaI(Tl) detectors. Various related technical arguments are briefly addressed.