The second DEM was produced interpolating contour lines digitized from the same topographic maps used to produce the soil map of \citet{MiguelEtAl2012}, the geologic map of \citet{MacielFilho1990}, and the land use map of (Figure \ref{fig:contour-dem}). Interpolation of the raster surface with 5-m pixel size was performed using the same hydrologically correct interpolation method employed by \citet{JarvisEtAl2008}. The algorithm is implemented in ArcGIS software by ESRI and includes an interpolation method based on the ANUDEM program developed by \cite{Hutchinson1989}. Vector files of contour lines at 10 meters equidistance, drainage network, lakes, and peaks, were used to generate the DEM and define the boundary of the study area.