O objetivo deste trabalho foi descrever um procedimento de modelagem de fertilidade do solo que integra propriedades químicas do solo utilizando-se do método de Monte Carlo. A espacialização das propriedades químicas foi obtida por procedimento geoestatístico de simulação estocástica, com modelagem das incertezas associadas às estimativas. As incertezas das propriedades químicas foram propagadas para o modelo de fertilidade resultante, possibilitando a geração de mapas de fertilidade condicionados a níveis de risco prédefinidos. O método aqui apresentado é ilustrado por um estudo de caso de fertilidade para cultura de soja, no Estado de Santa Catarina, considerando as seguintes propriedades químicas do solo: alumínio trocável, potássio e capacidade de troca catiônica.2004