Next Article in Journal
Transferable Principles for Managing the Nexus: Lessons from Historical Global Water Modelling of Central Asia
Previous Article in Journal
Comparison of the Physical Properties of Showers that the Satisfaction of Shower Feeling among Users in Three Asian Countries
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Reconstruction of Lake Level Changes of Groundwater-Fed Lakes in Northeastern Germany Using RapidEye Time Series

1
Helmholtz Centre Potsdam, GFZ German Research Centre for Geosciences, Telegrafenberg, Potsdam 14473, Germany
2
Geoinformation in Environmental Planning Lab, Technische Universität Berlin, Straße des 17. Juni 145, Berlin 10623, Germany
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Water 2015, 7(8), 4175-4199; https://doi.org/10.3390/w7084175
Submission received: 2 April 2015 / Revised: 21 July 2015 / Accepted: 23 July 2015 / Published: 31 July 2015

Abstract

:
Groundwater-fed lakes in northeastern Germany are characterized by significant lake level changes, but for only a few lakes are in situ water level measurements available. In this study, we test the potential of RapidEye satellite images for indirectly reconstructing lake level changes. The lake levels are derived by intersecting water-land borders with a high-resolution digital elevation model (DEM). Based on Lake Fürstenseer (LF), we define requirements and limitations of the method. Water-land borders were extracted automatically from the 37 RapidEye images available for the period between 2009 and 2014. Otsu’s threshold was used for the NIR band and for the normalized difference water index (NDWI). The results were validated with in situ gauging, contour lines from the DEM, and in situ Differential Global Positioning System (DGPS) measurements of the shoreline. Using an ideal shoreline subset, the lake levels could be reconstructed with decimeter accuracy using the NIR water-land border, but the levels were systematically underestimated by 0–20 cm. The accuracy of the reconstructed lake level retrieval strongly depends on the precision of the water-land border retrieval, on the accuracy of the DEM, and on the lake level itself. A clear shift of the water-land border with increasing lake level is also essential for the unambiguous reconstruction of different levels. This shift needs to be several times larger than the pixel size. The biggest challenges for lake level reconstruction are the presence of vegetation at the shorelines, the quality of the topographic data in the underwater area, the slope of the shoreline, and shadows in combination with low solar angles.

1. Introduction

The glacial landscape of the northeastern Central European lowlands in Germany and Poland is characterized by a high number of natural lakes. Most are groundwater lakes whose water levels are linked to the local uncovered aquifer. Groundwater and precipitation determine the inflow, evapotranspiration the outflow [1,2]. These processes affect the natural dynamic of the lakes and lake level changes. Previous studies showed that groundwater-fed lakes are characterized by small seasonal fluctuations and large perennial changes in lake levels [1,2]. In addition to natural level fluctuations, human interference such as land use and hydromelioration, as well as climate change, influence lake levels [3,4]. Lake level changes influence the lakes’ physics, chemistry and biology. Lakes in the northeastern Central European lowlands play an important role in ecosystem functions [1]. The most important services of lakes are the protection of biodiversity, water and carbon storage as well as their recreational value (tourism) [4].
In the last few years, field measurements at selected lakes in Germany and Poland showed a significant increase in their lake levels after a long period of continuous sinking [1,5]. A lake level decrease leads to the drying out of the shallow water areas at the shoreline, up to the appearance of sandy sediments. As a reaction, macrophytes, reed, and macrozoobenthos move to deeper water areas. If the change is abrupt the existing flora and fauna vanish [6]. The (temporal) loss of typical underwater and shallow water vegetation is a serious problem because of its contribution to water filtering processes and its function as habitat for fauna [7,8]. Large macrophyte stocks help to stabilize clear water conditions [1,9,10]. Their loss can cause a transition from macrophyte-dominated states of lakes with clear water conditions to murky plankton-dominated states [11]. This affects specifically shallow lakes as the smaller lake volume and depth accelerate biological processes and increase eutrophication [12]. Additionally, the loss of vegetation and the drying out of shorelines lead to an increased mineralization of organic matter. If the water level rises and those sediments are flooded again, the mineral substance again dissolves and increases the nutrient content of the water [13].
So far, lake level and shoreline changes have only been monitored at selected lakes. As field monitoring is time consuming and expensive, most lakes in northeastern Germany and Poland are not gauged and their lake level fluctuations could only be estimated. In future, the aim should be to monitor a large number of lakes to build a database for trend analysis. Remote sensing images can be used for lake monitoring and reconstructing past oscillations. Multispectral remote sensing data have been successfully used in the past to monitor and map lakes [14,15,16,17] and to monitor coastline changes [18,19,20,21,22,23]. Lake levels can be measured directly by radar altimetry [24,25,26,27] or indirectly by combining topographic data and water–land boundaries [16,28,29,30,31,32,33]. Indirect reconstruction of lake volumes or levels has been accomplished using aerial photos [33] and medium‑resolution data [16,28]. It has also been used to reconstruct water surface elevations during inundations using various types of remote sensing data [29,30,31,32].
As radar altimeters are profiling tools and their orbital spacing misses too many lakes [31], we use in this study an indirect method to reconstruct lake levels: we test the potential of high-resolution multi-spectral satellite images (RapidEye) for the indirect measurement and reconstruction of small and short-term lake levels changes. The short-term analysis of small lake level changes and high-resolution satellite image data leads to new methodological challenges: in this study, we discuss in detail the influence of shorelines slopes, vegetation and shadows on the accuracy of lake level reconstruction. The large data set for LF, including daily gauging and a shoreline measured in situ, enables absolute validation, whereas previous studies could only compare water-land borders or levels extracted at different times [21,22]. In order to assess the transferability of this lake level reconstruction approach, we examine the attainable accuracy and define requirements and limitations.

2. Study Area

Our test site, LF, is a well-researched groundwater-fed lake located near Neustrelitz in Mecklenburg-Vorpommern, Germany (Figure 1). The lake is part of the Mecklenburg Lake District and the geology of its catchment is sandur from the last ice age [34]. The mean annual precipitation at the closest climate station, Carpin-Serrahn [35], is 656 mm (reference period: 1988–2014). The distribution of the annual precipitation is bimodal with maxima in July (mean: 75 mm) and December (mean: 58 mm) and minima in April (mean: 35 mm) and October (mean: 48 mm).
LF is representative of a large number of groundwater-fed lakes in northeastern Central Europe and it has a heterogeneous shoreline, which is essential for analyzing shoreline slope and vegetation.
The lake does not have an active inlet or outlet and is characterized by significant seasonal and inter‑annual lake level changes [1,34]. The range of annual lake level oscillation is in the range of 0.2 m (reference period: 1988–2014) with a local maximum level in April and a local minimum in October in a typical annual cycle. Periods with high lake levels over 83.8 m a.s.l. (1988/1988, 1995/1996 and 2012/2013) are followed by low water periods with an absolute minimum in 2007 (63.2 m a.s.l.). Consequently, the water surface area of LF varies between 2.1 am 2.6 km2 (GIS analysis based on DEM and in situ level measurements). Between 2009 and 2014, there is a positive trend in lake level due to extraordinarily high precipitation. The gauge data since 1988 and the decomposition of the time series are illustrated in detail in Figure S1.
The average depth of the lake is 8.0 m; the maximum depth is 24.5 m at a lake level of 63.5 m a.s.l. (bathymetric survey). The lake is a mesotrophic clear water lake characterized by water stratification [36]. The maximum sight depth, measured via secchi disk between 2004 and 2012, ranges from 4 to 8 m [37]. In winter, the lake is on average completely covered with ice on 51 days and partly covered on 14 days (reference period: winter 2009/2010–winter 2013/2014) [38].
Figure 1. The northern Central European lowlands are characterized by a large number of water bodies. Lake Fürstenseer (LF) is located near Neustrelitz in Mecklenburg-Vorpommern, Germany. The close-up shows a RapidEye satellite image of LF in September 2009 (low lake level). The orange rectangle marks the region of interest for calculating the Otsu thresholds. The red circles mark different subsets of the shorelines (A–E), and the yellow triangle shows the position of the Fürstensee gauging station.
Figure 1. The northern Central European lowlands are characterized by a large number of water bodies. Lake Fürstenseer (LF) is located near Neustrelitz in Mecklenburg-Vorpommern, Germany. The close-up shows a RapidEye satellite image of LF in September 2009 (low lake level). The orange rectangle marks the region of interest for calculating the Otsu thresholds. The red circles mark different subsets of the shorelines (A–E), and the yellow triangle shows the position of the Fürstensee gauging station.
Water 07 04175 g001

3. Materials and Methods

3.1. Data Acquisition and Preparation

