The interaction of aerosol-cloud-precipitation has an important impact on the global climate. The understanding of this issue is related to the uncertainty of climate change prediction. The traditional indirect effect of aerosols suggests that when the number of aerosols increases, it will act to suppress precipitation. However, recent studies on satellite observations have found that aerosols are positively correlated with precipitation, which is contrary to conventional views. This study attempts to use the A-Train satellite product to verify the correlation between aerosol and precipitation, and further reveals the possible physical mechanism of the positive relationship between aerosols and precipitation in satellite observations through the three-dimensional structure of clouds. The study found that the precipitation intensity is positively correlated with the aerosol optical depth, while the relationship between the cloud droplet concentration and the precipitation intensity is related to the liquid water path; and when the number of aerosols increases, the radar reflectance spectrum is widened, and the precipitation increases, while the cloud droplet concentration shows the opposite phenomenon. It can be concluded that the possible cause of positive correlation between aerosol and precipitation is the negative correlation between CDNC and AOD.

The $^{1}S_{0}$ nucleonic superfluids are investigated within the relativistic mean-field model and Bardeen-Cooper-Schrieffer theory in hyperonic neutron stars. The $^{1}S_{0}$ pairing gaps of neutrons and protons are calculated based on the Reid soft-core interaction as the nucleon-nucleon interaction. We have studied particularly the influence of hyperons degrees of freedom on the $^{1}S_{0}$ nucleonic pairing gap in neutron star matter. It is found that the appearance of hyperons has little impact on baryonic density range and size for the $^{1}S_{0}$ neutronic pairing gap, the $^{1}S_{0}$ protonic pairing gap also decreases slightly in this region $\rho_B=0.0-0.393$ fm$^{-3}$. However, if baryonic density becomes greater than 0.393 fm${^{-3}}$, the $^{1}S_{0}$ protonic pairing gap obviously increases. In addition, the protonic superfluid range is obviously enlarged due to the presence of hyperons. In our results, the hyperons change the $^{1}S_{0}$ protonic pairing gap which must change the cooling properties of neutron stars.

The paper investigates the overall and detailed features of cosmic ray (CR) spectra in the knee region using the scenario of nuclei-photon interactions around the acceleration sources. Young supernova remnants can be the physical realities of such kind of CR acceleration sites. The results show that the model can well explain the following problems simultaneously with one set of source parameters: the knee of CR spectra and the sharpness of the knee, the detailed irregular structures of CR spectra, the so-called "component B" of Galactic CRs, and the electron/positron excesses reported by recent observations. The coherent explanation serves as evidence that at least a portion of CRs might be accelerated at the sources similar to young supernova remnants, and one set of source parameters indicates that this portion mainly comes from standard sources or from a single source.

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.

Supernova remnants have long been regarded as sources of the Galactic cosmic rays up to petaelectronvolts, but convincing evidence is still lacking. In this work we explore the common origin of the subtle features of the cosmic ray spectra, such as the knee of cosmic ray spectra and the excesses of electron/positron fluxes recently observed by ATIC, H.E.S.S., Fermi-LAT and PAMELA. Numerical calculation shows that those features of cosmic ray spectra can be well reproduced in a scenario with e$^+$e$^-$ pair production by interactions between high energy cosmic rays and background photons in an environment similar to the young supernova remnant. The success of such a coherent explanation serves in turn as an evidence that at least a portion of cosmic rays might be accelerated at young supernova remnants.

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.