The nature of instability and disturbance amplification in the laminar natural convection boundary layer over a vertical flat surface with uniform heat flux has been theoretically investigated. The coupled Orr-Sommerfeld equation has been numerically integrated for a Prandtl number of 6ยท7, with the boundary condition that the disturbance heat flux be zero at the surface. The spatial amplification characteristics of disturbances convected downstream were analyzed, and constant amplification rate contours were determined. The relative amplification has been calculated from these contours and is presented in the form of amplitude ratio contours. An important feature of these results is that the low frequency disturbances, which become unstable first, amplify very slowly and also have wavelengths which are much longer than the distance to the leading edge. The higher frequency, shorter wavelength, disturbances amplify much faster and are, therefore, presumed to be the dominant ones in stability considerations. The nature of the velocity and temperature amplitudes and phase profiles across the boundary layer has also been examined.