Next Article in Journal
Spatio-Temporal Hydrological Model Structure and Parametrization Analysis
Next Article in Special Issue
Essential Ocean Variables for Marine Environment Monitoring: Metrological Case Studies
Previous Article in Journal
How Good Is the STW Sensor? An Account from a Larger Shipping Company
Previous Article in Special Issue
Sea State Monitoring by Ship Motion Measurements Onboard a Research Ship in the Antarctic Waters
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Improved Current Estimates from Spar Buoy-Mounted ADCP Measurement Station: A Case Study in the Ligurian Sea

1
ENEA-Marine Environment Research Centre of S. Teresa, 19032 Pozzuolo di Lerici, Italy
2
CNR-IAS, Via De Marini, 6-16149 Genova, Italy
3
INRiM-Istituto Nazionale di Ricerca Metrologica, Strada delle Cacce 91, 10135 Torino, Italy
4
IIM-Istituto Idrografico della Marina Militare, Passo dell’Osservatorio 4, 16134 Genova, Italy
*
Author to whom correspondence should be addressed.
J. Mar. Sci. Eng. 2021, 9(5), 466; https://doi.org/10.3390/jmse9050466
Submission received: 30 March 2021 / Revised: 22 April 2021 / Accepted: 23 April 2021 / Published: 26 April 2021
(This article belongs to the Special Issue Marine Metrology and Oceanographic Measurements 2020)

Abstract

:
Current measurements in the open sea are generally acquired by Acoustic Doppler Current Profilers (ADCPs). In the case of ADCPs mounted on spar buoy, current profiles require to be post-processed, to properly take into account the buoy influence: in fact, ADCP compass may reflect alterations induced by the metal structure of the buoy and apparent currents can occur due to the large displacement of the platform. Uncertainty analysis is finally required to properly consider both these effects and to compute robust velocity estimates. A new methodology is tested for a measurement station in the Ligurian Sea, where an ADCP was mounted on the surface buoy of the W1-M3A (Western 1 Mediterranean Moored Multisensor Array) oceanographic observatory, facing upwards at the depth of about 40 m. Marine current numerical models and historical data in the area have been used as a basis for comparison to test the consistency of the proposed method. A very good agreement is obtained. Only minor discrepancies are reported (e.g., monthly averages from the reference model slightly underestimate the west-east current component along the entire profile), but, in general, the application of the proposed methodology ensures that the spar buoy-mounted ADCP system is able to provide reliable measurements for oceanographic studies and validation of 3D hydrodynamic models.

1. Introduction

Current data are a major asset for oceanographic research both for validation and assimilation in hydrodynamic models. Current monitoring is also an essential tool for understanding ocean marine ecosystems, planning and securing maritime operations or designing offshore structures. The measurement of current is the basis for any resource assessment related to offshore renewable energy devices. Accurate current data can help in tracking sources of pollution, like chemical or sewage spills. Same data can be used to efficiently and reliably forecast oil spill trajectories or assess coastal water quality in general. Reliable information on current conditions and consequent monitoring are also required for habitat restoration or coastal protection projects, especially considering the future spread of nature-based solutions (NBS).
Many efforts have been made in the last decades to monitor ocean currents using Lagrangian and Eulerian platforms. Classical Lagrangian platforms (i.e., drifters and floats) could infer current measurements over a wide area only through indirect methods, although Argo floats have been recently upgraded to host high frequency acoustic current meters. Modern autonomous moving platforms (i.e., gliders and other types of marine drones) can further increase spatial and temporal sampling of ocean current since most of them can submerge for long periods to get measurements at a particular depth, but cannot provide long term time-series at fixed point. High frequency coastal radar can ensure continuous monitoring of surface ocean currents at proper spatial and time scale (range up to 200 km and time intervals typically of 1 h), but cannot provide information about the current profile along the water column. Eulerian platforms (i.e., buoys and moorings) may collect in-situ current observations at a high temporal resolution along the water column, if their scientific payload includes array of current meters or Acoustic Doppler Current Profilers (ADCPs).
The issue of temporal and depth coverage is commonly tackled by ADCPs. They are based on the principle of Doppler Shift: they emit an acoustic signal and receive it backscattered by the particles moving at various depths in the water column due to current action [1]. The profile of speed and direction of the current at different depths is then calculated from the frequency shift of the backscattered signal. ADCPs are normally deployed on the seafloor, fixed along mooring lines, or mounted inside the keel of a boat (vessel-mounted), but can also be mounted on spar buoy. In that case, current profiles require to be post-processed to properly model the buoy influence. The following two effects occur, both attributable to the buoy action on the ADCP:
  • alterations of the ADCP compass due to the local distortion of the Earth’s magnetic field induced by the metal structure of the buoy: such distortion is not constant over time due to the rotation of the buoy itself caused by winds, waves and currents;
  • apparent currents measured by ADCP due to displacements of the buoy following the action of currents and winds: in case of long mooring line, the buoy is in fact free to move around the anchor.