The LF database includes in situ lake levels, a digital surface model of the terrestrial area adjacent to the lake, bathymetric point data for the lake, and Differential Global Positioning System (DGPS) measurements of the shoreline in August 2014 (Table 1). A workflow of the preparation of the topographic data and the RapidEye images is illustrated in Figure S2.
The Fürstensee gauging site is located at the new public beach (Figure 1, subset A). The processing of the lake level data includes the transfer of the measured relative lake level (in cm) to absolute level values in m above sea level (a.s.l.). The transfer is based on height information provided by the Staatliches Amt für Landwirtschaft und Umwelt Mecklenburgische Seenplatte [38].
Digital surface models with 1 m resolution, locally only 2 m resolution, are available nationwide from German government agencies. The digital surface model with 1 m resolution (ATKIS-DGM1) covering LF is based on two LIDAR acquisitions in winter 2010.
The bathymetric survey took place on 7 October 2002. Despite the time since the survey, the bathymetry is assumed to be still representative because of the closed basin and thus very low levels of erosion and sedimentation at the lake. The survey was conducted by a surveying company on behalf of Landesamt für innere Verwaltung Mecklenburg Vorpommern [39]. In a boat, the surveyors navigated in rows over the lake and measured the relative depth (z) from the lake bottom to the water surface. For the transformation to absolute height values, the z values were subtracted from the in situ gauged lake level in October 2002 (63.53 m). Ca. 46,000 points were measured on the 2.4 km2 lake surface. The distance between rows was 50 m and the shoreline was not mapped in situ. Based on the in situ measured level, the shoreline was reconstructed as a contour from the ATKIS-DGM1. The shoreline was added to the points of the bathymetric survey as points with zero depth. Then, the underwater surface model was calculated as a Triangulated Irregular Network (TIN). The TIN was resampled to a raster with 1 m pixel grid size. The accuracy of the underwater surface model is highly dependent on the density of measured points.
Table 1. Overview of the available data.
Table 1. Overview of the available data.
DataAbbreviationMethod/DeviceAcquisition DateResolution (Accuracy)Source
In situ measured lake levelManual reading at gauging sitemonthly since 1987, daily since 20061 cmStaatliches Amt für Landwirtschaft und Umwelt Mecklenburgische Seenplatte (MS) [38]
Digital surface modelATKIS-DGM1Pre-processed LiDAR1 November 2010 and 15 December 20101 m (vertical: 0.15–0.2 m)Landesamt für innere Verwaltung Mecklenburg-Vorpommern [40]
Bathymetric point dataSonar via SIMRAD-Echolot7 October 2002(Horizontal: 1 m, vertical: 0.1 m)Ministerium für Landwirtschaft, Umwelt und Verbraucherschutz M-V [37]
In situ measured water-land borderDifferential Global Positioning System via Triumph-VS Receiver12 August 2014(Horizontal RMSE < 1 m)
Finally, an over- and under-water digital elevation model (DEM) of LF was generated by merging the ATKIS-DGM1 and the bathymetric data. Below 63.53 m, the DEM was based on the bathymetric survey and above 63.53 m on the ATKIS-DGM1 (cf. Figure S4 and Figure 6). The DEM has a 1 m resolution.
The water–land borders of LF were extracted from RapidEye (RE) satellite images. The data were provided by BlackBridge AG, Berlin, Germany (www.blackbridge.com). RapidEye images were chosen because of their high spatial resolution of 6.5 m, their high temporal resolution, the large size of the archive including data collected since 2009, and the free-of-charge availability for German research projects using the RapidEye Research Archive. Five micro-satellites have provided products since February 2009 with a theoretical temporal resolution of one day (off-nadir) or 5.5 days (at nadir). The sensors on the satellites acquire five spectral bands at wavelengths between 440 and 850 nm, including the red edge band (band 1: 440–510 nm, band 2: 520–590 nm, band 3: 630–685 nm, band 4: 690–730 nm, and band 5: 760–850 nm) [39]. For this study the images were provided in processing level 3A, delivered as rectangular image tiles. Level 3A processing includes orthorectification and resampling to 5 m × 5 m pixels, as well as radiometric, geometric, and terrain correction of the images [41].
A total of 108 scenes are available showing LF between April 2009 and October 2014. However, 71 scenes were excluded from further analysis due to ice cover, snow or a closed cloud cover above the lake. The remaining 37 suitable RapidEye images collected between April 2009 and October 2014 are listed in Table 2. The further preprocessing of the images included the calculation of the top of atmosphere reflectance (TOAR) [41], the co-registration of the stack of images, and cloud masking.
Table 2. Available RapidEye archive data showing LF between April 2009 and October 2014, including the metadata for the acquisitions. The lack of data for 2013 is due to the scarcity of acquisitions because of cloud cover and ice.
Table 2. Available RapidEye archive data showing LF between April 2009 and October 2014, including the metadata for the acquisitions. The lack of data for 2013 is due to the scarcity of acquisitions because of cloud cover and ice.
Acquisition DateSatelliteSensorSensor Viewing Angle (°)Sun Elevation Angle (SEA) (°)Lake Levels (m a.s.l.) Measured in situ
2009-04-04RE-1MSI0.0642.5963.37
2009-04-13RE-1MSI13.4345.8163.36
2009-04-21RE-4MSI6.6948.7263.34
2009-08-31RE-2MSI−3.0545.3163.18
2009-09-20RE-3MSI−2.9037.7563.15
2010-06-03RE-2MSI17.0359.1263.4
2010-06-17RE-2MSI20.4660.0663.39
2010-07-03RE-3MSI3.6859.7463.34
2010-07-19RE-5MSI10.3857.5663.29
2010-09-22RE-3MSI3.5337.0563.32
2010-10-04RE-1MSI13.1632.4063.31
2011-04-20RE-3MSI−2.9448.2663.54
2011-05-07RE-1MSI3.4953.5463.51
2011-05-11RE-5MSI−3.1854.5963.5
2011-05-30RE-5MSI−2.9558.5063.52
2011-06-04RE-5MSI−9.7959.1163.51
2011-06-27RE-5MSI13.6060.1063.53
2011-09-24RE-3MSI6.8836.3463.7
2011-10-02RE‑1MSI−2.9633.0563.7
2011-10-13RE-3MSI7.0729.0063.7
2011-10-17RE-2MSI3.9127.4263.7
2011-10-22RE-2MSI−6.1125.5263.7
2011-11-13RE-5MSI−6.2118.6363.68
2012-04-05RE-1MSI−9.6143.0563.92
2012-05-01RE-4MSI10.2552.0863.94
2012-05-23RE-2MSI10.2457.4863.91
2012-06-18RE-4MSI−2.9560.2163.84
2012-07-24RE-2MSI7.0356.5363.83
2012-10-12RE-1MSI10.0029.0763.77
2012-11-14RE-5MSI−6.2518.1463.77
2014-03-10RE-5MSI−5.8732.7763.86
2014-03-20RE-1MSI6.6936.7063.86
2014-04-25RE-3MSI−13.1149.9563.87
2014-05-01RE-5MSI0.3151.9363.86
2014-05-20RE‑5MSI0.3456.8263.84
2014-08-06RE-2MSI3.5453.4463.76
2014-09-05RE-3MSI−9.7743.5263.72
To analyze changing water–land borders, precise locating of the RapidEye images is required. We chose a well-positioned master RapidEye image based on a comparison with digital orthophotos. All RapidEye slave images were co-registered to this master image using the script of Behling et al. [42]. Co-registration was based on 100 points per image tile. The maximum root mean squared error (RMSE) of the co-registered slave images to the master image ranges from 0.04 to 1.46 m and is on average 0.75 m. Finally, clouds were masked using two thresholds in band 1 (blue) and band 3 (red) improved by an additional buffer area of 50 m.
Finally, we measured in situ the shoreline in August 2014 via DGPS. The DGPS measurement of the water–land border took place on 12 August 2014, close to the RapidEye acquisition on 6 August 2014. Additionally, a description of the vegetation at the shoreline and in the water was acquired. The correction of the GPS signal was performed online by the Landesamt für innere Verwaltung Mecklenburg-Vorpommern via SAPOS-HEPS during the field measurement [43]. However, high vegetation at the northern part of the lake led to gaps in the mobile phone network and to inaccuracies in the online correction of the DGPS signal. We removed all points with a horizontal root mean square error >1 m to maintain high accuracy.

3.2. Method

