The skin concentration of a substance after a topical application or exposure determines both local treatment outcomes and the dermal toxicity assessment of various products. However, quantifying the time course of those concentrations at skin effect sites, such as the viable epidermal, superficial dermis and appendages in humans is especially problematic in vivo, making physiologically based mathematical modelling an essential tool to meet this need. This work further develops our published physiologically based pharmacokinetic and COMSOL based dermal transport modelling by considering the impact of the superficial subpapillary dermal plexus, which we represent as two well stirred compartments. The work also studied the impact on dermal concentrations of subpapillary plexus size, depth, blood velocity and density of subpapillary plexus vessels. Sensitivity analyses are used to define the most important transport determinants of skin concentrations after topical application of a substance, with previously published results used to validate the resulting analyses. This resulting model describes the available experimental data better than previous models, especially at deeper dermal depths.