Twenty-four lakes (Table 2) <16 m deep in the Kangerlussuaq (67° 00^′ N, 50° 43^′20^′′ W) region located 0.4 to 48.3 km from the Greenland ice sheet (GrIS; Fig. 1) were sampled in July 2021. The transect characterised expected spatial trends in environmental and limnological conditions (e.g. glacial inflow versus aridity, depth, size) in the region. The bedrock of the area is Precambrian gneiss. The study lakes are predominantly hydrologically isolated with no inflows and negligible input from groundwater or snowmelt (Anderson and Brodersen, 2001; Johansson et al., 2015), except those close to the GrIS, which are, in some circumstances, glacially fed (Burpee et al., 2018). The lakes furthest from the GrIS in Kellyville (Fig. 1) are oligohaline due to low precipitation, high rates of evaporation and locally derived salts (primarily an aeolian input). Annual precipitation is low (<250 mm y⁻¹; Mernild et al., 2018; Box et al., 2023) and the mean annual temperature between 1961 to 1990 was -5.7 °C (DMI, 2026). For the period July 2020 to July 2021, mean annual temperature was -2.2 °C with a maximum of 21 °C and a minimum of -34.6 °C. Mean January air temperatures in 2021 were -16.5 °C and mean July temperatures were 11.1 °C (DMI, 2026). Conditions close to the ice sheet, at higher elevations, are on average 2–3 °C cooler and with at least 10 mm more precipitation per month (Fowler et al., 2020).

Water column profiles of chlorophyll-a, dissolved oxygen (DO), pH, temperature and electrical conductivity (EC) at the deepest point of the lakes were measured using a YSI EXO2 multiparameter probe, which was set to log at 2 s intervals. The presence of an oxycline was assessed by subtracting the DO of bottom waters from the DO of surface waters. For pH, temperature, EC and chlorophyll-a, an average value for the top 2 m was calculated to represent values for the surface mixed layer and where ostracods were sampled. Secchi depth to determine water clarity was measured using a Secchi disk. Water samples for the determination of nutrient concentrations were collected from the same location using HCl-washed bottles. Samples were filtered with Whatman GF/F 0.7 µm filters, or with 0.4 µm polypropylene filters for silicate. Ostracods were collected in a 250 µm mesh zooplankton net from the littoral zone by sampling the top ∼1 cm of sediment. Where submerged macrophytes were present, samples were collected with the net from amongst the vegetation and included sampling the top ∼1 cm of sediment. Water colour, substrate and the dominant vegetation were noted from the littoral zone. Benthic images and videos were taken using an underwater camera and used to identify if submerged macrophyte cover was evident across the littoral zone. These videos were also used to visually determine water colour as clear, clear/green, green/brown and brown. According to the Nutrient Colour Paradigm, nutrient concentrations in the lakes are too low for colour to be classified as green. However, to visually distinguish between low nutrient clear water and more turbid water of “green” colour, these lakes have been visually classified as “clear/green”. Lake area was calculated at a later date using the polygon function in Google Earth.

Spatial patterns in species composition were determined by grouping the 24 lakes into clusters according to their similarity in ostracod assemblages based on Ward's method hierarchical clustering using the Pheatmap package in R version 4.4.1 (Kolde, 2025). The Corrplot package in R version 4.4.1 (Wei and Simko, 2024) was used to test for any statistically significant correlations between environmental variables. Any variables that were significantly correlated were excluded from further analysis. Multivariate correlations between ostracod species and environmental variables were examined using Redundancy Analysis (RDA). Variance partitioning analysis (VPA) was undertaken using the Vegan package in R version 4.4.1 (Oksanen et al., 2025). For inclusion in statistical tests, the dominant submerged macrophytes, presence of an oxycline, water colour and substrate were converted from categorical data to numeric.

For the 24 lakes, the dominant substrate was organic lake mud with low growing isoetid form macrophytes being the most common macrophytes. In seven lakes, no aquatic plants were observed. Water colour was categorised through observations as clear in four lakes, clear/green in nine lakes, brown in seven lakes and brown/green in three lakes (Table 2).

