The WILLI detector, built in IFIN-HH Bucharest, in collaboration with KIT Karlsruhe, is a rotatable modular detector for measuring charge ratio for cosmic muons with energy $<$ 1 GeV. It is under construction a mini-array for measuring the muon charge ratio in Extensive Air Showers. The EAS simulations have been performed with CORSIKA code. The values of the muon flux, calculated with semi-analytical formula, and simulated with CORSIKA code, based on DPMJET and QGSJET models for the hadronic interactions, are compared with the experimental data determined with WILLI detector. No significant differences between the two models and experimental data are observed. The measurements of the muon charge ratio for different angles-of-incidence, (performed with WILLI detector) shows an asymmetry due to the influence of magnetic field on muons trajectory; the values are in agreement with the simulations based on DPMJET hadronic interaction model. The simulations of muon charge ratio in EAS performed with CORSIKA code based on three hadronic interaction models (QGSJET2, EPOS and SYBILL) show relative small difference between models for H and for the Fe showers; the effect is more pronounced at higher inclination of WILLI detector. The future measurements should indicate which model is suitable.

In view of the current efforts to extend the KASCADE experiment (KASCADE-Grande) for observations of Extensive Air Showers (EAS) of primary energies up to 1 EeV, the features of muon arrival time distributions and their correlations with other observable EAS quantities have been scrutinised on basis of high-energy EAS, simulated with the Monte Carlo code CORSIKA and using in general the QGSJET model as generator. Methodically various correlations of adequately defined arrival time parameters with other EAS parameters have been investigated by invoking non-parametric methods for the analysis of multivariate distributions, studying the classification and misclassification probabilities of various observable sets. It turns out that adding the arrival time information and the multiplicity of muons spanning the observed time distributions has distinct effects improving the mass discrimination. A further outcome of the studies is the feature that for the considered ranges of primary energies and of distances from the shower axis the discrimination power of global arrival time distributions referring to the arrival time of the shower core is only marginally enhanced as compared to local distributions referring to the arrival of the locally first muon.