Open Access

This article has an erratum: []

Aquat. Living Resour.
Volume 36, 2023
Article Number 21
Number of page(s) 11
Published online 25 July 2023

© K. Papadopoulo et al., Published by EDP Sciences 2023

Licence Creative CommonsThis is an Open Access article distributed under the terms of the Creative Commons Attribution License (, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

1 Introduction

Animal movement is fundamental to life and shapes many ecological and evolutionary processes (Nathan et al., 2008; Cooke et al., 2022). For instance, spatial behaviour affects population dynamics, connectivity, exposure to threats and the ability to cope with environmental changes (Morales et al., 2010). At the same time, the movement of organisms is influenced by a myriad of biotic and abiotic factors at multiple temporal and spatial scales (Lédée 2015; Cooke et al., 2022). Understanding the causes and consequences of movement is therefore of high interest in ecological and evolutionary research (Shaw, 2020).

From a practical standpoint, a comprehensive understanding of the spatial ecology of marine organisms is essential to assess and inform marine conservation policies and management plans (Crossin et al., 2017; Hays et al., 2019). Combining information on the movement ecology of highly migratory species with the spatial distribution of fishing efforts has revealed a high and concerning overlap between fish abundance and fishing pressure (Queiroz et al., 2019). A growing number of marine protected areas (MPAs) are now being designed and evaluated based on the extent of species home ranges and habitat selection (MacKeracher et al., 2018; Gilmour et al., 2022). Spatial ecology studies have revealed essential habitats for marine animals such as foraging grounds (Warwick-Evans et al., 2018) or spawning and nursing areas (Hays et al., 2019). Finally, innovative research using movement monitoring succeeded in estimating key demographic parameters such as mortality rates for neonates and juveniles (Heupel and Simpfendorfer, 2002) as well as survival rates from fisheries discard (Morfin et al., 2019; Alonso-Fernández et al., 2022).

Over one-third of species of the class Chondrichthyes are threatened by the loss and degradation of habitat, climate change, pollution and overfishing (Dulvy et al., 2021). Skate (Rajidae) populations have been experiencing a severe decline worldwide as a result of habitat loss and overexploitation (McPhie and Campana, 2009, Simpson et al., 2020, Dulvy et al., 2021). Most skate species have low population growth due to their slow development rates and late sexual maturity (Licandeo et al., 2006; Pardo et al., 2016). Consequently, the majority of skate species are unable to withstand high levels of fishing pressure (McPhie and Campana, 2009; Dulvy et al., 2014) and are considered among the most vulnerable groups of fish (Dulvy et al., 2000). Current population declines of skates have resulted in the inclusion of many Rajidae species on the IUCN (International Union for the Conservation of Nature) Red List. As of today, 158 species of Rajidae are present on the IUCN Red List with 36% of species (ray and skates) classified as threatened (n = 220 of 611 species) (Dulvy et al., 2021; IUCN, 2022). Despite a few studies focused on skate ecology (Neat et al., 2014; Sousa et al., 2019; Simpson et al., 2020), many aspects of their life cycle still remains unknown. In fact, 13.3% of all skate species present on the IUCN Red List are still considered “data deficient” (IUCN, 2022) and the populations of skates in Europe are considered data-limited stocks without analytical stock assessments (ICES, 2022). Filling these ecological knowledge gaps is a pressing need if we want to counteract the global diversity loss (Joppa et al., 2016).

Technological advancements have made it possible to track marine animals in conditions that would otherwise be impossible to achieve, thus fostering the rapid development of the field of movement ecology (Lennox et al., 2017; Lowerre-Barbieri et al., 2019). In systems ranging from lakes and rivers to the open ocean, acoustic telemetry is the most used method to track submerged aquatic organisms (Hussey et al., 2015). The use of underwater acoustic telemetry allows for an in-depth understanding of fish movement ecology, such as home range (Leeb et al., 2021), activity (Bohaboy et al., 2022) as well as diel and seasonal differences in movement patterns during the year (Williams-Grove and Szedlmayer, 2016). This represents a particularly suitable technique to meet the actions of the “Sustainable Development Goal 14” of the United Nations (Alós et al., 2022). In a context of ever-growing anthropogenic disturbances such as habitat degradation, migration barriers and climate change, furthering our understanding of the drivers of animal behaviour is of the utmost importance to increase the effectiveness of conservation efforts (Hays et al., 2019; Lowerre-Barbieri et al., 2019).

The thornback skate, Raja clavata (Linnaeus 1758), is amongst the most common and widespread skates in the northeast Atlantic and Mediterranean Sea (Ellis, 2016). This bottom-dwelling and coastal species can be found from the South of the Arctic Circle (Iceland, Norway) to the east Atlantic coast of south Africa (Stehmann, 1995). Like most skates, the thornback skate spends the majority of its time buried in fine sediment (Albert et al., 2022). Often occurring at depths ranging from 0 to 60 meters, thornback skates can inhabit areas as deep as the upper limit of the continental slope (300 m) (Whitehead et al., 1986; Last et al., 2016; Trenkel et al., 2022). Characterised as an opportunistic, mobile and active predator, thornback skates feed on a wide variety of prey (mainly invertebrates and some fish), thus granting it a wide foraging area (Ellis et al., 1996; Farias et al., 2006). As with all Rajidae species, female thornback skates spawn egg capsules which they bury or attach to the substrate (Maia et al., 2015). Classified globally as “Near Threatened” by the IUCN Red List (Ellis, 2016), the thornback skate is among the most frequently captured skates by commercial fisheries in northwest Europe (Santos et al., 2021). In 2019, the global catch of thornback skate reached 6,874 tonnes, most of them being captured as bycatch of trawls and gillnets (FAO, 2021). The thornback skate is a common species in European multispecies and multi-gear fisheries, including the small-scale fishing sector (Bañón Díaz et al., 2008; Figueiredo et al., 2020), in partially protected areas (Di Lorenzo et al., 2022). Thornback skates rarely escape trawl nets because of their large size and thorns; this, coupled with the slow growth rate and low fertility of skates, makes overfishing a great threat to their populations (Ellis, 2016). It is important to note that, the conservation status of R. clavata appears to show signs of improvement in recent years, with increased biomass and indications of low exploitation levels in the northeast Atlantic (ICES, 2022). However, variations and uncertainties in different assessments, emphasize the need for continued monitoring and management measures to ensure the long-term sustainability of the species. Although the general latitudinal and bathymetric distribution of the species is understood, there have been virtually no studies assessing the fine-scale behaviour of the species in which may however be relevant for understanding the applicability of conservation measures such as marine protected areas. In this study, we filled this knowledge gap by using fine-scale positioning based on acoustic telemetry to track the spatial ecology of thornback skates in the Cíes Islands, a small MPA in the northwest of the Iberian Peninsula. Our objectives were to (1) identify the residence patterns in the study area; (2) assess the timing, duration and direction of the excursions out of the study area; (3) quantify the temporal variation in space use. This study allows inferences about thornback skates in the Cíes archipelago that were previously not possible.

2 Methods

2.1 Study array

This study was carried out in the Cíes islands, at the mouth of the Ría de Vigo (Galicia, northwest Spain), between October 2020 and June 2022. The Cíes islands are part of a partially protected area called the Parque Nacional Marítimo-Terrestre das Illas Atlánticas de Galicia (PNMTIAG). This archipelago is made up of several islands and islets that together cover an area of 31 km2 (Figs. 1a–c). Inside the PNMTIAG recreational fishing is prohibited and different uses are regulated (Xunta de Galicia, Conselleira de Medio Ambiente, Territorio y Vivienda, 27/12/2018). Waters surrounding the Cíes islands are subject to upwelling events, rendering them particularly productive and therefore valued as important fishing grounds by the small-scale local fishers (Arístegui et al., 2004; Broullón et al., 2023).

An array of 22 ©Innovasea (formerly Vemco) VR2W acoustic omni-directional receivers was deployed in the study area. Receivers were mounted at the top of auger anchors (140 cm high and screwed 60 cm deep into the sandy substrate), at depths ranging from 3.3 to 13.1 m (Villegas-Ríos et al., 2013) (Fig. S1). Two ©Innovasea reference transmitters (V13 and V16) were placed at fixed positions within the receiver array to assess potential environmental effects on the detection patterns and the error associated to skate positions (Payne et al., 2010). Following range tests (Leeb et al., 2021), the distance between stations was set at ∼150 m to ensure adequate coverage and overlap between receiver's detection ranges. Assuming an average detection range of ≃150 m (corresponding to 50% of the proportion of the detections received during the range test (Leeb et al., 2021)) the area covered by the array would be 0.58 km2.

thumbnail Fig. 1

(a) Location of the study array in the Iberian Peninsula (red square). (b) Position of the study area (red shaded area) within the national park (green polygons). (c) Map of the study area showing the location of the ©Innovasea acoustic receivers (coloured dots), reference tags (red triangles) with the temperature data logger (green circles and triangles). (d) Detailed map of the study area displaying ©Innovasea acoustic receivers divided into four sectors, reference tags and bathymetry.

2.2 Skate tagging

All thornback skates (n = 14) were caught by a small commercial boat using “palangrillo” (Galician local name for small bottom long-lines) baited with pilchard or squid near or within the acoustic array. The fishing gear used in the study comprised sets of four and six pieces containing approximately 160 hooks per piece. Longline sets covered an average distance of 2.26 ± 0.88 km and the soak time of the experimental fishing trips varied between 26 min and 4 h and 15 min. Eleven thornback skates were externally tagged (anchored with monofilament) in the pectoral fin (Fig. S1) with ©Innovasea V13P-1x (lifetime: 855 days; signal transmission delay: 80–160 s) and three with ©Innovasea V9AP-2x (lifetime: 453 days; delay: 80–160 s) transmitters equipped with pressure sensors. Fish were tagged on two dates: 08th of October and 19th of November 2020. Each tagged skate was sexed and the disc length (DL) measured to the nearest cm; individuals were subsequently returned to the sea as quickly as possible (<5 min). External T-bar tags (©Floy Tag) were attached to the pectoral fin of the individuals to prevent possible repeated tagging in subsequent samplings and to enable fishers and divers to report their recaptures (Fig. S1).

A dead thorback skate was equipped with an acoustic transmitter (©Innovasea V13P) and released inside the telemetry array to determine the effects of currents, waves, or scavengers on the movement of a dead individual. As a result, we were able to identify a “dead” pattern, which could be used to assess whether a tagged skate had died inside the array (Villegas‐Ríos et al., 2020; Alonso-Fernández et al., 2022) (Fig. S2). This is important to make sure that the behavioural variables are only estimated when the fish are alive (Villegas‐Ríos et al., 2020).

This study complied with animal welfare regulations of the regional government (Xunta de Galicia) starting on the 14th of November 2019. It follows the Experimental Animal Project Authorization: ES360570202001/19/FUN01/BIOL AN.08/AAF01.

2.3 Data treatment and analysis

Detection data were offloaded from acoustic receivers every 6 months starting in September 2020 until June 2022. The data was stripped from single detections occurring within 24 hours at receivers as they were considered false detections (Meyer et al., 2007). The fate and “fate date” for each tagged skate were assigned by examining plots of position and depth over the study time based on Centres of Activity (COAs) (Villegas‐Ríos et al., 2020; Alonso-Fernández et al., 2022) (Fig. S3). Behavioural analyses only included skates that were alive and within the acoustic array for at least three days post-release to exclude individuals which died after release. Moreover, any detections post the assigned fate date were filtered from the analysis (Fig. S3).

2.4 Residency

The residence index (RI) for each skate was computed as follows:

where DD corresponds to the number of days an individual was detected and TD to the total number of days between the tagging date and the end of the study. RI range of values goes from 0 (never detected in the array) to 1 (continuous residency in the array) (Papastamatiou et al., 2010; Espinoza et al., 2011).

2.5 Characterizing excursions

To determine possible preferred routes used for entering and exiting the study area, we divided the array into four sectors (Fig. 1d). Excursions out of the study area were identified as instances when an individual was absent for more than 24 h. For all excursions, we retained: (i) the sector where the last detection of the skate leaving the array was recorded, (ii) the sector which recorded the first detection of the skate returning into the array and (iii) the duration of the excursion. We defined four sectors mirroring the four cardinal points: ‘north’ − suggesting a movement towards the northern side of the archipelago; ‘east’ − suggesting a movement towards the inner part of the Ría de Vigo; ‘south’ − which suggests that the skate was heading towards the southernmost island; and ‘west’ − which is linked with the channel that separates the two main islands and ultimately leads to the continental shelf. Excursions were classified based on their duration as: (1) ‘short-term’ − when the time spent outside the study area was less than seven days; (2) ‘medium-term’ − when it was comprised between seven days and one month; and (3) ‘long-term’ − when it was longer than one month. Lastly, the degree of consistency with which fish used different sectors to exit and re-enter the study area was examined by calculating the number of excursions that had the same exit and entry sector.

2.6 Activity space

Activity space was defined as the 95% Kernel Utilization Distribution (KUDs) of each skate (i.e. a central area or volume within which an animal is 95% likely to be found) (Powell, 2000). KUDs were computed on a weekly basis using “adehabitatHR” packages in R (Calenge, 2006). We did not estimate KUDs for weeks with positions in less than four days (consecutive or non-consecutive) to eliminate bias from activity space computation based on weeks with few detection days (Leeb et al., 2021).

2.7 Environmental variables

Oceanographic conditions within and near the acoustic telemetry arrays were monitored throughout the study period. Sea levels were acquired from the tide gauge “Vigo 2” located at 42°1424N; 8°4348W (Prontuario instalaciones, 2022). Sea bottom temperature (‘Temp’) was recorded every half hour using a combination of data loggers (©Star: ODDI DSTcenti-T) and Thelma acoustic receivers only used as temperature loggers and covering different depths and areas of the study array (Fig. 1d).

2.8 Statistical analysis

Along with observed environmental conditions, the day of the year (DOY: 1-366) was included in the models to explain seasonal effects. Similarly, week of the year (WOY: 1-53) was added to describe the yearly cycles in skate behaviour. Day and night times based on sunrise and sunset at 42° 12 47.6634 N; −8° 543.9522 W were calculated using the “suncalc” package in R (Thieurmel and Elmarhraoui, 2022).

We used Generalized Additive Mixed-effects Models (GAMMs) implemented in the R package “mgcv” to assess the relationship between independent variables and skate probability of presence as well as activity space (Wood et al., 2014, 2017; Li and Wood, 2019). In all models, individuals were handled as a random effect (via random intercept), an autoregressive term of order one (corAR1) was used since observations were made throughout time in a sequential manner (Dormann, 2007). The addition of the autoregressive term when describing weekly integrated activity space was based on a previous study demonstrating that employing weekly replicates of behavior yielded less biased repeatability estimates (Villegas-Ríos et al., 2017). Common fixed effects in the model included sex (male, female) and disc length (cm). Maximum likelihood criteria were used to fit GAMMs and a backward (decreasing number of variables) selection method was used to construct them. Non-parametric smoothing functions s were fitted with four knots to model the non-linear effect of sea bottom temperature on a day t (‘Tempt’) on the residency. Day of the year (‘DOY’) and week of the year (‘WOY’) were fitted as non-parametric smoothing functions s with four knots and cyclic cubic splines.

3 Results

A total of 381,884 detections (29,344 post filtration) were retrieved from the receiver array at the end of the study period, on the 14th of June 2022. After filtering out the detections of two dead individuals that remained dead inside the array during most of the study (DESTAC-SPP-20-03 and DESTAC-SPP-20-04), we retained 29,344 detections for analysis. The number of tracking days varied greatly among skates ranging from just 3 to 160 days (Tabs. S1 and S2).

3.1 Residency

Of the 14 tagged individuals (mean disc length = 34.79 cm, range: 22–49 cm), only nine were present in the array for three or more days after the tagging date and thus included in the analyses (Tab. S1). Amongst those nine skates, two individuals were classified as dead during the study period (Fig. 2). The residence index was low, ranging from 0.005 to 0.260 (mean RI = 0.073) (Tab. S2).

The probability of presence of R. clavata in the study area was influenced by both biotic and abiotic factors. R. clavata residency in the array followed a bell curve (Fig. 3a), reaching its maximum in summer (DOY = 159). Females had a higher probability of presence in the study area compared to males (Tab. 1; Fig. 3b). Finally, there was some support (Tab. 1) for a higher probability of presence at temperatures around 15 °C with a small decrease towards higher and lower temperatures (Fig. 3c).

thumbnail Fig. 2

Abacus plot showing the daily presence of thornback skate (Raja clavata) in the study area. Days when an individual was present, are coloured in beige. The two black lozenges at the start and end of each time series represent respectively the tagging and end date of the study. Daily presence is displayed for the reference tag (red) and the dead skate (blue). Green lozenges represent death events: DESTAC-SPP-20-03 & DESTAC-SPP-20-04 on 19/10/2020.

thumbnail Fig. 3

Predicted probability of the presence of Raja clavata in the study area as a function of day of the year (a), sex (b) and sea bottom temperature (c). Black bars (b) and grey shaded areas (a, c) represent the 95% confidence interval. Black dots correspond to the raw data of probability of presence. Values used for predictions: sex = female, sea bottom temperature = 14 °C, day of the year = 260.

Table 1

Summary of the optimal generalized additive mixed-effects models investigating the (i) probability of presence and (ii) activity space of Raja clavata in the study area†.

3.2 Excursions

Over the course of the study, we identified 70 excursions. Among them, six were long-term (average duration = 175 ± 54 days), 12 were medium term (average duration = 15 ± 7 days) and 52 were short-term (average duration = 3 ± 2 days). The number of excursions per individual ranged from 0 to 26, with a mean of 7.8 excursions per individual (Tab. S2). In the 615 days of the study, individuals exited and re-entered the study area 55 times in the direction of the Ría de Vigo (East). Interestingly, from the four sectors composing the array, the sector connecting to the West (continental shelf) was never used. High consistency was observed in the sector used to exit and re-enter the study area by R. clavata during excursions, with 82% (43 out of 52) of the short terms excursions and 66% (8 out of 12) of the medium term excursions having the same exit and re-entry sector. Long-term excursions had less consistency with only 50% (3 out of 6) having the same exit and re-entry sector.

3.3 Activity space

R. clavata activity space in the study area ranged from 0.12 to 0.50 km2 (mean AS = 0.27 ± 0.13 km2) (Tab. S2, Fig. S4). Observation of raw data based on COAs indicated a sexual variation in activity space with females displaying a higher occupied space (0.33 ± 0.11 km2) than males (0.21 ± 0.08 km2) (Tab. S2). While some individuals may have utilized the entire depth range of the study array, it appeared that R. clavata had a preference for the deeper waters within the array (Fig. S4). The activity space of R. clavata had no significant diel variation, with daytime activity space (0.25 ± 0.11 km2) mirroring the activity space at night (0.26 ± 0.13 km2). Finally, results of the GAMM (Tab. 1) confirmed seasonal variation of the activity space of R. clavata, reaching its maximum in summer (0.25 ± 0.03 km2 at WOY = 29) (Fig. 5). It is important to note that, having very few individuals present in the array in summer, this estimate is likely influenced by this small sample size.

thumbnail Fig. 4

Overview of all sectors, showing the number of exits (a) and entry (b) routes taken during the study period by Raja clavata.

thumbnail Fig. 5

Predicted activity space of Raja clavata in the study area as a function of the week of the year. Grey-shaded areas represent the 95% confidence interval. Black dots correspond to the raw data of activity space. Values used for predictions: sex = female, sea bottom temperature = 14 °C.

4 Discussion

The analysis of data retrieved from nine acoustically tagged R. clavata showed significant variations in their spatial behaviour at different timescales, influenced by both biotic and abiotic factors. On average, the probability of presence in the study area for R. clavata was very low (mean RI = 0.073). The likelihood of presence was mostly affected by sex, with females being more likely to be present than males. The excursion patterns indicated a strong connection to inshore waters, specifically the Ría de Vigo. Our results showed a major increase in activity space in summer. This study provides an important baseline for understanding and linking the ecology of the vulnerable thornback skate with possible conservation actions.

The first key result of our study is that R. clavata spent very little time inside the study area. Most of the individuals left after five days or less and were not detected again. This alone suggests that small MPAs of just a few square kilometres are probably not enough to protect this species over long periods of time. It is important to point out the pronounced difference in residency between Scyliorhinus canicula RI (0.27) and R. clavata RI (0.073) in the same study area (Papadopoulo et al., 2023). An alternative explanation for the limited use of the study area by R. clavata may indicate a higher level of residence in a nearby area or to sporadic dispersal movements. For the few individuals that were detected for more than 30 days (n = 4), a seasonal pattern of presence was observed, with a higher presence in summer. In fact, two individuals were detected in two consecutive spring-summer seasons after being absent in winter. Although based on a small number of individuals, this may suggest a seasonal use of the study area for at least part of the population, matching the behaviour observed in another skate species in the same study area (Leeb et al., 2021). The reason why individuals of these species spend more time in the study area in spring and summer remains unknown, but it might be related to their life cycle (Chevolot, 2006). Prior studies have indicated that skates may exhibit seasonal movements, with individuals moving from deeper waters during winter, to shallower waters during spring and summer, where they are believed to mate and spawn (Walker et al., 1997; Hunter et al., 2005a). Although no egg cases of R. clavata have been reported inside the array, these seasonal patterns in skate movement could account for their migration in and out of the study area.

The second main result is the observation that space use varied seasonally. R. clavata had the highest activity space inside the array in summer and the smallest in spring. Seasonal variation in the home range of R. clavata has been observed in different studies, however, they focused on the variation in its vertical range (Hunter et al., 2005a, 2005b; Cabral, 2014). In these studies, the breeding cycle appeared to be the main driver of variation in R. clavata vertical range, thus affecting its home range (Hunter et al., 2005a, 2005b). Within the study area, another species of skate (Raja undulata) displayed a seasonal variation in activity space, although with a peak in spring and a steady decrease towards winter (Leeb et al., 2021). The main drivers of this variation were also associated with breeding cycles with R. undulata using sheltered, shallow habitats as nursery grounds. The ultimate drivers of variation of activity space in R. clavata could be related to biological needs and cycles (e.g. feeding or mating) (Hunter et al., 2005b) but compelling evidence is not available. It is worth mentioning that this study did not investigate any evidence of sex/size-related segregation patterns that could account for the observed variation in space use. Therefore, any association between seasonal variation in space use and the reproductive cycle of R. clavata should be approached with caution.

The third key result is about the patterns of excursions from and to the study area. The telemetry data revealed the many instances (70) during when R. clavata individuals travelled outside the array to later return. These observations combined with the overall low residency index suggest that R. clavata true home range extends far beyond the limit of the study area and is in agreement with previous studies conducted on two different elasmobranchs species at this site (Leeb et al., 2021; Papadopoulo et al., 2023). A previous mark-recapture study found that most adult R. clavata were recaptured within 37 km of the release site (Chevolot et al., 2006). Interestingly, R. clavata demonstrated high consistency in the sector taken to exit and re-enter the array with the majority of excursions heading to the Ría de Vigo. This result mirrors the excursion patterns of S. canicula at the same site (Papadopoulo et al., 2023) and suggests importance of the inshore waters of the Ría de Vigo for coastal elasmobranchs. Notably, although R. clavata appears to follow a specific route for entering and exiting the array, it is important to acknowledge that our approach points to the most probable direction taken after leaving the array, but it does not provide confirmation of the final destination of the excursions.

Our findings establish a significant relationship between sea bottom temperature, sex, day of the year and the probability of presence, thus contributing significantly to our understanding of animal movements and the ecological drivers of their behaviour. However, the amount of variance explained by all those factors was low (adjusted R2 = 0.132). This suggests that other unaccounted biotic factors such as life cycle and abiotic factors such as habitat type, salinity or tides are also likely drivers of the presence of animals in the study area, as it has been observed in other studies (Collins et al., 2008; Heupel and Simpfendorfer, 2014; Elston et al., 2022).

Despite being limited in sample size, temporal and spatial scope, this study has provided a preliminary insight into the spatial ecology of the thornback skate in Galicia. In order to effectively manage mobile elasmobranchs, it is crucial to gather information on the factors that influence their use of space and length of stay in a certain area. In this study, we have identified two distinct patterns of movement exhibited by the thornback skate. (1) A peak of presence in the study area in summer, when their space use is at its maximum. Similar to R. clavata, R. undulata displayed a peak in its probability of presence during the summer, while its highest activity space occurred in late spring (Leeb et al., 2021). In contrast, S. canicula had the lowest probability of presence in summer and did not exhibit significant variations in its activity space throughout the year (Papadopoulo et al., 2023). (2) Mirroring S. canicula, R. clavata prefered sector for exiting and reentering the array was the sector pointing towards the Ría de Vigo. These findings highlight distinct patterns of presence and activity space among three elasmobranch species in the same area, emphasizing the importance of considering species-specific behaviours and environmental factors in understanding their distribution and ecology. This knowledge can be useful to implement conservation strategies to better protect the thornback skate.

Further research is needed to explore the migratory patterns of thornback skates once they leave the array in the direction of the Ría de Vigo. It is important to understand their foraging habits and reproductive cycle in order to identify other areas and environmental conditions that are crucial for the completion of their life cycle. Continued research into the behavioural ecology of thornback skates is therefore necessary to allow the implementation of effective management plans around this commercially important and near threatened species of skate.

Declaration of interest statement

The authors report there are no competing interests to declare.


This work was supported by European Union's Erasmus program, Fundación Biodiversidad, DESTAC and CONECTEE.


K. Papadopoulo was supported by the European Union's Erasmus+ programme. This study was funded by the projects DESTAC and CONECTEE with the collaboration of the Fundación Biodiversidad, from the Ministerio para la Transición Ecológica y el Reto Demográfico (Spanish Government), through the Pleamar program, cofounded by the FEMP and the research program of Red Parques Nacionales. The acoustic telemetry array was supported by the ATLAZUL project (0755_ATLAZUL_6_E) co-funded by the European Regional Development Fund (ERDF) through the Interreg V-A Spain-Portugal Program (POCTEP) 2014-2020. DVR has received funding from the Ramón y Cajal programme from the Spanish Ministry of Science and Innovation (RYC2021-032594-I).


  • Albert L, Olivier F, Jolivet A, Chauvaud L, Chauvaud S. 2022. Insights into the behavioural responses of juvenile thornback ray Raja clavata to alternating and direct current magnetic fields. J Fish Biol 100: 645–659. [CrossRef] [PubMed] [Google Scholar]
  • Alonso-Fernández A, Mucientes G, Villegas-Ríos D. 2022. Discard survival of coastal elasmobranchs in a small-scale fishery using acoustic telemetry and recapture data. Estuarine Coast Shelf Sci 276: 108037. [CrossRef] [Google Scholar]
  • Alós J, Aarestrup K, Abecasis D, Afonso P, Alonso‐Fernandez A, Aspillaga E, Barcelo‐Serra M, Bolland J, Cabanellas‐Reboredo M, Lennox R, McGill R, Özgül A, Reubens J, Villegas‐Ríos D. 2022. Toward a decade of ocean science for sustainable development through acoustic animal tracking. Glob Change Biol 28: 5630–5653. [CrossRef] [PubMed] [Google Scholar]
  • Arístegui J, Alvarez-Salgado XA, Barton ED, Figueiras FG, Hernández-León S, Roy C, Santos AMP. 2004. Oceanography and fisheries of the Canary Current/Iberian region of the Eastern North Atlantic (18a, E). The global coastal ocean: Interdisciplinary regional studies and syntheses 14: 879. [Google Scholar]
  • Bañón Díaz R, Quinteiro Fernández R, Juncal Caldas LM,Campelos Álvarez JM, Gancedo Baranda A, Lamas Rodríguez F, Morales de la Fuente C, Ribó Landín J. 2008. Composición, distribución y abundancia de rayas (Elasmobranchii: Rajidae) en aguas costeras de Galicia. Foro Dos Recursos Mariños e Da Acuicultura Das Rías Galegas 10: 325–331. [Google Scholar]
  • Bohaboy EC, Cass-Calay SL, Patterson WF. 2022. Fine-scale movement of northern Gulf of Mexico red snapper and gray triggerfish estimated with three-dimensional acoustic telemetry. Sci Rep 12: 14274. [CrossRef] [PubMed] [Google Scholar]
  • Broullón E., Franks PJS, Fernández Castro B, Gilcoto M, Fuentes‐Lema A, Pérez‐Lorenzo M, Fernández E, Mouriño‐Carballido B. 2023. Rapid phytoplankton response to wind forcing influences productivity in upwelling bays. Limnol Oceanogr Lett. [Google Scholar]
  • Cabral S. 2014. Raja spp. spatial dynamics in the Arrábida Marine Park (Msc Thesis). University of Algarve. [Google Scholar]
  • Calenge C. 2006. The package “adehabitat” for the R software: a tool for the analysis of space and habitat use by animals. Ecol Model 197: 516–519. [CrossRef] [Google Scholar]
  • Chevolot M, Ellis JR, Hoarau G, Rijnsdorp AD, Stam WT, Olsen JL. 2006. Population structure of the thornback ray (Raja clavata L.) in British waters. J Sea Res 56: 305–316. [CrossRef] [Google Scholar]
  • Chevolot MSCOM. 2006. Assessing genetic structure of thornback ray, Raja clavata: a thorny situation? (PhD Thesis). University of Groningen. [Google Scholar]
  • Collins AB, Heupel MR, Simpfendorfer CA. 2008. Spatial distribution and long-term movement patterns of Cownose Rays Rhinoptera bonasus within an estuarine river. Estuar Coasts 31: 1174–1183. [CrossRef] [Google Scholar]
  • Cooke SJ, Bergman JN, Twardek WM, Piczak ML, Casselberry GA, Lutek K, Dahlmo LS, Birnie‐Gauvin K, Griffin LP, Brownscombe JW, Raby GD, Standen EM, Horodysky AZ, Johnsen S, Danylchuk AJ, Furey NB, Gallagher AJ, Lédée EJI., Midwood JD, Gutowsky LFG. 2022. The movement ecology of fishes. J Fish Biol 101: 756–779. [CrossRef] [PubMed] [Google Scholar]
  • Crossin GT, Heupel MR, Holbrook CM, Hussey NE, Lowerre-Barbieri SK, Nguyen VM, Raby GD, Cooke SJ. 2017. Acoustic telemetry and fisheries management. Ecol Appl 27: 1031–1049. [CrossRef] [PubMed] [Google Scholar]
  • Di Lorenzo M, Calò A, Di Franco A, Milisenda G, Aglieri G, Cattano C, Milazzo M, Guidetti P. 2022. Small-scale fisheries catch more threatened elasmobranchs inside partially protected areas than in unprotected areas. Nat Commun 13. [PubMed] [Google Scholar]
  • Dormann CF. 2007. Assessing the validity of autologistic regression. Ecol Modell 207: 234–242. [CrossRef] [Google Scholar]
  • Dulvy NK, Fowler SL, Musick JA, Cavanagh RD, Kyne PM, Harrison LR, Carlson JK, Davidson LN, Fordham SV, Francis MP, Pollock CM, Simpfendorfer CA, Burgess GH, Carpenter KE, Compagno LJ, Ebert DA, Gibson C, Heupel MR, Livingstone SR, Sanciangco JC. 2014. Extinction risk and conservation of the world's sharks and rays. eLife 3: e 00590. [CrossRef] [Google Scholar]
  • Dulvy NK, Metcalfe JD, Glanville J, Pawson MG, Reynolds JD. 2000. Fishery stability, local extinctions, and shifts in community structure in skates. Conserv Biol 14: 283–293. [CrossRef] [Google Scholar]
  • Dulvy NK, Pacoureau N, Rigby CL, Pollom RA, Jabado RW, Ebert DA, Finucci B, Pollock CM, Cheok J, Derrick DH, Herman KB, Sherman CS, VanderWright WJ, Lawson JM, Walls RHL, Carlson JK, Charvet P, Bineesh KK, Fernando D, Ralph GM. 2021. Overfishing drives over one-third of all sharks and rays toward a global extinction crisis. Curr Biol 31. [Google Scholar]
  • Ellis J. 2016. Raja clavata. IUCN Red List of Threatened Species. [Google Scholar]
  • Ellis JR, Pawson MG, Shackley SE. 1996. The comparative feeding ecology of six species of shark and four species of ray (Elasmobranchii) in the north-east atlantic. J Mar Biol Assoc UK 76: 89–106. [Google Scholar]
  • Espinoza M, Farrugia TJ, Lowe CG. 2011. Habitat use, movements and site fidelity of the gray smooth-hound shark (Mustelus californicus Gill 1863) in a newly restored southern California estuary. J Exp Mar Biol Ecol 401: 63–74. [CrossRef] [Google Scholar]
  • Elston C, Cowley PD, von Brandis RG, Lea J. 2022. Stingray habitat use is dynamically influenced by temperature and tides. Front Mar Sci 8. [CrossRef] [Google Scholar]
  • FAO. 2021. Fishery and Aquaculture Statistics . Global Capture Production 1950–2019 (FishstatJ). Quebec City, QC, Canada: 1; FAO. [Google Scholar]
  • Farias I, FigueiredoI, Moura T, Serrano Gordo L, Neves A, Serra-Pereira B. 2006. Diet comparison of four ray species (Raja clavata, Raja brachyura, Raja montaguiand, Leucoraja naevus) caught along the Portuguese continental shelf. Aquat Liv Resour 19: 105–114. [CrossRef] [EDP Sciences] [Google Scholar]
  • Figueiredo I, Maia C, Lagarto N, Serra-Pereira B. 2020. Bycatch estimation of Rajiformes in multispecies and multigear fisheries. Fish Res 232: 105727. [CrossRef] [Google Scholar]
  • Gilmour ME, Adams J, Block BA, Caselle JE, Friedlander AM, Game ET, Hazen EL, Holmes ND, Lafferty KD, Maxwell SM, McCauley DJ, Oleson EM, Pollock K, Shaffer SA, Wolff NH, Wegmann A. 2022. Evaluation of MPA designs that protect highly mobile megafauna now and under climate change scenarios. Glob Ecol Conserv 35: e02070. [CrossRef] [Google Scholar]
  • Hays GC, Bailey H, Bograd SJ, Bowen WD, Campagna C, Carmichael RH, Casale P, Chiaradia A, Costa DP, Cuevas E, Nico de Bruyn PJ, Dias MP, Duarte CM, Dunn DC, Dutton PH, Esteban N, Friedlaender A,Goetz KT, Godley BJ, Halpin PN. 2019. Translating marine animal tracking data into conservation policy and management. Trends Ecol Evolut 34: 459–473. [CrossRef] [Google Scholar]
  • Heupel M, Simpfendorfer C. 2014. Importance of environmental and biological drivers in the presence and space use of a reef‑associated shark. Mar Ecol Progr Ser 496: 47–57. [CrossRef] [Google Scholar]
  • Heupel MR, Simpfendorfer CA. 2002. Estimation of mortality of juvenile blacktip sharks, Carcharhinus limbatus, within a nursery area using telemetry data. Can J Fish Aquat Sci 59: 624–632. [CrossRef] [Google Scholar]
  • Hunter E, Buckley AA, Stewart C, Metcalfe JD, 2005. Repeated seasonal migration by a thornback ray in the southern north sea. J Mar Biol Assoc UK 85: 1199–1200. [CrossRef] [Google Scholar]
  • Hunter E, Buckley AA, Stewart C, Metcalfe JD. 2005. Migratory behaviour of the thornback ray, Raja clavata, in the southern north sea. J Mar Biol Assoc UK 85: 1095–1105. [CrossRef] [Google Scholar]
  • Hussey NE, Kessel ST, Aarestrup K, Cooke SJ, Cowley PD, Fisk AT, Harcourt RG, Holland KN, Iverson SJ, Kocik JF, Flemming JEM, Whoriskey FG. 2015. Aquatic animal telemetry: a panoramic window into the underwater world. Science 348. [CrossRef] [Google Scholar]
  • ICES. 2022. Working Group on Elasmobranch Fishes (WGEF). ICES Scientific Reports. Report. [Google Scholar]
  • IUCN Red List of Threatened Species [WWW Document]. 2023. IUCN Red List of Threatened Species. URL https://www. & searchType=species (accessed 4.14.22). [Google Scholar]
  • Joppa LN, O'Connor B, Visconti P, Smith C, Geldmann J, Hoffmann M, Watson JEM, Butchart SHM, Virah-Sawmy M, Halpern BS, Ahmed SE, Balmford A, Sutherland WJ, Harfoot M, Hilton-Taylor C, Foden W, Minin ED, Pagad S, Genovesi P, Hutton J. 2016. Filling in biodiversity threat gaps. Science 352: 416–418. [CrossRef] [PubMed] [Google Scholar]
  • Last P, Naylor G, Séret B, White W, Stehmann M, de Carvalho M. 2016. Rays of the World. CSIRO Publishing. [Google Scholar]
  • Lédée EJI. 2015. Improving understanding of spatial ecology through network analysis of acoustic monitoring data (PhD Thesis). James Cook University. [Google Scholar]
  • Leeb K, Villegas‐Ríos D, Mucientes G, Garci M, Gilcoto M, Alonso‐Fernández A. 2021. Drivers of spatial behaviour of the endangered undulate skate, Raja undulata. Aquat Conserv 31: 3466–3479. [CrossRef] [Google Scholar]
  • Lennox RJ, Aarestrup K, Cooke SJ, Cowley PD, Deng ZD, Fisk AT, Harcourt RG, Heupel M, Hinch SG, Holland KN, Hussey NE, Iverson SJ, Kessel ST, Kocik JF, Lucas MC, Flemming JM, Nguyen VM, Stokesbury MJW, Vagle S, VanderZwaag DL. 2017. Envisioning the future of aquatic animal tracking: technology, science, and application. BioScience 67: 884–896. [CrossRef] [Google Scholar]
  • Li Z, Wood SN. 2019. Faster model matrix crossproducts for large generalized linear models with discretized covariates. Stat Comput 30: 19–25. [Google Scholar]
  • Licandeo RR, Lamilla JG, Rubilar PG, Vega RM. 2006. Age, growth, and sexual maturity of the yellownose skate Dipturus chilensis in the south-eastern Pacific. J Fish Biol 68: 488–506. [CrossRef] [Google Scholar]
  • Lowerre-Barbieri SK, Kays R, Thorson JT, Wikelski M. 2019. The ocean's movescape: fisheries management in the bio-logging decade (2018- 2028). ICES J Mar Sci 76: 477–488. [CrossRef] [Google Scholar]
  • MacKeracher T, Diedrich A, Simpfendorfer CA. 2018. Sharks, rays and marine protected areas: a critical evaluation of current perspectives. Fish Fisheries 20: 255–267. [Google Scholar]
  • Maia C, Serra-Pereira B, Erzini K, Figueiredo I. 2015. How is the morphology of the oviducal gland and of the resulting egg capsule associated with the egg laying habitats of Rajidae species? Environ Biol Fishes 98: 2037–2048. [CrossRef] [Google Scholar]
  • Mcphie RP, Campana SE. 2009. Reproductive characteristics and population decline of four species of skate (Rajidae) off the eastern coast of Canada. J Fish Biol 75: 223–246. [CrossRef] [PubMed] [Google Scholar]
  • Meyer C, Holland K, Papastamatiou Y. 2007. Seasonal and diel movements of giant trevally Caranx ignobilis at remote Hawaiian atolls: implications for the design of Marine Protected Areas. Mar Ecol Progr Ser 333: 13–25. [CrossRef] [Google Scholar]
  • Morales JM, Moorcroft PR, Matthiopoulos J, Frair JL,Kie JG, Powell RA, Merrill EH, Haydon DT. 2010. Building the bridge between animal movement and population dynamics. Philos Trans Royal Soc B: Biol Sci 365: 2289–2301. [CrossRef] [PubMed] [Google Scholar]
  • Morfin M, Simon J, Morandeau F, Baulier L, Méhault S, Kopp D. 2019. Using acoustic telemetry to estimate post-release survival of undulate ray Raja undulata (Rajidae) in northeast Altantic. Ocean Coast Manag 178: 104848. [CrossRef] [Google Scholar]
  • Nathan R, Getz WM, Revilla E, Holyoak M, Kadmon R, Saltz D, Smouse PE. 2008. A movement ecology paradigm for unifying organismal movement research. Proc Natl Acad Sci 105: 19052–19059. [CrossRef] [PubMed] [Google Scholar]
  • Neat F, Pinto C, Burrett I, Cowie L, Travis J, Thorburn J, Gibb F, Wright PJ. 2014. Site fidelity, survival and conservation options for the threatened flapper skate (Dipturus cf. intermedia). Aquat Conserv 25: 6–20. [Google Scholar]
  • Papadopoulo K, Villegas-Ríos D, Mucientes G, Hillinger A, Alonso-Fernández A. 2023. Drivers of behaviour and spatial ecology of the small spotted catshark (Scyliorhinus canicula). Aquat Conserv 1–15. [Google Scholar]
  • Papastamatiou YP, Friedlander AM, Caselle JE, Lowe CG. 2010. Long-term movement patterns and trophic ecology of blacktip reef sharks (Carcharhinus melanopterus) at Palmyra Atoll. J Exp Mar Biol Ecol 386: 94–102. [CrossRef] [Google Scholar]
  • Pardo SA, Kindsvater HK, Reynolds JD, Dulvy NK. 2016. Maximum intrinsic rate of population increase in sharks, rays, and chimaeras: the importance of survival to maturity. Can J Fish Aquat Sci 73: 1159–1163. [Google Scholar]
  • Payne N, Gillanders B, Webber D, Semmens J. 2010. Interpreting diel activity patterns from acoustic telemetry: the need for controls. Mar Ecol Progr Ser 419: 295–301. [CrossRef] [Google Scholar]
  • Powell RA. 2000 Animal home ranges and territories and home range estimators. Research techniques in animal ecology: controversies and consequences. New York, NY: Columbia University Press, pp. 65−103. [Google Scholar]
  • Prontuario instalaciones [WWW Document]. 2022. URL [Google Scholar]
  • Queiroz N, Humphries NE, Couto A, Vedor M, da Costa I, Sequeira AMM, Mucientes G, Santos AM, Abascal FJ, Abercrombie DL, Abrantes K, Acuña-Marrero D, Afonso AS, Afonso P, Anders D, Araujo G, Arauz R, Bach P, Barnett A, Bernal D. 2019. Global spatial risk assessment of sharks under the footprint of fisheries. Nature 572: 461–466. [CrossRef] [PubMed] [Google Scholar]
  • Santos R, Medeiros-Leal W, Novoa-Pabon A, Crespo O, Pinho M. 2021. Biological Knowledge of Thornback Ray (Raja clavata) from the Azores: improving scientific information for the effectiveness of species-specific management measures. Biology 10: 676. [CrossRef] [PubMed] [Google Scholar]
  • Shaw AK. 2020. Causes and consequences of individual variation in animal movement. Movement Ecol 8: 12. [CrossRef] [Google Scholar]
  • Simpson SJ, Humphries NE, Sims DW. 2020. The spatial ecology of Rajidae from mark-recapture tagging and its implications for assessing fishery interactions and efficacy of Marine Protected Areas. Fish Res 228: 105569. [CrossRef] [Google Scholar]
  • Sousa I, Baeyaert J, Gonçalves JMS, Erzini K. 2019. Preliminary insights into the spatial ecology and movement patterns of a regionally critically endangered skate (Rostroraja alba) associated with a marine protected area. Mar Freshw Behav Physiol 52: 283–299. [CrossRef] [Google Scholar]
  • Stehmann M. 1995. A record of Raja clavata, the Eastern Atlantic Thornback Skate, from the southern Madagascar Ridge at Walters Shoal (Elasmobranchii, Rajidae). J Ichthyol 35: 63–73. [Google Scholar]
  • Thieurmel B, Elmarhraoui A. 2022. Package "suncalc" Title Compute Sun Position, Sunlight Phases, Moon Position and Lunar Phase. Available from: [accessed 16 April 2023]. [Google Scholar]
  • Trenkel VM, Charrier G, Lorance P, Bravington MV. 2022. Close-kin mark-recapture abundance estimation: practical insights and lessons learned. ICES J Mar Sci 79: 413–422. [CrossRef] [Google Scholar]
  • Villegas-Ríos D, Alós J, March D, Palmer M, Mucientes G, Saborido-Rey F. 2013. Home range and diel behavior of the ballan wrasse, Labrus bergylta, determined by acoustic telemetry. J Sea Res 80: 61–71. [CrossRef] [Google Scholar]
  • Villegas‐Ríos D, Freitas C, Moland E, Thorbjørnsen SH, Olsen EM. 2020. Inferring individual fate from aquatic acoustic telemetry data. Methods Ecol Evol 11: 1186–1198. [CrossRef] [Google Scholar]
  • Villegas-Ríos D, Réale D, Freitas C, Moland E, Olsen EM. 2017. Individual level consistency and correlations of fish spatial behaviour assessed from aquatic animal telemetry. Animal Behav 124: 83–94. [CrossRef] [Google Scholar]
  • Walker P, Howlett G, Millner R. 1997. Distribution, movement and stock structure of three ray species in the North Sea and eastern English Channel. ICES J Mar Sci 54: 797–808. [CrossRef] [Google Scholar]
  • Warwick-Evans V, Ratcliffe N, Lowther AD, Manco F, Ireland L, Clewlow HL, Trathan PN. 2018. Using habitat models for chinstrap penguins Pygoscelis antarcticato advise krill fisheries management during the penguin breeding season. Divers Distrib 24: 1756–1771. [CrossRef] [Google Scholar]
  • Whitehead PJP, Baucho ML, Hureau JC, Nielsen J, Tortonese E. 1986. Fishes of the North-eastern Atlantic and the Mediterranean. First edition, UNESCO. [Google Scholar]
  • Williams-Grove L, Szedlmayer S. 2016. Acoustic positioning and movement patterns of red snapper, Lutjanus campechanus, around artificial reefs in the northern Gulf of Mexico. Mar Ecol Progr Ser 553: 233–251. [CrossRef] [Google Scholar]
  • Wood SN, Goude Y, Shaw S. 2014. Generalized additive models for large data sets. J Royal Stat Soc: Ser C (Applied Statistics) 64: 139–155. [Google Scholar]
  • Wood SN, Li Z, Shaddick G, Augustin NH. 2017. Generalized additive models for Gigadata: modeling the U.K. black smoke network daily data. J Am Stat Assoc 112: 1199–1210. [CrossRef] [Google Scholar]
  • Xunta de Galicia, Conselleira de Medio Ambiente,Territorio y Vivienda. (27/12/2018). DECRETO 177/2018, de 27 de diciembre, por el que se aprueba el Plan rector de uso y gesti on del Parque Nacional MarítimoTerrestre de las Islas Atlánticas de Galicia. Diario Oficial de Galicia, Num. 34: 9778. Available at: [Accessed 23 May 2022]. [Google Scholar]

Cite this article as: Papadopoulo K, Villegas-Ríos D, Mucientes G, Hillinger A, Alonso-Fernández A. 2023. Drivers of the spatial behaviour of the threatened thornback skate (Raja clavata). Aquat. Living Resour. 36: 21

Supplementary Material

Table S1. Information of all Raja clavata tagged (including the dead individual used as control) displaying their biometric parameters and tagging information. Individuals were excluded from the analysis if they had less than three detection days.

Table S2. Behavioural metrics of the nine individuals of Raja clavata considered in the analysis. Individual fish ID, FD: number of filtered detections, DD: detected days, TD: total tracked days, RI: residence index, AS Tot: total activity space (km2), AS day: activity space (km2) during the day, AS night: activity space (km2) during the night, Excursion: number of times an individual left and returned to the study area,/: not enough data to estimate parameters.

Figure S1. Panel displaying: (a) ©Innovasea omnidirectional acoustic receiver attached to an auger anchor, screwed in the seabed (b) acoustic transmitter V13P-1x externally attached with double zip tie fixation system, (c) Raja clavata on measuring ruler marked with both T-bar tags (©Floy Tag) (left) and acoustic transmitter (right).

Figure S2. Time series of latitudinal and longitudinal position as well as distance travelled between consecutive centres of activity (see definition in the main text) and depth records of a discarded control dead individual of Raja clavata. The depth variation in the control individual corresponds to the tidal range

Figure S3. Workflow displaying the steps taken to obtain the different parameters of Raja clavata considered in this study.

Figure S4. Total, day and night activity space areas for all Raja clavata based on centre of activities. Red dots symbolise ©Innovasea acoustic telemetry receiver locations and the shaded area of the activity space (KUD = kernel utilization distribution).

Figure S5. Total, day and night activity space areas for all Raja clavata based on centre of activities. Red dots symbolise ©Innovasea acoustic telemetry receiver locations and the shaded area of the activity space (KUD = kernel utilization distribution).

Access here

All Tables

Table 1

Summary of the optimal generalized additive mixed-effects models investigating the (i) probability of presence and (ii) activity space of Raja clavata in the study area†.

All Figures

thumbnail Fig. 1

(a) Location of the study array in the Iberian Peninsula (red square). (b) Position of the study area (red shaded area) within the national park (green polygons). (c) Map of the study area showing the location of the ©Innovasea acoustic receivers (coloured dots), reference tags (red triangles) with the temperature data logger (green circles and triangles). (d) Detailed map of the study area displaying ©Innovasea acoustic receivers divided into four sectors, reference tags and bathymetry.

In the text
thumbnail Fig. 2

Abacus plot showing the daily presence of thornback skate (Raja clavata) in the study area. Days when an individual was present, are coloured in beige. The two black lozenges at the start and end of each time series represent respectively the tagging and end date of the study. Daily presence is displayed for the reference tag (red) and the dead skate (blue). Green lozenges represent death events: DESTAC-SPP-20-03 & DESTAC-SPP-20-04 on 19/10/2020.

In the text
thumbnail Fig. 3

Predicted probability of the presence of Raja clavata in the study area as a function of day of the year (a), sex (b) and sea bottom temperature (c). Black bars (b) and grey shaded areas (a, c) represent the 95% confidence interval. Black dots correspond to the raw data of probability of presence. Values used for predictions: sex = female, sea bottom temperature = 14 °C, day of the year = 260.

In the text
thumbnail Fig. 4

Overview of all sectors, showing the number of exits (a) and entry (b) routes taken during the study period by Raja clavata.

In the text
thumbnail Fig. 5

Predicted activity space of Raja clavata in the study area as a function of the week of the year. Grey-shaded areas represent the 95% confidence interval. Black dots correspond to the raw data of activity space. Values used for predictions: sex = female, sea bottom temperature = 14 °C.

In the text

Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.

Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.

Initial download of the metrics may take a while.