The electrical conductivity (EC) of the lakes ranged from 0.01 to 4.1 mS cm⁻¹, with the highest EC in lake 7 and lowest in lake 15 (Table 3). Water temperature ranged from 9.5 to 10.5 °C at lake 1 and lake 24. Total Alkalinity ranged from 0.2 meq L⁻¹ in lake 15 to 20.9 meq L⁻¹ in lake 7. pH ranged from 7.7 in lake 13 to 9.7 in lake 19. All lakes were ultra-oligotrophic to mesotrophic; soluble reactive phosphorus (SRP) ranged 1.9 to 49.7 µg L⁻¹ from lake 9 to lake 3 and nitrate (NO₃-N) concentrations were highest in lake 12 at 12.3 µg L⁻¹. The concentrations were below detection limit in lakes 1, 2, 7, 14 and 16. Chlorophyll-a ranged from 0.06 µg L⁻¹ in lake 22 to 2.65 µg L⁻¹ in lake 13.

Secchi depth, altitude, silicate, total dissolved inorganic nitrogen, nitrite, ammonium, water temperature and the macrophyte cover were all significantly correlated (at p< 0.001 or p< 0.01 and an r> 0.42) with at least two other variables (Table S1 and Fig. S1 in the Supplement).

Thirteen species of ostracods were recorded across the study lakes (Figs. 2, 3). Other than Limnocythere friabilis (Benson and MacDonald, 1963) and species of Ilyocypris, some individuals with soft parts were collected (Table S2) and so were assumed to be living at the time of collection. The most abundant species was Cypris pubera (O.F. Müller, 1776), with a maximum abundance of 87 valves per gram in lake 17. Candona candida (O.F. Müller, 1776) was present in the most sites (n=20) with Ilyocypris gibba (Ramdohr, 1808) present in the fewest (n=1). Shells (carapaces and valves) of Candona juveniles were present in ten sites, all of which also had adult C. candida individuals present. Cypridopsis vidua (O.F. Müller, 1776), Ilyocypris bradyi (Sars, 1890), Ilyocypris gibba, Sarscypridopsis aculeata (Costa, 1847) and Potamocypris parva (Schmidt, 1976) were rarely recorded, occurring in one to six sites at an abundance of no more than 3, 2, 2 and 2 valves per gram, respectively. Arctocypris sp., Heterocypris incongruens (Ramdohr, 1808), L. sanctipatricii (Brady and Robertson, 1869) and L. inopinata (Baird, 1843) were somewhat abundant with maximum abundances of 7, 15, 15 and 11 valves per gram, and occurring in four to seven sites. Ostracods were not recorded in three sites (lake 1, 4 and 16). The highest diversity of eight species was recorded in lake 6. Lowest diversity of one species was recorded in lakes 13, 14, 20 and 24. At each of these lakes, the only species recorded was C. candida.

All ostracod species except Arctocypris sp. were most abundant at low total alkalinity (≤5.4 meq L⁻¹; Fig. 4a). Cypridopsis sp. and S. aculeata were not present at total alkalinity values above this. Cypris pubera and C. candida were the only ostracod species abundant above chlorophyll-a concentrations of 2 µg L⁻¹, with only Candona juveniles also present (Fig. 4b). Heterocypris incongruens and Candona juveniles were abundant between 0.2 and 0.8 µg Chl-a L⁻¹. Heterocypris incongruens, Cypridopsis vidua, Cypris pubera, S. aculeata, I. bradyi, I. gibba, P. parva and Candona juveniles were not present below Chlorophyll-a concentrations of 0.2 µg L⁻¹. Candona candida and C. pubera were most abundant at NO₃-N concentration between 1.3 and 4.2 µg L⁻¹ (Fig. 4c). Only C. pubera, H. incongruens, C. candida and Candona juveniles were present at concentrations above this. Other than L. friabilis, which was most abundant at concentrations of 14.2 µg L⁻¹, all species were most abundant at SRP concentrations of <6.8 µg L⁻¹ (Fig. 4d). Only three species were present at the highest SRP concentration of 49.7 µg L⁻¹ namely C. pubera, which was very abundant (at 75 valves per gram), C. candida and H. incongruens. Potamocypris parva, L. friabilis, L. inopinata L. sanctipatricii, I. bradyi, I. gibba., C. candida and Candona juveniles are most abundant at depths >9 m (Fig. 4e). Cypridopsis vidua was only present at depths <6 m and C. pubera was most abundant at depths <4 m. All ostracod species were most abundant at EC <1 mS cm⁻¹, other than Arctocypris sp. (Fig. 4f). Cypris pubera, Cypridopsis vidua and S. aculeata were not present at ECs above this. Cypridopsis vidua was only present at ECs between 0.15 and 0.16 mS cm⁻¹. Lakes 5, 6, 7, 8, 9 and 10 are oligohaline and nine species were present at EC >3 mS cm⁻¹; of these L. friabilis L. sanctipatricii, L. inopinata, I. bradyi, I. gibba, C. candida and Arctocypris sp. were present at the EC of 4.09 mS cm⁻¹. At the lowest EC of 0.01 mS cm⁻¹, only C. candida, Candona juveniles and Arctocypris sp. were present. Both species of Ilyocypris were only present at ECs above 0.54 mS cm⁻¹ and L. inopinata and L. sanctipatricii at concentrations >0.32 mS cm⁻¹ (Fig. 4f).

