alphaFeedback

WorksPhysics and Astronomy

Short spatial mooring-tilt variations from deep Mediterranean observations

AI contributions
Not recorded for this paper
View details

Abstract

Interaction between energy-abundant mesoscale eddies and internal waves can lead to convection-turbulence generation and may prove important for deep-sea life and circulation. However, the size of scales of interacting flows is not well known. In this paper, a diagnostic tool of tilt is tested near the single top-buoyancy of 40 mooring lines 9.5 m apart horizontally and compared with 50-m scale relative vorticity and waterflow above a 2500-m deep flat Northwestern-Mediterranean seafloor. Whilst tilt relates to first order with flow-speed squared induced by mooring-line drag, considerable deviations from this relationship and larger tilt occur when the amplitude of relative vorticity attains values O(f), f the inertial frequency of planetary vorticity. The sign of relative vorticity is of no importance. Such larger-tilt events occur during most intense convection turbulence via warm-water slanting from above. During these events, variations in tilt-angle magnitude are as large as the average tilt O(0.1)degree, thereby reducing variational scales from 50 to 9.5 m. Thus, in a deep-sea environment where flow speeds are <0.07 m s^-1, O(0.01) m s^-1 flow-speed variations provide important turbulent mixing, without deep dense-water formation.

Introduction

In a stably stratified environment like the sun-heated ocean, downward pulses of warm water seem impossible in terms of irreversible turbulent convection. In the deep weakly stratified Northwestern Mediterranean Sea however, such pulses are observed to dominate local turbulent mixing about 40% of time, governed by sub-mesoscale eddies [22] and inertio-gravity waves IGW [29], underneath a large-scale boundary flow that varies at mesoscales [4]. The convection-type turbulence is found on average three times more intense than geothermal heating through the seafloor, and not associated with deep dense-water formation that was absent [29].

The region also shows extensive generation of near-inertial waves, which under conditions of weakly stratified waters can lead, like eddies, to slantwise convection, as has been proposed in general [13, 21] and inferred from shipborne profiling observations in the open Western Mediterranean [26]. The convection may develop to irreversible 3D turbulence, while vertical (opposite-to-gravity) density profiles appear stably stratified [21, 20]. Under such conditions, the traditional internal-wave bounds [f,N][f, N] of inertial f=2Ωsin⁡φf = 2\Omega\sin\varphi and buoyancy NN frequencies extend to IGW bounds [ωmin⁡ ωmax⁡][\omega_{\min}\ \omega_{\max}] [12].

Here, φ\varphi denotes latitude, and Ω\Omega the Earth rotational frequency. The ωmin⁡<f\omega_{\min}<f and ωmax⁡>2Ω\omega_{\max}>2\Omega, for N<2ΩN < 2\Omega, are functions of NN, or shear magnitude ∣S∣∼N|S| \sim N in regularly observed marginal stability conditions (e.g., [27]), φ\varphi, and of direction of wave propagation [12, 8].

ωmax⁡, ωmin⁡=(A±(A2−B2)1/2)1/22,(1)\omega_{\max},\ \omega_{\min} = \frac{\left(A \pm\left(A^2-B^2\right)^{1/2}\right)^{1/2}}{\sqrt{2}}, \tag*{(1)}

in which A=N2+f2+fs2A = N^2 + f^2 + f_s^2, B=2fNB = 2fN, and fs=fhsin⁡αφf_s = f_h\sin\alpha_\varphi, fh=2Ωcos⁡φf_h = 2\Omega\cos\varphi, and αφ\alpha_\varphi the angle to zonal direction.

In the deep Mediterranean N=O(f)N = O(f). The associated IGW band (1) thus extends free-propagating wave frequencies to sub-inertial, sub-mesoscales. In theory, the bounds may be further extended locally under horizontal flow differences generating the vertical component of relative vorticity (e.g., [11]),

ζ=∂V∂x−∂U∂y,(2)\zeta= \frac{\partial V}{\partial x} - \frac{\partial U}{\partial y}, \tag*{(2)}

for horizontal waterflow components (UU, VV) in East and North directions, respectively. Depending on its sign, this relative vorticity may add to or subtract from planetary vorticity component ff. In the Northwestern Mediterranean values of ∣ζ∣=f/2|\zeta| = f/2 are reported from mid-depth drifter observations [22]. The question is what its value amounts to in the deep, very weakly stratified Mediterranean and what its scales are compared to that of the motions it might trap or reflect.

In this paper, observational investigations of convective water pulses are further pursued that occur frequently in the deep Northwestern Mediterranean. An unusual tool of vertical mooring-line tilt is tested for its usefulness to study vibrations that may associate with convection turbulence. For this purpose, an additional tool was installed at two levels in a 3D mooring array [28]. At the top and bottom of each of the closely spaced 45 mooring lines a tilt sensor was installed. Although typical tilt-angle variations are only ±0.1∘\pm0.1^\circ to the vertical in waterflow with speeds <0.07 m s−1< 0.07\ \mathrm{m}\ \mathrm{s}^{-1}, they reflect large-scale dynamical variations that are driven by mesoscale processes. As will be demonstrated, between the lines at 9.5 m distance however, variations are similar as at large scales, which reflects the importance of small-scale variations.

Materials and Methods

A ‘large-ring mooring’ array was deployed at the <1∘<1^\circ flat and 2458-m deep seafloor of 42∘ 49.50′N42^\circ\ 49.50^\prime\mathrm{N}, 006∘ 11.78′E006^\circ\ 11.78^\prime\mathrm{E} just 5 km south of the foot of the steep continental slope of the Northwestern Mediterranean Sea (fig:1a), in October 2020. Details of mooring layout and deployment are given in van Haren et al. (2021) and of data processing in van Haren et al. (2026). In summary, 65 high-resolution NIOZ4 temperature sensors (van Haren, 2018) were taped at 2-m intervals to each of 45 vertical lines of 0.006-m diameter and 125-m tall that were tensioned to 1.3 kN by a single 0.4-m-diameter cylinder buoy on top. The lines were attached at 9.5-m horizontal intervals to a steel-cable grid that was tensioned inside a 70-m diameter steel-tube ring (fig:1b), which functioned as a float prior to deployment and as an anchor after landing at the seafloor.

