Environmental determinants of fish abundance in the littoral zone of gravel pit lakes

The type and extent of habitats along the shoreline specify the distribution of fish in the littoral zone of lakes, but effects are likely species and size-specific and might be overwhelmed by lake-level environmental factors that drive fish abundance (e.g. trophic state). We applied a replicated transect-sampling design by electrofishing assessing fish abundance and distribution along the banks of 20 gravel pit lakes in Lower Saxony (Germany). Boosted regression trees were used to analyse the impact of different characteristic habitat types (e.g. vegetated, woody or open water zones), shoreline water depth and lake-level environmental variables on species-specific fish abundances. In contrast to earlier studies, lake-level environment and transect-level habitat type similarly influenced the abundances of differently sized fish species in the littoral zone of gravel pit lakes. The abundance of almost all fish species increased with lake productivity and extent of structured littoral habitats, mostly following non-linear relationships. Our work suggests that investments into the quality of littoral habitat, and not merely the control of nutrient inputs or other lake-level environmental factors, can promote abundance of most gravel pit lake fish species, in particular those who depend on the littoral zone for at least part of their life-cycle.

Most lake fish use the littoral zone on diel, seasonal or ontogenetic scales (Hofmann & Fischer, 2001;Amundsen et al., 2003;Westrelin et al., 2018) for feeding, spawning, nursery or as refuge habitat (Hölker et al., 2002;Lewin et al., 2004;Winfield, 2004). Yet, comparative studies across lakes, and time series analysis from individual lakes, have shown that lake-level environmental characteristics, in particular trophic state, morphology or water clarity, have stronger structuring effects on most lake fish communities than the quality and quantity of the habitats present in the littoral zone (Persson et al., 1991;Jeppesen et al., 2000;Diekmann et al., 2005;Mehner et al., 2005;Lewin et al., 2014). However, the characteristics of littoral habitats have been found to co-determine the abundance of selected fish species in the littoral zone (Fischer & Eckmann, 1997;Brosse et al., 1999;Helmus & Sass, 2008;Lewin et al., 2014), because availability and quality of littoral habitats may fundamentally constrain certain life-history stages (Scheuerell & Schindler, 2004;Ahrenstorff et al., 2009).
In cultural landscapes, artificially created lentic water bodies, in particular gravel pit lakes, have become common landscape elements (Soni et al., 2014;Blanchette & Lund, 2016). These artificial water bodies are often characterized by steep slopes, sandy habitats and high water depth (Gee, 1978;Blanchette & Lund, 2016;Nikolaus et al. 2020). Correspondingly, the littoral zone of gravel pit lakes appears functionally simplified relative to those of natural lakes (Emmrich et al., 2014). An increased relative importance of littoral zones has been reported to affect the abundance of selected fish species in lakes with limited littoral zones (Gasith, 1991;Hampton et al., 2011;Lewin et al., 2014). Hence, the relevance of littoral habitat features for driving the fish abundance in the limited littoral zone of gravel pit lakes might overrule the effect of environmental lake-level characteristics. However, this has not been quantified, yet.
Littoral habitat structures mainly encompass wood, coarse woody debris, reeds and submerged macrophytes, which are known to affect the abundances of fish species in the littoral zone of lakes (e.g. Okun & Mehner, 2005;Sass et al., 2012;Lewin et al., 2014). For example, after controlling for lake-level impacts rudd (Scardinius erythrophthalmus (Linnaeus, 1758)), tench (Tinca tinca (Linnaeus, 1758)) and northern pike (Esox lucius Linnaeus, 1758) have been found to be more abundant in complex habitats formed by submerged and emerged macrophytes, while European eel (Anguilla anguilla (Linnaeus, 1758)), perch (Perca fluviatilis Linnaeus, 1758) and roach (Rutilus rutilus (Linnaeus, 1758)) have been preferentially found in woody habitats (Perrow et al., 1996;Lewin et al., 2004Lewin et al., , 2014. Ecological processes such as preferences for spawning substrate and foraging habitat might explain these findings. At the same time, too dense structures can also limit foraging success (Savino & Stein, 1982, 1989Diehl, 1988) and, thus, non-linear relationships with a peak at intermediate levels can be expected for the relationship of the extent of specific habitat structures and the abundance of selected species in the littoral zone of lakes.
In particular small fish face a trade-off between foraging and shelter seeking to avoid predation, which is at the cornerstone of classical ecological theory explaining habitat choice as a function of a growthsurvival trade-off (Werner & Hall, 1988;Ahrens et al., 2012). Risk-sensitive foraging is under strong natural selection, and fish have evolved a sensitive repertoire to balance shelter seeking with foraging (Ahrens et al., 2012). Highly structured habitat is usually safe as is very shallow water for small-bodied fish, and these habitats are thus preferred by larval and juvenile fishes (e.g. Brosse & Lek, 2002;Okun & Mehner, 2005). As fish grow, risk of predation by gape-limited predators decreases (Lorenzen, 2000;Gaeta et al., 2018), which might loosen the attachment of larger-bodied fish to shelter habitat (Ř íha et al., 2015). Therefore, in particular larval and juvenile fish of several species should benefit from functional littoral zones with shallow areas and structural complexity to cope with the growth-survival trade-off mentioned above by switching among foraging and refuge habitats (Brosse & Lek, 2002).
In lentic waters comparative studies of microhabitat use by fish have mainly focused on the effects of littoral substrates (e.g. sand or rocks) rather than littoral structures (e.g. submerged macrophytes or deadwood; Fischer & Eckmann, 1997;Brosse et al., 1999;Brosse & Lek, 2002;Š mejkal et al., 2014;Ř íha et al., 2015). Lewin et al. (2014) examined the relative importance of lake-level and shoreline-level environmental determinants of the abundance of fish species across German lakes. However, this study did not consider how the littoral zone might affect different size classes of fish. Abundances of small individuals should be more dependent on littoral habitat and shoreline characters (Poizat & Pont, 1996;Grift et al., 2003;Pierce & Tomcko, 2005), whereas the abundances of large individuals should be better predicted by lake-level environmental factors (Persson et al., 1991;Jeppesen et al., 2000;Mehner et al., 2005). To test this proposition, a species-and size-specific analysis was performed of fish distributions in the littoral zone (here the nearshore margin up to a maximum depth of 2 m) of small gravel pit lakes, one of the most abundant and at the same time most understudied water body type in many cultural landscapes (Saulnier-Talbot & Lavoie, 2018;Søndergaard et al., 2018).
We hypothesized that (I) the relative importance of lake-level vs. littoral-level environmental descriptors in driving fish abundance differs between species and size-classes, (II) increasing amounts of complex, littoral habitats and shallow areas positively affect abundances of small fish, and (III) the impact of different littoral structures on abundances of large fish is non-linear and species-specific, with thresholds expected in the relationship of littoral habitat and fish abundance.

Sampling location and data collection
We sampled the littoral fish abundance in 20 gravel pit lakes in the lowlands of Lower Saxony, North Western Germany, located in the Central Plains ecoregion (Fig. 1). For littoral electrofishing, we divided the entire shoreline of each lake into transects. The number of transects varied between 4 and 27 per lake. The individual transect length ranged from 30 to 244 m. All transects were sampled by boat electrofishing (FEG 8000 electrofishing device; 8 kW; 150-300 V/300-600 V; EFKO Fischfanggeräte GmbH; www.efko-gmbh.de) with one anodic hand net (40 cm diameter and mesh size 6 mm) once per year in autumn from 2016 to 2019. This configuration enabled an effective electric fishing field of about 5 m diameter. Accordingly, along each transect the activated anode was swiftly immersed every three to five meter and all immobilized fish were netted. In rare cases, the number of transects varied between sampling years for some lakes for logistical reasons, varying water levels and due to ongoing habitat management actions. In each transect, all fish were identified and total length was measured to the nearest mm.
The transect-based sampling design was chosen to produce robust density estimates and to avoid zeroinflation of the data (e.g. Reid, 2011). Transect-based sampling further allowed the estimation of the relative composition of certain habitat types, e.g. the relative fraction of reed vs. woody habitat for a given transect. Electrofishing data were processed, and the catch per unit effort (CPUE) on transect-level was quantified for each year, lake and transect as individuals caught per 50 m electrofishing. Lake-level CPUE values were calculated as individuals caught per 50 m electrofishing per year and lake by summing up all catches and dividing them by the sum of all transect lengths. We then studied the effect of littoral structure on fish abundance (using transect-level CPUE in fish per 50 m as an abundance index), with transect-based CPUE values by size class and species nested as samples within lakes.
We confined our analysis of the relative impact of littoral structures, littoral water depth and lake-level environmental factors to six fish species, which regularly use and partly depend on the littoral zone in temperate European lakes: European eel, a relevant fisheries resource, the predators perch and pike, the benthivorous tench and the smaller-bodied ''forage fish'' roach and rudd. These species are common in German gravel pit lakes and naturally reproduce there, except for the stocked eel (Emmrich et al., 2014;Matern et al., 2019). To study the size-structured utilization of the littoral zone by fishes, all species except pike were separated into ''small'' and ''large'' size classes. For perch, roach, rudd and tench the threshold was set to 100 mm total length (TL).
Fish B 100 mm TL are particularly vulnerable to predators (Gaeta et al., 2018) and should especially use structured habitats in the littoral zone. For the larger-bodied eel, the threshold was set at 300 mm TL. Pike was the only species analysed in three size classes: small (B 200 mm TL), medium ([ 200 and B 400 mm TL), and large ([ 400 mm TL) to account for size-related habitat choice as a consequence of decreasing risk of cannibalism with increasing length (Nilsson, 2006) and size-dependent reliance on submerged macrophytes previously reported in the literature (Casselman, 1996;Casselman & Lewis, 1996).
Environmental lake-level variables Lake area and shoreline length of the gravel pit lakes were calculated using QGIS (QGIS Development Team, 2019). The shoreline development factor (SDF) was calculated after Hutchinson (1957) as an index of shoreline complexity and extent. Mean and maximum lake depth and the percentage of shallow lake areas (less than 3 m depth) were extracted from depth contour maps that were calculated by ordinary kriging in R following Monk and Arlinghaus (2017). Conductivity and Secchi depth were measured at each sampling event above the deepest point of the lake. Furthermore, water samples of the epilimnion were taken to analyse the total phosphorus concentration (TP) and the chlorophyll a concentration (Chl a) of the lakes to indicate trophic state. TP determination was conducted following the molybdenum blue method (Murphy & Riley, 1962;ISO, 2004), and the Chl a concentration was determined using high performance liquid chromatography (Mantoura & Llewellyn, 1983;Wright et al., 1991). To control for annual variation, Secchi depth, Chl a and TP values from two additional summer samplings (2017 and 2018) and one additional sampling in early spring (2017) were used to calculate a more robust mean for each lake. The extent of different shoreline structures and average transect depth were visually determined for each transect at every sampling event. We differentiated the extent (in percentage of transect length) of reeds, wood, deadwood, submerged macrophytes and the absence of structure, termed ''open littoral''. Reed habitats were mainly created by common reed (Phragmites australis) and in minor fractions by cattail (Typha sp.) and rush (Juncus sp.). The category ''wood'' was mainly represented by branches of living trees (various species) that extended from above the surface into the water column. Deadwood and roots were assigned to the category ''deadwood'' as they create a wooden structure in the water column and on the lake bottom. The category ''submerged macrophytes'' was represented by various species. For simplicity, some rare floating-leaved macrophytes, mainly water lilies (Nymphaea sp.), were also included in the category ''submerged macrophytes''. The fifth category ''open littoral'' represented littoral areas without any structures and open, mostly soft bottom or sandy substrates such as sandy beaches or unstructured angling sites. The average fished depth in each transect was estimated and noted as ''littoral depth''. Transects with a littoral depth of more than 2 m were removed from the analysis due to a reduced catchability of the electrofishing device (Zalewski & Cowx, 1990). The variable littoral depth contained three levels: very shallow (VS; \ 0.5 m); shallow (S; 0.5-1.0 m) and deep (D; 1.0-2.0 m).

Data analysis
We first conducted a principal component analysis (PCA) of the z-transformed lake-level variables mean lake depth, maximum lake depth, share of shallow lake area (0-3 m depth), TP, Chl a, Secchi depth, SDF and conductivity to reduce the dimensionality of the environmental data (Table S1 and Fig. S1 in Supplementary Information). The broken stick method was chosen as stopping rule in the PCA (Jackson, 1993). Only one relevant PC axis was retained, which displayed information about the lake's productivity. SDF and conductivity loaded on separate axes (Table S1 and Fig. S1 in Supplementary Information). PC scores from PC 1 were extracted and used for further analysis (Table 1). Two gravel pit lakes showed ''unusual'' productivities (one very high, the other very low), however, both lakes were kept for further analysis to cover a larger productivity gradient.
To model the influence of environmental predictors (lake-level and transect-level) on the abundances of the different size classes of the six fish species, we used boosted regression trees (BRTs), similar to Lewin et al. (2014). BRTs are a machine-learning technique that produces several simple models (trees) and combines (boosts) them to produce a model with an improved predictive performance (Elith et al., 2008;Buston & Elith, 2011). BRTs are superior to other statistical methods (e.g. generalized mixed models and generalized additive models) in identifying the relative importance of different predictor variables when relationships are non-linear (Elith et al., 2006;Leathwick et al., 2006). Importantly, they are able to detect non-linear effects of the predictor variables on the dependent variable (Elith et al., 2008). BRTs were modelled separately by species and size class. We also used linear mixed effects models and found qualitatively similar results, and for the sake of space and comparability with Lewin et al. (2014) opted to only present the BRT results here.
Lake-level predictor variables included PC scores representing productivity and raw values of SDF and conductivity. SDF was used in the BRTs as it represents shoreline quantity. Conductivity was integrated in the analysis as it corrects for electrofishing gear effects, i.e. lower efficiency at low conductivity, and species-specific reactivity to electrofishing at different conductivity. Transect variables describing the littoral habitat comprised the extent (in percentage of each transect) of different shoreline structures as well as average transect depth (as categorical variable with three levels). The BRTs calculated a species-and size-class specific relative importance for each variable for fish abundance. These relative importance values were scaled so that all variables summed to 100 (Elith et al., 2008) and can be compared between all lake-level and all transect-level variables to reveal the overall importance of lake vs. transect effects for the abundance of each species and size-class.
The bag fraction determines the proportion of the dataset that is used to build a single tree. It should be set between 0.5 and 0.75 to receive robust results (Elith et al., 2008). In our analysis, the bag fraction for all BRTs was set at 0.75 to ensure highest accuracy  even for datasets with a low number of observations (e.g. large pike). The tree complexity was set to 5 to account for potential interaction effects in all BRT models. The learning rate was adjusted by species and size class to ensure at least 1000 trees per BRT model as recommended by Elith et al. (2008). The BRT output was displayed by using partial dependence plot to show the effect size of a variable on the response after accounting for the average effect of all other variables (Elith et al., 2008). Fish abundances in neighbouring transects within a lake might be spatially autocorrelated violating assumptions of independence of data within lakes. Therefore, the spatial autocorrelation index Moran's I was separately estimated by fish species, lake and sampling event. Due to multiple comparisons, Pvalues were Š idák-corrected (Šidák, 1967). We found no evidence for spatial autocorrelation with just three detected events out of a total of 338 analyses (Table S2 in Supplementary Information) and thus rejected the assumption of autocorrelation. Each lake-year was assumed an independent sample as we were not interested in specific lake-level effects, similar to Lewin et al. (2014).
All statistical analyses were performed in R version 3.6.1 (R Core Team, 2019). The PCR was conducted using the package vegan (Oksanen et al., 2019), Moran's I was calculated using the package ape (Paradis & Schliep, 2018), and the BRTs were modelled using the packages gbm (Greenwell et al., 2019) and dismo (Hijmans et al., 2017). BRT bootstrapping and visualization of the BRT results was conducted using modified commands of the ggBRT package (Jouffray et al., 2019).

Lake environment and descriptive information on sampling outcomes
The gravel pit lakes were on average 7.1 ha in size with the smallest lake being 1 ha and the largest lake being 19.5 ha ( Table 1). The SDF varied between 1.1 and 2.2 with an average of 1.5. Mean lake depth ranged between 1.7 and 11.9 m with an average of 4.4 m, and the shallow areas comprised between 8.6 and 67% of the total lake area. The TP concentration ranged between 9 and 183 lg l -1 (average 30 lg l -1 ), and the Chl a concentration ranged between 2 and 59 lg l -1 (average of 13 lg l -1 ). The Secchi depth varied between 0.5 and 4.4 m (average 2.2 m) and the conductivity between 138 and 1004 lS cm -1 (average 457 lS cm -1 ).
Open littoral was the dominant shoreline structure with a share of 35.1% over all lakes and sampling events (Table 2). Reeds and wood were also common in the littoral zone with relative frequencies of 22.9% and 24.2%, respectively. Submerged macrophytes were only detected on 13.2% of the fished shoreline length, and deadwood represented the scarcest structure covering only 4.5% of the shorelines.
In total 65,261 individuals of 25 fish species were caught during 80 electrofishing surveys. The catch of all six fish species of interest summed up to a total of 53,853 fishes (Table 3). Small rudd (\ 100 mm TL) was the most abundant species and size class with 25,293 individuals, while large pike ([ 400 mm TL) were least abundant with 130 individuals, respectively. In eight gravel pit lakes, all six fish species were caught, while in three lakes only one or two of these species occurred (Table 4). Tench were caught in most of the lakes (90%), followed by perch (85%) and eel (80%). Pike and roach were each caught in 15 gravel pit lakes (75%), while rudd occurred least frequently (70%). Accordingly, the number of lake samples used in the models varied by species and size class with a minimum sample size (lake sampling events) of 31 for large roach and a maximum of 64 for small perch (Table 3).
Highest median CPUE values (number of fish per 50 m) were revealed for small perch with 4.8 (range 0.1-36.1) on lake level and 2.6 (range 0-123.8) on transect level (Table 3). The highest single species CPUE values were recorded for small rudd (lake level CPUE: 317.1; transect level CPUE: 779.1). Large pike were least abundant in the electrofishing catches with a median CPUE of 0.01 (range 0.01-0.56) on lake level and 0 (range 0-1.9) on transect level. In general, the species CPUEs for larger size classes were lower than for smaller size classes, except for stocked eel.
Importance of lake-level vs. transect-level variables Aggregated lake-and transect-habitat variables both showed similar relative importance for most species  for each year and are presented on lake level and transect level and size classes (Table 5), with few exceptions: The relative importance of lake-and transect-variables was 46.7% and 53.3%, respectively, for small roach, while for large roach the lake-variables (81.6%) were much more influential for the littoral abundance than transect-variables (18.4%). This pattern was strongly   driven by the high relevance of lake productivity (61.7%; Fig. 2) for large roach. For stocked eel, the relative importance of lake-variables also exceeded those of transect-variables with 70.4% for small eels and 61.8% for large eels. For pike the relative importance of the transect-variables compared to lake-level variables increased with size classes (small: 46.8%; medium; 49.8% and large: 57.3%), and the importance of lake-level variables decreased accordingly (small: 53.2%; medium: 50.2% and large: 42.7%).

Lake-level effects
In each of the species and size classes, one of the lake-level variables had the highest relative importance in affecting fish abundance (Figs. 2, 3 and 4). In all cases, except for small perch, productivity or conductivity most strongly explained the observed abundances. The abundance of small perch was best explained by the SDF, a variable measured on lakelevel, but describing the quantity of littoral habitat (Fig. 2). With the exception of large rudd, productivity was of high relative importance for the three cyprinid species roach (small: 25.4%; large: 61.7%), rudd (small: 25.5%) and tench (small: 22.4%; large: 38.7%). In all cases, abundances of these three species increased at higher productivity levels. Productivity was also of highest relative importance for large perch (18.8%), small pike (28.9%) and medium sized pike (26.6%). For large perch and mediumsized pike increasing productivity led to increases in abundance, while abundances of small pike and also small perch peaked at an intermediate level of productivity. Conductivity positively affected CPUE and appeared of high relative importance for abundances of eel (small: 20.9%; large: 28.6%), small tench (30.7%), large rudd (34.3%), and large pike (26.4%). A relatively high importance of conductivity was also detected for perch (small: 16.7%; large: 13.3%), but without a clear positive or negative effect on the fish catches.

Importance of specific littoral habitats
Relative importance, effect size and direction of effects of the different habitat structures varied for fish species and size classes. At transect level open littoral best described the abundance of small perch (19%), small rudd (24.3%), eel of both size classes (small: 13.1%, large: 14.4%) and large pike (17.7%). Increasing shares of unstructured, open littoral correlated with decreasing abundances of all species, except small eel. Wood was the most important habitat variable for explaining the abundance of small tench (12.9%), small and large roach (17.6% and 7.6%). High shares of over 50% wood in the littoral zone were positively related to the abundance of small and large roach, while abundances of small tench decreased at over 70% wood within a transect. A high relative importance of deadwood was found for explaining the abundance of large perch (14.9%) and medium sized pike (12%). Shares of 40% deadwood on the shoreline had the largest effects, but the number of transects with high deadwood abundances was low (Figs. 2 and 4). Reeds were found an important habitat variable for the abundance of small roach (15.7%), large eel (10.3%), large rudd (27.5%), as well as all size classes of pike (small: 10.9%, medium: 14.5% and large: 16.2%). Except for small roach, reeds influenced the fish abundances positively. Finally, submerged macrophytes were found highly important for the abundance of large perch (14.8%) and small pike (17.6%). Littoral depth had the highest relative importance for explaining the abundances of small roach (8.3%) and large pike (8.7%). In both cases highest abundances were detected in shallow water.

Non-linearities and threshold effects
There were several non-linear responses of fish abundance to habitat features and lake-level variables. In particular, non-linear relationships among a variable and fish abundance were detected for productivity impacts on tench of both size classes (Fig. 3). Abundances of large perch also increased with productivity in a non-linear pattern, with a threshold value that was smaller than the one for tench.
Non-linear pattern and threshold effects were also detected for habitat variables measured at transect level. For large perch and small pike, abundances peaked at an intermediate submerged macrophytes cover of approximately 70% and declined thereafter ( Figs. 2 and 4). By contrast, deadwood was an important shoreline structure for medium sized pike and large perch reaching highest abundances at about 40% of deadwood, while even larger amounts did not further affect the abundances. For medium and large pike, abundances strongly increased in reed-dominated transects at a threshold of at least 50% reed cover.

General findings
We studied the influence of specific shoreline structures on littoral fish species by comparing the importance of littoral habitats to those of environmental lake-level variables on species-specific abundances. In support of H1, the importance of littoral habitats was found to be species-and size-specific, and in some cases (e.g. both size classes of perch), the transect variables better explained the abundances than the lake-level environmental variables. This finding shows that the littoral fish abundance was strongly influenced by the littoral environment in gravel pit lakes. We further hypothesized a positive effect of specific littoral structures especially on the abundances of small fishes (H2), which was confirmed only for small and medium sized pike. Furthermore, we found non-linear patterns with species-and size-class specific threshold effects for certain shoreline structures on fish abundance (H3), in particular vegetation, woody structures and the extent of fully unstructured habitats. Littoral depth played only a minor role for fish abundances in gravel pit lakes, indicating that in our study lakes the habitat types were dominant factors affecting transect-level fish abundance.

Lake-level vs. transect-level impacts on fish abundance
Our results revealed a high importance of both, lake environment and shoreline habitats, for the local abundance of most species and size-classes. Lakelevel environmental variables, such as morphology and nutrient level, are known to be important drivers of lake fish communities and abundances (Persson et al., 1991;Jeppesen et al., 2000;de Leeuw et al., 2003;Mehner et al., 2005;Lewin et al., 2014). However, in our work the combined effect of all littoral variables was of similar or sometimes even higher importance than that of the pooled lake environment for driving local fish abundance in gravel pit lakes. Littoral structures are known to be crucial for the lifecycle of certain lentic fish species (e.g. pike; Casselman & Lewis, 1996;Nilsson et al., 2014) and hence changes in the shoreline habitat strongly affect fishes on species and community levels as well as abundance (Whitfield, 1986;Sass et al., 2006Sass et al., , 2012Helmus & Sass, 2008;Ziegler et al., 2017). Using a similar methodological approach to ours, Lewin et al. (2014) found substantially higher importance of lakelevel variables compared to the littoral characteristics in natural German lakes (compare Table S3). Our results are different, by revealing a higher relative importance of littoral variables compared to lake-level variables in gravel pit lakes. Compared to natural lakes, man-made gravel pit lakes are often characterized by steep slopes and higher littoral depth (Gee, 1978;Emmrich et al., 2014) and thus the distribution of submerged macrophytes and reeds is often more limited (Duarte & Kalf, 1986). The amount of deadwood is also lower in gravel pit lakes than in natural lakes, due to their young age (Robichon, unpublished data). Because of these deficits in littoral structures, their relative increases may have a stronger impact in gravel pit lakes compared to natural lakes with natural shorelines studied by Lewin et al. (2014). Nevertheless, the littoral fish community structure does not differ between gravel pit lakes and natural lakes as found in previous work in the same study region (Emmrich et al., 2014).

Lake-level environmental effects
We found lake productivity as an important, but not outstanding factor influencing littoral fish abundances in small gravel pit lakes. The relevance of the lake's trophic state for fish abundance is well established in the fish ecological literature: productivity and carrying capacity for fish biomass are strongly related to nutrient levels (Hanson & Leggett, 1982;Downing et al., 1990). Abundances of cyprinid species peak in eutrophic to polytrophic lakes, because these species benefit from increasing algal biomass and related zooplankton as food source (Persson et al., 1991;Jeppesen et al., 2000;de Leeuw et al., 2003;Mehner et al., 2005). Our study mostly included oligotrophic to mesotrophic lakes and this might have downplayed the relevance of trophic state as key lake environmental factor, relative to Lewin et al. (2014). Nevertheless and in line with literature, we also saw a strong relevance of our productivity index for abundances of cyprinids, particularly roach and tench, in the littoral zone. Abundances of pike and perch have previously found to peak at a meso-to slightly eutrophic state and decrease afterwards (Persson et al., 1991;Jeppesen et al., 2000;de Leeuw et al., 2003;Mehner et al., 2005;Lewin et al., 2014). Similarly, Haugen & Vøllestad (2018) reported highest pike abundance at intermediate phosphorus concentrations of 15 lg l -1 in shallow lakes, likely because of the loss of macrophytes with increasing trophic state. We detected only small perch and small pike abundances to peak at intermediate productivity levels, while medium sized pike and large pike benefited from higher productivity levels, suggesting that the trophic states we observed were not limiting to the abundance of pike. Relatedly, Jeppesen et al. (2000) failed to find inverse relationships of productivity and abundance of pike in shallow lakes in Denmark, and Haugen & Vøllestad (2018) reported no decline of pike abundance with nutrient levels in deep lakes, indicating that pike can tolerate a wide range of productivities and might even benefit from increasing trophic state through higher prey availability (Nilsson et al., 2009). Overall, our lakes did not consistently cover the full productivity gradient from oligotrophic to polytrophic states and, thus, we might not have covered the threshold of decreasing abundances with elevated trophic state.
We found SDF to be of major importance for perch abundances. The SDF broadly describes the availability of littoral habitats. The littoral zone displays an important habitat for fishes, especially in deep lakes, by providing food (Hampton et al., 2011) and shelter from predators (Beauchamp et al., 1994;Stoll et al., 2008). Perch, the dominating fish species in mesotrophic natural lakes (Persson et al., 1991;Mehner et al., 2005) and gravel pit lakes (Matern et al., 2019), highly rely on the littoral zone during their ontogeny (Amundsen et al., 2003). Furthermore, SDF is an important factor determining growth of perch as revealed for some of our study lakes (Höhne et al., 2020). Hence, our findings are in accordance with perch literature and show the SDF as good surrogate to explain not only growth (Höhne et al., 2020), but also abundance of perch.
Conductivity was identified as important lake-level variable especially for the catch of large pike, rudd and eel, but also of small tench. Generally, the catches of all species and size classes were positively related to conductivity. Conductivity is known to influence the size of the electric field created through the electrofishing device, with larger efficient field sizes at higher conductivity (Bohlin et al., 1989;Zalewski & Cowx, 1990) and hence increased catches. However, at low conductivity catches of larger individuals were disproportionally lower and probably the abundances were also underestimated in low conductivity lakes.

Effect of littoral structures at the transect level
We found a strong positive effect of specific littoral structures on abundances of small and medium sized pike. The observed abundances of small pike peaked in highly, but not completely macrophyte-dominated transects, which is in agreement with previous studies (Grimm, 1989;Casselman & Lewis, 1996). While low or entirely unvegetated habitats were avoided due to increased risk of cannibalism (Grimm & Klinge, 1996;Skov et al., 2003;Skov & Koed, 2004), too dense habitat structures might limit foraging success (Savino & Stein, 1989;Eklöv, 1997), explaining the detected patterns. Abundances of medium sized pike increased in less complex reeds and deadwood, which corresponds with previously reported raising independence from vegetation cover with increasing body size in this species (Chapman & Mackay, 1984;Rosell & MacOscar, 2002;Kobler et al., 2008Kobler et al., , 2009. The expected highest abundances of large fish in medium-structured habitats were only observed for large perch. Perch are superior competitors to roach in medium structured littoral habitats (Persson & Eklöv, 1995), however, highly complex and dense aquatic vegetation negatively influences foraging efficiency of percids (Savino & Stein, 1982;Diehl, 1988;Gotceitas & Colgan, 1989) and, thus, abundances (Brosse & Lek, 2002). Our results suggest that similarly to pike, large perch benefit from intermediate macrophytes coverage.
Submerged macrophytes were also highly relevant for large tench and this pattern was rather linear. Tench select for densely vegetated habitats of floating and submerged macrophytes, as well as reeds (Perrow et al., 1996;Gallardo et al., 2006;Lewin et al., 2014). Hence, highest tench abundances can be expected in habitats with highest structural complexity, which is in agreement with our results.
A further positive linear effect was detected for littoral reed stands on large rudd abundances. Aquatic vegetation has always been stated as typical rudd habitats (Eklöv & Hamrin, 1989), but Lewin et al. (2014) already revealed a higher relevance of reed habitats compared to submerged macrophytes for rudd abundance. Our results confirm this finding for rudd larger 100 mm.
Open littoral was found to negatively affect the abundance of small rudd and large eels. Lewin et al. (2014) found similar eel abundances in structured and unstructured habitats, but did not investigate size effects. River studies on Anguilla species showed an increased importance of riparian cover for eels larger 300 mm, while smaller individuals selected for diverse sediments as shelter habitat (Jellyman et al., 2002;Glova et al., 2010). Hence, the absence of littoral cover negatively affects the abundances of larger eels, while small eel abundances might be more influenced by the diversity of the sediment, a variable not measured here. Furthermore, our findings indicate a general dependence of small rudd and large eel on littoral habitat complexity, but not specifically on a distinct habitat type.
Importance of littoral water depth was generally low throughout all species and size classes compared to littoral structures. Generally, shallow water zones are often stated as habitat for fishes because they warm quicker and offer refuge from larger bodied predators (e. g. Paterson & Whitfield, 2000), and Brosse et al. (1999) detected a higher species-specific importance of littoral depth. However, Brosse et al. (1999) conducted electrofishing from June to mid-September, while our sampling was conducted from end of August to mid-October. Hence, littoral depth might be particularly important for larval fishes, and the importance of littoral depth decreases with increasing fish size.

Special case of stocked eel
Pattern and variances in eel catches were best explained by lake-level rather than transect-level variables. Lewin et al. (2014) already detected a high importance of the lake-level variables ''trophic index'' and ''mean lake depth'' on eel abundances in shallow lakes. In shallow lakes eel abundances and growth rates are usually highest at high trophic levels (Anwand, 1982;Lewin et al., 2014). In comparison to that, in deep lakes shoreline development factor and shoreline structure become more important for eel . However, eel populations in unconnected gravel pit lakes completely rely on stocking through recreational-fisheries management (Emmrich et al., 2014;Matern et al., 2019), which can positively impact their abundances if the environment offers suitable conditions (Simon & Dörner, 2014;Arlinghaus et al., 2015). Our analysis did not consider stocking effects on lake-level and, therefore, does neither allow for an accurate interpretation of the lake variable effects nor for a comparison of lake environment vs. littoral habitat effects as these results may be biased by individual stocking actions per lake. Nevertheless, the comparison of the different shoreline structures and their specific effects on the eel abundances remain valid.

Study limitations and future research needs
In our study, all data were collected by electrofishing, which is known for being size and species-selective and exhibiting lower efficiency in deep water and open habitats. However, it is still the most broadly applicable standard method for fish sampling in littoral habitats (Mueller et al., 2017). We analysed conductivity effects to account for lake-specific catchability. Nevertheless, the abundances of large fish in general and of all fish in the unstructured, open littoral habitats might have been underestimated. However, even if absolutely underestimated, at higher stock densities relatively more large fish get captured and also more fish in the less structured habitats should have been caught, such that the relative differences among habitats should broadly hold.
In our study, we exclusively used data from daytime electrofishing samplings in autumn (end of August until mid-October). Littoral habitat utilization of fishes differs between day and night (Copp, 2010;Ř íha et al., 2015) and between seasons (Kobler et al., 2008;Nakayama et al., 2018;Westrelin et al., 2018). Hence, fish samplings over all seasons and both times of the day might reveal other patterns than we presented here. Future research focusing on these aspects would complete the knowledge of the littoral zone to verify its importance for the fish abundance in gravel pit lakes. As a final limitation, we have to mention that we did not fully randomly select study lakes such that our inferences are confined to the environmental gradients characteristic of our study lakes. However, in the study region, similar to many other cases in central Europe, most gravel pit lakes are mesotrophic (Emmrich et al., 2014;Zhao et al., 2016;Søndergaard et al., 2018) and match the environmental characteristics of our study lakes. Therefore, our results remain transferable to most gravel pit lakes in Europe and might also hold for natural lakes of similar ecological conditions.

Conclusions and implications
Our results suggest that the availability and characteristics of littoral habitats are highly relevant and shape the abundance of littoral fish in gravel pit lakes, and in certain cases (e.g. perch of both size classes) their importance can even exceed those of environmental lake-level effects. We further demonstrated that fish-habitat associations are dependent on the species, size-class and the type of littoral habitat, often revealing non-linearities and threshold effects. Hence, we recommend that fisheries managers should not only rely on stocking and harvest regulations, but also consider habitat enhancement in the littoral zone of lakes as a suitable means to manage fish abundance in lakes (Sass et al., 2017). Generally, increasing amounts of littoral structure lead to increasing fish abundances for most species, with exceptions of small pike and large perch, which benefit from intermediate coverage of submerged macrophytes. We can conclude that diverse shoreline structures and the availability of various habitat types support a diverse fish community (e.g. Werner et al., 1977;Dustin & Vondracek, 2017) and elevate fish abundance in gravel pit lakes. These results are likely to hold for natural water bodies with degraded shorelines as well.