The spatial and temporal development of shear-induced, overturning billows associated with breaking internal solitary waves is studied by means of a combined laboratory and numerical investigation. The waves are generated in the laboratory by a lock exchange mechanism and they are simulated numerically via a contour-advective, semi-Lagrangian method. The properties of individual billows (maximum height attained, time of collapse, growth rate, speed, wave length, Thorpe scale) are determined in each case and the billow interaction processes are studied and classified. For broad, flat waves, similar characteristics are seen to those in parallel shear flow but, for waves not at the conjugate flow limit, billow characteristics are affected by the spatially-varying, wave-induced shear flow. Wave steepness and wave amplitude are shown to have a crucial influence on determining the type of interaction that occurs between billows and whether billow overturning can be arrested. Examples are given in which billows (i) evolve independently of one another (ii) pair with one another, (iii) engulf/entrain one another and (iv) fail to completely overturn. It is shown that the vertical extent a billow can attain (and the associated Thorpe scales of the billows) are dependent on wave amplitude but that their values saturate once a given amplitude is reached. It is interesting to note that this amplitude is less than the conjugate flow limit amplitude. The number of billows that form on a wave are shown to be dependent on wavelength; shorter waves support fewer but bigger billows than their long wave counterparts for a given stratification.Key words: Authors should not enter keywords on the manuscript, as these must be chosen by the author during the online submission process and will then be added during the typesetting process (see http://journals.cambridge.org/data/relatedlink/jfmkeywords.pdf for the full list)