O estado do Piauí é um dos maiores produtores de commodities agrícolas da região Nordeste do Brasil. O crescimento da produção de grãos é estimulado pelo manejo de fertilizantes que contém elementos terras raras (ETRs) na sua composição. Diante dos riscos ambientais e para a saúde, as concentrações de ETRs em solos devem ser monitoradas. Técnicas tradicionais empregadas na determinação de ETRs são onerosas, demoradas e requerem intensa preparação das amostras. Neste estudo, investigamos o potencial da espectroscopia de infravermelho próximo na predição de ETRs e a contribuição desse método no mapeamento destes elementos em grandes territórios. Foram coletadas 243 amostras compostas de solos em áreas sob nenhuma ou mínima influência antrópica. Os ETRs foram quantificados em ICP-OES. Os espectros foram obtidos na faixa do infravermelho próximo (NIR - 1000 a 2500 nm). A regressão de mínimos quadrados parciais (PLSR) foi selecionada para a calibração dos espectros. O filtro de transformação normal padrão (SNV) foi aplicado no pré-tratamento dos espectros. Os resultados mostraram que o modelo apresentou boas predições para o La, Ce, Nd, Dy e ƩETRPs e predições razoáveis para o Sm, Eu, Gd, Tb, Dy, Er, Yb, Lu, ƩETRLs e ƩETRs. A maior precisão (R2 ajust.) do PLSR foi observado para La, Ce, Nd e ƩETRPs. As concentrações estimadas por espectroscopia de reflectância difusa no NIR foram utilizadas na distribuição espacial. Dessa forma, observou-se que a espectroscopia de reflectância no NIR é um método rápido e seguro para a determinação de ETRs no solo, além de facilitar o trabalho de mapeamento de concentrações de ETRs, tarefa importante para o monitoramento da qualidade dos solos.