In this paper, we present here a detailed study of several observables as a function of pseudorapidity from several hadron production models commonly used for extended air shower simulations and compare these predictions with [Formula: see text] collisions measured by the CMS experiment at [Formula: see text] TeV. To facilitate the model/data comparisons, the pseudorapidity binning used in the energy density spectra data was also applied to the Monte Carlo analysis ([Formula: see text] and [Formula: see text]). Given that no model currently describes the full breadth of available hadronic collision data, four independent models were used in the comparison study, each with differences in the handling of soft quantum chromodynamics processes, in particular. The analysis results show that predictions from EPOS-LHC provide good agreement with the data, in contrast with the other models used in the study (QGSJETII-04, DPMJET-3.19, Sibyll2.3d). Significant modifications may be required in the rest of the model-based generators for them to provide a better description of the experimental observations.