Radiometric Top-of-Atmosphere Reﬂectance Consistency Assessment for Landsat 8 / OLI, Sentinel-2 / MSI, PROBA-V, and DEIMOS-1 over Libya-4 and RadCalNet Calibration Sites

: There is a clear trend toward the use of higher spatial resolution satellite sensors. Due to the low revisit time of these sensors and frequent cloud coverage, many applications require data from di ﬀ erent sensors to be combined in order to have more frequent observations. This raises concerns regarding data interoperability and consistency. The initial pre-requisite is that there are no radiometric di ﬀ erences in top-of-atmosphere (TOA) observations. This paper aims to quantitatively assess di ﬀ erences in the TOA signal provided by PROBA-V, Sentinel-2A and Sentinel-2B, Landsat-8, and Deimos-1 by using observations over both the Libya-4 desert calibration site and the RadCalNet sites. The results obtained over the Libya-4 site indicate that for all sensors investigated, the inter-sensor deviations are negligible, i.e., within ± 2% for comparable spectral bands, with the exception of the Deimos-1 Green band. Clear BRDF (bi-directional reﬂectance distribution function) e ﬀ ects were observed over the RadCalNet sites, thereby preventing consistent conclusions on inter-sensor deviations from being made. In order to fully explore the potential of the RadCalNet sites, it is recommended that BRDF characterizations be additionally incorporated into the RadCalNet simulations and made publicly available through the distribution portal.


