In this article we derive a macroscopic model for the time evolution of root density, starting from a discrete mesh of roots, using homogenization techniques. In the microscopic model each root grows vertically according to an ordinary differential equation. The roots growth rates depend on the spatial distribution of nutrient in the soil, which also evolves in time, leading to a fully coupled non-linear problem. We derive an effective partial differential equation for the root tip surface and for the nutrient density.