Based on ostracod assemblages, three clusters of lakes were identified (Fig. 5). Cluster 1 (4 lakes – 17, 19, 3, 10) is characterised by high abundance of C. pubera. Cluster 2 (twelve lakes) is characterised by high abundance of C. candida. Cluster 3 (five lakes – 13, 15, 21, 6, 9) is characterised by a more diverse ostracod fauna with high abundances of Candona juv, H. incongruens, Arctocypris sp., L. inopinata, L. sanctipatricii and P. parva with intermediate abundance of C. candida. Cluster 3 lakes also have the lowest abundances of C. pubera and Cypridopsis vidua. A fourth cluster (three lakes – 1, 4 and 16), not depicted on Fig. 5, is characterised by the absence of ostracods.

RDA axis 1 is negatively correlated with EC and positively correlated with lake area and substrate and explains 16.23 % of variation (Fig. 6). Axis 2 is positively correlated with presence of an oxycline and nutrients and is negatively correlated with the presence of submerged macrophytes and depth. Axes 1 and 2 together explain 28.94 % of variation. Cluster 1 lakes are characterised by the lowest depths (mean 4.6 m, ranging from 1.4 to 10.7 m), higher SRP (mean 14.7 µg L⁻¹, ranging from 2.7 to 49.7 µg L⁻¹), higher NO₃-N (mean 3.6 µg L⁻¹, ranging from 1.3 to 7.8 µg L⁻¹) and higher chlorophyll-a (mean 0.82 µg L⁻¹, ranging from 0.51 to 1.27 µg L⁻¹). Cluster 2 lakes, on the other hand, are characterised by lowest mean chlorophyll-a (0.50 µg L⁻¹, ranging from 0.09 to 2.36 µg L⁻¹) and, although still high for lakes, the lowest pH (8.7, ranging from 8.0 to 9.3). Cluster 3 lakes are the deepest (mean 9.1 m, ranging from 2.8 to 15.7 m), encompass the oligohaline lakes 6 and 9 so have the highest EC (0.72 mS cm⁻¹, ranging from 0.01 to 3.00 mS cm⁻¹), lowest SRP (4.1 µg L⁻¹, ranging from 1.9 to 5.0 µg L⁻¹) and the majority of lakes (3 of 5 lakes) are within 10 km of the GrIS. In contrast, the majority of lakes in cluster 4 (2 of 3) are ≤40 km from the GrIS and have the largest average area (4.4 km², ranging from 0.01 to 12.8 km²). Lakes in cluster 4 also have the lowest average NO₃-N (1.29 µg L⁻¹, ranging from BDL to 3.88 µg L⁻¹) and lowest total alkalinity (3.1 meq L⁻¹, ranging from 1.1 to 3.0 meq L⁻¹). There were also no macrophytes present in cluster 4 lakes other than the presence of filamentous algae in lake 3.

Our results indicate a more complex pattern of nutrient distribution with lakes close to the GrIS and a spatial cluster of lakes (3 and 4) having higher concentrations of SRP and NO₃-N. Total phosphorus (TP) and NO₃-N concentrations in glacially fed lakes have been shown to be three times higher than in snowmelt fed lakes (Grider et al., 2025). Lower bioavailable nutrient concentrations in Kellyville lakes (those located furthest from the GrIS; Fig. 1) could be related to being further from the source of dust and to the reduced wind speed at >10 km from the GrIS (Heinemann, 1999) therefore decreasing dust derived P (Burpee et al., 2016; Prater et al., 2022). As expected, therefore, highest NO₃-N concentrations (12.27 µg L⁻¹) were in lake 12, which is glacially fed and located 0.37 km from the GrIS. After lakes 3 and 4, lake 12 had the third highest SRP concentrations of 14.43 µg L⁻¹. Lake 15 (GL6 in Grider et al., 2025) is located 3.71 km from the GrIS but has relatively low NO₃-N concentrations (1.75 µg L⁻¹) with high SRP concentrations (4.96 µg L⁻¹). Meltwater is therefore likely a dominant source of nutrients with N concentrations derived from atmospheric deposition on the ice sheet and P derived from geological weathering of the glacial bed (Hawkings et al., 2016). However, previous work has suggested that most of this mineralogically-derived P is not biologically available (Burpee et al., 2018).

