We present preliminary results obtained by a novel difference method for the study of the nature of the knee in the energy spectrum of the primary cosmic radiation. We have applied this method to data from the GAMMA experiment in Armenia. The analysis provides evidence for the possible existence of a nearby source of primary cosmic rays in the Southern hemisphere.

The present status of the GAMMA facility consisting of an enlarged surface EAS array (116 of 1 m^2 scintillation detectors) and underground muon carpet (150 m^2 detectors) is described. The recent results on mass composition and energy spectrum at the energy region above the knee obtained on the basis of the GAMMA experimental data are presented. It is shown that the power law after the knee is not invariable like -3.1. The slope of the energy spectrum becomes more flat at E0>20 PeV. The strong irregularities of the energy spectrum at about 70-80 PeV are discussed in comparison with other experiments. The bump can be described by a two-component model of primary cosmic ray origin, where additional (pulsar) Fe components are included with a very flat power law energy spectrum.

We present all-particle primary cosmic-ray energy spectrum in the 3-200 PeV energy range obtained by a multi-parametric event-by-event evaluation of the primary energy. The results are obtained on the basis of an expanded EAS data set detected at mountain level (700 g/cm^2) by the GAMMA experiment. The energy evaluation method has been developed using the EAS simulation with the SIBYLL interaction model taking into account the response of GAMMA detectors and reconstruction uncertainties of EAS parameters. Nearly unbiased (<5%) energy estimations regardless of a primary nuclear mass with an accuracy of about 15-10% in the 3-200 PeV energy range respectively are attained. An irregularity ('bump') in the spectrum is observed at primary energies of ~74 PeV. This bump exceeds a smooth power-law fit to the data by about 4 standard deviations. Not rejecting stochastic nature of the bump completely, we examined the systematic uncertainties of our methods and conclude that they cannot be responsible for the observed feature.

On the basis of the extensive air shower (EAS) data obtained by the GAMMA experiment, the energy spectra and elemental composition of the primary cosmic rays are derived in the 1-100 PeV energy range. The reconstruction of the primary energy spectra is carried out using an EAS inverse approach in the framework of the SIBYLL2.1 and QGSJET01 interaction models and the hypothesis of power-law primary energy spectra with rigidity-dependent knees. The energy spectra of primary H, He, O-like and Fe-like nuclei obtained with the SIBYLL interaction model agree with corresponding extrapolations of the balloon and satellite data to ~1 PeV energies. The energy spectra obtained from the QGSJET model show a predominantly proton composition in the knee region. The rigidity-dependent knee feature of the primary energy spectra for each interaction model is displayed at the following rigidities: ~2.5+/-0.2 PV (SIBYLL) and ~3.1-4.2 PV (QGSJET). All the results presented are derived taking into account the detector response, the reconstruction uncertainties of the EAS parameters, and fluctuations in the EAS development.

On the basis of the Extensive Air Shower data observed by the GAMMA experiment the energy spectra and elemental composition of the primary cosmic rays have been derived in the 1-100 PeV energy range. Reconstruction of the primary energy spectra is carried out in the framework of the SIBYLL and QGSJET interaction models and the hypothesis of the power-law steepening primary energy spectra. All presented results are derived taking into account the detector response, reconstruction uncertainties of EAS parameters and fluctuation of EAS development.

The CELESTE experiment will be an Atmospheric Cherenkov detector designed to bridge the gap in energy sensitivity between current satellite and ground-based gamma-ray telescopes, 20 to 300 GeV. We present test results made at the former solar power plant, Themis, in the French Pyrenees. The tests confirm the viability of using a central tower heliostat array for Cherenkov wavefront sampling.