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.