The origin of Galactic CRs up the knee energy remains unanswered and provides strong motivation for the study of gamma-ray sources at energies above 10 TeV. We discuss recent results from ground-based gamma-ray Cherenkov imaging systems at these energies as well as future observational efforts in this direction. The exciting results of H.E.S.S. give clues as to the nature of Galactic CR accelerators, and suggest that there is a population of Galactic gamma-ray sources with emission extending beyond 10 TeV. A dedicated system of Cherenkov imaging telescopes optimised for higher energies appears to be a promising way to study the multi-TeV gamma-ray sky.

We investigate properties of astrophysical electromagnetic cascades developed in matter, photon gas and magnetic fields, and discuss similarities and differences between characteristics of electron-photon showers developed in these 3 substances.

We suggest a new approach to study the cosmis ray (CR) mass composition in the energy region from 30 TeV/nucleus up to the "knee" region, i.e. up to a few PeV/nucleus, using an array of imaging atmospheric Cherenkov telescopes (IACTs) of a special architecture. This array consists of telescopes with a relatively small mirror size (~10 square meters) separated from each other by large distances (~500 meters) and equipped by multichannel cameras with a modest pixel size (0.3-0.5 degree) and a sufficiently large viewing angle (6-7 degree). Compared to traditional IACT systems (like HEGRA, HESS or VERITAS) the IACT array considered here could provide a very large detection area (several square kilometers or more). At the same time, it allows an accurate measurement of the energy of CR induced air showers (the energy resolution ranges within 25-35%) and an effective separation of air showers created by different nuclei. Particularly, it is possible to enrich air showers belonging to the nucleus group assigned for selection up to ~90% purity at a detection efficiency of 15-20% of such showers.

The spectra of very high energy gamma-radiation from distant extragalactic objects suffer significant deformations during the passage of primary gamma-rays through the intergalactic medium. The recently reported fluxes of diffuse infrared background radiation indicate that we detect, most probably, heavily absorbed TeV radiation from Mkn 421 and Mkn 501. This implies that the absorption-corrected spectrum of Mkn 501 may contain a sharp pile-up which contradicts to the predictions of the conventional models of TeV blazars, and thus may leads to the so-called "IR background-TeV gamma-ray crisis". To overcome this difficulty, in this paper we propose two independent hypotheses assuming that (i) the TeV radiation from Mkn 501 has a secondary origin, i.e. it is formed during the development of electron-photon cascades in the intergalactic medium initiated by primary gamma-rays; (ii) the pile-up in the source spectrum is a result of comptonization (in deep Klein-Nishina regime) of ambient optical radiation by an ultrarelativistic conical cold outflow (jet) with bulk motion Lorentz factor Gamma_0 >= 3 10^7. Within the uncertainties caused by the limited energy resolution of spectral measurements, the observed TeV radiation of Mkn 501 formally can be explained by the intergalactic cascade gamma-rays, assuming however an extremely low intergalactic magnetic field in the direction to the source at the level of < 10^{-18} G. We also demonstrate that the "bulk motion comptonization" scenario can quite naturally reproduce the unusual spectral features in the absorption-corrected TeV spectrum of Mkn 501, and briefly discuss the astrophysical implications of this hypothesis.

The technique of adjoint cascade equations has been applied to calculate the properties of extremely high energy gamma-rays in the energy range 10^18--10^22 eV with taking into account the LPM effect and interactions of gamma-rays with the geomagnetic field. Such characteristics are analysed as the electron and muon contents at the observation level, the electron cascade curves, the lateral distribution functions of photoproduced muons.

The sensitivity of ground-based imaging atmospheric Cherenkov gamma-ray observatories depends critically on the primary particle identification methods which are used to retain photon-initiated events and suppress the spurious background produced by cosmic rays. We suggest a new discrimination technique which utilizes differences in the fluctuations of the light intensity in the images of showers initiated by photons and those initiated by protons or heavier nuclei. The database of simulated events for the proposed VERITAS observatory has been used to evaluate the efficiency of the new technique. Analysis has been performed for both a single VERITAS imaging telescope, and a system of these telescopes. We demonstrate that a discrimination efficiency of > 1.5 - 2.0 can be achieved in addition to traditional background rejection methods based on image shape parameters.

Ground-based atmospheric Cherenkov telescopes are proven to be effective instruments for observations of very high energy (VHE) gamma-radiation from celestial objects. For effective use of such technique one needs detailed Monte Carlo simulations of gamma-ray- and proton/nuclei-induced air showers in Earth atmosphere. Here we discuss in detail the algorithms used in the numerical code ALTAI developed particularly for the simulations of Cherenkov light emission from air showers of energy below 50 TeV. The specific scheme of sampling the charged particle transport in the atmosphere allows the performance of very fast and accurate simulations used for interpretation of the VHE gamma-ray observations.