The lake’s water level was reconstructed by intersecting the water–land border with a digital elevation model (DEM), which results in the absolute height of the lake surface (=the lake level) [16,32,33]. The principle is demonstrated in Figure 2. Figure 2b also illustrates the challenges caused by vegetation and steep slopes that reduce the distances between contour lines. In practical application, inaccuracies in the extraction of the water-land borders or in the DEM must also be considered. While in the ideal case the water-land border lies on one contour line, in practice the water-land border cuts through several contour lines. Thus, the intersection of the water-land border and DEM results in a range of height values for each date. After a statistical analysis of the height values, the mode value (=most frequent value) for each date is taken as the according reconstructed lake level.
Figure 2. (a) Scheme of the side view of a lake including its topography; (b) scheme of the top view of a lake including its topography (contour lines). The top view illustrates a satellite image in which the water-land border is delineated. In an ideal case, the water-land border lies on one contour line, in this example at 64 m. This height of the water-land border (=contour line) equals the lake level.
Figure 2. (a) Scheme of the side view of a lake including its topography; (b) scheme of the top view of a lake including its topography (contour lines). The top view illustrates a satellite image in which the water-land border is delineated. In an ideal case, the water-land border lies on one contour line, in this example at 64 m. This height of the water-land border (=contour line) equals the lake level.
Water 07 04175 g002
Two state‑of‑the‑art approaches for distinguishing water-land borders in multispectral images are thresholds in the NIR band [15,44,45] and the use of the Normalized Difference Water Index (NDWI) by McFeeters [46,47,48,49]. The formula for calculating NDWI according to McFeeters (1996) [50] is:
N D W I = R e f ( g r e e n ) R e f ( N I R ) R e f ( g r e e n ) + R e f ( N I R )
where Ref is Reflectance.
To delineate the water-land border in the RE images and reconstruct former levels we tested and compared both approaches, the thresholding of the NIR band and the use of NDWI. To avoid subjectivity in the choice of the threshold and to maximize the level of automation, we used the automatic Otsu threshold [49,51]. In a bimodal histogram, the Otsu algorithm calculates the optimum threshold separating the two classes by minimizing the intra‑class variance or maximizing the inter‑class variance [52]. To calculate the Otsu thresholds, we selected a subset of the southern part of LF and its surroundings (Figure 1). The histograms of NIR and NDWI for this region of interest showed a clear bimodal distribution with the two classes water and land. Based on this threshold, the entire RapidEye image was classified and the shoreline of LF converted to vectors. To reconstruct the lake level, only the longest water-land border of LF was selected and small misclassified polygons were removed.
The extracted water-land borders (NIR and NDWI) from the RapidEye image of 6 August 2014 were validated by the water-land border determined in situ with the DGPS signal. As a degree of precision, we calculated the minimum distances between each DGPS point and the extracted RE water-land borders. Additionally, in the time series we compared all extracted RapidEye water-land borders to the corresponding contour line from the DEM and the in situ gauged lake level. For the evaluation, we separated each RapidEye water-land border into a row of points with 1 m distance and based on this separation calculated from the DEM the minimum distance to the corresponding contour lines.
The final step was the intersection of the extracted RapidEye water-land border with the high-resolution DEM for the lake level reconstruction. As the shoreline is very heterogeneous with respect to land cover and shoreline topography, we selected and compared different subsets, as recommend by Kaiser et al., 2015 [33]. The subsets were selected based on shoreline slope and vegetation. The slope of the shoreline of LF is defined as the slope (in percent) between 63 and 64 m a.s.l. The range from 63 to 64 m a.s.l. represents the rounded range of lake level changes in the last decade. The selected subsets were approximately 40 m wide to have a minimum size for a statistical analysis of the reconstructed level values.
The entire processing chain of water classification and level reconstruction was implemented using the free software R (version 3.0.1) and is illustrated as a workflow in Figure S2.

4. Results

4.1. Automatic Water-Land Border Extraction

The comparison of measured DGPS points and the extracted water-land borders (NDWI and NIR) in August 2014 are illustrated in Figure 3a,b. The boxplots in Figure 3a represent the calculated minimum distances between the DGPS points and the extracted water-land borders using NDWI and NIR. Negative distances imply that the extracted water-land border is inside the in situ determined boundary (underestimation of the water extent); positive distances imply that the extracted water-land border is outside the in situ determined boundary (overestimation of the water extent).
Figure 3. (a) shows the minimum distances between the DGPS points and the extracted water-land borders (NDWI and NIR) at different subsets (cf. Figure 1). Negative distances show an underestimation of the extracted water-land border, meaning that the extracted water-land border is inside the in situ determined boundary. Positive distances show an overestimation, meaning that the extracted water-land border is outside the in situ determined boundary. The black dots represent the median value; (b) illustrates the reconstructed and in situ determined water-land borders at the southeastern part of LF. The reconstructed water-land borders were mapped in the RapidEye image using Otsu’s threshold in the NIR (pink line) and the NDWI (yellow line) images. The in situ determined water-land border was measured via DGPS (light blue dots).
Figure 3. (a) shows the minimum distances between the DGPS points and the extracted water-land borders (NDWI and NIR) at different subsets (cf. Figure 1). Negative distances show an underestimation of the extracted water-land border, meaning that the extracted water-land border is inside the in situ determined boundary. Positive distances show an overestimation, meaning that the extracted water-land border is outside the in situ determined boundary. The black dots represent the median value; (b) illustrates the reconstructed and in situ determined water-land borders at the southeastern part of LF. The reconstructed water-land borders were mapped in the RapidEye image using Otsu’s threshold in the NIR (pink line) and the NDWI (yellow line) images. The in situ determined water-land border was measured via DGPS (light blue dots).
Water 07 04175 g003
Eye-catching in Figure 3a are the ranges of negative and positive deviations. High deviations for NIR in subsets B and D result from misclassifications caused by the dark cloud shadows in the RapidEye image (6 August 2014). The NIR water-land border spreads massively into the land area with deviations up to 80 m. The NDWI water-land border also shows misclassifications in the land area, but the lake shoreline remains intact. Large negative deviations are caused by dense reed and overhanging trees (subset E). Subset B is not completely free of vegetation as reed grows there sparsely (mapped in summer 2014). There are beech trees growing around the lake. In summer 2014, at a level of 63.77 m, the water was close to, but not yet under the trees. Here, the deviations between DGPS and NIR/NDWI water-land borders are caused by inaccuracies in water extraction. The overestimation of the NIR water-land border at subset B is mainly caused by cloud shadow, but also in the smaller subset without direct cloud shadow the NIR water-land border is overestimated on average by 4 m (range: 0 m–7 m). The NDWI water-land border at the subset with 2.8% slope differs on average by 1 m from the DGPS-derived shoreline. The height differences between the DGPS-derived shoreline and the RapidEye water-land borders range from 0.05 to 0.25 m (NIR) and 0 to 0.20 m (NDWI) at the subset with 2.8% slope.
The new public beach (subset A) is only 50 m wide, with a very small sandy stripe (ca. 2.5 m in August 2014) and lawns behind. Here, the extracted water-land borders using NDWI and NIR are identical, with an underestimation of 0 m–10 m. The height difference between the DGPS shoreline and the NDWI/NIR water-land borders is between 0 and 0.25 m.
Not only cloud shadow, but also shadows behind rows of trees and buildings caused misclassifications using the NIR image. These small misclassified polygons were removed for the reconstruction of the lake level and only the longest line was taken as the water-land border.

4.2. Shoreline Changes