Introduction
High spatial resolution satellite missions enable local operational earth observation (EO) services to be established. However, due to frequent cloud coverage, combined with the lower revisit time of higher spatial resolution sensors, operational services have to combine data from different sensors in order to collect observations of sufficient frequency. A seamless combination of EO data products obtained from different missions is not straightforward. The joint use of data from different sensors raises some clear concerns regarding data consistency. One of these is that there might be radiometric differences at TOA (top-of-atmosphere) level. Furthermore, even for a perfectly cross-calibrated constellation of instruments, intrinsic differences in the relative spectral response functions (RSRF) of spectrally comparable bands might cause discrepancies in the final products. Biases in the data products may also be introduced through the use of different processing chains and algorithms (e.g., for the atmospheric correction). One option is to normalize time series at product level, following an empirical approach. The drawback of such an approach is that the root cause of an observed bias is not identified. Furthermore, a non-linear empirical relationship could be identified with associated high uncertainty. Lastly, such an empirical relationship is only applicable to the sites and/or periods for which it has been derived.
In the harmonization of multi-mission high-resolution time series: Application to BELAIR (BELHARMONY) project [1], the aim is to obtain harmonized time series from Level 1 (L1) through value-added Level 3 (L3) bio-geophysical products for Landsat 8 (LS8) operational land imager (OLI), Sentinel-2 A/B (S2A/S2B) multi-spectral instrument (MSI), Project for On-Board Autonomy-Vegetation (PROBA-V), and Deimos-1 within the VNIR-SWIR spectral range. In BELHARMONY, a bottom-up approach (L1-L3) is used to assess and improve the consistency of multi-mission high-resolution time series. This bottom-up approach starts from the L1 TOA observations. This is important, as differences at L1 can be strongly amplified into L2 and L3 data products and should therefore already be corrected for at L1. Vicarious radiometric calibration methods are used to assess the consistency of the L1 TOA observations. Following the CEOS (Committee on Earth Observation Satellites) definition, "Vicarious calibration" refers to "techniques that make use of natural or artificial sites on the surface of the Earth for the post-launch calibration of sensors. These targets are imaged in near-coincident fashion by the sensor to be calibrated and by one or more well-calibrated sensors from satellite or aircraft platforms or on the ground".
The simultaneous nadir overpass (SNO) method [2] enables inter-sensor L1 cross-calibration. As highlighted in [3], an important limitation of the SNO approach is the need for simultaneous nadir observations of the sensors over the same target and strict constraints on the solar zenith and viewing zenith angles (SZA and VZA, respectively). Therefore, only a very low number of eligible image pairs are obtained, which becomes an even greater issue if more than two sensors are to be cross-calibrated, as in this case.
Another method of inter-sensor cross-calibration relies on the comparison of the sensor-measured TOA spectral reflectance ("MEAS") against modeled TOA ("REF") reflectance for a number of observations.
Ground-based reference surface reflectance spectra combined with atmospheric information, measured simultaneously with the satellite data acquisition, can be used as input to radiative transfer models (RTMs) in order to simulate/model the reference TOA reflectance values ("REF"), to which the sensor-measured values ("MEAS") are compared. Since the same reference approach is used for all instruments under investigation, the difference in the "MEAS/REF" ratio between the instruments is a direct measure of the inter-sensor calibration.
To increase the number of matchups between in-situ measurements and space sensor observations and to reduce the overall uncertainties, land equipped sites (LES) have been established. In-situ measurements are fully automated at these LES, thereby negating the presence of a field crew. To ease data access from LES, the RadCalNet, a radiometric calibration network of automated instruments, was established by a dedicated working group, with members from various international (space) organizations. RadCalNet  . The La Crau site, among others, has been used for vicarious calibration/validation of the modular optoelectronic multispectral scanner (MOMS)-2P-camera onboard the Russian MIR space station, the high-resolution visible (HRV) cameras on the Satellite Pour l'Observation de la Terre-Végétation (SPOT) satellites and Sentinel-2. In addition, the Railroad Valley site has been used for the calibration of various sensors, such as Landsat 8 OLI, Landsat 5 thematic mapper (TM), Landsat 7 enhanced thematic mapper plus (ETM+), the moderate-resolution imaging spectroradiometer (MODIS) onboard the Terra and Aqua satellites, the advanced spaceborne thermal emission and reflection radiometer (ASTER), the multiangle imaging spectroradiometer (MISR), and the constellation of RapidEye satellites [4].
In addition to the use of LES, 'pseudo invariant calibration sites' or PICS, usually corresponding to stable desert sites exhibiting high reflectance, can be used for inter-sensor cross-calibration [5,6]. A list of CEOS-endorsed PICS is included in [7]. Of these, the Libya-4 site has been shown to be one of the most stable sites and hence it is widely used.
The objective of this paper is to assess differences in the TOA signal provided by the various instruments (i.e., PROBA-V, S2A/S2B, L8, and Deimos-1) by using observations over both the Libya-4 desert site and the RadCalNet sites. The optical sensor calibration with simulated radiances (OSCAR) desert calibration approach [5] is used over Libya-4. As part of this approach, simulated TOA bidirectional reflectance factors (BRFs) define an absolute reference against which optical sensors can be cross-calibrated. In [5], the OSCAR Libya-4 site was used for the cross-calibration of low-resolution satellite missions, such as SPOT-VGT-2, MODIS-Aqua, the medium-resolution imaging spectrometer (MERIS), and polarization and anisotropy of reflectances for atmospheric sciences coupled with observations from a lidar (PARASOL). The same method is applied here. For the RadCalNet sites, the objective is to evaluate the extent to which the modeled nadir TOA hyperspectral reflectance data, which are publicly available through the RadCalNet portal, allow for inter-sensor cross-calibration.
We start with a brief overview of the sensors considered in the study (Section 2). Section 3 describes the OSCAR Libya-4 and RadCalNet methodology and data used. The results for both approaches are given in Section 3 and are discussed further in Section 4. Finally, conclusions are drawn in Section 5.

PROBA-V
The project for on-Board autonomy-VEGETATION (PROBA-V) is a small satellite, less than a cubic meter in size, designed to monitor global vegetation. It was launched on 6 May 2013 and was designed to bridge the gap in space-borne vegetation measurements between SPOT-VGT (March 1998-May 2014) and the two Sentinel-3 satellites, the first of which has been in orbit since 16 February 2016 and the second since 25 April 2018.
The VEGETATION instrument has a field of view of 102 • , resulting in a swath width of 2295 km. This swath width ensures near-global daily coverage (90%), whereas full global coverage is achieved every two days. The optical design of PROBA-V consists of three cameras. The central camera (considered in this study) observes at 100 m nominal resolution and covers a swath of approximately 517 km, which ensures global coverage every five days. PROBA-V observes in four spectral bands: BLUE (centered at 0.463 µm), RED (0.655 µm), near infrared (NIR, 0.837 µm), and short-wave Infrared (SWIR, 1.603 µm, see Figure 1). Due to the lack of propellant onboard PROBA-V, the overpass time (10:45 at launch) decreases as a result of increasing atmospheric drag [8] and was 09:46 as of 24 September 2019. Due to the absence of on-board calibration devices, the radiometric calibration and stability monitoring of the PROBA-V instrument relies on vicarious and lunar calibration [9,10].