Consequently, cluster 1 lakes, which have high SRP (mean 14.7 µg L⁻¹, ranging from 2.7 to 49.7 µg L⁻¹), NO₃-N (mean 3.57 µg L⁻¹, ranging from 1.3 to 7.8 µg L⁻¹), and chlorophyll-a (mean 0.82 µg L⁻¹, ranging from 0.51 to 1.27 µg L⁻¹) are likely to become more dominant in the Kangerlussuaq landscape in the future. Higher nitrate concentration in lakes is associated with visual water colour (brown and brown/green) and higher chlorophyll-a concentrations (0.51 to 1.27 µg L⁻¹; Fig. 6). Cluster 1 lakes are characterised by a high abundance of C. pubera (Fig. 5), which is most abundant at depths <4 m but present and still relatively abundant between 8 and 12 m. Here the species is only present in lakes with an EC <1 mS cm⁻¹ but reportedly found at salinity up to 4 ‰ (∼7.3 mS cm⁻¹; Stephanides, 1948). In paleolimnological records of higher salinity lakes in the region (SS6, Lille Saltsø and Store Saltsø), C. pubera has been considered rare, being only recorded in two lake basins across Greenland (Bennike, 2000; Bennike et al., 2000, 2010). Whilst C. pubera is considered abundant in our study, C. pubera was not present in the Kellyville “salt” lakes, suggesting high salinity significantly limits distribution in the Kangerlussuaq region.

In general, information on nutrient status of lakes is often not included when documenting ostracod species presence and abundance. In a study of three ponds in Patagonia, however, C. pubera was most abundant in the pond with highest TP of 121.4 µg L⁻¹ (Coviaga et al., 2018), suggesting an ecological preference for higher nutrient availability. A preference for waters with higher nutrient concentrations may be related to food source. Cypris pubera is omnivorous, feeding on algae, bacteria and Daphnia (Meisch, 2000; Coviaga et al., 2018). In lakes 3, 10 and 19, C. pubera was present in very high numbers (75, 38 and 67 valves per gram respectively). These lakes have high coverage of filamentous algae or large Nostoc cyanobacteria balls, colloquially named sea tomatoes. It is likely, therefore, that C. pubera is present in large numbers in these lakes due to a dietary preference.

Cluster 1 lakes are typically shallower (<4.6 m) than those belonging to other clusters. Most ostracod species that are abundant in deeper lakes are not present above NO₃-N concentrations of 4.24 µg L⁻¹ (Fig. 4c, e). Cluster 3 lakes are on average the deepest and are characterised by a diverse ostracod fauna including L. inopinata, H. incongruens, Arctocypris sp., P. parva, L. sanctipatricii, I. bradyi and C. candida. Limnocythere inopinata and C. candida are known to inhabit deep lakes, but both are also present in shallow lakes across Europe (Meisch, 2000). Both are also found in lakes across a range of salinities. Depth and EC are therefore likely not the controls on distribution for L. inopinata in this region. Total alkalinity has also been suggested as a control on L. inopinata abundance (Löffler, 1959). Indeed lakes 6 and 7 have the highest alkalinities of 20.9 and 15.6 meq L⁻¹. As suggested by Jungwirth (1979), L. inopinata is also not present in lakes with clay or gravel substrates. In the Canadian Arctic, the abundance of L. inopinata is negatively correlated with chlorophyll-a concentrations (Viehberg and Pienitz, 2017). Our results suggest L. inopinata is not present at concentrations above 0.8 µg L⁻¹ and is most abundant in Cluster 2 lake 8.

Cluster 2 lakes are characterised by the lowest mean chlorophyll-a concentrations (0.50 µg L⁻¹, ranging from 0.09 to 2.36 µg L⁻¹) and an abundance of C. candida. Candona candida is considered to be oligothermophilic (Vesper, 1975), preferring low nutrient concentrations and adults are present throughout the year in waters where the water temperature does not exceed 18 °C in the summer (Hartmann and Hiller, 1977), suggesting an upper temperature limit on adult life stage persistence. The species has a known Holarctic distribution and its life cycle preference for cooler summers would suggest abundance of the species in the Arctic; indeed, it has been shown to be abundant at higher latitudes (Alkalaj et al., 2019). Increasing temperatures are, however, likely to affect the life cycle and abundance, but not presence of this species.

