Voyager 2 observations of electrostatic electron and ion cyclotron harmonic waves in Neptune's magnetosphere are addressed. A model of electron Bernstein modes generated by a loss cone distribution of superthermal electrons is scaled to Neptune parameters and a comparison of theory with the observed electron flux shows good agreement. A model of proton Bernstein modes generated by a ring distribution of Tritonogenic nitrogen ions is also investigated and satisfactory agreement with the data are obtained compatible with known properties of the magnetosphere. The success of the model in accounting for electrostatic emissions observed by Voyager over a wide range of sampled parameters recommends its general applicability to planetary magnetospheres.