DEIMOS-1
The DEIMOS-1 mission is fully owned and operated by Deimos Imaging (DMI, Spain).

Landsat-8
Landsat-8 (LS8) was launched on 11 February 2013. There are two imagers onboard: The operational land imager (OLI) and the thermal infrared sensor (TIRS). The OLI has nine spectral bands (Figure 1) covering the BLUE to SWIR wavelengths with a spatial resolution of 30 m for bands 1-7 and 9. The spatial resolution for band 8 (panchromatic) is 15 m. Landsat-8/OLI has a 16-day repeat cycle, with an equatorial crossing time of 10:00 ± 15 min. The OLI has three primary radiometric calibration devices built into its hardware: (1) A shutter, (2) lamps, and (3) solar diffusers. These are used in combination with lunar observations and terrestrial targets to provide the radiometric calibration of the OLI [12].

OSCAR Libya-4
For the L1 TOA inter-comparison over the Libya-4 PICS desert site, we follow the OSCAR desert calibration approach described in [5], where it was successfully applied for the cross-calibration of low-resolution satellite missions. The OSCAR Libya-4 method relies on the comparison between cloud-free TOA reflectance as measured by the satellite sensor ("MEAS") and the modeled TOA reflectance values ("REF"). These modeled values serve as the 'absolute' calibration. The modeled TOA reflectance values are calculated with the Second Simulation of the Satellite Signal in the Solar Spectrum-Vector version (6SV) RTM using as input: Improved Rahman-Pinty-Verstraete (RPV) bi-directional reflectance distribution factor (BRDF) model [13] parameters derived from MODIS and MISR data by [5], following the approach originally described in [14]. Furthermore, European Centre for Medium-Range Weather Forecasts (ECMWF) numerical weather prediction (NWP) data (atmospheric pressure, total column ozone, and total column water vapor), a specific desert aerosol model derived for the Sahara region, a monthly variable aerosol optical thickness (AOT), based on an analysis of aerosol robotic network (AERONET [15]), and the terrain altitude are used as input. The modeled TOA reflectance values are simulated for the actual illumination and observation geometry and account for the instruments' actual spectral response. 6SV generates spectral TOA reflectance at 2.5 nm resolution. These are spectrally resampled to the S2A, S2B, LS8, Deimos-1, and PROBA-V spectral bands, considering the actual response curves of the sensors.
Among the six desert PICS sites endorsed by the CEOS/WGCV/IVOS subgroup (Committee on Earth Observation Satellites/Working Group on Calibration and Validation/Infrared and Visible Optical Sensors), the Libya-4 site has been shown one of the most stable sites [6]. However, the Libya-4 site is composed of long sand dune ridges that impact the surface BRDF as a function of the solar azimuth angle [16]. In order to account for these surface directional effects, BRDF models are traditionally used in the vicarious calibration methods [17]. In this study, the RPV BRDF model is used as described above.
The OSCAR Libya-4 calibration method has been applied to the S2A, S2B, LS8, Deimos-1, and PROBA-V TOA extractions over the Libya-4 ROI. Table 1 gives the number of cloud-free extractions used. All scenes were visually inspected from red-green-blue (RGB) quick looks to ensure that there were no undetected clouds or cloud-shadows at the site or in its vicinity (within~5 km).

