Thermodynamic modelling of (1) osumilite solid solutions and (2) dehydration melting in pelitic compositions within the KFMASH system is quite successful in reproducing the invariant and univariant reactions determined in experimental studies. Even though rather preliminary, such melt thermodynamic models may be very useful in interpolating and extrapolating the limited information available from a small number of experimental runs. These methods allow the compositions of all phases to be monitored as a function of pressure, temperature and equilibrium phase assemblage for any desired bulk composition. Locating the higher variance phase fields (e.g. quadrivariant, quinivariant) is often difficult or impossible by inspection, but is made relatively easy using thermodynamic software such as thermocalc. In the KFMASH system the calculated partition of Fe and Mg between osumilite, garnet, cordierite, orthopyroxene and biotite are shown to be in good agreement with experimental and natural data and allow reliable calculation of mineral compositions coexisting with quartz-saturated and H2O-undersaturated melts for a variety of bulk compositions. These phase diagram calculations allow quite tight limits to be placed on the pressure, temperature and water activity conditions which accompanied metamorphism of natural osumilite occurrences in Nain, Namaqualand, and Rogaland. At fixed bulk composition, the initial melting of pelites by dehydration of biotite can occur via univariant, divariant or trivariant equilibria depending upon pressure of metamorphism. Of particular interest is that, for low pressures or more magnesian bulk compositions, fluid-absent melting begins by generating liquid from the high-variance assemblage biotite+cordierite+K-feldspar+ quartz. This type of modelling allows investigation, at least qualitatively, of the fine scale details of melt production as a function of changes in pressure, temperature and bulk composition.