Burrow emergence rhythms of deep-water Mediterranean Norway lobsters (Nephrops norvegicus) revealed by acoustic telemetry

N. norvegicus supports one of the most commercially-important fisheries in the Mediterranean, and there is considerable interest in developing non-invasive sampling stock assessment methods. Currently, stock assessments are conducted via trawling or by UnderWater TeleVision (UWTV) surveys with limited capacity to provide direct population data due to the burrowing behavior of the species. Here, we used acoustic telemetry to characterize the burrow emergence and movement patterns of N. norvegicus in relation to internal tides and inertial currents in deep-sea habitats of the northwestern Mediterranean. We deployed acoustic tags on 25 adults between May and June 2019, and collected up to 4 months of detection data from each tag. Tagged lobsters spent approximately 12% of their time in door-keeping (i.e., detections at burrow entrance with no displacements) but differences existed among the different behavioral rhythms identified. We observed that mixed day-night, tidal and inertial rhythms in field settings were similar to those observed in laboratory tests. The presence of mixed day-night and tidal periodicity poses the question of why N. norvegicus exhibits flexible responses to environmental cues other than photoperiod. It is possible that tidal regulation of locomotor activity could reduce energy expenditure in relation to hydrodynamic drag. Inertial periodicity occurs in animals with no clear burrowing-oriented activity (always present on the seabed). Possibly, inertial-related movements are the result of a disruption of the biological clock. Our results are discussed in the context of how burrow emergence rhythms may bias UWTV surveys and how novel in situ monitoring approaches address these biases.


Introduction
Managing fishery resources in Marine Protected Areas (MPAs) requires the development of non-invasive sampling methods that, unlike trawling, do not damage habitats or remove biomass (Belgrano et al.Vol:.( 1234567890) 2021; Fernández-Chacón et al. 2021;Grinyó et al. 2022).For this reason, there is a general shift away from using trawl-based stock-assessment methods to those based on imaging (Aguzzi et al. 2022a).However, such methods will still remain limited in their capacity to provide accurate population data unless they also incorporate knowledge of how species' spatio-temporal movements can alter their detectability (e.g., Aguzzi et al. 2022a).One way to generate this data is to use acoustic telemetry devices to track the long-term movements and behavior of species of interest (Aspillaga et al. 2019;Matley et al. 2022).
The Norway lobster (Nephrops norvegicus) supports one of the most commercially-important fisheries in the EU (STECF 2020; Aguzzi et al. 2022b;Issifu et al. 2022).As such, managing its populations is of strategic importance, and it requires tools for the efficient and reliable stock assessments (Aguzzi et al. 2022a).Currently, stock assessments in the Atlantic and the Mediterranean are conducted either via trawling or by UnderWater TeleVision (UWTV) surveys, which are cameras mounted on towed sledges (Dobby et al. 2021).However, both approaches are unable to provide direct population counts of animals due to the burrowing behavior of N. norvegicus (Aguzzi et al. 2022a).As animals can only be captured in trawls or filmed by UWTV surveys when they are outside their burrows (Aguzzi et al. 2022a), the population-level patterns of burrow emergence have a massive impact on detection of animals (Bell et al. 2006;Aguzzi and Sardà 2008;Sardà and Aguzzi 2012).Typically, individuals emerge between nighttime and crepuscular hours on upper and lower shelf and at midday on slope (e.g., Chapman and Howard 1979;Oakley 1979;Aguzzi et al. 2003).UWTVs surveys address this issue by counting burrows instead of individuals under the assumption that one animal occupies each burrow complex (Dobby et al. 2021).However, burrows may be empty or occupied by dead individuals or other opportunistic species (e.g., Aguzzi et al. 2009a;Sardà and Aguzzi 2012).
While burrow emergence behavior appears largely conserved across populations (Aguzzi et al. 2022a), there is still a lack of high-resolution temporal data (i.e., at hourly frequency, over consecutive days) of this phenomenon.To date, most information comes from laboratory experiments on the locomotor activity rhythms that sustain burrow emergence.In controlled conditions, tests on north-western Mediterranean slope (400 m depth) animals showed displacement rhythms of mixed day-night, tidal-like, and inertial-currents oriented periodicity, plus a nonnegligible quantity of arrhythmic individuals (Aguzzi et al. 2004b, a;2011;Chiesa et al. 2010;Sbragalia et al. 2015).Also, animals in the same tank showed aggressive territoriality, with dominant individuals attempting to evict subordinates from their burrows (Sbragaglia et al. 2017).
Acoustic tracking technology can be deployed for long time into an MPA, staying safe from the threat of removal from trawling, being hence a promising tool for spatio-temporally monitoring of the burrowing behavior of N. norvegicus (Masmitja et al. 2020;Vigo et al. 2021).This method can reveal the reaction of individuals to light-intensity and water current speed cycles (when moored current sensors are also provided) over consecutive days, disclosing strategic information on species life-traits.Here, our objective was to characterize the burrow emergence activity in terms of a) detection and b) displaced distances in freely motile, acoustically-tracked N. norvegicus and their relationship with internal tides and inertial currents in a north-western Mediterranean deep-sea MPA.To achieve this, we revisited and reanalyzed acoustic detection data by Vigo et al. (2021) with an individual tracking focus, related to chronobiology, to describe burrow emergence behavior as a combination of stationary door-keeping and locomotor free epibenthic wandering.We aimed at relating the periodicity of this latter behavior (which has direct impacts on stock assessments) to the environmental modulation.In particular, we focused on the occurrence of peaks in day-night and inertial (current-related) periodicity in-situ, hence providing a comparative and more reliable dynamic for two decades of laboratory experiments (with animals collected in the same depth range and geographic area).