An ad hoc uncertainty analysis is then required to properly consider both these effects and to compute robust velocity estimates. As regards metal body interference, Munchow et al. [2] tested a submerged towed vehicle with both an upward- and a downward-facing ADCP, designed to have a minor exposure to surface disturbance and wave wake. They stressed for the first time the need for a careful in-situ calibration of magnetic compass, to remove the magnetic variation induced by the ship’s hull and the tow platform and to improve the accuracy of the data. The correction of errors on direction measurements brought by the current profilers installed on mooring cages is introduced by Le Menn et al. [3], who developed a calibration platform for magnetic compasses and tilt sensors, following a method perfected in 2007. In combination with this platform, [4] suggested an approach based on the measurement of the frequency of pulses and on the simulation of received echoes, allowing the logging and calibration of most of the current meters manufactured by Nortek Instruments.
The main focus of the present work is illustrating the details of an experimental activity aimed at measuring current profiles in open sea using an ADCP mounted on a spar buoy: the proposed methodology, from the set-up of the experiment to its implementation, from recovery to data post-processing is described in the following.
The paper is organized as follows: a brief description of both the W1-M3A (Western 1 Mediterranean Moored Multisensor Array) oceanographic infrastructure and of the used ADCP is illustrated in Section 2. In Section 3, the sea-current dataset is analyzed as an applicative case study of the post-processing correction methodology, that is discussed in detail: actually, the application of the method must be understood as a result, rather than the interpretation of the data from an oceanographic point of view. Discussion is covered by Section 4, where comparisons with historical data and numerical models in the area are presented to test the consistency of the proposed method, together with its uncertainty analysis performed as an example on a sub-set of data. Conclusions, together with directions of future work, are drawn in Section 5. The collected dataset, together with the proposed post-processing techniques, allowed to achieve a basic knowledge, aimed at implementing a monitoring system able to operate reliable ADCP measurements even from this kind of moving platforms, or from profiling floats equipped with ADCPs ([5]).

2. Experimental Setting: The W1-M3A Spar Buoy and the ADCP

The W1-M3A is located about 40 NM offshore Genoa, on a 1270 m deep sea bed [6,7] (43°48.90′ N, 09°06.80′ E) (Figure 1a). The observatory is composed of two sub-systems, a large spar buoy (known as “ODAS Italia 1”, 51 m long of which about 37 m immersed and 12 tons weight, shown in Figure 1b: ODAS, or more simply buoy, in the following) and a sub-surface mooring acquiring oceanographic data down to the ocean interior. The Buoy is a completely autonomous laboratory that is part of the of European Multidisciplinary Seafloor and water-column Observatory and of the Integrated Carbon Observation System-European Research Infrastructure Consortiums (the EMSO and ICOS-ERIC respectively). Its power supply is guaranteed using solar and wind energy, while the communication is ensured through satellite connection. Monitoring the meteorological conditions of the basin and the physical and biogeochemical status of the area, the Buoy is a useful support for operational oceanography and for meteorology centers, contributing to weather and marine modeling and forecasting, and allowing us to develop comparisons between the observations from satellites and in situ measurements.
A three-beam ADCP, manufactured by Nortek (model Aquapro, [8]), was selected in order to avoid the possible interference of a fourth beam with the structure of the Buoy. The general characteristics of the instrument are reported hereafter:
  • acoustic frequency: 400 kHz;
  • max working depth: 200 m;
  • max profiling range: 60–90 m;
  • cell size: 2–8 m;
  • beam width: 3.7°;
  • min blanking distance: 1 m;
  • max number of cells: 128;
  • velocity range: ±10 m s−1;
  • accuracy declared by the manufacturer: 1% of measured values ±0.005 m s−1;
  • max sampling rate: 1 Hz.
Before the installation on the buoy, the Nortek ADCP was verified against a four-beam RDI ADCP (Figure 2), at water depth of about 37 m, comparable with the damping disk depth of the Buoy. Same general settings and the ENU reference system (East-North-Up coordinates) were configured for both the ADCPs, with automatic correction of roll, pitch and heading. The two ADCPs were installed on the frame of a CTD probe, with reference beams facing the same direction (beam X for the Nortek, beam no. 3 for the RDI). Mean current measures of the two ADCPs, along the water column, resulted reasonably comparable inside their standard uncertainties, calculated as the standard deviation of the collected datasets (e.g., at the depth of 18 m, Nortek and RDI measured a mean north component of the current equal to (6.7 ± 1.5) cm s−1 and (5.6 ± 1.7) cm s−1, respectively).
The ADCP was then set to collect current signals at 1800 s intervals (average measure as a result of 900 s acquisition period, with 10 s compass update rate), in ten 4 m-long depth cells, with a (nominal) horizontal velocity precision of 0.006 m s−1. The ADCP mechanical configuration on the Buoy was properly designed to avoid any interferences of the acoustic beams with the Buoy structure: ADCP was placed on the damping disk of the Buoy at the depth of about 37 m and facing upwards (Figure 3). Specifically, the frame of the ADCP was clamped by divers to pre-existing bolts that were fixed to the damping disk when the buoy was in the dry dock. By this way, it was assured that any frame attached to those bolts was perfectly aligned to the bow of the buoy. Current measures were continuously collected for a period of 5 months (from April to August 2017, with a brief interruption due to battery change, for a total of 7256 current profiles).