In addition to the validation using DGPS measurements, the extracted water-land-borders (NDWI and NIR) were compared to the according contour lines at different lake levels. Figure 4 illustrates the water-land borders and the contour lines at minimum level (63.15 m a.s.l.) on 20 September 2009, two intermediate levels, and the near-maximum level (63.91 m a.s.l.) on 23 May 2012 (cf. Figure S1). The changes of the water-land border at the southeastern beach of LF are illustrated in smaller scale in Figure S3. The lake levels on 1 May 2012 (63.94 m a.s.l., maximum level) and 5 April 2012 (63.92 m a.s.l.) were slightly higher, but clouds covered large parts of LF.
The contour lines from the DEM (Figure 4c) show clear changes of the shoreline with increasing lake level. In the southwestern very shallow area (slope: 0%–3%), the contour lines move up to 150 m between the minimum and maximum extent. The sandbank that is visible in the RapidEye image is, however, not detected. This error is caused by the large gap of bathymetric points in the shallow water area and the resulting significant underestimation of the height in this area (cf. Figure S4).
On the vegetation free southeastern beach, the NIR water-land borders (Figure 4a) clearly show the changes in lake level. The NIR water-land border moves up to 35 m between the minimum and maximum extent, whereas the NDWI water-land borders (Figure 4b) are clearly overestimated at the low and intermediate lake levels (cf. the close-up in Figure S3). The position of the NDWI water-land border at low water levels (20 September 2009) fits neither the GIS analysis nor photos from 2009 that show a wide sand beach (Figure S5).
The sandbank in the middle of the southern part of LF is detected via NIR and NDWI. In the NIR image it has a size of ca. 3300 m2 at minimum water level (20 September 2009) and shrinks with increasing levels to ca. 300 m2 on 21 April 2009. In the NDWI image it is at ca. 300 m2 also much smaller than in the NIR water-land border and only detected at the minimum level.
Due to the importance of topographic data for lake level reconstruction, Figure 5 illustrates the NIR water-land borders and the according contour lines at low and high levels as well as the topographic data. The underwater surface model is based on bathymetric points and the contour line on the day of the bathymetric survey. The gap between the point data and the contour line ranges in the shallow water area from ca. 15 m to ca. 40 m, depending on the slope of the shoreline. It is unavoidably larger in the shallower subsets.
Whereas the low lake level falls in the area of the underwater surface model, the high level falls in the area covered by ATKIS-DGM1. Thus, these two quality levels of the topography need to be kept in mind for the further analysis. The distance between the contour line and the NIR water-land border at low lake levels is larger than the distance between both at high lake levels (cf. Figure 6). In the subset with 2.7% slope the high deviation between the extracted NIR water-land border and the contour line sticks out. This deviation is caused by flooded trees. It must also be considered that there is a deviation between the NIR water-land borders and the contour lines because of the different resolutions of the RapidEye image and the DEM.
Figure 4. Illustration of the water-land borders on 21 April 2009, 20 September 2009, 20 April 2011, and 23 May 2012. (a) shows the water-land borders extracted from the RapidEye images using the NDWI; (b) shows the water-land borders extracted from the RapidEye images using the NIR; (c) represents the GIS analysis with contour lines: The levels measured in situ on the four dates are delineated as contour lines in the DEM. The minimum level (63.15 m a.s.l.) between 2009 and 2014 was measured on 20 September 2009, the maximum level on a cloud free image (63.91 m a.s.l.) on 23 May 2012.
Figure 4. Illustration of the water-land borders on 21 April 2009, 20 September 2009, 20 April 2011, and 23 May 2012. (a) shows the water-land borders extracted from the RapidEye images using the NDWI; (b) shows the water-land borders extracted from the RapidEye images using the NIR; (c) represents the GIS analysis with contour lines: The levels measured in situ on the four dates are delineated as contour lines in the DEM. The minimum level (63.15 m a.s.l.) between 2009 and 2014 was measured on 20 September 2009, the maximum level on a cloud free image (63.91 m a.s.l.) on 23 May 2012.
Water 07 04175 g004
Figure 5. Digital elevation model, NIR water-land borders and corresponding contour lines at the southeastern part of LF (subsets B and C). The red and blue lines illustrate the deviation between the reconstructed water-land borders and the contour lines at the minimum lake level (20 September 2009, 63.15 m a.s.l.) and its maximum level in a cloud free image (23 May 2012, 63.91 m a.s.l.). The toothed shorelines (dark red and dark blue) were extracted from the RapidEye NIR image. The smoother lines (bright red and bright blue) are the contour lines from the DEM based on levels measured in situ. The green line illustrates the contour line on the day of the bathymetric survey (63.53 m a.s.l.). The measured bathymetric points are the black dots. Different subsets and their slopes (in percentage) are noted.
Figure 5. Digital elevation model, NIR water-land borders and corresponding contour lines at the southeastern part of LF (subsets B and C). The red and blue lines illustrate the deviation between the reconstructed water-land borders and the contour lines at the minimum lake level (20 September 2009, 63.15 m a.s.l.) and its maximum level in a cloud free image (23 May 2012, 63.91 m a.s.l.). The toothed shorelines (dark red and dark blue) were extracted from the RapidEye NIR image. The smoother lines (bright red and bright blue) are the contour lines from the DEM based on levels measured in situ. The green line illustrates the contour line on the day of the bathymetric survey (63.53 m a.s.l.). The measured bathymetric points are the black dots. Different subsets and their slopes (in percentage) are noted.
Water 07 04175 g005
Figure 6. Boxplot of the minimum distances (in m) between the RapidEye water-land borders (NIR) and the according DEM contours at subset B. The black dots represent the median. The root mean squared error (RMSE) for all dates is 21.1 m. The RMSE is reduced to 10.8 m when dates with a very low sun elevation angle (SEA) below 30° and dates with a SEA angle below 45° in cases of high lake levels (≥63.7 m a.s.l.) are excluded.
Figure 6. Boxplot of the minimum distances (in m) between the RapidEye water-land borders (NIR) and the according DEM contours at subset B. The black dots represent the median. The root mean squared error (RMSE) for all dates is 21.1 m. The RMSE is reduced to 10.8 m when dates with a very low sun elevation angle (SEA) below 30° and dates with a SEA angle below 45° in cases of high lake levels (≥63.7 m a.s.l.) are excluded.
Water 07 04175 g006
A comparison of the contour lines from the DEM with the extracted NIR water-land borders of the entire time series is illustrated in Figure 6. It is based on subset B. On 24 of the 37 dates the RapidEye water-land border (NIR) compared with the median of the distances to the according DEM contour line is negative, in the range of −2 m to −10 m. At near-maximum lake level on 23 May 2012 the NIR water-land border is underestimated on average by −4.1 m, at the minimum lake level on 20 September 2009 by −9.7 m. Thirteen dates of the time series are characterized by (high) positive median distances up to several tens of meters. One of those outliers is 6 August 2014, the date of the field campaign with misclassifications due to cloud shadow. Data collected on 5 April 2012 is heavily disturbed by cloud shadow, data for 19 July 2010 only slightly. The six images between 12 October and 15 November are all characterized by a low sun elevation angle (SEA) <30°. Due to the higher noise in those images, the automatic classification fails and the classification spreads into a dark clearing in the surrounding forest (cf. Figure S6). The remaining images that overestimate the water surface area are dated 24 September 2011, 2 October 2011, 5 April 2012, 10 March 2014, 20 March 2014, and 24 September 2014. These dates have in common a SEA of 30° and 45° and a high lake level ≥63.7 m a.s.l. At this high lake level, the sandy beach is completely flooded. By contrast, images with a sandy beach (4 April 2009, 20 September 2009, 22 September 2010, and 4 October 2010) and a SEA between 30° and 45° underestimate the water surface area. A comparison of the contour lines with the extracted NIR water-land borders using the subset with 2.8% slope (ideal subset) can be found in Figure S7. Using the ideal subset the root mean squared error (RMSE) for all dates is 13.5 m. After excluding the 12 images with low SEA, the RMSE sinks to 5.0 m.

4.3. Subset Selection for Lake Level Reconstruction

A precise extraction of the water-land border is essential for an accurate reconstruction of the lake level. Another requirement for lake level reconstruction is the visibility of a water-land border change. Table 3 illustrates the shift of the water-land border at the minimum lake level to the water-land border at the maximum lake level in subsets with different slopes. The position of the subsets at the southeastern shoreline is shown in Figure 5. The subset with 8.7% slope is located at the new public beach (cf. Figure 1). With shallower slope the distance between the minimum and maximum water-land border increases for the extracted water-land borders and for the contour lines from the DEM. The small shift of 5 m (=1 pixel) is too small for a reliable result, as intermediate lake extents cannot be delineated unambiguously. A shift of 10 m (=2 pixels) at least is theoretically required for the reconstruction of intermediate levels. At LF this means that only subsets with slopes ≤7.1% are suitable for lake level reconstruction.
Table 3. Average distances (in m) between the water-land borders at the minimum lake level (20 September 2009, 63.15 m a.s.l.) and the maximum level of a cloud-free image (23 May 2012, 63.91 m a.s.l.) in subsets with various slopes (in percent).
Table 3. Average distances (in m) between the water-land borders at the minimum lake level (20 September 2009, 63.15 m a.s.l.) and the maximum level of a cloud-free image (23 May 2012, 63.91 m a.s.l.) in subsets with various slopes (in percent).
Slope (%)Distance between Minimum and Maximum NIR Water-Land Border (m)Distance between Minimum and Maximum Contour Lines (m)Location (cf. Figure 1 and Figure 3)Notes
2.730.140.2Subset CTrees flooded in August 2014, high lake level
2.829.629.0Subset B
3.729.121.6Subset B
4.921.719.1Subset B
5.122.216.3Subset B
7.110.210.8Subset B
8.75.38.9Subset A
As this study aims at defining the most attainable accuracy, we conducted the reconstructed lake levels based on an “ideal” subset. The largest distances between the minimum and maximum water surface area are in the two shallowest subsets with 2.7% and 2.8% slopes (cf. Table 3). However, only the subset with 2.8% slope is more or less free of vegetation. Thus, this subset was selected as the ideal subset.

4.4. Lake Level Reconstruction

