Resumen: Remote sensing was used as an early alert tool for water clarity changes in five Araucanian Lakes in South-Central Chile. Turbidity records are scarce or unavailable over large and remote areas needed to fully understand the factors associated with turbidity, and their spatial-temporal representation remains a limitation. This work aimed to develop and validate empirical models to estimate values of turbidity from Landsat images and determine the spatial distribution of estimated turbidity in the selected Araucanian Lakes. Secchi disk depth measurements were linked with turbidity measurements to obtain a turbidity dataset. This in turn was used to develop and validate a set of empirical models to predict turbidity based on four single bands and 16 combination bands from 15 multispectral Landsat images. The best empirical models predicted turbidity over the range of 0.3–12.3 NTUs with RMSE values around 0.31–1.03 NTU, R2 (Index of Agreement IA) around 0.93–0.99 (0.85–0.97) and mean bias error (MBE) around (−0.36–0.44 NTU). Estimation maps to analyze the temporal-spatial turbidity variation in the lakes were constructed. Finally, it was found that the meteorological conditions may affect the variation of turbidity, mainly precipitation and wind speed. The data indicate that the turbidity has slightly increased in winter–spring. These models will be used in the future to reconstruct large datasets that allow analyzing transparency trends in those lakes.