Daily Monitoring of Shallow and Fine-Grained Water Patterns in Wet Grasslands Combining Aerial LiDAR Data and In Situ Piezometric Measurements

: The real-time monitoring of hydrodynamics in wetlands at ﬁne spatial and temporal scales is crucial for understanding ecological and hydrological processes. The key interest of light detection and ranging (LiDAR) data is its ability to accurately detect microtopography. However, how such data may account for subtle wetland ﬂooding changes in both space and time still needs to be tested, even though the degree to which these changes impact biodiversity patterns is of upmost importance. This study assesses the use of 1 m × 1 m resolution aerial LiDAR data in combination with in situ piezometric measurements in order to predict the ﬂooded areas at a daily scale along a one-year hydrological period. The simulation was applied over 663 ha of wet grasslands distributed on six sites across the Marais Poitevin (France). A set of seven remote sensing images was used as the reference data in order to validate the simulation and provide a high overall accuracy (76–94%). The best results were observed in areas where the ditch density was low, whereas the highly drained sites showed a discrepancy with the predicted ﬂooded areas. The landscape proportion index was calculated for the daily steps. The results highlighted the spatiotemporal dynamics of the shallow ﬂooded areas. We showed that the differences in the ﬂooding durations among the years were mainly related to a narrow contrast in topography (40 cm), and occurred over a short period of time (two months).


Introduction
In wetlands, hydrological processes affect biogeochemical and ecological functions together with hydrological hazards [1]. As an example, runoff from rivers to wetlands prevent flood risks in downstream urbanized areas [2]; hydrological pulses transport phosphorus into the floodplain drainage network [3]; the flood duration drives plant diversity and functional traits in wet grasslands [4]; and the water level in grasslands affects bird diversity [5]. However, the flooding pattern and its dynamics along space and time are still poorly understood, and as a consequence they are also poorly predicted [6]. Some of the main challenges include producing: (1) reliable predictions while taking the subtle spatial heterogeneity of marshlands into account, and (2) predictions that may fit the large variations that are observed for the flooding areas over a short period of time. In order to grasslands. The advances as well as the limitations of this model for the understanding of ecohydrology processes are discussed.

Study Area
The study area comprises six sites for a total of 663 hectares of wet grassland distributed across the Marais Poitevin, a marshland area near the French Atlantic coast (46.4°N, 1.2°W) (Figure 1). The climate is oceanic temperate, with a mean monthly minimum/maximum temperature ranging from 2/10 °C in winter to 12/24 °C in summer. The mean annual precipitation ranges from 700 mm to 900 mm, with a summer water deficit. Deriving from salt marshes after embankment, the studied marshland shows an elevation gradient of limited extent, with depressions that are 50-80 cm deep, and upper flats. The elevation ranges between 1.4-3.1 m above sea level. A high-density ditch network efficiently drained the marsh for agriculture activities (Table 1). However, the grasslands may be flooded from winter to spring, when the soils are waterlogged and there is significant precipitation. Therefore, the study area is affected with biodiversity [32], and is part of the NATURA 2000 European network [33]. As a consequence, management of the water level is subjected to a trade-off between economic and ecological issues.

