In this study, combinations of the global arrival time, ( τ global ), pseudorapidity, and lateral density distribution (ρ μ ) of muons, which are three mass-sensitive observables of cosmic ray induced extensive air showers, have been used as new parameters to study the primary mass discrimination around the knee energies (100 TeV-10 PeV). This is a simulation-based study and the simulations have been performed for the KASCADE array at Karlsruhe and the Alborz-I array at Tehran to study the effect of the altitude on the quality of the primary mass discrimination. The merit factors of the single and combined three mass-sensitive observables have been calculated to compare the discrimination power of combined and single observables. We have used the CORSIKA 7.4 code to simulate the extensive air showers (EASs) sample sets. Considering all aspects of our study, it is found that the ratio of the global time to the lateral density distribution of the muons gives better results than other ratios; also in the case of single observables, the muon density gives better results compared with the other observables. Also it is shown that below 1 PeV primary energies, the ratio of the muon global time to the muon density ( τ global /ρ μ ) results in a better mass discrimination relative to the muon density only.