3. Results

In the following, post-processing of ADCP data is described, step-by-step.

3.1. Estimated Average Depth of the ADCP during the Deployment Period

From the pressure measurements recorded on board of the ADCP, the average depth of the head of the acoustic sensor was derived, in order to correctly scale the depths associated with the centers of the various cells into which the current profile was divided. The procedure can be described in the following points:
  • average deployment pressure of the ADCP throughout the entire deployment period (and associated standard deviation):
p ¯ A D C P = ( 37.46 ± 0.20 )   dbar
  • conversion of pressure in depth, taking into account the following parameters:
    -
    average position of the buoy, in terms of latitude: 43.8430° N,
    -
    average temperature measured by the ADCP: 15.98 °C,
    -
    average salinity (from historical data, a priori set in the ADCP configuration): 38,
the average depth to which the ADCP was deployed is the following:
d ¯ A D C P = ( 37.16 ± 0.20 )   m
  • from the calibration certificate provided by Nortek, the pressure measurements on board the ADCP were declared with a tolerance of ±0.5% of the full scale (200 m). The contributions due to drift over time, resolution and uncertainty associated with the single pressure measurement were compassed in this range of variability. By combining the various contributions of standard uncertainty, the following measurement of the average depth of immersion of the ADCP was finally obtained:
d ¯ A D C P = ( 37.2 ± 0.6 )   m
  • the valid measuring range of the ADCP could be estimated from the well-known empirical formula [8] (where θ = 25° indicates the angle of the single beam with respect to the vertical):
Δ L A D C P = d ¯ A D C P × cos θ c e l l   s i z e   ~   30   m
The measurement range is therefore valid for depths between about 37 m and 7 m: considering that there is 1 m of blank, cell num. 8, whose center is located about 4 m deep, is the first of the critical cells in terms of accuracy (and will therefore be excluded from the analysis).
In conclusion, the actual depths of the cells constituting the current profiles measured by the ADCP are shown in Table 1.
As additional information, Figure 4 shows the trends over time of the pressure (and temperature) values recorded by the sensors on board the ADCP: pressure oscillations, thanks to the damping action of the Buoy disk, are always much lower than the bin size of 4 m.

3.2. Intensity Check on the Three Acoustic Beams of ADCP

To verify that the acoustic signals recorded by ADCP show no anomalies, the average values (vs time) of the counts on the three acoustic beams were compared, for each cell depth (Figure 5a).
It can be noticed how the three beams behaved in a homogeneous manner (within the respective standard deviations): in fact, acoustic signals correctly decreased towards the sea surface (due to an increase of the distance from the emitting/receiving sensor). Greater consistency was found for the cells in the range 2 to 7. For what concerns cells no. 10, 9 and 8 their exclusion from the analysis of the current profiles was confirmed: positions of cells no 10 and 9 were directly in air and close to the water–air interface, respectively. Cell no. 8 was characterized by an excessively high number of counts (if compared to deeper cells): this may be due again to the water–air interface effect.
Finally, the signal-to-noise ratio (SNR) of the echo signal was maintained at acceptable levels during the whole measurement period, for all the three acoustic beams (Figure 5b).

3.3. Buoy Heading and GPS Data Matched with ADCP Data: Correction for Heading Error and for the Velocity Vector Associated with the Buoy Displacement

