Current Monte Carlo simulations do not provide a good description of the muon component extensive air showers. Many shower experiments report discrepancies between their data and predictions, ranging from TeV scale up to highest energies. In these proceedings, we address seasonal variation multi-muon events observed by NOvA Near Detector (ND). For our studies, use general-purpose code FLUKA tre...