This paper presents an efficient channel estimation method for near-field massive multi-input multi-output (MIMO) systems. To estimate the channel parameters in terms of angles and ranges to the scatterers at the transmitter and receiver sides, we formulate the signal model for massive MIMO systems as a function of angle and range parameters. The signal model is then used to construct the fourth-order cumulant matrix, from which 2-dimensional (2D) multiple signal classification (MUSIC) pseudo-spectrums for the angles and ranges are derived. The effectiveness of the proposed algorithm is verified through simulation by comparing with the conventional channel estimation techniques.