Although validation via DGPS showed more accurate results for the NDWI then for the NIR water-land border extraction, the time series analysis revealed an overestimation of the NDWI lake extent at low lake levels (Chapters 4.1 and 4.2). Thus, the following lake level reconstruction is based on the NIR water-land-borders and the ideal subset (cf. Chapter 4.3).
Even in the case of the ideal subset, the intersection of the water-land border and the DEM results in a range of height values for every image/date. This is due to the different spatial resolutions of DEM and RapidEye images and small inaccuracies in the delineation of the water-land border. The range of height values per date is illustrated in Figure 7. A comparison of the mean, median and mode height values shows that the mode value is the most robust parameter. The mode value best represents the lake level, even in cases of misclassification of the water-land-border (cf. dates 13 October 2011, 5 April 2012 and 10 March 2014). Thus, the mode value of the height values for one date is considered to be its lake level.
The reconstructed lake level is underestimated for 28 of 37 dates. The underestimate ranges from −20 to −2 cm, the average being −11 cm. For three dates the precise in situ lake level is reconstructed (19 July 2010, 2 October 2011, and 6 August 2014); however, two of those images are disturbed by cloud shadow (yellow highlighting in Figure 7). For six dates (17 October 2011, 22 October 2011, 13 November 2011, 12 October 2012, 14 November 2012, and 20 March 2014) the reconstructed lake levels are overestimated in the range of 6 cm to 79 cm. The overestimations are caused by misclassifications in the extraction of the water-land border. As discussed before, the automatic detection of the water-land border has a reduced accuracy for images acquired at very low SEA <30° and for images with a SEA <45° in cases of high lake level (cf. Table 2). In Figure 7, all these images are highlighted in red. With removal of the 12 images, the RMSE between in situ measured and reconstructed levels is reduced from 22 to 12 cm, and the absolute mean error from 14 to 11 cm.
In a second step, we analyzed the relative level changes Δ h = h ( t ) h ( t + 1 ) of dates without misclassifications in the ideal subset in Figure 8. For the monitoring period, the absolute mean of difference of Δ h ( i n   s i t u ) and Δ h ( r e c o n s t r u c t e d ) is 5 cm, the RMSE 6 cm. In 2009 and 2010, the lake level is so low that there is a wide sandy area between the lake and surrounding vegetation. At those times, the in situ measured and the reconstructed lake level changes are very similar: the difference between Δ h ( i n   s i t u ) and Δ h ( r e c o n s t r u c t e d ) ranges only about ±5 cm (absolute average = 3 cm, RMSE = 3 cm). The only larger deviations on 19 July 2010 and 22 September 2010 can be explained by cloud shadows on image 19 July 2010 and the resulting overestimation of the water level. In 2011 the lake level increases significantly and the sandy beach is flooded. Also, in 2011 the deviations between Δ h ( i n   s i t u ) and Δ h ( r e c o n s t r u c t e d ) grow larger (absolute average 5 cm, RMSE = 6 cm), with a maximum delta of −13 cm. However still, the reconstructed level changes are close to the trends of the changes determined by in situ measurement. Since 2012, the year with the highest lake level, the accordance of reconstructed and in situ measured level changes decreases significantly. The difference of Δ h ( i n   s i t u ) and Δ h ( r e c o n s t r u c t e d ) ranges from −7 cm to +19 cm (absolute average 6 cm, RMSE = 8 cm). Several times the reconstructed and the in situ measured lake levels show opposite lake level changes. The most conspicuous outlier (6 August 2014) is caused by cloud shadow.
Figure 7. Illustration of the in situ measured and reconstructed lake levels at LF. The boxplots show the range of height values based on the intersection of the water-land border in the ideal subset with 2.8% slope with the high-resolution DEM. The mode value of a boxplot is considered to be the lake level at the according date. The blue line shows the in situ measured lake level ranges between January 2009 and November 2014. The red dots are the mode value, the light blue dots the median value for each date. Red highlighted are all dates with a very low sun elevation angle (SEA) below <30° and dates with a SEA angle below <45° in cases of high lake levels (≥63.7 m a.s.l.). Yellow highlighted RE images are disturbed by significant cloud shadows.
Figure 7. Illustration of the in situ measured and reconstructed lake levels at LF. The boxplots show the range of height values based on the intersection of the water-land border in the ideal subset with 2.8% slope with the high-resolution DEM. The mode value of a boxplot is considered to be the lake level at the according date. The blue line shows the in situ measured lake level ranges between January 2009 and November 2014. The red dots are the mode value, the light blue dots the median value for each date. Red highlighted are all dates with a very low sun elevation angle (SEA) below <30° and dates with a SEA angle below <45° in cases of high lake levels (≥63.7 m a.s.l.). Yellow highlighted RE images are disturbed by significant cloud shadows.
Water 07 04175 g007
Figure 8. The black line shows the relative changes Δ h ( i n   s i t u ) of the in situ measured level, the red line the relative changes Δ h ( r e c o n s t r u c t e d ) of the reconstructed “mode” levels. Twelve dates were removed because of low sun elevation angles.
Figure 8. The black line shows the relative changes Δ h ( i n   s i t u ) of the in situ measured level, the red line the relative changes Δ h ( r e c o n s t r u c t e d ) of the reconstructed “mode” levels. Twelve dates were removed because of low sun elevation angles.
Water 07 04175 g008
Finally, we compare the relative changes of the ideal subset with the results of the five other subsets with various slopes (cf. Figure 5). The results are depicted in Figure 9.
Figure 9. Relative changes ( Δ h = h ( t ) h ( t + 1 ) ) of the reconstructed lake levels (“mode values”). The black dashed line is the relative change of the level measured in situ. The colored lines represent the relative changes of the reconstructed lake levels.
Figure 9. Relative changes ( Δ h = h ( t ) h ( t + 1 ) ) of the reconstructed lake levels (“mode values”). The black dashed line is the relative change of the level measured in situ. The colored lines represent the relative changes of the reconstructed lake levels.
Water 07 04175 g009
As expected, the shallower the slopes, the higher the precision in the reconstruction of the lake level changes: the level changes using the two shallowest subsets are very similar until the end of 2011. After 2011, the shallowest subset (2.7% slope) shows larger deviations and partly contrary changes. As stated above, the inaccuracies in the shallowest subset are most likely caused by vegetation that covers the water surface. The subset with 3.7% slope still shows the same ups and downs as the lake level change measured in situ, but generally has larger deviations and outliers in 2014. The three steepest subsets (slopes of 4.9%, 5.1% and 7.1%) are all characterized by large deviations and outliers, most eye-catching at 20 April 2011.

5. Discussion

This study suggests that high‑resolution RapidEye archive data are a valuable source for reconstructing former lake levels. It also demonstrates the importance of the shoreline slope, vegetation, and shadows (and thus solar angle) for accurately reconstructing the lake level.
Under ideal conditions, the lake levels can be measured indirectly with the NIR water-land borders in decimeter resolution. This is a significant improvement compared to previous lake level retrievals based on medium resolution data with only one-meter precision [16]. Annual and seasonal fluctuations as well as short-term (within two weeks) jumps in lake level were detected. However, for a very high-resolution temporal monitoring of daily changes or changes within one week, the changes in lake level were too small. The maximum level changes within one week are only in the range of ±5 cm, with small changes in the water-land border. The decimeter accuracy (RMSE 12 cm, absolute mean error 11 cm) of the indirect lake level reconstruction is in line with the accuracy of flood level estimations based on high‑resolution data [53,54] and even better than a previous study in northeastern Germany using aerial photos [33]. However, the levels were systematically underestimated in the range of −20–0 cm (excluding images with sever misclassification of the water-land border). The underestimation of the lake levels is caused by the underestimation of the water extent in the automatic extraction of the water-land border. Due to this systematic underestimation of absolute lake levels, an analysis of relative changes is preferred. The accuracy of the relative changes is generally high (RMSE = 6 cm), but difference between times with sandy beaches and times with flooded beaches became visible: whereas the accuracy in 2009 and 2010 with a sandy beach was very high with (RMSE = 3 cm), it declined massively to approx. ±20 cm (RMSE 6 cm) with higher lake levels and a flawless transition of water to vegetation. This shows that the accuracy of the lake level retrieval depends on terrain and vegetation cover on the lake shore. Thus, the accuracy of the lake level retrieval can change over the monitoring period due to changing shoreline conditions at different dates and lake levels.
Essential for an accurate lake level reconstruction is the precise delineation of the water-land border and precise topographic data. Other studies just compare their extracted coastline at different times without absolute validation of the shoreline precision [21,22]. We validated the delineated water-land borders with DGPS measurements and with contour lines from the DEM. The water-land borders were extracted from the RE images automatically using Otsu’s threshold on the NIR bands and the NDWI images. Whereas the NDWI water-land border had a higher accuracy in the DGPS validation, the time series revealed significate overestimations of the water extent, which fits to the findings of Li et al., 2013 [49]. The extraction of the water-land borders via NIR showed a higher sensitivity to (cloud) shadows, but was also more precise in the delineation of water at low water levels. Thus, we decided to use NIR for the lake level reconstruction.
The sensitivity of NIR and NDWI to shadows is a well-known problem and causes misclassifications [14,16,47,55,56,57]. Shadows of trees and buildings caused small misclassifications that were easily filtered by only taking the longest line as the water-land border. Analysis of the NIR water-land borders and the reconstructed lake levels showed that cloud shadows, such as in the RapidEye image contemporary to the field campaign with the DGPS measurements, could lead to an overestimation of water extent and thus to a relatively higher reconstructed lake level (smaller or no underestimation). This finding confirms the importance of validation with contour lines.
The accuracy of the NIR water-land borders (without misclassified images due to low solar angles) in comparison to the according contour lines at the ideal subset with 2.8% slope is RMSE = 5.0 m; for the subset B the RMSE = 10.8 m. A deviation between the NIR water-land border and the contour line is expected due to the different spatial resolutions of the RapidEye images and the DEM. In a previous study with a very high resolution IKONOS image, the RMSE between the delineation of sandy coastlines and the in situ measured water-land border was approx. 6 m [20]. Another study with medium-resolution Landsat data estimates its accuracy as 1.3 pixels (pixel size = 30 m) [18]. Systematic underestimation of the NIR water-land border is caused by the inherent sensitivity of NIR reflectance in littoral zones with very shallow water columns and by the spectral mixture of water and land [17,58]. Underestimation of the water-land border via NIR band is confirmed by Wang et al., 2014 [17]. However, it needs to be considered that the contour lines in the GIS analysis in this study probably differ from the in situ determined water-land borders due to inaccuracies in the DEM, specifically in the area of the underwater surface model and because of the different resolutions of the RapidEye image and the DEM.
Essential for the unambiguous reconstruction of different levels is also a clear shift of the water-land border with increasing lake level. This shift needs to be several times larger than the pixel size. This study finds that for RapidEye images a 30 m (=6 pixels) shift between minimum and maximum levels (80 cm level difference) is necessary for reconstructing lake levels with decimeter accuracy. A 30 m shift is reached at shoreline slopes <4%. A smaller shift significantly reduced the accuracy of lake level reconstruction, so only a trend in lake level changes could be estimated. The ideal subset for lake level reconstruction at LF has a 2.8% slope.
In summary, the requirements for lake level reconstruction with decimeter accuracy are:
  • The water-land border needs to be delineated precisely. The automatic Otsu threshold on the NIR band showed the best results, but the approach is sensitive to shadows and thus to low solar angles. Dense vegetation cover hinders the accurate retrieval of water-land borders.
  • The topographic data needs to be accurate as the lake level reconstruction is only as good as the underlying DEM. This is specifically true for the underwater surface model, which is not of very high quality due to sparse bathymetric point measurements. Due to the different quality levels of the DEM, the accuracy of lake level reconstruction depends here also on the lake level itself.
  • The shoreline subset that is used for the retrieval of the lake level needs to be very shallow, so that the shift of the water-land border with changing lake level is maximized. Using RapidEye images, a decimeter accuracy of lake level reconstruction is only feasible if the shoreline slope is less than 3%.
