European Radiometry Buoy and Infrastructure (EURYBIA): A Contribution to the Design of the European Copernicus Infrastructure for Ocean Colour System Vicarious Calibration

: In the context of the Copernicus Program, EUMETSAT prioritizes the creation of an ocean color infrastructure for system vicarious calibration (OC-SVC). This work aims to reply to this need by proposing the European Radiometry Buoy and Infrastructure (EURYBIA). EURYBIA is designed as an autonomous European infrastructure operating within the Marine Optical Network (MarONet) established by University of Miami (Miami, FL, USA) based on the Marine Optical Buoy (MOBY) experience and NASA support. MarONet addresses SVC requirements in di ﬀ erent sites, consistently and in a traceable way. The selected EURYBIA installation is close to the Lampedusa Island in the central Mediterranean Sea. This area is widely studied and hosts an Atmospheric and Oceanographic Observatory for long-term climate monitoring. The EURYBIA ﬁeld segment comprises o ﬀ -shore and on-shore infrastructures to manage the observation system and perform routine sensors calibrations. The ground segment includes the telemetry center for data communication and the processing center to compute data products and uncertainty budgets. The study shows that the overall uncertainty of EURYBIA SVC gains computed for the Sentinel-3 OLCI mission under EUMETSAT protocols is of about 0.05% in the blue-green wavelengths after a decade of measurements, similar to that of the reference site in Hawaii and in compliance with requirements for climate studies.


