Offshore wind farms and the attraction–production hypothesis: insights from a combination of stomach content and stable isotope analyses

Offshore wind farms (OWFs) act as artificial reefs, attracting high abundances of fish, which could potentially increase their local production. This study investigates the feeding ecology of fish species that abundantly occur at artificial habitats, such as OWFs, by examining the short- and the long-term dietary composition of five species: the benthopelagic Gadus morhua and Trisopterus luscus, the pelagic Scomber scombrus and Trachurus trachurus, and the benthic Myoxocephalus scorpioides. We conducted combined stomach content and stable isotope analyses to examine the short- and the time-integrated dietary composition, respectively. Our results indicated that benthopelagic and benthic species utilize artificial reefs, such as OWFs, as feeding grounds for a prolonged period, since both analyses indicated that they exploit fouling organisms occurring exclusively on artificial hard substrates. Trachurus trachurus only occasionally uses artificial reefs as oases of highly abundant resources. Scomber scombrus does not feed on fouling fauna and therefore its augmented presence in OWFs is probably related to reasons other than the enhanced food availability. The long-termed feeding preferences of benthic and benthopelagic species contribute to the hypothesis that the artificial reefs of OWFs could potentially increase the fish production in the area. However, this was not supported for the pelagic species.


Introduction
Worldwide, the number of offshore renewable energy installations is increasing fast as a measure to mitigate the effects of climate change. In the southern part of the North Sea, offshore wind farms (OWFs) are being continuously constructed (Soma et al., 2019), adding artificial hard substrates (the foundations and the associated scour protection layers) to usually naturally soft bottom areas. These wind turbines act as artificial reefs, drastically altering the local habitat (Petersen & Malm, 2006) and adding hard bottom communities to the original environment (Mangi, 2013). These communities consist of species exclusively occurring on artificial hard substrates and rarely observed on natural hard habitats (Zintzen, 2007). In the southern part of the North Sea, natural reefs are rare and are located on the Hinder Banks (gravel fields), the Cleaver Bank (gravel fields), the Texel Rough (boulder clay), and the Borkum Reef Grounds (rocky outcrops) (Veenstra, 1969;Dörjes, 1977;Coolen et al., 2015). Natural hard substrates (gravel beds) in the Belgian part of the North Sea are located at * 15 km distance from the OWFs (Houziaux et al., 2008).
Artificial habitats usually accommodate higher densities of fish compared to natural reefs (Bohnsack & Sutherland, 1985), although the design and location of artificial reefs may influence these densities (Strelcheck et al., 2005). While OWFs act as fish aggregation devices (Wilhelmsson et al., 2006), the exploitation of the locally attracted fish is either not allowed within the OWF areas (e.g. Belgium) or it is partially allowed (e.g. Netherlands) or, finally, it is fully allowed but not largely exploited by fishers (e.g. UK). Two different hypotheses have been proposed to explain the increased fish abundances in artificial reefs (Brickhill et al., 2005). The attraction hypothesis proposes that fish move from the surrounding environment towards the artificial reefs, where their local production does not increase (Brickhill et al., 2005). The attraction to artificial reefs can be explained by the increased structural complexity resulting in an enhanced protection against predators, shelter from currents (Bohnsack & Sutherland, 1985) and/or an increased food availability (Reubens et al., 2011). The production hypothesis suggests that the carrying capacity of the original environment increases due to the creation of the new habitats, favouring attraction and enhanced growth (Lindberg, 1997;Powers et al., 2003;Brickhill et al., 2005). These two hypotheses were put forward as the two extremes of a wide range of ecological processes (Svane & Petersen, 2001;Brickhill et al., 2005). Attraction could for example, potentially facilitate production, since these two actions are not necessarily mutually exclusive (Svane & Petersen, 2001;Brickhill et al., 2005;Reubens et al., 2013c;Cresson et al., 2019). For example, if attraction towards the OWFs increases fish landings, in combination with fisheries prohibition/limitation (Samples & Sproul, 1985), it could potentially lead to production.
Understanding the attraction-production hypotheses with respect to the artificial reefs at OWFs is crucial both for their efficient management and for the associated fisheries in the area (Cresson et al., 2014). So far, the attraction-production hypotheses have mainly been investigated by using stomach content analysis (Leitão et al., 2007;Reubens et al., 2011) and fish movement in close proximity to artificial reefs (Reubens et al., 2014b). However, approaching the attraction-production hypotheses by investigating the trophic relationships on artificial reefs and the transfer of organic matter from producers to consumers with the use of tracers has only been recently studied in oil and gas platforms the Gulf of Mexico (Daigle et al., 2013;Dance et al., 2018;Reeves et al., 2019).
Research on the short-term feeding ecology is often reliant on stomach content analysis (SCA), providing a snapshot of the recently ingested (up to 10 h) food items (Hyslop, 1980). SCA can provide a high taxonomic resolution (e.g. at the species level) of the ingested food items that may be difficult to attain with other methods (Lin et al., 2007). However, it also has multiple limitations, such as the requirement of high sampling frequencies in order to obtain a timeintegrated overview of the feeding habits of a species (St John, 1999;Parkyn et al., 2001;Daly et al., 2013) and a minimal indication of the assimilated dietary items (Parkyn et al., 2001;Jensen et al., 2012). Hence, stomach contents mainly consisting of hard substrate fouling organisms can be an indication of attraction towards the artificial structures for foraging activities.
The majority of the limitations of the SCA can be overcome by combining this method with stable isotope analysis (SIA). SIA has considerable advantages, showing a time-(up to 6 months) and space-integrated dietary estimation (Bearhop et al., 2004;Fry, 2006;Newsome et al., 2007), providing long-term information of a species feeding ecology. The time-integrated information can be acquired with a few samples (absolute minimum three replicates- Jackson et al., 2011) per species. A diet based on fouling fauna for a prolonged period could indicate a long-term residency close to artificial reefs, which is a clear indication of attraction, but could potentially also lead to increased local production. SIA provides information on the truly assimilated food items (Boecklen et al., 2011) and thus on the nutritionally important food sources (Melville & Connolly, 2003;Daly et al., 2013). Moreover, it allows for investigating the isotopic niche of a species, which is a proxy for its trophic niche (Jackson et al., 2011). However, this method does not provide the high taxonomic resolution as provided by SCA (Daly et al., 2013). Thus, a combination of SIA and SCA is a powerful approach that can provide a complete picture of the feeding ecology of fish species and can improve the interpretation of aquatic food webs (Parkyn et al., 2001).
In this paper, we assessed the feeding ecology of the five fish species by using a combination of SCA and SIA. We aimed at investigating whether fish use artificial reefs as feeding grounds for prolonged periods, according to three assumptions: (a) fish species use artificial reefs as feeding grounds for a prolonged period when both the SCA and the SIA show similar results and the fish' diet relies on fouling fauna that only occurs on artificial hard substrates in the southern North Sea; (b) fish occasionally use artificial reefs as feeding grounds if SCA reflects diet composition on food items derived from artificial structures (i.e. fouling organisms) which is not confirmed by SIA; and (c) fish probably do not use the artificial reefs as feeding grounds if both methods reveal a diet based on food items that are not associated with artificial hard substrates.

