The detection of gravitational waves and neutrinos from astrophysical sources with gamma-ray counterparts officially started the era of Multi-Messenger Astronomy. Their transient and extreme nature implies that monitoring the VHE sky is fundamental to investigate the non-electromagnetic signals. However, the limited effective area of space-borne instruments prevents observations above a few hundred GeV, while the small field of view and low duty cycle of IACTs make them unsuited for extensive monitoring activities and prompt response to transients. Extensive Air Shower arrays (EAS) can provide a large field of view, a wide effective area and a very high duty cycle. Their main difficulty is the distinction between gamma-ray and cosmic-ray initiated air showers, especially below the TeV range. Here we present some case studies stressing the importance that a new EAS array in the Southern Hemisphere will be able to survey the sky from below 100 GeV up to several TeV. In the energy domain between 100 and 400 GeV we expect the strongest electromagnetic signatures of the acceleration of ultra-relativistic particles in sources like SNRs, blazar jets and gamma-ray bursts, as recently proved by IACT observations. This spectral window is also crucial to understand the Universe opacity to high energy radiation, thus providing constraints on the cosmological parameters. We will discuss the implications of VHE radiation on the mechanisms at work and we will focus on the advantages resulting from the ability to monitor the energy window lying between the domain of space-borne detectors and ground-based facilities.

A new concept for the direct measurement of muons in air showers is presented. The concept is based on resistive plate chambers (RPCs), which can directly measure muons with very good space and time resolution. The muon detector is shielded by placing it under another detector able to absorb and measure the electromagnetic component of the showers such as a water-Cherenkov detector, commonly used in air shower arrays. The combination of the two detectors in a single, compact detector unit provides a unique measurement that opens rich possibilities in the study of air showers.

We focus on the primary composition of cosmic rays with the highest energies that cause extensive air showers in the Earth's atmosphere. A way of examining the two lowest order moments of the sample distribution of the depth of shower maximum is presented. The aim is to show that useful information about the composition of the primary beam can be inferred with limited knowledge we have about processes underlying these observations. In order to describe how the moments of the depth of shower maximum depend on the type of primary particles and their energies, we utilize a superposition model. Using the principle of maximum entropy, we are able to determine what trends in the primary composition are consistent with the input data, while relying on a limited amount of information from shower physics. Some capabilities and limitations of the proposed method are discussed. In order to achieve a realistic description of the primary mass composition, we pay special attention to the choice of the parameters of the superposition model. We present two examples that demonstrate what consequences can be drawn for energy dependent changes in the primary composition.

We consider a hypothetical observatory of ultra-high energy cosmic rays consisting of two surface detector arrays that measure independently electromagnetic and muon signals induced by air showers. Using the constant intensity cut method, sets of events ordered according to each of both signal sizes are compared giving the number of matched events. Based on its dependence on the zenith angle, a parameter sensitive to the dispersion of the distribution of the logarithmic mass of cosmic rays is introduced. The results obtained using two post-LHC models of hadronic interactions are very similar and indicate a weak dependence on details of these interactions.