O estudo objetiva modelar e mapear a variabilidade espacial do volume de um povoamento de Pinus sp. e avaliar as zonas de incertezas das inferências geradas pelo melhor interpolador geoestatístico testado. Foi realizada a krigagem ordinária testando seis modelos para estimativa do volume, sendo: Circular, Exponencial, Pentaesférico, Esférico, Tetraesférico e Gaussiano. O desempenho dos modelos foi avaliado por meio do grau de dependência espacial (GD), raiz do erro médio quadrático (RMS), erro padrão da estimativa (Syx%) e coeficiente de determinação (R²) da validação cruzada. Após a escolha do melhor modelo testado, foram geradas as zonas de incerteza das estimativas visando proporcionar o controle nos intervalos de incerteza nas estimativas preditas de forma variável espacialmente, diferente do que tradicionalmente vem sendo aplicado, com a consideração de uma incerteza média para toda a área. O modelo Circular proporcionou os melhores resultados, apresentando o maior grau de dependência espacial e menor RMS e Syx%, sendo de 1,7 m³ e 11,2% respectivamente, assim como o maior R² (0,9). A zona de certeza representou 77,79%, e as de incerteza de 22,21% da área, o que torna a distribuição espacial das parcelas amostradas confiáveis para interpolação do volume comercial.