Data Collection and Pre-Processing
A 1 m × 1 m resolution digital terrain model (DTM)-the Litto3D product-was delivered by the French Mapping Agency (IGN, Saint Mandé, France). It was derived from airborne LiDAR (light detection and ranging) data acquired in the summer of 2012 with a 2 pts/m 2 point density. In grasslands, the altimetric accuracy of the LiDAR data usually ranges from 5 cm to 30 cm, depending on the vegetation height and density cover [15]. In order to assess the altimetric accuracy of the LiDAR-based DTM, differential global positioning system (DGPS) plots (n = 395) with a 2-cm altimetry accuracy were recorded across the six sites along topographic transects. The DTM values are expressed relative to sea level in the French national georeferencing system (NGF-93). In order to reduce the raster size, the DTMs were clipped and analyzed according to the boundaries of the sites.
A piezometric probe was installed in October 2014 in the lower part of each site. The water level elevations were recorded hourly with an accuracy of 2 cm. The probe values are also expressed in NGF-93. The database for each probe was collected during the field visits (the last in March 2016), and the hourly water levels were averaged daily.
A set of seven earth observation images, including two aerial photographs and five SPOT-6/7 multispectral images, was collected in order to validate the simulated flood maps, and used as the reference data ( Table 2). The two aerial images acquired by the microlight aircraft during the winter flooding periods were used as the reference data, with a first acquisition campaign on 4 March 2015 for sites 2-6, and a second campaign on 15 January 2016 for site 1. The aerial images were orthorectified and georeferenced using the spline transformation from ArcGIS (ESRI, Redlands, CA, USA). Approximately 50 to 60 reference points were chosen on each site based on cross ditches that were easily identifiable on both the aerial image and the LiDAR-based DTM. The SPOT-6/7 images were delivered free of charge from the Geosud archive catalog (ids.equipex-geosud.fr) as an orthorectified level product. A total of five cloud-free SPOT-6/7 images were found to be usable for all or part of the area of the six sites. The SPOT-6/7 sensors featured panchromatic (450-745 nm), blue (450-525 nm), green (530-590 nm), red (625-695 nm), and near-infrared (760-890 nm) spectral bands, and had a spatial resolution of 6.0 m and 1.5 m in multispectral and panchromatic mode, respectively. In order to improve the spatial resolution of the multispectral bands to 1.5 m, a pan-sharpening process was applied by Gram-Schmidt fusion [34]. Each site was covered by at least three remote sensing images.

Model Calibration and Validation
To model the flood extent at a daily step interval, the LiDAR-based DTM dataset was used and transformed into a binary raster (0 = unflooded; 1 = flooded) that was specific to each site, according to the daily water level measured by the piezometric probe. The model was iteratively performed using the R software [27] with the Raster package [28]. Since the accuracy of the LiDAR-based DTM altimetry values range between 5-30 cm in grasslands [29], a site-specific DTM calibration was necessary. For each DGPS plot, the elevation value was compared with the corresponding DTM pixel value. Within each site, the median elevation difference value measured between the DGPS plots and the DTM was assigned as the offset value. We assumed that the altimetry offset is constant across a site.
Consequently, the same offset value was added to the DTM pixel value, and resulted in site-specific calibrated DTMs.
In order to validate the model, we performed a cross-validation for each site between the flooding areas, which was measured from the remote sensing reference images with the areas simulated in the LiDAR-based DTM. The flood extents observed on each remote sensing image were semi-automatically extracted by an object-based approach. Specifically, a multiresolution segmentation was applied to merge each adjacent pixel into objects with similar values [35] using the eCognition software (Definiens Imaging, Munich, Germany). As the spectral response of the flooded areas showed large variability depending on the sensor type (pan-sharpened multispectral imagery, microlight aircraft) and environmental conditions (water depth, vegetation height, wind force, etc.), the objects were manually sorted into unflooded or flooded categories. Considering that the studied sites were covered by herbaceous vegetation, we assumed that a minimum water depth of 15 cm was required in order to be detectable on the aerial image. Simulated flooded pixels with a water depth of less than 15 cm were thereafter reclassified as unflooded for the validation.
The comparison between the observed and simulated flood extents produced overall accuracy (OA) and Kappa location (K loc ) indexes. These analyses were carried out using the Idrisi GIS software (Clark labs, Worcester, MA, USA). The popular OA index provides information on the proportion of pixels that are correctly simulated. The K loc index is a variation of the standard Kappa index of agreement, but takes the spatial location of the errors into account [36]. It was calculated according to the following Equation (1): where Po is the observed proportion of cells that are correctly classified, Pc is the expected proportion that is correct due to chance, and Pp is the proportion that is correct when the classification is perfect. Pc was expressed as follows in Equation (2): where Pc is the expected proportion that is correct due to chance, S j is the proportion of cells in category j in the simulation, and R j is the proportion of cells in category j in the reference.
Pp was expressed as follows in Equation (3): where Pp is the proportion that is correct when the classification is perfect, S j is the proportion of cells in category j in the simulation, and R j is the proportion of cells in category j in the reference. Regarding the K loc values, we considered the simulation performances as poor when <0.4; as moderate for values between 0.41-0.60; and as good or excellent when respectively between 0.61-0.80 and >0.81.

