The patterns of spatial variation of diatom assemblages from surface sediments in Lake Lama were quantified using a combined approach of ordination and geostatistics. The aims were (i) to estimate the amount of variation between diatom assemblages within the lake, (ii) to model the spatial variability of the diatom assemblages and their diversity, and (iii) to map the diatom distributions in the lake. A correspondence analysis (CA) separated the diatom assemblages into a planktonic and a periphytic group. Rheophilic taxa were found within the periphytic group. Variogram analysis showed that only the sample scores of the first CA axis and the Shannon diversity index were spatially structured. The range of spatial correlation was estimated to be 55 km for both variables. The diversity and, to a lesser extent, the sample scores had considerable small-scale variability of about 20 and 3%, respectively. Estimates of the first component of the CA and the Shannon index were derived using block-kriging. The maps of the estimates provided a basis for partitioning Lake Lama according to the spatial structures into an eastern and a western basin, a north–south connection between the basins, and a north–south directed tip at the far eastern end. It was shown that variation in diatom assemblages is mainly spatially structured at the catchment scale and that there is a considerable amount of variation at smaller scales. According to the modeled spatial distribution, the assemblages are most likely affected by the lake size, morphology, and the water and nutrient input introduced by rivers. This has to be taken into account when paleolimnological interpretations are drawn from records of complex lake systems like Lake Lama.
