Geologic maps were registered and geocoded with the prediction grid using the nearest neighbor resampling method. This method was used to avoid changing raster values. Systematic positional errors \citep{Samuel-RosaEtAl2014} were corrected using affine transformation (R-package vec2dtransf \citep{Carrillo2012}).

SG-I

Inferior Sequence of the Serra Geral Formation. Composed mainly of basic igneous rocks (tholeiitic basalt and andesite);

SG-S

Superior Sequence of the Serra Geral Formation. Composed mainly of acid igneous rocks (granophyric rhyolite and rhyodacite);

BT

Botucatu Formation. Composed mainly of aeolian sandstones;

CT

Caturrita Formation, is which composed mainly by fluvial sandstones;

QD

Quaternary deposits of fluvial, alluvial, and colluvial origin.

Land use maps

The first land use map available was produced by manually digitizing land use data included in the topographic maps (scale of 1:) used to produce the more detailed pedological and geologic maps.

The most up-to-date land use map was prepared using high resolution orbital images (Quick Bird satellite). It covers the years of 2008 and 2009, and was prepared at a scale of 1: \citep{SamuelRosaEtAl2011a}.

Land use maps were registered and geocoded with the prediction grid using the nearest neighbor resampling method. This method was used to avoid changing raster values. Systematic positional errors \citep{Samuel-RosaEtAl2014} were corrected using affine transformation (R-package vec2dtransf \citep{Carrillo2012}).