Quantifying Abundance and Mapping Distribution of Loggerhead Turtles in the Mediterranean Sea Using Aerial Surveys: Implications for Conservation

: In the Mediterranean, incidental captures in fishing gear contribute to the high mortality of loggerhead turtles ( Caretta caretta ). Understanding the effects of bycatch is complex and requires robust knowledge of baseline population parameters such as abundance and density, as well as an understanding of animals’ distribution in relation to commercial fishing efforts. Based on data collected during multi-species line transect aerial surveys conducted between 2009 and 2017, we present density and abundance estimates, corrected for availability bias, for a large sector of the central Mediterranean, discuss temporal and spatial patterns and provide Potential Biological Removal (PBR) values for the monitored areas. Sightings data were also used to evaluate the spatial and temporal usage areas. Strong latitudinal and longitudinal gradients in density, abundance and area usage emerged from the analysis, with turtles occurring in higher numbers in the deeper pelagic waters of the Tunisian Plateau, the Ionian Sea and the Gulf of Taranto, irrespective of the season. PBR values derived from this study are likely unsustainable. This paper investigates the implications of commercial fisheries for Mediterranean loggerhead turtles across an area rarely included in sea turtle monitoring and has the potential to be relevant towards informed management and conservation of this species and highlights the necessity of international collaborative efforts in the region.


