# Monitoring the Snowpack Volume in a Sinkhole on Mount Lebanon using Time Lapse Photogrammetry

C. Abou Chakra, Simon Gascoin, J. Somma, Pascal Fanise, Laurent Drapeau

► To cite this version:

C. Abou Chakra, Simon Gascoin, J. Somma, Pascal Fanise, Laurent Drapeau. Monitoring the Snow-pack Volume in a Sinkhole on Mount Lebanon using Time Lapse Photogrammetry. Sensors, 2019, 19 (18), pp.1-16. 10.3390/s19183890. [hal-02351618](https://hal.science/hal-02351618)

HAL Id: hal-02351618

https://hal.science/hal-02351618

Submitted on 25 May 2020

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

Article

# Monitoring the Snowpack Volume in a Sinkhole on Mount Lebanon using Time Lapse Photogrammetry

Charbel Abou Chakra<sup>1,2,3,*</sup>, Simon Gascoin<sup>4,*</sup>①, Janine Somma<sup>1</sup>, Pascal Fanise<sup>4</sup> and Laurent Drapeau<sup>4</sup>

Received: 15 July 2019; Accepted: 6 September 2019; Published: 9 September 2019

<figure><img src="image_3.png" /></figure>

check for
updates

**Abstract:** Lebanon has experienced serious water scarcity issues recently, despite being one of the wealthiest countries in the Middle East for water resources. A large fraction of the water resources originates from the melting of the seasonal snow on Mount Lebanon. Therefore, continuous and systematic monitoring of the Lebanese snowpack is becoming crucial. The top of Mount Lebanon is punctuated by karstic hollows named sinkholes, which play a key role in the hydrological regime as natural snow reservoirs. However, monitoring these natural snow reservoirs remains challenging using traditional in situ and remote sensing techniques. Here, we present a new system in monitoring the evolution of the snowpack volume in a pilot sinkhole located in Mount Lebanon. The system uses three compact time-lapse cameras and photogrammetric software to reconstruct the elevation of the snow surface within the sinkhole. The approach is validated by standard topographic surveys. The results indicate that the snow height can be retrieved with an accuracy between 20 and 60 cm (residuals standard deviation) and a low bias of 50 cm after co-registration of the digital elevation models. This system can be used to derive the snowpack volume in the sinkhole on a daily basis at low cost.

**Keywords:** snow; time-lapse camera; photogrammetry; structure-from-motion; water resource

## 1. Introduction

Lebanon is one of the wealthiest countries in the Middle East for water resources. Despite this, the people in Lebanon are frequently affected by water restrictions, especially during the dry summer season. In addition, the water demand is increasing, due to the growth of the native population and the large number of Middle Eastern refugees received in recent years [1].

In Lebanon, precipitation is enhanced by the Lebanon and Anti-Lebanon mountain ranges, which run parallel to the Mediterranean coast. An important fraction of the precipitation falls as snow on the Lebanese mountains, and the snow covers around 25% of the total Lebanese territory in winter [2–5]. Furthermore, it is estimated that approximately 30% of the national water resource originates from seasonal snow melt [4]. However, the seasonal evolution of the snowpack remains poorly known, as there is currently no operational snow national monitoring program [6,7]. Continuous and systematic snow observations are urgently required to improve water management planning [6,7].

The upper area of Mount Lebanon is a large upper cretaceous carbonate plateau [8]. A unique aspect of Mount Lebanon snow hydrology is related to its karstic features like sinkholes [9]. These

<footnote n="1">¹ Laboratoire de Télédétection, Centre de Recherche en Environnement-Espace Méditerranée Orientale, Université Saint-Joseph, Beirut BP 17-5208, Lebanon</footnote>

<footnote n="2">² Arab Union of Surveyors, Beirut BP 9300, Lebanon</footnote>

<footnote n="3">³ Faculty of Engineering, Université Libano-Canadienne, Aintoura BP 32, Lebanon</footnote>

<footnote n="4">⁴ Centre d'Etudes Spatiales de la Biosphère (CESBIO), Université de Toulouse, CNES/CNRS/INRA/IRD/UPS, 31401 Toulouse, France</footnote>

<footnote>* Correspondence: charbel.abouchakra@gmail.com (C.A.C.); simon.gascoin@cesbio.cnes.fr (S.G.); Tel.: +961-7090-7922 (C.A.C.)</footnote>
basin-shaped hollows typically have a diameter of a few tens of meters or more, and a depth of a few meters. These sinkholes are closed and form natural collectors for rainfall and meltwater for groundwater recharge. In addition, due to their hollow shape, sinkholes on Mount Lebanon form natural traps for snow for several reasons [3,6]. First, sinkholes are preferential areas of wind-drifted snow deposition, due to the localized drop in wind speed on their lee side. Second, they isolate the snowpack from the surface atmospheric turbulence, which reduces the melting and sublimation rates. As a result, sinkholes on Mount Lebanon help delay the low flow period in downstream areas by providing meltwater late in summer [6]. This significant hydrological contribution has motivated the implementation of an observation system that monitors the evolution of the seasonal snowpack in a sinkhole.

From a water management perspective, the snow water equivalent (SWE) is the key variable to monitoring [10]. However, existing automatic SWE measurement devices are costly or difficult to install in rugged terrain. On the other hand, SWE variability is closely related to the variability of snow depth, therefore, snow depth monitoring is useful in hydrological studies [11,12]. In addition, bulk snow density can be simulated with reasonable uncertainty using empirical formulas or a snowpack model [13]. However, snow depth can be measured using sonic gauges in automatic weather stations. As sonic snow gauges provide only a point-scale measurement, given the limited footprint of the reflected echoes, this kind of device is not adapted to snow monitoring in sinkholes over Mount Lebanon where snow depth can vary from 0 to 15 m within a few tens of meters, due to their specific morphology. Spatially-distributed and accurate observations of snow depth can be obtained at low cost from drone photogrammetry [14]. However, the temporal frequency of drone surveys in Lebanon is often limited by the meteorological conditions or military restrictions. Another option is to use very high resolution stereo satellite images [15]. However, the acquisition of stereoscopic images has proven to be difficult due to the frequent cloud cover in winter and the fact that very high resolution satellites (in this case, Pléiades) are often tasked by the military, who have priority over the region. Finally, terrestrial lidar scanners are very efficient to determine the snow surface topography, but these devices are very costly and cannot be deployed in the field for continuous monitoring [16].

Terrestrial time-lapse photography has been increasingly used in environmental studies, including for snow cover monitoring in small- and medium-sized catchments [17–19]. This new technique is not expensive and can be easily implemented in mountain environments for snow studies [20,21]. However to the best of our knowledge, this method has only used to determine the extent of the snow cover and other derived metrics, like snow cover duration or the snow melt-out date. Revuelto et al. [21] used the snow cover area retrieved from time-lapse camera images to perform backward snow depth reconstruction. However, this approach relied on a degree-day model that was subject to uncertainties and required local weather data. On the other hand, photogrammetry is now widely used to generate three-dimensional (3D) models of the land surface, using multi-view images from digital cameras, in particular due to the development of structure-from-motion (SfM) algorithms [22–26]. For instance, Mertes et al. [27] used SfM from oblique terrestrial and aerial images to map glacier elevation changes.

Here, we present a new system using low cost digital cameras to monitor the snowpack volume in a sinkhole in Mount Lebanon. Three time lapse cameras were employed to derive a 3D model of the snow surface on a daily basis. We also evaluate the method using topographic surveys with a total station.

## 2. Study Area

The study area is located on the high plateau of Mount Lebanon (Figure 1). This karstic Cenomanian geological formation is punctuated with hundreds of sinkholes at elevation, ranging between 1700 and 3000 m. In this area, the seasonal snow can persist over six months [4]. For this study we focus on a pilot site at 2300 m in the Jabal el Dib region (Ouyoun el Simane-Kfardebian). This site is a sinkhole of about 16 m depth, 100 m wide, and 170 m long, located near the ski resort of Mzaar Kfardebian
(35°51′28″ E; 33′58′53″ N). The bottom of the hole has an oval shape, with a diameter ranging from 30 to 36 m. There is a very sparse vegetation cover in the studied sinkhole (herbaceous and small shrubs), but most of the surface is covered with rocks and bare soil, including the floor and the walls of the sinkhole. Since 2010, the site has been equipped with an automatic weather station [7]. There is also a geodetic point very close to the sinkhole, which, together with the buildings of the ski center, facilitates geodetic surveys.

<figure><img src="image_4.png" /><figcaption>Figure 1. Location map.</figcaption></figure>

## 3. Method

### 3.1. *System Setup*

For this study, three cameras were installed to image a surface area of about 4500 m², with stereoscopic overlap around 80 % (Figure 2). The cameras were installed facing the dominant wind direction, in order to face the snow cornice, which is formed every year by wind drift (Figure 2). Each camera was mounted on 2 m height galvanized pole, fixed in a concrete block, which was poured within the soil (Figure 2). The system was designed to resist harsh conditions, including strong gusts in winter and to require minimal maintenance.

The spacing distance between each camera was determined empirically by taking a dozen photographs of the site, from different location, and testing different combinations. We found that an approximate spacing of 5 m between each camera allowed a good resolution of the surface topography. The cameras were set up with a downward inclination of the focal plane of about 30° from the zenith. We used three Canon PowerShot A1400 cameras (Konka Group Co., Ltd., Shenzhen, China), with a resolution of 4608 by 3456 pixels and a focal length of 5 mm. The cameras were connected to a Campbell CR200 data-logger and powered by a 12 V lithium battery. The CR200 was only used to switch on and off the cameras at fixed times. A DC-to-DC converter was installed within each camera case to convert the 12 V voltage from the CR200 to 3.6 V. The Canon Hack Development Kit (https://chdk.fandom.com/) was used to trigger photograph acquisition in automatic mode when the camera was switched on. The photographs were saved in the 32 Gb SD cards of each camera. The system was programmed to take three triplets of photographs per day at 08:00, 10:00, and 12:00 local time (UTC + 2 h). Morning acquisitions were preferred given the westward orientation of the cameras. The acquisition period started on 15 March 2015, and the system worked smoothly until complete snow ablation on 11 July 2015. However for this study, we only processed the images that were acquired on six different days during this period, when field validation surveys were conducted (see below).
<figure><img src="image_5.png" /><figcaption>**Figure 2.** Setup: (a) Terrestrial view of the setup; (b) Drone ortho-image of the study area, taken on 24 April 2018 with 100 m as flight altitude.</figcaption></figure>

## 3.2. Photogrammetric Processing

The snow volume in the sinkhole was computed by differencing a digital elevation models (DEM) of the snow-covered scene (snow-on DEM) with a DEM of same area without snow (snow-off DEM). To complete this operation, it was important that both DEMs are well aligned (co-registered). Hence, we defined two experiments, which differ by the method in co-registering the snow-on with the snow-off DEM. In addition, we also evaluated the impact of the source imagery to compute the snow-off DEM. The three experiments (Exp. 1, Exp. 2 and Exp. 3) are detailed below and summarized in Table 1.

<table><caption>**Table 1.** Characteristics of the three experiments to compute the snow-on surface elevation (GCP: ground control point, ICP: iterative closest point).</caption><thead><tr><td></td><td>Snow-off DEM Source Imagery</td><td>Snow-on DEM Source Imagery</td><td>Method to Register Snow-on DEM</td></tr></thead><tbody><tr><td>Exp. 1</td><td>cameras</td><td>cameras</td><td>GCP</td></tr><tr><td>Exp. 2</td><td>cameras</td><td>cameras</td><td>ICP</td></tr><tr><td>Exp. 3</td><td>drone</td><td>cameras</td><td>ICP</td></tr></tbody></table>

To compute the DEMs, we used Agisoft Photoscan Professional to generate a time-stamped DEM, from three simultaneous photographs, taken by the time-lapse cameras. We selected Agisoft Photoscan because it is a state-of-the-art photogrammetric. The software can also run in batch mode, which integrates the SfM workflow with a user-friendly interface, but also because it can be used in batch mode through a Python application programming interface (API) to automate future processing. Numerous studies have successfully utilized Agisoft Photoscan in this manner [22–26,28–30]. We used the images that were acquired at 08:00 local time. Agisoft was configured to produce a DEM at 10 cm horizontal resolution. The coordinates of the cameras were surveyed with a dGPS and provided as input data to Agisoft. For each date, this process also allowed the generation of an ortho-rectified image of the scene from the photographs of the three cameras using the maximum of intensity at each pixel among the three images. This ortho-image was used in the co-registration step below.

We used six ground control points (GCPs) in the scene to generate the snow-off DEM with Agisoft for all the dates as it is generally done to allow a good registration of the DEM to the spatial reference system of the study area (Exp. 1). The GCPs (Figure 3) were measured by total station (Pentax R-1505N)
The computation of the differences between the co-registered snow-on DEMs and snow-off DEM was done using the geo-diff tool in ASP (step 4). This gave us a raster map of the elevation differences by pixel (dDEM), which can be compared to field measurements of the snow depth. Finally, the snow volume in the sinkhole was computed by spatial integration of the dDEM within the snow cover area polygon. Table 1 summarizes the different experiments that were conducted to compute the snow-on surface elevation for six validation dates.

### 3.3. Ground Truth Measurements

The elevation of the snow surface in the study area was surveyed on 20 May 2015, 4 June 2015, 18 June 2015, 29 June 2018, and 4 July 2015 (Figure 4) using a total station (Leica TS06). The total stations have been, until now, one of the main instruments used for obtaining high accuracy surveying measurements, including topography, hydrography, construction, geography, hydrology, and many other fields [35,36]. As an example, the total stations have also validated the accuracy of the United States geological digital elevation models (DEMs) [35].

<figure><img src="image_6.png" /><figcaption>**Figure 4.** Profiles of the snow surface extracted from the DEM produced from the total station data total station.</figcaption></figure>

Although, measurements with a total station do not provide a dense cloud of the surveyed area, it is possible to obtain centimeter to millimeter accuracy of any observed point position. Hence, the points obtained by the total station were used as reference altitudes to validate the DEM that were obtained from the terrestrial photogrammetry. On each date, the random points were measured on the snow cover. The criteria for choosing the points locations was a space of around 5 m for a flat area; the points were densified to approximately 1 m, when a change of slope was noticed. In other words, the locations of the points were randomly selected, but the spacing between he points varied, depending on the topography (closer points where the snow surface was more variable). The number of measured points for each date is shown in Table 2. For the comparison of the total station DEM with the camera DEM, the total station points were converted to a raster, with the same spatial resolution of 10 cm. The points were first converted to triangular irregular network (TIN) and then to a regular elevation grid in ArcMap software (Figure 5).
<table><caption>**Table 2.** Number of points measured by the total station for each surveyed day.</caption><thead><tr><td>Date</td><td>Number of Points</td></tr></thead><tbody><tr><td>8 May 2019</td><td>210</td></tr><tr><td>20 May 2019</td><td>250</td></tr><tr><td>4 June 2019</td><td>219</td></tr><tr><td>18 June 2019</td><td>140</td></tr><tr><td>29 June 2019</td><td>168</td></tr><tr><td>4 July 2019</td><td>92</td></tr></tbody></table>

<figure><img src="image_7.png" /><figcaption>**Figure 5.** Steps to obtain a raster DEM from total station points.</figcaption></figure>

During the same topographic surveys, the limit of the snow cover was marked using a specific code when acquiring the point coordinates on the edge of the snow patch. These points were assembled as a single polygon to delineate the snow limit.

## 4. Results

We present here the results of the validation for each experiment (Table 1). The snow-on DEMs were compared to the total station measurements. Figure 6 shows the spatial variability in the elevation residuals and Table 3 provides the statistical description of the residuals. The results indicate that the snow depth was retrieved with accuracy (as expressed by the standard deviation of the residuals) between 0.08 and 1.4 m, depending on the date and the method (median value 0.31 m). Except for the first two dates, the absolute bias is generally below 1 m for all experiments (median of the absolute bias values: 0.35 m). The bias is below 50 cm in 10 out of 17 cases. The root mean squared error (RMSE) values on the snow surface elevation range between 2.63 m for Exp. 2 DEM, that was conducted on 5 May 2015, to 0.09 m for the Exp. 1, which was conducted on 4 July 2015. Similarly, in relation to bias, the large range of RMSE is mainly explained by the date of the experiment (Figure 7). In particular, low RMSE values were obtained at the end of the snow season (4-Jun-2015, 29-Jun-2015, and 4-July-2015). For Exp. 2 and Exp. 3, a reduction of the RMSE through the snow season was observed on both, the stable area elevation and the snow surface elevation (Figure 8). Figure 6 shows that the spatial distribution of the error is consistent between Exp. 2 and Exp. 3, which suggests that the elevation errors are primarily due to the snow-on DEMs, since only the snow-off DEMs differ between both experiments.
<figure><img src="image_8.png" /><figcaption>**Figure 6.** Maps of the elevation difference between total station data and camera DEM for each date and experiment.</figcaption></figure>
<table><caption>**Table 3.** Statistics of the residuals between the snow-on DEMs and the field measurements. The RMSE over the stable area is provided as an indicator of the performance of the ICP adjustment. The results are not available for Exp. 1 because the GCPs were covered by the snow and thus not visible in the photographs.</caption><thead><tr><th>Date</th><th>Exp.</th><th>Min (m)</th><th>Max (m)</th><th>Mean (m)</th><th>Std (m)</th><th>RMSE (m)</th><th>Stable Area RMSE (m)</th></tr></thead><tbody><tr><td rowspan="3">8-May-2015</td><td>1</td><td>N/A</td><td>N/A</td><td>N/A</td><td>N/A</td><td>N/A</td><td>N/A</td></tr><tr><td>2</td><td>-6.78</td><td>0.09</td><td>-2.23</td><td>1.40</td><td>2.63</td><td>0.57</td></tr><tr><td>3</td><td>-3.32</td><td>0.53</td><td>-0.79</td><td>0.63</td><td>1.01</td><td>1.02</td></tr><tr><td rowspan="3">20-May-2015</td><td>1</td><td>-0.21</td><td>3.39</td><td>1.50</td><td>0.81</td><td>1.70</td><td></td></tr><tr><td>2</td><td>-3.10</td><td>3.35</td><td>0.66</td><td>1.06</td><td>1.25</td><td>0.82</td></tr><tr><td>3</td><td>-0.06</td><td>4.45</td><td>2.05</td><td>1.03</td><td>2.29</td><td>0.87</td></tr><tr><td rowspan="3">4-Jun-2015</td><td>1</td><td>-0.57</td><td>0.81</td><td>0.17</td><td>0.20</td><td>0.26</td><td></td></tr><tr><td>2</td><td>-1.86</td><td>0.32</td><td>-0.35</td><td>0.41</td><td>0.54</td><td>0.44</td></tr><tr><td>3</td><td>-3.12</td><td>0.11</td><td>-0.73</td><td>0.54</td><td>0.91</td><td>0.53</td></tr><tr><td rowspan="3">18-June-2015</td><td>1</td><td>-0.46</td><td>1.44</td><td>0.28</td><td>0.27</td><td>0.39</td><td></td></tr><tr><td>2</td><td>-0.54</td><td>2.2</td><td>0.74</td><td>0.63</td><td>0.97</td><td>0.34</td></tr><tr><td>3</td><td>-0.94</td><td>1.00</td><td>0.17</td><td>0.31</td><td>0.35</td><td>0.46</td></tr><tr><td rowspan="3">29-June-2015</td><td>1</td><td>-0.03</td><td>0.77</td><td>0.30</td><td>0.12</td><td>0.32</td><td></td></tr><tr><td>2</td><td>-0.10</td><td>0.91</td><td>0.43</td><td>0.21</td><td>0.48</td><td>0.29</td></tr><tr><td>3</td><td>-0.48</td><td>0.33</td><td>0.02</td><td>0.14</td><td>0.14</td><td>0.30</td></tr><tr><td rowspan="3">4-July-2015</td><td>1</td><td>-0.18</td><td>0.33</td><td>0.03</td><td>0.08</td><td>0.09</td><td></td></tr><tr><td>2</td><td>-0.22</td><td>0.68</td><td>0.27</td><td>0.19</td><td>0.33</td><td>0.29</td></tr><tr><td>3</td><td>-0.96</td><td>0.27</td><td>-0.15</td><td>0.25</td><td>0.29</td><td>0.32</td></tr></tbody></table>

<figure><img src="image_9.png" /><figcaption>**Figure 7.** (Left) scatterplot of the RMSE (root mean square error) of the elevation on the stable areas (snow-free) and the RMSE on the snow surface, (Right) evolution of the RMSE on the snow surface with time.</figcaption></figure>
<figure><img src="image_10.png" /><figcaption>**Figure 8.** Comparison of the RMSE on the surface elevation from Exp.2 and Exp.3 over the snow surface and the stable area.</figcaption></figure>

The results also indicate that the accuracy on stable area elevation is slightly better with Exp. 3 than Exp. 2 (Figure 8). However, this does not mean there was better accuracy for the snow surface elevation, since the RMSE on the snow surface elevation is not systematically higher with Exp. 2.

**5. Discussion**

The fact that the largest biases were found at the beginning of the season is clearly linked to the co-registration method (step 3, Section 3.2). Indeed, for Exp. 2 and Exp. 3, the RMSE on the stable area elevation also decreases when the RMSE on the snow surface elevation decreases (Figure 7). This can be understood by the fact that the snow-free areas, used in step (3), are located near the edges of the camera photographs where the deformation is expected to be maximal (Figure 6). Hence, at the beginning of the snow season, the stable terrain area is too limited to allow a proper alignment of the snow-on and snow-off DEMs. Our workflow provides acceptable results after 20 May 2015, when the snow covered area is not too large to allow a proper alignment between the snow-on and snow-off DEMs. However, despite a smaller snow covered area, there was still considerable snow after that date, since the maximum snow height on 04 Jun-2015, was 5.6 m (Figure A1, Figure A2, Figure A3).

The fact that Exp. 3 provided a better representation of the stable terrain was expected, since the drone images are less oblique than the cameras images. Interestingly, the results suggest that this does not necessarily improve the accuracy on the snow surface elevation, meaning that a drone survey to create a snow-off DEM might not be necessary for another installation. This could be due to the fact that the snow-off camera DEM was generated with photographs having the same acquisition geometry than as the snow-on DEMs, looking toward the main slope of the sinkhole. In this region, the acquisition geometry yields better DEM resolution, due to the steep slope facing the cameras. Since this is the area where most of the snow accumulates, the error on the snow elevation values is not affected by the source of the snow-off DEM.
Several factors could further affect the accuracy of the estimated snow depth, including the lens distortion and camera orientation, which cannot be fully corrected by the structure-from-motion algorithm, given the low number of photographs. The theory of structure-from-motion suggest that a higher accuracy could be achieved using more cameras, but this would increase the cost and the work-load to maintain the system. It is also likely that, more expensive cameras could provide better results especially where digital cameras, with higher dynamic range, may better resolve the texture variations over the snow surface. This would help the image correlation step within the structure-from-motion processing. However, it should be noted that the present accuracy between 20 and 60 cm (standard deviation of the residuals) is already acceptable, given that the observed snow depth ranges between 0 and 8 m in May, and between 0 and 5 m in early June. The error only becomes relatively significant in late June/July when the snow depth ranges between 0 and 2 m (Figures A1–A3 in Appendix A).

**6. Conclusions**

We have studied the accuracy of a new low-cost system to monitor the snow depth in a sinkhole, in the snow-dominated area of Mount-Lebanon. The system is based on three standard compact digital camera photographs, which typically cost USD $100 each. The cost of the system can be further reduced by using a single-board microcontroller instead of the CR200 data-logger, in order to trigger image acquisition.

The results indicate that a simple workflow, using Agisoft Photoscan and Ames Stereo Pipeline software, can estimate the snow surface elevation, with an acceptable accuracy between 20 and 60 cm (standard deviation of the residuals), and a low bias (below 50 cm) when there is sufficient stable terrain around the snow patch. We observed that the co-registration worked well with a ratio of stable area to snow area of roughly one third, but this is probably highly-dependent on the distribution of the stable terrain across the scene. The Agisoft Photoscan software is commercial, and was used in this study because it is efficient and user-friendly, however, open-source alternatives exist to perform the computation of the 3D model with similar accuracy (e.g., MicMac [37]).

The system is autonomous and requires only minimal maintenance. It can be fully automated to generate a time series of snow surface DEMs, as the processing is not based on field data, apart from the coordinates of the cameras that need to be surveyed once. In addition, it was possible to setup a near real-time system by transmitting the cameras images using mobile telecommunications networks. This could be useful in supporting the monitoring of the water resource, since the snow volume in the sinkhole can be considered a proxy of the state of the snow resources in the Mzaar area. Our future work will focus on the generation and analysis of the daily time series of snow volume in the sinkhole over three complete snow seasons from these time-lapse cameras.

**Author Contributions:** Conceptualization, S.G. and J.S.; methodology, C.A.C., S.G., J.S., L.D., and P.F.; software, C.A.C., S.G., and P.F.; validation, C.A.C., S.G., and J.S.; formal analysis, C.A.C., S.G., and J.S.; data curation, C.A.C.; resources (instrumentation), C.A.C., P.F., and L.D., writing—original draft preparation, C.A.C., S.G., and J.S.; writing—review and editing, C.A.C., S.G., and J.S.; supervision, J.S. and L.D.; project administration, J.S. and L.D.; funding acquisition, J.S. and L.D.

**Funding:** This research was funded by the Programme National de Télédétection Spatiale (PNTS, http://www.insu.cnrs.fr/pnts), grant no. PNTS-2018-4.

**Acknowledgments:** This work has been supported by the Research Council at USJ and the United Nations Children's Fund (UNICEF). It was carried out in the frame of the ERANET-MED03-62 CHAAMS project and the Observatoire Franco-Libanais pour l'Environnement (O-LIFE). We warmly thank the two anonymous reviewers for their comments and suggestions, which helped improve the paper.

**Conflicts of Interest:** The authors declare no conflict of interest.
Appendix A

<figure><img src="image_11.png" /><figcaption>Figure A1. Snow depth maps of Exp.1.</figcaption></figure>
<figure><img src="image_12.png" /><figcaption>Figure A2. Snow depth maps of Exp.2.</figcaption></figure>
<figure><img src="image_13.png" /><figcaption>Figure A3. Snow depth maps of Exp.3.</figcaption></figure>
## References

1. Baylouny, A.M.; Klingseis, S.J. Water thieves or Political catalysts? syrian refugees in Jordan and Lebanon. *Middle east Policy* **2018**, *25*, 104–123. [CrossRef]

2. Shaban, A.; Faour, G.; Khawlie, M.; Abdallah, C. Remote sensing application to estimate the volume of water in the form of snow on Mount Lebanon/Application de la télédétection à l'estimation du volume d'eau sous forme de neige sur le Mont Liban. *Hydrol. Sci. J.* **2004**, *49*, 37–41. [CrossRef]

3. Shaban, A.; Dariwich, T. Investigating snow characteristics on mountain chains of Lebanon. In Proceedings of the First CIGR Inter-regional Conference on Land and Water Challenges, Bari, Italy, 10–14 September 2013.

4. Telesca, L.; Shaban, A.; Gascoin, S.; Darwich, T.; Drapeau, L.; El Hage, M.; Faour, G. Characterization of the time dynamics of monthly satellite snow cover data on Mountain Chains in Lebanon. *J. Hydrol.* **2014**, *519*, 3214–3222. [CrossRef]

5. Mhawej, M.; Faour, G.; Fayad, A.; Shaban, A. Towards an enhanced method to map snow cover areas and derive snow-water equivalent in Lebanon. *J. Hydrol.* **2014**, *513*, 274–282. [CrossRef]

6. Somma, J.; Drapeau, L.; Abou Chakra, C. Caractérisation de la capacité de rétention de neige d'un plateau karstique. Essai méthodologique pour le plateau du Jabal Jraid (Liban). *Revue internationale de géomatique* **2014**, *24*, 87–99. [CrossRef]

7. Fayad, A.; Gascoin, S.; Faour, G.; Fanise, P.; Drapeau, L.; Somma, J.; Fadel, A.; Al Bitar, A.; Escadafal, R. Snow observations in Mount Lebanon (2011–2016). *Earth Syst. Sci. Data* **2017**, *9*, 573–587. [CrossRef]

8. Dubertret, L. Carte Géologique du Liban. Available online: https://www.worldcat.org/title/carte-geologique-du-liban/oclc/21747724 (accessed on 9 September 2019).

9. Edgell, H.S. Karst and hydrogeology of Lebanon. *Carbonates Evaporites* **1997**, *12*, 220. [CrossRef]

10. Fayad, A.; Gascoin, S.; Faour, G.; López-Moreno, J.I.; Drapeau, L.; Le Page, M.; Escadafal, R. Snow hydrology in Mediterranean mountain regions: A review. *J. Hydrol.* **2017**, *551*, 374–396. [CrossRef]

11. Jonas, T.; Marty, C.; Magnusson, J. Estimating the snow water equivalent from snow depth measurements in the Swiss Alps. *J. Hydrol.* **2009**, *378*, 161–167. [CrossRef]

12. López-Moreno, J.I.; Fassnacht, S.R.; Heath, J.T.; Musselman, K.N.; Revuelto, J.; Latron, J.; Morán-Tejeda, E.; Jonas, T. Small scale spatial variability of snow density and depth over complex alpine terrain: Implications for estimating snow water equivalent. *Adv. Water Resour.* **2013**, *55*, 40–52. [CrossRef]

13. Sturm, M.; Taras, B.; Liston, G.E.; Derksen, C.; Jonas, T.; Lea, J.; Sturm, M.; Taras, B.; Liston, G.E.; Derksen, C.; et al. Estimating Snow Water Equivalent Using Snow Depth Data and Climate Classes. *J. Hydrol.* **2010**, *11*, 1380–1394. [CrossRef]

14. Bühler, Y.; Adams, M.S.; Bösch, R.; Stoffel, A. Mapping snow depth in alpine terrain with unmanned aerial systems (UASs): Potential and limitations. *Cryosphere* **2016**, *10*, 1075–1088. [CrossRef]

15. Marti, R.; Gascoin, S.; Berthier, E.; De Pinel, M.; Houet, T.; Laffly, D. Mapping snow depth in open alpine terrain from stereo satellite imagery. *Cryosphere* **2016**, *10*, 1361–1380. [CrossRef]

16. Deems, J.S.; Painter, T.H.; Finnegan, D.C. Lidar measurement of snow depth: A review. *J. Glaciol.* **2013**, *59*, 467–479. [CrossRef]

17. Garvelmann, J.; Pohl, S.; Weiler, M. From observation to the quantification of snow processes with a time-lapse camera network. *Hydrol. Earth Syst. Sci.* **2013**, *17*, 1415–1429. [CrossRef]

18. Lui, J.; Chen, R.; Wang, G. Snowline and snow cover monitoring at high spatial resolution in a mountainous river basin based on a time-lapse camera at a daily scale. *J. Mt Sci.* **2015**, *12*, 60–69.

19. Kępski, D.; Luks, B.; Migała, K.; Wawrzyniak, T.; Westermann, S.; Wojtuń, B. Terrestrial Remote Sensing of Snowmelt in a Diverse High-Arctic Tundra Environment Using Time-Lapse Imagery. *Remote Sens.* **2017**, *9*, 733. [CrossRef]

20. Parajka, J.; Haas, P.; Kirnbauer, R.; Jansa, J.; Blöschl, G. Potential of time-lapse photography of snow for hydrological purposes at the small catchment scale. *Hydrol. Processes* **2012**, *26*, 3327–3337. [CrossRef]

21. Revuelto, J.; Jonas, T.; López-Moreno, J.I. Backward snow depth reconstruction at high spatial resolution based on time-lapse photography. *Hydrol. Processes* **2016**, *30*, 2976–2990. [CrossRef]

22. Westoby, M.J.; Brasington, J.; Glasser, N.F.; Hambrey, M.J.; Reynolds, J.M. 'Structure-from-Motion' photogrammetry: A low-cost, effective tool for geoscience applications. *Geomorphology* **2012**, *179*, 300–314. [CrossRef]
23. Fonstad, M.A.; Dietrich, J.T.; Courville, B.C.; Jensen, J.L.; Carbonneau, P.E. Topographic structure from motion: A new development in photogrammetric measurement. *Earth Surf. Processes Landforms* **2013**, *38*, 421–430. [CrossRef]

24. Micheletti, N.; Chandler, J.H.; Lane, S.N. Structure from Motion (SfM) Photogrammetry. In *Geomorphological Techniques (Online Edition)*; Clarke, L.E., Nield, J.M., Eds.; British Society for Geomorphology: London, UK, 2015; Chap. 2.

25. Morgan, J.A.; Brogan, D.J.; Nelson, P.A. Application of Structure-from-Motion photogrammetry in laboratory flumes. *Geomorphology* **2017**, *276*, 125–143. [CrossRef]

26. Balaguer-Puig, M.; Marqués-Mateu, Á.; Lerma, J.L.; Ibáñez-Asensio, S. Estimation of small-scale soil erosion in laboratory experiments with Structure from Motion photogrammetry. *Geomorphology* **2017**, *295*, 285–296. [CrossRef]

27. Mertes, J.R.; Gulley, J.D.; Benn, D.I.; Thompson, S.S.; Nicholson, L.I. Using structure-from-motion to create glacier DEMs and orthoimagery from historical terrestrial and oblique aerial imagery. *Earth Surf. Processes Landforms* **2017**, *42*, 2350–2364. [CrossRef]

28. Cimoli, E. Determining Snow Depth Distribution from Unmanned Aerial Vehicles and Digital Photogrammetry. Master's Thesis, Technical University of Denmark, Kongens Lyngby, Denmark, August 2015.

29. Brunier, G.; Fleury, J.; Anthony, E.J.; Gardel, A.; Dussouillez, P. Close-range airborne Structure-from-Motion Photogrammetry for high-resolution beach morphometric surveys: Examples from an embayed rotating beach. *Geomorphology* **2016**, *261*, 76–88. [CrossRef]

30. Gindraux, S.; Boesch, R.; Farinotti, D. Accuracy assessment of digital surface models from Unmanned Aerial Vehicles' imagery on glaciers. *Remote Sens.* **2017**, *9*, 186. [CrossRef]

31. Nissen, E.; Krishnan, A.K.; Arrowsmith, J.R.; Saripalli, S. Three-dimensional surface displacements and rotations from differencing pre-and post-earthquake LiDAR point clouds. *Geophys. Res. Lett.* **2012**, *39*, 1–6. [CrossRef]

32. He, Y.; Liang, B.; Yang, J.; Li, S.; He, J. An Iterative Closest Points Algorithm for Registration of 3D Laser Scanner Point Clouds with Geometric Features. *Sensors* **2017**, *17*, 1862. [CrossRef]

33. Teza, G.; Galgaro, A.; Zaltron, N.; Genevois, R. Terrestrial laser scanner to detect landslide displacement fields: A new approach. *Int. J. Remote Sens.* **2007**, *28*, 3425–3446. [CrossRef]

34. Shean, D.; Alexandrov, O.; Moratto, Z.; Smith, B.; Joughin, I. An automated, open-source pipeline for mass production of digital elevation models (DEMs) from very-high-resolution commercial stereo satellite imagery. *ISPRS J. Photogramm. Remote Sens.* **2016**, *116*, 101–117. [CrossRef]

35. Lavine, A.; Gardner, J.N.; Reneau, S.L. Total station geologic mapping: An innovative approach to analyzing surface-faulting hazards. *Eng. Geol.* **2003**, *70*, 71–91. [CrossRef]

36. Wolf, D.G.P.R. *Elementary Surveying*; Pearson: Essex, UK, 2014.

37. Jaud, M.; Passot, S.; Le Bivic, R.; Delacourt, C.; Grandjean, P.; Le Dantec, N. Assessing the accuracy of high resolution digital surface models computed by PhotoScan® and MicMac® in sub-optimal survey conditions. *Remote Sens.* **2016**, *8*, 465. [CrossRef]

<figure><img src="image_14.png" /></figure>