We revisit neutrino oscillations in matter considering the open quantum system framework, which allows to introduce possible decoherence effects generated by New Physics in a phenomenological manner. We assume that the decoherence parameters γ i j may depend on the neutrino energy, asThe case of non-uniform matter is studied in detail and, in particular, we develop a consistent formalism to study the nonadiabatic case dividing the matter profile into an arbitrary number of layers of constant densities. This formalism is then applied to explore the sensitivity of IceCube and DeepCore to this type of effects. Our study is the first atmospheric neutrino analysis where a consistent treatment of the matter effects in the three-neutrino case is performed in presence of decoherence. We show that matter effects are indeed extremely relevant in this context. We find that IceCube is able to considerably improve over current bounds in the solar sector (γ 21 ) and in the atmospheric sector (γ 31 and γ 32 ) for n = 0, 1, 2 and, in particular, by several orders of magnitude (between 3 and 9) for the n = 1, 2 cases. For n = 0 we find γ 32 , γ 31 < 4.0 × 10 −24 (1.3 × 10 −24 ) GeV and γ 21 < 1.3 × 10 −24 (4.1 × 10 −24 ) GeV at the 95% CL, for normal (inverted) mass ordering.