Introduction
The European Copernicus Program for Earth monitoring (http://www.copernicus.eu/) aims to deliver remote sensing capabilities, field measurements, and data processing services that provide users with reliable and up-to-date information related to environmental, climate, and security issues. As part of the Copernicus space component, the European Commission is developing, operating, and planning a series of Sentinel missions for the next 15 years, together with ESA and EUMETSAT. Among these missions, the Sentinel-3 Ocean and Land Colour Instrument (S3 OLCI) is a key component for monitoring the biogeochemical status of the ocean. OLCI is a primary source of ocean color (OC) data for the Copernicus Marine Environment Monitoring Service (CMEMS). For instance, OLCI data are operationally used to advance our understanding of living marine environments and assimilated in the CMEMS models to forecast the evolution of marine ecosystems, together with the other OC third party missions. The Sentinel-2 mission, although mainly developed for land applications, already demonstrates a great potential for coastal applications like sediment and chlorophyll retrieval (e.g., reference [1]), and it is expected to become an operational complement of the CMEMS next phases. OC is also an important dataset for Copernicus Climate Change Service (C3S) to monitor the impact of climate change on marine ecosystems and the ocean carbon cycle [2].
OC science is making an increasing contribution to Ocean Observing systems thanks to continuous global monitoring of biogeochemical variables, to the point that OC missions have become fundamental to the success of the Copernicus Marine and Climate services. This has also raised expectations on the quality of data products. The current threshold requirement for calibration uncertainty of Copernicus satellite OC sensors is 0.5% in the blue-green spectral region, with a further goal of 0.3% [3][4][5]. OC System Vicarious Calibration (SVC) [4,6] complements pre-launch and on-board calibrations through highly accurate in-situ measurements of water-leaving radiances. These measurements are indeed the principal source for the vicarious calibration of space-born radiometric data. By reducing residual biases in water-leaving radiances, SVC is currently the only way to attain the target OC product uncertainties [7], satisfy OC mission requirements [8], and enable marine and climate data services.
Presently, OC-SVC relies on two sites-the Marine Optical BuoY (MOBY; https://www.mlml. calstate.edu/moby/) [9] and the BOUée pour l'acquiSition d'une Série Optique à Long termE (BOUSSOLE; http://www.obs-vlfr.fr/Boussole/) [10] -established in the last two decades with different instrumental designs, respectively, in the Pacific Ocean and in the Mediterranean Sea. None of these systems is, however, fully compliant with the operational requirements of the Copernicus Programme. The creation of an OC-SVC infrastructure dedicated to the Copernicus Programme is a key point to improve the performance of OC missions that meet the long-term and high-quality standards underpinning marine bio-geochemistry data products delivered by Climate and Marine Services [5]. EUMETSAT has identified a series of steps to achieve this goal. The first step accomplished in 2017 defined the requirements for the Copernicus OC-SVC Infrastructure [3]. These requirements represent the baseline for the present second step: preliminary design of the Copernicus OC-SVC Infrastructure. The overall objective of this paper is to contribute to this European effort by designing the OC-SVC infrastructure able to provide highly accurate, state-of-the-art Fiducial Reference Measurements (FRM) in the very tight development schedule (36 months from the start development) requested by Copernicus Programme. Our design is to commit to the highest radiometry standards and fast operational readiness. This work will then contribute to the third step: technical definition, specifications, and engineering design that are expected to start in 2020.
The main objective of this paper is to describe the proposed European Radiometry Buoy and Infrastructure (EURYBIA), named for the goddess of power over and mastery of the sea (Eὐρυβία: wide-force; who presided over external forces that influence the main such as the rise of the constellations, seasonal weather, and the power of the winds; Theoi Greek Mythology, www.theoi.com.) for SVC of current and future Copernicus missions (Sentinel-3 A, B, C, D, Sentinel-2 C, D), as well as of third party missions (e.g., NOAA/VIIRS, NASA/PACE). The paper will demonstrate that the proposed OC-SVC Infrastructure ( Figure 1) meets these requirements in terms of core optical system and field Remote Sens. 2020, 12,1178 3 of 41 deployment structure hosting it, data processing, and delivery. The deployment site selection will be justified on the basis of environmental conditions, logistic factors, and a complete estimate of the uncertainty budget of the products for the proposed OC-SVC infrastructure.
Remote Sens. 2020, 12, x FOR PEER REVIEW 3 of 41 deployment site selection will be justified on the basis of environmental conditions, logistic factors, and a complete estimate of the uncertainty budget of the products for the proposed OC-SVC infrastructure. Based on OC-SVC requirements, this paper is structured as follows. Section 2 describes the core EURYBIA optical system, the characterization and calibration procedures, as well as acquisition protocol. The optical buoy and the mooring buoy for the deployment of the optical system and ancillary equipment are presented in Section 3. Solutions featured by the Ground Segment for data telemetry, processing, quality control, dissemination, and the scheduling of maintenance operations are discussed in Section 4. Regarding the location of the site, previous studies [4,[11][12][13] investigated the suitability of a set of candidate regions using, as reference metrics, the number of valid matchups resulting from a set of single threshold criteria applied to a restricted set of variables. With respect to the previous studies, the OC-SCV requirements [3] encompass a larger number of variables and corresponding constraints taken into account by this paper. This environment assessment includes oceanographic and atmospheric properties and adjacency effect of the proposed selected sites and is summarized in Section 5. Section 6 summarizes the suitability of the site by developing a model to estimate the uncertainty on the vicarious calibration gains and by applying it to OLCI observations [14] and to the results of the characterization study. Final conclusions are addressed in Section 7. Based on OC-SVC requirements, this paper is structured as follows. Section 2 describes the core EURYBIA optical system, the characterization and calibration procedures, as well as acquisition protocol. The optical buoy and the mooring buoy for the deployment of the optical system and ancillary equipment are presented in Section 3. Solutions featured by the Ground Segment for data telemetry, processing, quality control, dissemination, and the scheduling of maintenance operations are discussed in Section 4. Regarding the location of the site, previous studies [4,[11][12][13] investigated the suitability of a set of candidate regions using, as reference metrics, the number of valid match-ups resulting from a set of single threshold criteria applied to a restricted set of variables. With respect to the previous studies, the OC-SCV requirements [3] encompass a larger number of variables and corresponding constraints taken into account by this paper. This environment assessment includes oceanographic and atmospheric properties and adjacency effect of the proposed selected sites and is summarized in Section 5. Section 6 summarizes the suitability of the site by developing a model to estimate the uncertainty on the vicarious calibration gains and by applying it to OLCI observations [14] and to the results of the characterization study. Final conclusions are addressed in Section 7.

The Optical System
This section provides a summary of the proposed optical system (2.1) and its characterization and calibration procedures (2.2). Additional information can be found in the existing literature and technical documents [15][16][17][18].

Optical System Description
The core instrumentation of EURYBIA is the optical system of the Marine Optical Network (MarONet), developed by the MOBY team on the basis of the long-term experience acquired in the last three decades. In order to fulfil spectral range and resolution requirements, the system actually consists of a common optical path from the radiation collectors to a device that directs the measured signals into two similar single spectrographs covering different spectral regions and named, on the basis of the range the "Blue" (BSG: 350-690 nm) and the "Red" (RSG: 550-900 nm). Figure 2 shows a scheme of the optical system. The optical signal is collected by the optical heads deployed into the water at the end of each of the three arms (upwelling radiance and downward irradiance measurements, Lu and Ed respectively, ) and on the top of the mast of the Optical Buoy (above-water irradiance measurement, Es) (see Section 3). The structure of the optical head depends on the radiometric variable (radiance vs. irradiance) to be acquired, as detailed by Zibordi et al. [12]. To mitigate bio-fouling effects on the optical windows, they are surrounded by a copper bezel and equipped with a UV anti-biofouling module consisting of a fiber optic coupled source of sterilizing UV light (285 nm). The exposure scheme will be programmable to be optimized for the actual site characteristics.
Optical fibers transfer the radiation from each collector to the bottom of the buoy where the two spectrometers are mounted. The signal is split with a polka-dot mirror and transmitted to the two spectrometers. In the following, we describe a single spectrometer. The optical fibers from each collector are connected to a shutter block consisting of 14 single shutters that can also host neutral density filters to adjust the signal level. These devices are also connected, by means of optical fibers, the reference lamps for system monitoring during deployment. Each shutter can be activated independently, allowing the sensor to be illuminated with all possible combinations of inputs.
Optical fibers are aligned on the entrance slit of the spectrometer to allow for the projection on the spectral dispersion system. The alignment order is determined to minimize stray light between different inputs. Band pass filters and short and long wave cutoff filters are mounted at the entrance The optical signal is collected by the optical heads deployed into the water at the end of each of the three arms (upwelling radiance and downward irradiance measurements, Lu and Ed, respectively) and on the top of the mast of the Optical Buoy (above-water irradiance measurement, Es) (see Section 3). The structure of the optical head depends on the radiometric variable (radiance vs. irradiance) to be acquired, as detailed by Zibordi et al. [12]. To mitigate bio-fouling effects on the optical windows, they are surrounded by a copper bezel and equipped with a UV anti-biofouling module consisting of a fiber optic coupled source of sterilizing UV light (285 nm). The exposure scheme will be programmable to be optimized for the actual site characteristics.
Optical fibers transfer the radiation from each collector to the bottom of the buoy where the two spectrometers are mounted. The signal is split with a polka-dot mirror and transmitted to the two spectrometers. In the following, we describe a single spectrometer. The optical fibers from each collector are connected to a shutter block consisting of 14 single shutters that can also host neutral density filters to adjust the signal level. These devices are also connected, by means of optical fibers, the reference lamps for system monitoring during deployment. Each shutter can be activated independently, allowing the sensor to be illuminated with all possible combinations of inputs. Optical fibers are aligned on the entrance slit of the spectrometer to allow for the projection on the spectral dispersion system. The alignment order is determined to minimize stray light between different inputs. Band pass filters and short and long wave cutoff filters are mounted at the entrance of the spectral dispersion system to reduce the effects of radiation outside the respective spectral range of the BSG or SRG spectrometer. The spectral dispersion system is a transmission Volume Phase Holographic spectrograph (VPH) distributed by Resonon Inc., Bozeman, MT 59715, USA. The VPH reduces the stray light by more than one order of magnitude with respect to the concave reflective gratings adopted in MOBY Heritage [19]. The VPH allows imaging all the environmental optical inputs at once: downwelling sky irradiance (Es), downwelling irradiance measurements from the three arms (Ed Top , Ed Mid , and Ed Bot ), and upwelling irradiance from the three arms (Lu Top , Lu Mid , and Lu Bot ).
The radiation sensor consists of a scientific Andor Ikon M934 camera, Oxford, UK. The camera is a back illuminated Charged Coupled Device (CCD) with a thicker active photosensitive area (i.e., deep depletion CCD) and fringe suppression to prevent etalon effects. The camera is retrofitted with wedged windows to minimize straylight effect (UV ghosting). Different base settings allow the definition of data acquisition characteristics and the sensor working temperature to optimize the camera performances. System synchronization capabilities are used to control the integrated C-mount shutter that is closed during readout to avoid vertical smear. The system is thermo-electrically cooled to reduce the dark current. For MarONet, the cooling temperature has been set to −60 • C as trade-off between noise reduction, power consumption, long term stability, and lowest drift. Cooling the CCD array provides stability to the camera system eliminating bright and dark pixels while also reducing thermal noise and dark current. An internal reference system is incorporated into the optical system to monitor the stability of the radiometric detectors, the electronics, and internal optics. These measurements are critical for establishing confidence in the observations acquired during a deployment cycle. The reference system consists of one optical feedback-regulated incandescent lamp and two current regulated, temperature stabilized LEDs. The diodes emit at a wavelength centered at 465 nm (Blue LED) and 705 nm (Red LED), with a bandwidth of approximately 100 nm. The lamps are run with current controlled circuitry and the temperature of the lamp holder block is monitored. Each source output is channeled through a one-inch Spectralon sphere to split between the RSG and BSG. The diffuser is illuminated in sequence by the incandescent lamp, the blue LED, and the red LED. These lamps are observed at the end of each data acquisition set. The internal reference sources give no information on changes of optical properties of the single collectors (for example due to bio-fouling) since these components are not included in the optical path of the reference lamps signals.
The system is designed to allow for fully (also remotely) programmable measurement protocols. This functionality, which can be done remotely, will be exploited mostly in the implementation phase to optimize the characteristics of the measurement protocol scenarios needed for SVC. Table 1 summarizes the characteristics of the proposed optical system and compares them against the corresponding requirements.
A set of additional sensors are associated with the optical system to: • monitor the status of the system: e.g., temperature of different sub-components, power supply, coolant flow, etc.; • provide ancillary/auxiliary variables for data processing and/or data quality assessment: e.g., sea water temperature and salinity, top arm depth sensor, tilt along two horizontal axes, system heading.

System Characterization and Calibration
The main objectives of the system characterization are to Remote Sens. 2020, 12, 1178 6 of 41 • allow the identification of interactions and dependencies between optical and electronic components; • identify, evaluate, and correct any systematic disturbances; • allow the design of the Level-0 to Level-1 data processing, estimate input variables, and verify the radiometric model (see Section 4); • allow the estimate of the contribution of each component to the uncertainty budget associated with the radiometric measurement (see Section 6). A detailed description of the characterization activities carried out in the period 2014-2018 can be found in Voss et al. (2018) [15]. The complete system characterization should be done for each new unit when the optical system is built and fully integrated. Leaving the details of the characterization to dedicated existing or future literature, we list the following relevant issues that are planned for the system characterization: • CCD tracks binning factor; In addition to the complete system-level characterization when fully integrated, extensive partial characterizations are planned to be performed pre-and post-deployment (see Section 3). The aspects that need to be periodically characterized will depend on the results of the pre-operational phase to be carried out after the infrastructure deployment. For this purpose, a set of tests that can monitor the status of the system can be introduced in the measurement protocol and/or data processing.
The objective of the radiometric calibration is to obtain absolute values of radiometric measurements compliant with the International System of Units (SI) and directly traceable to National Metrology Institute (NMI) primary standards. This is obtained by absolutely calibrating each system before and after optical buoy deployment and by monitoring the radiometric calibration during each deployment [18].

EURYBIA Field Segment
Based on the MarONet system design [5], the EURYBIA Infrastructure consists of a twin buoy system: an optical buoy tethered to a mooring buoy anchored to the sea bottom ( Figure 3). The first buoy carries the optical system to perform the radiometric measurements, whereas the second buoy hosts additional scientific meteorological and oceanographic instrumentation and collects ancillary data. In addition to the complete system-level characterization when fully integrated, extensive partial characterizations are planned to be performed pre-and post-deployment (see Section 3). The aspects that need to be periodically characterized will depend on the results of the pre-operational phase to be carried out after the infrastructure deployment. For this purpose, a set of tests that can monitor the status of the system can be introduced in the measurement protocol and/or data processing.
The objective of the radiometric calibration is to obtain absolute values of radiometric measurements compliant with the International System of Units (SI) and directly traceable to National Metrology Institute (NMI) primary standards. This is obtained by absolutely calibrating each system before and after optical buoy deployment and by monitoring the radiometric calibration during each deployment [18].

EURYBIA Field Segment
Based on the MarONet system design [5], the EURYBIA Infrastructure consists of a twin buoy system: an optical buoy tethered to a mooring buoy anchored to the sea bottom ( Figure 3). The first buoy carries the optical system to perform the radiometric measurements, whereas the second buoy hosts additional scientific meteorological and oceanographic instrumentation and collects ancillary data.  Indeed, collecting measurements useful for OC-SVC implies strict constraints on the motion of the buoy that, on one hand, must follow the wavy surface to keep each optical collector at a constant depth and, on the other, must maintain tilting oscillations to within a few degrees in order not to bias the optical measurements.
In order to guarantee the operational measurements, we propose to design and build two identical infrastructures in order to enable rotation in the field, always having one system operational at sea and the other ashore for maintenance operations to make it ready for the next deployment. This design allows the system to be compliant with the critical requirements of the measurement traceability achievable only with frequent periodic on-site and on-land maintenance and calibration operations while ensuring the continuity of the long-term observations.

Optical Buoy
The optical buoy, hosting the radiometric system, consists of a medium-size spar buoy [9] whose actual dimensions, distribution of masses along the structure, and tether connection to the mooring buoy are properly designed to optimize the buoy stability due to local environmental conditions (i.e., reduced tilting with respect to the sea state conditions) and self-shading effects. Figure 3 shows a preliminary scheme of the optical buoy; the reported dimensions and depths of the three arms hosting the in-water optical collectors are based on MarONet design. The final design, optimized to the environmental conditions of the selected site, might require minor adjustments. Indeed, a spar buoy is characterized by a reduced section at the water surface, a deep draft, a long vertical hull, a disk on the deepest level, and a ballast. These characteristics, relevant to reduce heave (i.e., the buoy vertical motion) and tilting also allow the deployment of instruments along the long vertical hull. For the EURYBIA application, the buoy is also designed to follow the surface in order to keep the optical instruments at constant depths with respect to the sea surface, which is achieved by ensuring a substantial capacity of floating reaction.
The structure of the optical buoy can be divided into the following main parts, each one contributing to the overall structure buoyancy while housing instruments and electronic components essential for the system operational activity: a top floatation element, a hull, and a lower bay.
The top floatation element hosts a payload consisting of an Es irradiance cosine collector, antennas for the positioning and communication system, a radar reflector, a lantern, photovoltaic panels and batteries to power up all electronic systems, an industrial controller collecting data and operating the whole system through A/D converters and/or serial channels, a modem for the communication system, and controllers. Along the hull, at different fixed depths, three arms support the in-water optical collectors-one for downwelling irradiance and one for upwelling radiance. The arms have different lengths, the longer being closest to the surface and the shorter at the deepest position in order to avoid any interference on upwelling radiance measurements. On the shallower arm, a pressure sensor is also installed in order to monitor the depth of the upper arm measurement.
Two industrial data loggers, connected through an Ethernet connection, are deployed on the optical spar buoy: one data logger collects the optical measurements, located underwater in the watertight lower bay, and a second data logger, which is installed in the watertight compartment on the main floatation element, collects all other data and monitoring parameters (i.e., orientation of the arms, temperatures inside the watertight compartment, depth of the shallower arm, trim of the buoy, current produced by the photovoltaic panels, voltage of the batteries, etc.). At the bottom of the buoy, an orientation module consisting of a pressure sensor and an attitude and heading reference system are installed in order to monitor the trim of the collector and to provide information for any compensation of the acquired irradiance. All data acquired on board are packed into a data cluster, timestamped, and geographically annotated by the upper controller, locally stored on the hard disk while constituting the stream of bytes to be transferred ashore with a communication (telemetry) system.
All the above-mentioned electronic systems and sensors are powered by the photovoltaic panels and batteries installed on the infrastructure itself, which are self-sustained for the entire period of activity, considering up to four acquisition cycles per day, also including the payload required by the necessary alert and monitoring services.

Mooring Buoy
A mooring buoy consisting in a standard disc-shaped surface buoy is anchored to the sea bottom and is tethered to the spar buoy to keep the optical buoy in place. It also hosts additional equipment to support the radiometric data collection, processing, quality assessment, and results interpretation.
The following parts of the payload of the mooring buoy are similar to the those installed in the optical buoy: antennas for the positioning and communication system, a radar reflector, a lantern, photovoltaic panels and batteries to power all electrical systems, an industrial controller collecting data and operating the whole system through A/D converters and/or serial channels, a modem for the communication, and controllers for the photovoltaic panels.
Additionally, the mooring buoy hosts an attitude and heading reference system to measure the motion of the buoy and instruments for collecting ancillary measurements above the sea surface and underwater. A meteorological station is required to process and interpret the radiometric observations. Optionally, above-water instrumentation may include radiometric measurements with the main objective of estimating the presence of clouds, and in-water standard bio-optical instrumentation could be installed to support main radiometric data processing and quality assessment.
As for the optical buoy, all above-mentioned electronic systems and sensors get their electrical power supply from a system of photovoltaic panels, charge regulators, and batteries.

Ground Segment
Ground segment (GS) tasks include telemetry operations, processing of the in-situ measurements to derive radiometric products and related uncertainties for the computation of SVC gains, dissemination of results, documentation of processing code, and scheduling of maintenance operations. The telemetry center and the processing center described in the next sub-sections are the two main infrastructural components that allow for GS activities, as detailed hereafter.

The Telemetry Center
The telemetry center executes direct and independent data transfers with the mooring and optical buoy by means of GSM communication. Its specific tasks are (1) uploading configuration files with instrumentation settings and acquisition protocols, (2) uploading firmware programs for the EURYBIA equipment, and software modules for the onboard controllers, (3) executing preliminary data processing and checks, (4) downloading in-situ measurements, and (5) verifying the status of buoys at sea and alerting in the case of out-of-order conditions. This latter assignment relies on a continuous (24/7) acquisition of monitoring parameters, even when radiometric and auxiliary measurements are not collected. Integrity checks are performed to detect transfer errors (checksum), missing values in data records, or out-of-range conditions (e.g., temperature, position and power supply). Redundant communication and storage resources, as well as functionalities to remotely supervise and restart telemetry operations, ensure continuous control of the buoys system at sea. In case of anomalies, including missing communication signal, alerts are immediately released to start recovery actions.

The Processing Centre
In-situ measurements are transferred from the telemetry center to the processing center to compute data products. The data processing is structured in levels to control the complexity of code, simplify its maintenance and evolution, facilitate the assessment of intermediate results, and allow for starting data reprocessing on demand from the specified processing level ( Figure 4). The Level-0 processing accounts for the raw data conversion into a standard format. The Level-1 processing applies corrections for the non-ideal instrument response and computes calibrated values in physical units. The Level-2A proceeds to the minimization of environmental perturbations, performs depth propagation of in-water radiometric measurements and computes full-resolution data products such as the water-leaving radiance. The Level-2B finally derives data products at center wavelengths of specific space-borne sensors (spectral integration). Intermediate results and final products are archived as independent output files. Processing settings and metadata for the identification of input measurements and the version of the processing model are also recorded.

Level-0
The format of raw data produced by the radiometric system might change through time due to instrument updates or revisions of the EURYBIA data acquisition program. The Level-0 processing is then scoped to convert raw data into a standard format.

Level-1
Data recorded at the CCD of the Resonon spectrogram are expressed in Analog Digital Units (ADU). The Level-1 processing corrects these input data accounting for the non-ideal radiometer response and generates results in physical units. An additional task is a spectral calibration to relate CCD column-wise pixels to wavelengths.
The formalism adopted in this manuscript is derived from the "Protocols for Satellite Ocean Colour Data Validation: In-situ Optical Radiometry" report of the International Ocean-Colour Coordinating Group (IOCCG) [7]. Denoting with the symbol ℑ(λ) either (λ) or (λ) in physical units (W • cm • sr • nm and W • cm • nm , respectively), the Level-1 conversion equation is The Level-0 processing accounts for the raw data conversion into a standard format. The Level-1 processing applies corrections for the non-ideal instrument response and computes calibrated values in physical units. The Level-2A proceeds to the minimization of environmental perturbations, performs depth propagation of in-water radiometric measurements and computes full-resolution data products such as the water-leaving radiance. The Level-2B finally derives data products at center wavelengths of specific space-borne sensors (spectral integration). Intermediate results and final products are archived as independent output files. Processing settings and metadata for the identification of input measurements and the version of the processing model are also recorded.

Level-0
The format of raw data produced by the radiometric system might change through time due to instrument updates or revisions of the EURYBIA data acquisition program. The Level-0 processing is then scoped to convert raw data into a standard format.

Level-1
Data recorded at the CCD of the Resonon spectrogram are expressed in Analog Digital Units (ADU). The Level-1 processing corrects these input data accounting for the non-ideal radiometer response and generates results in physical units. An additional task is a spectral calibration to relate CCD column-wise pixels to wavelengths.
The formalism adopted in this manuscript is derived from the "Protocols for Satellite Ocean Colour Data Validation: In-situ Optical Radiometry" report of the International Ocean-Colour Coordinating Group (IOCCG) [7]. Denoting with the symbol (λ) either L(λ) or E(λ) in physical units (W·cm −2 ·sr −1 ·nm −1 and W·cm −2 ·nm −1 , respectively), the Level-1 conversion equation is where: 1.
(λ) represents net radiometric data obtained by subtracting the dark contribution from signal data after scaling the raw CCD readouts by the measurement integration time (sec) and the binning factor (unit of a pixel, pix)-(λ) is thus in units of ADU·sec −1 ·pix −1 , see for instance [20]; 2.
F I (λ) is the immersion factor accounting for the change in responsivity of the sensor when immersed in water with respect to air (dimensionless).
Results obtained from measurements performed by the blue and red spectrographs are finally merged into a unique spectrum (for each sampling acquisition and radiometric sensor). This processing step requires as input the predefined cut-off wavelength at which the spectra are joined.

Level-2A
The first processing step is the minimization of shading effects due to the deployment structure and the radiometers housing by means of look-up tables (LUTs) produced with Monte Carlo (MC) simulations of the radiative transfer process for given boundary conditions [29][30][31]. The scheme to determine the shading correction factors executes, for the same environmental case, MC ray tracing with and without accounting for the presence of the optical buoy. The difference between the two results defines the shading correction coefficients [32].
The biofouling effect is evaluated by means of specific quality checks performed during the periodical cleaning of the radiometers. These monitoring results are used to evaluate changes in the instrument performance and compile information to assess the data products quality.
The sub-surface value of the upwelling radiance is derived from L u (z i , λ) measurements performed at three fixed depths z i (i = 1, 2, 3), with z 1 > z 2 > z 3 and assuming positive depth increasing towards the bottom [28,33]. Precisely, the value of the upwelling radiance just below the sea-surface L u (0 − , λ) is derived from the L u (z 1 , λ) measurement at the shallowest depth z 1 as where K L (0, 1, λ) is the attenuation coefficient of the diffuse upwelling radiance for the depths from 0 to z 1 -i.e., K 01 L [33]. The diffuse attenuation coefficient The value of K L (0, 1, λ) cannot be directly measured and is determined from L u (z 1 , λ), L u (z 2 , λ) and L u (z 3 , λ) values [33]. To this end, the data processing accounts for the increasing fraction of light generated by inelastic radiative processes (Raman effect and chlorophyll fluorescence) from the blue region of the spectrum, where energy is abundant, to the red region, where the incoming light is rapidly attenuated (this implies a K L tendency to decrease with depth for λ > 575 nm). MC simulations are employed to compute K L lookup tables accounting for Raman scattering [33]. The effect of chlorophyll fluorescence is instead estimated from in-situ measurements in the [0, z 1 ] interval performed with ancillary instrumentation of the moored buoy. The K 01 L estimate is finally modelled as a function of the K 12 L and K 13 L [33]. The water-leaving radiance [6,34,35], L w θ, φ, λ j , is computed as where φ is the azimuth angle, θ is the zenith angle in air, θ is the corresponding refracted nadir angle in water, ρ(θ, θ ) is the internal Fresnel reflectance, and n λ j , T, S is the refractive index of seawater, with T and S indicating the temperature and salinity, respectively [36]. The fully normalized water-leaving radiance, L wN (λ i ), is the radiance leaving the sea surface from the water body (i.e., excluding the contribution of incident light reflected above the sea surface) as if the sun was at the zenith, the Earth was at its mean distance from the sun and in the absence of any atmospheric loss [37]. L wN (λ i ) is determined as where (1) E s (λ i ) is the incident solar irradiance; (2) C f /Q θ 0 , λ j , Chla is the BRDF correction term obtained from look-up-tables as a function of the solar zenith angle θ 0 , the wavelength λ j , and the Chlorophyll-a concentration Chla [38]; and 3)

Level-2B
Full resolution data are determined in the 350-900 nm interval. Spectral integration is addressed to compute radiometric values for the channel identified by the central wavelength λ i of selected space-borne sensors as indicating with SRF i the spectral response function of the i − th band of the space-borne sensor. Note that the numerical integral corresponding to Equation (7) is performed after interpolation of SRF discrete values to the central wavelengths of the in-situ measurement. The numerical integration scheme is then designed to minimize the uncertainty budget. Selected ocean color space sensors of interest include

Quality Assurance and Quality Control
The determination of data products is complemented by the quality assurance and control (QA/QC) of results. This task is jointly based on objective parameters (e.g., coefficient of variation of CCD data at Level-1, stability of pre-and post-deployment calibration coefficients, presence of optical stratifications limiting the capability to perform depth propagation of the upwelling radiance at Level-2), as well as the visual inspection of results (e.g., abnormal presence of spikes in Level-1 data) and complementary information (e.g., cloud-cover conditions). Additional QA/QC is performed checking the consistency of individual records with historical measurements once the time series permits (in analogy to the procedure adopted by the MOBY Team). For a final assessment, the stability of calibration lamps within their life cycle (i.e., determined based on traceability with respect to standards performed by the National Metrology Institute) is also verified.
QA/QC results are finally summarized through the following labels: • GOOD indicates the highest quality data for SVC applications (reprocessed data only) or for validation and monitoring (these include near-real-time products release); • QUESTIONABLE refers to data products derived upon applying significant corrections (e.g., excessive sensor tilting or instability of the instrument response); and • BAD is applied to data that are unreliable and cannot be used even if the current set of correction schemes are applied (usually caused by unanticipated problems during data acquisition or clouds).

Data Products Dissemination
The data processing cycle accounts for the following three main release versions: 1.
Near Real-Time (NRT) results are delivered once the following operations have been executed: visual data inspection, data correction and calibration, computation of data products such as the normalized water-leaving radiance, and preliminary data quality assurance and control (QA/QC. The timely dissemination of high-quality daily data products derived from in-situ measurements is scoped to release measurements results for the routine assessment and monitoring of orbiting OC remote sensing sensors.

2.
Post-Deployment Calibration (PDC) results are computed by performing a data re-processing when the optical system has been recovered to shore and re-calibrated. Namely, pre-and post-deployment calibration coefficients are used to account for instruments responsivity change.
Any instrument re-characterization for the non-ideal instrument response is also to be taken into consideration. 3.
Finally, Traceable to Calibration Standards (TCS) of NMI supporting OC-SVC activities are determined by means of an additional data re-processing at the end of the operational cycle of calibration lamps (which is in the order of 50 working hours). At this stage, in fact, the lamps are sent to NMI for traceability. Data products are then recomputed when reference calibration sources are replaced (note that the source stability is continuously monitored at the calibration facilities to allow for prompt actions in the case of out-of-order variations).
All GS dissemination services are performed within the Copernicus infrastructure and in compliance with the Copernicus principles of offering free, full, and open access to data and processing models. An interactive data browser will be designed to allow end-users to look up field measurements and derived products for specified search conditions and generate a quick statistical summary of the selected data. Deliverables include details on the pre-and post-deployment calibration coefficients, as well as the stability of calibration sources. Additional services are (1) a web service for the documentation of processing modules and (2) a forum for registered users to provide general information, answer questions, and receive comments and suggestions.

Maintenance Operations
A web interface is finally dedicated to facilitate the maintenance of the ground segment (e.g., software revisions, hardware replacement) and the field segment (e.g., infrastructure turn-over, instrumentation calibration, cleaning, and replacement). Provided on-line services include the scheduling of maintenance operations and the posting of services logs.

Geography and Illumination
As discussed in Section 1, the area around the island of Lampedusa, south of the Strait of Sicily, is proposed as the EURYBIA candidate site for the OC-SVC infrastructure. Lampedusa, part of the small Pelagian archipelago, is a small island (about 20 km 2 surface area) in the open Mediterranean Sea, at 35.5 • latitude North. About 6000 inhabitants live in Lampedusa and Linosa, the two major islands of the archipelago. The Lampedusa area has been extensively studied with respect to geomorphology, marine and terrestrial ecology, and many environmental parameters. The archipelago is located on the African tectonic plate, at the southern shoulder of a Plio-Quaternary rift zone that shows deep fault-controlled structural depressions oriented along the North-West to South East direction (e.g., the grabens of Pantelleria, Linosa, and Malta). The associated deep depressions, just North of Lampedusa, are more than 1500 m deep and separate the African continental shelf from that of Sicily. The rift zone is characterized by a moderate seismicity, mostly located in the Linosa graben.
Geologically, the island of Lampedusa is formed by Upper Miocene neritic limestones [40]. The island is about 10 km wide, with a tabular morphology and a surface sloping toward the south and south-east. The maximum altitude is about 130 m.
Calculations of the solar illumination and occurrence of sun glint conditions have been carried out. The illumination conditions at Lampedusa and MOBY are similar in summer, while in winter, the available solar energy at Lampedusa is a minimum of 70% of that available at MOBY. Higher solar zenith angles associated with higher latitude at Lampedusa than at MOBY are favorable for minimizing the influence of sun glint conditions.

Existing Infrastructures
Relevant logistical and scientific infrastructures, and several natural protection areas are present in and around Lampedusa (http://www.ampisolepelagie.it/). An airport and a primary and a secondary port exist on the island. Regular daily flights and ferry boats connect Lampedusa and Sicily throughout the year.
With respect to scientific infrastructures, Lampedusa hosts a climate observatory (http://www. lampedusa.enea.it), which has been operational since 1997 and is composed of an atmospheric observatory (35.52 • N, 12.63 • E) [41] and an oceanographic observatory (35.49 • N, 12.47 • E) [42]. A large number of atmospheric and oceanographic parameters are routinely measured at Lampedusa; these data, together with satellite observations and model results, have been used for the characterization of the Lampedusa area for OC-SVC. In addition, they offer relevant and reliable information for the definition of the design requirements of the buoys structure in order to meet the necessary stability for a higher number of quality measurements throughout the year. The Lampedusa Climate Observatory contributes to global measurement networks and to European environmental research infrastructures, such as the Integrated Carbon Observation System (ICOS) and the Aerosol, Clouds, and Trace Gases (ACTRIS). Lampedusa also contributes to the national network of the European Multidisciplinary Seafloor and water-column Observatory (EMSO), research infrastructure, with a number of oceanographic measurements. The operational instruments include sun photometers, also part of the Aerosol Robotic Network (AERONET), a ceilometer, an aerosol and water vapor lidar, radiometers and spectrometers for downwelling and upwelling radiation and atmospheric remote sensing, oceanographic sensors for temperature, salinity, chlorophyll, CDOM, backscattering, and underwater radiation.
As part of the Lampedusa Climate Observatory, instruments and capabilities in radiometers and spectrometers characterization and calibration (see e.g., references [43][44][45]) are available on site.
Scientific data acquired at the Lampedusa Climate Observatory have been extensively used, and more than 200 international scientific papers referring to these data have been published so far.
Many studies make advantage of these data for inter-comparison, validation, and combined analyses of satellite observations.

Bathymetry
Two regions with bottom depth larger than 300 m at 35.5248 • N-12.7667 • E and 35.7430 • N-12.3579 • E, LMP1 and LMP2 respectively, have been identified as possible sites for OC-SVC Infrastructure ( Figure 5). These two sites are characterized by a negligible impact from fishing activities and ship traffic and are located at a reasonable distance from Lampedusa (less than 20 km), thus enabling technical operations and allowing for timely interventions, while avoiding any possible adjacency effect (see Section 5.4) and possible influence from emissions occurring on the island.
Remote Sens. 2020, 12, x FOR PEER REVIEW 15 of 41 studies make advantage of these data for inter-comparison, validation, and combined analyses of satellite observations.

Bathymetry
Two regions with bottom depth larger than 300 m at 35.5248°N-12.7667°E and 35.7430°N-12.3579°E, LMP1 and LMP2 respectively, have been identified as possible sites for OC-SVC Infrastructure ( Figure 5). These two sites are characterized by a negligible impact from fishing activities and ship traffic and are located at a reasonable distance from Lampedusa (less than 20 km), thus enabling technical operations and allowing for timely interventions, while avoiding any possible adjacency effect (see Section 5.4) and possible influence from emissions occurring on the island. The bottom effect has been analyzed by performing radiative transfer MC simulations [29,30]. The study acknowledges that the selected region of the Mediterranean Sea is characterized by a relatively high CDOM concentration in comparison with oceanic oligotrophic waters. Light penetration in the water column at the lower end of the spectral region is then significantly limited by the exponential increase of CDOM absorption as wavelength reduces. On the other hand, water absorption lessens the red and near-infrared light transmission. On this basis, 470 nm has been considered as a test wavelength for the bottom effect analysis (specific simulation parameters in Table  2). The bottom contribution has been estimated as percent difference of and values assuming a sandy (S) and a black (B) bottom: where and indicate the bottom and the sensor depth, respectively. The value of has been varied between 50 and 200 m. The value adopted for is 14 m, which corresponds to the depth of the sensor that is mostly affected by bottom reflectance.  The bottom effect has been analyzed by performing radiative transfer MC simulations [29,30]. The study acknowledges that the selected region of the Mediterranean Sea is characterized by a relatively high CDOM concentration in comparison with oceanic oligotrophic waters. Light penetration in the water column at the lower end of the spectral region is then significantly limited by the exponential increase of CDOM absorption as wavelength reduces. On the other hand, water absorption lessens the red and near-infrared light transmission. On this basis, 470 nm has been considered as a test wavelength for the bottom effect analysis (specific simulation parameters in Table 2). The bottom contribution has been estimated as percent difference of L S u and L B u values assuming a sandy (S) and a black (B) bottom: where z btm and z sns indicate the bottom and the sensor depth, respectively. The value of z btm has been varied between 50 and 200 m. The value adopted for z sns is 14 m, which corresponds to the depth of the sensor that is mostly affected by bottom reflectance.  The Harrison-Coombes model has been employed for the sky radiance [46,47]. The values of total irradiance and the diffuse-to-direct irradiance ratio have been determined according to [48]. Inherent Optical Properties (IOPs) to model light transfer in the water column have been selected considering the lowest values of absorption, attenuation and backscattering fraction collected in the field [49]. It is highlighted that the minimum measured b b /b value is 0.021. However, MC simulations have been executed considering a lower b b /b value, namely employing the Volume Scattering Function of Petzold (i.e., b b /b = 0.0183) [50]. The simulated bottom effect is hence expected to overestimate reality due to larger forward scattering. A Lambertian surface with the reflectance value of coral sand is considered for the bottom [51]. Raman scattering [33] has been neglected because it only accounts for a second-order contribution. Examples of optical profiles for a black and sandy bottom are presented in Figure 6 for water depths of 50 m and 200 m. The Harrison-Coombes model has been employed for the sky radiance [46,47]. The values of total irradiance and the diffuse-to-direct irradiance ratio have been determined according to [48]. Inherent Optical Properties (IOPs) to model light transfer in the water column have been selected considering the lowest values of absorption, attenuation and backscattering fraction collected in the field [49]. It is highlighted that the minimum measured bb/b value is 0.021. However, MC simulations have been executed considering a lower bb/b value, namely employing the Volume Scattering Function of Petzold (i.e., bb/b = 0.0183) [50]. The simulated bottom effect is hence expected to overestimate reality due to larger forward scattering. A Lambertian surface with the reflectance value of coral sand is considered for the bottom [51]. Raman scattering [33] has been neglected because it only accounts for a second-order contribution. Examples of optical profiles for a black and sandy bottom are presented in Figure 6 for water depths of 50 m and 200 m.  Table 3 and Figure 7 suggest that the impact of bottom reflectance is negligible for a water column of 200 m or deeper, making the selected site with 300 m depth suitable considering the radiometric accuracy required for SVC (i.e., assuming the representativeness of the parameters employed for the MC simulations).    The characterization of the ocean temperature and salinity in the Lampedusa area has been based on in-situ and satellite observations. Data from a conductivity temperature depth (CTD) at 18 m depth, and from a temperature and pressure sensor at a 1 m depth at the Lampedusa Oceanographic Observatory have been used for the reference period 2018-2019. The annual evolution of the temperature and salinity monthly mean values at 18 m depth is shown in Figure 8. By comparison, the range of temperature and salinity at the MOBY site, derived from Feinholz et al. [16], are also shown in Figure 8. The annual temperature cycle is larger at Lampedusa than at MOBY; conversely, the annual salinity changes are significantly smaller at Lampedusa. These variations constitute essential information for the determination of the water refractive index.

Results presented in
Temperature measurements at a 1 m depth show a slightly larger annual cycle than at 18 m, with similar temperatures in winter and somewhat larger temperatures in summer. The temperature daily cycle is generally within ± 3% of the average daily value.
The Sea Surface Temperature (SST) spatial and temporal variability in the region around Lampedusa has been studied based on the daily Level-4 (SST_MED_SST_L4_NRT_OBSERVATIONS_010_004) product from the Copernicus Marine Environment Monitoring Service (CMEMS) for the reference year 2013. The annual average temperature in the Lampedusa area is about 21 °C. Average day-to-day SST anomalies are generally smaller than 1.5% in the whole area around Lampedusa.

Temperature and Salinity
The characterization of the ocean temperature and salinity in the Lampedusa area has been based on in-situ and satellite observations. Data from a conductivity temperature depth (CTD) at 18 m depth, and from a temperature and pressure sensor at a 1 m depth at the Lampedusa Oceanographic Observatory have been used for the reference period 2018-2019. The annual evolution of the temperature and salinity monthly mean values at 18 m depth is shown in Figure 8. By comparison, the range of temperature and salinity at the MOBY site, derived from Feinholz et al. [16], are also shown in Figure 8. The annual temperature cycle is larger at Lampedusa than at MOBY; conversely, the annual salinity changes are significantly smaller at Lampedusa. These variations constitute essential information for the determination of the water refractive index.
Temperature measurements at a 1 m depth show a slightly larger annual cycle than at 18 m, with similar temperatures in winter and somewhat larger temperatures in summer. The temperature daily cycle is generally within ± 3% of the average daily value.
The Sea Surface Temperature (SST) spatial and temporal variability in the region around Lampedusa has been studied based on the daily Level-4 (SST_MED_SST_L4_NRT_OBSERVATIONS_010_004) product from the Copernicus Marine Environment Monitoring Service (CMEMS) for the reference year 2013. The annual average temperature in the Lampedusa area is about 21 • C. Average day-to-day SST anomalies are generally smaller than 1.5% in the whole area around Lampedusa.

Currents
Data for currents at different depths at the proposed sites have been extracted from the CMEMS Med Sea Physical Reanalysis dataset [52]. Daily data for the period 1987-2013 have been used. The information on currents at the Lampedusa area has been compared with current data for the MOBY site retrieved from the CMEMS GLORYS12V1 global reanalysis (http://resources.marine.copernicus.eu/documents/PUM/CMEMS-GLO-PUM-001-030.pdf).
The reference period that has been considered is 1993-2013. Figure 9 shows profiles of the median, 99th percentile, and maximum values of currents at LMP1, LMP2, and MOBY. For comparison, also values at the Boussole area (Ligurian Sea), also derived from the CMEMS Med Sea Physical Reanalysis dataset, are displayed. Currents are generally weaker at the Lampedusa area than at MOBY and in the Ligurian sea, although higher current may occur in specific events. All sites show a marked seasonal variability for the currents, with median smaller values at Lampedusa in all seasons. Dominant currents are from

Currents
Data for currents at different depths at the proposed sites have been extracted from the CMEMS Med Sea Physical Reanalysis dataset [52]. Daily data for the period 1987-2013 have been used. The information on currents at the Lampedusa area has been compared with current data for the MOBY site retrieved from the CMEMS GLORYS12V1 global reanalysis (http://resources.marine.copernicus.eu/ documents/PUM/CMEMS-GLO-PUM-001-030.pdf). The reference period that has been considered is 1993-2013. Figure 9 shows profiles of the median, 99th percentile, and maximum values of currents at LMP1, LMP2, and MOBY. For comparison, also values at the Boussole area (Ligurian Sea), also derived from the CMEMS Med Sea Physical Reanalysis dataset, are displayed.

Currents
Data for currents at different depths at the proposed sites have been extracted from the CMEMS Med Sea Physical Reanalysis dataset [52]. Daily data for the period 1987-2013 have been used. The information on currents at the Lampedusa area has been compared with current data for the MOBY site retrieved from the CMEMS GLORYS12V1 global reanalysis (http://resources.marine.copernicus.eu/documents/PUM/CMEMS-GLO-PUM-001-030.pdf).
The reference period that has been considered is 1993-2013. Figure 9 shows profiles of the median, 99th percentile, and maximum values of currents at LMP1, LMP2, and MOBY. For comparison, also values at the Boussole area (Ligurian Sea), also derived from the CMEMS Med Sea Physical Reanalysis dataset, are displayed. Currents are generally weaker at the Lampedusa area than at MOBY and in the Ligurian sea, although higher current may occur in specific events. All sites show a marked seasonal variability for the currents, with median smaller values at Lampedusa in all seasons. Dominant currents are from Currents are generally weaker at the Lampedusa area than at MOBY and in the Ligurian sea, although higher current may occur in specific events. All sites show a marked seasonal variability for the currents, with median smaller values at Lampedusa in all seasons. Dominant currents are from the north-west [53] and are associated with the Atlantic-Tunisian current flowing on the Tunisian shelf.

Waves
Data on wave characteristics have been derived from the ECMWF ERA-5 dataset [54]. Data from the period 1979-2019 have been used in the analysis of wave characteristics for both the Lampedusa and MOBY areas. The median of the significant wave height at Lampedusa is about 0.9 m, while it is about 1.1 m in the MOBY area. The significant wave height shows a very limited spatial variability in the region around Lampedusa, much smaller than in the MOBY region. For Lampedusa, the significant wave height is expected to be > 1 m for about 40% of the time. Moreover, relatively high waves are mainly concentrated in winter. Waves from the North-North-West sector are the most frequent, with other smaller contributing sectors from East and South-East.
Significant differences in the wave period exist between the Lampedusa and the MOBY areas. Wave periods are generally shorter at Lampedusa than at MOBY, although a wider range of periods may be found in Lampedusa. The wave characteristics are expected to impact the stability of the buoy attitude; hence, the full wave climate characterization provided with this study aims to offer the backbone of the necessary requirements for the buoys structure design, which is essential for good quality measurements in the long-term.

Inherent Optical Properties
The Inherent Optical Properties (IOPs) were investigated using different the CMEMS datasets, for the January 2017-September 2019 period in the Mediterranean, and for the July 2018-November 2019 for the MOBY area; in-situ measurements of chlorophyll, CDOM and backscatter made at the Lampedusa Oceanographic Observatory at 4.5 m depth in the period June 2018-July 2019 [49].
Overall, a very good agreement between satellite and in-situ determinations of chlorophyll is found and confirms that the central Mediterranean is an oligotrophic region. The chlorophyll amount shows a clear annual cycle, with a maximum in winter and a minimum in summer. The annual average chlorophyll amount is <0.08 mg m −3 in both the Lampedusa and MOBY regions. The day-to-day variability, as retrieved from satellite observations, is <10% at both sites. The chlorophyll spatial variability is also small (<10%) in the Lampedusa area. Figure 10 shows the chlorophyll seasonal cumulative distribution as derived from in-situ observations. The chlorophyll amount exceeds 0.2 mg m −3 only during winter. In-situ continuous observations also show that CDOM and volume backscattering coefficients display an annual cycle, with maxima in winter. Maxima of CDOM and volume backscattering coefficient are 1.1 ppb and 6 × 10 −4 m −1 sr −1 , respectively. Minima are around 0.5 ppb for CDOM and 2.5 × 10 −4 m −1 sr −1 for the backscattering coefficient. Remote Sens. 2020, 12, x FOR PEER REVIEW 20 of 41

Atmospheric properties
In this section, the main atmospheric properties that may affect SVC are presented and analyzed considering ground-based, satellite, and in-situ measurements.

Meteorology and Local Scale Effects
Because of its reduced surface extension and low and smooth topography, the meteorological conditions at Lampedusa are dominated by synoptic-scale phenomena, while local effects are negligible. Weather patterns display a significant seasonal cycle, with a larger variability during autumn and winter, and generally stable and dry conditions during spring and, particularly, in summer. Precipitation is concentrated in autumn and winter, with a maximum in October. Highest temperatures are registered in August (about 26 °C) and minima in January, February, and March (13.5-14 °C).
The annual cycle in weather variability is also evidenced by tropospheric temperature profiles, showing a temperature inversion from April to August, due to the differential warming of the air and sea surface. From September to March, the absence of temperature inversions is linked to more unstable conditions [55].
The comparison of the temporal evolution of meteorological parameters measured at the atmospheric observatory and at the oceanographic observatory confirms that the island of Lampedusa produces a very small perturbation to the atmospheric dynamics, and local effects are negligible.

Wind
Wind speed is linked to wave formation and with the sea surface roughening and foam, thus limiting the availability of optimal conditions for SVC.
Wind speed patterns are driven by the main synoptic structures, showing prevailing directions from NW, with limited interannual variability. The most frequent wind speed is between 4 m/s and 8 m/s (annual average of about 6 m/s); in about 67% of the cases the wind speed values are <8 m/s ( Figure 11).
The wind speed and direction display a seasonal variability, also connected with seasonal changes of the mesoscale Mediterranean atmospheric circulation, with a larger variability during autumn and winter.

Atmospheric Properties
In this section, the main atmospheric properties that may affect SVC are presented and analyzed considering ground-based, satellite, and in-situ measurements.

Meteorology and Local Scale Effects
Because of its reduced surface extension and low and smooth topography, the meteorological conditions at Lampedusa are dominated by synoptic-scale phenomena, while local effects are negligible. Weather patterns display a significant seasonal cycle, with a larger variability during autumn and winter, and generally stable and dry conditions during spring and, particularly, in summer. Precipitation is concentrated in autumn and winter, with a maximum in October. Highest temperatures are registered in August (about 26 • C) and minima in January, February, and March (13.5-14 • C).
The annual cycle in weather variability is also evidenced by tropospheric temperature profiles, showing a temperature inversion from April to August, due to the differential warming of the air and sea surface. From September to March, the absence of temperature inversions is linked to more unstable conditions [55].
The comparison of the temporal evolution of meteorological parameters measured at the atmospheric observatory and at the oceanographic observatory confirms that the island of Lampedusa produces a very small perturbation to the atmospheric dynamics, and local effects are negligible.

Wind
Wind speed is linked to wave formation and with the sea surface roughening and foam, thus limiting the availability of optimal conditions for SVC.
Wind speed patterns are driven by the main synoptic structures, showing prevailing directions from NW, with limited interannual variability. The most frequent wind speed is between 4 m/s and 8 m/s (annual average of about 6 m/s); in about 67% of the cases the wind speed values are <8 m/s ( Figure 11).
The wind speed and direction display a seasonal variability, also connected with seasonal changes of the mesoscale Mediterranean atmospheric circulation, with a larger variability during autumn and winter.

Cloud Cover
The Sicily channel is an area with relatively low cloud cover occurrence, a moderate seasonal cycle ranging from about 10% in summer to 40% in winter and a relatively low (<5%) interannual variability [56].
The detection of cloud-free periods at Lampedusa is based on ground-based measurements collected by all sky imagers and radiometers. A specific algorithm has been developed to identify cloud-free cases (cloud cover <2 oktas, sun unobstructed) from global and diffuse irradiance measurements made with a Multi Filter Rotating Shadowband Radiometer (MFRSR) [57]. The algorithm is based on the scheme proposed by Long and Ackerman [58], which has been modified to identify periods devoid of clouds also during cases with elevated aerosol amounts, in particular Saharan dust. The method is capable to identify cloud-free cases and has been tested against all-sky imager data. In the paper by [59] (Figure 4 of the cited reference) the frequency of occurrence of cloudfree conditions during the central part of the day (UTC time between 09:00 and 14:00), as derived from 15 years of MFRSR measurements, is shown. Cloud free conditions occur in more than 60% of the cases at Lampedusa during June, July, and August, with peaks of 80%. The frequency is <20% in February, November and December. On average, the area has a 37% overall probability of cloud free occurrence; the interannual variability is relatively low.
Cloud properties retrieved from the MODIS sensor on-board Aqua have been investigated over a region of 2° × 2° around Lampedusa analyzing Level-3 monthly mean cloud optical thickness, cloud fraction (CF), and cloud top pressure (CTP) from July 2002 to December 2016 (see Figure 10 [59]).
The annual evolution of cloud fraction, with summer minima, is consistent with that of the occurrence of cloud-free conditions. Moreover, covariance of CF and CTP is observed, indicating lowest clouds and smallest CF in July, when low cloud optical thickness values are also observed; this may be associated to the formation of thin low-level cloud at the top of the atmospheric inversion characterizing the marine boundary layer in summer [55]. The spatial distribution of the cloud cover occurrence is highly homogeneous around Lampedusa throughout the year.
Recently, an algorithm to derive cloud cover from measurements of diffuse irradiance at two wavelengths was developed. This algorithm is based on the method developed by Min et al. [60], which uses measurements of diffuse irradiances at 415 nm and 870 nm obtained with an MFRSR instrument. An improved version of the algorithm was implemented, taking into account the effect of aerosol and varying solar zenith angle and leading to an improved determination of cloud cover and to a more general applicability of the method. High time resolution cloud cover was obtained throughout the year 2016 and compared with ground-based sky imager observations, and cloud cover derived from MODIS satellite data. The analysis of those data has shown the absence of a daily cycle in cloud cover.

Cloud Cover
The Sicily channel is an area with relatively low cloud cover occurrence, a moderate seasonal cycle ranging from about 10% in summer to 40% in winter and a relatively low (<5%) interannual variability [56].
The detection of cloud-free periods at Lampedusa is based on ground-based measurements collected by all sky imagers and radiometers. A specific algorithm has been developed to identify cloud-free cases (cloud cover <2 oktas, sun unobstructed) from global and diffuse irradiance measurements made with a Multi Filter Rotating Shadowband Radiometer (MFRSR) [57]. The algorithm is based on the scheme proposed by Long and Ackerman [58], which has been modified to identify periods devoid of clouds also during cases with elevated aerosol amounts, in particular Saharan dust. The method is capable to identify cloud-free cases and has been tested against all-sky imager data. In the paper by [59] (Figure 4 of the cited reference) the frequency of occurrence of cloud-free conditions during the central part of the day (UTC time between 09:00 and 14:00), as derived from 15 years of MFRSR measurements, is shown. Cloud free conditions occur in more than 60% of the cases at Lampedusa during June, July, and August, with peaks of 80%. The frequency is <20% in February, November and December. On average, the area has a 37% overall probability of cloud free occurrence; the interannual variability is relatively low.
Cloud properties retrieved from the MODIS sensor on-board Aqua have been investigated over a region of 2 • × 2 • around Lampedusa analyzing Level-3 monthly mean cloud optical thickness, cloud fraction (CF), and cloud top pressure (CTP) from July 2002 to December 2016 (see Figure 10 [59]).
The annual evolution of cloud fraction, with summer minima, is consistent with that of the occurrence of cloud-free conditions. Moreover, covariance of CF and CTP is observed, indicating lowest clouds and smallest CF in July, when low cloud optical thickness values are also observed; this may be associated to the formation of thin low-level cloud at the top of the atmospheric inversion characterizing the marine boundary layer in summer [55]. The spatial distribution of the cloud cover occurrence is highly homogeneous around Lampedusa throughout the year.
Recently, an algorithm to derive cloud cover from measurements of diffuse irradiance at two wavelengths was developed. This algorithm is based on the method developed by Min et al. [60], which uses measurements of diffuse irradiances at 415 nm and 870 nm obtained with an MFRSR instrument. An improved version of the algorithm was implemented, taking into account the effect of aerosol and varying solar zenith angle and leading to an improved determination of cloud cover and to a more general applicability of the method. High time resolution cloud cover was obtained throughout the year 2016 and compared with ground-based sky imager observations, and cloud cover derived from MODIS satellite data. The analysis of those data has shown the absence of a daily cycle in cloud cover.

Aerosols
Aerosol transport patterns and optical properties over Lampedusa have been extensively studied using ground-based measurements and satellite data (e.g., references [57,[61][62][63]). Measurements of the aerosol optical properties were started at Lampedusa in 1999 using an MFRSR instrument; a CIMEL sunphotometer contributing to AERONET was set up in 2000, although it worked with some discontinuity during the first years of operation. An extended dataset combining these two datasets has been implemented, and a correction scheme has been developed to take into account differences in the field of view of the two instruments [63].
Several studies show that, except for the cases directly influenced by Saharan dust transport events, Lampedusa is characterized by low aerosol optical depth (AOD). When such events are not considered in the statistics the annual mean aerosol optical depth at 870 nm is about 0.09 (as derived from the data by [61]). The distribution of the AOD daily average values at 870 nm for the period 1999-2018 (> 4500 days with measurements) is shown in Figure 12: about 75% of all AOD values at 870 nm are below 0.15, while larger AOD values are measured during sporadic events of Saharan dust transport that have a seasonal cycle and typically last for one or two days [57]. The distribution of Ångström exponent values, which is linked to the particles' dimension, is shown in Figure 12. Cases of pollution transport from the European continent, characterized by small particles, i.e., larger Ångström exponent values [61], are infrequent; values of Ångström exponent below 1 represents 78% of all values.

Aerosols
Aerosol transport patterns and optical properties over Lampedusa have been extensively studied using ground-based measurements and satellite data (e.g., references [57,[61][62][63]). Measurements of the aerosol optical properties were started at Lampedusa in 1999 using an MFRSR instrument; a CIMEL sunphotometer contributing to AERONET was set up in 2000, although it worked with some discontinuity during the first years of operation. An extended dataset combining these two datasets has been implemented, and a correction scheme has been developed to take into account differences in the field of view of the two instruments [63].
Several studies show that, except for the cases directly influenced by Saharan dust transport events, Lampedusa is characterized by low aerosol optical depth (AOD). When such events are not considered in the statistics the annual mean aerosol optical depth at 870 nm is about 0.09 (as derived from the data by [61]). The distribution of the AOD daily average values at 870 nm for the period 1999-2018 (> 4500 days with measurements) is shown in Figure 12: about 75% of all AOD values at 870 nm are below 0.15, while larger AOD values are measured during sporadic events of Saharan dust transport that have a seasonal cycle and typically last for one or two days [57]. The distribution of Ångström exponent values, which is linked to the particles' dimension, is shown in Figure 12. Cases of pollution transport from the European continent, characterized by small particles, i.e., larger Ångström exponent values [61], are infrequent; values of Ångström exponent below 1 represents 78% of all values.
Chemical analyses of in-situ aerosol samples [64] show that natural sources (sea salt, mineral dust, and biogenic emissions) give the largest contribution to PM10 all year round. Sea salt represents the largest absolute and relative contribution in all seasons, with a maximum (54%) in winter and minima (33%) in summer, while desert dust contributes by 17%-37% (largest in autumn).

Interfering Gases
Among all the atmospheric gases with absorption in the UV-visible-near infrared spectral intervals, water vapor, ozone, nitrogen dioxide, and sulphur dioxide have the largest absorption.
Total ozone column (TOC) has been measured at Lampedusa since 1998 by means of the Brewer spectrophotometer #123. TOC time series in the period 2000-2016 has a clear annual cycle with maxima in spring and minima in autumn, which is the result of a balance between transport and photochemical loss, whereby the former is dominant in winter, while photochemical loss dominates Chemical analyses of in-situ aerosol samples [64] show that natural sources (sea salt, mineral dust, and biogenic emissions) give the largest contribution to PM10 all year round. Sea salt represents the largest absolute and relative contribution in all seasons, with a maximum (54%) in winter and minima (33%) in summer, while desert dust contributes by 17%-37% (largest in autumn).

Interfering Gases
Among all the atmospheric gases with absorption in the UV-visible-near infrared spectral intervals, water vapor, ozone, nitrogen dioxide, and sulphur dioxide have the largest absorption.
Total ozone column (TOC) has been measured at Lampedusa since 1998 by means of the Brewer spectrophotometer #123. TOC time series in the period 2000-2016 has a clear annual cycle with maxima in spring and minima in autumn, which is the result of a balance between transport and photochemical loss, whereby the former is dominant in winter, while photochemical loss dominates in summer, leading to a minimum in autumn. The TOC mean annual cycle at Lampedusa obtained from daily average has an amplitude of about 60 DU. The TOC diurnal variability is mainly caused by photochemical and transport processes occurring in the troposphere. At Lampedusa, diurnal variability is limited, below 2.5% under cloud-free conditions.
The total NO 2 column is not measured at Lampedusa, so data from the Ozone Monitoring Instrument satellite (OMI) onboard Aura satellite have been examined. The temporal evolution of NO 2 column density provided by OMI measurements have been derived in the period 2010-2018 from the NASA Giovanni web site (https://giovanni.gsfc.nasa.gov/giovanni/). Maxima are reached in summer and minima in winter, with an amplitude of the annual cycle of about 2 × 10 15 cm −2 . The day-to-day column NO 2 variability is <17%.
Finally, in-situ measurements of O 3 and NO 2 mixing ratio show that their diurnal variability is negligible.

Adjacency Effect
The contribution of the diffuse radiance reflected by the presence of land to the radiance measured by the satellite for an ocean scene, called the adjacency effect (AE), perturbs the spectral satellite data, thus leading to uncertainties in derived primary products. The AE has to be minimized by considering a minimum distance from the coast (i.e., [12]).
The adjacency radiance contribution depends on the land surface radiative properties (spectral and geometric), on the extension and morphology. The adjacency effect field around Lampedusa island has been studied by [65] simulating OLCI measurements at the central wavelengths of four OLCI channels (λ = 490 nm, 555 nm, 670 nm, and 865 nm), accounting for typical surface and atmospheric conditions derived from satellite or remote sensing surface observations.
According to the results of this study, AE is avoided at distances larger than 15 km from the coast of the island.
Since land surface albedo has a key role in the estimation of the AE, spectrally detailed measurements of this quantity should be made at the same time and at the same wavelength as OC measurements [65]. The estimation of the AE will take advantage of surface albedo measurements carried out at the Atmospheric Observatory. While spectral downward irradiance has been measured by spectroradiometers since 2013, measurements of reflected spectral irradiance are about to be started.
Measurements of broadband solar albedo at the Atmospheric Observatory started in 2016. The annual distribution of the surface albedo monthly mean values presents maxima of 0.19-0.21 from July to January, and a decrease to 0.16-0.17 in April and May, when the occurrence of rain favors the growth of vegetation, decreasing the reflecting properties of the surface.

Uncertainty of the Radiometric Measurement
The estimated uncertainty budget (EUB) of the infrastructure is firstly assessed for the radiometric measurements. For each radiometric correction, the measurement equation is implemented in the ground segment to provide end-to-end uncertainties of every individual measurement ( Figure 13). Uncertainty propagation is then theoretically achievable with this formalism (Type A uncertainty, defined in a statistical sense), using either the Guide to the Expression of Uncertainty in Measurement (GUM) first-order assumption [66] or Monte Carlo simulations [67].
For the quantified EUB of the optical system, we prefer, however, to rely on the actual uncertainties recently assessed by the MOBY team with the MarONet optical system [17,18]. This refers to Type B uncertainties, i.e., uncertainties not derived by a statistical analysis, but field experiences, calibration reports, and background knowledge on the system. It shall be noted that uncertainties of MarONet are reduced compared to the previous MOBY system [68]. The EUB of the optical system is split into three main parts:
Uncertainty of the downwelling irradiance (Ed) covering the optical system and environmental effects ( Table 6); 3.
Remote Sens. 2020, 12, x FOR PEER REVIEW 24 of 41 1. Uncertainty of water-leaving radiance (Lw(0+)) covering the optical system (Table 4), deployment structure and ground segment processing (Table 5); 2. Uncertainty of the downwelling irradiance (Ed) covering the optical system and environmental effects (Table 6); 3. Uncertainty of the resulting remote-sensing reflectance, Rrs (Table 7).  A basic principle in the uncertainty budget is that corrections are multiplicative. The formalism adopted in the study follows the protocols of the IOCCG report "Protocols for Satellite Ocean Color Data Validation: In-situ Optical Radiometry" [28]. Denoting with the symbol (λ) either L(λ) or E(λ) in physical units (W·cm −2 ·sr −1 ·nm −1 and W·cm −2 ·nm −1 , respectively), the conversion equation for a given correction, ℵ i (λ), from measurement i−1 to measurement i , is The term ℵ i (λ) accounts for the non-ideal instrument response of the instrument and implicitly depends on various input parameters (possibly other wavelengths, like for straylight correction). The final radiometry is thus given by the product of all corrections: where 0 represents the net radiometric data obtained by subtracting the dark contribution from signal data after scaling the raw CCD readouts by the measurement integration time and the bin factor. Applications of the GUM gives This justifies why the quadratic uncertainty budget is additive, in relative unit, in the next tables. Furthermore, the tables tentatively distribute the uncertainties sources between their random and systematic parts. This has an important consequence on the SVC gains, further discussed in the next section.  Table 4 gives the EUB of the optical system (Lu), and Table 5 gives the EUB of the deployment structure (Lu in the field) and of the depth extrapolation and post-processing up to Lw(0+). Most of the terms are from Voss et al. 2017 and Johnson et al. 2017 [17,18]. Uncertainty due to spectral resolution is covered in the next section. The five sources of uncertainty related to radiometric calibration refer, respectively, to the sphere calibration, sphere drift, non-uniformity, calibration measurement uncertainty and spectral interpolation at the radiometer bands. The effect of angular responses is neglected thanks to very limited FOV (1.73 • full angle). The uncertainty of the immersion factor comes from Feinholz et al. 2017 [16]. Thermal stability refers to terms "instantaneous temperature" in Johnson et al. 2017 [18]. Dark current is implicitly taken into account in the temporal averaging and random component of noise calibration below. Uncertainty related to polarization comes from Voss et al. 2017 [17] and is very conservative, considering the use of a depolarizer. Noise characterization covers MOBY integration time and reproducibility of Johnson et al. 2017 [18]. The effect of environmental conditions on Lu is assessed with a very conservative assumption as the sum of uncertainties related to MOBY stability during deployment-system response, in-water internal calibration, and wavelength stability [68].
The quantified uncertainty of the depth propagation at MOBY of 0.4% [17] is kept for EURYBIA. Indeed, with chlorophyll concentration below 0.15 mg/m 3 most of the time relevant for SVC (only up to 0.2 in January), the attenuation coefficient at Lampedusa is below the maximum considered in Voss et al. 2017 [33]; real acquisitions at Lampedusa done with a Satlantic profiler (with OCR507 radiometers) in 2017 confirms values of K L below 0.04 m −1 from 412 nm to 510 nm, of about 0.06 m −1 at 555 nm and 0.22 m −1 at 665 nm, fully compatible with the analysis of Voss et al. 2017 [33] at MOBY (see their Figure 6). A conservative uncertainty of 1% is assigned to the biofouling correction [68].
The total uncertainty of Lw(0+) ranges from 3% to 5.5%, with spectral variation essentially due to spectral emission of the calibration lamp and noise level.
EUB of solar irradiance (Es) is given in Table 6. The uncertainty induced by the cosine correction depends on the solar zenith angle and the wavelength. For relevant conditions used in the SVC (most of the time θ 0 < 60 • ), the uncertainties range from below 0.1% in the blue and green bands to 0.2% in the red. Uncertainty of radiometric calibration accounts for these five sources respectively [18]: FEL lamp calibration, drift, bench effect, calibration measurement uncertainty, and spectral interpolation. The EUB includes also an estimated 2% environmental uncertainty, yielding a total uncertainty between 2% and 3% at all bands. Table 7 eventually gives the uncertainty of remote-sensing reflectance (Rrs), through simple combination between Lw and Es. We emphasize that the requirement to conduct SVC with Lw rather than Rrs = Lw/Es is not agreed by all agencies; using Rrs may increase the total EUB and is here considered as a conservative assumption. The total uncertainty remains below 5% from 412 nm to 560 nm and below 6% in the red. These numbers come from various conservative assumptions and provide more an upper threshold than optimal estimates for the EURYBIA concept. Slightly smaller uncertainties can likely be expected during the detailed design of the infrastructure.

Uncertainty of the SVC Gains
The second level of the EUB refers to the SVC process, i.e., uncertainty of the vicarious gains. Although gain computation is not part of the SVC infrastructure, we consider that it is of highest importance to quantify the end-to-end uncertainty of the system. This computation is done for the identified locations around Lampedusa as well as for the MOBY location in Hawaii (MOBY), and using OLCI data and effective criteria defined by EUMETSAT, detailed below. Because both Lampedusa candidate locations yield very similar numbers, in the following we only report the analysis for location LMP1 (East of Lampedusa). The EUB contains two parts, both based on realistic assumption on the number of match-ups and radiometry at the calibration site.

1.
Ground segment post-processing for in-situ and satellite data merging in the gain computation, i.e., match-ups, spectral integration (

2.
Gain computation, up to final mission average gains over a decade ( Figure 14).
We cover two years of Sentinel-3 OLCI data (Full Resolution products) from July 2017 to July 2019. We apply the strict criteria used by EUMETSAT for SVC gain generation, assuming that the criteria related to the in-situ measurement should be met with an operational infrastructure: 1.
For each pixel, sensor zenith angle should be below 56 • and the sun zenith angle below 70 • ; 3.
All pixels within the 5 × 5 region of interest (ROI) should be valid; 5.
Mean chlorophyll within the ROI should be below 0.2 mg/m 3 ; 6.
Mean aerosol optical depth (AOD) for the channel at 865 nm should be below 0.15; 7.
Coefficient of Variation CV = σ * µ * should be below 0.15 for remote sensing reflectances at bands between 412 and 560 nm, where µ* and σ* are mean and standard deviation, respectively, computed after the removal of pixel outliers. A pixel is considered as an outlier for variable x i if |x i -µ| < 1.5 σ where µ is the mean and σ is the standard deviation of that variable over the 5 × 5 box.
In two years, the number of valid data is 36 at Lampedusa and 64 at Hawaii ( Table 9). The difference is essentially explained by the threshold on AOD, and to a lesser extent on chlorophyll concentration, noticing that the number of match-ups after the flag and geometry screening only is exactly similar (70 match-up in two years for both sites). However, recent analyses suggest that OLCI may overestimate AOD [69], possibly reducing the number of cases for valid matchups at Lampedusa.
In the uncertainty of the ground segment post-processing, our analysis includes the spatial variability of the satellite observation (the first item of Table 8). This uncertainty is quantified for every match-up in each relative unit by the Coefficient of Variation (CV) and strongly grows toward red bands for clear waters. This term appears to be crucial in the EUB; neglecting it would drastically reduce the final uncertainty on the SVC gains in the NIR. For a given band λ in the visible part of the solar spectrum, the SVC of OLCI consists first in computing individual gains, g(λ), for every match-up and pixels [70]: where ρ t w (λ) is the in-situ marine reflectance (de-normalise for BRDF effect, i.e., in the observation and solar geometry of the actual pixel), ρ path (λ) and t(λ) are the atmospheric path reflectance and diffuse total transmittance, respectively, accounting for the Rayleigh and aerosol components reflectance and ρ gc (λ) is the actual TOA reflectance pre-corrected for gaseous absorption and sun glint effect. Note that computation and application of gains directly on the TOA Level-1 radiometry would not change the present analysis. When dealing with standard atmospheric correction [37,71], the aerosol amount and type are identified by the near-infrared (NIR) bands only, assumed to be calibrated beforehand by other means [70]. Application of the GUM first-order assumption [66] gives the uncertainty of individual gain: This equation implicitly discards the uncertainty of the satellite radiometry before SVC, i.e., u ρ gc , because it essentially depends on the satellite characteristics (radiometric noise) and not on the OC-SVC infrastructure nor aerosol or marine conditions. In term of covariances, we do not expect either any correlation with ρ gc , insofar as the standard atmospheric correction does not detect aerosol from the VIS domain: r ρ gc , ρ path =0 and r ρ gc , t =0. On the contrary, covariance between scattering functions exists and is included through u ρ path , t . The multiplicative factor tρ t w /ρ gc decreases in absorbing waters and theoretically minimizes the uncertainty of g. However, complex waters may counterbalance this decrease by an increased uncertainty of the additive components in the square brackets-less reliable in-situ marine reflectance (temporal variability, radiometric measurement) and possibly more complex atmosphere. Also, a lower ratio gives more importance to the atmospheric uncertainty through the factor ρ path /tρ t w . The balance between all terms is important to get the realistic uncertainty estimate. Eventually, when the uncertainties of individual gains are computed, the uncertainty on the mission average gains g is simply given by This formula assumes a simple average of the individual gains. More evolved averaging techniques are also possible, such as the mean of the semi-interquartile range to remove outliers, or a weighted average. In the present context we have no access to the distribution of the gains, only to their uncertainties, so the arithmetic averaging is kept as a conservative estimate. The value of u(g) is affected by the number of years considered in the time-series, here Y = 2. A more objective metric is the Relative Standard Error of the Mean (RSEM), introduced by Zibordi et al. 2015 [4] to study the stability requirements of SVC gains per decade. With an average SVC gain of the order of unity, RSEM, can be related to our computed u(g) by simply scaling the number of match-ups obtained in Y years to a decade.
The relative uncertainty of in-situ marine reflectance, , is assumed to be of 5% at all bands and all calibration sites (Lampedusa, Hawaii). This assumption, previously used in Zibordi et al. 2015 [4], is consistent with the overall quantification detailed in the previous sections. It has obviously a direct impact on the final uncertainties and could be refined in a later stage. The purpose of the present analysis is to simulate the gain uncertainties that would be achieved with OLCI over these locations assuming high-quality and daily radiometric measurements. All other terms involved in the computation of u(g) are assessed by satellite data, notably the uncertainties of the atmospheric scattering functions, u ρ path and u(t). This allows conducting a harmonized uncertainty budget for all sites under consideration, while only Lampedusa is fully characterized by long-term AERONET measurements. The choice is made here to fully rely on OLCI aerosol products, as it would be done in real SVC operation. The actual uncertainties u ρ path and u(t) for a given match-up used in the SVC process may be due to many sources, such as the aerosol modelling in the atmospheric correction, the aerosol identification in the NIR, the propagation in the visible, and include both systematic and random components. Because the purpose of SVC is to remove the systematic effects, the present analysis determines u ρ path and u(t) by the local variation of ρ path and t around the calibration site. This approach considers that heterogeneity in the atmospheric functions degrades the confidence in the gains. The method is implemented as follows: 1.
Consider, for each potential match-up, a ROI of 5 × 5 OLCI pixels (full resolution), consistent with EUMETSAT protocols for SVC as defined previously; 2.
Consider all OLCI metadata and ancillary data at pixel level (wind speed, pressure, geometry); 3.
Identify, for each pixel, the best matching aerosol model in the OLCI LUTs; 5.
Given geometry and ancillary data, compute ρ path and t by the OLCI LUTs (Rayleigh + aerosol), as done in OLCI atmospheric correction; 6.
Compute u ρ path and u(t) by their local standard-deviation over the ROI.
Final values of RSEM are given in Figure 14-the top-left for Lampedusa (LMP1 location) and the top-right for the Hawaii location. RSEM is of about 0.05% at all bands and both sites, with a peak closer to 0.08% in the Mediterranean at 490 nm. These numbers are consistent with the stability requirement over a decade in the blue-green bands but do not reach the requirement of 0.005% in the red. It is worth noting that this target is hardly reached by any existing SVC infrastructure. Indeed, the RSEM values computed by Zibordi et al. 2015 [4] based on real SVC gains with real in-situ measurements shows that only MOBY gains reach the required RSEM of 0.05% in the blue bands, and not in the red (their Figure 4, reprinted here in Figure 14, bottom-right). Despite the consideration of spatial variability in our uncertainty estimate (more than 50% at 674 nm, Table 8), this strongly suggests that the actual uncertainties of marine reflectance are, for many infrastructures, largely above 5%, especially in the red bands for clear waters.
Another analysis has been conducted in parallel following [4], by simplifying the uncertainty computation of individual gain by u(g) = tρ t w /ρ gc * u ρ t w /ρ t w , where the overline stands for a temporal average. At Lampedusa, this average factor is of 8.4% at 412 nm, 10.7% at 443 nm, 12.4% at 490%, 5.6% at 560 nm, and 1.2% at 674 nm. While neglecting the atmospheric uncertainties, this formulation allows easy computation of the uncertainty of gains, hence of RSEM, for testing various distributions of in-situ uncertainty sources between their random and systematic parts. Here two cases have been investigated:

1.
A purely random case: the uncertainty sources are all assigned to the random component.

2.
A mix between random and systematic sources (Figure 14, bottom-left, triangles): in Tables 4-8, an attempt is made to assign the uncertainty of various sources to the systematic component, only (i.e., without any random counterpart). Note that the global uncertainty of a given source is identical to the previous case (i.e., the quadratic sum of the random and systematic components equals the square random component of the previous case). The first case is an ideal condition, that reduces the uncertainty of mission-average gains as a function of the square root of the number of match-ups. Such an assumption is generally retained in the SVC literature: Antoine et al 1999 [70] and Zibordi et al. 2015 [4] have estimated the final SVC uncertainty through the standard-error of the mean or the relative standard error of the mean of SVC gains, what amounts to dividing the standard-deviation of individual gains by the square root of number of match-ups. Under this purely random assumption, the EURYBIA assessment retrieves an uncertainty of about 0.05% over the full spectral range (Figure 14, bottom-left with circles) which is The first case is an ideal condition, that reduces the uncertainty of mission-average gains as a function of the square root of the number of match-ups. Such an assumption is generally retained in the SVC literature: Antoine et al 1999 [70] and Zibordi et al. 2015 [4] have estimated the final SVC uncertainty through the standard-error of the mean or the relative standard error of the mean of SVC gains, what amounts to dividing the standard-deviation of individual gains by the square root of number of match-ups. Under this purely random assumption, the EURYBIA assessment retrieves an uncertainty of about 0.05% over the full spectral range (Figure 14, bottom-left with circles) which is extremely similar to the assessment of Zibordi et al. 2015 [4] for MOBY (reprint in Figure 14, bottom-right, blue curve for MOBY). The perfect consistency between these two assessments, made under completely different datasets and methods, is noteworthy-the EURYBIA computation relies on the uncertainty budget of the optical system (MarONet quantification) over Lampedusa conditions observed by OLCI, while the Zibordi et al. 2015 [4] estimate relies on the statistical analysis of real SVC gains computed at MOBY for the SeaWiFS sensor. These results demonstrate the relevance of EURYBIA, based on both the state-of-the-art optical system and favorable conditions, while other systems cannot meet such performance (Figure 14, bottom-right, other systems than MOBY).
This purely random assumption may, however, not be strictly true for various sources of uncertainty. For instance, it is likely that the radiometric calibration, stray light correction, depth extrapolation, have systematic uncertainties. The resulting uncertainty of average gains is then not driven by the number of match-ups, but by these systematic contributions which do not cancel out in the averaging process. Uncertainties can now reach nearly 0.25% in the blue and 0.1% in the red (Figure 14, bottom-left with triangles). Importantly, although these numbers show some similarity with the curves of Zibordi et al. 2015 [4] for non-MOBY systems, they are not comparable in the sense that the relatively poor performance of these latter systems still benefits from the ideal random assumption. We also emphasize that the "rand.+syst." case refers to the worst case scenario, since all sources should in reality have a random part. The distinction between both components is currently not addressed [18,33] and should be investigated in the next phase. The reality probably lies between the ideal random assumption and the worst-case scenario. Quantifying the detailed structure of the uncertainties is a key scientific challenge that should be addressed by coordinated effort among SVC teams and metrology institutes part of the future OC-SVC network.
In summary, the present analysis based on two years of OLCI data demonstrates that the Lampedusa location (either east or north-west, the latter not shown here) is theoretically suitable for an SVC site, in comparison with the existing site in Hawaii. Clearly, the screening criteria currently defined for an ideal location such as Hawaii removes a lot of data in the Mediterranean Sea (about 48% removal at Lampedusa versus less than 9% at Hawaii, starting from the number of valid data after flag and geometrical screening). However, the amplitude of the marine signal at Lampedusa is below that of the oligotrophic waters of Hawaii (while not being as low as for mesotrophic waters), what advantageously drives the uncertainty propagation from sea level to the SVC gains (term tρ t w /ρ gc ).

Discussion and Conclusions
In this paper, we presented the overall design of the EURYBIA OC-SVC infrastructure, taking into account the Copernicus long-term and multi-mission perspective to deliver OC products and services.
Our strategy, to reach the goal to deliver state-of-the-art Fiducial Reference Measurements (FRM) in a very tight development schedule (expected operational target after 36 months from the beginning of the development), is to operate independently but in compliance with MarONet specification [72]. This solution brings to a network concept (e.g., AERONET-OC, [73]) and takes full advantage of more than 30 years of MOBY experience in operation and technological developments. EURYBIA will be compatible with the MarONet network in terms of radiometric data collection, traceability, and protocols, but autonomous regarding detailed architecture, development, funding, and operations.
Our proposed solution for the European OC-SVC Infrastructure (Figure 1) offers the possibility to develop and deploy a technologically proven system within a realistic timeframe for the Copernicus needs. Moreover, it allows to avoid the compatibility issues experienced in the past two decades due to the adoption of different instrumentation and processing approaches. The MarONet optical system, constituting the core instrumentation of the infrastructure, has been verified to be fully compliant with the radiometric and spectral requirements ([3], see Table 1). The EURYBIA Field Segment accounts for the instrumentation, infrastructure and operations required to acquire radiometric FRM in the field. The field segment includes also an on-shore laboratory for preand post-deployment instrument calibration, characterization, and verification of the calibration sources' stability. Complete SI traceability will be insured by a Metrology Institute. The EURYBIA Ground Segment communicates with the buoys to transmit/receive data and to monitor the infrastructure status ensuring operations on 24/7/365. The processing unit computes data products with their associated uncertainties and performs the quality assurance of results. Data products are distributed within the official Copernicus mechanism through the dissemination unit, which also gives access to the open-source processing code and documentation.
The EURYBIA infrastructure is proposed to be settled in the Strait of Sicily offshore of Lampedusa Island, one of the potential European OC-SVC sites [11][12][13]. Lampedusa is a small flat island in the central Mediterranean Sea far away from the mainland and main traffic routes, but very well connected to Europe. The island is surrounded by oligotrophic waters and hosts a Marine Protected Area.
Our analysis, based on satellite and in-situ observations, revealed the occurrence of consistent and homogeneous cloud free conditions (> 60% during summer). The median value of the AOD at 870 nm is <0.1, hence below the currently adopted threshold value (0.15) for acceptance of measurements for vicarious calibration. Sporadic larger values of AOT are observed only during desert dust transport events characterized by relatively short time duration (1-2 days) over the whole Mediterranean region [74].
Being Lampedusa a very small island characterized by low and smooth topography, its meteorological conditions are dominated by synoptic scale phenomena, and local effects are negligible, as demonstrated by the absence of a diurnal cycle in cloudiness and aerosol load. Wind patterns are driven by the main synoptic structures, showing prevailing directions from NW, with a limited interannual variability and most frequent wind speed between 4 m/s and 8 m/s. Wave and current conditions are suitable to ensure the required stability of the optical buoy.
The analysis of the optical characteristics of the waters surrounding the Island were compared with those observed at the MOBY site. This reveals that in Lampedusa the annual average chlorophyll is <0.08 mg m −3 , with a day-to-day and local spatial variability below 10%. An end-to-end uncertainty budget based on OLCI data adopting EUMETSAT screening protocols for SVC and site characterization demonstrates the relevance of this site. The overall uncertainty of EURYBIA SVC gains is estimated to be about 0.05% per decade in the blue-green wavelengths. This value is comparable with those estimated for the reference site in Hawaii and in compliance with stability requirements for climate studies in this spectral domain.
We propose to settle the location site at more than 15 km off the island of Lampedusa at suitable sea depth conditions (> 300m) to avoid adjacency and bathymetry effects. Two candidate locations ( Figure 5), satisfying these geophysical criteria and relatively distant from the main marine traffic routes, have been identified-LMP1 (35.5248 • N-12.7667 • E) and LMP2 (35.7430 • N-12.3579 • E). Both candidate sites have been evaluated to be fully compliant with the environmental (e.g., cloudiness, aerosols, surface wind, currents, and waves) and logistical requirements for vicarious calibration identified by Mazeran et al. [3]. Recently, Bulgarelli and Zibordi [13] carried out a complete study on the adjacency effect of the Lampedusa Island to address implications on a hypothetical nearby system vicarious calibration infrastructure for satellite ocean color sensors. Their results, as expected, confirmed that the site should be located at distances larger than approximately 14 km from the coast (requirement satisfied by both LMP1 and LMP2). They also suggested that the region in the direction of the reflected sunbeam should be avoided. Based on this latter recommendation LMP1 location should be preferred. The final selection of the exact location where to deploy the optical buoy will require further studies on logistics and environmental characterization, including analyses of ad-hoc in-situ observation to be carried out around the Island. This will be part of future studies to be carried out in phase 3 of the EUMETSAT Programme.
It is important to underline that the EURYBIA design is modular-any update of individual subsystems will be possible without the need to reconsider the whole system architecture. The proposed infrastructure will be able to host, on the moored buoy, additional instrumentation that can be used also for CAL/VAL of other satellite missions. The scheduled cruises for infrastructure maintenance will give an additional opportunity to collect in-situ measurements needed to improve the current OC product algorithms as well as to develop new ones.
Our EURYBIA design is independent from the proposed installation site, so other potential locations for its deployment could be considered. A preliminary analysis of other potential sites was carried out during the EUMETSAT project (not shown) revealing that Lampedusa is a very good candidate for the installation of the European OC-SVC Infrastructure, also in comparison to other locations. A detailed evaluation of other potential European sites will require a deep inter-comparison of environmental characteristics, including a complete evaluation of end-to-end uncertainty budget as we carried out for Lampedusa and MOBY. This will be part of our future studies.
In conclusion, the logistics, the capability and required expertise to operate the infrastructure will be ensured by the fact that Lampedusa has been hosting a long-term well-established Climate Observatory for more than twenty years. In addition, the presence of the Climate Observatory has the advantage to make available relevant observations on atmospherics variables, including aerosols profiles that could be used to improve the accuracy of system vicarious calibration. Lampedusa, as proposed site, represents a trade-off between ideal conditions reached in Hawaii, for the reference MOBY infrastructure, and the priority to host the EURYBIA system in European waters. Marullo are thanked for making available temperature and salinity data from the Lampedusa Oceanographic Observatory. Analyses of OMI NO 2 data used in this paper were produced with the Giovanni online data system (https://giovanni.gsfc.nasa.gov/giovanni/). developed and maintained by the NASA GES DISC. This study has been conducted using EU Copernicus Marine Service Information.

Conflicts of Interest:
The authors declare no conflict of interest.