Next Article in Journal
New Method for Automated Quantification of Vertical Spatio-Temporal Changes within Gully Cross-Sections Based on Very-High-Resolution Models
Previous Article in Journal
A Continental-Scale Assessment of Density, Size, Distribution and Historical Trends of Farm Dams Using Deep Learning Convolutional Neural Networks
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Using Sentinel-2 Images to Estimate Topography, Tidal-Stage Lags and Exposure Periods over Large Intertidal Areas

1
CESAM, Faculdade de Ciências, Universidade de Lisboa, 1749-016 Lisboa, Portugal
2
IDL, Faculdade de Ciências, Universidade de Lisboa, 1749-016 Lisboa, Portugal
*
Author to whom correspondence should be addressed.
The first two authors contributed equally to this study.
Submission received: 18 December 2020 / Revised: 13 January 2021 / Accepted: 15 January 2021 / Published: 19 January 2021
(This article belongs to the Section Environmental Remote Sensing)

Abstract

:
Intertidal areas provide key ecosystem services but are declining worldwide. Digital elevation models (DEMs) are important tools to monitor the evolution of such areas. In this study, we aim at (i) estimating the intertidal topography based on an established pixel-wise algorithm, from Sentinel-2 MultiSpectral Instrument scenes, (ii) implementing a set of procedures to improve the quality of such estimation, and (iii) estimating the exposure period of the intertidal area of the Bijagós Archipelago, Guinea-Bissau. We first propose a four-parameter logistic regression to estimate intertidal topography. Afterwards, we develop a novel method to estimate tide-stage lags in the area covered by a Sentinel-2 scene to correct for geographical bias in topographic estimation resulting from differences in water height within each image. Our method searches for the minimum differences in height estimates obtained from rising and ebbing tides separately, enabling the estimation of cotidal lines. Tidal-stage differences estimated closely matched those published by official authorities. We re-estimated pixel heights from which we produced a model of intertidal exposure period. We obtained a high correlation between predicted and in-situ measurements of exposure period. We highlight the importance of remote sensing to deliver large-scale intertidal DEM and tide-stage data, with relevance for coastal safety, ecology and biodiversity conservation.

Graphical Abstract

1. Introduction

Intertidal flats represent transitional ecosystems linking marine and terrestrial areas and are characterized by vast extensions of sediments undergoing regular cycles of inundation by seawater and exposure to air [1]. These are very productive areas, intensively used by underwater consumers such as fish during high tide and sustaining a large variety and abundance of air-breeding predators such as shorebirds during low tide. Benthic producers and invertebrates that live in the sediments are pivotal in these ecosystems, as they constitute the main source of food for these two disparate predator assemblages. Intertidal systems and associated margins also provide for key ecosystem services, acting as a nursery and habitat for economically important resources (e.g., crustaceans, molluscs and fishes, [2]), and also contributing to coastal protection, nutrient cycling and carbon storage [3,4]. Still, these ecosystems are threatened at a global scale by land reclamation [5], sediment dredging [6] and sea-level rise [7]. Indeed, a large proportion of the area covered by intertidal systems declined globally by ca. 16% between 1984 and 2016 [7] and 68% of their extent is currently subject to moderate to very high human pressure [8]. Hence, characterizing and monitoring changes in such ecologically valuable systems is a priority in order to delineate more effective conservation actions.
Over the past three decades, several techniques have been used to map intertidal areas and also to develop digital elevation models (DEM) for these areas. Some studies have involved the use of acoustic sounders installed on ships which provide reliable measurements of topography, but these techniques are very expensive and time-consuming when large areas are involved [9]. Airborne Light Detection and Ranging (LIDAR) is a very promising technique, being able to accurately generate bathymetry estimates [10] at depths of up to 70 m in clear water [11]. However, the reliability of such techniques can be affected by water turbidity, and its application is usually constrained to local and regional scales due to very high acquisition costs [12].
Due to these difficulties, intertidal topography remains poorly mapped worldwide [13,14], despite its importance to predict flood impacts and other foreseeable changes on many coastal ecosystems due to sea-level rise [15,16]. In fact, accurately generated DEMs allow for the determination of the coastal slope and for the identification of areas with a low tidal range, which are significant parameters for calculating coastal vulnerability indices, and for the assessment of the impacts of the climate change-induced sea-level rise on these areas [17,18]. DEMs have also other applications, being useful for instance, in estimating inundation/exposure periods of intertidal areas [19]. This is a meaningful exercise from an ecological perspective as one important factor affecting biota zonation in intertidal areas is the exposure period [20,21], an attribute difficult to obtain over large scales by traditional sampling or observation [19].
The “waterline method” constitutes a simple and cost-effective alternative to derive the intertidal topography over large geographic regions [22]. This multi-temporal approach estimates intertidal topography from a set of topographic contours obtained from satellite images, depicting the position of the waterline at the widest possible range of tidal stages [14]. These contours can be assigned to a specific elevation by using hydrodynamic models to calculate the water height at the time of the acquisition of each image. Synthetic Aperture Radar (SAR) or Multispectral Instruments (MSI), such as Sentinel-2 sensors, have been used to this end and may be the only practical way to extract elevation values in remote and inaccessible places [23,24,25]. Recently, the authors of [26] proposed a variation to the waterline method, consisting of a pixel-based technique using the backscatter of a time-series of SAR Sentinel-1 images to generate a DEM of intertidal areas. The technique was later generalized for multispectral images [27]. This technique explores the relationship between the near infrared (NIR) reflectance and pixel elevation, which can be described by a logistic regression, relating a set of values of the water height at the time of image acquisition with the corresponding NIR values of each pixel. In this relationship, NIR reflectance values at a given point will suddenly shift from high values, when the pixel is exposed, to low values, when the pixel is covered by water (which strongly absorbs IR radiation). The water height at this shift point is assumed to represent the height of that particular pixel [26,27].
The application of this method requires precise information concerning the height of the water at the particular time of image acquisition. Such data are often calculated from tidal models developed for a given region, based on longitudinal observations collected by a network of sea-level gauges which are used to derive the regions’ tidal constituents. However, in remote locations the number of available gauges is often very limited and consequently tidal models are unable to capture spatial differences in the stage of the tide over large extents. In such situations, one often assumes that the water height measured at one reference point is the same across the entire region, a simplification that can lead to significant errors in the estimates of elevation values when tide-stage differences are significant [27]. An alternative approach is to use tidal prediction models, such as the Oregon State University Tidal Prediction Software (OTPS) [28], which had been used for continental scale time-series tide modelling [25,29]. However, the accuracy of tide heights estimation again decreases with increasing distances from the gauge location, most particularly in estuarine and insular areas where river flows and landmass configuration can alter tidal currents and cause considerable tide-stage lags [25].
In this study, we propose a set of procedures to address the issues outlined above and further develop the quality of DEM estimation of large intertidal flats using images sensed by Sentinel-2 Multi Spectral Instrument (S2A&B, hereafter). Using a time-series of images of the Bijagós Archipelago off the coast of Guinea-Bissau, West Africa, we (i) propose a four-parameter alternative for the pixel-based logistic algorithm proposed by the authors of [27] to produce a DEM of the intertidal area, and (ii) implement a pre-calibration procedure for all bands of the images in the time-series to reduce temporal variability at the pixel level. Furthermore, we (iii) propose a novel method to estimate cotidal lines in the entire extent of the Sentinel-2 scene and (iv) produce a final DEM of the area, after correcting for differences in tidal-stage across the region. Finally, (v) we use the final DEM to produce a detailed (10 m resolution) map of exposure periods, and (vi) assess the quality of the prediction using measurements of the exposure period made in the field.