Fieldwork procedures
We used data from 25 acoustically-tagged N. norvegicus adults that were originally sampled as a part of the study by Vigo et al. (2021)   a north-western Mediterranean slope (315-475 m depth) MPA (Fig. 1).Creels were recovered at night, with deck operations conducted in dim-red light to avoid retinal damage to captured individuals (Shelton et al. 1985).Individuals were transferred in complete darkness to separate 20 × 30 cm cells in a compartmented aquarium with water at 13 °C (i.e., the temperature of the Mediterranean Sea below the 200 m thermocline; Hopkins et al. 1985).
Carapace Length (CL ± 0.05 mm) and sex (determined via the thickness of the pleopods) were recorded for all acoustically-tagged individuals (see Table 1).Acoustic transmitters were deployed just prior to their release the day after, the 8th of June 2019 (see Vigo et al. 2021 for additional details).Briefly, acoustic emitters (VEMCO V7-L1-69 k; Innovasea.18 mm length and a weight of 0.7 gr) were glued to the ventral side of the cephalothorax.Each The same analysis is reported for water column pressure, water salinity, and current velocity, as markers for tides and inertial currents.Total displacements in bold black represent the maximum and the minimum for the detected range of displacements per day, while bold black periods correspond to significant periods based on the Lomb-Scargle periodogram analysis transmitter spent a specific frequency signal (i.e., the ping) as the Identification Number (ID), every 120 s.Individuals were released at night at 1-2 m above the seabed, with a release-system consisting of a PVC receptacle with a dissolving latch closure on the door (see Vigo et al. 2021 for additional details).The release site was at the center of a monitoring area composed of 4 vertically moored hydrophones (VEMCO models: VR2W-69 k and VR2AR-69 k) (see Fig. 1).The 4 hydrophones were attached to their individual moorings at 20 m above the sea bottom and 150 m apart from each other, hence creating a detection radius of ~ 350 m.
During the tracking period, we measured several environmental parameters including: salinity, temperature, and pressure (CTD, SBE 37) on one of the hydrophone mooring lines (see Fig. 1).Pressure was used as a proxy for the tidal cycle; e.g., Aguzzi et al. 2020c).In addition to the CTD, a single point current meter (Aquadopp Nortek) was attached to the mooring line.