Nutrient concentration has also been suggested as a control on L. sanctipatricii, which shows a preference for oligotrophic habitats (Scharf, 1981). The species is documented to have disappeared from Lake Mondsee, Austria, following anthropogenically-derived eutrophication (Danielopol et al., 1993). It is considered to be a cold-water indicator and has been found previously in Greenland (Table 1) as well as in Arctic Siberia (Wetterich et al., 2008). The presence of L. sanctipatricii is characteristic of Cluster 3 lakes 9 and 21 (Fig. 5), which are relatively deep, large lakes with relatively low nutrient concentrations (Fig. 6). Future increased water temperatures and nutrient concentrations in the region, may therefore limit the abundance and distribution of L. sanctipatricii.

It may be considered surprising that variance partitioning analysis (VPA) of nutrients (SRP, NO₃-N, Chl-a), EC, and habitat (the dominant submerged macrophytes and macrophyte cover), explained little of the overall variation in ostracod species composition (∼2.5 %). EC contributed the largest unique contribution (adjusted R2=0.035) but, overall, the variation is not explained by these three categories of variables (residuals =1.01; Fig. S2). However, for some species that are nektobenthic and large (e.g. C. pubera), provision of food and protection from predation offered by macrophyte cover may be a larger contributor to its presence and abundance than can be determined from this dataset. Due to the sampling strategy, it is also likely that variables such as pH, Chl-a, macrophyte cover and bioavailable nutrients vary within and between seasons, particularly in the late summer with longer ice-free periods (McGowan et al., 2018). Our results also provide a time-averaged “present day” living ostracod fauna, which is not an unreasonable approach in these remote environments, but it is likely that different ostracod species will be abundant in different seasons and that the measured parameters do not reflect the full range of habitat and environmental preferences.

Our record of Limnocythere friabilis is the only published occurrence of recent to living individuals outside North America. The North American species Limnocythere friabilis, considered to be a senior synonym of the extinct European species Limnocythere suessenbornensis (Horne et al., 2023), is unique to cluster 3 and only recorded in lake 7. Limnocythere friabilis is common in the Great Lakes region, which shares limnological features, which may favour L. friabilis, with the Kangerlussuaq region such as seasonal ice cover, increasing anthropogenic N and P enrichment since the 1970s CE (Nelligan et al., 2021) and are deep with an average depth of 19 m in Lake Erie. Candona candida, L. inopinata, and L. suessenbornensis occur together in interglacial records in Europe (e.g. Benardout, 2015; Marchegiano et al., 2020; Horne et al., 2023), with L. suessenbornensis regarded as a cold-water species. Limnocythere suessenbornensis is present during warm interglacial periods when, at least, UK summer temperatures are suggested to be similar or slightly warmer than today, but winter temperatures were up to 10 °C cooler (Benardout, 2015; Horne et al., 2023). There is, therefore, likely a significant winter temperature control on its life cycle and hence distribution. Our record would corroborate the requirement for significantly cooler average winter temperatures. Increased temperature and earlier ice out in western Greenland may therefore have an adverse impact on the distribution and abundance of L. friabilis.

