A variabilidade dos elementos do solo vem sendo estudada por técnicas de modelagem ambiental por meio da geoestatística. Essa pode ser uma ferramenta essencial para o desenvolvimento de projetos relacionados a métodos de interpolação e metodologias estatísticas para validar a correlação entre as características referentes aos solos. Nessa perspectiva, esse trabalho tem como objetivo modelar as características químicas naturais dos elementos maiores P, K, Mg, Ca, e Na dos solos de Feira de Santana-BA. Para tal, realizou-se a coleta das amostras de solos na área de estudo e foram realizadas as análises químicas em laboratório. Posteriormente foram removidos os outliers e efetuou-se a estatística descritiva dos atributos, avaliando também o coeficiente de correlação linear de Pearson entre os elementos. Calculou-se o tamanho do pixel, definiu-se o interpolador Krigagem por meio do cálculo dos resíduos e selecionou-se os critérios para classificação das propriedades. Os semivariogramas experimentais foram ajustados e em seguida produziu-se a Krigagem Ordinária, sendo também gerados mapas 2,5D. Observou-se que na maioria das vezes, as unidades de alta distribuição de elementos maiores estão ao oeste, em contrapartida, as áreas de baixa ao leste do município.