Detection and positioning of tagged individuals
The detection of each tagged individual (i.e., the presence/absence as a proxy for burrowing behavior) was reported when an acoustic ID signal was identified by at least one hydrophone.Conversely, non-detection indicated that either the animal was concealed in its burrow (with the sediment shielding the sound propagation toward the hydrophones) or that it had moved outside the detection range of the hydrophones.This method allows the high-frequency and continuous monitoring of animals' presence and displacement by day and night, over consecutive months (see below for the achieved date ranges), until tags and hydrophone batteries are depleted.The displacement of acoustic-tagged individuals (i.e., as a proxy for the expression of spatially-oriented locomotor activity) was computed by triangulating simultaneous detections by a minimum of three hydrophones.Each detected acoustic tag signal had a different timestamp (in terms of its position and the time of flight).Thus, the exact position of the sound emitting source was calculated by applying least squares regression to the time difference of the arrival of signal at the different hydrophones (see Masmitja et al. 2020, for additional details).To do so, we had to treat the internal clock-drift effect of each hydrophone (i.e., the same ping being received at different hydrophones at slightly different recorded time stamps, with the signal erroneously interpreted as a change in position).This was solved by attaching an additional acoustic tag to each of the 4 moorings (i.e., 1 m above each hydrophone) for synchronization.

Data processing
We obtained tracking data for just under 4 months: 8th of June at 23:00 to 23rd of September 2020 at 17:00.Time series for detections and displacement (as proxy for burrowing behavior) were compiled by summing 120 s readings by one hour.To compile the displacement time series, we discarded all the triangulated positions with velocities above a maximum threshold of 0.5 m/s, since they were biologically unrealistic (see Vigo et al. 2021 for details).For environmental readings, we extracted hourly averages from 30-min frequency readings.
Given the open nature of the monitoring area, with animals free to enter and exit the hydrophones' detection ranges, the detection, and displacement time series were interspersed with data gaps of variable length.Therefore, the longest time series segment, corresponding to each individual, was selected for further analysis.In order to do so, we established a length threshold of minimum three continuous days (i.e., the minimum reliable data segment length to statistically assess the occurrence of periodicity in behavior; Dunlap et al. 2003).
We considered a door-keeping behavior when individuals were detected and, at the same time, they presented less than 0 m or 0.5 m of displacement during one hour.Detection and movements time series were aligned at corresponding hours for each individual, having detection data for every 2 min.When the displacement was between 0 and 0.5 m and detections were present (no 0 s in every 2 min), the hours were counted.We could then estimate the relative time of door-keeping behavior, by calculating the percentage over the all the time series length.Data with 0 detections were also counted and estimated in relative percentages, to describe the duration individuals spent inside their burrow (no door-keeping behavior nor emergence).We also estimated the relative time of activity, when individuals were outside the burrow and presenting movement.Obtained temporal measurements were also related to each type of rhythm that the individuals presented in periodogram analysis (see below), to look for potential behavioral differences.
To avoid redundancy, the occurrence of significant periodicity was screened only in time series of displacement data.We used the Lomb-Scargle Periodogram of the "Lomb" package of the R statistical software (Ruf 1999; R Core Team).Periodicity was screened within 10 and 27-h intervals (i.e., accounting for day-night photic and tidal or inertial hydrographic cycles).Given the high variability in behavioral data usually reported in field studies, the classification of periodicity was as follows: significant (i.e., at α = 0.001) periodogram peaks between 22.5 and 24.4-h indicated day-night rhythms; peaks above 24.5-htidal rhythms; and finally, peaks between 16 and 20-h inertial periods (based on the latitude of the studied area; reviewed by Aguzzi et al. 2009b).
Secondly, a waveform analysis was conducted to assess the phase of detection and displacement rhythms.That analysis was performed at the 24-h (i.e., day-night) cycle scale, since fluctuations in photoperiod and natural light intensity largely modulate the burrow emergence behavior of N. norvegicus (Aguzzi et al. 2003;2009b;2015).Resulting waveforms for detections and displacements were superimposed to evidence the occurrence of door-keeping behavior (animals were detected at their tunnel entrance, but no displacement registered; Fig. 2).
For each individual, the selected detection and displacement time series were divided into 24-h sub-segments (Aguzzi et al. 2006).All values corresponding to each sub-segment (i.e., measurements at the same time of the day over the entire length of the selected period) were averaged.The resulting mean (± SE) values were plotted to determine the waveform and depict the mean fluctuation over a standard period of 24 h.The phase was then assessed according to the Midline Estimating Statistic of Rhythm (MESOR) method (Aguzzi et al. 2006).The MESOR value was estimated by re-averaging all waveform values, and the resulting value was represented as a horizontal threshold line onto the waveform itself.The Onset and Offset of waveform peaks were estimated by considering the first and the last value above the MESOR, respectively, indicating the peak's temporal limits.
Periodogram analyses were also applied to pressure, temperature, salinity, and water speed time series.Those time series showing significant periodicity at α = 0.05 were considered for further waveform analysis.
An integrated chart of waveform peaks (Aguzzi et al. 2012;2015) was used to visualize the temporal relationships of detection and displacement rhythms among all individuals and in relation to the hydrographic and day-night cycles.The duration of the waveforms' peaks as delimited by Onset and Offset values (see above) is presented in the form of continuous horizontal bars depicting the peak temporal amplitude for an individual.The temporal amplitudes for all individuals were represented together and in association with the results of the same waveform analysis for oceanographic data (including the day hours as proxy for the day-night cycle).

