We model the unsteady evolution of turbulent buoyant plumes following temporal changes to the source conditions. The integral model is derived from radial integration of the governing equations expressing the evolution of mass, axial momentum and buoyancy in the plume. The non-uniform radial profiles of the axial velocity and density deficit in the plume are explicitly captured by shape factors in the integral equations; the commonly assumed top-hat profiles lead to shape factors equal to unity. The resultant model for unsteady plumes is hyperbolic when the momentum shape factor, determined from the radial profile of the mean axial velocity in the plume, differs from unity. The solutions of the model when source conditions are maintained at constant values are shown to retain the form of the well-established steady plume solutions. We demonstrate through a linear stability analysis of these steady solutions that the inclusion of a momentum shape factor in the governing equations that differs from unity leads to a well-posed integral model. Therefore, our model does not exhibit the mathematical pathologies that appear in previously proposed unsteady integral models of turbulent plumes. A stability threshold for the value of the shape factor is also identified, resulting in a range of its values where the amplitudes of small perturbations to the steady solutions decay with distance from the source. The hyperbolic character of the system of equations allows the formation of discontinuities in the fields describing the plume properties during the unsteady evolution, and we compute numerical solutions to illustrate the transient development of a plume following an abrupt change in the source conditions. The adjustment of the plume to the new source conditions occurs through the propagation of a pulse of fluid through the plume. The dynamics of this pulse is described by a similarity solution and, through the construction of this new similarity solution, we identify three regimes in which the evolution of the transient pulse following adjustment of the source qualitatively differs.