These requirements are very strict and only very few lakes can fulfill the demands. Suitable lakes with shallow shorelines are, for example, Lake Wittwesee near Reihnsberg, Lake Redernswalder See near Parlow-Glambeck, Lake Briesensee near Poritz, and Lake Stechlin near Neuglobsow. The above-mentioned suitable lake sites are characterized by having little or no vegetation on sandy beaches during low level periods in the time range of the RapidEye acquisitions. However, at high lake levels reed vegetation often colonizes the shallow shorelines.
Another limitation on regular lake level monitoring is the actual temporal resolution of the RapidEye time series, with often less than one image per month and no suitable images from September/October until April. Frequent clouds coverage as well as cloud shadows reduced the temporal resolution of the RapidEye time series massively. Additionally, low solar angles during wintertime result in larger shadows and a worse signal-to-noise-ratio so that land is mistaken for water in the automatic classification. As a consequence, images with a very low SEA <30° and images with a low SEA between 30° and 45° during high lake levels needed to be removed. Thus, at LS only 75% of the RapidEye images in the archive could be finally used for the lake level reconstruction.
Nevertheless, monitoring lake level changes remains a hot topic as only ca. 500 of the ca. 5000 lakes in northeastern Germany have regular gauging stations [59,60] and because the RapidEye time series analysis in this study clearly showed the significance of lake level and shoreline changes for ecologic and economic reasons: Within one season, between April 2009 and September 2009, the size of the sandbank decreased by 90% and within one year the sandbank and shallow sandy beaches were completely flooded. This rapid change has high ecological impact on LF. For tourism the vanishing old public beach (subset B) with its wide sandy beach was a loss of attraction and at the new public beach (subset A) a swimming pier needed to be relocated. Also important was the flooding of the land connection of the peninsula in the southwest and the drowning of trees (subset C). The increase in lake level also influenced the submerged and aquatic vegetation: whereas the drying out of the shallow areas at low lake levels led to a significant reduction in the amount of pondweed and stonewort [1]; the shallow water areas were again colonized with reed in 2014. It is hoped that with a continuous high lake level the former underwater macrophytes will recolonize to stabilized and secure clear water conditions.

6. Conclusions

In this study, we measured lake levels indirectly by combining information on the water-land borders derived from high-resolution multi-spectral satellite images (RapidEye) and high-resolution topographic data. The results showed that the reconstruction of lake level changes worked well even for small and short-term lake levels changes. The RMSE between in situ measured lake level changes and reconstructed lake level changes was 6 cm. The absolute lake levels were reconstructed with decimeter accuracy, but they showed a systematic underestimation by −20–0 cm. The RMSE between in situ measured lake levels and reconstructed lake levels was 12 cm.
We tested two approaches for retrieval of water-land borders: automatic thresholds (Otsu) of the NIR bands and of the NDWI images. This study showed that NIR water-land borders were more precise than NDWI water-land borders, but also more sensitive to shadows. The RMSE between NIR water-land borders and according contour lines was 5.0 m (excluding winter images). Limitations in the accuracy of the delineation of the water-land border are caused by the characteristics of the RapidEye sensors and by the interfering shoreline vegetation. Mixed water and land pixels lead to an underestimation of the NIR water-land border. Due to changing shoreline conditions at different dates and lake levels, the accuracy of the lake level delineation can change over the monitoring period. The sensitivity of the NIR band to shadows causes severe misclassifications in the NIR water-land borders, especially in winter with low SEA angles. Clouds and shadows reduce significantly the repeat interval of usable RapidEye data acquisitions. Possible alternatives to tackle the limitations of the RapidEye optical system could be very high resolution data or Synthetic Aperture Radar (SAR). However, both alternatives bring new challenges. The use of very high resolution data (e.g., Worldview-2, Quickbird, Spot) for higher precision of water-land border determination comes with high cost of the images. SAR sensors are sunlight independent and able to penetrate clouds and light rain [61]. However, the side-looking geometry of the sensors is a disadvantage for spatially precise monitoring of shorelines [5].
The time series analysis of the NIR water-land borders at LF showed significant shifts of the shoreline with increasing lake levels. This is essential for the unambiguous reconstruction of different lake levels. The shift needs to be several times larger than the pixel size. For RapidEye images, at least a 30 m (6 pixels) shift between minimum and maximum lake levels (80 cm level difference) is necessary for reconstructing lake levels with decimeter accuracy. Such a 30 m shift results from shoreline slopes <4%. A smaller shift significantly reduces the accuracy of lake level reconstruction, so only a trend in lake level changes could be estimated.
Our results show that the reconstruction of small lake level changes can only be achieved using high resolution topographic data. For better accuracy, LIDAR acquisitions at low lake levels could substitute medium-resolution bathymetric surveys in shallow parts of the lake. As long as erosion and sedimentation rates are low, one-time acquisition of additional topographic data would be sufficient.
For the natural lakes in northeastern Germany and Poland, the transferability of the approach is limited, because of the specific requirements of the method. The very shallow shoreline slopes of the natural lakes often exhibit dense shoreline vegetation. Our experience indicates that the approach can be transferred to lakes with shorelines affected by human activity which keeps the shores free of vegetation such as public beaches or watering places. Additionally, the use of the lakes as reservoirs and fish ponds often introduces significant lake level changes.

Supplementary Materials

Supplementary materials can be found at: https://www.mdpi.com/2073-4441/7/8/4175/s1.

Acknowledgments

We would like to thank Julian Oeser and Robert Behling for their help in processing the RapidEye images and Christian Hohmann for his help in the field campaign. We are also grateful to our anonymous reviewers for their helpful and constructive comments.
This study was funded by the “Helmholtz Association of German Research Centres Initiative—Networking Fund for funding a Helmholtz Virtual Institute” (VH-VI-415). Satellite imagery was provided by ESA (proposal no. 14611) and BlackBridge (RESA Projekt-ID A1274).

Author Contributions