Spatiotemporal Monitoring and Characterization of the Flood Patterns
The flood extents were modeled daily for each site during one hydrological year, from 1 November 2014 to 31 October 2015. Each daily simulated map of the flood extent was iteratively added in order to produce the annual flood duration map. Transitional maps representing the cumulative flood duration were produced at a bi-monthly interval. All of the output maps were modeled with a 1 m × 1 m resolution.
In order to quantify the temporal changes regarding the flooded area, the "landscape proportion index" (LPI) was calculated daily and for each site, from the simulated flood maps using the SDMTools R-package [32]. The landscape analysis was performed separately on both the flooded and unflooded categories. The LPI describes the proportional abundance of each patch category in each site. The LPI is calculated as follows in Equation (4): where LPI is the landscape proportion index; a ij is the area of patch ij in m 2 ; and A is the total area in m 2 in each site. In order to compare the hydrological performance between the sites, a hierarchical clustering of the LPI time series was applied using a shape-based distance [37] that was implemented in the dtwclust R-package [38].

Accuracy Assessment
The median LiDAR-based DTM altimetry errors (i.e., offset) ranged from 6 cm to 26 cm, depending on the site ( Figure 2). The minimum and maximum altimetry errors were found on sites 5 and 1, respectively. For sites 2, 3, 4, and 6, the values for the median errors were similar, and ranged from 16 cm to 19 cm. The standard deviation errors ranged from 6.5 (site 5) to 14.6 (site 6) cm, and revealed a moderate error variability within the sites. For site 1, the higher altimetry errors could be explained by the vegetation cover (tall and dense grasslands,~80 cm in height). The offset (i.e., median error) values applied for the calibration of the LiDAR-based DTMs are listed in Table 2. and unflooded categories. The LPI describes the proportional abundance of each patch category in each site. The LPI is calculated as follows in Equation (4): where LPI is the landscape proportion index; aij is the area of patch ij in m²; and A is the total area in m² in each site. In order to compare the hydrological performance between the sites, a hierarchical clustering of the LPI time series was applied using a shape-based distance [37] that was implemented in the dtwclust R-package [38].