In accordance with procedure described in detail in [9], to ensure the correct matching between ADCP and buoy data (acquired at half an hour and every hour, respectively) it is necessary to interleave the data of the latter, in terms of its heading and GPS longitude–latitude values. Velocity vector (in terms of modulus and bearing) associated with buoy displacement can then be calculated by applying the well-known spherical law of cosines formula [9,10].
First of all, since ADCP compass values could be potentially influenced by the presence of the metal body of the buoy, the ODAS heading measures were adopted as a reference and consequently the correction to be adopted on the ADCP measures is evaluated. The heading values for both the ODAS and the ADCP are shown in Figure 6 (where the value of 0° indicates the north direction, as per convention).
The analogy with the information obtainable from the GPS position of the buoy (Figure 7) can be noticed: the buoy, in fact, remained usually in a stationary position for most of the time in a direction approximately equal to the North-North-West, evidently put in traction by a current and/or by the action of winds coming from South-South-East. Being the mooring attachment point the mechanical reference of the buoy used to define its orientation (or heading), it can be noticed that the buoy heading values were those expected, that is, on average positioned in the 2nd quadrant (South-South-East, Figure 6).
At the same time, since the ADCP was mounted at the antipodes with respect to the buoy mooring, it can be observed that the ADCP heading values were reversed compared to those of the buoy, unless an error was to be evaluated and then corrected for each ADCP data acquired.
A mean evaluation of ADCP heading error can be obtained as follows: averaging the values shown in Figure 6, the average heading measured for the buoy and ADCP, respectively, is obtained in the same period. From the first of the two values, adding an angle of 180°, the theoretical average direction of the ADCP heading is then obtained. It can be seen (Figure 8) how this value differs from the average value measured by the buoy by about 14°: it can be concluded that this is on average the value of the phase shift induced by the metal body of the buoy on the compass of the ADCP.
Actually, the correction of the ADCP heading by means of buoy heading was carried out for each individual profile acquired: the difference between the two values (which is an estimate of the error made by the ADCP compass for the measurements taken every half hour) was subtracted from the various angles of the current velocity vectors cell by cell, profile by profile. As an example, Figure 9 shows the effect of the heading correction for the current measurements detected at cell no. 4 (in the middle of the useful profile, at a depth of about 20 m).
Finally, Figure 10 shows the overall trend of the ADCP compass error as a function of the buoy compass values, showing how this effect depended on the orientation of the buoy itself.
After the correction for heading error, ADCP data were ready to be corrected for the relative buoy movement. Since the buoy, in addition to rotation on its own axis, moved in the plane under the action of current, waves and/or winds, the ADCP was affected by apparent currents for which the measured current profiles had to be corrected. If the buoy moved with velocity v, the ADCP measured an equal and opposite apparent current -v: the correction for the ODAS movement could then be implemented by vectorially composing the buoy and ADCP velocity vectors cell by cell, profile by profile.
Figure 11 shows in polar coordinates an overview of all the velocity vectors measured by the ADCP in the seven considered cells (already corrected for the buoy heading), together with the velocity vectors of the buoy derived from GPS measurements.
Finally, Figure 12 shows the polar diagrams relating to all corrected currents values in the two cells at the extremes of the profile range (cells 1 and 7, at depths of about 32 m and 8 m, respectively), during the period considered: these were the final values of the current profiles, corrected both for the buoy heading and for the displacement of the buoy itself. The information contained in Figure 12 is part of final, corrected experimental dataset: it had to be intended as the best estimate provided by the implemented method.

4. Discussion

