In this study, a methodology to estimate the intertidal bathymetry from multispectral remote sensing images is presented. The technique is based on the temporal variability of the water and the intertidal zone reflectance and their correlation with the tidal height. The water spectral behavior is characterized by high absorption at the infrared (IR) band or radiation with higher wavelengths. Due to tidal cycles, pixels on the intertidal zone have higher temporal variability on the near IR spectral reflectance. The variability of IR reflectivity in time is modeled through a sigmoid function of three parameters, where the inflection parameter corresponds to the pixel elevation. The methodology was tested at the Tagus river estuary in Lisbon, Portugal, and at the Bijagós archipelago, in the West African nation of Guinea-Bissau. Multispectral images from Sentinel-2 satellites were used, after atmospheric corrections from ACOLITE processor and the derived bathymetric model validated with in situ data. The presented method does not require additional depth data for calibration, and the output can generate intertidal digital elevation models at 10 m spatial resolution, without any manual editing by the operator. The results show a standard deviation of 0.34 m at the Tagus tidal zone, with −0.50 m bias, performing better than the Stumpf ratio transform algorithm, also applied to the test areas to derive intertidal bathymetry. This methodology can be used to update intertidal elevation models with clear benefits to monitoring of intertidal dynamics, morphodynamic modeling, and cartographic update.
The skills of CODAR SeaSonde coastal high-frequency radars (HFR) located in the West Iberian Peninsula on measuring wave parameters are compared to in situ (buoy) and satellite altimeters (SA) wave observations. Significant wave heights (SWH), wave periods, and wave directions are compared over a time window of 36-months, from January 2017 to December 2019. The ability of HFR systems to capture extreme wave events is also assessed by comparing SWH measurements during the Emma storm, which hit the Iberian Peninsula in March 2018. The analysis presented in this study shows a slight overestimation of the SWH by the HFR systems. Comparisons with in situ observations revealed correlation coefficients (R) of the order of 0.69–0.87, biases below 0.60 m, root-mean-squared errors (RMSE) between 0.89 m to 1.18 m, and a slope regression between 1.01 and 1.26. Using buoy observations as reference ground truth, the comparisons with SA revealed Rs higher than 0.94, biases under 0.19 m, and RMSEs between 0.17 m and 0.42 m. Since in situ observations do not overlap all the HFR range cells (RC), and its correlation coefficients with SA have shown good agreement (R > 0.94), Sentinel-3 SA (SRAL) SWH measurements are further used for the validation of the HFR systems SWH observations. The comparison between the HFR and the SA collocated SWH observations allowed the evaluation of the ability of the radars to retrieve wave data as a function of the distance to the coast, particularly during extreme wave events. The comparison of the lower frequency (4.86 MHz) HFR coastal radars with the SA measurements showed an R of 0.94–0.99, a negative but reduced bias (−0.37), and an RMSE of 0.53 m. The higher frequency HFR systems (12–13.5 MHz) showed R between 0.53 and 0.82, and a clear overestimation of the SWH by the HFR sites.