Every line held two T-sensors, one just below the buoy (Figure 2) and another near the cable-grid, of which the tilt sensor was also activated. These T-tilt sensors effectively sampled at a rate of once per 2.4 s. This allowed for the monitoring of the mooring-line deviation from vertical at two levels, a lower at about h=1.5h = 1.5 m above seafloor and an upper at about h=125.5h = 125.5 m. Due to unknown causes all T-sensors switched off unintentionally when the file-size on the memory card reached 30 MB after 7500 startups for writing a 4-kB data block. It implied that a maximum of about three months of temperature-tilt data was obtained. (The temperature-only sensors recorded about 20 months of data.) In total 8 out of 90 tilt sensors failed for various electronic reasons. The analysis in this paper focuses on the upper-level T-sensor data, and their short variations with time and space. Besides time series and spectral analyses, 2D and quasi-3D movies are presented, for their construction see [29].

Large-ring mooring after unfolding underwater

Figure 2. Large-ring mooring after unfolding underwater. (a) Three buoys on top of mooring lines. Underneath the left buoy in front, three small instruments are seen taped to the mooring line, up-down: anode for synchronization, uppermost temperature T-only sensor, upper T-tilt sensor. Strapped to the right buoy a downlooking current meter can be seen. Buoyancy and drag forces to tilt the mooring line from rest over angle α\alpha. Many small suspended creatures and matter are visible as dust-ticks. Video still by ROV Victor (Ifremer, France). (b) Schematic of full array with correct aspect ratio (artist impression by J. van Bennekom).

Although tilt is recorded along three axes to provide a vector, only its magnitude, more precisely the relative magnitude of its angle to the vertical $|\alpha|1, is considered here as a proxy for mooring vibrations and the various dynamic waterflow processes causing them. While this follows from practical reasons since a recent calibration is not available and electronic drift is about 0.1^\circ$ per 3 months, the tilt-vector direction is of no use since the mooring line can rotate freely due to a swivel. The T-tilt data are referenced to current-meter tilt data obtained simultaneously at the upper level only (Appendix A).

Three buoys held a 2-MHz single-point Nortek AquaDopp acoustic current meter recording at a rate of once per 600 s waterflow and tilt (fig:1b; fig:2a). The triangular setup in the same horizontal plane affords the possibility to compute ζ\zeta in (2). While the current-meter data are considerably noisier and less resolved than the T-sensor data, their tilt sensor shows less bias from electronic drift. While the nominal instrumental error in horizontal waterflow speed measurements is 0.0090.009 m s−1^{-1} for each 600-s sample, the acoustic deep-sea measurements require noise-filtering with cut-off at 4 cpd (cycles per day), which reduces the precision or one standard deviation error to 0.00150.0015 m s−1^{-1}. Similarly, the noise-reduced precision in tilt-angle magnitude amounts about 0.004∘0.004^\circ for current-meter measurements and 0.002∘0.002^\circ for T-sensor measurements. (T-sensor temperature is a very low-noise measurement, and only requires filter cut-offs at O(1000)O(1000) cpd.)

General tilt expectations for large-ring mooring array

A vertical-line mooring with a single buoy on top is not displaced exactly like a stiff pendulum, but in a slight J-shape with largest line-tilts near the anchor. In rest with zero waterflow, the upright position is fully due to buoyancy force FbF_{\mathrm{b}} in negative direction of gravity. For non-zero waterflows that are assumed most horizontal, the drag force FdF_{\mathrm{d}} by resistance of the mooring line and elements is in disbalance with FbF_{\mathrm{b}} because they are initially perpendicular to each other from rest (fig:2a). A balance of force (amplitudes) is obtained when, to first order [14, 23, 5],

∣Fbtan⁡α∣=∣Fd∣=CρA∣U∣2,(3)\lvert F_{\mathrm{b}}\tan\alpha\rvert= \lvert F_{\mathrm{d}}\rvert= C\rho A\lvert U\rvert^2, \tag*{(3)}

