Introduction

Transitional waters (estuaries and coastal lagoons) are highly productive ecosystems that support a variety of fishery activities, most of which are practised at small-scale (Granzotto et al. 2001; Crespi 2002; Gomez et al. 2006; Pérez-Ruzafa and Marcos 2012). In Mediterranean coastal lagoons, fixed gears and barriers at sea inlets are some of the most widely used techniques (Cataudella et al. 2015). They exploit faunal movements either within the lagoon basin or between the lagoon and the adjacent marine and freshwater ecosystems, and target a variety of species including fish, decapods and cephalopods (Cataudella and Ferlin 1984; Quignard 1984; Ardizzone et al. 1988; Chauvet 1988; Pérez-Ruzafa and Marcos 2012; Pranovi et al. 2013; Zucchetta et al. 2016). Nekton assemblages found in lagoons represented for millennia one of the major sources of livelihood and income for local populations, making fishery a primary ecosystem service provided by the high productivity and resilience of transitional waters (Pranovi et al. 2013; Cataudella et al. 2015; Rova et al. 2015). Fishery productivity in these ecosystems has been maintained by means of a variety of hydraulic and morphological interventions, which evolved through time along with the development of fishing techniques themselves. Lagoon fisheries therefore not only play an important role in food provisioning for local communities, but also indirectly contribute to the conservation of physical and biological diversity (Pérez-Ruzafa and Marcos 2012; Cataudella et al. 2015). In the Venice Lagoon, the historical preservation of shallow water morphologies and management of hydraulic patterns fostered the development of a particularly complex and diversified fishery, which is still one of the most important in the Mediterranean basin in terms of both fishing effort and yields (Provincia di Venezia 2000; Provincia di Venezia 2009; Provincia di Venezia 2015). The regulating and supporting role of morphological heterogeneity in the Venice Lagoon has also motivated the inclusion of habitat restoration practices among the management options for the lagoon basin in more recent years (Curiel et al. 2005; Facca et al. 2014; Feola et al. 2018; Scapin et al. 2019b, c), integrating the traditional approaches based on water quality and hydraulic management (Pérez-Ruzafa et al. 2011; Cataudella et al. 2015).

Spatial and temporal dynamics in physico-chemical parameters play a major role in driving the distribution of nekton in transitional water ecosystems (Elliott and Hemingway 2002; Harrison and Whitfield 2006). Species found within estuaries and coastal lagoons are usually well adapted to such variability; however, temperature, salinity, dissolved oxygen and turbidity levels (among others) can determine physiological and behavioural responses in nekton species, including those that are targeted by local fisheries (Blaber and Blaber 1980; Whitfield et al. 1981; Marshall and Elliott 1998; Maes et al. 2004). For instance, water temperature is known to affect both reproductive investment and fishery catches of grass goby (Zosterisessor ophiocephalus) (Zucchetta et al. 2012; Pranovi et al. 2013), one of the primary fishery resources in the Venice Lagoon (Provincia di Venezia 2009; Provincia di Venezia 2015).

Transitional water ecosystems typically feature a high morphological heterogeneity, resulting in a variety of habitats that are used by nekton species as reproductive, feeding and nursery grounds (Pihl et al. 2002; Pérez-Ruzafa et al. 2007; Whitfield 2017). Morphology significantly influences not only nekton assemblages, but also fishery yields in coastal lagoons. In particular, higher shoreline complexity and amount of shallow-water areas can improve yields by intensifying physico-chemical gradients and nutrient inputs (Pérez-Ruzafa et al. 2007; Pérez-Ruzafa and Marcos 2012). In recent years, seascape ecology emerged as a new approach to effectively investigate the role of habitat mosaic on faunal distribution in coastal ecosystems (Boström et al. 2011). The habitat mosaic approach is successfully contributing in advancing transitional water ecology, by demonstrating the influence of habitat composition (i.e. typology, extent and diversity) and configuration (i.e. spatial arrangement and geometrical complexity) on nekton distribution in estuaries and coastal lagoons (Irlandi and Crawford 1997; Pittman et al. 2007; Staveley et al. 2017; Scapin et al. 2018). Seascape studies have major implications for ecosystem management and biodiversity conservation, tackling ecological issues from a large-scale perspective and taking into account spatial heterogeneity (Mumby 2006; Engelhard et al. 2016; Betzabeth and de los Ángeles 2017). Similarly, understanding how seascape structure determines nekton distribution in transitional water ecosystems may support fishery resource management, for instance by providing a context under which the different options for preserving stocks can be compared (Mizerek et al. 2011), or by identifying areas with the highest potential for recovery from fishing impacts (Stamoulis et al. 2018).