Iris Heine developed the methodological framework, performed programming, and conducted the analysis. Birgit Kleinschmit and Sibylle Itzerott were involved in formulating the research questions, preparing the manuscript, and contributing to critical discussions. Iris Heine, Birgit Kleinschmit and Sibylle Itzerott were involved in the general paper review. Peter Stüve supported the study by providing in situ data as well as expert knowledge about the gauging and the study area.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Kaiser, K.; Friedrich, J.; Oldorff, S.; Germer, S.; Mauersberger, R.; Natkhin, M.; Hupfer, M.; Pingel, A.; Schönfelder, J.; Spicher, V.; et al. Aktuelle hydrologische Veränderungen von Seen in Nordostdeutschland: Wasserspiegeltrends, ökologische Konsequenzen, Handlungsmöglichkeiten. In Wasserbezogene Anpassungsmaßnahmen an den Landschafts- und Klimawandel; Grünewald, U., Bens, O., Fischer, H., Hüttl, R.F., Kaiser, K., Eds.; Schweizerbart Science Publishers: Stuttgart, Germany, 2012; pp. 148–170. [Google Scholar]
  2. Mauersberger, R. Klassifikation der Seen für die Naturraumerkundung des nordostdeutschen Tieflandes. Arch. Naturschutz Landschaftsforsch. 2006, 3, 51–90. [Google Scholar]
  3. Germer, S.; Kaiser, K.; Bens, O.; Hüttl, R.F. Water Balance Changes and Responses of Ecosystems and Society in the Berlin-Brandenburg Region—A Review. DIE ERDE J. Geogr. Soc. Berlin 2011, 142, 65–95. [Google Scholar]
  4. Kaiser, K.; Koch, P.J.; Mauersberger, R.; Stüve, P.; Dreibrodt, J.; Bens, O. Detection and attribution of lake-level dynamics in north-eastern central Europe in recent decades. Reg. Environ. Chang. 2014, 14, 1587–1600. [Google Scholar] [CrossRef]
  5. Heine, I.; Francke, T.; Rogass, C.; Medeiros, P.H.A.; Bronstert, A.; Foerster, S. Monitoring Seasonal Changes in the Water Surface Areas of Reservoirs Using TerraSAR-X Time Series Data in Semiarid Northeastern Brazil. J. Sel. Top. Appl. Earth Obs. Remote Sens. 2014, 7, 3190–3199. [Google Scholar] [CrossRef]
  6. Van de Weyer, K.; Päzolt, J.; Tigges, P.; Raape, C.; Oldorff, S. Flächenbilanzierungen submerser Pflanzenbestände—Dargestellt am Beispiel des Großen Stechlinsees (Brandenburg) im Zeitraum von 1962–2008. Naturschutz Landschaftspfl. Brand. 2009, 18, 1–16. [Google Scholar]
  7. Brauns, M.; Garcia, X.F.; Pusch, M.T. Potential effects of water-level fluctuations on littoral invertebrates in lowland lakes. Hydrobiologia 2008, 613, 5–12. [Google Scholar] [CrossRef]
  8. Landesamt Brandenburg. Ökologische Charakterisierung der wichtigsten Brutgebiete für Wasservögel in Brandenburg. Stud. Tagungsberichte 2008, 57, 1–181. [Google Scholar]
  9. Landesamt Brandenburg (Ed.) Leitfaden zur Renaturierung von Feuchtgebieten in Brandenburg. Studien und Tagungsberichte des Landesumweltamtes, 2004; No. 50. Available online: http://www.lugv.brandenburg.de/cms/media.php/lbm1.a.3310.de/lua_bd50.pdf (accessed on 29 July 2015).
  10. Schmieder, K.; Dienst, M.; Ostendorp, W.; Joehnk, K. Effects of water level variations on the dynamics of the reed belts of Lake Constance. Ecohydrol. Hydrobiol. 2004, 4, 469–480. [Google Scholar]
  11. Hilt, S.; Henschke, I.; Rücker, J.; Nixdorf, B. Can submerged macrophytes influence turbidity and trophic state in deep lakes? Suggestions from a case study. J. Environ. Qual. 2010, 39, 725–733. [Google Scholar] [CrossRef] [PubMed]
  12. Hupfer, M.; Nixdorf, B. Zustand und Entwicklung von Seen in Berlin und Brandenburg. Materialien der Interdisziplinären Arbeitsgruppen (IAG Globaler Wandel—Regionale Entwicklung); Berlin-Brandenburgische Akademie der Wissenschaften: Leipzig, Germany, 2011; Diskussionspapier, No. 11. [Google Scholar]
  13. Ulrich, K.U. Vergleichende Untersuchungen zur Auswirkungen des Sediments auf die Wasserbeschaffenheit in Trinkwassertalsperren unter Berücksichtigung von Stauspiegelschwankungen; Cuvillier Verlag: Göttingen, Germany, 1998. [Google Scholar]
  14. Verpoorter, C.; Kutser, T.; Tranvik, L. Automated mapping of water bodies using Landsat multispectral data. Limnol. Oceanogr. Methods 2012, 10, 1037–1050. [Google Scholar] [CrossRef]
  15. Muster, S.; Heim, B.; Abnizova, A.; Boike, J. Water Body Distributions Across Scales: A Remote Sensing Based Comparison of Three Arctic TundraWetlands. Remote Sens. 2013, 5, 1498–1523. [Google Scholar] [CrossRef] [Green Version]
  16. Maillard, P.; Pivari, M.O.; Luis, C.H.P. Remote Sensing for Mapping and Monitoring Wetlands and Small Lakes in Southeast Brazil. In Remote Sensing of Planet Earth; Chemin, Y., Ed.; InTech: Rijeka, Croatia, 2012; pp. 23–46. [Google Scholar]
  17. Wang, J.; Sheng, Y.; Tong, T.S.D. Monitoring decadal lake dynamics across the Yangtze Basin downstream of Three Gorges Dam. Remote Sens. Environ. 2014, 152, 251–269. [Google Scholar] [CrossRef]
  18. Alesheikh, A.; Ghorbanali, A.; Nouri, N. Coastline change detection using remote sensing. Int. J. Environ. Sci. Technol. 2007, 4, 61–66. [Google Scholar] [CrossRef]
  19. White, K.; El Asmar, H.M. Monitoring changing position of coastlines using Thematic Mapper imagery, an example from the Nile Delta. Geomorphology 1999, 29, 93–105. [Google Scholar] [CrossRef]
  20. Lipakis, M.; Chrysoulakis, N. Shoreline extraction using satellite imagery. In Beach Erosion Monitoring: Results from BEACHMED-e/OpTIMAL Project; Pranzini, E., Wetzel, L., Eds.; Nuova Grafica Fiorentina: Florence, Italy, 2008; pp. 81–96. [Google Scholar]
  21. Li, X.; Damen, M.C.J. Coastline change detection with satellite remote sensing for environmental management of the Pearl River Estuary, China. J. Mar. Syst. 2010, 82, S54–S61. [Google Scholar] [CrossRef]
  22. Gens, R. Remote sensing of coastlines: Detection, extraction and monitoring. Int. J. Remote Sens. 2010, 31, 1819–1836. [Google Scholar] [CrossRef]
  23. Pardo-Pascual, J.E.; Almonacid-Caballer, J.; Ruiz, L.A.; Palomar-Vázquez, J. Automatic extraction of shorelines from Landsat TM and ETM+ multi-temporal images with subpixel precision. Remote Sens. Environ. 2012, 123, 1–11. [Google Scholar] [CrossRef]
  24. Baup, F.; Frappart, F.; Maubant, J. Combining high-resolution satellite images and altimetry to estimate the volume of small lakes. Hydrol. Earth Syst. Sci. 2014, 18, 2007–2020. [Google Scholar] [CrossRef]
  25. Maillard, P.; Bercher, N.; Calmant, S. New processing approaches on the retrieval of water levels in Envisat and SARAL radar altimetry over rivers: A case study of the São Francisco River, Brazil. Remote Sens. Environ. 2015, 156, 226–241. [Google Scholar] [CrossRef]
  26. Abarca-Del-Rio, R.; CrÉtaux, J.F.; Berge-Nguyen, M.; Maisongrande, P. Does Lake Titicaca still control the Lake Poopó system water levels? An investigation using satellite altimetry and MODIS data (2000–2009). Remote Sens. Lett. 2012, 3, 707–714. [Google Scholar] [CrossRef]
  27. Da Silva, J.S.; Seyler, F.; Calmant, S.; Rotunno Filho, O.C.; Roux, E.; Araújo, A.A.M.; Guyot, J.L. Water level dynamics of Amazon wetlands at the watershed scale by satellite altimetry. Int. J. Remote Sens. 2012, 33, 3323–3353. [Google Scholar] [CrossRef]
  28. Gupta, R.; Banerji, S. Monitoring of reservoir volume using LANDSAT data. J. Hydrol. 1985, 77, 159–170. [Google Scholar] [CrossRef]
  29. Smith, L.C. Satellite remote sensing of river inundation area, stage, and discharge: A review. Hydrol. Process. 1997, 11, 1427–1439. [Google Scholar] [CrossRef]
  30. Brakenridge, G.R.; Nghiem, S.V.; Anderson, E.; Chien, S. Space-based measurement of river runoff. Eos Trans. Am. Geophys. Union 2005, 86, 185. [Google Scholar] [CrossRef]
  31. Alsdorf, D.E.; Rodriguez, E.; Lettenmaier, D.P. Measuring surface water from space. Rev. Geophys. 2007, 45, 1–24. [Google Scholar] [CrossRef]
  32. Hostache, R.; Matgen, P.; Schumann, G.; Member, S.; Puech, C.; Hoffmann, L.; Pfister, L. Model Calibration Uncertainties Using Satellite SAR Images of Floods. Trans. Geosci. Remote Sens. 2009, 47, 431–441. [Google Scholar]
  33. Kaiser, K.; Heinrich, I.; Heine, I.; Natkhin, M.; Dannowski, R.; Lischeid, G.; Schneider, T.; Henkel, J.; Küster, M.; Heussner, K.; et al. Multi-decadal lake-level dynamics in north-eastern Germany as derived by a combination of gauging, proxy-data and modelling. J. Hydrol. 2015, in press. [Google Scholar] [CrossRef]
  34. Germer, S.; Kaiser, K.; Mauersberger, R. Sinkende Seespiegel in Nordostdeutschland: Vielzahl hydrologischer Spezialfälle oder Gruppen von Ähnlichen Seesystemen? In Aktuelle Probleme im Wasserhaushalt von Nordostdeutschland: Trends, Ursachen, Lösungen; Scientific Technical Report 10/10; Deutsches GeoForschungsZentrum: Potsdam, Germany, 2010. [Google Scholar]
  35. Data from Deutscher Wetterdienst. Available online: http://www.dwd.de/ (accessed on 10 March 2015).
  36. Landesamt für Umwelt Naturschutz und Geologie Mecklenburg-Vorpommern. Gutachtlicher Landschaftsrahmenplan Mecklenburgische Seenplatte. 2011. Available online: http://www.lung.mv-regierung.de/dateien/glrp_ms_06_2011.pdf (accessed on 29 July 2015).
  37. Ministerium für Landwirtschaft, Umwelt und Verbraucherschutz Mecklenburg-Vorpommern. Available online: http://www.regierung-mv.de/cms2/Regierungsportal_prod/Regierungsportal/de/lm/index.jsp (accessed on 4 December 2014).
  38. Staatliches Amt für Landwirtschaft und Umwelt Mecklenburgische Seenplatte (MS). Available online: http://www.stalu-mv.de/cms2/StALU_prod/StALU/de/ms/index.jsp (accessed on 8 January 2015).
  39. Sandau, R. Status and trends of small satellite missions for Earth observation. Acta Astronaut. 2010, 66, 1–12. [Google Scholar] [CrossRef]
  40. Landesamt für innere Verwaltung Mecklenburg Vorpommern. Available online: http://www.laiv-mv.de/land-mv/LAiV_prod/LAiV/AfGVK/index.jsp (accessed on 3 May 2013).
  41. RapidEye—Satellite Imagery Product Specifications, Version 6.1, April 2015. Available online: http://blackbridge.com/rapideye/upload/RE_Product_Specifications_ENG.pdf (accessed on 27 July 2015).
  42. Behling, R.; Roessner, S.; Segl, K.; Kleinschmit, B.; Kaufmann, H. Robust automated image co-registration of optical multi-sensor time series data: Database generation for multi-temporal landslide detection. Remote Sens. 2014, 6, 2572–2600. [Google Scholar] [CrossRef]
  43. Amt für Geoinformation, V.K. Satellitenpositionierungsdienst-SAPOS®. Available online: http://www.laiv-mv.de/land-mv/LAiV_prod/LAiV/AfGVK/Festpunkte,_SAPOS/SAPOS/index.jsp (accessed on 27 July 2015).
  44. Arbor, A.; Gilmer, D.S.; Fish, U.S.; Service, W. Utilization of Satellite Data for Inventorying Prairie Ponds and Lakes. Photogramm. Eng. Remote Sens. 1976, 42, 685–694. [Google Scholar]
  45. Kropáček, J.; Braun, A.; Kang, S.; Feng, C.; Ye, Q.; Hochschild, V. Analysis of lake level changes in Nam Co in central Tibet utilizing synergistic satellite altimetry and optical imagery. Int. J. Appl. Earth Obs. Geoinf. 2012, 17, 3–11. [Google Scholar] [CrossRef]
  46. Roessler, S.; Wolf, P.; Schneider, T. Multispectral Remote Sensing of Invasive Aquatic Plants Using RapidEye. In Earth Observation of Global Changes (EOGC); Krisp, J.M., Meng, L., Pail, R., Stilla, U., Eds.; Lecture Notes in Geoinformation and Cartography; Springer: Berlin, Germany, 2013; pp. 109–123. [Google Scholar]
  47. Ji, L.; Zhang, L.; Wylie, B. Analysis of Dynamic Thresholds for the Normalized Difference Water Index. Photogramm. Eng. Remote Sens. 2009, 75, 1307–1317. [Google Scholar] [CrossRef]
  48. Lu, S.; Ouyang, N.; Wu, B.; Wei, Y.; Tesemma, Z. Lake water volume calculation with time series remote-sensing images. Int. J. Remote Sens. 2013, 34, 7962–7973. [Google Scholar] [CrossRef]
  49. Li, W.; Du, Z.; Ling, F.; Zhou, D.; Wang, H.; Gui, Y.; Sun, B.; Zhang, X. A Comparison of Land Surface Water Mapping Using the Normalized Difference Water Index from TM, ETM+ and ALI. Remote Sens. 2013, 5, 5530–5549. [Google Scholar] [CrossRef]
  50. McFeeters, S.K. The use of the Normalized Difference Water Index (NDWI) in the delineation of open water features. Int. J. Remote Sens. 1996, 17, 1425–1432. [Google Scholar] [CrossRef]
  51. Lu, S.; Wu, B.; Yan, N.; Wang, H. Water body mapping method with HJ-1A/B satellite imagery. Int. J. Appl. Earth Obs. Geoinf. 2011, 13, 428–434. [Google Scholar] [CrossRef]
  52. Otsu, N. A threshold selection method from Gray-level. IEEE Trans. Syst. Man Cybern. 1979, 9, 62–66. [Google Scholar]
  53. Puech, C.; Raclot, D. Using geographical information systems and aerial photographs to determine water levels during floods. Hydrol. Process. 2002, 16, 1593–1602. [Google Scholar] [CrossRef]
  54. Hostache, R.; Matgen, P.; Schumann, G.; Puech, C.; Hoffmann, L.; Pfiste, L. Water level estimation and reduction of hydraulic model calibration uncertainties using satellite SAR images of floods. Trans. Geosci. Remote Sens. 2009, 47, 431–441. [Google Scholar] [CrossRef]
  55. Bochow, M.; Heim, B.; Küster, T.; Rogaß, C.; Bartsch, I.; Segl, K.; Reigber, S.; Kaufmann, H. On the Use of Airborne Imaging Spectroscopy Data for the Automatic Detection and Delineation of Surface Water Bodies. In Remote Sensing of Planet Earth; Chemin, Y., Ed.; InTech: Rijeka, Croatia, 2012; pp. 3–22. [Google Scholar]
  56. Lira, J. Segmentation and morphology of open water bodies from multispectral images. Int. J. Remote Sens. 2006, 27, 4015–4038. [Google Scholar] [CrossRef]
  57. Qiao, C.; Luo, J.; Sheng, Y.; Shen, Z.; Zhu, Z.; Ming, D. An Adaptive Water Extraction Method from Remote Sensing Image Based on NDWI. J. Indian Soc. Remote Sens. 2012, 40, 421–433. [Google Scholar] [CrossRef]
  58. Davaasuren, N.; Meesters, H.W.G. Extent and Health of Mangroves Mangroves in Lac Bay Bonaire Using Satellite Data; IMARES Wageningen UR: Haringkade, The Netherland, 2012. [Google Scholar]
  59. Päzolt, J.; Landesamt für Umwelt, Gesundheit und Verbraucherschutz, Brandenburg, Germany. Personal Communication, 2015.
  60. Mathes, J.; Ministerium für Landwirtschaft, Umwelt und Geologie, Mecklenburg-Vorpommers, Germany. Personal Communication, 2015.
  61. Lillesand, T.M.; Kiefer, R.W.; Chipman, J.W. Remote Sensing and Image Interpretation, 6th ed.; Wiley: Hoboken, NJ, USA, 2008. [Google Scholar]

