Aquatic systems have traditionally played a key role in the development of human life, providing multiple ecosystem services to society and being a reservoir for a wide biodiversity of organisms. Among them, bacteria belonging to Legionella stand out, mainly because they are of great interest both in the field of microbial ecology and public health, since some of them turn out to be pathogenic for humans. The aim of this work was to study the monthly temporal dynamics of Legionella spp. and its relationship with the environmental variables measured in two Pampean shallow lakes (Gómez and Carpincho, Buenos Aires Province, Argentina). The analysis was carried out using a quantitative approach by real-time Polymerase Chain Reaction (qPCR) and a non-quantitative approach using bacterial diversity data obtained by Next Generation Sequencing (NGS), using the Illumina MiSeq platform. Our results showed that the overall Legionella abundance was very high in the studied Pampean shallow lakes. Notably, fluctuations in dissolved organic carbon and temperature influenced the dynamics shifts in Legionella abundances. Correlation analyses between Legionella reads from NGS and copy numbers obtained through qPCR revealed positive relationships, unveiling distinctions attributable to the diverse sequence processing algorithms employed in the analysis of NGS data.