Despite the increasing amount of literature dealing with habitat mosaic influence on commercially important nekton, only very few studies have directly investigated whether seascape attributes influence fishery yields in estuarine and coastal ecosystems (Meynecke et al. 2008; De la Torre-Castro et al. 2014). While previous studies emphasised the role of broad morphological categories in influencing fishery yields in coastal lagoons, such as lagoon size, degree of complexity, shallowness (Pérez-Ruzafa et al. 2007), there is a general lack of knowledge about the contribution of habitat mosaic structure in supporting local fisheries. Going beyond the traditional approach of fishery studies, which focused on water quality and productivity as the main drivers of fishery production in transitional waters, this paper takes into account an additional level of complexity and investigates the potential influence of seascape structure on fishery yields in the Venice Lagoon. Following the habitat mosaic approach and including two different spatial scales in the analysis, it disentangles the role of morphology, habitat distribution and habitat mosaic properties from that of water physico-chemical parameters. The outcomes of this study could help improve our understanding of how the conservation and management of coastal lagoon seascapes may benefit local fisheries in the next decades.

Methods

Venice Lagoon and Study Areas

The Venice Lagoon is the largest coastal lagoon in the Mediterranean Sea (approximately 550 km2). It is connected to the northern Adriatic Sea by three inlets, and experiences microtidal conditions with a tidal range of ± 0.50 m during spring tides (Umgiesser et al. 2004). Most of the lagoon basin is composed of shallow water areas (average lagoon depth 1.2 m) that are intersected by a network of channels that could exceed 15 m of depth (Solidoro et al. 2004; Molinaroli et al. 2007). Environmental conditions in the lagoon are characterised by strong spatial and temporal gradients that are influenced by both natural processes and anthropogenic pressures (Solidoro et al. 2010). The high degree of environmental heterogeneity and the geomorphological variability contribute to creating a highly diverse seascape of islands, salt marshes, tidal creeks, mud- and sandflats with or without transient macroalgal beds, seagrass meadows, channels and man-made structures (Franzoi et al. 2010; Solidoro et al. 2010; Scapin et al. 2018).

For the purposes of this study, following the approach of Zucchetta et al. (2016) six broad macroareas were identified in the lagoon based on degree of confinement (Cucco et al. 2009), salinity, sediment characteristics and overall seascape features: Chioggia (CH), Ca’ Zane (CZ), Lago dei Teneri (LT), Ponte della Libertà (PL), Poveglia (PO) and San Giuliano (SG) (Fig. 1a).

Fig. 1
figure 1

Six macroareas in the Venice Lagoon (a) and an example of subdivision in fishing areas and spatial scales considered in the study (b). The location of fishing areas (numbered) can change over time due to changes in the position of fishing gears within macroareas (see text for details). A stylisation of the fishing gear considered in the study is shown. Habitats quantified at the broad scale are also indicated

Fishery Monitoring

Fyke nets are the most important fishing gear currently used by artisanal fishery in the Venice Lagoon (Granzotto et al. 2001, 2004; Pranovi et al. 2013; Zucchetta et al. 2016). Fyke nets consist of funnel-shaped, unbaited traps arranged in series along barrier nets that are deployed on lagoon shallow bottoms (Provincia di Venezia 2009; Provincia di Venezia 2015). Following Zucchetta et al. (2016), artisanal fishery activities were monitored within the six broad areas previously described. Monitoring took place monthly during the main fishing periods, namely spring (April to June) and autumn (September to November) (Provincia di Venezia 2009). Data from 8 years (2010 to 2017) were collected. Within each broad area, fishermen often change fyke net locations depending on the season, weather and tide conditions and on the specific targets pursued (Provincia di Venezia 2015). As a result, location and extent of actual fishing areas could change among monthly surveys, while always remaining within the boundaries of each broad area as defined in the previous paragraph (Fig. 1b). Hence, also the associated morphological and seascape features could vary over time. This allowed the seascape dataset to better represent the morphological and habitat variability of the lagoon.

For each study area, an average number of 40.9 traps (standard deviation 18.2) were inspected during each sampling occasion. Trap content was inspected on board with the collaboration of professional fishermen, and samples were collected to confirm identification in the laboratory. Catches were identified at species level and weighed (± 1 g). A total of 2574 samples were then included in the analysis.

Nekton fauna (fish, swimming decapods and cephalopods) were considered in the analysis. Also, total catches of green crab Carcinus aestuarii were included in the analysis, due to their importance for artisanal fishery in the Venice Lagoon (Provincia di Venezia 2009; Provincia di Venezia 2015). Green crab was not included among target species since only specific life stages (i.e. moulting specimens and egg-carrying females) have commercial value and are retained by fishermen. In the present work, it was not possible to quantify proportions of commercial stages since selection of individuals is usually performed by fishermen during the days after the capture (Matozzo et al. 2013; Pranovi et al. 2013). Cumulative biomass per species was then standardised per trap and days of gear deployment, hence expressed as catch per unit effort (cpue). Only for fish species, abundance was also recorded, and the average specimen size was estimated as weight/abundance ratio per species and sample.

Definition and Characterisation of Spatial Scales

In order to investigate the influence of lagoon morphology and habitat mosaic on fishery yields, data were analysed at a finer fyke-net scale, and at a broader fishing-area scale (see also Fig. 1b). Each spatial scale was characterised with a set of environmental parameters.