Share and Cite

MDPI and ACS Style

Heine, I.; Stüve, P.; Kleinschmit, B.; Itzerott, S. Reconstruction of Lake Level Changes of Groundwater-Fed Lakes in Northeastern Germany Using RapidEye Time Series. Water 2015, 7, 4175-4199. https://doi.org/10.3390/w7084175

AMA Style

Heine I, Stüve P, Kleinschmit B, Itzerott S. Reconstruction of Lake Level Changes of Groundwater-Fed Lakes in Northeastern Germany Using RapidEye Time Series. Water. 2015; 7(8):4175-4199. https://doi.org/10.3390/w7084175

Chicago/Turabian Style

Heine, Iris, Peter Stüve, Birgit Kleinschmit, and Sibylle Itzerott. 2015. "Reconstruction of Lake Level Changes of Groundwater-Fed Lakes in Northeastern Germany Using RapidEye Time Series" Water 7, no. 8: 4175-4199. https://doi.org/10.3390/w7084175

APA Style

Heine, I., Stüve, P., Kleinschmit, B., & Itzerott, S. (2015). Reconstruction of Lake Level Changes of Groundwater-Fed Lakes in Northeastern Germany Using RapidEye Time Series. Water, 7(8), 4175-4199. https://doi.org/10.3390/w7084175

Article Metrics

Back to TopTop