Control over cell growth by mobile regulators underlies much of eukaryotic morphogenesis. In plant roots, cell division and elongation are separated into distinct longitudinal zones and both division and elongation are influenced by the growth regulatory hormone gibberellin (GA). Previously, a multicellular mathematical model predicted a GA maximum at the border of the meristematic and elongation zones. However, GA in roots was recently measured using a genetically encoded fluorescent biosensor, nlsGPS1, and found to be low in the meristematic zone grading to a maximum at the end of the elongation zone. Furthermore, the accumulation rate of exogenous GA was also found to be higher in the elongation zone. It was still unknown which biochemical activities were responsible for these mobile small molecule gradients and whether the spatiotemporal correlation between GA levels and cell length is important for root cell division and elongation patterns. Using a mathematical modeling approach in combination with high-resolution GA measurements in vivo, we now show how differentials in several biosynthetic enzyme steps contribute to the endogenous GA gradient and how differential cellular permeability contributes to an accumulation gradient of exogenous GA. We also analyzed the effects of altered GA distribution in roots and did not find significant phenotypes resulting from increased GA levels or signaling. We did find a substantial temporal delay between complementation of GA distribution and cell division and elongation phenotypes in a GA deficient mutant. Together, our results provide models of how GA gradients are directed and in turn direct root growth.