where ∣U∣=(U2+V2)1/2\lvert U\rvert= (U^2 + V^2)^{1/2} denotes the waterflow amplitude, AA the surface of flow-obstructing objects, ρ≈1027 kg m−3\rho\approx1027\ \mathrm{kg\ m}^{-3} the density of seawater and C≈1C \approx1 the drag coefficient for a cylinder object. (More about distributed drag on mooring lines can be found in the mooring design software by R. Dewey; https://rkdewey.github.io/Mooring-Design-and-Dynamics/mdd.html, last accessed 18 August 2026.)

For ∣Fb∣=1300 N\lvert F_{\mathrm{b}}\rvert= 1300\ \mathrm{N} of a single mooring line [28], A≈125×0.006+1.4×0.4=1.3 m2A \approx125 \times0.006 + 1.4 \times0.4 = 1.3\ \mathrm{m}^2 surface of mooring line of height h=125 mh = 125\ \mathrm{m} plus buoy, and typical flow speed of ∣U∣=0.05 m s−1\lvert U\rvert= 0.05\ \mathrm{m\ s}^{-1}, (2) provides an expected tilt angle-magnitude of ∣α∣=3.3/1300=0.15∘\lvert\alpha\rvert= 3.3/1300 = 0.15^\circ. Hereby, ∣U∣\lvert U\rvert is taken constant over the entire mooring line. A more realistic flow includes amplitude and direction changes in the vertical, also over h=125 mh = 125\ \mathrm{m} in the deep Mediterranean. A change of ΔU=±0.01 m s−1\Delta U = \pm0.01\ \mathrm{m\ s}^{-1} gives Δα≈±0.05∘\Delta\alpha\approx\pm0.05^\circ around above expected value of ∣α∣\lvert\alpha\rvert.

The displaced height of the buoy is then h(1−cos⁡∣α∣)<0.01 mh(1 - \cos\lvert\alpha\rvert) < 0.01\ \mathrm{m} and thus negligible for typical values of waterflow and induced tilt. (Typical horizontal excursions are hsin⁡∣α∣≈0.3 mh\sin\lvert\alpha\rvert\approx0.3\ \mathrm{m}.) These height displacements are much smaller than the height variations Δh=±0.4 m\Delta h = \pm0.4\ \mathrm{m} due to the doming of the large-ring mooring array’s cable grid, with largest heights in the centre of the grid [25].

The orientation of the mooring array, pointing towards the north-northwest (fig:1b), demonstrates that the three vertical lines equipped with current meters on top form a triangle in the horizontal plain of which the base is closely aligned to the west-east direction. Line 14 is located approximately due north, its perpendicular intersecting the baseline approximately halfway between lines 57 and 35. The doming of the steel-cable grid suggests that line 35 is highest, as it is furthest from the large ring, about 0.5 m above line 14 (lowest) and about 0.2 m above line 57.

Results

The general overview of three months of stored tilt and other parameter data from h≈125h \approx125 m above seafloor is presented in Figure 3. It demonstrates time variability in waterflow which is dominated by mesoscale (∼\simmonthly), sub-mesoscale (∼\sim1--10 days) and near-inertial variations with time (fig:3a). Waterflow speeds seldom exceed 0.06 m s−1^{-1}. The sub-inertial <1<1-cpd low-pass filtered dominant U-component data from the three current meters demonstrate consistently similar waterflow speeds when ∣U∣≲0.03|U| \lesssim0.03 m s−1^{-1} (fig:3a), near-inertial flow amplitudes are largest (fig:3a), temperature is near minimum value of Θ=12.825∘C\Theta= 12.825^{\circ}\mathrm{C} (fig:3b) and when tilt is close to its relative base-value of about 0.07∘0.07^{\circ} (fig:3c).

Three months timeseries of waterflow and tilt at the top of the mooring array

Figure 3. Three months timeseries of waterflow and tilt at the top of the mooring array. (a) Total noise-filtered waterflow speed (black) measured at line 57, and sub-inertial 1-cpd low-pass filtered East-West component for lines 57 (blue), 35 (cyan) and 14 (magenta). (b) Conservative Temperature ([9] 2010) measured at upper T-sensor. (c) As a., but for relative tilt-angle amplitude from line 57, in comparison with upper T-sensor tilt data that are not-detrended noise filtered (green) and its sub-inertial data (yellow). (d) Effective inertial frequency of combined vertical components of planetary and local relative vorticity for sub-inertial data, scaled with ff. (e) As c., but for tilt difference between lines 57 and 14 from T-sensor (green; yellow) and current meter (blue; red), in comparison with current-meter tilt difference between lines 57 and 35 (cyan; magenta). The record pairs are arbitrarily shifted along the yy-axis for display purposes.

Larger tilt generally associates with larger ∣U∣2|U|^2 as in (3), but exceptions occur, visibly e.g. on days 330 and 386. Larger tilt generally associates with higher temperatures (fig:3b), including on days 330 and 386, and which relates with larger turbulence events [29]. Mostly when ∣U∣≳0.05|U| \gtrsim0.05 m s−1^{-1}, the three waterflow speed measurements diverge, with the one from line 14 being often higher than one or both others by up to 0.01 m s−1^{-1} during the presented 3-month period. (Later episodes without T-tilt measurements showed also higher speeds at other instruments.) As line 14 is north of the other two that are closely aligned to the west-east axis, this may associate with dominant zonal flow (fig:3a). The fact that line 14 is about 0.5 m lower than line 35 is not expected to cause great effect, as it would imply a vertical current shear of up to 0.02 s−1^{-1}, or 200 times N∼fN \sim f, which is unrealistic.

The array’s horizontal waterflow differences result in non-negligible ζ\zeta, which adds to the planetary vorticity component ff for motions at horizontal scales smaller than 50 m (fig:3d). The effective Coriolis frequency, feff=(f2+fζ−∂U/∂x⋅∂V/∂y)1/2f_{\mathrm{eff}} = (f^2 + f\zeta- \partial U/\partial x \cdot\partial V/\partial y)^{1/2} [11], varies by about ±f\pm f, and has values <f<f when waterflows are generally eastward, and >f>f when waterflows are westward. Tilt-enhancement of 0.10.1–0.3∘0.3^{\circ} above base-value (fig:3c) occurs under both fefff_{\mathrm{eff}}-conditions and thus relates with ∣ζ∣|\zeta|, seemingly independent of the sign of ζ\zeta. However, the horizontal difference in tilt (fig:3e) changes sign with fefff_{\mathrm{eff}} and thus relates with ζ\zeta especially consistently in the west-east direction between lines 57 and 35. As tilt is a magnitude without direction sensitivity, this implies that waterflows to the east with negative relative vorticity induce more tilt on the southwest-side of the array compared to its east-side, and vice-versa. Considering also the northern line 14, results are more ambiguous. However, it is unclear what precisely causes the tilt differences, as the wakes from the vertical lines are expected to induce relatively weak turbulence due to the small diameters of line and buoy so that vortex shedding is not expected to be noticeable over scales O(0.01)O(0.01) m for the line and O(1)O(1) m for the buoy.

Although tilt-angle amplitude associates with temperature variations, no direct quantitative connections can be made with turbulence processes and possibly their intensity. This is because the tilt measurements are limited in resolution of the full turbulence range, capturing only the largest eddies at best (Appendix B; Figure A2). Most of registrations is resolving the IGW and sub-mesoscale motions governing the turbulence processes. However, the white-noise part of tilt demonstrates qualitative information on turbulence intensity (Figure A3).

Spectrum showing low-pass filter cutoffs for T-sensor tilt measurements

Figure A2. Spectrum showing low-pass filter cutoffs for T-sensor tilt measurements, applied for those of line 57. Original tilt data (black), noise filtered data with cut-off at 10 cpd (blue), and sub-inertial data with cut-off at 1 cpd (green). (The noise-filter cut-off for current meter tilt is at 4 cpd). The vertical dashed lines indicate the IGW bounds for N=1.0fN = 1.0f, 2Ω2\Omega indicates the semidiurnal frequency.

“Noise” high-pass filtered data and spectral power plots

Figure A3. “Noise”, 1000-cpd high-pass filtered ‘hf’ data from upper T-sensor tilt measurements of line 57. (a) Three-month time series of standard deviation of noise. Periods are indicated from which data in panels b.–d. are extracted. (b) Magnification of typical 24000-s time series from period between days 318–322. (c) As b., but for days 322.5–326.5. (d) As b., but for days 327–331. (e) Spectral power from high-frequency parts of periods in b. (black), c. (blue) and d. (green). The records are not deliberately off-set vertically.

Spectra in Figure 4 confirm the association of tilt with flow amplitude, as these show a distinct peak around ff, and a particularly steep slope p<−2p < -2 of power-law ωp\omega^p in the sub-mesoscale frequency range 0.1<ω<0.40.1 < \omega< 0.4 cpd, which is part of the sub-inertial range ω<f\omega< f. (This slope extends to the range 0.1<ω<0.80.1 < \omega< 0.8 cpd for the kinetic energy spectrum, not shown, which has a sharper ff-peak than ∣U∣|U| that is slightly blue-shifted from ff rather than ∣U∣|U|'s red-shift.) That sub-inertial range is at frequencies ω<ωmin⁡(N=f)\omega< \omega_{\min}(N=f) largely below those of the IGW band, for N=fN = f, and considered to largely contain energy from sub-mesoscale eddy motions. However, given the variability in effective inertial frequency, hence relative vorticity, and regularly occurring feff→0f_{\mathrm{eff}} \to0 (fig:3d), it may include freely propagating small-scale IGW at sub-inertial frequencies as ωmin⁡→0\omega_{\min} \to0.

Three-month mean spectra comparing waterflow speed, relative vorticity, current-meter tilt-angle magnitude and T-sensor temperature

Figure 4. Three-month mean spectra from data at h≈125h \approx125 m above seafloor for waterflow speed squared (blue) and 10-cpd low-pass filtered vertical component of relative vorticity (magenta), are compared with current-meter tilt-angle magnitude (green) and its horizontal difference (red; vertically shifted by a factor of 0.45), and with T-sensor temperature (cyan). Indicated are inertio-gravity wave ‘IGW’ frequency bounds [ωmin⁡,ωmax⁡][\omega_{\min},\omega_{\max}] (for buoyancy frequency N=fN=f, the inertial frequency) and semidiurnal frequency 2Ω2\Omega. Several slopes pp of spectral-model power-laws ωp\omega^p are shown (see text). The reader may use a transparent ruler or set square with parallel lines to verify model slopes with observed spectral slopes.

The best-fit slope for ∣α∣|\alpha|- and ∣U∣|U|-spectra is around p=−11/5p = -11/5, which complies with the power-law model of buoyancy-driven motions by [3] and [16], henceforth BO, but errors are considerable. Nonetheless, deviations tend away from other power-law model-slopes like p=−2p = -2 for finestructure contamination ([17]; [18]) and internal waves established for open seas away from boundaries at f≪ω≪Nf \ll\omega\ll N ([7]), p=−5/3p = -5/3 established for the inertial

@bibliography keys subrange of passive-scalar shear-induced turbulence ([10]; [15]), and, not shown, less steep p=−1p=-1 for intermittency of a deterministic chaos process ([19]; [2]).

In contrast, the tilt’s horizontal difference shows a distinctly different spectrum, with a weaker near-inertial peak (Figure 4). Except for an unknown sub-peak at 0.5 cpd, its spectral slope in the sub-inertial, sub-mesoscale frequency range matches p=−7/5p=-7/5, is significantly different from that of tilt and kinetic energy, and complies with the BO-model for active scalars (and which thus have a different slope from that of the associated advective-flow kinetic energy). However, it matches the sub-mesoscale slope of relative vorticity, and that of its magnitude, and it matches the slope of scalar quantity temperature, which suggests dominant buoyancy-driven motions.

The limited spectral information from only three current meter data records can be extended by exploring variations in tilt and tilt difference between various T-sensors at different lines (Figure 5). Tilt is spectrally similarly varying vertically at the top and bottom of a given line, which however does differ in variance by a constant value at all frequencies, the larger values at the bottom, as expected for a single top-buoy mooring line forming a J-shape. The spectra for horizontal tilt-difference signals are indeed less featureless, and show considerable variation in variance at all frequency ranges, including the sub-mesoscale range. Nevertheless the mesoscale bulge and the spectral slope remain statistically significant features.

Three-month mean spectra for T-sensor tilt-angle magnitude and top-sensor tilt differences

Figure 5. Three-month mean spectra for T-sensor tilt-angle magnitude of line 57 top (black) and bottom (cyan), and top-sensor tilt difference between line 57 and lines 12 (red), 14 (blue), 17 (green) and 68 (magenta). The spectra are a factor of two less smoothed than in Figure 4, focusing on the entire subinertial and IGW bands.

Comparing tilt and vorticity at sub-inertial frequencies

Although to first order mooring-line tilt-angle magnitude is imposed by ∣U∣2|\mathrm{U}|^2 due to the drag force, without considering a particular direction, a tight constant relationship as in (3) is not achieved in practice. This is partially due to the mooring construction that deviates from a stiff pendulum, but probably more due to flow variations. While vertical flow variations are unknown, because not measured, at the deep Mediterranean site, horizontal variations have been captured in ζ\zeta.

The observed small-scale ∣ζ∣=O⁡(f)|\zeta|=\operatorname{O}(f) are very large compared to reported open-ocean values of <0.2f<0.2f ([11]), but they are commensurate with, albeit larger than, mid-depth Mediterranean drifter observations of about $0.5f ([22]). Associated with peak-values of |\zeta|=4fareare|\Delta U| =∣ΔV∣=0.02 m s−1=|\Delta V|=0.02\ \mathrm{m}\ \mathrm{s}^{-1} over horizontal scales of ∣Δx∣|\Delta x|, ∣Δy∣=50 m|\Delta y|=50\ \mathrm{m}, i.e. a flow difference of 0.004 m s−10.004\ \mathrm{m}\ \mathrm{s}^{-1} over 10 m10\ \mathrm{m}, at sub-inertial scales. Using an aspect ratio of 1, this would yield a shear magnitude ∣S∣=4f|S|=4f that complies with observed small-scale stratification layering so that N=4fN=4f, in marginal stability with gradient Richardson number Ri=N2/∣S∣2≈1\mathrm{Ri}=N^{2}/|S|^{2}\approx1. As a result, the relative vorticity magnitude potentially reflects vertical flow differences that seriously affect mooring-line tilt-angle magnitude in addition to large-scale flow speeds.

Investigation of a one-month detail in Figure 6 shows the ambiguity of relationship between ∣α∣|\alpha| and ζ\zeta. It changes from out-of-phase between days 320–330 to in-phase between days 340–350. A more consistent relationship is found between ∣α∣|\alpha| and ∣ζ∣|\zeta|, which implies that the sign of relative vorticity does not matter. In general, correlation is good, although fails over short periods like around day 330 when relatively intense turbulence ([29]) and warm waters were observed (Figure 3b). In comparison for this monthly period, a positive relationship is found between ∣α∣|\alpha| and ∣U∣|U| for ∣α∣≲0.13∘|\alpha|\lesssim0.13^\circ (the dotted line in Figure 6), and a negative relationship for larger ∣α∣|\alpha|.

Examining the relationship between tilt and relative vorticity

Figure 6. Examining the relationship between tilt and relative vorticity. (a) One month time-series with scale of relative tilt-angle magnitude, from current meter at line 57, given to the right. For reference, ∣U∣|U| in green (scale far left). (b) Relative vorticity scaled with ff versus tilt-angle magnitude for first 10 days of a. when ζ<0\zeta< 0 and ∣α∣∼−ζ|\alpha| \sim-\zeta (dotted red), while average flow U>0U > 0 is eastward (Figure 3a). In background-grey, the entire three months of data. (c) As b., but for middle 10 days. (d) As b., but for last 10 days when ζ>0\zeta> 0 and ∣α∣∼ζ|\alpha| \sim\zeta, while U<0U < 0 (Figure 3a).

Details of a variable warming period

A six-day detail governed by easterly waterflow, slantwise warming from above/sideways, including relatively large ∣α∣|\alpha| and ∣ζ∣|\zeta| is investigated in more detail in Figs 7–9. In Figure 7 time series are shown of waterflow, relative tilt-angle magnitude, negative relative vorticity, and temperature with depth.

Six-day timeseries of waterflow and tilt variations at the top of the mooring array, in comparison with vertical temperature distribution

Figure 7. Six-day timeseries of waterflow and tilt variations at the top of the mooring array, in comparison with vertical temperature distribution. All data are noise-filtered unless indicated. (a) Waterflow speed from lines 57 (blue; 1-cpd low-pass filtered black), 35 (cyan) and 14 (magenta). (b) Relative tilt from line 57 using current meter data (blue; 1-cpd low-pass filtered black) and using top T-sensor (green; 1-cpd low-pass filtered yellow). (c) Relative sub-inertial vorticity scaled with inertial frequency ff. (d) Time-depth series of Conservative Temperature from line 25. Black contours are drawn every 0.001 ∘C0.001\,^\circ\mathrm{C}.

When slantwise warming is still relatively weak up to day 327, waterflow speeds are often consistent between the three current meters (fig:7a), except that the (east-component of the) inertial/IGW motion of line 14 wanders in- and out-of-phase with one of the other two. After day 327 more differences occur between the current meters 50 m apart horizontally with typical values of 0.010.01–0.02 m s−10.02\ \mathrm{m}\ \mathrm{s}^{-1}, also when the total waterflow speed is about 0.02 m s−10.02\ \mathrm{m}\ \mathrm{s}^{-1}, at day 330, when general relative mooring-line tilt is largest (fig:7b).

Thus, tilt is not only dependent on local waterflow speed, but, occasionally, also on that varying elsewhere. Tilt measured by current meter and upper T-sensor are generally found consistent, also at near-inertial scales. Some deviations occur at super-inertial scales, because the T-sensor has better tilt resolution, and also at sub-inertial scales at which the T-sensor tilt is smaller than the current meter’s by 0.010.01–0.02∘0.02^\circ when ∣U∣>0.04 m s−1|U| > 0.04\ \mathrm{m\,s^{-1}}, and larger by about the same amount when ∣U∣<0.04 m s−1|U| < 0.04\ \mathrm{m\,s^{-1}}.

Only during the first 3 days of the record, sub-mesoscale tilt follows negative relative vorticity, with a >1>1 day delay towards the end of the six-day period (Figure 7c). On the other hand, a consistent correspondence exists between temperature time-depth variations (Figure 7d) and tilt, with largest, turbulent warm-water pulses following tilt by a few hours up to one inertial period. This is confirmed by the movies associated with Figs 8 and 9, which also show the small-scale variability between different lines.

In Figure 8, the mooring array is viewed from above in the lower panel and the time-line is visible in the upper panel. It shows relative tilt from 10-cpd noise-filtered T-sensor data from 40 upper-level instruments for the period in Figure 7. The 72-s movie is accelerated so that 2 hours of real-time appear in 1-s movie-time. The associated movie starts with hardly any tilt above base-value, and very little variation between the lines. Soon though some light variations are visible, starting at the southwestern side, and notably appearing near the edges of the large ring at about 50-m horizontal distances. This pattern is followed for three consecutive near-inertial periods before intensifying and shorter-scale varying down to 9.5 m. Most intense values and variations occur around the entire ring-interior, around day 330. The movie ends like it started, although some emphasis is now seen on the eastern side of the ring.

Still from 2D movie showing low-pass filtered tilt data and sensor array

Figure 8. Still from 2D movie of 10-cpd low-pass filtered tilt data from 40 upper T-sensors during the period in Figure 7. In the lower panel the array is viewed from above, and each sensor is represented by a small filled circle of which the colour represents relative tilt magnitude in the scale to the right. In the movie’s upper panel, a white time-line progresses in a 6-d timeseries of tilt from all 40 lines. The 72-s movie is accelerated by a factor of 7200 with respect to real-time.

In Figure 9, a quasi-3D cube presents the temperature sensor array in the lower panel, with its colour scheme next to the upper-panel time line. The small panel to the right shows the upper-level flow speed by stick diagram for the three current meters. Waterflow speed is noise filtered at 4 cpd, whilst temperature at 3000 cpd. The associated movie has exactly the same duration and acceleration as that of Figure 8. Nevertheless it appears must faster and more vigorous, with frequent slashes of warm water from above. If the movie would be made slower (e.g. for half-day episodes in [29]), it could be seen that the warming also comes sideways. The waterflow is consistently to the east, but varies considerably with time and between the current meters. All flickering in temperature reflects turbulent motions, as noise has been filtered out, partially similar in intensity at all lines, but after more detailed inspection also varying between lines. Most intense is the period between days 330.0 and 331.0, about half an inertial period after the most intense tilt period of Figure 8.

Still from quasi-3D movie of 3000-cpd low-pass filtered temperature data

Figure 9. Still from quasi-3D movie of 3000-cpd low-pass filtered temperature data from about 2800 T-sensors under stratified water conditions during the period in Figure 7. In the lower panel’s cube, which is vertically depressed by a factor of two relative to horizontal scales, each sensor is represented by a small filled circle of which the colour represents Conservative Temperature in the scale above. In the movie’s upper panel, a white time-line progresses in a 6-d/124-m time/depth image from line 24 on the east-side of the cube. To the right of the cube, the noise-filtered waterflow from the three current meters is indicated by vector sticks. The 72-s movie is accelerated by a factor of 7200 with respect to real-time, like Figure 8.

Figure 10 shows small-scale variations develop over the 3D large-ring mooring array during a 24-minutes detail period. The movie-time is still 20 times faster than real-time, and reflects the slow but persistent character of deep-sea turbulence (N. Oakey, pers. comm.).

Temperature field and three-dimensional measurement detail

Figure 10. Still from a 24-minutes long slow-motion detail movie, still 20 times faster than real-time.

Discussion

Waterflows are best registered by dedicated instrumentation like current meters. If such rather costly instrumentation is not available, one may rely on proxies from instruments like mooring-line tilt sensor.

To first order, tilt-angle magnitude is proportional to water-flow speed squared. Such a proxy works best when the flow is uniform in space and purely horizontal, perpendicular to the direction of buoyancy (-gravity). However, waterflows vary over many scales and are thus not uniform also not in the deep sea, as is demonstrated in this paper.

Like any proxy, the effectiveness of tilt mimicking waterflows depends on the instrumental setup.

Under uniform flow, a mooring line with buoyancy evenly distributed over its entire length will behave closely like a stiff pendulum, with identical tilt along the line, whereas a single top-buoy mooring line has larger tilt near the anchor. The latter-case J-shape mooring effect is shown here for tilt-angle magnitudes as small as <0.5∘<0.5^\circ. It is therefore not surprising that tilt is observed only partially corresponding with flow-speed squared, although it remains to be investigated why such correspondence is found for small tilt only and even though direct measurements of vertical shear are missing.

Indirectly, the observed relatively large 50-m horizontal flow differences, having values of up to 100% of the mean flow value, demonstrate the presence of high flow variability at the deep Mediterranean site. These differences, and the associated relative vorticity magnitude, not only affect tilt, with a good correspondence with tilt variations over scales O(10)O(10) m, but also the energy transport to motions at other scales and frequencies, varying from 2D subinertial eddies to 3D turbulence, which is important for deep-sea life.

The lack of correspondence between tilt and relative vorticity suggests a reversal of cause and effect.

Instead of anticipated relative vorticity widening (for negative values) or shifting to higher frequencies (for positive values) of the IGW band, so that for instance near-inertial motions may be trapped or reflected thereby causing variations in shear and stratification, the actual process may be in reversed order. This is brought about by the observation that tilt relates with relative-vorticity magnitude, and probably shear, and precedes convection from above. Slantwise convection turbulence plumes may increase local vertical stratification, which becomes in balance with shear magnitude (of horizontal flow) to the level of marginal stability. This is diagnosed by relative-vorticity magnitude and, as a proxy, tilt magnitude, which reflects the arrival of an ‘underwater squall’ by its small-scale motion variations. Paradoxically, in the slantwise direction stratification may become near-homogeneous N∼0.3fN \sim0.3f, so that the IGW bounds may extend their subinertial range thereby explaining some of the observed tilt and tilt-difference spectral peaks around 0.5–0.7 cpd independently from the sign of relative vorticity, indirectly. As the convection turbulence is mainly driven by mesoscale variations, this also indirectly relates tilt with atmospheric disturbances [29].

While the observations have been made in weak flows and weak stratification, it may be worthwhile to test tilt as a proxy for strong water flows O⁡(0.1\operatorname{O}(0.1–1)1) m s−1^{-1} under more stratified conditions, so that convection turbulence may be less dominant. It would be especially interesting to see the variability in tilt between closely spaced mooring lines, whether these would be limited to O⁡(10)\operatorname{O}(10) m scales, as presumed for turbulence in the density-stratified vertical [6]. The present observations of weak-flow tilt suggest such a limit under the largest observed convection turbulence with variations O⁡(1)\operatorname{O}(1) times mean tilt, and which may also play a role in a forest under wind forcing.

Conclusions

Tightly tensioned, single-buoyancy mooring lines tilt <0.5∘<0.5^\circ under weak <0.07<0.07 m s−1^{-1} flow speeds at a site in deep Northwestern Mediterranean Sea. The tilt-angle magnitude can be used as a proxy for flow speed squared, but to first order only when tilt <0.13∘<0.13^\circ. Given up to half-inertial period delay from relative vorticity magnitude, mooring-line tilt seems a reasonably proxy for warm-water convection turbulence.

Horizontal, and presumed vertical, flow speed variations O(0.01)O(0.01) m s−1^{-1} cause tilt-angle magnitude variations O(0.1)∘O(0.1)^\circ over 50-m scales during weak turbulence and over O(10)O(10) m scales during more turbulent episodes, as is demonstrated in short movies.

The observed 50-m scale relative vorticity does not relate with tilt and hence does not affect the effective local vorticity by enlarging the IGW bounds.

The proxy for convection turbulence is also found in very high frequency tilt variations, which appear as quasi-random white noise in spectra. Mooring-line vibrations thus extend across a frequency range from subinertial O(0.1–1)O(0.1\text{--}1) cpd to Strouhal O(105)O(10^5) cpd, which may reflect energy transfer in an unknown manner.

Conflict of Interest

The author declares no conflict of interest relevant to this study.

Data availability

The movies to figures 8–10 can be found in [29], “Movies to: Short spatial mooring-tilt variations from deep Mediterranean observations”, Mendeley Data, V1, https://doi.org/10.17632/9bg7z7gz7p.1. Current meter data are available from [?], “Large-ring mooring current meter and CTD data”, Mendeley Data, V1, https://doi.org/10.17632/f8kfwcvtdn.1.

Acknowledgements Captains and crews of R/V Pelagia are thanked for the very pleasant cooperation. I also thank the team of ROV Holland I for the well-performed underwater mission to recover the instrumentation of the large ring. NIOZ colleagues notably from NMF department are thanked for their indispensable contributions during the long preparatory and construction phases to make the unique sea-operation successful. I highly appreciated working with colleagues within the KM3NeT collaboration. I has been a pleasure and a privilege. The author acknowledges the financial support of Nederlandse organisatie voor Wetenschappelijk Onderzoek (NWO), the Netherlands.

Determination of relative tilt-angle magnitude

Upper T-sensor tilt-angle magnitude can be referenced at the three vertical lines holding a current meter at the buoy (Figure A1), with the notion that the different tilt measurements may slightly differ because the buoy is swiveled and has a shackle-connection with the mooring-line. Also, no exact absolute reference was made after the mounting of the current meters to the buoy, so that a slight non-zero tilt-angle magnitude of about 0.05–0.1∘0.05\text{--}0.1^\circ was accepted.

Test for reference of line 14 upper T-sensor tilt measurements

Figure A1. Test for reference of line 14 upper T-sensor tilt measurements, noise filtered, as compared with current meter tilt data, noise filtered (yellow) and sub-inertial low-pass filtered (black). T-sensor reference dates with small waterflow <0.01 m s−1< 0.01\ \mathrm{m}\ \mathrm{s}^{-1} are: days 305.25 (blue record), 316.5 (red) and 317.7038 (green). (a) Detrended noise-filtered tilt-angle amplitude relative to an arbitrary offset. (b) As a., but for detrended tilt difference between lines 57 and 14.

Because of the larger electronic drift in T-sensor tilt compared to current-meter tilt, its data are first referenced to the mean value over a short period in time when the waterflow speed ∣U∣<0.01 m s−1|U| < 0.01\ \mathrm{m}\,\mathrm{s}^{-1}, so that the mooring line is expected to be close to vertical. (Unfortunately, no waterflow measurements were made at other levels than at h=125.5 mh = 125.5\ \mathrm{m} above seafloor. Thus, non-negligible waterflow speeds closer to the seafloor attributing to vertical current shear cannot be ruled out during these short periods.)

Three short periods from the start of the record are identified, on days 305.25, 316.5 and 317.7038. These tune the T-sensor tilt by about ±0.02∘\pm0.02^\circ.

In general after fitting a non-zero reference base-value of 0.07∘0.07^\circ and choosing the day-317 reference, a reasonable comparison is obtained between current-meter and T-sensor tilt at line 14. All excursions up from base-value are identified in both types of instrumentation including near-inertial periodic and sub-mesoscale ones (fig:A1a), although the amplitude of T-sensor tilt in general exceeds that of current-meter tilt during largest excursions such as around day 330, and vice-versa during smaller excursions such as around day 380. It is not clear what causes these different tilt-angle amplitudes between the instrumentation, but it is likely due to the non-stiff connection, and about 1 m different vertical position, between buoy and vertical line and thus a different response upon forcing. The shear may also be larger during largest excursions, affecting the line’s tilt more than that of the buoy. Such a different response is also found across horizontal distance of about 50 m, with tilt-angle differences ±0.1∘\pm0.1^\circ, well exceeding noise levels by more than one order of magnitude, and with similar magnitude for the two instrumentation systems, albeit not providing identical variations.

Filters and noise

Instrumental noise level for T-sensor tilt-angle magnitude is reached at about 10 cpd (Figure A2), which gives half an order of magnitude larger range than for current-meter tilt (Figure 4). Noise-filtering thus implies low-pass filtering with presented cut-off at 10 cpd (Figure A2), and 4 cpd for current-meter tilt. Sub-inertial filtering is achieved by applying a low-pass filter with cut-off close to 1 cpd (Figure A2), which effectively removes most of inertio-gravity waves and retains daily and longer motions, including at (sub-)mesoscales.

The flatness of the T-sensor tilt spectrum in Figure A2 for ω>10\omega> 10 cpd suggests random white noise, or a pure instrumental-determined cause. However, it may also contain a quasi-random response of the instrumentation to environmental changes, as mooring-line tilt is strongly dependent on flow magnitude that causes line-oscillations. This is visualized in Figure A3, in which the three-month record of tilt is inspected in terms of its very high-1000-cpd-pass filtered data. In fig:A3a, the full record’s standard-deviation is shown, which mimics all of the tilt variations in fig:3c. Dissemination of tilt from three episodes varying from calm, around day 320 (fig:A3b), via moderately turbulent around day 325 (fig:A3c), to relatively intensely turbulent around day 329 (fig:A3d), demonstrates a non-continuous amplitude variation with similar intermittency albeit at different levels, for all of these. The 1000–1800-cpd spectra show half–one order of magnitude change in variance, besides the general horizontal flatness, mimicking white noise well within error, and yet some, insignificant, variation with frequency, such as a common local high around 1300 cpd and local lows around 1100 and 1800 cpd, for all three episodes. The response of the instrumental setup thus qualitatively reflects the convection turbulence at all scales of tilt-angle magnitude.

References

Mooring array location and layout

Figure 1. Mooring array location and layout. (a) Location named "Temp Array" (orange dot) on map off southern France. The mooring is well east of main neutrino-telescope site “ORCA” of KM3NeT [1] and just northeast of the former “ANTARES” neutrino-telescope site. Isobaths are drawn every 100 m. The grey lines indicate Toulon-harbour approach sectors. (b) Mooring-array orientation at seabed and layout, with steel-cable grid and small rings holding the vertical lines at 9.5 m intervals. Lines are numbered in six synchronization groups. Single synchronizer S is at ring 51. Waterflow ‘WF’ instruments are at buoys on top of lines 14 (indicated in magenta ellipse), 35 (cyan) and 57 (blue). The insert shows part of the large ring just prior to deployment at sea, with free-fall drag-parachute in front.

References

  1. [1]Adrián-Martinez, S., Ageron, M., Aharonian, F., Aiello, S., Albert, A., Ameli, F., et al. (2016). Letter of intent for KM3NeT 2.0. Journal of Physics G, 43(8), 084001. https://doi.org/10.1088/0954-3899/43/8/084001
  2. [2]Bak, P., Tang, C., & Wiesenfeld, K. (1987). Self-organized criticality: An explanation of the 1/f noise. Physical Review Letters, 59, 381-384. https://doi.org/10.1103/PhysRevLett.59.381
  3. [3]Bolgiano, R. (1959). Turbulent spectra in a stably stratified atmosphere. Journal of Geophysical Research, 64, 2226-2229. https://doi.org/10.1029/JZ064i012p02226
  4. [4]Crepon, M., Wald, L., & Monget, J. M. (1982). Low-frequency waves in the Ligurian Sea during December 1977. Journal of Geophysical Research, 87, 595-600. https://doi.org/10.1029/jc087ic01p00595
  5. [5]Dewey, R. K. (1999). Mooring design & dynamics—a Matlab® package for designing and analyzing oceanographic moorings. Marine Models, 1, 103-157. https://doi.org/10.1016/S1369-9350(00)00002-X
  6. [6]Gargett, A. E., Hendricks, P. J., Sanford, T. B., Osborn, T. R. & Williams, A.J. (1981). A composite spectrum of vertical shear in the upper ocean. Journal of Physical Oceanography, 11, 1258-1271. https://doi.org/10.1175/1520-0485(1981)011<1258:ACSOVS>2.0.CO;2
  7. [7]Garrett, C., & Munk, W. (1972). Space-time scales of internal waves. Geophysical Fluid Dynamics, 3, 225-264. https://doi.org/10.1080/03091927208236082
  8. [8]Gerkema, T., Zimmerman, J. T. F., Maas, L. R. M., & van Haren, H. (2008). Geophysical and astrophysical fluid dynamics beyond the traditional approximation. Reviews of Geophysics, 46, RG2004. doi:10.1029/2006RG000220DOI
  9. [9]IOC, SCOR, & IAPSO (2010). The international thermodynamic equation of seawater—2010: Calculation and use of thermodynamic properties. In Intergovernmental Oceanographic Commission, Manuals and Guides No. 56. UNESCO.
  10. [10]Kolmogorov, A. N. (1941). The local structure of turbulence in incompressible viscous fluid for very large Reynolds numbers. Doklady Akademii nauk SSSR, 30, 301-305.
  11. [11]Kunze, E. (1985). Near-inertial wave propagation in geostrophic shear. Journal of Physical Oceanography, 15(5), 544–565. https://doi.org/10.1175/1520-0485(1985)015<0544:niwpig>2.0.co;2
  12. [12]LeBlond, P. H., & Mysak, L. A. (1978). Waves in the ocean. Elsevier.
  13. [13]Marshall, J., & Schott, F. (1999). Open-ocean convection: Observations, theory, and models. Reviews of Geophysics, 37, 1-64. https://doi.org/10.1029/98rg02739
  14. [14]Moller, D. A. (1976). A Computer Program for the design and static analysis of single-point subsurface mooring systems: NOYFB. WHOI-76-59. Woods Hole, MA: Woods Hole Oceanographic Institution, 106 pp. https://darchive.mblwhoilibrary.org/server/api/core/bitstreams/0f41541c-7db6-5641-8412-02f68276b439/contentDOI
  15. [15]Obukhov, A. M. (1949). Structure of the temperature field in a turbulent flow. Izvestia Akademii nauk SSSR, Seriya Geograficheskaya i Geofizika, 13, 58-69.
  16. [16]Obukhov, A. M. (1959). Effect of buoyancy forces on the structure of temperature field in a turbulent flow. Doklady Akademii nauk SSSR, 125, 1246-1248.
  17. [17]Phillips, O. M. (1971). On spectra measured in an undulating layered medium. Journal of Physical Oceanography, 1, 1-6. https://doi.org/10.1175/1520-0485(1971)001<0001:OSMIAU>2.0.CO;2
  18. [18]Reid, R. O. (1971). A special case of Phillips’ general theory of sampling statistics for a layered medium. Journal of Physical Oceanography, 1, 61-62.
  19. [19]Schuster, H. G. (1984). Deterministic chaos: An introduction. Physik-Verlag.
  20. [20]Sheremet, V. A. (2004). Laboratory experiments with tilted convective plumes on a centrifuge: A finite angle between the buoyancy force and the axis of rotation. Journal of Fluid Mechanics, 506, 217-244. doi:10.1017/S0022112004008572DOI
  21. [21]Straneo, F., Kawase, M., & Riser, S. C. (2002). Idealized models of slantwise convection in a baroclinic flow. Journal of Physical Oceanography, 32, 558-572. https://doi.org/10.1175/1520-0485(2002)032<0558:imosci>2.0.co;2
  22. [22]Testor, P., & Gascard, J.-C. (2006). Post-convection spreading phase in the Northwestern Mediterranean Sea. Deep-Sea Research I, 53, 869-893. https://doi.org/10.1016/j.dsr.2006.02.004
  23. [23]van Haren, H. (1996). Correcting false thermistor string data. Continental Shelf Research, 16, 1871-1883. https://doi.org/10.1016/0278-4343(96)00020-9
  24. [25]van Haren, H., (2026). Technical note: Determination of the height of deep-sea mooring lines above seafloor using turbulence measurements. Ocean Science, 22, 1501-1513. https://doi.org/10.5194/os-22-1501-2026
  25. [26]van Haren, H., & Millot, C. (2009). Slantwise convection: a candidate for homogenization of deep newly formed dense waters. Geophysical Research Letters, 36, L12604. https://doi.org/10.1029/2009GL038736
  26. [27]van Haren, H., Maas, L., Zimmerman, J. T. F., Ridderinkhof, H. & Malschaert, H. (1999). Strong inertial currents and marginal internal wave stability in the central North Sea. Geophysical Research Letters, 26, 2993-2996. https://doi.org/10.1029/1999GL002352
  27. [28]van Haren, H., Bakker, R., Witte, Y., Laan, M., & van Heerwaarden, J. (2021). Half a cubic hectometer mooring array 3D-T of 3000 temperature sensors in the deep sea. Journal of Atmospheric and Oceanic Technology, 38, 1585-1597. https://doi.org/10.1175/JTECH-D-21-0045.1
  28. [29]van Haren, H., Adriani, O., Albert, A., Alhebsi, A.R., Alshalloudi, S., et al. (2026). Whipped and mixed warm clouds in the deep sea. Geophysical Research Letters, 53, e2025GL119998. https://doi.org/10.1029/2025GL119998

Paper details

vlnc:W_ynurkdfv
v1 · latestHistory (1)

published ·

added ·

Usage

views · 0

downloads · 0

A view counts once per browser per day, a download each time the PDF is requested.

CC BY 4.0

Original source

  • Physics and Astronomy

Not recorded for this paper · View statement →

Looking for related papers…
Contents