The residence time distribution of noneliminated solutes in the liver can be represented by a variety of stochastic models. The dispersion model (closed and mixed boundary conditions), gamma distribution, log normal distribution and normal distribution models were used to describe output concentration-time profiles after bolus injections into the liver of labeled erythrocytes and albumin. The dispersion model and log normal distribution model provide the best representation of the data and give similar estimates of relative dispersion and availability for varying hepatocellular enzyme activity. The availability of solutes eliminated from the liver by first-order kinetics is determined by the residence time distribution of the solute in the liver and not on events occurring in the liver when a uniform enzyme distribution is assumed. Both enzyme heterogeneity (axial or transverse) and hepatocyte permeability may affect solute availability. A more complex model accounting for enzyme distribution and the micromixing of solute within the liver is required for solutes undergoing saturable kinetics.