A single power law flux spectrum of high energy neutrinos does not adequately explain the entire 60 TeV to 10 PeV event spectrum observed at IceCube, specially the excess of PeV events and the lack of Glashow resonance events expected at 6.3 PeV cannot be simultaneously explained by a single power law source neutrino flux. Here we consider a model of neutrino splitting ν → 3ν over cosmological distances. Starting from a single power-law spectrum expected from the astrophysical sources, we show that by adjusting the decay length and spectral index one can give a better fit to the observed IceCube events over the entire 1 TeV-6 PeV, compared to that from a single power spectrum. For ν → 3ν neutrino splitting, the flavor ratios of the daughter neutrinos are different from the standard oscillation or invisible decay cases and can be used as a test of this scenario. We propose a ν2HDM where a light Higgs (∼0.1 eV) mediates neutrino splitting via a one-loop box diagram. The split in the masses of the scalars in the doublet gives a large contribution to the oblique T parameter which is severely constrained. This constraint from the S,T,U oblique parameters can be evaded by the introduction of an extra vector lepton doublet and with mass ∼ 200 GeV.