Accuracy Assessment
The median LiDAR-based DTM altimetry errors (i.e., offset) ranged from 6 cm to 26 cm, depending on the site ( Figure 2). The minimum and maximum altimetry errors were found on sites 5 and 1, respectively. For sites 2, 3, 4, and 6, the values for the median errors were similar, and ranged from 16 cm to 19 cm. The standard deviation errors ranged from 6.5 (site 5) to 14.6 (site 6) cm, and revealed a moderate error variability within the sites. For site 1, the higher altimetry errors could be explained by the vegetation cover (tall and dense grasslands, ~80 cm in height). The offset (i.e., median error) values applied for the calibration of the LiDAR-based DTMs are listed in Table 2. Statistical comparisons between the maps of the flooded areas observed from the seven remote sensing images and the pattern simulated from the LiDAR data are shown in Table 3. In general, the mean OA values ranged from 82.3% to 91.5%, and thus show that there is a good agreement between the simulations and the remotely-sensed observations, regardless of the site ( Table 3). The analysis Statistical comparisons between the maps of the flooded areas observed from the seven remote sensing images and the pattern simulated from the LiDAR data are shown in Table 3. In general, the mean OA values ranged from 82.3% to 91.5%, and thus show that there is a good agreement between the simulations and the remotely-sensed observations, regardless of the site ( Table 3). The analysis of the mean K loc index also highlights the differences among the sites in terms of the error location. Mean K loc values higher than 0.6 were found on sites 3 and 4-revealing a spatial concentration for the error location-whereas K loc values lower than 0.6 were measured on sites 1, 2, 5, and 6-expressing a high spatial dispersion of the error. More specifically, a site-per-site interpretation highlights variations in the flood simulation performance among the seasons. On site 1, the OA values were excellent in both June (88.9-92.2%) and January (84.7%), and the K loc index ranged from poor/moderate in June  As a general trend, the observed and simulated flood extents showed a high level of similarity, in accordance with the OA and K loc values (Figure 3). More specifically, two situations were recorded. Firstly, during the period when grasslands were dry and water was only present in ditches and permanent ponds on the reference maps (e.g., 28 June 2015 for sites 1-3 or 12 April 2015 for site 5), the flooded simulated extent was nearly null in the grasslands, as well as in the ditches and permanent ponds. These errors can be explained by the presence of water in the ditches and permanent ponds at the time of the LiDAR acquisition, which is an obstacle to LiDAR signal penetration, and therefore to microtopographic characterization. Secondly, when the sites were widely flooded (e.g., 15 January 2016 for site 1, or 14 March 2016 for sites [4][5], then, the fine-scale spatial patterns of the flooded areas were found to be reliably simulated, regardless of whether the flooding occurred in large depressions or in narrow ditches. In this situation, the simulation performs well, even several hundred meters away from the piezometric probe. However, some underdetection errors could be noticed (e.g., 4 March 2015 on site 2, and in the eastern part of site 3), which was probably related to the altimetric error of the LiDAR DTM in relation to the vegetation cover (tall and dense grasslands,~80 cm in height).

Spatiotemporal Characterization of the Flood Patterns
The LPI analyses showed that the flooding events generally occurred between November 2014 and June 2015 (Figure 4). Within this general pattern, a fine temporal variability of the LPI can be noticed, particularly on sites 2-5. The highest LPI values can be seen on sites 4-6. The flood extents

Spatiotemporal Characterization of the Flood Patterns
The LPI analyses showed that the flooding events generally occurred between November 2014 and June 2015 (Figure 4). Within this general pattern, a fine temporal variability of the LPI can be noticed, particularly on sites 2-5. The highest LPI values can be seen on sites 4-6. The flood extents became restricted in summer and autumn, except on site 1, where the LPI increased again from September onward.  The comparison of the hydrological behavior between the sites is presented in Figure 5. Based on the visual inspection of the dendrogram, we chose to set the threshold distance value to 0.15. Sites 1, 2, and 5 are classified into three different clusters, while sites 3, 4, and 6 are grouped within a common cluster.  The comparison of the hydrological behavior between the sites is presented in Figure 5. Based on the visual inspection of the dendrogram, we chose to set the threshold distance value to 0.15. Sites 1, 2, and 5 are classified into three different clusters, while sites 3, 4, and 6 are grouped within a common cluster.  The comparison of the hydrological behavior between the sites is presented in Figure 5. Based on the visual inspection of the dendrogram, we chose to set the threshold distance value to 0.15. Sites 1, 2, and 5 are classified into three different clusters, while sites 3, 4, and 6 are grouped within a common cluster. The map of the annual flooded duration was produced for all of the sites (Figure 6). Even within the limited altitudinal gradient (<1 m), large variations in the annual flood duration were found, ranging from 0 for the upper locations, and up to 365 days for the lower areas. On sites 1-3, the The map of the annual flooded duration was produced for all of the sites ( Figure 6). Even within the limited altitudinal gradient (<1 m), large variations in the annual flood duration were found, ranging from 0 for the upper locations, and up to 365 days for the lower areas. On sites 1-3, the longest flooded areas (>180 days a year) were located in paleochannels with pronounced depressions (>50 cm). However, the situation was quite different on sites 4-6, where the longest flooded areas occurred and where the water was lengthily retained by managements that may not only have been in the deeper depressions. Overall, the annual flood durations were lower on site 2 (60-120 days), and especially on site 5 (0-30 days), where a dense ditch network may have efficiently drained the water.  In order to illustrate the fine-grained and sinuous shape of the flooded areas in the studied coastal marshland, Figure 7 focused on site 4, which shows the typical patterns in the region. Each pixel is characterized by its temporal flooding profile. As an example, three temporal profiles were extracted from three plots (i.e., pixels) chosen along a fine-i.e., 40-cm large-elevation gradient. These results showed that the differences in the flood duration among the sites mainly occurred during late spring and early summer. In addition, a high variability in the flood duration can be observed between plots even only a few meters away. Plot 1, which was situated on a higher elevation, was flooded for 92 days from January to April; whereas plot 2, which was located at an intermediate elevation, was flooded for 143 days, from January to June; and Plot 3, located at a lower elevation, was flooded for 216 days, from November to July. In order to illustrate the fine-grained and sinuous shape of the flooded areas in the studied coastal marshland, Figure 7 focused on site 4, which shows the typical patterns in the region. Each pixel is characterized by its temporal flooding profile. As an example, three temporal profiles were extracted from three plots (i.e., pixels) chosen along a fine-i.e., 40-cm large-elevation gradient. These results showed that the differences in the flood duration among the sites mainly occurred during late spring and early summer. In addition, a high variability in the flood duration can be observed between plots even only a few meters away. Plot 1, which was situated on a higher elevation, was flooded for 92 days from January to April; whereas plot 2, which was located at an intermediate elevation, was flooded for 143 days, from January to June; and Plot 3, located at a lower elevation, was flooded for 216 days, from November to July.

Model Performance and Limitations
The calibration of the predictive model of the flooding extent was based on a constant offset value that was related to the median altimetry error of the LiDAR-based DTM. The comparison between the DGPS measurements and DTM pixels showed strong similarities in terms of the elevation values (a variation of 10 centimeters in the altimetry errors). Although some differences higher than 10 cm were observed in the dense herbaceous vegetation and shallow water bodies, this finding corroborates with previous studies [15,28]. To minimize the altimetric error of the LiDAR-based DTM, we tested the effect of applying a plant species-specific offset; this correction by plant height was shown to be efficient, although it may not have fully balanced the errors: the post-correction accuracy assessments reached a mean altimetric error of ~20cm in the salt marshes [39]. This remaining altimetric error may partially explain part of the errors in the flooding simulations we produced, while no statistically significant relationship was found between the LiDAR altimetry error and the accuracy of the simulated flooding areas. As an example, sites 1 and 5 respectively had high (26 cm) and low (6 cm) median LiDAR altimetry errors, but similar OA values (88.8% and 87.6% respectively). Variations in the Kloc values appeared to be strongly correlated to the ditch density. On sites 2, 3, and 5, where the Kloc values were the lowest in the winter period, the ditch density was the highest. On the contrary, sites 1, 4, and 6, where the Kloc values were the highest in the winter period, showed the lowest density of ditches. These findings confirm previous studies that highlighted the ditch network density as an important control variable in hydrological processes in wet grasslands [1,26].
The validation of the model was based on seven remote sensing images acquired between late winter and early summer. This is a critical period in temperate wetlands in terms of the large

Model Performance and Limitations
The calibration of the predictive model of the flooding extent was based on a constant offset value that was related to the median altimetry error of the LiDAR-based DTM. The comparison between the DGPS measurements and DTM pixels showed strong similarities in terms of the elevation values (a variation of 10 centimeters in the altimetry errors). Although some differences higher than 10 cm were observed in the dense herbaceous vegetation and shallow water bodies, this finding corroborates with previous studies [15,28]. To minimize the altimetric error of the LiDAR-based DTM, we tested the effect of applying a plant species-specific offset; this correction by plant height was shown to be efficient, although it may not have fully balanced the errors: the post-correction accuracy assessments reached a mean altimetric error of~20cm in the salt marshes [39]. This remaining altimetric error may partially explain part of the errors in the flooding simulations we produced, while no statistically significant relationship was found between the LiDAR altimetry error and the accuracy of the simulated flooding areas. As an example, sites 1 and 5 respectively had high (26 cm) and low (6 cm) median LiDAR altimetry errors, but similar OA values (88.8% and 87.6% respectively). Variations in the K loc values appeared to be strongly correlated to the ditch density. On sites 2, 3, and 5, where the K loc values were the lowest in the winter period, the ditch density was the highest. On the contrary, sites 1, 4, and 6, where the K loc values were the highest in the winter period, showed the lowest density of ditches. These findings confirm previous studies that highlighted the ditch network density as an important control variable in hydrological processes in wet grasslands [1,26].
The validation of the model was based on seven remote sensing images acquired between late winter and early summer. This is a critical period in temperate wetlands in terms of the large variations occurring in the flood extent [10] and plant growth dynamics [4]. The results showed the good quality of the predictions made by the model we used, as the OA values ranged between 76-94%. Compared to the in situ validation, this approach presents some pros and cons. The first limitation concerns the use of three-date reference data per site: as a consequence, we assumed that this relationship between the piezometric altimetry value and flooded areas is globally constant throughout the year. With respect to the complex hydrological processes occurring in both the surface water and groundwater [1,40], this assumption could be considered simplistic. Additional data, spread out throughout the year, are needed in order to better calibrate the model and take the variation in the piezometric altimetry values and the preduction of flooded areas into account. Recently, Keys et al. [41] developed a cost-effective approach that was based on the use of a camera placed on top of a 10-m high field tower, which records hourly time-step images. Another large-scale efficient alternative could be the use of a new SAR time series. Recently, Cazals et al. [10] showed that flooded areas may be readily detectable in the Marais Poitevin (100,000 ha) using a 10-day revisit frequency for the Sentinel-1 satellite time series, at the regional scale. However, this work highlighted the limits of the radar data: flooded areas lower than 0.8 ha and/or shallow water with emergent vegetation were poorly detected, whereas they are actually very frequent in wetlands. On the contrary, the LiDAR-based simulations were able to successfully produce a daily fine-grained pattern characterization, but were constrained on sites measuring only a few hectares. In a nutshell, the observed and simulated approaches appeared to be complementary, since they were suitable for regional/monthly and local/daily scales, respectively.
The second limitation lies in the possible visual misinterpretations of the flood extents of the remotely sensed data that were used to establish the reference flooding maps. Hence, we assumed that the flooded areas were detected with a minimum water depth of 15 cm, which is relevant based on our field knowledge; however, this water depth threshold may locally underdetect or overdetect the flooding areas that were covered by vegetation taller than 30 cm or smaller than 10 cm.
Conversely, the main advantage of the approach proposed here is its significant spatial sampling based on a raster reference map, including hundreds of pixels (i.e., validation samples). The reference map offers the opportunity to calculate a K loc index, which assesses the spatial dispersion of the error, and is well adapted (better than the Kappa index) for the assessment of the simulated map [36]. Another interest of the method used here is that the errors on the simulated maps can be spatialized: the comparison between the reference maps showed that the simulated maps were correctly modeled, and not only close to the piezometric probe, but also hundreds of meters away.
In contrast to hydrological models that require a high level of expertise and intensive data acquisition [40], the model we developed, which was based on LiDAR and piezometric data, is simple and suitable for mapping both surface and subsurface water with spatially fine-grained precision across many hectares. Our results underlined how microtopography data may be useful and help predict flooding in wetlands [22,42].
It should also be noticed that while many studies have pointed out the importance of a hydrological network [19], vegetation resistance [18,20,21,42], and infiltration [23] in flow modeling, these factors were not considered here, and could explain some errors in our simulation. However, the integration of these factors in hydrological modeling remains challenging for three reasons. Firstly, we acknowledge that the hydrological connectivity between wet grasslands and the hydrological network is crucial for assessing wetland functions [19], but this connectivity is difficult, perhaps even impossible, to be reliably assessed since floodplains are connected through small intermittent creeks or waterflows that are not captured by LiDAR data [25]. Recently, McDonough et al. [43] suggested monitoring the surface hydrologic connectivity between wetlands and streams with a float switch data logger system. The coupling of this system with the use of both LiDAR and piezometric data, as done in this study, would provide a great way to monitor the effective connection/disconnection pattern in time and space. Alternatively, this limitation could be circumvented by using new bathymetric LiDAR data, which are suitable for capturing channel beds [44]. Secondly, the vegetation resistance is highly variable in time in relation to the growing season, and the somewhat unreliable inputs from these parameters could be another source of errors, rather than a significant improvement [42]. Thirdly, the infiltration parameter is difficult to spatialize at the fine scale with sufficient certainty. The integration of remotely-sensed SAR (such as TerraSAR-X) time series in flood inundation models could provide promising avenues for mapping flood hazards and overland flows [10,13].

Toward a Better Knowledge of Ecohydrological Processes
The daily simulation of fine-grained flood patterns in wet grasslands offers a relevant technical and methodological solution for a real-time ecohydrology [45]. For example, this approach is able to detect rapid changes in the characteristics of the signal "pulse" based on the LPI value, as in site 6 in May (Figure 4). Considering the flood duration map (Figure 7), areas with significant changes compared to their surroundings can be identified, and considered "hot spots" [45]. Some "thresholds" can be highlighted in May, when small changes in the water management led to an abrupt decrease in the extent of the flooded areas ( Figure 4). Such abrupt variations in spring, when reproduction is occurring for most of the flora and fauna, may present strong consequences for the biodiversity patterns. Interestingly, the hierarchical clustering of the LPI time series highlighted a common cluster including all three sites concerned by conservation planning for waterbirds, which are generally accompanied by a longer flooding period in the spring ( Figure 5).
In this study, the water levels measured by the piezometric probes result from both natural (precipitation and evapotranspiration) and anthropogenic drivers (water management). Although the daily-step flooded maps provide interesting information for local managers, the real impacts of water management on flooded areas remain unclear with regard to climatic conditions. The recent development of an interactive simulation tool to promote adaptation to climate variability [6] constitutes an interesting track to follow in order to discriminate the respective contributions of natural and anthropogenic variables on the water level dynamics.
There will be many applications for the inundation map that has been produced, in particular in the fields of plant and animal ecology. Regarding plant biodiversity and communities, it will provide needed information for predicting the distribution of plant species and communities [46], as well as plant functional traits and species richness, as shown by Violle et al. [4]. Our simulated map of flooding at a daily step could explore this relationship with a larger environmental characterization of the flood patterns and a more in-depth characterization of the key time period and interannual variations. Within each site, the differences in the annual flood duration were mainly related to contrasts occurring in narrow intervals both in time (two months) and topography (40 cm).
The modeling of the daily flooded areas over a one-year period, over 661 ha, with a 1 m × 1 m spatial resolution produces a huge amount of data (>7 GB). It could be interesting to carry out advanced analyses, such as the classification of pixels based on their temporal profiles. However, the spatial and temporal analyses of these large datasets are still challenging, and require further algorithmic and computational developments [47].

Conclusions
In this study, we modeled the flooding pattern in wet and highly drained grasslands with both fine temporal precision (one day) and on a fine spatial scale (one m). The model was built considering both LiDAR and in situ piezometric data. The comparison with the remote sensing-based reference maps showed that the model produces a good quality prediction. Better results were observed on the sites with a low ditch density (54-87 m/ha); whereas the highly drained sites (163-233 m/ha) showed a discrepancy with the model predictions. Further modeling developments that couple LiDAR data and drainage networks will be needed in order to improve the accuracy of the flood simulation in the high drained wet grasslands. The maps of the model outputs provide fine-grained maps showing the flooding and drying-out pattern of the shallow water patches. On this basis, the relationships between biodiversity and flooding dynamics may be studied, together with their respective spatial patterns.