Hidrocarbonetos de petróleo são encontrados em áreas contaminadas na forma de líquido leve de fase não aquosa (light non-aqueous phase liquids – lnapl). o processo de biodegradação natural desses contaminantes é conhecido como depleção natural na zona da fonte, do inglês natural source zone depletion (nszd). a determinação das taxas de depleção natural é fundamental para avaliar a biodegradação natural do lnapl como uma solução eficaz, econômica e sustentável na gestão de passivos. estudos recentes apresentam relações diretas entre a temperatura do solo como evidência da ocorrência da nszd. nesse contexto, o objetivo deste trabalho é apresentar um novo algoritmo para determinar a depleção da massa de lnapl em função da variação de temperatura na fonte de contaminação. o algoritmo foi testado com dados sintéticos gerados a partir de um modelo analítico. para avaliar a sensibilidade do algoritmo, foi simulado o efeito de aumentos de espaçamento das sondas e intervalos de medição, bem como os efeitos de erros na medição de temperatura e umidade. para concepções com espaçamentos entre sondas de até dois metros, e com erros de medição de temperatura e umidade distribuídos normalmente com desvio padrão de 0,5ºc e 20%, respectivamente, observou-se que os resíduos de medição de calor nos 2 anos se limitaram a uma superestimativa em 10% do calor liberado. finalmente, foram comparados os resultados do algoritmo desenvolvido com os resultados dos algoritmos propostos por warren e bekins (2015) e stockwell (2015). a análise comparativa revelou que os algoritmos de warren e bekins (2015) e stockwell (2015) captaram apenas 30% e 51% do calor na zona da fonte de contaminação, respectivamente. já o novo algoritmo apresentou melhor desempenho, uma vez que considera a variação vertical dos parâmetros térmicos do meio poroso na zona da fonte. estudos futuros preveem a validação do algoritmo a partir da medição da temperatura com instrumentação de campo em áreas experimentais contaminadas.