The methodical properties of the original difference method for the search of the anisotropy at the knee region of the primary cosmic radiation energy spectrum are analyzed. The main feature of the suggested method is a study of the difference in the EAS characteristics in different directions but not their intensity. It is shown that the method is stable to the random experimental errors and allows to separate the anomalies related to the laboratory coordinate system from the anomalies in the celestial coordinates. The method uses multiple scattering of the charge particles in the Galaxy magnetic fields to study the whole celestial sphere including the regions outside of the line of sight of the installation.

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 energy spectrum of the primary cosmic radiation in the energy range 1 - 100 PeV and the extensive air shower (EAS) characteristics obtained on the basis of the expanded data bank of the GAMMA experiment (Mt. Aragats, Armenia) are presented. With increased statistics we confirm our previous results on the energy spectrum. The spectral index above the knee is about -3.1, but at energies beyond 20 PeV a flattening of the spectrum is observed. The existence of the 'bump' at about 70 PeV is confirmed with a significance of more than 4{\sigma}. In the energy range of 10 - 100 PeV the shower age becomes energy independent and we observe a direct proportionality of the EAS size to the primary energy. This suggests an approximately constant depth of the EAS maximum in this energy range. This is evidence in favour of an increasing average mass of primary particles at energies above 20 PeV. The additional source scenario, which is a possible explanation of the 'bump' in the spectrum, also leads to the conclusion of increasing mass of the primary cosmic rays. A comparison with the data of other experiments is presented.

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 an upper limit of galactic diffuse gamma-ray flux, measured with the GAMMA experiment at energy about 175 TeV. The results were obtained using selection of muon poor extensive air showers at mountain level (700 g/cm^2, Mt. Aragats, Armenia) for 5 GeV muon energy threshold. An upper limit for the differential flux at energy E=175+/-25 TeV is equal to (5.8-7.0)x10^-12 [erg m^2 s sr]^-1 for 95% confidence level.

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.