2. Materials and Methods

2.1. Study Area

This study took place in the Bijagós Archipelago, located off the coast of Guinea-Bissau (11°52′N, 15°36′W, Figure 1). The archipelago comprises 88 islands and islets, and it is classified as a Biosphere Reserve by UNESCO (2011) and as a Wetland of International Importance under the Ramsar Convention. The total area of this archipelago is ca. 10,470 km2, of which 1200 km2 are occupied by intertidal mudflats and sandbanks [30] and ca. 650 km2 are covered by mangrove forests [31]. The Bijagós Archipelago has a semi-diurnal tide regime, with strongest spring tides ranging from 0.3 m to 4.8 m according to the Portuguese Hydrographic Institute (IH—Instituto Hidrográfico) [32].

2.2. Image Selection and Pre-Processing

A total of 35 Level-1C-processed images from S2A&B were downloaded from Sentinel’s Scientific Data Hub (https://scihub.copernicus.eu/), georeferenced in WGS84/UTM 28 N time zone [33,34] (the specifications of the MultiSpectral Instrument onboard Sentinel-2 can be found at https://sentinel.esa.int/web/sentinel/missions/sentinel-2/instrument-payload/resolution-and-swath). The images were sensed in two different periods: between December 2017 and May 2018 and November 2019 and April 2020 (Table S1, Supplementary Material). In addition to presenting a cloud cover lower than 10%, images were selected to depict the most diverse situations of tidal stages in order to get the most accurate representation of the spatial progression of the water across the intertidal area.
Atmospheric correction is a crucial step in the analysis of multispectral satellite images, namely to improve the accuracy of bathymetric models [23,35,36]. Hence, we converted all Level-1C top of atmosphere reflectance products (ρtoa) into Level-2A (L2A) bottom of atmosphere surface reflectance products (ρs), using the ACOLITE software (version 20190326.0) developed by the Royal Belgian Institute of Natural Sciences (RBINS) [37]. The ACOLITE processor was especially designed for aquatic applications and it uses the “dark spectrum fitting” approach, described in [38,39,40]. For the purposes of this study, we only used the green (560 nm) and NIR (~833 nm) bands as processed by ACOLITE, both at 10 m resolution. All other analyses were carried out in R software v.4.0.0 [41], mostly using the library raster [42] and its dependencies.
Despite this atmospheric correction, scenes from the same area obtained in different dates will often show some differences in the range of band reflectance, even if representing regions that remained almost unchanged. Such differences can result from slight differences in light conditions at the corresponding acquisition dates and can be detected by simply plotting the reflectance of a given band obtained in two different periods for the same area (Figure 2). To calibrate the bands of each scene to a common dynamic range, we used one of the scenes (informally selected to have an average reflectance range in the green and NIR bands) as a reference and then we calculated band-wise major axis regressions with all other images, which were then corrected using the corresponding regression coefficients (Figure 2). Major axis regressions minimize the sum of squared distances between the points and the regression line, measured perpendicularly to the line, and were estimated using lmodel2 package [43] in R software. While fitting these regressions, we only included pixels with very low (<0.05) and high (0.2) NIR reflectance, representing seawater and land areas, respectively, and excluded intertidal areas which vary considerably among scenes, according to the tidal-stage (Figure 2).

2.3. Estimation of Water Heights in Sentinel-2 Scenes

We compiled information on the time of high- and low-water at the Bijagós Archipelago, as well as the corresponding water heights for each scene used in this study, from the Portuguese Hydrographic Institute (https://www.hidrografico.pt/). The data available refer to a single location in the archipelago, Bubaque (Longitude 15°50.0′W; Latitude 11°18.1′N). Based on these values, we estimated the water height at the exact time of image acquisition (hsat) as:
h s a t = h H i g h W a t e r ( h H i g h W a t e r h L o w W a t e r ) × ( c o s ( π × T s a t T L o w W a t e r T H i g h W a t e r T L o w W a t e r ) + 1 ) 2
where hsat corresponds to the water height predicted at the satellite sensing time (Tsat), in meters (m), hLowWater and hHighWater correspond to the low- and high-water marks, respectively, and TLowWater and THighWater to the time of low and high water in hours, respectively (water heights are presented in Table S1, Supplementary Material).

2.4. Identification of the Intertidal Area

The method developed by authors in [27] is computationally intensive, since it requires a set of operations carried out at each pixel of the image. To avoid unnecessary calculations, the process starts by identifying the zone of interest, comprising only pixels in the intertidal area. To achieve this, we first calculated the Normalized Difference Water Index (NDWI, proposed by authors in [44] for each image), to maximize the contrast between water and other land cover types. This index is calculated as a standardized ratio between the green and NIR reflectance ( ρ ), as given by equation:
N D W I = ρ G r e e n ρ N I R ρ G r e e n + ρ N I R
Pixels corresponding to permanent water bodies and to land have a low temporal variability in their NDWI values, since they represent relatively constant conditions. In contrast, pixels in the intertidal zone typically show a higher variability in NDWI, as they are subject to cycles of exposure to air and coverage by water. We followed the authors of [27] to identify areas of land, water and intertidal areas using increasing thresholds of temporal variability in NDWI. Using the NDWI of all scenes, we then calculated the temporal standard deviation in NDWI of each pixel, based on all scenes, as:
σ ( x ,   y ) = 1 M i = 1 M ( N D W I i N D W I ¯ ) 2
where N D W I ¯ is the temporal mean of the NDWI of each pixel (x, y) across all images (M).

2.5. Estimation of the Height of Intertidal Pixels

The elevation of each pixel can be estimated by the height of the water at which its NIR reflectance suffers a sudden change in magnitude. Emerged areas of the intertidal typically show high NIR reflectance, which will suddenly drop to very low values when the water line reaches that area, due to the high NIR absorption. Therefore, as proposed by [27], one expects a sigmoid relationship between the height of the water and the NIR reflectance of a given pixel, given by:
ρ i ( x , y ) = ( t o p a s y m p t o t e b o t t o m a s y m p t o t e ) 1 + e a ( h i h p i x e l ) + b o t t o m a s y m p t o t e , i = 1 , , M
where ρi is the standardized NIR reflectance of the pixel with coordinates (x, y), hpixel corresponds to the elevation (in meters) of the pixel (x, y) and hi is the water height associated with each scene (in meters, Figure 3). Besides the inflection point (hpixel), which is our main parameter of interest, the logistic function is characterized by three other parameters: top asymptote, bottom asymptote and steepness (a). These four parameters were estimated for each intertidal pixel using the nplr package [45], running in R software [41]. The algorithm uses a Newton-Raphson method to minimize a weighted sum of squared errors [45], and whenever model convergence was not reached, the pixel was marked as being outside of the intertidal zone. This procedure was run for each pixel of the intertidal area, resulting in a preliminary DEM of the archipelago.

2.6. Estimation of Differences in Tidal-Stage across the Archipelago

As mentioned above, the times of high- and low-tides available for the Bijagós Archipelago have been calculated from measurements at a single location, in the Island of Bubaque (see Figure 1). However, the archipelago extends over ca. 100 km in latitude and in longitude, and therefore considerable differences in tidal stages are expected to occur across the region. Such differences will obviously impact the estimation of the heights of each pixel, since the height of the water at a particular location at the time of acquisition of the satellite images (hi in Equation (4)) will differ from that of the reference location (in Bubaque). The difference between tidal stages will obviously increase with an increasing distance between the locations along the axis of tide-wave progression, introducing a systematic geographic bias in the estimation of intertidal DEM. To account for this problem, we estimated the tide lags for the Bijagós Archipelago by examining the differences between the predictions of pixel height obtained using satellite scenes representing rising tides only, as compared with those obtained with images acquired during ebbing tides. In fact, if the local tidal-stage is correctly determined, the pixel height estimated from images depicting ebbing tides should be fairly similar to that obtained with data representing rising tides. However, if the predicted water height at a given location occurs e.g., 1 h after the reference value used in the model, one would expect that predictions of pixel height based on rising tides will be underestimated (as the water is yet to reach the height required to cover the pixel), while predictions based on ebbing tide deliver overestimated heights (because the water will still be covering that pixel for 1 h, see Figure 4c). Likewise, when the tidal-stage at a place occurs, for example, 1 h before the time in the reference location, the opposite will occur (Figure 4a).
We explored these differences among tide stages to estimate the time-lags of each pixel. To do this, we run two separate logistic models, one based on images representing ebbing tides and another one using rising tides. We then defined a sequence of potential tide time-lags ranging from −90 min to +90 min (at 5 min intervals) in relation to the tidal time of our reference point (at Bubaque), and for each of these values we recalculated the water height of all images within each dataset (one for scenes acquired during ebbing tides and another for rising tides). We fit a set of models for each time-lag and dataset, and finally selected the time-lag that resulted in the smallest difference between the estimated pixel height of ebbing- and rising-tide models (see Figure 5). This computationally intensive procedure was repeated for 50,000 random pixels in the intertidal area, selected among those previously estimated to have an elevation within 2.47 ± 0.25 m (see above). This value corresponds to the mean water height of all satellite images, ensuring that the distribution of points around the inflection point is more balanced and therefore parameter estimation is more accurate.
To estimate the time differences in tidal-stage for all intertidal pixels of our study region, we fit a generalized additive model (GAM) using the time differences in tidal-stage estimated at the 50,000 random pixels as dependent variable and the longitude and latitude as predictors, using mgcv package [46]. GAM fit involved the use of penalized thin-plate regression splines, with regression parameters (i.e., amount of smoothing) selected by generalized cross-validation [46]. We choose this type of model because tide-stage delays are expected to show a very strong spatial (latitudinal and longitudinal) correlation, hence a very smooth spatial variation which will be well captured by the use of thin-plate splines.

2.7. Production of the Final Bathymetric Map

Based on the previous procedure, we produced a map representing the spatial differences in tidal-stage across the archipelago. Using this information, we re-estimated the heights of all intertidal pixels using the logistic function as explained above, but now calculating the water heights of all images according to the time-lag in the tide stage, in relation to the reference point (Bubaque).

2.8. Estimation and Validation of the Intertidal Exposure Period

We used the corrected pixel elevation to estimate the exposure period of all pixels in the intertidal zone, assuming a simple sinusoidal function and a tide cycle lasting 12.40 h. This value corresponds to the mean cycle duration recorded in our Sentinel-2 scenes dataset (mean = 12.40 ± 0.13 h, n = 35 days). We used this value in a simple sinusoidal function to estimate the exposure period (in hours) of each pixel:
E x p o s u r e ( h ) = 12.40 × ( 1 a c o s ( 2 × ( h p i x e l h L o w W a t e r h H i g h W a t e r h L o w W a t e r ) 1 ) π )
where hpixel is the height of the pixel (m), and hHighWater and hLowWater are the high- and low-water marks (m), respectively. We used this formula to produce a map representing the exposure periods of intertidal sediments for an average tide, using a mean high- and low-water mark of 1.05 m (±0.55 m SD) and 3.90 m (±0.54 m SD), respectively, calculated as averages of the values obtained in the Sentinel-2 scenes dataset (n = 35).
To validate the map of the exposure period, we measured the exposure/inundation periods at 66 sites using time-lapse cameras in different islands of the archipelago. Between February and March 2020, we set several cameras in the intertidal area, fixed them to a stick and facing at its base to record the exact time the water reached each site during ebbing and rising tides (Figure S1 in Supplementary Material). The location of each camera was georeferenced with a portable GPS (ca. 4 m accuracy), and recorded footage at 1-min intervals which were subsequently examined to determine the exposure period of each site. These field measurements were compared with exposure values estimated by the models at the same pixels, as determined by their GPS positions. To calculate the predicted height of each pixel, we used the tide parameters (hHighWater and hLowWater) for the corresponding camera sampling dates. A schematic representation of all procedures is presented in Figure 6.

3. Results

3.1. Identification of the Intertidal Area

The inter-calibration of the images was the first step of all the analysis described below and had a visible effect on the dynamic range of the NIR reflectance in our image set (see Figure 2), most especially by reducing the variability at low reflectance levels which in turn affected the pixel height estimation (Figure S2 in Supplementary Material). The distribution of the standard deviations in NDWI calculated for all pixels showed three clear peaks (Figure 7), with the lowest values associated with land (highest stability in NDWI values), and the highest temporal variability associated with intertidal areas. The highest peak, corresponding to areas of water, showed intermediate variability in NDWI. We defined a threshold of SD(NDWI) = 0.2, above which pixels were considered to represent intertidal areas (see Figure S3 in Supplementary Material). This conservative value was defined as to minimize the omission error for pixels in the intertidal area, and any pixel misclassified as water will be excluded in subsequent analytical steps.

3.2. DEM of the Bijagós Archipelago

The automated application of the logistic function (Equation (4)) to each pixel provided a map of the elevation of the intertidal flats of the archipelago with a high spatial detail, capturing the slopes associated with small channels in the intertidal area (Figure 8). The estimated elevation values ranged from 1.04 to 4.69 m, limited by the lowest and highest tides recorded in our scene dataset. According to the Portuguese Hydrographic Institute, the hydrographic reference at Bubaque is 2.54 m below mean sea level.

3.3. Estimation of Spatial Differences in Tidal-Stage at the Archipelago Scale

The representation of differences in tidal-stage across the archipelago is shown in Figure 9. The orientation of the cotidal lines indicates that the tide progresses along a SW-NE axis, and involves a difference of ca. 1 h in the tidal-stage between the most distant regions in that axis. The stage differences provided by the Hydrographic Institute at four locations were relatively close to those estimated by our model (average absolute difference = 6.6 min), with a maximum discrepancy in Abu (Figure 1), where our prediction underestimates the reported tidal stage by 15 min (see Table S2, in Supplementary Material). A final DEM of the archipelago was calculated after correcting the water height of each image and accounting for the spatial differences in the tidal stage (Figure S4, in Supplementary Material).

3.4. Mapping and Validating the Intertidal Exposure Period

Figure 10 represents the exposure period estimated from the DEM, for an average tide in the archipelago (low water mark = 1.05 m, high water mark = 3.9 m, amplitude = 2.85 m). There was a strong agreement between the predicted exposure periods and those measured in the field using time-lapse cameras (r2 = 0.94; n = 66, p < 0.001), with low exposure values presenting lower prediction accuracy (Figure 11).

4. Discussion

In this study we present, for the first time, a set of automated remote sensing methods that allow the estimation of tidal-stage lags over large intertidal areas, the production of cotidal lines, and the predictions of the exposure time of intertidal sediments, which were validated by in-situ measurements of the position of the water line. These procedures can significantly increase the quality of intertidal DEM by removing a geographical bias in water level estimation. We also refined and expanded the application of logistic regressions in a time-series of remote sensed images to estimate the topography of intertidal areas, as proposed by [26,27]. These refinements also included a prior inter-calibration of NIR reflectance among all images in the time series, and a more flexible formulation of the logistic function (now with four parameters).
The quality of the topographic predictions of the intertidal areas using our approach is strongly contingent on the availability of time-series images for a given site. The lowest and highest water marks recorded in the series (generally during spring tides) determine the range of height predictions that the models can provide. On the other hand, the availability of images portraying intermediate water heights is critical to reduce the uncertainty associated with the inflection points of the logistic function, which corresponds to the estimate of the height of the pixel. In this study, we used data from 35 Sentinel-2 A&B scenes, almost twice the number previously used to derive the topography in this area [27]. In doing so, we were able to expand the range of tides from 1.39–4.16 m [27] to 1.04–4.69 m, and also to gather more data from intermediate tidal- stages (7 out of 35 images with water heights between 2 and 3 m). We could not find images with tidal heights lower than 1.04 m, thus preventing the identification of the lowest lying intertidal areas, which only occur in the strongest spring tides, exposing sediments at heights of 0.5 to 1.0 m. However, such strong spring tides are rather uncommon, which reduces the limitation of our approach.
The atmospheric correction of multispectral images is an important step prior to image analysis, as aerosols and gases have a strong influence in the nature of signal, particularly in water bodies [47], where they can constitute up to 90% of the signal [48,49]. This is particularly relevant for the method presented here, which relies to a very large extent on the large and rapid change in pixel reflectance in relation to the presence of water. Our study shows that, despite the atmospheric correction implemented by ACOLITE, some discrepancies persisted in the dynamic range of NIR reflectance of images obtained in different dates (see Figure 2). The application of a major axis regression visibly reduced this unwanted source of variability in the NIR values (see Figure 3), ultimately improving the estimation of the inflection point of the logistic function, from which the height of each pixel is calculated. This linear correction is very easy to implement and it is likely to be of use in several other remote sensing techniques involving the use of time-series of satellite images, such as crop monitoring [50] or annual monitoring of tree vitality [51].
Modelling the water level considering the tide-stage lags over large areas generally improves the accuracy of DEMs [52,53]. Some efforts have been developed to address this issue, but they usually involve continental-scale models (e.g., [25]), where the gravitational effects of the moon are much more significant than, e.g., the smaller scale physiography of the coast of land. Furthermore, at a larger scale there is higher availability from tide gauges to help modelling tide-stage differences. In this study we propose a new method to assess spatial differences in tidal-stage at a smaller scale based on the expected differences in pixel-height estimation obtained from using ebbing tides and rising tides separately. Our approach is still based on the NIR reflectivity of each used pixel and requires a single site as a reference height. This new development enables the estimation of cotidal lines from pixel height estimations made in intertidal areas, and such methodology can be applied in other regions whenever tidal areas are distributed more or less evenly across the region of interest. The Bijagós Archipelago represents an almost ideal setting for such application, since there are extensive intertidal areas associated with each island providing a good spatial coverage for tide lag estimation. Cotidal maps are by themselves useful products but are also a significant advance as a means of reducing the geographical bias in techniques that rely on exact water heights registration across the region of interest. This is particularly relevant in remote areas where tide gauges are scarce and, therefore, information on tide stage lags is totally absent.
The estimated cotidal lines for the Bijagós Archipelago indicate a tide progression roughly from SW to NE, being more compressed in areas with higher landmass concentration, and shallower waters both causing a stronger gradient in the tidal wave. The cotidal line corresponding to the absence of tide-stage lag lies very close to the eastern tip of Bubaque (see Figure 8), thus correctly matching the location of the reference tide measurements. This observation, together with the low mean absolute error (6.6 min, see Supplementary Material Table S1), is a good indication of the adequacy of the method, despite the somewhat larger error observed at the northern part of the areas (15 min in relation to the value reported by the Hydrographic Institute). These larger differences might be due to the fact that the reference data provided by IH is probably outdated, as they are based on tidal observations carried out in the 1960s. However, they may also steam from the inability of the method to distinguish between the small discrepancies in tide lags associated with rising and ebbing tides, with spring and neap tides, with shallower and deeper areas, therefore representing average conditions for the region.
To assess the quality of our final DEM, we validated our predictions of exposure periods (produced based on the corrected DEM) with a large set of direct measurements of the waterline position. There was a very good agreement between the observed exposure periods and those predicted by our model. As expected, pixels with lower elevation showed a higher-than-expected exposure period, producing a systematic bias, particularly at heights lower than 2.2 m. To some extent, this can occur because the intertidal sediments are predominantly muddy [27]. In fact, in these areas the intertidal sediments are exposed for only a short period, and therefore often retain a thin film of surface water or form small water pools that never dry completely. This circumstance ultimately affects the NIR reflectance (because NIR radiation is absorbed by this film of water) and results in a slight overestimation of the corresponding bathymetric values. Nonetheless, the final exposure product has a very strong correlation with the observed values showing that, with some subsequent calibration work, the remote sensing approach we present can deliver large scale and very detailed DEMs and average exposure data, with relevance for many areas, from coastal safety to ecology and conservation of these critically important habitats, which are globally threatened due to climate changes [7].

5. Conclusions

The method proposed in this study allows for the estimation of differences in water elevation due to the spatial variation in the tidal phase, using Sentinel-2 scenes. We were able to detect differences of over 1 h in the tide-phase within the Bijagós Archipelago, an area encompassed by a single Sentinel-2 scene. Although this study was carried out with Sentinel-2 scenes, our approach can probably be applied to scenes acquired by other sensors with similar characteristics, including the ones onboard Landsat 8. Landsat 8 scenes have an extent of 170 km by 185 km, and therefore spatial differences in water height are likely to be even more significant. The methodological refinements proposed in the study, and the possibility to correct for differences in tide-phases, allow for important improvements in the accuracy of intertidal DEMs. This is of significant relevance, since these coastal habitats are likely to undergo important topographical changes in the near future due to sediment mobilization caused by sea-level rise and by increased frequency of extreme weather events associated with the ongoing climatic changes. Accurate intertidal maps and DEMs are also an important tool to resource (including wildlife) managers, since the topography influences several important ecosystem processes (e.g., the settlement of specific life-stages of commercially relevant marine organisms), which in turn affects the spatial variation in the structure of communities, from benthic organisms to higher vertebrates.

Supplementary Materials

The following are available online at https://0-www-mdpi-com.brum.beds.ac.uk/2072-4292/13/2/320/s1, Table S1: Sentinel-2 scenes used in this study, with indication of the estimated height of the water at Bubaque at the time of image acquisition and of the corresponding tidal stage; Table S2: Estimated and measured tide delay (minutes). Measured delays were provided by Portuguese Hydrographic Institute at four locations: Abú, Bruce, Eguba and Uno, see Figure 1. The mean absolute difference is 6.6 min (n = 5); Figure S1: Time-lapse cameras pointing down and set to detect the exposure period of intertidal sampling points; Figure S2: Example of the application of the logistic function to reflectance of a sample pixel, using images that were (red dots) and were not inter-calibrated (black dots), and its effects on the estimation of the height of the pixel. Note the reduction in variability, particularly on lower values of NIR reflectance; Figure S3: Temporal variability of the NDWI of the 35 Sentinel-2 scenes in a small part of the study area (a), showing the clear distinction between land (darker regions, with low variability), water (medium variability) and intertidal areas (clear colors, corresponding to high variability). Panel (b) shows the areas with values of variability higher than 0.2 (see methods); Figure S4: DEM of the Bijagós archipelago obtained by recalculation of height of water at each pixel in all images, according to the tide-stage differences estimated for each pixel.

Author Contributions

Conceptualization: J.P.G., J.B., M.H., J.C. and T.C.; methodology: J.P.G., J.B., J.C.; validation: M.H., J.B.; formal analysis: J.P.G., J.B.; writing—original draft preparation: J.P.G., J.B.; writing—review & editing: M.H., J.C., T.C.; funding acquisition: J.P.G., T.C.; project administration: J.P.G., T.C. All authors have read and agreed to the submitted version of the manuscript.

Funding

This study was supported by MAVA Foundation, through the project “Waders of the Bijagós: Securing the ecological integrity of the Bijagós archipelago as a key site for waders along the East Atlantic Flyway”. We also thank Fundação para a Ciência e Tecnologia (FCT/MCTES) for providing financial support through project “MigraWebs: Migrants as a seasonal ecological force shaping communities and ecosystem functions in temperate and tropical coastal wetlands” (PTDC/BIA-ECO/28205/2017) and through the financial support to CESAM (UIDP/50017/2020 and UIDB/50017/2020), from national funds. M.H. and T.C. were supported by FCT by a doctoral grant (SFRH/BD/131148/2017) and a contract IF/00694/2015, respectively.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Acknowledgments

We thank the logistical support of Instituto da Biodiversidade e Áreas Protegidas Alfredo Simões da Silva in Guinea-Bissau, and particularly Justino Biai, Aissa Regalla, Samuel Pontes, Quintino Tchatchalam and Saído Gomes. We also thank the NGO Tiniguena-Esta Terra É Nossa, for facilitating our work in Formosa Island and for logistical support and Seco Leão and several other volunteers for help during fieldwork. We also thank the suggestion made by the referees, which helped to improve the ms.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript; or in the decision to publish the results.

References

  1. Rafaelli, D.; Hawkins, S.J. The shore environment: Major gradients. In Intertidal Ecology; Springer Science & Buisness Media: Berlin/Heidelberg, Germany, 2012. [Google Scholar]
  2. Chen, Y.; Dong, J.; Xiao, X.; Zhang, M.; Tian, B.; Zhou, Y.; Li, B.; Ma, Z. Land claim and loss of tidal flats in the Yangtze Estuary. Sci. Rep. 2016, 6, 1–10. [Google Scholar] [CrossRef]
  3. Chmura, G.L.; Anisfeld, S.C.; Cahoon, D.R.; Lynch, J.C. Global carbon sequestration in tidal, saline wetland soils. Glob. Biogeochem. Cycles 2003, 17. [Google Scholar] [CrossRef]
  4. Billerbeck, M.; Werner, U.; Bosselmann, K.; Walpersdorf, E.; Huettel, M. Nutrient release from an exposed intertidal sandflat. Mar. Ecol. Prog. Ser. 2006, 316, 35–51. [Google Scholar] [CrossRef]
  5. Day, J.W.; Crump, B.C.; Kemp, W.M.; Yáñez-Arancibia, A. Introduction to Estuarine Ecology. In Estuarine Ecology; Springer Science & Buisness Media: Berlin/Heidelberg, Germany, 2012; pp. 1–18. ISBN 9780471755678. [Google Scholar]
  6. Piersma, T.; Koolhaas, A.; Dekinga, A.; Beukema, J.J.; Dekker, R.; Essink, K. Long-term indirect effects of mechanical cockle-dredging on intertidal bivalve stocks in the Wadden Sea. J. Appl. Ecol. 2001, 38, 976–990. [Google Scholar] [CrossRef] [Green Version]
  7. Murray, N.J.; Phinn, S.R.; DeWitt, M.; Ferrari, R.; Johnston, R.; Lyons, M.B.; Clinton, N.; Thau, D.; Fuller, R.A. The global distribution and trajectory of tidal flats. Nature 2019, 565, 222–225. [Google Scholar] [CrossRef] [PubMed]
  8. Hill, N.K.; Woodworth, B.K.; Phinn, S.R.; Murray, N.J.; Fuller, R.A. Global protected-area coverage and human pressure on tidal flats. Conserv. Biol. 2020, 1–19. [Google Scholar] [CrossRef]
  9. Chénier, R.; Faucher, M.A.; Ahola, R. Satellite-derived bathymetry for improving Canadian Hydrographic Service charts. ISPRS Int. J. Geo Inf. 2018, 7, 306. [Google Scholar] [CrossRef] [Green Version]
  10. Bishop-Taylor, R.; Sagar, S.; Lymburner, L.; Beaman, R.J. Between the tides: Modelling the elevation of Australia’s exposed intertidal zone at continental scale. Estuar. Coast. Shelf Sci. 2019, 223, 115–128. [Google Scholar] [CrossRef]
  11. Su, H.; Liu, H.; Heyman, W. Automated derivation of bathymetric information from multi-spectral satellite imagery using a non-linear inversion model. Mar. Geod. 2008, 31, 281–298. [Google Scholar] [CrossRef]
  12. Klemas, V. Beach profiling and LIDAR bathymetry: An overview with case studies. J. Coast. Res. 2011, 27, 1019–1028. [Google Scholar] [CrossRef]
  13. Eakins, B.W.; Grothe, P.R. Challenges in building coastal digital elevation models. J. Coast. Res. 2014, 30, 942–953. [Google Scholar] [CrossRef]
  14. Tseng, K.H.; Kuo, C.Y.; Lin, T.H.; Huang, Z.C.; Lin, Y.C.; Liao, W.H.; Chen, C.F. Reconstruction of time-varying tidal flat topography using optical remote sensing imageries. ISPRS J. Photogramm. Remote Sens. 2017, 131, 92–103. [Google Scholar] [CrossRef]
  15. Bird, C.O.; Bell, P.S.; Plater, A.J. Application of marine radar to monitoring seasonal and event-based changes in intertidal morphology. Geomorphology 2017, 285, 1–15. [Google Scholar] [CrossRef] [Green Version]
  16. Silva, A.N.; Taborda, R.; Catalão, J.; Freire, P. DTM extraction using video-monitoring techniques: Application to a fetch limited beach. J. Coast. Res. 2009, 203–207. [Google Scholar]
  17. Boruff, B.J.; Emrich, C.; Cutter, S.L. Erosion hazard vulnerability of US coastal counties. J. Coast. Res. 2005, 21, 932–942. [Google Scholar] [CrossRef] [Green Version]
  18. Karymbalis, E.; Chalkias, C.; Ferentinou, M.; Chalkias, G.; Magklara, M. Assessment of the Sensitivity of Salamina (Saronic Gulf) and Elafonissos (Lakonic Gulf) islands to Sea-level Rise. J. Coast. Res. 2014, 70, 378–384. [Google Scholar] [CrossRef]
  19. Hickey, R. Tidal inundation modeling within GIS. J. Coast. Conserv. 2019, 23, 599–606. [Google Scholar] [CrossRef]
  20. Granadeiro, J.P.; Andrade, J.; Palmeirim, J.M. Modelling the distribution of shorebirds in estuarine areas using generalised additive models. J. Sea Res. 2004, 52, 227–240. [Google Scholar] [CrossRef]
  21. Peterson, C.H. Intertidal zonation of marine invertebrates in sand and mud. Am. Sci. 1991, 79, 236–249. [Google Scholar]
  22. Mason, D.C.; Davenport, L.J. Accurate and efficient determination of the shoreline in ERS-1 SAR images. IEEE Trans. Geosci. Remote Sens. 1996, 34, 1243–1253. [Google Scholar] [CrossRef]
  23. Casal, G.; Monteys, X.; Hedley, J.; Harris, P.; Cahalane, C.; McCarthy, T. Assessment of empirical algorithms for bathymetry extraction using Sentinel-2 data. Int. J. Remote Sens. 2019, 40, 2855–2879. [Google Scholar] [CrossRef]
  24. Pahlevan, N.; Sarkar, S.; Franz, B.A.; Balasubramanian, S.V.; He, J. Sentinel-2 MultiSpectral Instrument (MSI) data processing for aquatic science applications: Demonstrations and validations. Remote Sens. Environ. 2017, 201, 47–56. [Google Scholar] [CrossRef]
  25. Sagar, S.; Roberts, D.; Bala, B.; Lymburner, L. Extracting the intertidal extent and topography of the Australian coastline from a 28 year time series of Landsat observations. Remote Sens. Environ. 2017, 195, 153–169. [Google Scholar] [CrossRef]
  26. Catalao, J.; Nico, G. Multitemporal Backscattering Logistic Analysis for Intertidal Bathymetry. IEEE Trans. Geosci. Remote Sens. 2017, 55, 1066–1073. [Google Scholar] [CrossRef]
  27. Bué, I.; Catalão, J.; Semedo, Á. Intertidal bathymetry extraction with multispectral images: A logistic regression approach. Remote Sens. 2020, 12, 1311. [Google Scholar] [CrossRef] [Green Version]
  28. Egbert, G.D.; Erofeeva, S.Y. Efficient inverse modeling of barotropic ocean tides. J. Atmos. Ocean. Technol. 2002, 19, 183–204. [Google Scholar] [CrossRef] [Green Version]
  29. Murray, N.J.; Phinn, S.R.; Clemens, R.S.; Roelfsema, C.M.; Fuller, R.A. Continental scale mapping of tidal flats across east Asia using the landsat archive. Remote Sens. 2012, 4, 3417–3426. [Google Scholar] [CrossRef] [Green Version]
  30. Cuq, F.; Capedron, P.; Giraudet, J.; Giraudet, E.; Gourmelon, F.; Pennober, G.; Silva, A. Un système d’information géographique pour l’aide à la gestion intégrée de l’archipel des Bijagos (Guinée-Bissau). In Notice de la Carte, Constitution et Exploitation du SIG; Cuq, F., Ed.; Brest Géosystèmes: Brest, France, 2001. [Google Scholar]
  31. Temudo, M.P.; Cabral, A.I. The Social Dynamics of Mangrove Forests in Guinea-Bissau, West Africa. Hum. Ecol. 2017, 45, 307–320. [Google Scholar] [CrossRef]
  32. Instituto Hidrográfico. Tabelas de Marés; Vol II Países Africanos de Língua Oficial Portuguesa e Macau; Instituto Hidrográfico: Lisbon, Portugal, 2020. [Google Scholar]
  33. European Space Agency. Sentinel-2 User Handbook; ESA Standard Document; ESA: Paris, France, 2015; Available online: https://sentinel.esa.int/documents/247904/685211/Sentinel-2_User_Handbook (accessed on 29 October 2020).
  34. European Space Agency. Sentinel-2 MSI Technical Guide. 2017. Available online: https://earth.esa.int/web/sentinel/technical-guides/sentinel-2-msi (accessed on 29 October 2020).
  35. Hedley, J.D.; Roelfsema, C.; Brando, V.; Giardino, C.; Kutser, T.; Phinn, S.; Mumby, P.J.; Barrilero, O.; Laporte, J.; Koetz, B. Coral reef applications of Sentinel-2: Coverage, characteristics, bathymetry and benthic mapping with comparison to Landsat 8. Remote Sens. Environ. 2018, 216, 598–614. [Google Scholar] [CrossRef]
  36. Vanhellemont, Q.; Ruddick, K. Atmospheric correction of metre-scale optical satellite data for inland and coastal water applications. Remote Sens. Environ. 2018, 216, 586–597. [Google Scholar] [CrossRef]
  37. Royal Belgium Institute of Natural Sciences. ACOLITE Python User Manual. 2019. Available online: https://odnature.naturalsciences.be/remsem/software-and-data/acolite (accessed on 29 October 2020).
  38. Vanhellemont, Q.; Ruddick, K. Turbid wakes associated with offshore wind turbines observed with Landsat 8. Remote Sens. Environ. 2014, 145, 105–115. [Google Scholar] [CrossRef] [Green Version]
  39. Vanhellemont, Q.; Ruddick, K. Advantages of high quality SWIR bands for ocean colour processing: Examples from Landsat-8. Remote Sens. Environ. 2015, 161, 89–106. [Google Scholar] [CrossRef] [Green Version]
  40. Vanhellemont, Q.; Ruddick, K. Acolite for Sentinel-2: Aquatic applications of MSI imagery. ESA Spec. Publ. 2016, SP-740, 9–13. [Google Scholar]
  41. R Core Team. R: A Language and Environment for Statistical Computing; R Foundation for Statistical Computing: Vienna, Austria, 2020; Available online: https://www.R-project.org/ (accessed on 29 October 2020).
  42. Hijmans, R.J. Raster: Geographic Data Analysis and Modeling. R Package Version 3.1-5. 2020. Available online: https://CRAN.R-project.org/package=raster (accessed on 20 June 2020).
  43. Legendre, P. lmodel2: Model II Regression. R Package Version 1.7-3. 2018. Available online: https://CRAN.R-project.org/package=lmodel2 (accessed on 20 June 2020).
  44. McFeeters, S.K. The use of the Normalized Difference Water Index (NDWI) in the delineation of open water features. Int. J. Remote Sens. 1996, 1425–1432. [Google Scholar] [CrossRef]
  45. Commo, F.; Bot, B.M. nplr: N-Parameter Logistic Regression. R Package Version 0.1-7. 2016. Available online: https://CRAN.R-project.org/package=lmodel2 (accessed on 20 June 2020).
  46. Wood, S.N. R, Generalized Additive Models: An Introduction with R, 2nd ed.; Chapman and Hall/CRC: Boca Raton, FL, USA, 2017. [Google Scholar]
  47. IOCCG (International Ocean-Colour Coordinating Group). Atmospheric Correction for Remotely-Sensed Ocean Colour Products; Wang, M., Ed.; Technical Report; International Ocean Colour Coordinating Group (IOCCG): Dartmouth, NS, Canada, 2010. [Google Scholar]
  48. Gordon, H.R.; Morel, A.Y. Remote assessment of Ocean Color for Interpretation of Satellite Visible Imagery. In A Review Lecture Notes ib Coastal and Estuarine Studies; Springer: New York, NY, USA, 1983; p. 114. [Google Scholar]
  49. Pereira-Sandoval, M.; Ruescas, A.; Urrego, P.; Ruiz-Verdú, A.; Delegido, J.; Tenjo, C.; Soria-Perpinyà, X.; Vicente, E.; Soria, J.; Moreno, J. Evaluation of atmospheric correction algorithms over spanish inland waters for sentinel-2 multi spectral imagery data. Remote Sens. 2019, 11, 1469. [Google Scholar] [CrossRef] [Green Version]
  50. Defourny, P.; Bontemps, S.; Bellemans, N.; Cara, C.; Dedieu, G.; Guzzonato, E.; Hagolle, O.; Inglada, J.; Nicola, L.; Rabaute, T.; et al. Near real-time agriculture monitoring at national scale at parcel resolution: Performance assessment of the Sen2-Agri automated system in various cropping systems around the world. Remote Sens. Environ. 2019, 221, 551–568. [Google Scholar] [CrossRef]
  51. Navarro, A.; Catalao, J.; Calvao, J. Assessing the use of Sentinel-2 time series data for monitoring Cork Oak decline in Portugal. Remote Sens. 2019, 11, 2515. [Google Scholar] [CrossRef] [Green Version]
  52. Bell, P.S.; Bird, C.O.; Plater, A.J. A temporal waterline approach to mapping intertidal areas using X-band marine radar. Coast. Eng. 2016, 107, 84–101. [Google Scholar] [CrossRef] [Green Version]
  53. Liu, Y.; Li, M.; Zhou, M.; Yang, K.; Mao, L. Quantitative analysis of the waterline method for topographical mapping of tidal flats: A case study in the dongsha sandbank, China. Remote Sens. 2013, 5, 6138–6158. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Study area in the Bijagós Archipelago, Guinea-Bissau. RGB-composite after atmospheric correction of Sentinel-2B imagery on 4 February 2020. Tide gauge of Bubaque is represented as a red triangle and the locations of tide measurements provided by Instituto Hidrográfico (IH) are represented as pink dots (Abú, Eguba, Uno and Bruce).
Figure 1. Study area in the Bijagós Archipelago, Guinea-Bissau. RGB-composite after atmospheric correction of Sentinel-2B imagery on 4 February 2020. Tide gauge of Bubaque is represented as a red triangle and the locations of tide measurements provided by Instituto Hidrográfico (IH) are represented as pink dots (Abú, Eguba, Uno and Bruce).
Remotesensing 13 00320 g001
Figure 2. Representation of near infrared (NIR) reflectance obtained in a reference image (a) and in a second image obtained in a different date before (b) and after calibration ((c), see Methods). Bottom panels show the relationship between NIR of both images (dashed line x = y) before (d) and after (e) the calibration using the coefficients of a major axis regression. All intertidal pixels (NIR between 0.05 and 0.2) were excluded from the analysis since their value differed strongly among scenes and since they depend on the extent of the water.
Figure 2. Representation of near infrared (NIR) reflectance obtained in a reference image (a) and in a second image obtained in a different date before (b) and after calibration ((c), see Methods). Bottom panels show the relationship between NIR of both images (dashed line x = y) before (d) and after (e) the calibration using the coefficients of a major axis regression. All intertidal pixels (NIR between 0.05 and 0.2) were excluded from the analysis since their value differed strongly among scenes and since they depend on the extent of the water.
Remotesensing 13 00320 g002
Figure 3. Relationship between the standardized NIR reflectance and the height of the water estimated for one pixel in the 35 images (each one represented as a point). The height of this particular pixel is estimated at 2.72 m by the four-parameter logistic function.
Figure 3. Relationship between the standardized NIR reflectance and the height of the water estimated for one pixel in the 35 images (each one represented as a point). The height of this particular pixel is estimated at 2.72 m by the four-parameter logistic function.
Remotesensing 13 00320 g003
Figure 4. Illustration of the expected height of the water (represented as black) in the intertidal sediments (in grey) during rising and ebbing tides, assuming different local tide stage delays: (a) −1 h, (b) 0 h and (c) +1 h in relation to the time recorded at the reference position.
Figure 4. Illustration of the expected height of the water (represented as black) in the intertidal sediments (in grey) during rising and ebbing tides, assuming different local tide stage delays: (a) −1 h, (b) 0 h and (c) +1 h in relation to the time recorded at the reference position.
Remotesensing 13 00320 g004
Figure 5. Example of a set of logistic models fitted separately to data representing images with rising (solid lines and black triangles) and ebbing tides (dashed lines and grey triangles), assuming different tide delays (from −60 min to +15 min) in relation to the reference point (in the island of Bubaque). In this example, the phase-delay of the pixel was estimated at −30 min (when both models show the closest predictions (b)).
Figure 5. Example of a set of logistic models fitted separately to data representing images with rising (solid lines and black triangles) and ebbing tides (dashed lines and grey triangles), assuming different tide delays (from −60 min to +15 min) in relation to the reference point (in the island of Bubaque). In this example, the phase-delay of the pixel was estimated at −30 min (when both models show the closest predictions (b)).
Remotesensing 13 00320 g005
Figure 6. Flow diagram of the sequence of procedures followed in this study.
Figure 6. Flow diagram of the sequence of procedures followed in this study.
Remotesensing 13 00320 g006
Figure 7. Histogram of standard deviation of Normalized Difference Water Index (NDWI) layers calculated for each pixel, with peaks representing land (lowest variability), seawater and intertidal areas (highest variability). Pixels with variability larger than 0.2 were considered to be in the intertidal area and were therefore selected for bathymetric estimation.
Figure 7. Histogram of standard deviation of Normalized Difference Water Index (NDWI) layers calculated for each pixel, with peaks representing land (lowest variability), seawater and intertidal areas (highest variability). Pixels with variability larger than 0.2 were considered to be in the intertidal area and were therefore selected for bathymetric estimation.
Remotesensing 13 00320 g007
Figure 8. Detail of a Sentinel-2 RGB composite (a) and of the corresponding elevation profile (b), obtained using the 4-parameter logistic model.
Figure 8. Detail of a Sentinel-2 RGB composite (a) and of the corresponding elevation profile (b), obtained using the 4-parameter logistic model.
Remotesensing 13 00320 g008
Figure 9. Cotidal lines, showing the time-lag (in min) from the reference point in Bubaque, estimated using data from 50,000 random intertidal pixels as a function of longitude and latitude, using a generalized additive model (GAM) with thin-plate regression splines.
Figure 9. Cotidal lines, showing the time-lag (in min) from the reference point in Bubaque, estimated using data from 50,000 random intertidal pixels as a function of longitude and latitude, using a generalized additive model (GAM) with thin-plate regression splines.
Remotesensing 13 00320 g009
Figure 10. Map of the exposure period as predicted by our model for the Archipelago (a), and details of the islands of Urok (b) and Bubaque (c).
Figure 10. Map of the exposure period as predicted by our model for the Archipelago (a), and details of the islands of Urok (b) and Bubaque (c).
Remotesensing 13 00320 g010
Figure 11. Location of in situ measurements of exposure period using time-lapse cameras (red dots) in Urok (a) and Bubaque (b). Comparison between values of exposure period measured in the field and the values predicted by the DEM (c).
Figure 11. Location of in situ measurements of exposure period using time-lapse cameras (red dots) in Urok (a) and Bubaque (b). Comparison between values of exposure period measured in the field and the values predicted by the DEM (c).
Remotesensing 13 00320 g011
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Granadeiro, J.P.; Belo, J.; Henriques, M.; Catalão, J.; Catry, T. Using Sentinel-2 Images to Estimate Topography, Tidal-Stage Lags and Exposure Periods over Large Intertidal Areas. Remote Sens. 2021, 13, 320. https://0-doi-org.brum.beds.ac.uk/10.3390/rs13020320

AMA Style

Granadeiro JP, Belo J, Henriques M, Catalão J, Catry T. Using Sentinel-2 Images to Estimate Topography, Tidal-Stage Lags and Exposure Periods over Large Intertidal Areas. Remote Sensing. 2021; 13(2):320. https://0-doi-org.brum.beds.ac.uk/10.3390/rs13020320

Chicago/Turabian Style

Granadeiro, José P., João Belo, Mohamed Henriques, João Catalão, and Teresa Catry. 2021. "Using Sentinel-2 Images to Estimate Topography, Tidal-Stage Lags and Exposure Periods over Large Intertidal Areas" Remote Sensing 13, no. 2: 320. https://0-doi-org.brum.beds.ac.uk/10.3390/rs13020320

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

Article Metrics

Back to TopTop