Methodology
Sampling area and data collection Sampling was conducted next to a gravity-based wind turbine (D6, coordinates: 51°33.04 0 N-02°55.42 0 E) within the C-Power wind farm at the Thornton Bank at the Belgian part of the North Sea (BPNS). The gravitybased foundation has a diameter of 6 m at the sea surface increasing to 14 m at the sea floor (Reubens et al., 2011), where it is surrounded by a rock collar forming the scour protection layer with a width of 62.5 m (SPL- Peire et al., 2009). We chose this particular location because a wealth of data on the investigated fish species was available (Reubens et al., 2011;Kerckhof et al., 2018).
Individuals of five species were collected during day time in August 2016: juvenile cod (Gadus morhua; length (mean ± SD): 320 ± 35 mm; weight (mean ± SD): 324 ± 100 g), pouting (Trisopterus luscus; length (mean ± SD): 237 ± 30 mm; weight (mean ± SD): 188 ± 80 g), mackerel (Scomber scombrus; length (mean ± SD): 289 ± 34 mm; weight (mean ± SD): 221 ± 66 g) and horse mackerel (Trachurus trachurus; length (mean ± SD): 254 ± 24 mm; weight (mean ± SD): 172 ± 51 g) were sampled from an inflatable boat at about 5 m from the turbine using fishing lines [hooks: Arca size 4, bait: Arenicola marina (Linnaeus, 1758)]. Sculpin (Myoxocephalus scorpioides; length (mean ± SD): 213 ± 64 mm; weight (mean ± SD): 168 ± 133 g) was collected by spear-fishing by scientific divers since this benthic species lives in between the rocks of the scour protection layer. The benthopelagic species (cod and pouting) were collected above the scour protection layer, while the pelagic species (mackerel and horse mackerel) were sampled from ca. 10 m depth of the water column. Only cod juveniles, belonging to the I-group (i.e. \ 50 cm (Reubens et al., 2013c), were sampled (from now on referred as cod). This size category is known to be abundantly frequent in the OWFs (Reubens et al., 2013c), while the older age groups (i.e. II-group and beyond) move towards other locations (Reubens et al., 2014b). Fish were immediately killed by pithing with a metal rod (Kadye & Booth, 2012) and measured (total length in mm). Subsequently, fish stomachs were extracted and preserved in an 8% formaldehyde-seawater solution for the stomach content analysis (Reubens et al., 2011). A part of the dorsal (white) muscle without skin of at least five individuals (from which we also collected the stomachs) per fish species was isolated for stable isotope analysis using a scalpel. The scalpel was cleaned with ethanol before every use (both between different species and between different individuals) to avoid cross-contamination. This tissue yields the lowest d 15 N and d 13 C variability compared to other fish tissues, such as red muscle, liver and heart, and it is considered the best tissue for use in stable isotope analysis (Pinnegar & Polunin, 1999). The muscle tissue was stored in a clean vial with filtered sea water in a -20°C freezer until further analysis. Stomachs and muscles were extracted from similar sized individuals per species (see above) to exclude any size-related dietary shifts (Mittelbach & Persson, 1998).

Stomach content analysis
For the short-term diet analysis, we obtained stomach content samples from 62 individuals in total (pouting: 17 ind.; cod: 16 ind.; horse mackerel: 12 ind.; mackerel: 12 ind.; sculpin: 7 ind.). The number of stomachs analysed was relatively small and deriving from individuals collected only in summer due to the logistic effort of collecting all consumers and food resources in all the seasons. When the small number of stomachs is unavoidable, pairing SCA with stable isotope analysis (SIA- Davis et al., 2012) can help validating trophic relationships elucidated through SCA (Clarke et al., 2005) and offer additional information on trophic flows and food sources (Dalsgaard et al., 2003;Melville & Connolly, 2003;Abrantes et al., 2014).
In the laboratory, stomachs were dissected, cleaned three times with water to assure that they were empty, and the content was collected and rinsed thoroughly with milli-Q water. Stomach content samples were drained in a 90 lm mesh sieve to thoroughly collect both the small zooplanktonic organisms and the half-digested parts of bigger organisms (i.e. crabs and fish). Consequently, they were immediately examined under a stereomicroscope to identify prey items to the lowest taxonomic level possible, and counted (Reubens et al., 2011;Rumolo et al., 2017). Since many prey individuals were incomplete (i.e. stomachs containing only the scales of fish or a cheliped from a decapod), we did not estimate the biomass contribution of prey species. Prey items were allocated into ten broader taxonomic groups (''Stable isotope analysis'' section) to facilitate the analysis. The relative numerical abundance (%), of each prey item in all non-empty stomachs of one species (Jardas et al., 2004) was estimated.
Differences in stomach content composition between the different fish species were visualized using non-metric multi-dimensional scaling (MDS) and tested using the one-way analysis of similarity (ANOSIM) routine. As the quantity of the prey items in the stomach contents varied across the fish individuals, the underlying Bray-Curtis similarity matrix (Clarke & Gorley, 2006) was constructed based on relative prey abundance data per stomach (De Crespin De Billy et al., 2000). Following a significant ANOSIM result, pairwise tests were conducted to investigate possible diet overlap. We considered pairwise comparisons with R \ 0.25 as substantial dietary overlap and with R [ 0.5 indicating relatively minimal dietary overlap, while 0.25 \ R \ 0.5 was considered moderate dietary overlap (Creque & Czesny, 2012).

Stable isotope analysis
For the analysis of the more long-term trophic ecology of the five fish species, carbon (d 13 C) and nitrogen (d 15 N) stable isotope analyses were performed (5 individuals of cod, pouting, horse mackerel and sculpin, and 11 individuals of mackerel). The frozen dorsal muscle tissue was thawed, rinsed with milli-Q water to remove salts (Nilsen et al., 2008) and dried overnight at 60°C. The dried samples were grounded to a homogenous powder using a pestle and a mortar and approximately 1 mg of dried tissue per individual per species was encapsulated in a Sn capsule (8 9 5 mm, Elemental Microanalysis UK). The capsules were stored in multi-well microtitre plates until the isotopic analysis for d 13 C and d 15 N was conducted at the UC Davis Stable Isotope Facility (University of California, USA) by a PDZ Europa ANCA-GSL elemental analyser, interfaced to a PDZ Europa 20-20 isotope ratio mass spectrometer (Sercon Ltd., Cheshire, UK). d 13 C and d 15 N data are expressed in the standard delta (d) notation relative to standard reference material: Vienna Pee Dee Belemnite for carbon and atmospheric N for nitrogen. The results are expressed in parts per thousand (%) according to the following equation: where X is 13 C or 15 N and R is the corresponding ratio of 13 C/ 12 C or 15 N/ 14 N. The trophic niches of the five fish species were estimated using the stable isotope Bayesian ellipses package in R (SIBER package, Jackson et al., 2011). SIBER generates sample size-corrected standard ellipse areas (SEA C ) that were used to estimate the isotopic niche size and the trophic overlap between the five fish species. SEA C nullifies the bias created from small sample sizes (n [ 3) and reflects the isotopic niche of a population, which is considered a proxy of its trophic niche (Jackson et al., 2011). By estimating the SEA C overlap, we obtained information about the isotopic niche overlap between the different fish species. Results are expressed as mean ± standard deviation (SD).
The relative contribution of each prey item to the diet of the different fish species was estimated using stable isotope mixing models in R (SIMMR package, Parnell et al., 2013). SIMMR can incorporate the variability that is associated with consumers, sources and trophic enrichment factors (TEFs, Parnell et al., 2013). By using Markov Chain Monte Carlo simulations, this model can determine the potential contribution of different sources to the diet of a consumer (Parnell et al., 2010). As the isotopic composition of the fish species was measured by using muscle tissues, the diet estimates represent the proportion of the sources that have been assimilated by this tissue (Lamontagne et al., 2016). This tissue has slower turnover rates than other body parts, and thus, its isotopic signature reflects a longer-term diet history than tissues with faster turnover rates (Bruestle et al., 2018). SIMMR considers the mean (l prey ) and the standard deviation (r prey ) of the stable isotope ratios of the prey species for the estimation of the proportion of a given prey in the diet (Parnell et al., 2013;Lamontagne et al., 2016). However, this model is based on the assumption that all the food resources are included in the analysis (Phillips et al., 2014). In order to meet this assumption, we used the food items that were identified by the SCA in the stomachs of juvenile cod, pouting, horse mackerel and sculpin as prior information for the SIMMR (Jensen et al., 2012). For mackerel, though, we also included uncountable prey items that were not identified in their stomachs to meet the assumption of the model that the consumers' mixtures fall within the polygon of the prey items. Food resources representing \ 1% of the stomach contents per species were not included in the mixing models. The prey items were grouped into 10 broad taxonomic categories: Amphipoda, Bivalvia, Cumacea, Decapoda, Echinodermata, Gastropoda, Hydrozoa, Nematoda (i.e. soft sediment and hard substrate benthos), and the pelagic prey items Fish and Zooplankton. Partially digested prey individuals that could not be allocated to one of these groups were not further considered. For mackerel, no prior information deriving from the SCA was used for the mixing model to meet the assumption that the mixtures of the consumer fall within the polygon of the prey items. Therefore, more prey items, such as suspended particulate organic matter from surface and bottom waters (SPOM-surface and SPOM-bottom, respectively), were considered as potential prey items for this species.
The prey items included in the mixing models of all the fish species were already a mixture of isotopically similar sources. This grouping was performed because the discriminatory power of the mixing models decreases with the number of food sources (Phillips et al., 2014). Furthermore, an a posteriori combination of prey items was performed when a strong negative correlation indicated that the prey items were isotopically similar and, thus, were undistinguishable (Parnell et al., 2013). Stable isotope signatures of these broad taxonomic categories, zooplankton, SPOMsurface and SPOM-bottom, sampled during the same sampling event were obtained from Mavraki et al. (2020a, b), while nematodes, Cumacea and fish (Engraulidae) stable isotope data were acquired from Franco et al. (2008), Kürten et al. (2013) and Jennings and Cogan (2015), respectively ( Table 1). The isotopic signatures of these prey items did not show a large variation (nematodes: d 13 C: -18.0 ± 1.0% and d 15 N: 10.0 ± 0.9%; Cumacea: d 13 C: -13.18 ± 0.2% and d 15 N: 9.68 ± 0.1% and fish: d 13 C: -18.8 ± 1.1% and d 15 N: 11.7 ± 0.3%) assuring that the results would not be affected by using bulk food resources. We further constrained the mixing model with the known dietary proportions for juvenile cod, pouting, horse mackerel and sculpin acquired by the SCA using the function simmr_elicit (Parnell & Inger, 2019).
Trophic enrichment factors (TEFs) were added to the mixing models. Mixing models are sensitive to assumptions regarding TEFs (Phillips et al., 2014) because of the multiplicity of factors that can affect them (Phillips et al., 2014). Hence, the use of the most appropriate TEF is crucial to produce the most accurate results. In order to accomplish this, we used the function simmr_mcm_tdf that can estimate the most appropriate TEF values for a known set of dietary proportions (Parnell & Inger, 2019). Hence, every species had different mean and standard deviation values of TEFs (Table 2). The TEFs defined by the model were within the reported range of -1 to 6% for nitrogen and -0.8 to 4% for carbon (Caut et al., 2009). The use of these TEFs assured that the fish' isotopic ratios fitted within the polygon of the prey items (Fig. S1), which is one of the prerequisite of isotopic mixing models (Phillips et al., 2014). Furthermore, the use of the functions simmr_elicit and simmr_mcm_tdf assured that the density distributions were not constrained (Fig. S2). The isotopic niche and the dietary proportion estimations were performed in R (R Development Core Team, 2018). The diet proportions are reported as the mean and standard deviation, but also as the median (50th percentile) and the 95% credible intervals (2.5th and 97.5th percentiles) to avoid misinterpretations of unique (mean) values (Phillips & Gregg, 2003;Phillips et al., 2014).

Diet analysis of the fish species
The SCA showed that the stomachs of the individuals of sculpin contained the lowest, and pouting the highest number of food items. Pouting, juvenile cod and horse mackerel stomachs mainly had the amphipod Jassa herdmani (Walker, 1893) (87.6 ± 10.2%,   , 1769) were observed in smaller quantities. Mackerel stomach content was dominated by zooplankton (83.9 ± 8.3%), while sculpin stomachs mainly contained both Pisidia longicornis (Linnaeus, 1767) (40.0 ± 11.3%) and fish (30.0 ± 10.0%- Table 3). After allocating the prey items to broader taxonomic groups (see ''Stable isotope analysis'' section), the stomach contents of sculpin showed the lowest variety (3 prey groups); highest prey group variety (9 groups) was observed for juvenile cod (Fig. 1). Occasionally, empty stomachs were found. One empty stomach was found for cod and mackerel, while the stomachs of three sculpins were empty as well. ANOSIM revealed significant differences in the stomach contents of the different fish species (Fig. 2-Global R = 0.244, P = 0.001). The pairwise tests indicated dietary overlaps between cod, horse mackerel and pouting. Moreover, dietary overlap was also observed between cod, horse mackerel and sculpin (pairwise R values \ 0.25) ( Table 4). Moderate dietary overlap was found between mackerel, horse mackerel and cod (0.25 \ R \ 0.5), while low overlap was observed between mackerel and sculpin and pouting (pairwise R [ 0.5).

Isotopic niches
Five individuals of cod, horse mackerel, sculpin and pouting and eleven mackerels were analysed for C and N stable isotopes ( Standard ellipse areas corrected for small sample sizes (SEA C ) (a proxy of the trophic niches) indicated differences between the five fish species (Fig. 3b, Table 6). The smallest SEA C (0.22% 2 ) was observed for cod, followed by pouting (0.88% 2 ). The largest SEA C (24.64% 2 ) was observed for mackerel. Generally, high overlaps were observed between the isotopic niches of the species with the same ecology (Table 6, Fig. 3b). The isotopic niche of the benthic species sculpin was the only one that did not overlap with any of the isotopic niches of the other species.

Discussion
Our study showed that two benthopelagic (cod and pouting) and one benthic (sculpin) fish species use artificial hard substrates, such as those in the OWFs, as feeding grounds for a prolonged period, since both the short-(stomach content analysis) and the long-   SPOM-surface and SPOM-bottom), which suggests that it might not exploit the artificial hard structures as feeding habitats.
The stomach content analysis derived from a limited number of individuals that have been collected during day time in a single season (summer), hence not necessarily fully reflecting the natural dietary variability. Increasing the number of stomach samples would result in a more reliable representation of the actual diet (Ferry & Cauilliet, 1996). However, this was not possible in this study due to the logistic efforts mentioned above. While the number of replicates was small, the stomach content results pointed in the same direction as previous dietary studies for these five species. Furthermore, the stomach contents did not differ remarkably between individuals of the same species, hence illustrating the reliability of our results. The stable isotope analysis was based on a small number of replicates, which, however, was always higher than the absolute minimum sample size (n [ 3, Jackson et al., 2011). The SIA results suggested a clear dietary separation of the fish species according to their ecology. The isotopic signatures of the benthic Table 6 Standard ellipse areas corrected for small size samples (SEA C -% 2 ) and overlaps (% 2 ) between the isotopic niches of the different fish species collected at the offshore wind turbine Species 1 Species 2 SEA C 1 (% 2 ) SEA C 2 (% 2 ) Overlap (% 2 ) (sculpin) and benthopelagic (pouting and juvenile cod) species indicated a diet based on artificial hard substrate prey items, while that of mackerel and horse mackerel showed pelagic feeding preferences. The d 13 C value of sculpin is characteristic of a benthic diet, which is explained by the exploitation of the crab Pisidia longicornis that consumes sediment organic carbon when it occurs at the scour protection layer (Mavraki et al., 2020a). The d 13 C signatures of pouting and cod reflect the signature of the amphipod Jassa herdmani that mainly feeds on zooplankton when it is found on the turbine foundations (Mavraki et al., 2020a). The depleted d 13 C values of the pelagic fish species indicate that they consume prey items that rely their diets on primary producers. Finally, the use of the informed priors in the mixing models did not affect the results of this study compared to conducting the models without the prior information. Hence our results were not biased towards more abundant but lower in biomass items in comparison with less numerous and larger prey items. Therefore, we are certain that our models highlighted the actual dietary proportions of the examined consumers.

Exploitation of the artificial habitats
The results of our study indicated that the benthopelagic species (pouting and juvenile cod) were mainly feeding on fouling organisms, and specifically on Jassa herdmani, while the benthic species (sculpin) was predominantly feeding on decapods, and especially on the porcelain crab Pisidia longicornis. Previous stomach content analyses of pouting and juvenile cod in the area have indicated similar results to ours, with J. herdmani being the most important food resource for these two benthopelagic species (Reubens et al., 2011(Reubens et al., , 2014a. Both J. herdmani and P. longicornis are exclusively associated with subtidal artificial hard substrates, such as wind turbine foundations and shipwrecks, where they occur in high densities (Leonhard & Pedersen, 2006;Zintzen, 2007;De Mesel et al., 2015). They are, however, rather rare Table 7 Proportion of prey categories in the diet of the five fish species sampled at close proximity from the offshore wind turbine as estimated by the combination of stomach content analysis and stable isotope mixing models Bold values indicate the highest percentages of the food items as they were identified bythe mixing models in natural hard substrate communities (Zintzen, 2007). Jassa herdmani and P. longicornis have sufficient lipid levels (2.83 and 3.19%WW, respectively-Heirman, 2011) reflecting their essential energy content, which substantiates their importance as food items for highly mobile species, such as cod and pouting. Both SCA and SIA indicated remarkable dietary overlaps between cod and pouting, while SCA suggested a possible trophic competition between cod, sculpin and horse mackerel. However, the level of trophic competition depends on the degree of the foodresource availability (López-Jamar et al., 1984). Since J. herdmani and P. longicornis occur in the OWF area in high densities (Zintzen, 2007;De Mesel et al., 2015), consumers could exploit the same resources with no or limited competition (López-Jamar et al., 1984). This in combination with the sufficient energy levels that have been reported for both J. herdmani and P. longicornis (De Troch et al., 2013) supports the idea that there is no trophic competition between the fish species.
In addition to their abundance, J. herdmani and P. longicornis are energy-and fatty-acid-rich food resources (De Troch et al., 2013). Jassa herdmani contains substantial concentrations of the fatty acid docosahexaenoic acid (DHA-De Troch et al., 2013), which is an essential fatty acid for the eye development of many fish species, such as cod (O'Brien-MacDonald et al., 2006). Pisidia longicornis, on the other hand, has higher lipid levels than J. herdmani (3.19 and 2.83%WW, respectively-Heirman, 2011) and contains considerable amounts of the essential fatty-acid eicosapentaenoic acid (EPA), which is an important fatty-acid for growth and fish development (Sargent et al., 1999). These two species contain higher lipid levels than the fouling species Monocorophium acherusicum and the J. herdmani mats (Heirman, 2011;De Troch et al., 2013). The energy content of J. herdmani and P. longicornis is also higher than that of food sources of pelagic origin, such as copepods (Michaud & Taggart, 2007), and of soft substrate origin, such as brown shrimp (Crangon crangon (Linnaeus, 1758)) larvae (Urzúa & Anger, 2013).
The pelagic horse mackerel seems to feed only occasionally on fouling organisms. This follows from the contradicting dietary results between the SCA (diet mainly based on fouling amphipods) and the SIA, as revealed by its isotopic niche (diet mainly composed of pelagic prey items-since the overlap of the isotopic niches of horse mackerel and mackerel was large). This was the only species for which SCA and SIA did not reveal similar results, confirming that both methods can yield complementary information on dietary compositions (Kadye & Booth, 2012;Hertz et al., 2017). Discrepancy between SIA and SCA can be a result of the differences in assimilation efficiencies among the prey items (Hertz et al., 2017). This discrepancy could also derive from the similarities in the stable isotope ratios of zooplankton and Jassa herdmani (Mavraki et al., 2020a) resulting in a less reliable resolution of the SIA for these two main food resources and leading to the a posteriori combination of these two prey items in the mixing model. A final explanation could also be found in the feeding behaviour of horse mackerels. Horse mackerels are opportunistic feeders that modify their diets according to the availability of the food items (Rumolo et al., 2017), and they can feed on both benthic and pelagic prey items (Van Ginderdeuren et al., 2014). Some pelagic species, such as horse mackerels, may take advantage of areas with high food availability to acquire and store energy (Jardas et al., 2004). Fat and energy content of horse mackerels has been reported to be lowest during and after the spawning period in late spring and summer, when considerable energy is spent on reproduction, while it is highest in autumn (Abaunza et al., 2003). In the North Sea, the energy content of horse mackerels increases rapidly in August and September (Abaunza et al., 2003), which coincides with our sampling period (see ''Sampling area and data collection'' section). Hence, we suggest that horse mackerels occasionally use the artificial habitats of OWFs as feeding grounds, where they can benefit from the locally increased availability of energy-rich food sources (J. herdmani) to support their highly mobile life style. It is possible, though, that this is a seasonally recurring pattern that needs further research.
Mackerel was the only species that was not predominantly feeding on fouling organisms, as both the short-and the long-term analyses indicated a diet based on pelagic prey items, which are indeed known to be its main food sources (Van Ginderdeuren et al., 2014;Bachiller et al., 2018). Hence, the occurrence of the mackerel in the vicinity of the OWF installations may either be associated with the presence of other pelagic species in the area (Reubens pers. com.), supporting the meeting point hypothesis that has been suggested for pelagic fish species (Dagorn & Fréon, 1999), or be linked to the enhanced protection against predation and the favourable current patterns that are provided by these habitats (Petersen & Malm, 2006).
All the fish individuals in this study were caught in close proximity to an offshore wind foundation. However, in the Belgian part of the North Sea, different types of artificial hard substrates exist, providing high abundance of fouling organisms (Zintzen, 2007). More than 300 shipwrecks, mainly from the World Wars I and II (Zintzen, 2007), have been mapped in the area (Fig. S3), while it is possible that many more unmapped wrecks are present (Reubens, 2013). Some of the wrecks are situated within OWF areas (Fig. S3). Even though the catch rates of pouting and juvenile cod are significantly higher at OWFs compared to shipwrecks in the Belgian part of the North Sea (Reubens et al., 2013a), we cannot exclude that the fish were also foraging for food resources at shipwrecks instead of offshore wind foundations. The only exception could be cod, since research on the juvenile individuals' movement has shown strong residency towards the OWFs (Reubens et al., 2013b). Cod juveniles exhibit crepuscular movements towards the OWFs, increasing their foraging success and reducing the predation pressure (Reubens et al., 2014a) and, at the same time, they settle in the OWF area (Reubens et al., 2014b). The combination of these studies with our results could suggest that cods remain in the OWF area to forage for a prolonged period of time.
Overall, the results of this study indicate that there are two main pathways supporting the feeding ecology of the fish species associated with OWFs: the 'pelagic pathway' that was based on the consumption of prey items of pelagic origin, such as zooplankton and SPOM, and the 'benthic pathway' that was based on the local benthic production. Both these pathways could lead to local production, since the main benthic prey item, J. herdmani, is a suspension feeder that mainly consumes zooplankton (Mavraki et al., 2020a). Thus, organic matter of pelagic origin might significantly contribute as a food resource due to the predominance of suspension feeders on the artificial hard substrates. These findings are in line with earlier observations on the feeding ecology of fish species related to artificial habitats in the Mediterranean Sea (Cresson et al., 2014) and could provide an indication of the co-existence of attraction and production on the same artificial reef system (Cresson et al., 2019).

Attraction-production hypotheses
Offshore wind farms and other artificial reefs could potentially act as production sites for fish species (Leonhard & Pedersen, 2006) that feed on the abundant fouling fauna. The production hypothesis can be supported when the artificial habitat contributes to an increasing biomass and/or abundance of fish (Svane & Petersen, 2001) and the fishing pressure is limited (Samples & Sproul, 1985). This is due to the increased feeding efficiency, survival rates of larvae and protection against predation at the artificial habitats, which encourage fish to settle and concentrate here (Brickhill et al., 2005;Petersen & Malm, 2006). The feeding ecology of the benthic (sculpin) and benthopelagic (pouting and juvenile cod) species as revealed in this study, in combination with the prohibition of fisheries within the Belgian OWF area since 2009 (De Backer et al., 2019), could be an indication that these species could potentially increase their local production. This is due to their long-term feeding ecology, which suggests that they remain close to artificial habitats for a long period of time to feed on prey items exclusively occurring on manmade structures. Similar studies in the Gulf of Mexico have suggested that artificial reefs can increase fish production when fish exploit primary and secondary producers associated with the artificial habitat (Daigle et al., 2013;Dance et al., 2018). However, it should be noted that the artificial reefs might promote fish production but this does not necessarily lead to enhanced fish stocks. Increased fish production could take place not by adding new fish in the system but by protecting and enhancing growth of the already present individuals, leading to enhanced fish biomass (Powers et al., 2003). Support for the production hypothesis has already been indicated for pouting and cod in the Belgian OWF area (De Troch et al., 2013;Reubens et al., 2013c). Specifically, research has shown that juvenile cod individuals indicate high residency and site fidelity in offshore wind foundations (Reubens et al., 2013b) and that they settle in the area, suggesting increased production in terms of additional recruitment of individuals (Reubens et al., 2014b). Furthermore, larger lengths and better fitness proxies (i.e. growth, condition index and diet) have been observed for pouting individuals in close proximity to OWFs compared to individuals occurring in sandy areas (Reubens et al., 2013c). Cod and pouting exploit lipidand fatty-acid-rich food sources (J. herdmani and P. longicornis) that can provide them with sufficient energy to sustain their metabolism and growth, and could also contribute to pouting's reproduction (De Troch et al., 2013). The combination of our findings with these previous studies could corroborate the production hypothesis.
The findings of our study imply that sculpins could also increase their local production within the artificial reef area of OWFs. Indeed, artificial habitats have been reported to have a positive effect on the production of benthic carnivores, due to the increased food availability and the low fishing mortality (Cresson et al., 2019). It is, thus, possible that artificial reefs, such as OWFs, could enhance the production of benthic species.
The results of our study corroborate the hypothesis that both attraction and production can co-exist on the same artificial reef, depending on the habits of every species (Cresson et al., 2019). It is suggested that benthic and benthopelagic species could have an increased production within OWFs, but this cannot be demonstrated for the pelagic species. The opposite patterns observed between pelagic and benthic/benthopelagic species are typical in natural reefs but were seldomly noticed for artificial reefs (Morris et al., 2018). Recently, similar patterns at the functional group-level have been observed for fish species at artificial reefs in the bay of Marseilles (Cresson et al., 2019), where pelagic species do not feed on fauna associated with artificial reefs. In our study, mackerel seems not to be predominantly feeding on artificial reef fauna, while horse mackerel exploits them opportunistically. Visual observations have shown that horse mackerels approach the offshore wind foundations to predate, taking advantage of the shade that the foundations offer (Reubens pers. comm.). Hence, the pelagic species seem to be attracted towards the offshore wind turbines, but further research is needed concerning their production in the area. Previous studies have indicated that the effect of artificial reefs on the production of pelagic zooplanktivores (i.e. mackerel and horse mackerel) is negligible (Powers et al., 2003), since most of these reefs are too small to significantly affect the pelagic communities (Cresson et al., 2019). To conclude, further investigation is needed about whether and how the pelagic species exploit the artificial reefs.

Conclusions
Our results indicated that both benthopelagic (pouting and juvenile cod) and benthic species (sculpin) that occur abundantly at offshore wind farms (OWFs) use the artificial reefs as feeding grounds for a prolonged period, since stomach content analysis (SCA) and stable isotope analysis (SIA) reveal a diet that is based on fauna typically found on artificial hard substrates. These species mainly consume fouling amphipods and decapods, which are abundant on these substrates. High trophic niche overlap was observed between the benthopelagic species, while the trophic niche of sculpin was not overlapping with any of the trophic niches of the other species. The pelagic species in our study either consume fouling organisms occasionally but in large proportions, i.e. horse mackerel, or they do not exploit them predominantly, i.e. mackerel. In the case of horse mackerel, artificial reefs act as feeding oases, providing high availability of prey items. Mackerel is largely feeding on pelagic prey items even when it is found in close proximity to OWFs and hence, the reason of its attraction to the artificial hard substrates requires further investigation. These findings corroborate to the hypothesis that OWFs could potentially increase the local production of fish species and especially of benthopelagic and benthic species, while further research is necessary for the pelagic species.
original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.