In natural convection flows, the fluid velocities are highly dependent on the thermal field and property variations can have a strong effect on both the flow and thermal fields. To examine the effect of property variations, at first, numerical analyses covering wide range of the Prandtl number under the same Rayleigh numbers have been carried out. Next, we have modeled the viscosity and thermal conductivity as parabolic functions of temperature and a comprehensive set of numerical solutions have been obtained to understand the effect. The Prandtl number dependence of Nusselt number is fairly strong even though the effect is still weak compared to the Rayleigh number dependence. When thermophysical properties are dependent on temperature, the flow field showed a fairly weak variation except near boundaries, whereas the temperature field is strongly affected, especially by the temperature dependent thermal conductivity.