In this work, we developed a number of generalised skin diffusion based pharmacokinetic models to relate published in vivo urinary excretion data to matching experimentally generated in vitro human skin permeation test (IVPT) data for a series of topically applied salicylate esters. A simplified linear in vivo model was found to inadequately describe the time course of urinary excretion over the entire sampling period. We represented the skin barrier as both a one layer (stratum corneum) and a two-layer (stratum corneum with viable epidermis) diffusion model and convoluted their Laplace solutions with that for a single exponential disposition phase to describe the urinary excretion profiles in the Laplace domain. We also derived asymptotic approximations for the model and estimated the conditions under which they could be used. We then sought to develop in vitro - in vivo relationships (IVIVR) for topically applied methyl, ethyl and glycol salicylates using our experimental IVPT data and the literature urinary excretion data. Good linear IVIVRs for ethyl and glycol salicylates were obtained, but the IVIVR for methyl salicylate was poor, perhaps because of topical stimulation of local skin blood flow by methyl salicylate. The ratio of the hydrated to dehydrated skin permeation for all salicylate esters was the same in both the IVPT and in vivo studies. A diffusion based one compartment pharmacokinetic model was also developed to describe the urinary excretion of solutes after removal of topical products and to compare the methyl salicylate skin permeation for five different body sites. The work presented here is consistent with the development of skin IVIVRs, but suggests that different skin conditions, application sites and local skin effects may affect model predictions.