Lipid coated nanobubbles (NBs) have attracted a great level of interest as ultrasound (US) contrast agents due to their ability to extravagate through leaky tumor vasculature. Their linear resonance frequency is in the range of ∼50 M Hz–200 MHz, leading to confusion over their observed strong contrast in diagnostic US frequencies. By solving the Marmottant model, the dynamics of uncoated and lipid coated NBs and microbubbles (MBs) are studied over the frequency and pressure ranges (6–12 MHz, 0.1–1.2 MPa) generally used in diagnostic US. A novel bifurcation analysis in tandem with the analysis of the frequency component of the scattered pressure are conducted. Results show that despite the increased linear resonance frequency and viscous damping due to the lipid shell, buckling and rupture of the shell enhances the generation of the 2nd and 3rd harmonic resonances at pressures as low as 0.2 MPa, not observed with uncoated NBs. The generation of the harmonic resonances are concomitant with an abrupt increase in the 2nd and 3rd harmonic frequency component of the scattered pressure with their pressure threshold (PT) increasing with decreasing NBs size. For the same gas volume, and above the PT, the maximum non-destructive 2nd and 3rd harmonic powers of NBs can become higher than the 2–4 μm MBs. Similar to the lower subharmonic pressure threshold of MBs, the dynamic variation of the NBs effective surface tension due to buckling and rupture may be the potential reason behind the observed harmonic echogenicity.