RadCalNet
In order to evaluate biases of current and future satellite sensors in a harmonized manner, the RadCalNet has been established to facilitate traceability and sensor cross-comparison. RadCalNet comprises a set of four LES: Railroad Playa (USA, RVUS), Baotou (China), La Crau (France, LCFR), and Gobabeb (Namibia, GONA). The Baotou site was not included in this study, due to the small size of the representative area compared to the sensors' spatial resolution. GONA is situated in the gravel plains at the edge of the Namib Desert and was selected because of its high spatial homogeneity, low cloud coverage, relatively high surface reflectance stability, and stable atmospheric conditions [18]. LCFR is located in southeastern France, and its surface is uniformly covered by pebbles and dry grass-like vegetation. RVUS is located in Nevada (USA) in a dry lakebed or playa. Changes in the surface characteristics are generally due to rain and snow [18,19].
At the RadCalNet sites, surface and atmospheric data are automatically acquired. The RadCalNet data are publicly distributed through the RadCalNet portal [20]. Processed RadCalNet data distributed through the portal include modeled nadir TOA hyperspectral reflectance (at 10 nm spectral sampling resolution), including uncertainties at 30-min intervals from 09:00-15:00 local time (LT). The modeled nadir TOA hyperspectral reflectances are simulated using the moderate-resolution atmospheric transmission (MODTRAN-5) RTM and, after spectral resampling considering the sensors' RSRFs, these serve as a direct "REF" to compare the satellite-measured TOA signal with. It is important to highlight here that although BRDF measurements are performed at the sites, the RadCalNet portal only gives access to the nadir surface reflectance measurements and the corresponding nadir TOA reflectance simulations [21].
For each sensor, TOA reflectances were extracted or calculated for~100 ×100 m regions centered at the RadCalNet sites' geolocation, which means 3 × 3 pixels for LS8, 5 × 5 pixels for S2A/S2B and Deimos-1, and a single pixel for PROBA-V. Subsequently, the sensor-measured TOA reflectances were compared to the corresponding simulated TOA reflectances, with a maximum time difference of ±30 min. To allow for direct inter-comparison to the measured TOA reflectances, the RadCalNet-simulated TOA reflectances provided at 10 nm spectral resolution were first interpolated to 1 nm resolution and then resampled to the S2A, S2B, LS8, Deimos-1, and PROBA-V spectral bands, taking into account actual their spectral response curves. To exclude any cloud contamination to the observed TOA reflectance RGB quick looks for all sensor observations were subjected to visual inspection upon acceptance. Table 2 shows the number of cloud-free match-ups between the various TOA sensor observations and RadCalNet simulations over the three sites. It is clear that the sample sizes are relatively small and that, mainly over RVUS for S2B, Landsat-8, and Deimos-1, there are very few match-ups (< 10). For Deimos-1, only very limited cloud-free data acquisitions over LCFR and RVUS were available, while no data were available over GONA. The matchups between the sensor and simulated TOA reflectances were generally constrained by the availability of the RadCalNet simulations and the cloud occurrence over or in the respective sites' vicinity (<~5 km). According to [4,22], typically only 25% of the LS8 acquisitions over the RVUS site are successful. Reasons for unsuccessful collections are clouds over the site or water/snow on the surface. For RVUS, the presence of efflorescent salts (i.e., the migration of a salt to the surface of a porous material) after the site had dried out from a precipitation event also rendered images unsuitable for calibration purposes, due to the negative impact on the site's homogeneity.

Libya-4
For each sensor and spectral band (excluding designated cirrus detection and strong water vapor absorption bands), the mean ratios (over all observations) of the satellite-measured TOA reflectances to the 6SV TOA reflectance simulations in comparison with S2A are plotted (see Figure 2 and Table 3).   The satellite-measured and simulated TOA reflectances agree to within ±4% for all sensors and spectral bands, and even to within ±2% for most bands. The largest deviation is observed for the Deimos-1 Green band, which is approximately 3.5% lower than the 6SV reference TOA reflectance.
In view of the BELHARMONY harmonization objective, we are mainly interested in relative differences between the sensors, rather than in the absolute reference. Therefore, it is more relevant to compare the results for spectrally similar bands. This also excludes issues attributable to possible intrinsic differences in the TOA reflectance simulations. We use S2A as the reference in the inter-comparison. Please note that using S2A as reference should not be interpreted here as implying that S2A has the highest absolute radiometric accuracy. The inter-comparison results are given in Table 4. Positive values for a certain band indicate that the sensor is brighter than S2A, while the opposite applies for negative values.
For comparable spectral bands, LS8, PROBA-V, S2B, and Deimos-1 agree with S2A to within ±2%, with the exception of the Deimos-1 Green band. Deviations observed between S2A and S2B are of the same magnitude as those observed between S2A and the other missions. For most bands, S2A is slightly brighter than S2B. This is in line with the findings of [23].  Table 4. Cross-calibration results: Inter-comparison of the mean ratios for spectrally similar bands. The percentage difference of the mean ratios given in Table 3 Figure 3 shows the obtained average ratios and spanned ±1 standard deviation (σ) range for all sensors over the three RadCalNet sites. It follows that the ratios generally increase with wavelength over all sites, but that patterns for the various sensors deviate significantly over the three sites.