Potamocypris parva is, to our knowledge, only recorded in the Kangerlussuaq area. It is considered to be endemic to Greenland, and specifically the oligohaline lakes within the Kangerlussuaq region. The species was first described from lakes close to Kellyville (Schmidt, 1976) and since then, it has been recorded in other saline lakes including Store Saltsø and SS6 (Bennike, 2000; Bennike et al., 2010). Our results suggest that P. parva is also present in lakes with lower EC but is most abundant (40 valves per gram) at EC up to 4.09 mS cm⁻¹. In carapace morphology it is closely similar to an African species, Potamocypris paludum (Gauthier, 1939), which has been found in European Pleistocene ostracod assemblages (Fuhrmann, 2012; Marchegiano et al., 2018). However, while P. parva has long antennal swimming setae, those of P. paludum are relatively short, suggesting that these are two distinct species (Claude Meisch, Musée national d'histoire naturelle, Luxembourg, personal communication, 22 September 2025); we have not found any specimens with preserved antennae in our material.

Sarscypridopsis aculeata is considered an indicator of saline waters. Here, the species is not present in the higher conductivity waters but is found in lakes with an EC of <0.47 mS cm⁻¹. However, it has previously been collected from lake 6 with an EC of 3 mS cm⁻¹ (L. Roberts, unpublished data) and is therefore likely still an indicator of more saline waters in this region. Previously, the species was not thought to be extant in Greenland (Bennike, 2000) with the only previous record in isolation basins formed during the Early Holocene. Ilyocypris bradyi was also not considered to be extant in Greenland with the last recorded occurrence in Store Saltsø in Kangerlussuaq during the Holocene warm period (∼7000 years BP) before going extinct in the region (Bennike, 2000). The recording of these species in the modern fauna suggests that increasing average temperature are now placing the region within the temperature tolerance of these species and they will continue to thrive.

Predictions for the Arctic are for temperature and precipitation to continue to increase in the 21st century (Hu et al., 2021; McCrystall et al., 2021). Higher temperatures in the Arctic will increase P and N loading into lakes by reactivating hydrological flows transporting soil- and dust-borne nutrients into lakes. For the limited number of glacially-fed lakes (e.g lakes 12 and 15) meltwater discharge into lakes will likely increase P and N loading associated with the release of ice-locked atmospheric deposition of N and P from glacial bed erosion (Hawkings et al., 2016). Predicted increases in precipitation will also increase N in lakes. Previous δ15N of NO₃-N in lakes from the Kangerlussuaq region suggest an addition from direct atmospheric N deposition (Anderson et al., 2017). In ice core records this has been shown to have increased over the last 50–100 years (Hastings et al., 2009). Future trends of N deposition, however, are reliant on policies to control emissions. Conversely, increased precipitation and meltwater may reduce P derived from dust by increasing the fluvial area, consequently reducing the effectiveness of aeolian erosion and transportation.

Increased NO₃-N concentrations in lakes is likely to favour C. pubera, C. candida, and H. incongruens occurrence, abundance and distribution (Figs. 4c, 6). With warmer temperatures, earlier ice out will alter the timing of, and support longer, growing seasons for phytoplankton and macrophytes, which may reduce bioavailable N due to N-uptake and denitrification (McGowan et al., 2018). During this process, P concentrations may increase due to internal loading. Increases in SRP concentrations would favour the same ostracod species (C. pubera, C. candida, and H. incongruens; Fig. 4d). The results presented here do not suggest a strong association between dominant macrophyte taxa and ostracod species distribution or abundance (Fig. 6) despite previous studies linking macrophyte cover and ostracod species occurrences (e.g. Roca and Danielopol, 1991; Roca et al., 1993; Frenzel et al., 2005). Perhaps the macrophyte species present, including sparse low form isoetids, do not offer the same habitat structure as species with leaves distributed throughout the water column (e.g. Potamogeton). However, a more systematic macrophyte survey would be needed to verify this; notwithstanding a lower diversity and abundance in the oligohaline lakes, a more diverse flora, including Potamogeton has been described in the freshwater lakes in Kangerlussuaq (Reuss et al., 2014). The implications of macrophyte persistence and diversity for future ostracod distribution are, therefore, currently uncertain.

With increased temperatures, precipitation is more likely to occur as rainfall, rather than snow. Coupled Model Intercomparison Project Phase 6 (CMIP6) experiments suggest that, across most of the Arctic, precipitation in winter will continue to have snowfall as the dominant type but with some rainfall and increasing in amount. However, in summer and autumn the dominant precipitation will be rainfall (McCrystall et al., 2021). By 2100, relative to the year 2000, there is a 422 % increase in CMIP6 predicted rainfall in winter, 261 % in spring, 71 % in summer, and 268 % in autumn. Greenland is predicted to have rainfall-dominated precipitation with 1.5 °C warming in CMIP6 (McCrystall et al., 2021). Currently, precipitation-evapotranspiration (P-E) patterns in the region are paced by ice and snowmelt-derived freshwater pulses. Increased rainfall will alter this seasonal pattern and could result in lower evaporation from lakes in Kellyville, with increased seasonal outflow already reported since 2023 from lake 6. Potamocypris parva and L. friabilis are most abundant at high EC (Fig. 4f), and therefore alterations to the P-E balance may affect the future distribution of these species, potentially restricting their distribution to coastal lakes.