Temperature (°C), salinity, dissolved oxygen (saturation %) and turbidity (FNU) were recorded in situ during each fyke net survey by means of a multi-parameter probe. In addition, water residence time (days) and total dissolved nitrogen, total dissolved phosphorus and chlorophyll-a concentrations in water (µg L−1) were estimated for each fyke net observation using spatial data collected in previous and in regular monitoring programmes (Cucco et al. 2009; Mag. Acque—SAMA, 2013; Regional Agency for Environmental Prevention and Protection (ARPAV), https://www.arpa.veneto.it/dati-ambientali/open-data/idrosfera/acque-di-transizione, URL visited on 27/11/2018). Such set of variables was then employed to characterise the finer scale (see also Appendix 2 Table 3).

Average bottom elevation (metres) and average sediment grain size (% of sand particles) were calculated for each fishing area and survey month using spatial data collected in previous studies (Mag. Acque, 2002; Mag. Acque—Selc, 2005; Mag. Acque—Thetis, 2005) (see also Appendix 2 Table 3). Six habitat typologies were taken into account: natural salt marsh, lagoon channels, bare intertidal flats, macroalgal beds, seagrass meadows and artificial marshes/flats. Mainland, islands and man-made fixed structures (e.g. bridges) were also included in the analysis and classified as “no habitat”, assuming them to represent physical obstacles for nekton movements. Following Scapin et al. (2018), the temporal variation of some habitats (namely natural salt marshes, creeks, intertidal/subtidal flats, channels and land) was considered negligible within the study period, hence the same land cover map (Mag. Acque, 2002) was used over the entire time span of the analysis. By contrast, temporal variation was taken into account for more dynamic habitats (namely seagrass meadows, macroalgal beds and artificial habitats, the latter created in different occasions during the study period). Cover maps of dynamic habitats were obtained from classification of Landsat-7, Landsat-8 (spatial resolution 30 m) and Sentinel-2 (spatial resolution 10 m) satellite imagery. Proportion (% of area occupied) of each habitat was calculated for each fishing area and survey month.

Due to the high collinearity among many commonly-used landscape metrics (McGarigal et al. 2002), habitat diversity and mosaic complexity were chosen as measures of mosaic composition and spatial configuration, respectively (Boström et al. 2011; Scapin et al. 2018). Habitat diversity (Shannon’s index calculated on habitat surface) and whole mosaic complexity (landscape shape index) were calculated following McGarigal et al. (2002) for each fishing area and survey month. Morphological properties (average bottom elevation and sediment grain size), habitat distribution and whole mosaic properties (habitat diversity and mosaic complexity) were employed to characterise fishing areas at the broader scale.

Data Analysis

Preliminary Analysis

A principal component analysis (PCA) was performed separately on standardised water parameters at the fine scale and on standardised morphological variables and habitat proportions at the broad scale. This served two purposes: (i) to summarise the main spatial and temporal patterns of environmental heterogeneity; (ii) to extract un-correlated information (PC axes) from sets of highly correlated variables and include them in a model framework as predictors of fishery yields (Zucchetta et al. 2016; Scapin et al. 2018). Relevant PC axes of both PCAs were selected for the analysis on the basis of proportion of variance explained and by visually inspecting PCA bi-plots, in order to assess which axes could be interpreted as meaningful physico-chemical and morphological patterns (Scapin et al. 2018).

Response Variables

All nekton species recorded during monitoring activities were categorised according to their importance for artisanal fishery in the Venice Lagoon. Species actively pursued by the fishery were classified as target species; non-target species with a commercial value were classified as accessory species (named incidental in Pranovi et al. 2013); species with no commercial value were classified as discard (Pranovi et al. 2013). As previously described, green crab was not included in fishery categories, to separate the contribution of this benthic species from that of the nekton compartment. See Appendix 1 Table 2 for a complete list of species caught and their allocation in fishery categories.

The influence of lagoon morphology, habitat distribution and whole mosaic properties was investigated on cumulative cpue of each one of the three fishery categories (target, accessory and discard), on green crab cpue, and on cpue and average size of each target species.

Model Formulation and Selection

Generalised linear mixed models (GLMMs; Zuur et al. 2009) were fitted separately for each response variable (cpue of fishery categories, target species and green crab and average size of target species). After inspecting the distributions and the mean/variance relationships of response variables (Zuur et al. 2007; Warton 2008; Warton et al. 2014), the tweedie family within the R package glmmTMB (Brooks et al. 2017) was selected as the most appropriate to model both cpues and average sizes.

To take into account the nestedness of sampling design (i.e. a set of fyke net traps sampled within each fishery area), and assuming that the intercept of the regression would change among both traps and fishery areas, a random intercept (with trap nested into fishery area) was specified for all the models. The interaction between survey season and year was included as fixed effect representing the temporal factor.

Following Zucchetta et al. (2016) and Scapin et al. (2018), a hierarchical approach to model selection was followed, so as to hypothesise different contributions of predictors measured at different spatial scales, ultimately disentangling the relative influence of morphology, habitat distribution and mosaic properties quantified at the broad (fishing area) scale from the effect of water parameters measured at the fine (fyke net) scale. Four model categories were included in the analysis, investigating the following hypotheses (Table 1): the response variable is not affected by any of the predictors (null model, category m0); the response variable is affected by temporal factor only (category m1); the response variable responds to temporal factor and water parameters at the fine scale (the latter included progressively as each relevant PC axis; category m2); and the response variable is affected by temporal factor, water parameters at the fine scale and by predictors measured at the broad scale (category m3). To further discern between morphology/habitat distribution and whole habitat mosaic properties, two formulations were introduced in category m3: a first formulation hypothesised that the response variable is affected by temporal factor, water parameters at the fine scale and by morphology and habitat distribution (the latter included as each relevant PC axis); a second formulation hypothesised that the response variable is affected by all the previous predictors, as well as by habitat diversity and mosaic complexity.

Table 1 Structure of GLMMs employed in this study. The hypothesis and the spatial scale explored by each model formulation are also described. c, random intercept, equals 1|Area/Fyke net; random structure is the same for all the model formulations

Best candidate models were selected by means of the Akaike information criterion corrected for small samples (AICc; Burnham and Anderson 2002), choosing the model with the lower AICc value (Zuur et al. 2009). Models with a difference in AICc lower than 10 were considered undistinguishable (Burnham and Anderson 2002). Each hypothesis formulated was then explored, by verifying if progressive addition of predictors would improve the overall fit of the model (Scapin et al. 2018).

All data analyses were performed in the R environment (R Core Team 2014). Plots were produced using the ggplot2 package (Wickham 2016).

Results

Environmental Patterns

The first three axes of the PCA performed on fine-scale variables were selected as relevant axes summarising physico-chemical parameters, since they explained 57.4% of total variance and highlighted three major environmental gradients (Fig. 2). The first axis identifies a confinement gradient, being positively correlated to salinity and negatively to water residence time and nitrogen concentration. The second axis highlights the seasonal variation in water characteristics, with temperature and chlorophyll-a concentration being negatively correlated to this axis, opposed to dissolved oxygen saturation. Finally, the third axis is associated to a water transparency gradient, with turbidity being negatively correlated to this axis. Despite identifying major environmental gradients, the PCA does not allow for a clear separation of observations made within different fishing areas (Fig. 2).

Fig. 2
figure 2

Biplots of principal component analysis performed on water parameters measured at the fine scale. Variables are abbreviated as follows: TDN, total dissolved nitrogen; t-RES, water residence time; TDP, total dissolved phosphorus; CHL-A, chlorophyll-a; TEMP, temperature; TURB, turbidity; SAL, salinity; DO, dissolved oxygen

The first two axes of the PCA performed on broad-scale variables were selected as relevant axes summarising morphology and habitat attributes, since they explained 61% of total variance and allowed for a clear separation among fishing areas (Fig. 3). The first axis is associated to salt marsh and bare intertidal flats, being also positively correlated with bottom elevation and negatively with proportion of macroalgal beds. As a result, it separates fishing areas characterised by higher proportions of salt marsh, intertidal flats and overall shallower bottoms such as CH and LT from areas dominated by macroalgal beds and subtidal bottoms, such as PL. The second axis highlights a major pattern of seagrass vegetation, with proportion of seagrass meadows and lagoon channels and sand content in sediments being positively correlated with this axis. It then separates fishing areas with sandier bottoms and higher proportion of seagrass meadows and channels, such as PO, CH and PL (only some observations) from areas characterised by lower seagrass coverage and higher proportion of macroalgae and fixed structures, such as PL, CZ (most observations) and SG.

Fig. 3
figure 3

Biplot of principal component analysis performed on morphological properties and habitat proportion quantified at the broad scale. See text for variable definitions

Interpreting the PCAs allowed the use of each PC axis as predictor in the GLMM analysis. Model category m2, investigating the relative influence of fine-scale water parameters, resulted in three alternative formulations, including (in addition to temporal factor): (i) the confinement gradient (first axis of PCA calculated on fine-scale variables); (ii) the confinement gradient and the seasonal variation in water parameters (first and second axes); (iii) the confinement gradient, the seasonal variation in water parameters and the water transparency gradient (first, second and third axes). Conversely, the first and second axes extracted from the PCA calculated on morphology and habitat distribution were both included in the two formulations belonging to model category m3.

Fishery Catch Composition

Total fishery catches ranged between 704 ± 31 g trap−1 day−1 (average ± SE) in CH and 3222 ± 245 g trap−1 day−1 in SG. Green crab accounted for a relevant proportion of total catches in all fishing areas, ranging on average between 35% in PO and 90% in SG. Target species accounted for the great majority of nekton catches, while the contribution of both accessory and discard species was the smallest in all the areas. Target nekton catches varied on average between 236 ± 18 g trap−1 day−1 in SG and 726 ± 102 g trap−1 day−1 in LT (Fig. 4a). On average, big-scale sand smelt Atherina boyeri and grass goby Zosterisessor ophiocephalus accounted for more than 65% of target catches in all areas. The proportion of big-scale sand smelt ranged on average between 32% in SG and 78% in LT. The contribution of grass goby Zosterisessor ophiocephalus was the largest in PL (39% of target catches) and the smallest in LT (4%) (Fig. 4b).

Fig. 4
figure 4

Total fishery catches (cpue, g trap−1 day−1) (a) and composition of target nekton catches (b) averaged by fishing area. Contributions of fishery categories (including green crab) and target species, respectively are shown

Some seasonal differences were also found (not illustrated). Three species accounted for more than 90% of fishery target cpue during spring: big-scale sand smelt Atherina boyeri (49.6% of target cpue), grass goby Zosterisessor ophiocephalus (30.7%) and common cuttlefish Sepia officinalis (10.1%). Conversely, during autumn an average of five species made up more than 90% of target cpue: big-scale sand smelt (66.4%), grass goby (11.4%), brown shrimp Crangon crangon (5.6%), lagoon shrimp Palaemon spp. (4.6%) and European eel Anguilla anguilla (4.3%). Target species also included European flounder Platichthys flesus, golden grey mullet Chelon auratus, thicklip grey mullet Chelon labrosus, flathead grey mullet Mugil cephalus and Mediterranean sand smelt Atherina hepsetus. The latter three were excluded from the following analyses due to of the low number of observations in the dataset.

Model Selection

Temporal factor (season and year) was a relevant predictor for all the three fishery categories considered, while water parameters at the fine scale were relevant for accessory and discard species. In particular, accessory species responded to confinement gradient, while discard was influenced by confinement, seasonal pattern in water parameters and transparency. Morphology and seascape properties at the broad scale were also relevant for all the categories. In particular, fishery targets were influenced by both morphology/habitat distribution and whole mosaic properties, while accessory and discard species responded only to mosaic. Green crab cpue responded to transparency at the fine scale and whole mosaic properties at the broad scale, in addition to temporal factor (Fig. 5).

Fig. 5
figure 5

Summary of model formulations selected for each response variable via AICc. Predictors included in selected models are highlighted in shades of grey representing temporal factor, water characteristics measured at the fine scale and seascape features quantified at the broad scale

In terms of cpue, all target species except for European eel responded to temporal factor. Furthermore, all target species except for European eel and golden grey mullet were influenced by water parameters at the fine scale. In particular, confinement gradient was a relevant predictor for brown shrimp, lagoon shrimp and grass goby. Seasonal patterns in water parameters were relevant for big-scale sand smelt, brown shrimp, common cuttlefish and grass goby, while transparency was relevant for brown shrimp and European flounder. At the broad scale, morphology/habitat distribution had an influence on big-scale sand smelt, lagoon shrimp and grass goby, while whole mosaic properties were relevant for all the considered species except for European eel, lagoon shrimp and common cuttlefish (Fig. 5).

Regarding average sizes, all species except for European eel and grass goby responded to temporal factor. Confinement gradient was a relevant predictor only for the average size of European flounder, seasonal patterns in water parameters were relevant only for the average size of common cuttlefish, while transparency influenced only the average size of European flounder. While morphology/habitat distribution did not influence the average size of any of the target species, whole mosaic properties influenced the average size of big-scale sand smelt and European flounder (Fig. 5).

Influence of Fine- and Broad-scale Environmental Parameters on Fishery Catches

Biomass of target species was positively influenced by seagrass distribution and habitat diversity, while being negatively influenced by the proportion of salt marsh and intertidal flats and by mosaic complexity. The biomass of accessory species responded positively to higher salinities and shorter water residence times (i.e. larger cpue values in less confined areas) and to habitat diversity, while being negatively affected by mosaic complexity. Discard biomass was positively affected by higher salinities and shorter residence times, being also negatively associated to seasonal water parameters (i.e. larger cpue during warmer periods) and to the transparency gradient (i.e. larger cpue at higher turbidities). Discard biomass was also positively affected by habitat diversity, and negatively associated to mosaic complexity (Fig. 6).

Fig. 6
figure 6

Standardised coefficients and standard errors of fine- and broad-scale environmental predictors for cpue of fishery categories and cpue and average size of target species, as estimated by selected model formulations. Fine-scale predictors include the first three axes of PCA performed on water parameters; broad-scale predictors include the first two axes of PCA performed on morphology and habitat distribution, and whole mosaic properties. *Green crab Carcinus aestuarii was not classified as target species, but was included in the plot for illustration purposes

The inverse confinement gradient had a negative influence on many target species, with biomass of brown shrimp, lagoon shrimp and grass goby and average size of European flounder showing larger values at lower salinity and longer water residence times. In terms of seasonal water parameters, the biomass of big-scale sand smelt showed larger values during warmer periods, while the biomass of brown shrimp, common cuttlefish and grass goby, as well as the average size of common cuttlefish showed larger values during colder periods. Transparency gradient had a negative influence on the biomass of brown shrimp and European flounder, and on the average size of European flounder (larger values in more turbid conditions), while showing a positive influence on the biomass of green crab (larger values in clearer conditions).

At the broad scale, morphology and habitat distribution had an effect on three lagoon resident species: proportion of salt marsh and bare intertidal flats in the mosaic had a positive effect on the biomass of big-scale sand smelt, while it showed a negative effect on the biomass of lagoon shrimp and grass goby; proportion of seagrass vegetation positively influenced the biomass of big-scale sand smelt, lagoon shrimp and grass goby.

Conversely, the effect of whole mosaic properties was evident on a larger group of species, comprising big-scale sand smelt, golden grey mullet, brown shrimp, European flounder, grass goby and green crab. Habitat diversity had an overall positive effect, with only the biomass of brown shrimp showing a negative relationship with this parameter. In turn, mosaic complexity showed an overall negative influence, with only the biomass of brown shrimp being positively related to it (Fig. 6).

Discussion

The hydrological, morphological and ecological complexity characterising coastal lagoons is regarded as one of the key determinants of their productivity, supporting substantial fishery and aquaculture industries (Kapetsky and Lasserre 1984; Quignard 1984; Ardizzone et al. 1988; Chauvet 1988; Pérez-Ruzafa and Marcos 2012). The presence and intensity of sea-mainland gradients within lagoon basins and the extent of shallow areas, in particular, are strongly linked to fishing yields. In Mediterranean coastal lagoons, the latter are correlated to the amount of water exchanged between lagoons and the sea, as well as to the intensity of freshwater discharges and to temperature ranges (Pérez-Ruzafa et al. 2007; Cataudella et al. 2015). To a significant extent, such parameters themselves are dependent on morphological traits of lagoon basins (Pérez-Ruzafa et al. 2007; Pérez-Ruzafa and Marcos 2012; Umgiesser et al. 2014). In this light, the present paper takes a step further in investigating the role of lagoon morphology as a determinant of fishery production. Morphological traits measured at the broad scale (4 to 30 km2), such as local bottom elevation and sediment grain size, and habitat mosaic properties, such as habitat proportion, diversity and shape complexity indeed influenced both catches and sizes of commercial nekton, in addition to physico-chemical properties measured at the fine scale.

Overall downward trends in fishery yields have been recorded in Mediterranean coastal lagoons over the last decades (Béné 2011; Cataudella et al. 2015). A comparison of fishery yields reported in previous studies (Provincia di Venezia 2000; Granzotto et al. 2001; Libralato et al. 2004; Provincia di Venezia 2009; Provincia di Venezia 2015) with data presented in this paper attests a decreasing trend in catches of various commercial species also in the Venice Lagoon, at least since the 1980s. These include big-scale sand smelt, grass goby, brown shrimp, lagoon shrimp, European eel and European flounder. Variations in environmental conditions due to multiple anthropogenic drivers played a role, including the long-term increase in temperature (Pranovi et al. 2013) and decrease in nutrient inputs (Pastres et al. 2004; Solidoro et al. 2010) recorded in the last decades.

However, also the general degradation of hydrological and morphological features undergone by coastal lagoons must be included among the factors that contributed to the decline of lagoon fisheries. Previous works have reported how the presence of structured habitats and their heterogeneity are crucial in maintaining fishery productivity, and how habitat degradation and loss are among the primary anthropogenic pressures on commercial nekton (Cataudella et al. 2015; Zucchetta et al. 2016; Cavraro et al. 2019). Breber et al. (2008), for instance, demonstrated how past reclamation of shallow water habitats in two European lagoons contributed to the observed decline in local fishery productivity. Similarly, loss in habitat connectivity impaired migration of two commercially relevant mugilid species, subsequently affecting fishery production in coastal lagoons of Tunisia (Romdhane et al. 2019). In the Venice Lagoon, erosion of lagoon seascape structures due to channel dredging, increased navigation and mechanical exploitation of Manila clam Ruditapes philippinarum that occurred in the last decades determined an overall increase in lagoon average depth and the loss of morphological heterogeneity (Sarretta et al. 2010; Solidoro et al. 2010). Some critical habitats for nekton species including salt marsh, intertidal and subtidal flats and seagrass meadows were particularly impacted (Franco et al. 2009; Thetis S.p.A. (Ed.) 2010; Zucchetta et al. 2016; Cavraro et al. 2017; Scapin et al. 2018). The intensive harvesting of manila clam had an additional effect on artisanal fishery. Not only did it affect sediment stability and dynamics, but it also determined a relevant shift in fishing targets and practices among fishermen in the Venice Lagoon. Such shift towards manila clam as a very abundant, readily available and more marketable product contributed to the observed decrease of artisanal fishing effort and to the subsequent decrease of fishery production during the late 1980s and 1990s (Provincia di Venezia 2000; Libralato et al. 2004; Provincia di Venezia 2009).

In this paper, both accessory and discard species among fishery categories appeared to be more influenced by water characteristics at the fine scale. Conversely, target species showed the more evident association with morphological and habitat mosaic properties at the broad scale. Target species that exhibited a clear link to morphological properties (bottom elevation and sediment grain size) and specific habitats (proportion of natural and artificial salt marsh and intertidal flats, seagrass meadows, macroalgal beds, lagoon channels) were essentially lagoon residents (Franco et al. 2008; Potter et al. 2013; Scapin et al. 2019a) (see Appendix 1 Table 2). In the Venice Lagoon as well as in other estuarine ecosystems, these species spend most of their life cycle in transitional waters, often showing morphological, behavioural and reproductive adaptations to particular habitats (Malavasi et al. 2005; Franzoi et al. 2010; Scapin et al. 2018; Kara and Quignard 2019).

Artisanal fishery in the Venice Lagoon shares some common features with other small-scale fisheries in Mediterranean coastal lagoons, including the widespread use of passive gears employed by local fishermen by means of small (< 12-m-long) boats (Pranovi et al. 2013). Nevertheless, it retains some relevant peculiarities that differentiate it from the others. Many lagoon fisheries rely on fixed structures, usually located at the sea inlet (e.g. lavorieri, bordigues and encañizadas in Italy, France and Spain respectively), mainly targeting highly motile, migrant nekton (mainly Mugilidae, Sparidae, Anguillidae and Moronidae) (Cataudella and Ferlin 1984; Kapetsky and Lasserre 1984; Ardizzone et al. 1988; Pérez-Ruzafa and Marcos 2012). Conversely, fishery performed with fyke nets in the Venice Lagoon relies to a larger extent on smaller-sized nektonic resident species such as big-scale sand smelt and grass goby as well as on green crab (Pranovi et al. 2013; Provincia di Venezia 2015). Such resident species are strongly linked to the morphological and habitat features that characterise the Venice Lagoon. Grass goby is usually associated to seagrass vegetation, whose root system help in stabilising burrow nests dug into soft sediments for wintering, summering and reproduction (Gandolfi et al. 1991; Malavasi et al. 2005; Franco et al. 2012). Big-scale sand smelt, despite being ubiquitous in the Venice Lagoon (Franco et al. 2006; Franzoi et al. 2010), select shallow-water structured habitats including macroalgal and seagrass beds as spawning grounds, to which they exhibit a certain degree of site fidelity (Henderson et al. 1988; Koutrakis et al. 2004; Kara and Quignard 2019). Fishery in the Venice Lagoon (and northern Adriatic lagoons nearby, e.g. Grado-Marano Lagoon) therefore relies to a larger extent, compared to those practised in other Mediterranean lagoons, on the presence of vast portions of fine-sediment, intertidal and subtidal shallow areas representing the main reproductive, feeding and shelter habitats for the major commercial targets.

This work highlighted that whole mosaic properties (habitat diversity and mosaic shape complexity) influenced to a higher extent marine migrant target species compared to resident ones. Most of these species share marine reproduction and spawning strategies, and enter the Venice Lagoon and other transitional water ecosystems either to exploit nursery habitats as juveniles (marine estuarine-dependent), for trophic purposes (marine estuarine-opportunists) or during migrations between marine and freshwater habitats (catadromous) (Franzoi et al. 2010; Scapin et al. 2019a). Compared to resident species, marine migrants are highly motile, and many of them exhibit substantial changes in habitat preference during their life cycle, while being subjected to ontogenic shifts in morphology and diet preferences (Pita et al. 2002; Tancioni et al. 2003; Georgalas et al. 2007; Andolina et al. 2020). Moreover, most species perform movements among habitats on a daily basis by following tidal cycles, seeking shelter and forage (Irlandi and Crawford 1997; Hampel et al. 2003; Rountree and Able 2007). This paper shows an overall positive response of target species to habitat diversity, suggesting that these species would benefit from the presence of more heterogeneous lagoon seascapes. For instance golden grey mullet C. auratus and European flounder P. flesus, which are traditionally targeted by transitional water fisheries in the Mediterranean, exploit the estuarine seascape during juvenile and sub-adult stages, performing movements among habitats while undergoing ontogenic diet shifts (Franco et al. 2010; Zucchetta et al. 2010; Amorim et al. 2018).

The results of the present paper therefore emphasise the critical importance of multi-functional shallow water habitats that are used as feeding, shelter and reproductive grounds by resident nekton in the Venice Lagoon and similar transitional water bodies. In addition, they highlight how a certain degree of habitat diversity is needed by sub-adult individuals of marine migrant species that significantly contribute to sustain lagoon fisheries.

Among commercial species considered in this analysis, common cuttlefish Sepia officinalis showed a clear response only to temporal variability (i.e. temporal factor and seasonal patterns in water parameters). This marine migrant species is regularly found in the Venice Lagoon (Provincia di Venezia 2009; Provincia di Venezia 2015). In contrast to most marine migrant species considered in this work, which perform migrations to exploit shallow-water nursery and feeding grounds, common cuttlefish penetrate transitional waters mainly for reproductive reasons. In the Venice Lagoon this species often gathers in deeper and more dynamic lagoon channels, possibly being attracted by a variety of structured habitats that are used for spawning (e.g. rocky substrata, subtidal seagrass meadows and artificial structures) (Guerra 2006; Guerra et al. 2016; Grati et al. 2018). As a result, cuttlefish is traditionally pursued with a wide range of fishing gears deployed in deeper waters, including trawl nets, gill nets and traps in addition to fyke nets (Provincia di Venezia 2015). This may partially explain the lack of response to both spatial gradients (i.e. confinement and transparency) and seascape parameters of fyke net catches in shallow waters.

Similarly, also green crab Carcinus aestuarii showed an overall weak response to shallow water seascape structure; this is probably due to the benthic habits and relatively low motility that differentiate this species from the nektonic taxa considered in this work. Despite that, in the Venice Lagoon green crab represent the most important taxon in fyke net catches in terms of total biomass, as well as one of the most valuable fishery resources during certain life stages (Provincia di Venezia 2000; Provincia di Venezia 2009; Provincia di Venezia 2015). In Northern Adriatic coastal lagoons, only moulting individuals and ripe females are indeed retained and marketed, after a temporary relaying phase following catch and selection (Matozzo et al. 2013; Pranovi et al. 2013).

Size of commercial species is of particular importance for fisheries, since larger specimens are usually more valuable. Body condition measures in fish are commonly used as indicators of ecosystem health in transitional waters, being employed to evaluate productivity and anthropogenic pressures among habitats (Amara et al. 2007; Piazza and La Peyre 2010; Cavraro et al. 2019). This work emphasises that fish individual weight in transitional waters can be influenced by morphology and seascape structure, these being potential drivers of fish individual health. This confirms findings of previous authors, who pointed out that human alterations to morphology in coastal lagoons impact negatively on individual fish condition (Cavraro et al. 2019). This suggests that preserving and restoring habitat mosaics in transitional water ecosystems could benefit local fisheries not only in terms of overall yields, but also in terms of sizes of individual catches.

Ecological restoration is considered a strategic approach to revert loss of biodiversity and degradation of ecosystem functionality (Suding 2011; Menz et al. 2013; Perring et al. 2015). It is also regarded as a way to contribute towards preserving and enhancing ecosystem services, including provision of natural resources such as food (Jones 2013; Perring et al. 2015). Restoration of estuaries and coastal lagoons is particularly critical, since they support strategic fishery and aquaculture industries, while being at the same time subjected to high levels of anthropogenic pressures (Béné 2011; Pérez-Ruzafa and Marcos 2012). Pérez-Ruzafa et al. (2007; 2012) previously pointed out the importance of features such as lagoon size, shallowness, coastline development and degree of sea-lagoon connectivity, emphasising that resource management in coastal lagoons aims in particular at maintaining or modifying such morphological traits in order to maximise productivity. In this light, this paper not only confirms the relevance of overall morphological features for fishery yields, but also introduces habitat composition and diversity as additional, primary factors affecting fishery productivity. By highlighting how morphological and seascape features can positively affect fishery yields, this paper suggests that conservation and restoration of morphological traits including habitats and seascape complexity and diversity can contribute to the management of lagoon fishery resources in the future.

Previous studies reported how habitat restoration would benefit fishery productivity in transitional waters. Restoration measures aiming at increasing marsh edge and overall complexity in salt marsh seascapes, for instance, could support the standing crop of commercial species associated with structured habitats (Rozas et al. 2005). Also restoring habitat connectivity would benefit fishery species, as pointed out by Mizerek et al (2011) in patchy seagrass seascapes. Moreover, large-scale restoration of lagoon and estuarine hydraulic regimes may substantially enhance both recruitment and commercial yields of migrant nekton species by increasing tidal exchanges and sea-lagoon connectivity (Mohapatra et al. 2007; Tsihrintzis et al. 2007; Reese et al. 2008).

This paper emphasises how management actions towards the preservation of fishery resources in coastal lagoons should be taken at the seascape level. The preservation of gradients in water characteristics, including salinity, turbidity and trophic status, should therefore be accompanied by the conservation and restoration of whole habitat mosaics. A seascape approach in coastal lagoon management could benefit local fisheries not only directly, by enhancing habitat availability and quality for target species as demonstrated here, but also indirectly, by tackling ecosystem-scale, long-term changes that could negatively affect fishery productivity (e.g. variations in temperature and trophic status). Coastal lagoons are inherently dynamic ecosystems, and their provisional services have always been strongly linked to human intervention (Béné 2011; Cataudella et al. 2015). In order to preserve lagoons’ capability to support populations of species of commercial interest, in the next years management should focus on restoring not only habitat functionality, but also connectivity among habitats and the overall mosaic heterogeneity.