RadCalNet
Over LCFR, deviations between the ratios obtained for the various sensors are largest in the 0.6-0.8 µm range. The PROBA-V RED band exhibits a ratio of 0.87, while for the other sensors the ratios are~0.95-1.00. It should be noted that the LCFR site is less homogeneous than the RVUS and GONA sites, with sparse vegetation cover.
Over GONA, all sensors have ratios roughly between 0.95 and 1.00 from~0.45-0.85 µm and these gradually increase towards~1.02, with the PROBA-V SWIR band being slightly higher (1.04). Further, it can be seen that the ratio variation is much smaller than for LCFR. Similar to LCFR, large differences are observed over RVUS, as can be seen in Figure 3. The most significant differences to be noted are that Deimos-1 has ratios significantly lower than 1.0, while these ratios over LCFR are within the 0.99-1.01 range. Further, the ratios for S2B are significantly larger than those for S2A throughout the entire VNIR and SWIR spectral ranges. This is in contrast with the ratios obtained over GONA and LCFR, which exhibit similar values, S2A being slightly larger than S2B, respectively. Note that no Deimos-1 data were available for GONA.
In order to assess whether BRDF effects could play a role in the obtained results, polar plots of the viewing and illumination geometries for S2A, S2B, PROBA-V, and Deimos-1 were made for LCFR, RVUS, and GONA. These plots are shown in Figures 4-6, left columns. Note that over GONA no Deimos-1 data were available. The plots shown that only Landsat-8 observes the various sites under almost nadir viewing conditions (i.e., VZA < 1 • ), while for S2A and S2B the VZA is generally <10 • and has a low variation, similar to LS8. Figures 5 and 6 clearly show that, over LCFR and RVUS, the viewing azimuth angle (VAA) for S2A and S2B is either~135 • or~270 • . For Deimos-1, the VAA alternates between~80 • and~280 • , while for PROBA-V, the VAA range changes more continuously betweeñ 100 • and~280 • . The GONA site is always observed by S2A and S2B under similar viewing conditions.
To further analyze the impact of the changing viewing azimuth angle, we divided the ratios of observed/RadCalNet-simulated TOA reflectances into two groups based on the observation VAA: VAA < 180 • and VAA > 180 • . These subsets are presented as box-whisker plots for the RED spectral range in Figures 4-6, right columns. For all investigated sensors over LCFR and RVUS, it can clearly be seen that the median and the upper and lower quartiles of the ratio distribution are generally higher for VAA < 180 • than for VAA > 180 • . These differences contribute considerably to the rather large variation in the RadCalNet ratios for each sensor. As a result, the variation in the ratios obtained for one sensor (the intra-sensor variation) is greater than between the sensors (the inter-sensor variation). Note that no results for Deimos-1 could be established over GONA.
To further exemplify the large intra-sensor variation, Figure 7 shows the TOA reflectance measured by S2A and S2B B4 (0.665 µm) over the RVUS site. The plot includes all cloud-free observations over the site, including observations without matchups with RadCalNet. Similar to Figures 4-6, a distinction is made between observations with VAA < 180 • and with VAA > 180 • . A clear difference in the TOA reflectance depending on the VAA can be seen once again. This difference is significantly greater than the differences observed between S2A and S2B.

Discussion
The Libya-4 OSCAR method produces cross-calibration results that are more consistent and statistically robust compared to the those from the RadCalNet inter-comparison, for a number of reasons. As shown in Figures 4-7, BRDF effects strongly influence the RadCalNet results. These effects are accounted for in the OSCAR Libya-4 results through the incorporation of a BRDF model. Furthermore, no matchups with in-situ measurements are required for the OSCAR Libya-4 method. As a result, more observations were obtained in the analyses than for the individual RadCalNet sites, thereby producing results that are statistically more robust. Although overall the largest number of matchups are retrieved for the LCFR site, the results for all sensors show large error bars for this site. References [23,24] also reported on the difficulty of retrieving accurate results over the LCFR site. The large variation in the calibration results obtained over this site may be related to the site's heterogeneity, including sparse vegetation. Furthermore, due to this vegetation coverage, the convolution of the RadCalNet TOA reflectance simulations, provided at 10 nm resolution, to the (e) Figure 6. Similar to Figure 4, but for RVUS. Note that over this site subfigure (e), containing results for Deimos-1, is included.

