This work described a methodology of fertility modeling that integrates chemical properties of the soil through the Monte Carlo method. The chemical attributes are spatialized, along with their estimation uncertainties, by a geostatistical procedure of random simulation. Uncertainties of the chemical properties are propagated to the resulting fertility model allowing the generation of fertility maps conditioned to preset risk levels. The method presented here is illustrated by a case study of fertility for soy culture, at the State of Santa Catarina, Brazil, considering the following chemical properties of the soil: exchangeable aluminum, potassium and cation exchange capacity.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