Search tips
Search criteria 


Logo of scirepAboutEditorial BoardFor AuthorsScientific Reports
Sci Rep. 2017; 7: 44938.
Published online 2017 March 23. doi:  10.1038/srep44938
PMCID: PMC5362963

Inertial bioluminescence rhythms at the Capo Passero (KM3NeT-Italia) site, Central Mediterranean Sea

J. Aguzzi,a,1 E. Fanelli,2 T. Ciuffardi,2 A. Schirone,2 J. Craig,3 and KM3NeT-Italia/NEMO Collaboration


In the deep sea, the sense of time is dependent on geophysical fluctuations, such as internal tides and atmospheric-related inertial currents, rather than day-night rhythms. Deep-sea neutrino telescopes instrumented with light detecting Photo-Multiplier Tubes (PMT) can be used to describe the synchronization of bioluminescent activity of abyssopelagic organisms with hydrodynamic cycles. PMT readings at 8 different depths (from 3069 to 3349 m) of the NEMO Phase 2 prototype, deployed offshore Capo Passero (Sicily) at the KM3NeT-Italia site, were used to characterize rhythmic bioluminescence patterns in June 2013, in response to water mass movements. We found a significant (p < 0.05) 20.5 h periodicity in the bioluminescence signal, corresponding to inertial fluctuations. Waveform and Fourier analyses of PMT data and tower orientation were carried out to identify phases (i.e. the timing of peaks) by subdividing time series on the length of detected inertial periodicity. A phase overlap between rhythms and cycles suggests a mechanical stimulation of bioluminescence, as organisms carried by currents collide with the telescope infrastructure, resulting in the emission of light. A bathymetric shift in PMT phases indicated that organisms travelled in discontinuous deep-sea undular vortices consisting of chains of inertially pulsating mesoscale cyclones/anticyclones, which to date remain poorly known.

Light as an evolutionary force shapes marine ecosystem functioning, regulating the biological clock of predator and prey species, which react by performing rhythmic displacements into photic and disphotic realms1. However, beyond the reach of solar light, where virtually all ecologically-functional radiance is of biogenic origin2, the sense of time is dependent upon other geophysical fluctuations such as internal tides and inertial currents3,4. The deep sea forms the largest biome on Earth, yet the role of these forces on the temporal regulation of water-column and benthic communities’ activity rhythms is to date almost unknown5. The behavioural and physiological synchronization of organisms to flow changes6 is likely, as currents can form a key selective agent in relation to important ecological functions, such as low-energy dispersal and feeding7.