Discussion
The Libya-4 OSCAR method produces cross-calibration results that are more consistent and statistically robust compared to the those from the RadCalNet inter-comparison, for a number of reasons. As shown in Figures 4-7, BRDF effects strongly influence the RadCalNet results. These effects are accounted for in the OSCAR Libya-4 results through the incorporation of a BRDF model. Furthermore, no matchups with in-situ measurements are required for the OSCAR Libya-4 method. As a result, more observations were obtained in the analyses than for the individual RadCalNet sites, thereby producing results that are statistically more robust. Although overall the largest number of matchups are retrieved for the LCFR site, the results for all sensors show large error bars for this site. References [23,24] also reported on the difficulty of retrieving accurate results over the LCFR site. The large variation in the calibration results obtained over this site may be related to the site's heterogeneity, including sparse vegetation. Furthermore, due to this vegetation coverage, the convolution of the RadCalNet TOA reflectance simulations, provided at 10 nm resolution, to the (e)

Discussion
The Libya-4 OSCAR method produces cross-calibration results that are more consistent and statistically robust compared to the those from the RadCalNet inter-comparison, for a number of reasons. As shown in Figures 4-7, BRDF effects strongly influence the RadCalNet results. These effects are accounted for in the OSCAR Libya-4 results through the incorporation of a BRDF model. Furthermore, no matchups with in-situ measurements are required for the OSCAR Libya-4 method. As a result, more observations were obtained in the analyses than for the individual RadCalNet sites, thereby producing results that are statistically more robust. Although overall the largest number of matchups are retrieved for the LCFR site, the results for all sensors show large error bars for this site. References [23,24] also reported on the difficulty of retrieving accurate results over the LCFR site. The large variation in the calibration results obtained over this site may be related to the site's heterogeneity, including sparse vegetation. Furthermore, due to this vegetation coverage, the convolution of the RadCalNet TOA reflectance simulations, provided at 10 nm resolution, to the sensors' spectral bands might contain a larger uncertainty for sensors with broad spectral bands that partly cover the red edge spectral region. This might explain the deviating RadCalNet results for the PROBA-V RED band observed for the LCFR site.
The calibration results obtained over the Libya-4 site using the OSCAR methodology confirm the slight bias between S2B and S2A as observed by [23,25], with S2B having slightly lower TOA reflectance values than S2A. The strong consistency between S2 and LS8 as reported in [23,25,26], with a slightly greater difference for the LS8 CA band (Band1, 0.443 µm) is also confirmed here. To our knowledge, cross-calibration results for Deimos-1 and PROBA-V with respect to S2 and LS8 have not been published before. In [27], PROBA-V was compared against MERIS, and strong consistency (within 2.5%) was indicated.

Conclusions
This work focused on the assessment of radiometric differences in the TOA signal provided by PROBA-V, Sentinel-2A and Sentinel-2B, Landsat-8, and Deimos-1. The work has been performed in the frame of the BELHARMONY project, which aims to assess and to improve multi-sensor high resolution time series consistency.
Libya-4 OSCAR desert calibration results for LS8, PROBA-V, Deimos-1, S2A, and S2B agree to within ±2% for comparable spectral bands, with the exception of the Deimos-1 Green band, for which a slightly greater difference of approximately −3.5% is observed. S2A is slightly brighter than S2B for most bands. Results confirm that all sensors investigated are well calibrated and that inter-sensor differences are minor, at least over the bright Libya-4 site.
No consistent results could be obtained over the three RadCalNet sites for the sensors investigated. In addition to the number of matchups being rather limited to make statistically robust conclusions, BRDF effects significantly influence the observed results, as many of the satellite observations are made under non-nadir conditions. Even for slightly off-nadir VZAs of~7 • , a difference in the TOA measured reflectance values over RVUS and LCFR could be observed between Eastern and Western viewing directions. Due to these BRDF effects, the variation in the obtained ratios for one sensor (the intra-sensor variation) is larger than that between the sensors (the inter-sensor variation). In order to diminish this intra-sensor variation and to draw firmer conclusions, the BRDF characterizations for all RadCalNet sites must additionally be incorporated into the RadCalNet simulations and made available in the user distribution portal.
Author Contributions: S.S. designed the structure of the paper, performed the Libya-4 analyses, and supervised the RadCalNet analyses. E.W. carried out the RadCalNet analyses and co-wrote the paper.
Funding: This research was funded by the Belgian Science Policy Office (Belspo) through the Stereo program (Belharmony project) (Contract number: NR SR/00/356).