Deriving High Spatial-Resolution Coastal Topography From Sub-meter Satellite Stereo Imagery

High spatial resolution coastal Digital Elevation Models (DEMs) are crucial to assess coastal vulnerability and hazards such as beach erosion, sedimentation, or inundation due to storm surges and sea level rise. This paper explores the possibility to use high spatial-resolution Pleiades (pixel size = 0.7 m) stereoscopic satellite imagery to retrieve a DEM on sandy coastline. A 40-km coastal stretch in the Southwest of France was selected as a pilot-site to compare topographic measurements obtained from Pleiades satellite imagery, Real Time Kinematic GPS (RTK-GPS) and airborne Light Detection and Ranging System (LiDAR). The derived 2-m Pleiades DEM shows an overall good agreement with concurrent methods (RTK-GPS and LiDAR; correlation coefficient of 0.9), with a vertical Root Mean Squared Error (RMS error) that ranges from 0.35 to 0.48 m, after absolute coregistration to the LiDAR dataset. The largest errors (RMS error > 0.5 m) occurred in the steep dune faces, particularly at shadowed areas. This work shows that DEMs derived from sub-meter satellite imagery capture local morphological features (e.g., berm or dune shape) on a sandy beach, over a large spatial domain.


Introduction
Accurate topographic data are frequently needed for the assessment of rapid morphological changes and for the implementation of models that can predict coastal evolution. High spatial resolution coastal Digital Elevation Models (DEMs-defined here as the representation of the terrain surface elevations at regularly spaced intervals) are used to support vulnerability and risk assessment of a range of coastal hazards, such as beach erosion and sedimentation, storm surges, inundation, and sea level rise [1]. For such studies, the availability of a topographic dataset is fundamental, in particular for coastal systems characterized by a complex, rapidly evolving morphology.
Among topographic survey methods of suitable quality, those based on Global Navigation Satellite Systems (GPS), such as Real Time Kinematic GPS (RTK-GPS), have been used extensively to map and monitor coastal morphology [2]. Beach topographic surveys using RTK-GPS method can be performed either by walking and carrying a GPS receiver, or driving a mobile unit (e.g., quad bike). In both cases, the vertical precision is approximately 0.05 to 0.1 m, depending on the terrain relief [2]. This method typically requires an intense human effort, which normally is optimized by reducing the number of measurements to a limited number of cross-shore sections of the beach. Nevertheless, this limited spatial coverage results in an incomplete representation of topographic spatial patterns and evolving features, especially in the case of complex topographies such as steep and unconsolidated slopes. In such cases, interpolation methods are typically required, introducing additional uncertainty into the DEM [3].
Remote sensing techniques, such as airborne LiDAR (Light Detection and Ranging) and Unmanned Aerial Vehicle (UAV), emerge in this context as a solution to overcome the limited spatial coverage of the RTK-GPS method [4][5][6][7][8][9]. The use of airborne LiDAR to measure geomorphological changes in coastal areas is relatively new. This instrumentation acquires millions of x, y, z points per hour, with a horizontal spacing of typically 1 to 3 m. This high spatial resolution, together with the capacity to survey over large areas (from 10 1 to 10 5 m), allows overcoming traditional survey limitations found with RTK-GPS [2]. The vertical accuracy of LiDAR ranges from 0.05 m to 0.15 m [5], which is in the same order as RTK-GPS and appropriate for studying beach morphology. Nonetheless, LiDAR-based DEMs are costly [5,6], which limits the frequent (e.g., monthly or-post-storm) acquisition of large-scale topographic data adequate for the evaluation of coastal changes.
Airborne optical remote sensing and 3D-mapping have been serving the needs of regional-scale low-altitude imaging and geospatial information [10]. The enhanced usability of recent UAV equipment with onboard accurate positioning, such as off-the-shelf drones, has resulted in a large change in their practical application. The RTK-GPS positioning of the camera, combined with the large number of overlapping images, makes any additional ground surveys trivial. Moreover, the high degree of automation of UAVs and the absolute vertical precision, of approximately 0.2 m, achieved by the DEMs suggests possible uses in the fields of natural hazards, disaster response, and high-resolution terrain analysis [6]. Despite these advantages, a few disadvantages still remain such as the cost of the photogrammetric software and computer power that can be relatively high [7], the difficulty in removing dense vegetation to obtain bare earth elevation estimates [11], the need for electric batteries for longer flight duration, or the usage limitations related to weather conditions [12].
Sub-meter satellite imagery can potentially provide an alternative to these field-based techniques in order to collect high spatial resolution topographic data over large areas. The first civil satellite constellation that acquired stereoscopic imagery and applied DEM reconstruction over large areas was the French SPOT mission (Satellite Pour l'Observation de la Terre) in 1986 [12]. Since then, several very high spatial resolution satellites with stereo capabilities were launched in response to an increased demand [13]. Among them, the Pleiades constellation (built by the French Space Agency (CNES), commercialized by AIRBUS Defence & Space), consists of two high spatial resolution optical spacecrafts: Pleiades -1A and -1B. Both satellites fly over the same near-polar sun-synchronous orbits at an altitude of 694 km with a 180 • phase and descending node. The optical sensors of these satellites have the capability to obtain images with sub-meter image resolution (0.7 m pixel size, resampled to 0.5 m) over a maximum area of 350 km × 20 km (swath width of 20 km at nadir). An important aspect of Pleiades is the capacity to revisit any location in the world within 1 day, which is of great interest to monitor rapidly changing processes (e.g., coastal erosion due to storm events). Recent studies based on Pleiades-1A stereo-imagery include snow height mapping in mountainous areas [14], large landmass deformations due to earthquakes [15], surface reconstruction after landslides [16], and glacier topography [17,18].
The aim of the present work is to explore the use of Pleiades satellite stereo-imagery to develop a high resolution DEM of a 40-km-long sandy coastal section. The satellite-derived DEM is compared to RTK-GPS cross-shore profiles and an airborne LiDAR-derived DEM. The differences between the concurrent methods are quantified and the precision and accuracy of Pleiades-DEM analysed.

Study Site and Data Acquisition
A 40-km stretch of sandy coast in the South West of France was selected as the study site for the present work ( Figure 1). This section of the Aquitanian coast presents a relatively low shoreface bordered by aeolian dunes with an average crest elevation of about 15 m [19][20][21]. The sediment consists of fine to medium quartz, with mean grain sizes ranging from 200 to 400 µm [21]. This section of the coast is characterized by a macro-tidal regime, with an average tidal range of 3.2 m that can reach 5 m during spring tides [22]. The coast is exposed to high energy North Atlantic swells travelling mainly from the W-NW sector [23].

Pleiades Stereo-imagery Acquisition and DEM Generation
In the present experiment, stereo images were acquired when the angles between the line-of-sight of the satellite camera and the horizontal plane of the ground were 72° and 76°, for the first and second images respectively. The ground projection of the Pleiades position during the stereo-pair collection was 218 and 174 km from the coast respectively. Pleiades was overlooking the coastal area of interest from the sea side. This satellite setup indicates that both images were collected close to nadir angle (90° from the ground), with an azimuth of 19° for the first and −8° for the second image, resulting in an azimuth angle difference of 27° between the two images.
The Pleiades panchromatic band of the stereo pair was processed using the Ames Stereo Pipeline, ASP [24] to generate a DEM and ortho-images at 2 and 0.5 m resolution, respectively. The ASP uses the rational polynomial coefficient (RPC) camera model format for the DEM generation. The RPC model is provided in the imagery metadata (by AIRBUS) and gives a relationship between the image coordinates and the ground coordinates. No ground control points were initially used in the DEM generation. The planimetric coordinates were referenced to the WGS84 UTM 30N The Pleiades-HR 1A (hereinafter referred to as PL1A) stereo-pair was acquired on 14 November 2017 over a predefined area (orange box in Figure 1). The optical stereo-pair was obtained between 11:15 a.m. and 11:16 a.m. with a 40-s time-lag. The satellite orbits at 694 km altitude (base to height ratio of 0.36) and follows a descending orbit trajectory (North-South) in WGS84 decimal coordinates. A topographic DEM and ortho-image, covering the entire area of interest, was subsequently produced using NASA's AMES Stereo Pipeline [24].
In-situ RTK-GPS beach topographic measurements were used as ground-truth for inter-comparison with the Pleiades and LiDAR DEMs. The RTK-GPS survey was performed in the central section of the area of interest (coastal region of Capbreton) between the 7 th and 9 th of November 2017 ( Figure 1). The average tidal range during the beach surveys was 3.4 m with a moderate wave climate (wave height under 3 m). Beach profiles were surveyed in continuous mode (waypoint every 1 second using a position dilution of precision-PDOD-mask of 3) from the waterline to the back of the frontal dune ridge along discrete cross-shore profiles spaced by approximately 250 m ( Figure 1). Note that due to radio transmission shadowing (resulting in no real time correction) or ground obstacles (e.g., fences, walls, etc.), the survey coverage presents some spatial irregularities ( Figure 1). In addition to the RTK-GPS dataset, a high resolution airborne-LiDAR topographic survey and orthophoto map (aerial imagery was orthorectified with the LiDAR observations) of the entire SW coast of France (performed by the Institut National de L'Information Geographic et Forestiere-IGN; in cooperation with the Bureau de Recherches Géologiques et Minières-BRGM) were acquired in October 2017.

Pleiades Stereo-imagery Acquisition and DEM Generation
In the present experiment, stereo images were acquired when the angles between the line-of-sight of the satellite camera and the horizontal plane of the ground were 72 • and 76 • , for the first and second images respectively. The ground projection of the Pleiades position during the stereo-pair collection was 218 and 174 km from the coast respectively. Pleiades was overlooking the coastal area of interest from the sea side. This satellite setup indicates that both images were collected close to nadir angle (90 • from the ground), with an azimuth of 19 • for the first and −8 • for the second image, resulting in an azimuth angle difference of 27 • between the two images.
The Pleiades panchromatic band of the stereo pair was processed using the Ames Stereo Pipeline, ASP [24] to generate a DEM and ortho-images at 2 and 0.5 m resolution, respectively. The ASP uses the rational polynomial coefficient (RPC) camera model format for the DEM generation. The RPC model is provided in the imagery metadata (by AIRBUS) and gives a relationship between the image coordinates and the ground coordinates. No ground control points were initially used in the DEM generation. The planimetric coordinates were referenced to the WGS84 UTM 30N coordinate system and the heights were computed above the WGS84 ellipsoid. The Pleiades DEM and ortho-images were a posteriori coregistered by applying a first order polynomial transformation defined by 37 concomitant points manually identified in the Pleiades-ortho and DEM and in the IGN/BRGM orthophoto map (used as the reference). This process was performed with QGIS software and Georeferencer GDAL plugin, with an average planimetric error of 0.5 m. To convert Pleiades altimetric data from WGS1984 to the French NGF-IGN vertical datum, the Pleiades elevations were corrected from the average difference with the LiDAR elevations (Pleiades elevations 12.2 m higher) determined at the same point locations and then used for the planimetric correction.

RTK-DGPS Topographic Survey
The planimetric coordinates of the RTK-DGPS topographic survey were referenced to the World Geodetic System (WGS84) while the vertical datum was referenced to NGF-IGN 1969 datum. A GPS base station was installed in a local geodesic point (located near Capbreton - Figure 1) and provided, in real time, the corrections to the mobile GPS unit via radio-transmission. After the survey, the planimetric coordinates were converted to the same coordinate system as the Pleiades products (WGS 84 UTM 30N), and spikes in the data (erroneous measurements) were eliminated. The processed topographic measurements were subsequently divided into individual profiles and interpolated in the cross-shore direction with 2 m spacing (to match the resolution of the Pleiades DEM). Comparisons between RTK-GPS ground-truth and remotely sensed Pleiades and LiDAR DEMs were performed by extracting values from the DEMs at each profile location (using all RTK-GPS point measurements). This task was performed using QGIS software (Lyon version) and the function "sample raster maps at point location" from GRASS-GIS toolbox. The data comparisons included the calculation of statistical parameters such as the correlation coefficient (CC), root mean squared (RMS) error, and bias (BIAS) using all RTK-GPS topographic observations. Remote Sens. 2019, 11, 590 5 of 11

Airborne LiDAR 3D Topographic Survey
As part of a regional coastal monitoring program, the Aquitaine coastal zone is surveyed every year with airborne LiDAR by the IGN in cooperation with the BRGM. In the present work, the LiDAR topographic survey of the Aquitaine coast performed between 4 and 7 of October 2017 was used for comparison with the RTK-GPS and Pleiades observations. Note that the LiDAR survey was executed approximately a month before the field campaign and Pleiades acquisition. The airborne-LiDAR survey was performed with the Leica ALS70-HP LiDAR, mounted on an aircraft, and acquiring topographic information with a density of 8 points per square meter. The planimetric coordinates of the LiDAR point-cloud were referenced to the Lambert 93 coordinate system, with a precision of 30 cm, and terrain topography referenced to the NGF-IGN 1969 datum, with an accuracy of 15 cm (information provided by the LiDAR survey metadata). Simultaneously with the LiDAR acquisition, high resolution (10 cm) aerial orthophotography was obtained with an 8-head IGN V2 (focal length of 135 mm) camera. The position and orientation of the images were obtained from the GPS and inertial sensors embedded in the aircraft. Two products were obtained from this flight: a DEM of the study area, with a spatial resolution of 1 m, and an orthophoto map with 10 cm resolution. For comparison with the other datasets, the coordinate systems of these LiDAR products were converted to WGS 84 UTM 30N. Figure 2 shows the vertical difference between the three concurrent survey methods over 4138 points (i.e., the number of data points measured during the RTK-GPS survey). The differences between remote sensing methods (Pleiades and LiDAR) and RTK-GPS are normally distributed, with mean differences (BIAS) of 0.01 m and 0.03 m, and RMS errors of 0.35 m and 0.37 m for Pleiades and LiDAR respectively (both Pleiades and LiDAR elevations are slightly higher than the RTK measurements). The observed slight mean difference between the RTK-GPS and remote sensing methods are within the accuracy of the RTK-GPS, thus indicating that the different methods have similar accuracy. It is important to note that even though the LiDAR data have a larger RMS error than the Pleiades one, the error distribution is narrower and skewed for negative values. The 1:1 scatter-comparison of the surveys shown in Figure 2 indicates that the remotely-sensed beach topography is highly correlated with the RTK-GPS observations (slope = 1.01 and CC = 0.99 for both Pleiades and LiDAR). It is also observed that the correlation with Pleiades values does not vary with elevation (extending from the back of the dune to the top of the swash zone) while for LiDAR some scattering is observed in the lower part of the beach (beach face).

Comparison between RTK-GPS, Pleiades and LiDAR Topography
A close inspection of the datasets indicates the scatter (observed in the beach face) between RTK-GPS and LiDAR topography was due to morphological changes (berm erosion- Figure 3) over the beach profile that occurred between the two surveys. During this period of the year (winter season), this coastal area is under energetic waves and significant morphological changes (> 1 m) in the beach and dune face are likely to occur [23]. Considering that the RTK-GPS survey and Pleiades image acquisition were days apart, small morphological changes likely occurred in the beach face, resulting in minor differences between the two datasets ( Figure 3).
Remote Sens. 2019, 11, x FOR PEER REVIEW 6 of 12 elevation (extending from the back of the dune to the top of the swash zone) while for LiDAR some scattering is observed in the lower part of the beach (beach face). A close inspection of the datasets indicates the scatter (observed in the beach face) between RTK-GPS and LiDAR topography was due to morphological changes (berm erosion- Figure 3) over the beach profile that occurred between the two surveys. During this period of the year (winter season), this coastal area is under energetic waves and significant morphological changes (> 1 m) in the beach and dune face are likely to occur [23]. Considering that the RTK-GPS survey and Pleiades image acquisition were days apart, small morphological changes likely occurred in the beach face, resulting in minor differences between the two datasets ( Figure 3).  Figure 4 shows the two DEMs produced using Pleiades stereo imagery and LiDAR data, together with the difference between the two DEMs (Pleiades DEM minus LiDAR DEM). The northern part of the study area (North of Capbreton- Figure 4) presents higher dune fields than the South, and Pleiades DEM was able to capture this spatial variability with the same quality as the LiDAR. The difference between the two DEMs shows that over the full domain, 70% of the difference lies within ± 0.5 m. Areas with positive elevation difference (i.e., Pleiades higher than LiDAR) are located between the frontal dune face and swash zone, while a negative difference was found more often at the back of the dune. A preliminary inspection of the alongshore error distribution allowed identifying larger differences at the dune face region between Capbreton and  A close inspection of the datasets indicates the scatter (observed in the beach face) between RTK-GPS and LiDAR topography was due to morphological changes (berm erosion- Figure 3) over the beach profile that occurred between the two surveys. During this period of the year (winter season), this coastal area is under energetic waves and significant morphological changes (> 1 m) in the beach and dune face are likely to occur [23]. Considering that the RTK-GPS survey and Pleiades image acquisition were days apart, small morphological changes likely occurred in the beach face,  Figure 4 shows the two DEMs produced using Pleiades stereo imagery and LiDAR data, together with the difference between the two DEMs (Pleiades DEM minus LiDAR DEM). The northern part of the study area (North of Capbreton- Figure 4) presents higher dune fields than the South, and Pleiades DEM was able to capture this spatial variability with the same quality as the LiDAR. The difference between the two DEMs shows that over the full domain, 70% of the difference lies within ± 0.5 m. Areas with positive elevation difference (i.e., Pleiades higher than LiDAR) are located between the frontal dune face and swash zone, while a negative difference was found more often at the back of the dune. A preliminary inspection of the alongshore error distribution allowed identifying larger differences at the dune face region between Capbreton and  Figure 4 shows the two DEMs produced using Pleiades stereo imagery and LiDAR data, together with the difference between the two DEMs (Pleiades DEM minus LiDAR DEM). The northern part of the study area (North of Capbreton- Figure 4) presents higher dune fields than the South, and Pleiades DEM was able to capture this spatial variability with the same quality as the LiDAR. The difference between the two DEMs shows that over the full domain, 70% of the difference lies within ± 0.5 m. Areas with positive elevation difference (i.e., Pleiades higher than LiDAR) are located between the frontal dune face and swash zone, while a negative difference was found more often at the back of the dune. A preliminary inspection of the alongshore error distribution allowed identifying larger differences at the dune face region between Capbreton and Labenne ( Figure 4) in comparison to the rest of the domain. Within the section with these particular large differences, it was possible to identify the presence of shadows (for each transect, the length of the shadow was manually digitized from the orthophotomap) in the dune face that coincided with the areas where the largest errors were observed ( Figure 5).

3D Beach Topography Comparisons
Specific ground characteristics, such as the slope and aspect [25,26], can have an indirect impact on the remotely-sensed DEM quality in regions with high relief, such as dunes. The presence of shadows in the dune face is determined by the steepness (slope) and orientation (aspect) of the topography in relation to the Sun light (for a given Sun altitude and azimuth). Optical remote sensing images from shadowed areas have low reflectance and texture, which alters the calculation of the disparity (which is later converted into elevation) of corresponding points in the stereo-pair and leads to less accurate elevations [27,28].
In order to investigate the effect of the alongshore variability in the dune slope (front and back faces) and of the shadowed areas on the Pleiades DEM error, 266 cross-shore transects spaced 150 m apart were created. For each transect, the RMS error of the elevation difference (between Pleiades and LiDAR DEMs) over the dune face and back of the dune was computed. Figure 6 shows the RMS error variation along the 40 km of measured coastline together with the variations in the dune slope and shadow length in the dune face. The RMS error shows significant variations between transects, with values of similar magnitude on the dune face and back of the dune. The exception, as noted in Figure 4, is the coastal stretch between Capbreton and Labenne where a peak of RMS error (>0.5 m) is identified on the dune face. The alongshore location of this peak coincides with a relatively steep dune face region (large slope values) and the largest shadow regions (Figure 6). Steep dune face slopes were also present in the North section of the study area; however, the errors in this region were within the average, suggesting that the dune slope by itself cannot explain the largest errors observed. Shadowing at the dune face is more likely to represent a decisive quality factor of the produced DEM, considering the strong correlation between the RMS error and shadowing (CC = 0.77). the areas where the largest errors were observed ( Figure 5).
Specific ground characteristics, such as the slope and aspect [25,26], can have an indirect impact on the remotely-sensed DEM quality in regions with high relief, such as dunes. The presence of shadows in the dune face is determined by the steepness (slope) and orientation (aspect) of the topography in relation to the Sun light (for a given Sun altitude and azimuth). Optical remote sensing images from shadowed areas have low reflectance and texture, which alters the calculation of the disparity (which is later converted into elevation) of corresponding points in the stereo-pair and leads to less accurate elevations [27,28].
In order to investigate the effect of the alongshore variability in the dune slope (front and back faces) and of the shadowed areas on the Pleiades DEM error, 266 cross-shore transects spaced 150 m apart were created. For each transect, the RMS error of the elevation difference (between Pleiades and LiDAR DEMs) over the dune face and back of the dune was computed. Figure 6 shows the RMS error variation along the 40 km of measured coastline together with the variations in the dune slope and shadow length in the dune face. The RMS error shows significant variations between transects, with values of similar magnitude on the dune face and back of the dune. The exception, as noted in Figure 4, is the coastal stretch between Capbreton and Labenne where a peak of RMS error (> 0.5 m) is identified on the dune face. The alongshore location of this peak coincides with a relatively steep dune face region (large slope values) and the largest shadow regions (Figure 6). Steep dune face slopes were also present in the North section of the study area; however, the errors in this region were within the average, suggesting that the dune slope by itself cannot explain the largest errors observed. Shadowing at the dune face is more likely to represent a decisive quality factor of the produced DEM, considering the strong correlation between the RMS error and shadowing (CC = 0.77).    Figure 7 show the statistical evaluation of the comparison between the Pleiades and LiDAR elevations extracted over the 266 cross-shore transects ( Figure 6). Data from the beach face, containing natural morphological changes (not related to the method), were removed from this statistical comparison. The error distribution indicates that differences between the LiDAR and Pleiades DEMs over the full domain are normally distributed with a mean difference of -0.015 m and RMS error of 0.48 m (LiDAR DEM is slightly higher). The 1:1 comparison (Figure 7) shows that Pleiades and LiDAR beach topographies are highly correlated (slope = 1.01 and CC = 0.99).   Figure 7 show the statistical evaluation of the comparison between the Pleiades and LiDAR elevations extracted over the 266 cross-shore transects ( Figure 6). Data from the beach face, containing natural morphological changes (not related to the method), were removed from this statistical comparison. The error distribution indicates that differences between the LiDAR and Pleiades DEMs over the full domain are normally distributed with a mean difference of -0.015 m and RMS error of 0.48 m (LiDAR DEM is slightly higher). The 1:1 comparison (Figure 7) shows that Pleiades and LiDAR beach topographies are highly correlated (slope = 1.01 and CC = 0.99).

The Use of Pleiades to Survey and Monitor Coastal Areas
The comparisons presented in sections 4.1 and 4.2 indicate that Pleiades DEM reaches a vertical accuracy (BIAS) and precision (RMS error) similar to state-of-the-art survey methodologies used in coastal areas such as the RTK-GPS and LiDAR. Although in this study a comparison between Pleiades and UAVs stereo-topographic DEM was not presented, there are a few differences between these two methods that are important to take into consideration when deciding which of the methods to use. Both methods obtain the DEM based on photogrammetric techniques that estimate a 3D point-cloud of the ground surface using a large number of matching object and textural features automatically detected in the overlapping images. The low altitude and capacity to obtain images from many different view angles allows UAVs to obtain detailed 3D reconstruction of the ground surface with centimetric spatial resolution [7,29]. These characteristics represent an important advantage of the UAV technique when surveying highly complex coastal features (e.g., irregular rocky coastlines; [30]) or avoid the influence of sunlight exposure in the target features, such as the presence of shadows.
As it was observed in the present work, the presence of shadows can affect the quality of the DEM from stereo-satellite imagery ( Figure 5 and 6). A similar problem was identified by [31] when validating the Pleiades tri-stereo digital surface model on an urban area. The fact that Pleiades, stereo and tri-stereo, only obtains images along a fixed trajectory, this limits the view angles of the ground surface (one backward looking, one forward looking, plus a third near-nadir image, in the tri-stereo configuration), making difficult an appropriate DEM estimation from features affected by shadows.
Despite this limitation, the present results show that Pleiades stereo imagery has the capability to capture local beach features such as berms or dunes crests and troughs over large domains, highlighting its incomparable advantage to any other existing methodology. These unique skills allow the Pleiades constellation to potentially overcome traditional survey challenges in coastal areas, such as the acquisition stereo-imagery measurements in large coastal segments within a short period of time (i.e., minutes). For traditional survey methods, this would mean a gigantic logistical challenge and expensive exercise. Furthermore, Pleiades' capacity to acquire imagery anywhere on the globe within 1 day is ideal for rapid response assessment of changes in the coastal zone (e.g., to assess morphological changes after storms).
A potential drawback of the Pleiades constellation is the on-demand availability, which is different from other observation missions like Landsat (NASA) or Sentinel (ESA) that acquire optical imagery on a regular basis without any previous request. In addition to this, another potential limitation of our processing flow is the dependence of ground control points to correct the vertical

The Use of Pleiades to Survey and Monitor Coastal Areas
The comparisons presented in Sections 4.1 and 4.2 indicate that Pleiades DEM reaches a vertical accuracy (BIAS) and precision (RMS error) similar to state-of-the-art survey methodologies used in coastal areas such as the RTK-GPS and LiDAR. Although in this study a comparison between Pleiades and UAVs stereo-topographic DEM was not presented, there are a few differences between these two methods that are important to take into consideration when deciding which of the methods to use. Both methods obtain the DEM based on photogrammetric techniques that estimate a 3D point-cloud of the ground surface using a large number of matching object and textural features automatically detected in the overlapping images. The low altitude and capacity to obtain images from many different view angles allows UAVs to obtain detailed 3D reconstruction of the ground surface with centimetric spatial resolution [7,29]. These characteristics represent an important advantage of the UAV technique when surveying highly complex coastal features (e.g., irregular rocky coastlines; [30]) or avoid the influence of sunlight exposure in the target features, such as the presence of shadows.
As it was observed in the present work, the presence of shadows can affect the quality of the DEM from stereo-satellite imagery (Figures 5 and 6). A similar problem was identified by [13] when validating the Pleiades tri-stereo digital surface model on an urban area. The fact that Pleiades, stereo and tri-stereo, only obtains images along a fixed trajectory, this limits the view angles of the ground surface (one backward looking, one forward looking, plus a third near-nadir image, in the tri-stereo configuration), making difficult an appropriate DEM estimation from features affected by shadows.
Despite this limitation, the present results show that Pleiades stereo imagery has the capability to capture local beach features such as berms or dunes crests and troughs over large domains, highlighting its incomparable advantage to any other existing methodology. These unique skills allow the Pleiades constellation to potentially overcome traditional survey challenges in coastal areas, such as the acquisition stereo-imagery measurements in large coastal segments within a short period of time (i.e., minutes). For traditional survey methods, this would mean a gigantic logistical challenge and expensive exercise. Furthermore, Pleiades' capacity to acquire imagery anywhere on the globe within 1 day is ideal for rapid response assessment of changes in the coastal zone (e.g., to assess morphological changes after storms).
A potential drawback of the Pleiades constellation is the on-demand availability, which is different from other observation missions like Landsat (NASA) or Sentinel (ESA) that acquire optical imagery on a regular basis without any previous request. In addition to this, another potential limitation of our processing flow is the dependence of ground control points to correct the vertical offset of the DEM and geometric inaccuracies of the raw data. This aspect limits the use of Pleiades DEM for applications where the absolute error in relation to a specific datum is required. Nevertheless, for cases where the offset to the local vertical datum is secondary, the inter-comparison between consecutive Pleiades DEMs is likely to result in errors within the precision of the method.

Conclusions
The beach topography of a 40-km-long sandy coastal stretch, located in the southwest of France, was surveyed by satellite Pleiades stereo-imagery. The computed DEM was compared with those obtained from traditional survey methodologies (RTK-GPS and LiDAR). Present findings indicate that Pleiades stereo-imagery allows the acquisition of high resolution DEM with a RMS error that ranges from 0.35 to 0.48 m. The largest errors were observed at the dune face, in regions with large shadow patches. Near-perfect agreement between Pleiades and concurrent methods (all computed CC were above 0.9) provides strong indications that this method can be used as a surveying tool to monitor detailed coastal morphological changes over large spatial domains.