Introduction
Monitoring wildlife population parameters such as abundance and density over time and across space is pivotal to identify underlying drivers of changes and constitutes the backbone of evidence-based conservation (e.g., [1]).In this regard, assessing abundance and density for a given species is not only essential for managing populations and prioritizing conservation measures [2,3], but also a legal obligation under several regional was down-listed in 2015 from Endangered to Least concern, according to the International Union for the Conservation of Nature (IUCN) Red List criteria [8].Despite the downlisting, incidental captures in fishing gears still contribute to high mortality [9,10].In fact, the Mediterranean Sea shows the highest unsustainable bycatch rates of marine turtles worldwide [7,10,11].At the Basin level, Lucchetti and Sala [9] and Casale [12] report 132,000 sea turtles unintentionally caught during commercial fishing operations annually, with over 44,000 deaths per year.In Italian waters, Lucchetti et al., [13] estimated more than 52,000 bycatch events and up to 10,000 deaths in 2014 alone.More recently, in the north Adriatic Sea, the annual turtle bycatch due to bottom trawling has been estimated as 8600 individuals [14].In the region, despite growing efforts to monitor the conservation status of sea turtles, the difficulties in assessing the volume and the impacts of bycatch [15,16] prevent the effective evaluation of the extent and the real effects of this threat [17].
Assessing absolute, unbiased abundance and trends over time is a priority for sea turtle conservation [7]; however, evaluating these parameters and the effects of bycatch is challenging due to the inherent complexity of sea turtle life history [18].Traditionally, the monitoring of nesting beaches has been used to assess population parameters [19][20][21][22], but the drawbacks of this approach [18,[23][24][25] can result in drastic overestimates, making atsea sampling crucial to complement traditional approaches.
Dedicated aerial surveys are considered one of the most robust tools to gather information on sea turtle density, abundance and distribution [7,10,26], but while aerial line transect distance sampling surveys are regularly used worldwide to assess these population parameters for sea turtles [27][28][29][30][31][32], in the Mediterranean Sea, this approach has been used primarily to monitor Cetaceans [33][34][35][36][37][38].Only recently aerial surveys have been used to assess loggerhead turtles' abundance at a regional level [39], with relatively little effort allocated in the past [35,40,41].
Wide scale assessments of loggerhead turtle abundance are a pivotal first step to evaluate the sustainability of bycatch and eventually identify areas where the susceptibility of these species to commercial fisheries is high.In this context, this study presents original abundance/density estimates of loggerhead sea turtles for the central Mediterranean Sea and presents maps of inter-annual and seasonal use of space for those sectors monitored multiple times over the years and seasons.Finally, based on the obtained abundance estimates, this study presents Potential Biological Removal (PBR) values and discusses them considering known bycatch pressure in the area.

Study Area, Survey Design and Data Collection
The monitored area covers a large sector of the western and central Mediterranean Sea, extending from 5°42′ E to 19°42′ E and from 35°4′ N to 44°36′ N. It includes the seas of Corsica and Sardinia, the Ligurian Sea, the entire Tyrrhenian Sea, a large portion of the Sicily Strait and portions of the Ionian Sea and the Gulf of Taranto (Figure S1).Overall, the area extends for about 736,000 Km 2 .
The area has been divided into 6 strata and 12 blocks, within which survey transect lines were designed using the software Distance 5 and 6 [42].Except for the Ionian Sea and the Sicily Strait, all the other strata were subdivided into sub-blocks.In total, 298 transects were designed (Figure 1).Transects were spaced 15 km apart for the surveys carried out between 2010 and 2016 while for two surveys carried out in the Pelagos Sanctuary in winter and summer 2009, transects were spaced 10 km apart.A high-wing double-engine aircraft equipped with bubble windows was used as the research platform.Flight altitude and speed were kept constant at 750 feet (229 m) and 100 knots (185 km/h), respectively, for all surveyed areas and periods.Data collection followed standard line transect distance sampling procedures as described in [33,34,36,37,43].Table 1 summarizes the details for each survey.[35].The map inset with the extent of the monitored area is highlighted in red.
Table 1.Summary of the information for each surveyed area."Strata" = number of strata in which each surveyed block has been subdivided; "Extent" = extent of each survey block; "Transects" = total number of transects per block; "Spacing" = linear distance between parallel transects; "Length" = total length of planned survey tracks and "Coverage" = proportion of covered area for each stratum as resulting from survey design.

Abundance Estimates
Datasets arising from the different surveys were considered independently in the analysis and presented by season.The choice not to pool datasets was based on the fact that surveys took part in different years and months.Following DiMatteo et al. [39], two seasons were defined as winter months between November and April and summer months between May and October.For the 2013 and 2010 summer surveys, alongside pooled estimates for all the strata, separate estimates were calculated for the Pelagos Sanctuary and the northern Tyrrhenian Sea to allow for comparisons.
For density and abundance estimates of data collected between 2010 and 2016, Conventional Distance Sampling (CDS) was used and analyses were performed using the software Distance 7.3 [42].The methodological and analytical approaches for the surveys carried out in the Pelagos Sanctuary in winter and summer 2009 were presented in Lauriano et al. [35] and are not further detailed here.In CDS, no additional explanatory variables are considered along with perpendicular distance to the sightings in the estimation of the detection function [44].The choice of CDS models over other approaches has been based on the fact that all surveys were conducted in optimal conditions to maximize animals' visibility and, therefore, the influence of sea state and weather covariates were already minimized during data collection [30].Aircraft model, data collection protocols and observers were also kept constant during all surveys, minimizing potential biases introduced by changing sighting conditions.
One of the main assumptions of Distance Sampling is that all the animals on the transect line are recorded (G(0) = 1; [44]).This assumption is often violated when animals are at the surface, but the observer fails to detect them (perception bias) and/or when they are under water and cannot be detected (availability bias), although on the track line.These biases are not mutually exclusive and can negatively affect the estimates, unless accounted for [45,46].
Perception bias varies with the target species and the experience of the observer and can be minimized by using experienced observers and being accounted for in a doubleplatform framework [42,47].On the other hand, availability bias is linked to the activity periods and habitat types of the species of interest [48].Integrating tagging data (i.e., data on dive-surfacing patterns) with aerial survey datasets can help incorporate availability bias in the estimation of abundance [49].Due to logistic and economic limitations, a double platform configuration could not be set up during our surveys to quantify perception bias.However, we consider this bias to be low and constant because the same team of experienced observers collected the data during all surveys.Knowledge on sea turtles movement patterns and migration, although unevenly distributed across the Basin, is available [50][51][52][53][54][55][56][57][58].In this paper, we present both uncorrected estimates as well as estimates accounting for availability bias.The latter includes correction factors (Table 2) differing by habitat (neritic and pelagic) and season (winter and summer), presented in DiMatteo et al., [39].Specifically, for the summer survey and winter surveys, G(0) values of 0.48 and 0.43, and 0.43 and 0.27 were used for the pelagic and neritic regions, respectively.Table 2. Density and abundance estimates and abundance estimates corrected for availability bias for each survey.ER = encounter rate calculated as number of sightings per km of survey; D = uncorrected density; N = estimated uncorrected abundance; CV = confidence intervals for the estimated uncorrected abundance; 95% CI = lower and upper 95% confidence intervals for the estimated uncorrected abundance; G(0) = applied correction factor for availability bias; Nc = estimated abundance corrected for availability bias; models (best model based on AIC: HN/c = half normal with cosine adjustment; HR/c = hazard rate with cosine adjustment; U/c = uniform with cosine adjustment.Detection functions are provided in the Supplementary Material as Figure S3.For each survey and area, model selection was based on the Akaike's Information Criterion (AIC; Akaike 1974).

Turtles Usage Areas
Turtle usage areas (hereafter 'TUA') were identified via utilization distribution (UD) analyses representing the probability that an animal is found at a given point in space [59,60].We used the kernelUD function in the adehabitatHR package [61] for the software for statistical computing R (R version 4.3.0-www.cran.r-project.org;R Core Team 2021) to generate kernel UDs (KUD) for each survey on a 1000 m 2 grid using the reference smoothing parameter for the bandwidth.

Turtle Potential Biological Removal (PBR)
PBR is a reference point for human-induced mortality developed for marine mammals [62].PBR is based on the concept that the anthropogenic mortality of a population should not exceed 50% of its potential maximum productivity rate (Rmax), adjusted by a recovery factor (F) which can vary from 0.1 to 1 depending on the conservation status of the target population (0.1 for endangered, 0.5 for threatened, and 1 for good status; [63])).The PBR is well suited to a data-limited situation [64] like the one observed for Cheloniids in the Mediterranean Sea, and once the minimum population size (Nmin) is known, it can be calculated as: For the objectives of this study, PBR values were calculated considering both the corrected and uncorrected abundance estimates and, for the latter, its associated 95% CI for each of the surveys.Considering the species conservation status in the Mediterranean Sea and the fact that an Rmax of 0.064/year calculated at nesting sites is available for the species [10,65], PBR values were calculated using a recovery factor of 0.5 and 1.

Sightings and Effort
In total, eight surveys took place between 2009 and 2016.In 2009, the Pelagos Sanctuary was surveyed during the winter and the summer seasons.During 2010, the Ionian Sea and a large block including the Pelagos Sanctuary, the Seas of Corsica and Sardinia (hereafter referred to as the "Extended Pelagos Area") and the northern Tyrrhenian Sea were monitored in winter and summer, respectively.The southern Tyrrhenian Sea was surveyed twice during the winter season, in 2010-2011 and in 2014, while the Sicily Strait was monitored during winter 2016.
The results of the analysis of the 2009 winter and summer surveys in the Pelagos Sanctuary are presented in Lauriano et al. [35].Here, we refer and quote those results to allow for comparisons with surveys carried out between 2010 and 2016.
Overall, 5906 loggerhead turtle sightings (Figures 2 and 3) were recorded between 2009 and 2016, with the highest number of observations recorded in the Ionian Sea and the Sicily Strait in winter 2010 and 2016, respectively.The turtles' encounter rate, expressed as number of groups of turtles per kilometer of positive effort on transects, ranged between 0.01 and 0.11, with the highest values obtained from the Ionian Sea and the south Tyrrhenian Sea (Winter 2014) (Table 2).A latitudinal and longitudinal gradient is evident in the distribution of turtle sightings (Figure 4), with a generally higher occurrence of turtles at lower latitudes and higher longitudes.

Abundance Estimates
After excluding sightings with no declination angle and/or sightings with missing group size information, and after right truncation, 5177 sightings were considered in the CDS analysis.Details on model selection are summarized in Table 2.
For the summer surveys, the overall highest abundance was estimated from the northern Tyrrhenian Sea survey (n = 12,912; %CV = 13.In the Pelagos Sanctuary, the only area surveyed in both winter (2009) and summer (2009 and 2010), results show higher abundance for the summer surveys.Table 2 summarizes the results of the distance sampling analysis and presents both uncorrected and corrected estimates for each survey.

Turtles Usage Areas
Overall, high turtle usage areas were found in the central, southern and eastern sectors across the entire monitored area, primarily encompassing deep pelagic waters, and occurred at lower latitudes and higher longitudes, irrespective of the season or year.Winter maximum values of usage are double than the highest summer ones.In particular, the results of the UD analysis highlight two areas intensely used by sea turtles distributed across the central deepest portion of the Sicily Strait and the south-western sector of the Ionian Sea (Figures 2 and 3).Important areas of loggerhead presence also occur in the central, deeper portion of the southern Tyrrhenian Sea stratum, monitored in winter 2010-2011 and 2014, for both survey years, with a slight shift in the occurrence of turtles between the two surveys.In particular, for the 2014 survey, areas of intense usage emerge in the south-western portion of the area between south-eastern Sardinia and north-western Sicily and, to a lesser extent, in the deep northern portion of the stratum, as for the 2010-11 survey (Figure 3).
When considering the summer surveys, despite an overall lower presence of turtles in the monitored areas when compared with the winter ones, the areas with higher usage can be located in the deeper portion of the Ligurian Sea for the 2009 survey in the Pelagos Sanctuary, and in the deeper waters of the northern Tyrrhenian Sea, located in the centralsouthern portion of the stratum (Figure 2).For the Pelagos Sanctuary, surveyed in 2009 and in 2010, a stark difference emerged between the two years, with relatively fewer sightings recorded in 2010 than in 2009 (Figure 2) and a consequent absence of high usage area in this portion of the northwestern Mediterranean for the later survey.The northern Tyrrhenian Sea, also surveyed twice in 2010 and 2013, while presenting a relatively higher number of sightings for the 2010 survey, seems to host a stable turtle high usage area in its central and deeper portion.

Turtle Potential Biological Removal
The PBR calculated on the minimum population abundance (i.e., uncorrected estimates) ranged between 30.48 and 1609.14 individuals for a recovery factor of 0.5 and between 60.96 and 3218.27individuals for a recovery factor of 1.The highest values were obtained for the Sicily Strait monitored in winter 2016.As expected, being the PBR proportional to the estimated abundance, the same patterns emerge when considering the upper and lower confidence intervals for the uncorrected estimates as well as the adjusted estimates.Table 3 summarizes the calculated PBR values for each study area and recovery factor.Table 3. Summary of PBR calculations for each survey obtained using both uncorrected (N) and corrected (Nc) abundance estimates, lower (L95%CI) and upper (U95%CI) 95% confidence intervals of uncorrected estimates.F = recovery factor.Estimates are not presented for the Pelagos Sanctuary winter survey carried out in 2009 due to the small estimates.* Areas surveyed as part of the Extended Pelagos Area.* Surveyed as part of the Extended Pelagos Area.

Discussion
The conservation status of Mediterranean Sea loggerhead turtles was recently downlisted and the species is currently classified as Least concern [8].Nonetheless, high bycatch rates reported across the entire basin are still considered unsustainable and represent a serious threat to the persistence of the species in the region [10].While consistent research effort has been allocated to understand loggerhead turtle ecology in the region, most of this research has been carried out in coastal waters.Aerial surveys are useful for assessing large, highly mobile marine vertebrates.However, these surveys usually have relatively restricted temporal coverage, which limits the range of hypotheses that can be tested.Our surveys were conducted during the summer and winter months, with only one area, the Pelagos Sanctuary, surveyed in both seasons.This spatial and temporal uneven distribution of effort prevents meaningful comparisons of seasonal distribution of loggerhead turtles in the monitored area and the ability to evaluate shifts in the species occurrence and assess seasonal trends in abundance and density.However, our study provides new knowledge on the presence of loggerhead turtles in offshore areas which are usually scarcely monitored, although crucial to assess the bulk of the population.This information is particularly relevant to inform conservation and management activities and represent an important first contribution towards the MSFD criteria related to bycatch, abundance and range [66].

Loggerhead Turtle Abundance-Results and Limitations
Three species of sea turtles occur in the Mediterranean Sea.While the morphology, size and coloration of leatherbacks strongly differ from those of loggerhead turtles, thus preventing species misidentification, the characteristics of green turtles might render species identification difficult during aerial surveys.However, green turtles primarily occur in the Levantine Basin and, in general, in the eastern Mediterranean region [7].Accordingly, given the geographic extent of our surveys, we considered all sightings to be of loggerhead turtles.
The Mediterranean hosts turtles from different rookeries [67], and a mixture of surface currents and habitat selection [41] leads to a clear pattern of usage of foraging grounds by juvenile sea turtles from these rookeries.In particular, the areas monitored seem to be inhabited primarily by turtles of Mediterranean origin from Lebanese and Greek rookeries [67][68][69].We therefore assume that the estimates presented here concern turtles belonging to the Mediterranean Regional Management Unit.
This study highlights an overall north-to-south and west-to-east gradient in the density and abundance of loggerhead turtles across the entire monitored area.This gradient can be explained in part by the oceanographic characteristics of the study areas that, in turn, affect prey distribution and the occurrence of turtles [70], but also by the influence of water temperature which might push some animals to move towards or aggregate in warmer southern regions during the winter [51].The highest occurrence of turtles in the southern and eastern strata might also result from surface circulation distributing juvenile individuals from nesting areas to the feeding grounds [68].Animals originating from rookeries in Greece spread, in fact, in both the Sicily Strait and the Ionian Sea, which are described as important foraging grounds [57,71].The same pattern is recognized to occur for the individuals originating from Libya, which are then hosted in the Ionian and in the Tyrrhenian Seas [68] and by animals originating from Greece, Crete and western Turkey, which can eventually reach the Provençal coast and waters off France [69].Finally, it is important to consider that given the known presence of loggerhead turtles of Atlantic origin in both the Strait of Sicily and the Ionian Sea [72,73], we cannot exclude that a portion of the animals observed in these two areas do not belong to the Mediterranean RMU and have somewhat contributed to the higher abundance in these two sectors.
While most of our surveys were conducted in different years and seasons, some areas were monitored more than once in different seasons and/or years, allowing for comparisons.Estimates (Table 2) for the Pelagos Sanctuary show that, overall, this portion of the Mediterranean Sea hosts the lowest number of turtles across the monitored areas and that numbers increase in the summer months, with estimates being one order of magnitude higher in summer.For the summer, despite slightly higher estimates derived from the 2009 survey, numbers are comparable and show a constant presence of loggerhead turtles in the area, although at low densities.This is also true for the northern Tyrrhenian Sea, where no changes between the two surveys are apparent.Turtles also seem to occur consistently in the southern Tyrrhenian Sea, with an evident higher abundance in 2014.This difference could be due to the fact that the first survey carried out in the area in 2010-1011 was interrupted and then resumed after 16 weeks due to rough weather and sea conditions, which might have allowed for a redistribution of the animals in the area and affected their presence and the following estimates.
The abundance estimates have been corrected for availability bias only.Our choice of using distinct correction factors is based on the fact that in the Mediterranean, tagging effort has been unevenly allocated across the basin [50][51][52][55][56][57]] and a single overall correction factor could have led to strongly biased estimates.Considering this, we stress the need for further systematic research into the diving patterns of loggerhead turtles across the basin allowing for the development of robust correction factors for abundance estimates in this species.
Loggerhead turtles have complex life history traits associated with dramatic ontogenic changes [74,75] that can affect the detectability of animals and bias estimates.Overall, in this study, roughly 83% of the monitored waters extended beyond the continental shelf (depth> 200 m; Figure S2) and variability exists amongst study blocks.Accordingly, the imbalance in the time spent searching in each habitat might have led to an imbalance in the proportion of animals sampled in each life stage.In this light, although difficult during multi-species aerial surveys, future monitoring programs should explore the feasibility of a different design that takes into consideration the ecological aspects related to ontological stages of loggerhead turtles.
Another potential source of bias in our estimates is the fact that some smaller turtles could have been missed given the flight altitude and speed.Recent estimates of surface litter in the Mediterranean [76] used a minimum detectable size of 30 cm for litter items to derive estimates of macro-litter abundance.In our surveys, we estimate that loggerhead turtles of about 20 cm of Curve Carapace Length (CCL) could be spotted.This assumption is based on the fact that the probability to detect a turtle does not only depend on its carapace size but also on the fact that animals' flippers and head were usually spotted, increasing the overall observable size and on the fact that turtles were commonly observed in motion while surfacing or starting a dive, again increasing their detectability.Evidence suggests that Mediterranean loggerhead turtles of ≤30 cm CCL are within the first four years of life and specifically that turtles of about 20 cm CCL are about 1.5 years old [77,78].Based on this assumption and on the knowledge of loggerhead turtle population demographic characteristics and age structure [12], we consider missing a relatively small portion of the turtles.

Sea Turtle Usage Areas
This study shows a striking difference in the presence of loggerhead turtles across the studied area, with strong latitudinal and longitudinal gradients in the occurrence and distribution of sightings (Figures 2-4).This gradient reflects the general gradient of mean annual sea surface temperature in the Mediterranean [79] and shows that sea turtles are not geographically randomly distributed but tend to occupy preferred temperature ranges which vary seasonally [80].
The distribution of the sightings highlights a strong depth stratification, with a small fraction of the sightings occurring along the continental shelf (depth ≤ 200 m; n = 556, 9.1%) and most of the sightings occurring at depths between 200 and 1000 (n = 1733, 26.2%) and between 2000 and 3000 (n = 1761, 26.6%) meters (Figure 5).This imbalance in the occurrence of sea turtles in coastal and oceanic habitats strongly matches habitat availability in the monitored areas.However, the depth stratification within oceanic habitats might reflect in part the exchange of turtles between the coastal and oceanic habitats (peak of occurrence between 200 and 1000 m of depth) and in part the fact that hatchlings drift to deep pelagic waters after incubation and remain there for several years without coming to shore (peak at depth greater than 2000 m).While it is well known that turtles in the early juvenile oceanic life-stage represent most of the population and this can explain the high occurrence of sea turtles in the deepest portions of the monitored areas, little knowledge exists on the rates of movements between neritic and oceanic habitats and the rates of loggerhead turtles relocating to nearshore habitats following the pelagic stage.The location of high usage areas in the Sicily Strait and the south-western portion of the Ionian Sea strata, and to a lesser extent the southern Tyrrhenian Sea, reflects the distribution of both stable and sporadic nesting sites for the species [7,82,83] as well as the location of known foraging grounds for juvenile loggerhead turtles [7].However, while the Sicily Strait and the easternmost portion of the southern Tyrrhenian Sea are well known pelagic summer grounds, our results suggest a high year-round occurrence of loggerheads in these areas.This, in turn, further complements previous studies highlighting the Tunisian shelf as an important hot-spot for sea turtles [84].The highest presence of turtles in these areas could also be explained by the fact that the Sicily Strait, in part the Ionian Sea, and more in general the north African coast represent a preferred migratory corridor for turtles that flock from the eastern Mediterranean rockeries towards oceanic grounds in the central basin [51].This pattern of occurrence of high usage areas in the deeper portions of all the monitored areas, irrelevant of the year or season, strongly reflects the surface circulation of the Mediterranean [85][86][87][88][89], where the north and southern branches of the Algerian Current and the associated gyres in the Tyrrhenian Sea and the Sicily Strait not only aggregate food resources and therefore affect the presence of turtles, but also tend to push and aggregate smaller individuals which drift with the current [90][91][92].
Finally, when comparing our results with a previous Aqua Map model of sea turtle distribution in the Mediterranean Sea developed by [93], with the exception of the Sicily Strait area, very little overlap exists between the predicted distribution by Coll and colleagues [93] and our results.However, these differences should not be considered as a disagreement between the results of the two studies but most likely the effect of different methodological approaches.The predicted turtle distribution by Coll et al., [93] is primarily based on nesting beach data and does not reflect the actual presence of animals at sea.This further calls attention to the necessity of integrating different data sources to obtain a robust and overall picture of loggerhead turtle presence in the region.

Potential Biological Removal (PBR)
The calculation of PBR values is a way to identify management goals for a given population on the base of simple assumptions, provided that baseline abundance estimates of the target species are available.
The PBR values obtained in this study for both corrected and uncorrected estimates and associated confidence intervals are lower than the bycatch rates reported in portions of the monitored area [13], supporting previous concerns that the mortality of loggerhead turtles in commercial fisheries alone is likely unsustainable in the Mediterranean Sea [10].While these findings are relevant to management and conservation, caution is needed when interpreting PBR values.When comparing our results with the bycatch rates reported by [13], it is evident that the surveys were conducted in different seasons and the study areas do not fully overlap.While, to ease comparisons, further mathematical exercises are possible to extrapolate abundance estimates and PBRs to areas not directly surveyed, this approach is discouraged.Abundance estimates and thus PBR values extrapolated to non-surveyed areas can be strongly biased due to the fact that the correlates between species densities and covariates in extrapolated areas are unknown [94,95].It is also important to stress that bycatch rates are often biased (i.e., under-reported [13]) and this adds further difficulties in understanding the actual contribution of commercial fisheries to loggerhead turtle mortality in support of management and conservation.Finally, the presence in the Mediterranean Sea of turtles of different origins introduces a further level of complexity in understanding the actual size of the Mediterranean loggerhead turtles' RMU and how it is affected by bycatch.
All these elements once again stress the need for future surveys covering the full geographic extent of the Mediterranean loggerhead turtles' RMU and the integration of different approaches and methodologies to resolve uncertainty in abundance and obtain range-wide PBR estimates.

Conclusions
Regional sea conventions available on a Mediterranean scale are necessary tools to ensure that a robust legal framework exists to support management decisions.Nevertheless, it is important to emphasize that ad hoc efforts are urgently needed to bridge the gap between scientists and the fishing industry to ensure a rapid progression into policy and effective conservation.Despite the limitations acknowledged in the above paragraphs, the information provided here poses the basis for more targeted bycatch management efforts.Given the broad spatial scale of this study, our results could be considered in regional planning efforts such as the EBSA [96] and the EcAp processes, as well as the EU MSFD.Our results also support the Important Marine Turtle Areas (IMTAs) process for loggerheads in the Mediterranean region [97].
The abundance estimates, the derived PBR and the loggerhead turtle UDs may contribute to establish, at sub regional scale, the threshold values requested for the Criteria D1C1 (mortality rate from incidental bycatch), D1C2 (population abundances are not adversely affected by anthropogenic pressure), D1C4 (distributional range in line with ecological conditions) and D1C5 (necessary habitat extent), under the umbrella of the MSFD.

Supplementary Materials:
The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/d15121159/s1, Figure S1.Overall monitored area between 2009 and 2916 with highlighted overlapping portion of Marine Subregions as delineated under the Marine Strategy Framework Directive.Geodatabases used to produce the map available at https://www.eea.europa.eu/data-and-maps/data/europe-seas-1; Figure S2.Depth of the overall surveyed area classified in coastal (depth between 0 and 200 meters) and deep waters (depth greater than 200 meters).Depth reclassification was based on the ETOPO1 Global Relief Model [81]; Figure S3.Detection functions for Conventional Distance Sampling (CDS) analysis for each survey and area.Details on selected models are provided in Table 3. A= Extended Pelagos Area (2010); B= Sicily Channel (2016); C= Central Tyrrhenian Sea (2010); D= Central Tyrrhenian Sea (2013); E= Ionian sea (2010); F= Pelagos Sanctuary (2010); G= Southern Tyrrhenian Sea (2010-2011); H= Southern Tyrrhenian Sea (2014).For the detection functions and related information for the 2009 winter and summer surveys refer to Lauriano et al., [35].

Figure 1 .
Figure 1.The area monitored between 2009 and 2016 with the survey strata and transects resulting from the Distance software design.Survey design for the winter and summer 2009 surveys in the Pelagos Sanctuary are available from Lauriano et al. [35].The map inset with the extent of the monitored area is highlighted in red.

Figure 2 .
Figure 2. Summer surveys; for each study area and year, the map on the left shows loggerhead turtle sightings, the map in the center shows the obtained UDs and the map on the left highlights the UD's 25th and 75th percentile.

Figure 3 .
Figure 3. Winter surveys; for each study area and year, the map on the left shows loggerhead turtle sightings, the map in the center shows the obtained UDs and the map on the left highlights the UD's 25th and 75th percentile.

Figure 4 .
Figure 4. Violin plots with included boxplots of turtle sightings' latitude (left) and longitude (right) per surveyed area showing north-south and west-east gradient in the occurrence of sightings.Wider sections of the violin plots indicate a higher number of sightings at a given latitude or longitude while the skinnier sections indicate a lower occurrence.Boxplots show lower (Q1), median and upper (Q3) quartiles and outliers (black dots).
3), followed by the Extended Pelagos Area survey (n = 33,981; %CV = 15.4) and the Pelagos Sanctuary survey (n = 1905; %CV = 16.2).The latter was monitored in 2009 and 2010-as part of the Extended Pelagos Area-with an estimated abundance being higher in 2009.The northern Tyrrhenian Sea was surveyed in 2010-as part of the Extended Pelagos Area-and 2013; it presented the highest abundance in 2010.Amongst the winter surveys, the highest abundance was estimated for the Sicily Strait surveyed in 2016 (n = 100,571; CV = 4%), followed by the Ionian Sea (n = 35,434; CV = 16.7%)monitored in 2010.The Pelagos Sanctuary, monitored in 2009, conversely, shows the lowest abundance (n = 237; CV = 34.33%) of loggerhead turtles (Lauriano et al., 2011 [35]), with values being about 300 times lower than the Sicily Strait.In the southern Tyrrhenian Sea, surveyed twice in 2010-2011 and in 2014, the estimated abundance and densities are substantially higher in the latter period (2010-11 survey: n = 17,972; CV = 20.7%;2014 survey: n = 33,217; CV = 10.9%).

Figure 5 .
Figure 5. Percent frequency distribution of loggerhead turtle sightings at different depths.Bin widths are uneven to show the separation between turtles' neritic and oceanic habitats.The width of the last bin is based on the highest depth as extracted from the ETOPO1 Global Relief Model [81].