Data on the modulation imposed by cyclic flow changes on the biology of deep-sea organisms, although sparse, suggest the synchronization of behaviour and underlying physiology at internal-tidal periodicities (e.g. refs 6 and 7. However, evidence for the synchronization of biorhythms to the weaker atmospheric-related inertial currents is more scant; e.g. laboratory tests with the deep-water Norway lobster (Nephrops norvegicus) and deep-sea crab Geryon longipes8,9,10. Inertial currents are generated by wind drift on the water surface producing downward propagating whirlpools, perceived at the seabed as hydrodynamic oscillations. Their periodicity depends on the Coriolis force and therefore on latitude11. In the Mediterranean Sea, where tidal forces are very small (except for the North Adriatic)12, inertial currents dominate over internal tidal current periodicities, therefore providing a natural laboratory to identify the importance of these motions on marine communities.

Bioluminescence is widespread in the deep sea where it is found in most of the major marine phyla, from bacteria to fish. Luminous signals are an important aspect of deep-sea pelagic and benthic life, where they are used for a wide range of inter-and intraspecific interactions, permitting organisms to find or avoid each other for purposes of defense, offence, and reproduction2,4,13. In undisturbed water, the occurrence of bioluminescent flashes is reported to be rare14,15, presumably because of the energetic cost of light production16. However, defensive bioluminescent responses are provoked when organisms, advected by water currents, collide with submerged objects such as underwater vehicles or complex topography, e.g. rocky outcrops and their protruding sessile fauna14,17,18,19. Presently, little is known of the ecological relevance of temporal patterning in the bioluminescence emitted by deep pelagic fauna that results from the interaction of animals with cyclic water flow changes.

Neutrino telescopes deployed in the deep sea, instrumented with light detecting photomultipliers (PMTs) as well as oceanographic and acoustic sensors, are providing data beyond their primary physical purpose and are being used to describe biological phenomena in poorly monitored deep-sea areas. Datasets generated from this technology are enabling new discoveries about life in the abyssopelagic environment20,21,22,23 and PMT data are providing evidence of deep blooms in bioluminescence coupled to seasonal changes in passive carbon sinking, linking atmospheric-climatic and demersal components24,25,26.

The Capo Passero site, selected for the installation of the KM3NeT-Italy Neutrino Telescope (http://www.km3net.org27 presents appropriate properties for the study of bioluminescence by abyssopelagic organisms and their interaction with water flows of inertial periodicity from the bottom up to few hundreds of meters into the water column. The area is highly oligotrophic28,29 which results in very low levels of bioluminescent bacteria and bio-fouling on deployed structures at depths greater than 2500 m30,31. The location of this detector enables data collection on deep layer dynamics of the Eastern basin, which plays an important role in the whole Mediterranean circulation. In this basin, water mass formation and stratification has shown large variability in the last few decades, following the occurrence of the well-known Eastern Mediterranean Transient32,33,34. A recent study by Rubino el al.35 has shown that abyssal quasi-inertial vortices, with mesoscale feature, exist at this site, below 2500 m.

Here, we used high-frequency data retrieved from a prototype detection unit for an underwater neutrino telescope, with two different objectives. An ecological objective was to characterise rhythms in bioluminescence at the different depths of the telescope structure, as a proxy of the benthopelagic movement of organisms within the water column, relative to the seabed. Further, we observed how deep-sea mesoscale vortices may affect bioluminescence rhythms at different depths, generated by organisms in response to the flow of water masses.


We obtained a total of 32,970 bioluminescence readings of 10 min at the 8 telescope floors. Time series for these readings have been reported at all floors in Fig. 1 (for ease of visualization the time stamp is 24-h). Marked rhythmic fluctuation patterns were visible at all floors, although overall bioluminescent activity appeared to increase in the second half of the testing period. Few and short gaps in data acquisition were reported, with sparse outliers. These spike values were neither eliminated nor smoothed by filtering (e.g. moving average) as they were considered true biological signals not produced by instrument malfunctioning. Background light emitted by 40K decay was around 50 kHz, except for floors no. 1 and 5 where a minimum of 45 kHz was measured.

Figure 1
Time series of PMT readings for the 8 sampling floors of the telescope tower, as recorded during the entire testing period of June 2013.

The intensity of bioluminescence varied through floors as a function of depth (Fig. 2). Overall, a large majority of bursts (i.e. photons above a certain threshold that trigger a reading from the PMT) were reported within the 50–100 kHz range and a minority above 100 kHz. In particular, the presence of bursts within the 100–150 and >150 kHz ranges increased moving from floors close to the seabed, up within the water column.

Figure 2
Percentages of bioluminescence bursts for each defined intensity range at each floor depth: white 50–100 (white); 100–150 (grey); and >150 kHz (black).

A preliminary inspection of PMT time series (see Fig. 1) revealed differences in overall intensity from the first to the second half of the month, in agreement with changes in overall hydrodynamic regime (suggested by the increasing amplitude of floor displacement shown for floor no. 5). Accordingly, time series analysis was carried out by subdividing all biological and arm movement datasets into two halves of 13 days each, 3rd–16th and 17th–30th (16 and 17 full cycles of 1230 minutes, respectively); these were selected to have the maximum number of cycles for each sub-period while maintaining statistical representativeness. Evidence of the increase in burst intensity, from the 1st to the 2nd testing period, is further shown in Table 1, and confirmed by the reported amplitudes resulting from the Fourier analysis.

Table 1
PMTs bioluminescence reading performance per floor and time series analysis outputs for the 2 consecutive periods in June 2013 (3rd–16th and 17th–30th).

We detected a significant periodicity in bioluminescent intensity, corresponding to the rhythms of inertial currents at the latitude of the telescope site at all floors during the 1st and the 2nd testing period (periods between 19.1 and 20.8-h; Table 1). Rhythm strength as amplitude of the oscillation was generally low, with periodogram peak variance less than 20%. However, that periodicity progressively got stronger (as indicated by the overall increase in periodogram variance; see mean values in Table 1) when: i. moving from the first to the second half of the testing period; and ii. when moving away from the seafloor. In relation to this second finding, periodogram variance increased moving from deeper (i.e. lower) to shallower (i.e. higher) floors, as an indication of a depth-trend not only in intensity but also in the strength of the reported temporal patterning.

In our results, we simultaneously assessed the composite changes in bioluminescence intensity both over time (at each floor) and depth (through floors). Bioluminescence rhythms appeared reduced in peak amplitude (hence periodicity strength) in the first monitoring period compared to the second monitoring period. In addition to this, the rhythms appeared stronger at the shallower floors. Clearly, these trends were combined and the results should be interpreted cautiously. Outputs from the waveform and Fourier analyses (Table 1 and Fig. 3) confirmed the periodogram outputs, evidencing the presence of mean patterns of fluctuations in bioluminescence at each telescope floor, with lower peak amplitudes in the first period (Fig. 3A) compared to the second one (Fig. 3B). In both monitoring periods, a dampening of bioluminescence occurred through the floors, becoming more evident in the second period.

Figure 3
Waveform analysis outputs for data sets in bioluminescence for the 1st (3–16th) and 2nd (17–30th) period (A and B, respectively) of testing in June 2013, indicating the occurrence of temporally coherent peaking over the inertial day (i.e. ...

To assess the likelihood of stimulation of bioluminescence by currents, outputs of biological and tower orientation waveform analysis were plotted together with Fourier components at the defined period (i.e. 1230 minutes), for floor no. 7 (Fig. 4). To perform a comparison, the rate of floor displacement, calculated using the floor orientation, was assumed to be a proxy for current speed. The corresponding continuous lines, derived from the Fourier analyses, basically reproduce the results of the waveform analysis: the main difference appears for the 1st testing period, where the Fourier component of the rate of displacement (upward right panel in Fig. 4) shows a periodicity that did not appear in the waveform analysis. Furthermore, in the 1st testing period, the PMT output appears in phase with the floor orientation, whereas in the 2nd testing period it was the rate of displacement that noticeably overlapped with the bioluminescence signal. With respect to the angular data, it should be noted that the rotation of the structure is expected to be minimal as a result of its size, in line with the angular values shown in Fig. 4 (i.e. less than 5 degrees on the 1st testing period and about 10 degrees on the 2nd).

Figure 4
Waveform analysis output coupling of bioluminescence (black) phase (i.e. peak timing) with floor (grey) orientation (left) and rate of displacement (right).

Figure 5 shows the depth-oriented time shift in maximum bioluminescence (i.e. phase amplitudes, as computed by waveform analysis on PMT data at each floor; see Fig. 3) for the two testing periods. Peak onsets and offsets show a progressive depth delay along the tower axis, also confirmed by the phase values derived from the Fourier analysis (Table 1). The PMT antiphase peak is greater at shallower compared to deeper floors (i.e. maximum uncoupling at floor no. 8 and maximum phase synchronicity at floors no. 3 and 2). No similar conclusions can be made for the first testing period, where the rate of displacement shows a scattered distribution of PMT peaks (i.e. waveform phases) through the depth layers. Additionally, the Fourier-derived phases of bioluminescence fluctuate around the Fourier-derived phase of floor orientation (standard deviation of 60 min). This is in line with the progressive dampening of mean bioluminescence fluctuations (both intensity and strength) from shallower to deeper floors (see Fig. 2) and between testing periods A and B (see Fig. 1). The comparison of phases (both in the waveform and Fourier analyses) evidences how, in the second testing period, a temporal patterning in floor current-related displacement elicits a consequent temporal patterning in PMT recordings (see Fig. 5B).

Figure 5
Phases’ relationship chart showing the amplitude of all waveform peaks for PMT (one per floor) and floor information (orientation and rate of displacement), as proxies of speed and angular direction, plotted together in relation to the inertial ...


In the present study, conducted in the deep central Mediterranean basin, inertial biological rhythms in bioluminescence appear as a response to inertial hydrodynamic cycles. We detected an inertial rhythm in the bioluminescent signals of abyssopelagic organisms, equivalent to a 20-h periodicity at the latitude of South Sicily36, onto background water emissions in agreement with results from Aiello et al.37. When tides occur under the form of weak bulge displacements, the much weaker inertial wind-driven currents propagate into the deep sea as descending whirlpools, detectable in current data. This phenomenon is thoroughly explained within the framework of our previous knowledge on deep-sea biological rhythms in the Mediterranean3,5,8.

The inertial-related periodicity reveals a coherent relationship between currents and bioluminescence. This suggests mechanical stimulation of bioluminescence, as organisms transported by flowing water collide with the deployed telescope infrastructure, resulting in the defensive emission of bioluminescent light24,38,39. The presence of temporally structured peaks and troughs in average overall bioluminescence and currents suggests that drifting organisms travel in discontinuous deep-sea undular vortices. These vortices consist of chains of inertially pulsating mesoscale cyclones/anticyclones advected by a background flow, created by huge water masses that slowly propagate (at speeds of approximately 2–3 cm/s) and can change direction, as has been described by Rubino et al.35. These appear to result from a delicate and complex balance between the forces acting in the basin and the topographic constraints, which require further investigation. Our time series analysis on bioluminescence data highlighted the occurrence of two composite patterns (see waveform analysis outputs in Fig. 3), supported by the occurrence of a temporal shift in peak timing as a function of depth (see the phases’ relationship in Fig. 5 and Table 1). Taken together, these results suggest that bioluminescent plankton are carried by the downward propagation of lens-like undular vortices, the force of which increased in our testing area during the second half of the monitoring period. Unfortunately, a precise identification of the origin of the described vortices is beyond the scope of this study and would require more precise oceanographic characterisation. However, it would be interesting to explore this aspect in future work, to determine whether the mechanism of generation is local or remote and whether it is of thermohaline origin.

Deep-sea currents can elicit a bioluminescent signal through mechanical stimulation17,38, as advected organisms collide with static structures of both man-made40 and natural41 origin. In a natural setting, bioluminescence is predicted to be stimulated by impinging currents at steep geomorphological structures such as outcrops, seamount flanks and tops, cliffs and escarpments41, as well as into canyons. We would therefore also expect a periodicity in the bioluminescent signal where these geomorphologies are exposed to rhythmic flow patterns. The ecological consequence of rhythmic variations in levels of bioluminescent activity at these sites is a matter for conjecture at this point with bioluminescence serving as an attractant or a feeding cue (e.g. refs 42 and 43) as well as a repellant (e.g. ref. 44). However, the prevalence of bioluminescence in the deep sea combined with its importance as a sensory cue would certainly indicate its potential as a synchronizing ecological cue to affect a rhythmic response in benthic communities. Likewise, temporal changes in bioluminescence produced at man-made structures could locally alter the behaviour of organisms, affecting their survival within the framework of predator-prey relationships. In view of our results we would recommend more extended measurement of bioluminescence and its interactions where water flow is interrupted, such as neutrino telescopes as well as complex seabed morphologies (especially seamounts and canyons), including in those areas where cabled observatory networks have been deployed and a further technological development is envisaged (e.g. the NEPTUNE network in Barkley canyon from the Ocean Network Canada, ONC;

Endogenous rhythms are those measured in laboratory controlled conditions, when the synchronizing environmental stimulus is removed (hence rhythms appear as “free-running” with a phase that progressively shifts over consecutive day cycles away from the time window driven by the absent external synchronizer)45. Constant condition protocols are therefore used to distinguish between endogenous (persistent, free-running) and exogenous (progressively fading out) rhythms. To our knowledge, no laboratory observations of deep-sea bioluminescent organisms exist, whose activity has been tracked over consecutive days in constant conditions (i.e. with light and current cycle deprivation). In the field, the most proximate conditions simulating constant laboratory darkness as a proxy of our deep-sea results would be the Arctic in winter46. Interesting field observations on 24-h based DVMs and bioluminescence patterning has been provided at this latitude and season. In our case, the bioluminescent patterning was not spontaneous but associated to currents and therefore its endogenous component is still elusive.

The question of the identity of the organisms producing the bioluminescence remains to be resolved. Bioluminescent organisms each have a characteristic light emission in terms of its spectrum, intensity, and kinetics13. Spectral and intensity data have been reported for several deep-sea pelagic organisms18,47. However, the intensity and kinetics of the light produced is also dependent on the type and strength of the stimulus (e.g. refs 48 and 49, and in the absence of a standardized mode of stimulus, to try and identify the organisms in question from the PMT data alone is not likely to be possible. However, our results show that a greater proportion of organisms with brighter bioluminescent emissions were likely present higher within the water column than near the seafloor. This may result from different communities present at different depths. The presence of changes in community throughout the benthopelagic layer is well documented (refs 50 and 51 and literature herein). Deep-sea zooplankton with the brightest emissions include some species of copepod and decapod, as well as gelatinous organisms, including scyphozoans, siphonophores, and pyrosomes (reviewed in ref. 18, and these may have been present in greater concentrations at increasing altitude above the seafloor. In the future, additional information could be obtained from neutrino telescopes about individual bioluminescent flashes and their probable source if they were instrumented with cameras capable of imaging individual flashes, such as the ICDeep50, LuSEApher52 or sCMOS53 cameras.

Continual monitoring of deep-sea bioluminescence by neutrino telescopes at various sites within the Mediterranean Sea has the potential to make important spatial and temporal comparisons of productivity and its export to the deep sea. Supporting the findings in the present study, tidal rhythms were also observed at a deep telescope site on the coast of Greece (Pylos)39, which were not detected at Capo Passero. At the western Mediterranean ANTARES, light production is generally higher by an order of magnitude than at the NEMO site (i.e. approx. 102 here versus 103 kHz, as reported by 24 for the same month). The difference in bioluminescence signal between the two sites is similarly reflected in bioluminescence data collected in the deep sea using low light video cameras17,38,54, most likely resulting from a combination of differences in depth as well as higher levels of productivity in the western Mediterranean, fueling a greater abundance and biomass of deep-sea organisms55. Seasonal variations have been reported from the ANTARES site with observations of strong bioluminescence signals, generated from bioluminescent bacteria and other marine organisms, associated with deep-water formation in winter24,26. Seasonal variations in the abundance of bioluminescent zooplankton have been observed at the Capo Passero site from video data, with a deep peak forming in spring54. This could be attributed to the arrival of phytodetritus to the deep sea with a delay of two-three months55,56 with respect to the peak of primary production at the surface, typical of this area in late winter (e.g. see February-March satellite data from The greater food availability for deep-sea organisms likely enhances zooplankton biomass, possibly resulting in the increase in bioluminescent activity observed in the latter half of June in the present study.

Cabled observatory networks including branching neutrino telescopes open an important window for real-time, continuous, and long-lasting observation of the abyss57,58. Further work along these lines can be used to develop the use of biological data to support oceanographic studies and vice versa: the detection of bioluminescence rhythms can reveal the existence of deep-sea mesoscale features and their variability. In this context, automation of data processing per our multiple step analysis protocol inspired by chronobiology (i.e. periodogram, waveform and phase integrated chart analyses) are of relevance for the real-time and continuous monitoring of abyssal bioluminescence diel and seasonal dynamics.

Materials and Methods

The NEMO Phase 2 prototype

Cherenkov neutrino telescopes are research infrastructures whose primary goal is the detection of high-energy neutrinos of astrophysical origin, observable in the vast darkness of the deep sea59. An underwater neutrino telescope is an array of several thousand PMTs and oceanographic equipment (e.g. hydrophones and oceanographic sensors) displaced over a volume of water of about a cubic kilometre. Power and data transport is obtained using a network of electro-optical cables, connected to a main cable landing on-shore60. Thanks to their remote, real-time and multiparametric measurement capability, these structures can be opened to more multidisciplinary research and presently serve a diversified community of end-users such as oceanographers, biologists, and geophysicists61.

A prototype detection unit, the NEMO Phase 2 tower37,62 was deployed at 3450 m depth, about 90 km off-shore Capo Passero in March 2013 (Ionian Sea; Fig. 6A,B), within the activities of the NEMO (Neutrino Mediterranean Observatory) collaboration. The tower contains 32 optical modules, each comprised of a glass sphere containing a Hamamatsu R708110 Photo-Multiplier Tube (PMT), with single photon sensitivity within the 400–700 nm wavelength range, and its readout electronics63. The four optical modules are deployed at both extremities of an 8 m-long horizontal arm, two looking horizontally, two downward oriented, to optimise reconstruction of horizontal and up-going neutrino tracks. Here, we used data from pairs of horizontal-looking optical sensors for each floor. The tower rises 420 metres above the bottom and contains 8 arms representing measuring floors, separated by 40 m (i.e., from average CTD readings, 3349, 3309, 3269, 3229, 3189, 3149, 3109, and finally 3069 m depth, see Fig. 6B). PMT readout electronics provides two kinds of data: a) pulse shape of detected photon hits (signals exceeding 0.3 single photoelectrons threshold) sampled with 5 ns bin; b) average rate of photon hits recorded in 1 ms every second. The latter are used in the present work.

Figure 6
Location (red dot) and bathymetry (A) of the Capo Passero Site, hosting the KM3NeT-Italia neutrino telescope, in the Central Mediterranean (the red dot indicates the coordinates of the prototype NEMO Phase 2 tower, Ionian Sea South-East of Sicily at 36° ...

The PMT counting rate of an undersea neutrino detector is affected by constant background environmental light. The decay of radioactive elements in seawater, mainly potassium 40 (40K), causes a steady production of photons; Cherenkov radiation64,65. The emission spectra of most bioluminescent deep-sea organisms are blue (centered around 475 nm13, which also falls within the sensitivity of the PMTs, therefore also adding to the optical background in the data. Bioluminescent bacteria produce a steady glow when cell density exceeds a critical level66. When present at sufficient densities, bacterial bioluminescence adds to the Cherenkov radiation to produce a constant baseline (background) counting rate in the PMT data. Sharp peaks often several orders of magnitude higher than the steady baseline rate, are also found in the PMT data67. These peaks result from collisions between bioluminescent organisms and the telescope structure and can be correlated with current velocity64.

Data collection

We focused on bioluminescence readings for the month of June 2013 (from 03/06 at 6:57 to 30/06 at 23:57), as detected separately by PMT pairs at the 8 floors of the NEMO tower prototype (GPS coordinates: N 36° 17′48.12″, E 15° 58′45.06″, 3450 m depth). Light intensity readings obtained by both PMTs (that is the optical rate measured in 1 ms every second) were averaged together at 10 min frequency. We carried out time series pre-processing by eliminating acquisition gaps with a moving average, calculated using a variable number of precedent and successive values depending on gap duration. This time series reconstruction (in 10 min bins) did not alter the detection of significant periodicities at diel scale (i.e. encompassing tidal, inertial and day-night related rhythms) as it was carried out well below the Nyquist frequency law restrictions.

The angular movement of the tower arms has been used as a proxy for water flow direction and speed changes. The arm orientation was measured via the electronic compass installed at the fifth floor, and cross-checked with data of the acoustic positioning system: acoustic transceivers were anchored at fixed positions on the sea floor, whereas acoustic receivers (i.e. hydrophones) were fixed along the tower. The hydrophone positions were determined by triangulation, using the travel times of acoustic signals between emitters and receivers62,68. The rate of displacement of the arm orientation was also calculated.

Data analysis

Although the burst rate cannot be aligned with specific groups of organisms, a comparison was made of the bioluminescence data with respect to burst intensity to assess broad differences in the bioluminescent community at different altitudes above the seafloor. Bioluminescent emission can vary according to the strength of the stimulus and the degree of exhaustion of an individual specimen18,49,69, however a reasonable comparison of bioluminescent intensity can be made assuming similar water and current conditions at each tier of the tower. So PMT data were subdivided into ranges of intensity for each floor. Bioluminescent bursts were counted and grouped into 3 arbitrary intensity categories: 50–100, 100–150, and > 150 kHz. Values were then transformed into percentages and visually compared through the different sampling depths, to give an indication of the distribution of bioluminescent emissions with respect to intensity.

The occurrence of significant periodicities in PMT time series (i.e. all detections above 45 kHz), tower angular movements and displacements was assessed by using the chi-square periodogram procedure70, with El Temps package ( The periodogram analysis is an extension of the Buys-Ballot Table (Williams and Naylor, 197871. When a periodogram analysis is performed, the periodicity of the inherent time series is unknown. Thus, the period is screened by testing for a range of wavelengths selected to cover all rhythms of interest (e.g. tidal, inertial, and day-night related). Periodicities were screened within the interval 400–1600 min. In the periodogram output plots, the highest peak exceeding the significant threshold (p-value < 0.05) represents the maximum percentage of total data variance explained by the inherent dominant periodicity.

A peak is defined as a single increase in the time series fluctuation over a cycle time window. When this fluctuation is averaged over several cycles in order to obtain an approximated timing of time series increase, a phase can be computed via waveform analysis. A waveform analysis72 was carried out on biological and floor orientation data, by subdividing time series on the length of the dominant inertial periodicity T at the latitude ϕ of the study area, i.e. between 20 and 21-h62. Briefly, we subdivided time series into sub datasets of 20.5 h length, equivalent to 1230 min. Then a consensus curve (i.e. the waveform) was obtained by averaging values within all subset per corresponding 10 min bins. Phase duration and timing for bioluminescence were statistically assessed using the Midline Estimating Statistic Of Rhythm (MESOR)73 all waveform values were re-averaged and the resulting mean represented as a horizontal threshold line on the waveform plots. This threshold was used to identify waveform peaks as temporal amplitude of significant increments (i.e. the phase duration74. To further validate our methodology and corroborate the waveform using a different approach, an independent Fast Fourier Transform analysis (FFT) was also conducted on both biological and angular data, for the 1st and the 2nd testing period to produce corresponding results.

Outputs for waveform analysis of bioluminescence and floor orientation data were plotted together to assess the phase relationship between bioluminescence and water flow as a proxy of cause-effect relationships (i.e. water flowing against the telescope structure causing mechanical stimulation of bioluminescence in advected organisms). To do so, we presented data for the floor at which the largest biological and angular amplitude fluctuations were reported. The changes in peak timing through floors for bioluminescence rhythms were evidenced by plotting all waveform phases into a common plot57. In that chart, the temporal amplitude of each waveform (i.e. the phase made by all values above the MESOR; see above), was shown as a continuous horizontal line. To evidence any depth-related trend in the rhythms of peak timing, used as a proxy of vertical movement and downward propagating whirlpools, different phase durations were plotted together for all floors. Then, the “eye-fitting method” generally in use in chronobiology45, was applied to visually identify a depth-shift in phases. An observer marked eye-fitted straight lines, connecting the onset and offset of waveform phases through floors. This fitting approach is used to predict the timing of a significant increase in a rhythm (i.e. its peak), based on its temporal pattering during previous 24-h cycles. This is very useful when an observation must be gathered in the field or an experimental treatment must be administered in the laboratory at a certain moment of a rhythm, without a real-time knowledge of its fluctuation state. Here, we adapted this method to identify bioluminescence phase shifts through depth, as temporal difference in activity onsets and offsets at all floors.

Additional Information

How to cite this article: Aguzzi, J. et al. Inertial bioluminescence rhythms at the Capo Passero (KM3NeT-Italia) site, Central Mediterranean Sea. Sci. Rep. 7, 44938; doi: 10.1038/srep44938 (2017).

Publisher's note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.


The authors declare no competing financial interests.

Author Contributions J.A. performed time series analysis on PMT data, T.C. and A.S. provided time series analysis on angular movement of the tower arms, used as a proxy for hydrodynamic conditions, E.F. provided background on deep-sea ethology and ecology, while J.C. provided detailed knowledge on emitting organisms. A.C., G.R., M.G.P. and S.V. from KM3NeT-Italia/NEMO Collaboration provided data and important information for the development of the paper. All those authors contributed to the writing of the paper. All remaining authors were instead involved in the development, maintenance of infrastructures, as well as in data elaboration. Therefore, their co-authorship is due within the framework of established data diffusion policies of the KM3NeT-Italia/NEMO Collaboration. These authors are: S.A., F.A., M.A., G.B., E.B., N.B., S.B., A.B., B.B., C.B., G.C., M.C., C.C., A.C., F.C., S.C., A.C., T.C., M.C., R.C., R.C., M.C., G.C., C.D.A., A.D.A., G.D.B., V.D.L., N.D., C.D., L.S. D.M.,P.F., G.F., V.F., L.A.F., F.G., V.G., A.G., R.G., G.G., C.H., M.I., V.K.G.L., D.L., K.P., E.L., P.L., A.L., F.L., D.L.P., E.M., A.M., A.M., A.M., R.M., R.M., P.M., E.M., A.M., C.M.M., M.M., M.M., P.M., M.M., C.A.N., A.O., A.O., R.P., C.P., M.G.P., C.P., P.P., E.P.,C.P., S.P., F.R., N.R., G.R., A.R., M.S., P.S., V.S., I.S., F.S., V.S., F.S., A.S., M.S., S.M.S., M.T., G.T., L.T., A.T., F.V., P.V., S.V., D.V.


  • Robison B. H. Deep pelagic biology. J. Exp. Mar. Biol. Ecol. 300, 253–272 (2004).
  • Widder E. A. Bioluminescence and the pelagic visual environment. Mar. Freshwat. Behav. Physiol. 35, 1–26 (2002).
  • Aguzzi J. & Company J. B. Chronobiology of deep-water decapod crustaceans on continental margins. Adv. Mar. Biol. 58, 155–225 (2010). [PubMed]
  • Haddock S. H. D., Moline M. A. & Case J. F. Bioluminescence in the sea. Ann. Rev. Mar. Sci. 2, 443–493 (2010). [PubMed]
  • Aguzzi J. et al. . Activity rhythms in the deep-sea: a chronobiological approach. Front. Biosci. 16, 131–150 (2011). [PubMed]
  • Wagner H. J., Kemp K., Mattheus U. & Priede I. Rhythms at the bottom of the deep sea: Cyclic current flow changes and melatonin patterns in two species of demersal fish. Deep-Sea Res. I 54, 1944–1956 (2007).
  • Aguzzi J. et al. . Behavioral rhythms of hydrocarbon seep fauna in relation to internal tides. Mar. Ecol. Prog. Ser. 418, 47–56 (2010).
  • Aguzzi J., Puig P. & Company J. B. Hydrodynamic, non-photic modulation of biorhythms in the Norway lobster, Nephrops norvegicus (L.). Deep-Sea Res. I 56, 366–373 (2009).
  • Sbragaglia V., García J. A., Chiesa J. J. & Aguzzi J. The effect of simulated tidal currents on burrow emergence rhythms of the Norway lobster (Nephrops norvegicus, L.). Mar. Biol. 162, 2007–2016 (2015).
  • Duñez J. D., Sbragaglia V., García J. A., Company J. B. & Aguzzi J. First laboratory insight on the behavioral rhythms of the bathyal crab Geryon longipes. Deep-Sea Res. I 116, 165–173 (2016).
  • Webster F. Observations of inertial-period motions in the deep sea. Rev. Geoph. 6, 447–472 (1968).
  • Pugh D. T. Tides, surges, and mean sea-level (John Wiley and Sons, 1987).
  • Widder E. A. Bioluminescence in the ocean: origins of biological, chemical, and ecological diversity. Science 328, 704–708 (2010). [PubMed]
  • Widder E. A. et al. . Bioluminescence in the Monterey Submarine Canyon: image analysis of video recordings from a midwater submersible. Mar. Biol. 100, 541–551 (1989).
  • Priede I. G., Bagley P. M., Way S., Herring P. J. & Partridge J. C. Bioluminescence in the deep sea: Free-fall lander observations in the Atlantic Ocean off Cape Verde. Deep-Sea Res. I 53, 1272–1283 (2006).
  • Rivers T. J. & Morin J. G. The relative cost of using luminescence for sex and defense: Light budgets in cypridinid ostracods. J. Exp. Biol. 215, 2860–2868 (2012). [PubMed]
  • Priede I. G., Jamieson A., Heger A., Craig J. & Zuur A. F. The potential influence of bioluminescence from marine animals on a deep-sea underwater neutrino telescope array in the Mediterranean Sea. Deep-Sea Res. I 55, 1474–1483 (2008).
  • Craig J., Priede I. G., Aguzzi J., Company J. B. & Jamieson A. J. Abundant bioluminescent sources of low light intensity in the deep Mediterranean Sea and North Atlantic Ocean. Mar. Biol. 162, 1637–1649 (2015a).
  • Johnsen S., Frank T. M., Haddock S. H., Widder E. A. & Messing C. G. Light and vision in the deep-sea benthos: I. Bioluminescence at 500–1000 m depth in the Bahamian Islands. J. Exp. Biol. 215, 3335–3343 (2012). [PubMed]
  • Nosengo N., Riccobene G. & Pavan G. The neutrino and the whale. Nature 462 (2009). [PubMed]
  • Riccobene G. et al. . 2009, Long-term measurements of acoustic background noise in very deep sea. Nuc. Instr. Meth. Phys. Res. A 604, 149–157 (2009).
  • Sciacca V. et al. . Annual Acoustic Presence of Fin Whale (Balaenoptera physalus) Offshore Eastern Sicily, Central Mediterranean Sea. PLoS ONE 10, e0141838, doi: org/10.1371/journal.pone.0141838 (2015). [PMC free article] [PubMed]
  • Caruso F. et al. . Size distribution of sperm whales acoustically identified during long term deep sea monitoring in the Ionian Sea. PLoS ONE 10, e0144503 doi: (2015).10.1371/journal.pone.0144503 [PMC free article] [PubMed] [Cross Ref]
  • Tamburini C. et al. . Deep-sea bioluminescence blooms after dense water formation at the ocean surface. PLoS ONE 8, e67523, doi: org/10.1371/journal.pone.0067523 (2013). [PMC free article] [PubMed]
  • Martini S. et al. . Temporal survey, presence and activity of bioluminescent bacteria at the deep ANTARES station (Mediterranean Sea, 2200 m depth). Rapp. Comm. int. Mer Médit. 40, 422 (2013).
  • Martini S., Nerini D. & Tamburini C. Relation between deep bioluminescence and oceanographic variables: A statistical analysis using time-frequency decompositions. Prog. Oceanog. 127, 117–128 (2014).
  • Riccobene G. et al. . Deep seawater inherent optical properties in the Southern Ionian Sea, Astroph. Phys. 27, 1–9 (2007).
  • Krom M. D., Herut B. & Mantoura R. F. C. Nutrient budget for the Eastern Mediterranean: Implications for phosphorus limitation. Limnol. Oceanogr. 49, 1582–1592 (2004).
  • Bosc E., Bricaud A. & Antoine D. Seasonal and interannual variability in algal biomass and primary production in the Mediterranean Sea, as derived from 4 years of SeaWiFS observations. Glob. Biogeochem. Cy. 18, GB1005 (2004).
  • KM3NeT Collaboration. Conceptual design report for a Deep-Sea research infrastructure incorporating a very large volume Neutrino Telescope in the Mediterranean Sea. (2008).
  • KM3NeT Collaboration. Technical Design Report for a Deep-Sea Research Infrastructure in the Mediterranean Sea.;; (2011).
  • Roether W. et al. . Recent changes in eastern Mediterranean deep waters. Science 271, 333–335 (1996).
  • Klein B. et al. . The large deep water transient in the eastern Mediterranean. Deep-Sea Res. I 46, 371–414 (1999).
  • Manca B., Budillon G., Scarazzato P. & Ursella L. Evolution of dynamics in the eastern Mediterranean affecting water mass structures and properties in the Ionian and Adriatic Seas. J. Geoph. Res. 108, 8102 (2003).
  • Rubino A. et al. . Abyssal undular vortices in the Eastern Mediterranean basin. Nat. Comm. 3, 834, doi: (2012).10.1038/ncomms1836 [PubMed] [Cross Ref]
  • Gasparini G. P. et al. . Tidal and subtidal currents in the strait of Sicily. J. Geoph. Res. 109, C02011 (2004).
  • Aiello S. et al. . Measurement of the atmospheric muon depth intensity relation with the NEMO Phase-2 tower. Astroph. Phys. 66, 1–7 (2015).
  • Craig J., Jamieson A. J., Heger A. & Priede I. G. Distribution of bioluminescent organisms in the Mediterranean Sea and predicted effects on a deep-sea neutrino telescope. Nuc. Instr. Met. Phys. Res. A 602, 224–226 (2009).
  • van Haren H., de Jong M. & Kooijman P. Yearlong moored bioluminescence and current data at KM3NeT neutrino telescope sites in the deep Ionian Sea. Astroph. Phys. 67, 1–7 (2015).
  • Amram P. et al. . Background light in potential sites for the ANTARES undersea neutrino telescope. Astroph. Phys. 13, 127–136 (2000).
  • Craig J., Jamieson A. J., Bagley P. M. & Priede M. Naturally occurring bioluminescence on the deep-sea floor. J. Mar. Syst. 88, 563–567 (2011a).
  • Gordon J. D. M., Bergstad O. A. & Pascoe P. L. The influence of artificial light on the capture of deep-water demersal fish by bottom trawling. J. Mar. Biol. Assoc. UK 82, 1–6 (2002).
  • Nishida S., Ohtsuka S. & Parker A. Functional morphology and food habits of deep-sea copepods of the genus Cephalophanes (Calanoida: Phaennidae): perception of bioluminescence as a strategy for food detection. Mar. Ecol. Prog. Ser. 227, 157–171 (2002).
  • Heger A. et al. . Benthic bioluminescence in the bathyal North East Atlantic: luminescent responses of Vargula norvegica (Ostracoda: Myodocopida) to predation by the deep-water eel (Synaphobranchus kaupii). Mar. Biol. 151, 1471–1478 (2007).
  • Refinetti R. Circadian physiology (CRC Press, 2006).
  • Berge J. et al. . Bioluminescence in the high Arctic during the polar night. Mar. Biol. 159, 231–237 (2012). [PMC free article] [PubMed]
  • Widder E. A., Latz M. I. & Case J. F. Marine bioluminescence spectra measured with an optical multichannel detection system. Biol. Bull. 165, 791–810 (1983).
  • Bowlby M. R., Widder E. A. & Case J. F. Patterns of stimulated bioluminescence in two pyrosomes (Tunicata: Pyrosomatidae). Biol. Bull. 179, 340–350 (1990).
  • Bowlby M. R. & Case J. F. Flash kinetics and spatial patterns of bioluminescence in the copepod Gaussia princeps. Mar. Biol. 110, 329–336 (1991).
  • Cartes J. E., Fanelli E., Lopez-Perez C. & Lebrato M. The distribution of deep-sea macroplankton (over 400 to 2300 m) at intermediate and near bottom waters: relationships with hydrographic factors. J. Mar. Syst. 113-114, 75–87 (2013).
  • Craig J., Youngbluth M., Jamieson A. J. & Priede I. G. Near seafloor bioluminescence, macrozooplankton and macroparticles at the Mid-Atlantic Ridge. Deep-Sea Res. I 98, 62–75 (2015b).
  • Dominjon A. et al. . An ebCMOS camera system for marine bioluminescence observation: The LuSEApher prototype. Nuc. Instr. Met. Phys. Res. A 695, 172–178 (2012).
  • Phillips B. T. et al. . Observations of in situ deep-sea marine bioluminescence with a high-speed, high-resolution sCMOS camera. Deep-Sea Res. I 111, 102–109 (2016).
  • Craig J., Jamieson A. J., Bagley P. M. & Priede M. Seasonal variation of deep-sea bioluminescence in the Ionian Sea. Nuc. Instr. Met. Phys. Res. A 626, 115–117 (2011b).
  • Cartes J. E., Fanelli E., Papiol V. & Zucca L. Distribution and diversity of open-ocean, near-bottom macrozooplankton in the western Mediterranean: analysis at different spatio-temporal scales. Deep-Sea Res. I 57, 1485–1498 (2010).
  • Fanelli E., Cartes J. E. & Papiol V. Trophodynamics of zooplankton fauna on the Catalan slope (NW Mediterranean): insight from δ13C and δ15N analysis. J. Mar. Syst. 87, 79–89 (2011).
  • Aguzzi J. et al. . Challenges to assessment of benthic populations and biodiversity as a result of rhythmic behaviour: Video solutions from cabled observatories. Oceanogr. Mar. Biol. 50, 235–286 (2012).
  • Favali P., Beranzoli L. & De Santis A. Seafloor observatories: A new vision of the earth from the abyss (Springer, 2015).
  • Sapienza P. & Riccobene G. High-energy neutrino astronomy. Rivista del Nuovo Cimento 32, 591–649 (2009).
  • Ameli F. et al. . The data acquisition and transport design for NEMO phase 1. IEEE Trans. Nuc. Sci. 55, 233–240, doi: (2008).10.1109/TNS.2007.914035 [Cross Ref]
  • Migneco E., Piattelli P. & Riccobene G. Underwater neutrino telescopes: Detectors for astroparticle physics and a gateway for deep-sea laboratories in Seafloor Observatories: A New Vision of the Earth from the Abyss 23–57 (Springer, 2015).
  • Adrian Martinez S. et al. . Long term monitoring of the optical background in the Capo Passero deep-sea site with the NEMO tower prototype, Europ. Phys. J. C 76, 1–11 (2016).
  • Aiello S. et al. . The optical modules of the phase-2 of the NEMO project. J. Instrument. 8., 7001–7021 (2013).
  • Ameli F., Bonori M. & Massa F. Optical background measurement in a potential site for the NEMO KM undersea neutrino telescope. Europ. Phys. J. C 25, 7567–7570 (2002).
  • Massa F. Optical radiation background from 40K decays in undersea neutrino telescopes. Europ. Phys. J. C 22, 749–756 (2002).
  • Miller M. B. & Bassler B. L. Quorum sensing in bacteria. Ann. Rev. Microbiol. 55, 165–199 (2001). [PubMed]
  • van Haren H. et al. . Acoustic and optical variations during rapid downward motion episodes in the deep north-western Mediterranean Sea. Deep-Sea Res. I 58, 875–884 (2011).
  • Viola S. et al. . NEMO-SMO acoustic array: A deep-sea test of a novel acoustic positioning system for a km(3)-scale underwater neutrino telescope. Nuc. Instr. Met. Phys. Res. A 725, 207–210 (2013).
  • Latz M. I., Bowlby M. R. & Case J. F. Recovery and stimulation of copepod bioluminescence. J. Exp. Mar. Biol. Ecol. 136, 1–22 (1990).
  • Sokolove P. G. & Bushell W. N. The chi square periodogram: Its utility for analysis of circadian rhythms. J. Theoret. Biol. 72, 131–160 (1978). [PubMed]
  • Williams J. A. & Naylor E. A procedure for the assessment of significance of rhythmicity in time-series of data. Int. J. Chronobiol. 5, 435–444 (1978).
  • Chiesa J. J., Aguzzi J., Garcia J. A., Sardà F. & de la Iglesia H. Light intensity determines temporal niche switching of behavioral activity in deep water Nephrops norvegicus (Crustacea: Decapoda). J. Biol. Rhyth. 25, 277–287 (2010). [PubMed]
  • Aguzzi J., Bullock N. & Tosini G. Spontaneous internal desynchronization of locomotor activity and body temperature rhythms from plasma melatonin rhythm in rats exposed to constant dim light. J. Circad. Rhyth. 4, 6, doi: (2006).10.1186/1740-3391-4-6 [PMC free article] [PubMed] [Cross Ref]
  • Schlitzer R. Ocean Data View, (2015).

Articles from Scientific Reports are provided here courtesy of Nature Publishing Group