Comparisons of rhythmic behavior by sex and body size
Individuals were assigned one of four behavioral levels based on their dominant displayed periodicity: "arrhythmic", "day-night", "inertial" and "other" (which included short semi-diurnal and long tidal periods).The counts of male and female individuals displaying each behavior were treated as Poisson events (O'Hara and Kotze 2010) described by a discrete probability distribution, with each combination of sex and behavioral label having an expected rate λ (non-negative real number).
We selected a Bayesian approach instead of conventional, frequency-based analysis, as it provides a full distribution instead of a single value for the estimations of the λ parameter, after reallocating the credibility of probabilities across all possible values.In addition, the credible intervals which express the uncertainty do not depend on a large sample size (i.e., are N-independent).Finally, it allows quantifying the support in favor of the null hypothesis, instead of only against it (Kruschke 2021).
For count ratio comparisons (i.e., if the sex ratio of the entire sample was reflected in each behavioral category), we used a Bayesian analogue (i.e., Bayes Factor) of the classic contingency table analysis.This method tests the independence assumption in a contingency table under a Poisson sampling plan without a fixed total size N (Gunel and Dickey 1974).The analysis was performed in the R package "Bayes-Factor" (Morey and Rouder 2022).Then, probability density distributions for each λ were estimated from 10 4 posterior draws which were generated by a Markov Chain Monte Carlo (MCMC; McGrayne 2011) sampling process.This procedure simulates a parameter in sequential steps, and the value of step i depends only on the value at step i-1 (i.e., the previous state, independently of the initial step 0 and the process that led to i-1).The distributions of all λ rates were visualized in density plots alongside their 90% credible intervals and real observed counts.
Finally, the potential effect of CL on rhythmic behavior was analyzed using box-plots and comparing mean, median and variation of observed CL per behavioral label.

Detection and displacement time series
Detection and displacement time series for all animals are reported in Figure S1 and S2, with a representative section per individual considered for quantitative and time series analyses (see below).
Detection rates varied among the 25 tracked individuals, resulting in time series of different lengths (see Table 1): 8 animals ≤ 1 week; 5 animals ≤ 1 month; 2 animals ≤ 2 months; and finally, 10 individuals < 2 months (of those latter, 4 were present over the whole period of recording).The average displacement distance was 216.0 m/day (min: 90.7 m/day, max: 424.7 m/day).
The individuals showed a relative time of 12.92% and 11.81%, respectively, presenting a door-keeping behavior identified.The mean time spent inside the burrows, when detections and movement were 0, represented 38.31% of their time series, while the mean relative time in activity was about 49.88%.The values for each individual are shown in Table 2.The mean values, estimated for each type of rhythm (Table 3), indicated that the day-night periodicity was the rhythm in which individuals presented less activity outside the burrows (35.94%), with a similar doorkeeping behavior between the inertial individuals (16.18 and 15.25% respectively).
Periodogram analyses for displacement time series evidenced the occurrence of 4 major groups of animals, showing locomotor patterns matching frequency ranges within day-night, tidal, and inertial ranges, or no significant periodic activity (see Table 1).Assuming that day-night periodicities are in the range of 22.5 and 24.4-h, 12 animals (48%) followed that rhythmic pattern with a major activity phase at daylight hours (see Table 1).An example is described by the waveform analysis output for an individual (ID 14468) as representative of this group, in Fig. 2A.That animal shows a significant increment in displacement (i.e., an increase above the MESOR) from 6:00 to 19:00.Superimposing detection on the displacement waveform (see Fig. 2A), moments of stationary door-keeping behavior (i.e., animals waiting at the entrance of their tunnels) were revealed, when displacements were greatly reduced.This is evident for that individual, which showed a baseline level of detections around 19 per hour, when displacements were close to 0 (i.e., approximately from 21:00 to 04:00).
Detections also increased broadly at daytime, being present also during the night, when displacements dropped to zero from 22:00 to 04:00 (again, for the occurrence of stationary door-keeping behavior).
3 animals (12%) displayed inertial (i.e., periodicity between 16 and 20-h) displacement rhythmicity (see Table 1).An example is reported in waveform analysis for an individual (ID 14475) within this category in Fig. 2C.No temporally-structured increment of displacement activity appears in relation to the light cycle.The temporal trend of detections follows a similar pattern.Since the waveforms for both data sets are similar in that animal and the others for this same group, it appears that door-keeping (and hence burrowing) behavior does not occur at all.

Environmental time series
Time series for environmental data are reported in Figure S3.Water temperature fluctuated between a minimum and a maximum of 13.66 and 13.92 °C, respectively (with an averaged value and standard deviation of 13.82 ± 0.04 °C).Salinity varied between 38.55 and 38.65 as maximum and minimum (with an averaged value of 38.63 ± 0.01).Current speed occurred in the minimum and maximum range of 0.01 and 0.21 m/s (with an averaged value of 0.06 ± 0.03 m/s).Finally, the water column pressure fluctuated between 357.45 and 358.03 dbar (with an averaged value of 357.70 ± 0.10 dbar).
Periodogram analysis indicated the presence of significant periodicity for pressure, salinity, and water speed (see Table 1 for period values and Figure S4 for periodogram plots).Pressure showed a mixed semi-diurnal / diurnal periodicity pattern, with the other three occurring predominantly within the inertial range.Temperature did not show any significant periodicity and was therefore discarded from waveform analysis.
Waveform analysis for current speed and water pressure showed phase opposition (Fig. 3).The waveform for speed indicated the occurrence of a significant peak (mean values above the MESOR) from 23:00 to 10:00.Conversely, a significant increase in the fluctuation status for the water column pressure occurs from 17:00 to 01:00.
Table 2 The relative duration of "door-keeping" behavior for each acoustic tagged (ID) individual is represented as a relative percentage of time (from the total time of recording; see Table 1 for data sets lengths in days) in two ways: ≤ 0.5 (m), where individuals had detections and movement of less than 0.5 m in an hour; and 0 (m), where individuals had detections and 0

movement in an hour
There is also the proportion (in percentage) of individuals within their burrows ("in burrow") when 0 detections and 0 movements in an hour were met (0(Det.)).The activity outside the burrow ("full emergence") is represented by the remaining time series percentages values for the two methodologies mentioned before, considering full emergence when movement is over 0.5 m (> 0.5 (m)), and when it is over 0 m.This refers to the time in which the animal was moving further away from the burrow entrance

Comparisons by sex and body size
The Bayes Factor comparisons for differences in counts per behavioral category between male and female N. norvegicus (i.e., in relation to the expected ones based on the initial sample, roughly a ratio of 2:1) did not reveal a strongly pronounced deviation from the null hypothesis of independence, overall.In particular, the outcome of the analysis was a Bayes Factor of 1.33.
In order to get more insight of the pairwise comparisons within each behavioral category, the posterior distributions of the corresponding Poisson rates λ were derived and visualized.Most categories showed the same pattern except for the few individuals with inertial-based periodic activity, where more females were encountered (see Table 1) and there is a marked overlap of the distributions of λ for male and female (Fig. 4).In regard to the CL comparisons, the box-plots (Figure S5) showed very similar characteristics.

Integrated chart of waveforms peaks for biological and oceanographic data
The temporal phase relationship, obtained by comparing peaks of all waveforms (i.e., values above the MESOR) for detections and displacement time series, is presented in Fig. 5.For displacement time series, significant increases (constituting waveform peaks) are concentrated at daytime (i.e., from 6:00 to 19:00) with a marked reduction at night.Detections follow the same pattern at daytime, but significant waveform peak values occur also at nighttime (bimodal).Nocturnal detections without displacement indicate the occurrence of door-keeping behavior as a proxy for burrowing activity (i.e., because animals are detected but do not leave the entrance).That overall temporal dynamic highlights those differences reported at the level of individuals' waveforms (see Fig. 2).
The comparison of both behavioral time series with waveform peaks for current velocity and water column pressure (see Fig. 3) suggests under which Table 3 The mean percentage (change for ± SD) of individuals presenting door-keeping behavior for each type of rhythm described in Table 1 N indicates the number of individuals classified in each type of rhythm.This door-keeping behavior is represented in two ways: < 0.5 (m), when individuals presented detections but were moving less than 0.5 m in an hour; 0 (m), when individuals presented detections but no movement.We also indicated the mean percentage when individuals did not present movement nor detections (0 (D)) and when they were active, presenting moving and detections Fig. 5 The integrated analysis of phases depicts peaks (i.e., waveform values above the MESOR) in displacement activity (black), detections (red) and in the fluctuation of measured oceanographic variables (blue) Vol:. ( 1234567890) levels of hydrodynamic fluctuations burrow emergence behavior may be inhibited.If we consider the time interval 06:00-19:00 (see Fig. 5), as the time window where the majority of displacement peaks occurred, animals globally increased their displacement activity when current speed starts decreasing (between around 0.064 -0.058 m/s).The general daytime activity also coincides with the midday minimum in water column pressure at 12.00-13:00.That activity is reduced when pressure drops below 357.70 dbar (06:00) and rises above 357.74dbar (19:00).

Discussion
Long-lasting and high-frequency acoustic-based tracking of N. norvegicus in deep-sea benthic habitats allowed us to study its burrow emergence rhythms over a long period of time.These observations provide new information on how this species' biological clock in relation to burrow occupation and doorkeeping (Atkinson and Naylor 1976;Arechiga and Atkinson 1980;Moller and Naylor 1980;Chiesa et al. 2010;Aguzzi et al. 2011), as well as emergence and non-burrow occupation behavior, that might affect stock assessments (Aguzzi et al. 2022a).Simultaneously, these results also explain at individual level, the historically-reported population shelf and slope catchability patterns on a day-night and tidal basis (Oakley 1979;Bell et al. 2006;Aguzzi and Sardà 2008;Sardà and Aguzzi 2012).Tagged lobsters did not show long periods of concealment into their burrows (i.e., no long-lasting absence of detections in Figure S1 are evident), being instead engaged in enduring or repeating door-keeping activities (i.e., detections with no displacements).The lobsters patrol the seabed around their tunnel entrances (Aguzzi et al. 2006) as has been recently confirmed in the same MPA study area by ROV surveys carried out over the 24-h (Vigo et al. 2023).Here, variable ranges of displacement were detected over consecutive days (i.e., see the MESOR estimate as comparison among displacement time series of variable length in days).These observations support the idea of marked inter-individual variability in locomotor activity rhythms as observed in laboratory experiments (e.g., Aguzzi et al. 2004b, a;Katho et al. 2013;Sbragalia et al. 2013Sbragalia et al. , 2014;;García et al. 2019).
We observed that the proportion of individuals conducting mixed day-night, tidal and inertial rhythms in field settings were similar to those observed in laboratory tests (Aguzzi et al. 2013).Laboratory experiments also reveal the occurrence of mixed periodicity that points out the plasticity of the species' biological clock; i.e., able to adapt to a wide range of photoperiod lengths (for latitude and seasons) as well as intensities (for depth) (Aguzzi et al. 2011).Those rhythms are occurring in the Mediterranean disphotic slope environment, where even though the monochromatic blue-green light (i.e., wavelength between 400-700 nm) is dampened by several orders of magnitude, it still elicits a turnover of the community on a 24-h basis (Aguzzi et al. 2003;2015).Therefore, while the ecological basis of the daynight rhythmic expression can be more apparent, the interpretation of tidal and inertial fluctuations is less straightforward.A large proportion of animals were diurnal, while one individual displayed a bimodal (tidal-like) periodicity, which at the surveyed depth seems to be similar to crepuscular rhythms of shallower shelf populations.That pattern in displacement could be recognized in the periodogram analysis as tidally-oriented, being instead of borderline periodicity with the day-night one.Finally, inertial periodicity has been detected as per the laboratory results and field hydrographic comparisons by Aguzzi et al. (2009a), but the ecological reasons remain harder to explain.
The presence of mixed day-night and tidal periodicity poses the question of why animals exhibit flexible responses to different environmental cues other than photoperiod.If N. norvegicus followed an optimal interval in environmental light intensity to perform emergence (Chapman et al. 1975) as a visual predator and scavenger (e.g., Farmer 1975;Oakley 1979) or to avoid incoming predators (Richardson 1996;Chapman et al. 2000;Aguzzi et al. 2009b;Albalat et al. 2016), it should be less sensitive to currents.Alternatively, our data indicate that emergence behavior is influenced by tidal patterns (Newland et al. 1988).Weak tidal patterns influence sediment resuspension and transport in deep-water and deep-sea areas of the western Mediterranean Sea (Palanques et al. 2002;Grinyó et al. 2017), likely affecting the phenology of sessile megafauna (Grinyó et al. 2018).A tidal pattern in the species catchability was evidenced in the North Sea (Bell et al. 2008) and may be the product of animals' locomotor activation or inhibition upon predators or preys' odor cues within water flows.Odor cues play a relevant role in the feeding activity and success of deep-sea fauna (Aguzzi et al. 2012;Stuckless et al. 2021) and may contribute to activity rhythm regulation (e.g., Aguzzi et al. 2018).At the same time, a tidal regulation of locomotor activity also occurs in laboratory hydrodynamically controlled environments with no chemical cues (Sbragaglia et al. 2013).This indicates that animals may vary behavioral activity to compensate for energy expenditures in relation to hydrodynamic drag (e.g., Aguzzi et al. 2010).As observed in other deeper dwelling species such as the deep-sea crab Geryon longipes from the same area at 1000 m depth, locomotor patterns can be synchronized to both tidal and inertial currents in the laboratory (Nuñez et al. 2016).
Contrary to the tidal periodicity of the water column depth, current speed and salinity showed inertial rhythms, while temperature was arrhythmic.In a previous analysis comparing locomotor activity in laboratory tested animals of the same depth range and geographic location by Aguzzi et al. (2009a), inertial rhythms were explained based on near seafloor mooring data obtained at 1,000 m depth.Here, with our shallower mooring deployment (380 m depth), we observed a mixture of inertial currents and a tidal signal.However, water speed ranges were different.In Aguzzi et al. (2009b), water mean speed was approximately 0.035 m/s during 30 days (from 15th September to 14th October), while here from the 8th of June to 23rd of September, our averaged values are of 0.06 ± 0.03 m/s.In any case, we still have inertial related periodicity in the displacement behavior of some animals that show no clear burrowing-oriented activity.Possibly, rather than being the sign of synchronization of behavioral rhythms upon this hydrodynamic cue, inertial periodicity in displacements could be the results of a disruption of the biological clock according to a population of coupled oscillators (Strogatz & Stewart 1993;Aguzzi et al. 2011).For unknown conditions, some animals clock oscillators would split in different functional subgroups, hence producing, inertial-like rhythms, equivalent to ultradian periodicity in laboratory studies (Aguzzi et al. 2011).
Our data on water properties are in line with the generally accepted preferences of N. norvegicus (i.e., highly saline waters with small temperature fluctuations; Eriksson et al. 2013).However, there is no sufficient literature describing the influence of salinity on this species' rhythmic behavior (apart from survival studies showing higher mortality under low salinity regimes).Differently, a negative impact of high salinity on CPUE has been recently identified (Chiarini et al. 2022).It is possible that animals react to salinity fluctuations as an additional chemical cue indicative of the change of tides and hydrodynamic activity, as it has been observed in laboratory conditions of rhythmic behavioral entrainment in other benthic decapods (Naylor 1985).
Sex and size (i.e., CL) did significantly affect burrow emergence rhythms of tracked individuals.The value of the Bayes Factor comparing the sex ratio contingency table was within the "not worthy of mentioning" category (i.e., for the range of 1 to 3.2), according to Kass and Raftery (1995).A potential exception could be the case of inertial rhythms, for which the observed sex-ratio deviates from the observed 1:2 (females:males) of all other behavioral categories and across the entire sample.Similarly, the median CL of individuals displaying inertial periodicity is lower.Females are generally reported to be smaller than males (McQuaid et al. 2006), which was also the case in our study.In any case, these results should be treated with caution, as the number of individuals is low (i.e., N = 3), and should not be considered conclusive evidence suggesting dimorphism in behavioral rhythmic response.Reproductive behavior occurring in spring-summer at our latitude (Aguzzi et al. 2004b, a) could affect reproductively mature individuals of both sexes in different ways, but more field observations are needed.
The cease of the detection of animals within the monitored area could be the product of their dispersion outside the detection range area.As the MPA gets progressively depleted of individuals, and the releasing areas experience a transitory high density, animals seem to adjust their territoriality with different ranges of wandering (see Vigo et al. 2021 for a spatial representation of those data).This could occur for a current process of individuals' redistribution as caused by aggressive hierarchic social interactions (Sbragaglia et al. 2013).In fact, animal burrowing capability is influenced by a balance between aggressive territorial interactions and overall local demography (Chapman and Rice 1971;Johnson et al. 2013).Possibly, individuals with no defined detection and Vol:.( 1234567890) wandering rhythmic patterns are those with no burrowing capability due to their subordinate social ranking level.Future studies are needed, to monitor this enigmatic component of the population; e.g., acoustic tracking should be accompanied by a network of video-recording cameras from fixed and docked mobile platforms, branching off the MPA central release area (Aguzzi et al. 2022a).HD images could be used to identify those displacing animals, evidencing the reasons for their burrow digging incapability.

Conclusion
The use of acoustic telemetry to characterize the movements and burrowing behavior of N. norvegicus shows promising applications for ecological monitoring and recovery.Firstly, it allows the validation of UWTV and trawling surveys with an individual tracking technology, indicating that the same animals do emerge constantly over consecutive days (at least for those time series sections of continuous data, we analyzed).In animals with no burrows (for unknown reasons impairing the digging capability), behavioral activity can be measured under the form of a constant wandering on the seabed at a periodicity that matches inertial currents.So, when the burrow is present, it allows the temporal structuring of behavior into different stationary door-keeping locomotor wandering activities according to light cycles.Those observations indicate that trawl or video sampling for stock assessment purposes by ROV and towed sledges should not only count burrows but also animals both in door-keeping and in full emergence over consecutive 24-h cycles, in the same areas.Even if the bias (i.e., under-or overestimation of true abundance) can assumed to be constant, knowing the timing of maximum emergence of populations may assist in setting the best timings for video surveys, where both animals and tunnels should be counted.With repeated, temporally-oriented surveys, the difference between counted burrow systems and animals can be estimated in order to obtain a proportion coefficient to adjust the derived stock estimates.Finally, we need to extend the number of moored hydrophones in order to enlarge the MPA monitored area, to get an insight on animals' longer range displacements.In the end, we need to better assess if captured and reintroduced animals may contribute to the local stock refurbishing into the MPA or individuals wander across its borders, vanishing the conservation strategy.
Vol.: (0123456789) Commons licence, and indicate if changes were made.The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material.If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder.To view a copy of this licence, visit http:// creat iveco mmons.org/ licen ses/ by/4.0/.

Fig. 1
Fig. 1 Study area in the Catalan margin (NW Mediterranean, Spain) A on the slope between two large deep-sea submarine canyons, Cap de Creus and the Fonera-Palamós.Spatial distribution of N. norvegicus is indicated as Landings Per Unit Effort (LPUE) along the Catalan coast, obtained by combining vessel monitoring system information and official landing data in the annual average data from period 2008 − 2021.The study area is represented in B indicating the location of the telemetry

Fig. 2
Fig. 2 Waveform analysis outputs (mean values ± SE) on time series of detections (red) and displacements (black) for a group of 4 representative individuals per each group of activity rhythms: A (ID: 14,468) day-night; B (ID: 14,472) tidal; and finally, C (ID: 14,475) inertial.MESORs are the dashed horizontal lines reported only for those locomotor patterns with evident activity increase over the 24-h (see Table 1).The black arrows indicate the start and end of activity (i.e., both detection -individual outside the burrow-and displacement -individual ).Briefly, animals were captured between May 30th to June 7th, 2019 by creels on the muddy bottoms of Vol.: (0123456789)

Table 1
Tagged (ID) individuals with sex (M-Male, F-Female), size (Carapace Length; CL), and time series (TS) length in days, and walked distance (Total and per Day) speci-fications, plus the results of waveform and periodogram analyses (MESORs plus periods in minutes-min and hours-h)