Based on the methodological approach described in Section 3, an overall assessment of data-set consistency is hereafter discussed. To evaluate the results of the applied corrections, an analysis of the obtained currents temporal evolution and general characteristics was performed. In the lack of other independent and concurrent ADCP measurements, as well as of a reliable climatology to assess the consistency of the data set, a comparison was done with available ADCP measurements collected in the area during the past years, and with results from numerical model.
General characteristic of the observed currents are resumed by the progressive vector (Figure 13) and basic statistics. Currents in the surface layer are directed to the North-West at an average speed of about 19 cm s−1 over the whole period of measurements, in accordance with the general circulation scheme of the Ligurian Sea [11,12] and previous ADCP observations from the same area [13]. Vertical structure is quite barotropic: currents have almost the same pattern at each level, with some attenuation of the speed with the depth.
Complex correlation coefficients are lower than expected (0.6–0.9), considering the high vertical resolution (each cell is 4 m) of the measurements and the short water column examined, probably resulting from a few cases having more than 15 cm s−1 difference in magnitude of single measurements between two adjacent cells.
Rotary spectral analysis was performed by means of FFT on 11 subsamples 340 h long, that were then averaged. The clockwise spectrum is dominated by the inertial period, which is 17.3 h at this latitude; a small relative maximum appears on 12 h period only in the clockwise part, while in the anticlockwise part there are no relevant signals, apart from a small peak at about 17 h. Inertial rotations are also well evident in the progressive vectors at all depths; this is well consistent with the Ligurian Sea circulation, where inertial movements account for the great part of sub-daily variability [14].
Existing ADCP data-set collected in the area during the same season [13,15], were recovered and re-analyzed. Period of measurements and setting are resumed in Table 2. The ADCP was a 300 kHz RDI looking up, located on a fixed mooring deployed close to the buoy (43°47.77’ N; 9°02.85’ E). As the depth of the ADCP changed from one deployment to the others, only the time series of currents collected at the bin including the layer at 22–24 m depth, which is a common depth for all the samples, were selected for the comparison.
The two time series of 2005 were merged after hourly averaging the time series sampled at 30 min. The resulting time series (bc) covers the correspondent period of the 2017, apart from 38 h of missing data which were replaced by a constant value in order to preserve the temporal evolution. Basic statistics are reported in Table 3.
Despite 2017 summer currents are higher in respect to 2005, the comparison of the temporal evolution of daily mean currents (Figure 14) reveals relevant common features. In particular, an important mesoscale variability, which characterizes the early spring and summer circulation of the area [11,16], is found in all the three examined samples. A comparison with rotary spectral analysis on both a and bc samples (Figure 15) confirms the dominance of this variability and shows quite the same amplitude of that of 2004.
After the comparison with historical in-situ measurements in the area, experimental data have been compared with available results from numerical simulation in the area, to further evaluate the reliability of the corrections applied by the post-processing method. The E.U. Mediterranean Copernicus Marine Environmental Monitoring Service (CMEMS, http://marine.copernicus.eu, accessed on the 1 March 2021 ) physical system has been selected [17], and in particular monthly current data mean from the Mediterranean physical reanalysis product [18]. It is composed of a hydrodynamic model, supplied by the Nucleous for European Modelling of the Ocean (NEMO), and of a variational data assimilation scheme (OceanVAR) for temperature and salinity vertical profiles and satellite Sea Level Anomaly along track data. The horizontal resolution is approximately 5 km and there are 141 unevenly spaced vertical levels. The grid-model closest to the position of the buoy has been identified, corresponding to a distance of about 3 km from the average position of the buoy.
In Figure 16 East and North mean monthly current profiles are reported, showing a reasonable agreement especially with regard to the north components. On the other hand, minor discrepancies appear in relation to the profiles of the east components, closer to positive values in the case of the CMEMS model. Values are however comparable within the limits of the dispersion bands (calculated as standard deviation of measure dispersions in the considered period). A strong barotropic behavior is evident for both components of the model, result of the application of the computational algorithms.
In order to evaluate the differences found between experimental and model data, we focus on a significant phenomenological event, to highlight whether there is congruence between CMEMS outputs at higher temporal resolution and experimental data. Hourly data from the model were then compared with the measurements carried out every 30 min by the ADCP at the end of June 2017, when a significant increase in the East current component and a simultaneous reduction in the North component occurred. Results are reported in Figure 17, showing how CMEMS forecast was able to correctly identify in space and time the eastward rotation of the currents along the water column. On the contrary, the model does not represent well the significant increase in the northern component, clearly noted by the ADCP around the 23 June 2017.
Comparisons with historical data and numerical models need to be supported by a proper uncertainty analysis (section C in [9]). If on one hand the standard uncertainty associated to the heading correction can be calculated by applying the law of uncertainty propagation [19], on the other the spherical law of cosines used to calculate the buoy apparent current is clearly a nonlinear model: therefore, due to correction for buoy apparent currents, standard uncertainties of both net current velocities and their directions need to be calculated in accordance to Monte Carlo method [20], as a by-product of the propagation of the normal probability distributions modeling buoy latitude and longitude values (i.e., the input quantities by which buoy bearing and displacement, and consequently its velocity, are derived).
As an applicative example, evaluation of standard uncertainties associated with corrected current velocities is finally reported in Figure 18, computed over the same reference period of Figure 17. This simulation was performed by R software [21], randomly generating M = 106 data for each input quantities.

5. Conclusions

In oceanography, current measurements mostly rely on current profilers, but metrological controls are still missing or poorly applied.
The procedure here described allows a proper uncertainty analysis to be applied on current meters profilers installed on spar buoys. Both the magnetic influence (exerted by the buoy body on the ADCP compass) and the displacement of the buoy itself (origin of apparent currents detectable by the ADCP) are taken into account to assess the measurement reliability. Both heading error and apparent currents can be actually corrected by measuring buoy heading and GPS position
One of the main conclusion is that a standard correction factor could not be applied for every spar buoy equipped with an ADCP, as the effects are case sensitive. The W1-M3A buoy showed to have its own behavior and errors are different for different buoy orientations. The only calibration of the compass is not sufficient to have reliable current measurements from a spar buoy as effects as apparent currents need to be taken into account (and their correction is not linear).
The application of this methodology is also important in order to define an appropriate metric for the comparison of in-situ and models data. The comparison with past ADCP measurements and concurrent model results indicates an overall consistency of the ADCP time series with the dynamics of the area, further supporting the good quality of the applied correction.
The improvement of current estimates can be used to properly validate hydrodynamical models and consequently verify their applicability for specific geographical area. It also allowed us to demonstrate a strong qualitative accordance, in time and space, between the experimental data, literature and hourly model data. The estimation of the performance of an ADCP on a mobile platform is also important in the light of future applications on Argo floats where ADCP can be embedded [5].

Author Contributions

Conceptualization and methodology, A.B., T.C., G.R. and A.S.; operational activities A.B., G.R., R.B. and S.P.; writing—original draft preparation, writing—review and editing, T.C., G.R. and P.P.; overall assessment of data-set consistency and uncertainty analysis, P.P. and F.P. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

The data presented in this study are post-processed data of the RIMA project (owner), used as test case to described the analysis method.

Acknowledgments

This work was partly developed within the framework of the Project RIMA (Development of technologies and softwares for an Integrated Forecasting Mediterranean Network as a management tool for the marine and coastal environmen) funded by the Italian Ministry of Education, University and Research. Maps have been realized by using Ocean Data View (Schlitzer, Reiner, Ocean Data View, odv.awi.de, 2020).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Emery, W.J.; Thomson, R.E. Data Analysis Methods in Physical Oceanography; Pergamon Elsevier Science: New York, NY, USA, 1997. [Google Scholar]
  2. Munchow, A.; Coughran, C.S.; Hendershott, M.C.; Winant, C.D. Performance and calibration of an acoustic Doppler current profiler towed below the surface. J. Atmos. Ocean. Technol. 1995, 12, 435–444. [Google Scholar] [CrossRef] [Green Version]
  3. Le Menn, M.; Pacaud, L. Calibration of currentmeters in direction: Results obtained on a stock of instruments with a new calibration platform. In Proceedings of the 17th International Congress of Metrology, Paris, France, 21–24 September 2015. [Google Scholar]
  4. Le Menn, M.; Morvan, S. Velocity Calibration of Doppler Current Profiler Transducers. J. Mar. Sci. Eng. 2020, 8, 847. [Google Scholar] [CrossRef]
  5. André, X.; Le Traon, P.-Y.; Le Reste, S.; Dutreuil, V.; Leymarie, E.; Malardé, D.; Marec, C.; Sagot, J.; Amice, M.; Babin, M.; et al. Preparing the New Phase of Argo: Technological Developments on Profiling Floats in the NAOS Project. Front. Mar. Sci. 2020, 7. [Google Scholar] [CrossRef]
  6. Canepa, E.; Pensieri, S.; Bozzano, R.; Faimali, M.; Traverso, P.; Cavaleri, L. The ODAS Italia 1 buoy: More than forty years of activity in the Ligurian Sea. Prog. Oceanogr. 2015, 135, 48–63. [Google Scholar] [CrossRef]
  7. Bozzano, R.; Pensieri, S.; Pensieri, L.; Cardin, V.; Brunetti, F.; Bensi, M.; Petihakis, G.; Margo Tsagaraki, T.; Ntoumas, M.; Podaras, D.; et al. The M3A Network of Open Ocean Observatories in the Mediterranean Sea. In Proceedings of the MTS/IEEE OCEANS 2013, Bergen, Norway, 10–13 June 2013; pp. 1–10. [Google Scholar] [CrossRef]
  8. Aquapro Comprehensive Manual, Nortek AS, 2013. Available online: https://www.nortekgroup.com/assets/documents/ComprehensiveManual_Oct2017_compressed.pdf (accessed on 26 April 2021).
  9. Raiteri, G.; Bordone, A.; Ciuffardi, T.; Schirone, A.; Pennecchi, F.; Bozzano, R.; Pensieri, S. A spar buoy-mounted ADCP measurement station in the Ligurian Sea: A metrological approach to correct current measures for bias effects and evaluate uncertainties. In Proceedings of the 2018 IEEE International Workshop on Metrology for the Sea: Learning to Measure Sea Health Parameters (MetroSea), Bari, Italy, 8–10 October 2018; pp. 82–87. [Google Scholar] [CrossRef]
  10. Movable Type Scripts. Calculate Distance, Bearing and More between Latitude/Longitude Points. Available online: http://www.movable-type.co.uk/scripts/latlong.html (accessed on 1 March 2021).
  11. Ciuffardi, T.; Napolitano, E.; Iacono, R.; Reseghetti, F.; Raiteri, G.; Bordone, A. Analysis of surface circulation structures along a frequently repeated XBT transect crossing the Ligurian and Tyrrhenian seas. Ocean Dynam. 2016, 66, 767–783. [Google Scholar] [CrossRef]
  12. Prieur, L.; D’ortenzio, F.; Taillandier, V.; Testor, P. Physical Oceanography of the Ligurian Sea. In The Proceedings of The Mediterranean Sea in the Era of Global Change 1: 30 Years of Multidisciplinary Study of the Ligurian Sea; Migon, C., Nival, P., Sciandra, A., Eds.; John Wiley & Sons, Inc.: Hoboken, NJ, USA, 2020; pp. 139–164. ISBN 1-78630-428-7. [Google Scholar] [CrossRef]
  13. Picco, P.; Cappelletti, A.; Sparnocchia, S.; Schiano, M.E.; Pensieri, S.; Bozzano, R. Upper layer current variability in the Central Ligurian Sea. Ocean Sci. 2010, 6, 825–836. [Google Scholar] [CrossRef]
  14. Vandenbulcke, L.; Beckers, J.M.; Barth, A. Correction of inertial oscillations by assimilation of HF radar data in a model of the Ligurian Sea. Ocean Dyn. 2017, 67, 117–135. [Google Scholar] [CrossRef]
  15. Bozzano, R.; Fanelli, E.; Pensieri, S.; Picco, P.; Schiano, M.E. Temporal variations in zooplankton biomass in the Ligurian Sea inferred from long time series of ADCP data. Ocean Sci. 2014, 10, 93–105. [Google Scholar] [CrossRef] [Green Version]
  16. Iacono, R.; Napolitano, E. Aspects of the summer circulation in the eastern Ligurian Sea. Deep Sea Res. Part I Oceanogr. Res. Pap. 2020, 166, 03407. [Google Scholar] [CrossRef]
  17. Clementi, E.; Pistoia, J.; Escudier, R.; Delrosso, D.; Drudi, M.; Grandi, A.; Lecci, R.; Cretí, S.; Ciliberti, S.; Coppini, G.; et al. Mediterranean Sea Analysis and Forecast (CMEMS MED-Currents, EAS5 System). Copernicus Marine Environment Monitoring Service (CMEMS) 2019. Available online: https://www.earth-prints.org/bitstream/2122/13246/1/CMEMS-MED-PUM-006-013.pdf (accessed on 26 April 2021).
  18. Escudier, R.; Clementi, E.; Omar, M.; Cipollone, A.; Pistoia, J.; Aydogdu, A.; Drudi, M.; Grandi, A.; Lyubartsev, V.; Lecci, R.; et al. Mediterranean Sea Physical Reanalysis (CMEMS MED-Currents) (Version 1). Copernicus Monitoring Environment Marine Service (CMEMS) 2020. Available online: https://resources.marine.copernicus.eu/?option=com_csw&view=details&product_id=MEDSEA_MULTIYEAR_PHY_006_004 (accessed on 26 April 2021). [CrossRef]
  19. JCGM 100:2008. Evaluation of Measurement Data—Guide to the Expression of Uncertainty in Measurement (GUM 1995 with Minor Corrections); BIPM: Sèvres, France, September 2008. [Google Scholar]
  20. JCGM 101:2008. Evaluation of Measurement Data—Supplement 1 to the “Guide to the Expression of Uncertainty in Measurement”—Propagation of Distributions using a Monte Carlo Method; BIPM: Sèvres, France, 2008. [Google Scholar]
  21. The Comprehensive R Archive Network. Available online: https://cran.r-project.org/ (accessed on 1 March 2021).
Figure 1. The spar buoy-mounted ADCP measurement station in the Ligurian Sea: (a) the W1-M3A observatory location, shown as a red four points star; (b) the buoy on land and in operation at sea.
Figure 1. The spar buoy-mounted ADCP measurement station in the Ligurian Sea: (a) the W1-M3A observatory location, shown as a red four points star; (b) the buoy on land and in operation at sea.
Jmse 09 00466 g001
Figure 2. ADCPs comparison test: (a) mechanical mounting configuration; (b) working area (red dot: deployment point); (c) deployment.
Figure 2. ADCPs comparison test: (a) mechanical mounting configuration; (b) working area (red dot: deployment point); (c) deployment.
Jmse 09 00466 g002
Figure 3. ADCP installation on the buoy. Assembly on the buoy damping disk: (a) overall view; (b) detail view; (c) re-installation of the ADCP on the mounting system following the battery replacement.
Figure 3. ADCP installation on the buoy. Assembly on the buoy damping disk: (a) overall view; (b) detail view; (c) re-installation of the ADCP on the mounting system following the battery replacement.
Jmse 09 00466 g003
Figure 4. ADCP working pressure and temperature during the whole deployment period.
Figure 4. ADCP working pressure and temperature during the whole deployment period.
Jmse 09 00466 g004
Figure 5. ADCP acoustic signals: (a) average counts of the three beams vs. cell depths; (b) SNR values.
Figure 5. ADCP acoustic signals: (a) average counts of the three beams vs. cell depths; (b) SNR values.
Jmse 09 00466 g005
Figure 6. Heading values detected by buoy and ADCP in the period from 30 March to 29 August 2017.
Figure 6. Heading values detected by buoy and ADCP in the period from 30 March to 29 August 2017.
Jmse 09 00466 g006
Figure 7. Positions of the buoy during the considered period, detected by GPS signal.
Figure 7. Positions of the buoy during the considered period, detected by GPS signal.
Jmse 09 00466 g007
Figure 8. Average heading values detected by buoy and ADCP: the theoretical average value of the ADCP heading is obtained by overturning the average value of the ODAS heading by 180°.
Figure 8. Average heading values detected by buoy and ADCP: the theoretical average value of the ADCP heading is obtained by overturning the average value of the ODAS heading by 180°.
Jmse 09 00466 g008
Figure 9. Example of application of the buoy heading correction: current values measured at cell no. 4 before (on the left, in red) and after the correction (on the right, in blue).
Figure 9. Example of application of the buoy heading correction: current values measured at cell no. 4 before (on the left, in red) and after the correction (on the right, in blue).
Jmse 09 00466 g009
Figure 10. Trend of the error committed by the ADCP compass as a function of the rotation of the buoy (measured by the ODAS compass), computed for each individual profile acquired.
Figure 10. Trend of the error committed by the ADCP compass as a function of the rotation of the buoy (measured by the ODAS compass), computed for each individual profile acquired.
Jmse 09 00466 g010
Figure 11. ADCP currents measured cell by cell, already corrected for buoy heading but still affected by the apparent current values induced by the buoy displacements over time (obtained from GPS).
Figure 11. ADCP currents measured cell by cell, already corrected for buoy heading but still affected by the apparent current values induced by the buoy displacements over time (obtained from GPS).
Jmse 09 00466 g011
Figure 12. Corrected modules and directions of the currents measured by the ADCP in two cells of the profile.
Figure 12. Corrected modules and directions of the currents measured by the ADCP in two cells of the profile.
Jmse 09 00466 g012
Figure 13. Progressive vectors at each cell (a) and a zoom (b) on a sub-sample dominated by inertial variability.
Figure 13. Progressive vectors at each cell (a) and a zoom (b) on a sub-sample dominated by inertial variability.
Jmse 09 00466 g013
Figure 14. Comparison of daily mean values of currents magnitude at the selected cell.
Figure 14. Comparison of daily mean values of currents magnitude at the selected cell.
Jmse 09 00466 g014
Figure 15. Comparison of rotary spectra at the selected cell (clockwise and anticlockwise).
Figure 15. Comparison of rotary spectra at the selected cell (clockwise and anticlockwise).
Jmse 09 00466 g015
Figure 16. East and north components of monthly average current profiles, associated respectively with the experimental measurements (in blue) and the data provided by the model (in red).
Figure 16. East and north components of monthly average current profiles, associated respectively with the experimental measurements (in blue) and the data provided by the model (in red).
Jmse 09 00466 g016
Figure 17. East and north current components trend over time, with respect to depth (top row: experimental data measured by the ADCP-bottom row: data obtained from the CMEMS forecast model).
Figure 17. East and north current components trend over time, with respect to depth (top row: experimental data measured by the ADCP-bottom row: data obtained from the CMEMS forecast model).
Jmse 09 00466 g017
Figure 18. Corrected velocity modules and their standard uncertainty vs time (first row) and corrected velocity directions and their standard uncertainty vs time (second row).
Figure 18. Corrected velocity modules and their standard uncertainty vs time (first row) and corrected velocity directions and their standard uncertainty vs time (second row).
Jmse 09 00466 g018
Table 1. Actual depths at which the centers of the ADCP measurement cells are located: depths for which the current measurements were not considered valid are shown in gray background. Negative depth indicates the cell whose center is in the air, above the surface of the sea.
Table 1. Actual depths at which the centers of the ADCP measurement cells are located: depths for which the current measurements were not considered valid are shown in gray background. Negative depth indicates the cell whose center is in the air, above the surface of the sea.
Cell No.Cell Center Distance
from Sensor Head
Cell Center
Depth
mm
1532.2
2928.2
31324.2
41720.2
52116.2
62512.2
7298.2
8334.2
9370.2
1041−3.8
Table 2. Period of measurements and settings of the ADCP measurements used for the comparison.
Table 2. Period of measurements and settings of the ADCP measurements used for the comparison.
Periods of MeasurementsADCP Depthdzdt
(a) 30/3–24/5 200462 m8m1 h
(b) 30/3–16/4 2005100 m8 m1 h
(c) 19/4–29/8 200580 m8 m30 min
Table 3. Temporal averaged magnitude with standard deviation and direction of currents at the cell including 22–24 m depth. Letter a stands for current measurements collected in 2004, while bc indicates the two time series merged for year 2005. O indicates the ADCP measurements from the spar buoy.
Table 3. Temporal averaged magnitude with standard deviation and direction of currents at the cell including 22–24 m depth. Letter a stands for current measurements collected in 2004, while bc indicates the two time series merged for year 2005. O indicates the ADCP measurements from the spar buoy.
22–24 mabcO
Magnitude cm s−113.96.817.7
Magnitude std cm s−1 7.04.76.4
Direction °N330280294
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Bordone, A.; Ciuffardi, T.; Raiteri, G.; Schirone, A.; Bozzano, R.; Pensieri, S.; Pennecchi, F.; Picco, P. Improved Current Estimates from Spar Buoy-Mounted ADCP Measurement Station: A Case Study in the Ligurian Sea. J. Mar. Sci. Eng. 2021, 9, 466. https://doi.org/10.3390/jmse9050466

AMA Style

Bordone A, Ciuffardi T, Raiteri G, Schirone A, Bozzano R, Pensieri S, Pennecchi F, Picco P. Improved Current Estimates from Spar Buoy-Mounted ADCP Measurement Station: A Case Study in the Ligurian Sea. Journal of Marine Science and Engineering. 2021; 9(5):466. https://doi.org/10.3390/jmse9050466

Chicago/Turabian Style

Bordone, Andrea, Tiziana Ciuffardi, Giancarlo Raiteri, Antonio Schirone, Roberto Bozzano, Sara Pensieri, Francesca Pennecchi, and Paola Picco. 2021. "Improved Current Estimates from Spar Buoy-Mounted ADCP Measurement Station: A Case Study in the Ligurian Sea" Journal of Marine Science and Engineering 9, no. 5: 466. https://doi.org/10.3390/jmse9050466

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop