A geomagnetic storm is the result of sustained interaction between solar wind with a southward magnetic field and the magnetosphere. To investigate the influence of various solar wind parameters on the intensity of major geomagnetic storm, 67 major geomagnetic storms that occurred between 1998 and 2006 were used to calculate the correlation coefficients (CCs) between the intensities of major geomagnetic storms and the time integrals of southward interplanetary magnetic field $B_s$, solar wind electric field ($E_y$) and injection function (Q) during the main phase of the associated geomagnetic storms. SYM-H$_{min}$ was used to indicate the intensity of the associated major geomagnetic storm, while I($B_z$), I($E_y$) and I(Q) were used to indicate the time integrals of $B_z$, $E_y$ and Q during the main phase of associated major geomagnetic storm respectively. The derived CC between I($B_z$) and SYM-H$_{min}$ is 0.33, while the CC between I($E_y$) and SYM-H$_{min}$ is 0.57 and the CC between I(Q) and SYM-H$_{min}$ is 0.86. The results provide statistical evidence that solar wind dynamic pressure or solar wind density plays a significant role in transferring solar wind energy into the magnetosphere, in addition to the southward magnetic field and solar wind speed. Solar wind that has a strong geoeffectiveness requires solar wind dynamic pressure $>$3 nPa or solar wind density $>3$ nPa$/V_{sw}^2$. Large and long duration $B_s$ alone cannot ensure a major geomagnetic storm, especially if the solar wind dynamic pressure is very low, as large and long duration Bs is not a full condition, only a necessary condition to trigger a major geomagnetic storm.

We have studied the dynamic proton spectra for the two solar energetic particle (SEP) events on 2000 July 14 (hereafter GLE59) and 2005 January 20 (hereafter GLE69). The source locations of GLE59 and GLE69 are N22W07 and N12W58 respectively. Proton fluxes >30 MeV have been used to compute the dynamic spectral indices of the two SEP events. The results show that spectral indices of the two SEP events increased more swiftly at early times, suggesting that the proton fluxes >30 MeV might be accelerated particularly by the concurrent flares at early times for the two SEP events. For the GLE69 with source location at N12W58, both flare site and shock nose are well connected with the Earth at the earliest time. However, only the particles accelerated by the shock driven by eastern flank of the CME can propagate along the interplanetary magnetic field line to the Earth after the flare. For the GLE59 with source location at N22W07, only the particles accelerated by the shock driven by western flank of the associated CME can reach the Earth after the flare. Results show that there was slightly more than one hour during which the proton spectra for GLE69 are softer than that for GLE59 after the flares, suggesting that the shock driven by eastern flank of the CME associated with GLE69 is weaker than the shock driven by the western flank of the CME associated with GLE59. The results support that quasi-perpendicular shock has stronger potential in accelerating particles than the quasi-parallel shock. The results also suggest that only a small part of the shock driven by western flank of the CME associated with the GLE59 is quasi-perpendicular.

An X1.7 flare at 10:15 UT and a halo CME with a projected speed of 942 km/s erupted from NOAA solar active region 9393 located at N20W19, were observed on 2001 March 29. When the CME reached the Earth, it triggered a super geomagnetic storm (hereafter super storm). We find that the CME always moved towards the Earth according to the intensity-time profiles of protons with different energies. The solar wind parameters responsible for the main phase of the super storm occurred on March 31, 2001 is analyzed taking into account the delayed geomagnetic effect of solar wind at the L1 point and using the SYM-H index. According to the variation properties of SYM-H index during the main phase of the super storm, the main phase of the super storm is divided into two parts. A comparative study of solar wind parameters responsible for the two parts shows the evidence that the solar wind density plays a significant role in transferring solar wind energy into the magnetosphere, besides the southward magnetic field and solar wind speed.

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.

According to the solar protons' data observed by Geostationary Operational Environmental Satellites (GOES) and neutron monitors on the ground of the Earth, and near-relativistic electrons data measured by the ACE spacecraft, the onset times of protons with different energies and near-relativistic electrons have been estimated and compared with the time of solar soft and hard X-ray and radio burst data. The results show that first arriving relativistic and non-relativistic protons and electrons may be accelerated by the concurrent flare. The results also suggest that release times of protons with different energies may be different, and the protons with lower energy may release earlier than those with higher energy. Some protons accelerated by concurrent flares may be further accelerated by CME-driven shock.

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 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.

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 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.