The aim of this work is to introduce a mathematical model representing the evolution of the temperature in a vegetation cover and the ground underneath it. Vegetation, and its interaction with soil, plays a very important role in the protection of soil surface from the action of sun and precipitations. A reduction in the vegetated mass increase the risk of desertification, soil erosion or surface runoffs which which can give rise to soil loss and sediment retention. These processes can favour climate change and global warming, which are major concerns nowadays. The mathematical model presented takes into account the main processes involved in vegetation cover and the interaction with the soil, among which, we can mention the Leaf Area Index, which is a dimensionless quantity defined as the one-sided green leaf area per unit ground surface area, or albedo and co-albedo which are clearly influenced by the vegetation. It is also considered a nonlinear heat capacity in the soil which incorporates the latent heat of fusion, when the phase change takes place. The numerical technique used to solve the mathematical model is based on a finite volume scheme with Weighted Essentially Non Oscillatory technique for spatial reconstruction and the third order Runge-Kutta Total Variation Diminishing numerical scheme is used for time integration. Some numerical examples are solved to obtain the distribution of temperature both in the vegetation cover and the soil.