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, R² (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.

The diffuse attenuation coefficient of photosynthetically active radiation is an important inherent optical property of the subaquatic light field. This parameter, as a measure of the transparency of the medium, is a good indicator of water quality. Degradation of the optical properties of water due to anthropogenic disturbances is a common phenomenon in freshwater ecosystems. In this study, we used four algorithm-based Landsat 8 OLI and Sentinel-2A/B MSI images to estimate the diffuse attenuation coefficient of photosynthetically active radiation in Lake Villarrica located in south-central Chile. The algorithms’ estimated data from the ACOLITE module were validated with in situ measurements from six sampling stations. Seasonal and intralake variations of the light attenuation coefficient were studied. The relationship between the diffuse attenuation coefficient of photosynthetically active radiation, meteorological parameters, and an optical classification was also explored. The best results were obtained with QAA v6 KdPAR Nechad (R² = 0.931, MBE = 0.023 m⁻¹, RMSE = 0.088 m⁻¹, and MAPE = 35.9%) for spring and QAA v5 Kd490 algorithms (R² = 0.919, MBE = −0.064 m⁻¹, RMSE = −0.09 m⁻¹, and MAPE = 30.3%) for summer. High KdPAR values are associated with the strong wind and precipitation events suggest they are caused by sediment resuspension. Finally, an optical classification of freshwater ecosystems was proposed for this lake. The promising results of this study suggest that the combination of in situ data and observation satellites can be useful for assessing the bio-optical state of water and water quality dynamics in Chilean aquatic systems.