Next Article in Journal
Detection of Urban Damage Using Remote Sensing and Machine Learning Algorithms: Revisiting the 2010 Haiti Earthquake
Next Article in Special Issue
Hyperspectral Sensors as a Management Tool to Prevent the Invasion of the Exotic Cordgrass Spartina densiflora in the Doñana Wetlands
Previous Article in Journal
Using Google Earth Surface Metrics to Predict Plant Species Richness in a Complex Landscape
Previous Article in Special Issue
Estimating the Total Nitrogen Concentration of Reed Canopy with Hyperspectral Measurements Considering a Non-Uniform Vertical Nitrogen Distribution
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Effect of Protection Level in the Hydroperiod of Water Bodies on Doñana’s Aeolian Sands

1
Remote Sensing and GIS Laboratory (LAST-EBD), Estación Biológica de Doñana, CSIC, C/Américo Vespucio s/n, Isla de la Cartuja, 41092 Sevilla, Spain
2
Department of Wetland Ecology, Estación Biológica de Doñana, CSIC, C/Américo Vespucio s/n, Isla de la Cartuja, 41092 Sevilla, Spain
*
Author to whom correspondence should be addressed.
Submission received: 28 July 2016 / Revised: 11 October 2016 / Accepted: 13 October 2016 / Published: 20 October 2016
(This article belongs to the Special Issue What can Remote Sensing Do for the Conservation of Wetlands?)

Abstract

:
Mediterranean temporary ponds on Doñana’s aeolian sands form an extensive system of small dynamic water bodies, dependent on precipitation and groundwater, of considerable importance for biodiversity conservation. Different areas of the aeolian sands have received different levels of environmental protection since 1969, and this has influenced the degree of conservation and the flooding dynamic of these temporary surface waters. We use the Landsat series of satellite images from 1985 to 2014 to study the temporal dynamic of small temporary water bodies on the aeolian sands in relation to the protection level and to distance to water abstraction pressures from agriculture and residential areas. The results show that even with small and ephemeral water bodies optical remote sensing time-series are an effective way to study their flooding temporal dynamics. The protected areas of the aeolian sands hold a better preserved system of temporary ponds, with a flooding dynamic that fluctuates with precipitation. The unprotected area shows an increase in mean hydroperiod duration, and surface flooded, and a decline in hydroperiod variability. This seems to be due to the creation of irrigation ponds and the artificialization of the flooding regime of the natural temporary ponds, that either receive excess irrigation water or dry-up due to the lowering of the groundwater table level. Although a decline in hydroperiod duration of temporary ponds is seen as negative to the system, an increase in hydroperiod of surface waters due to artificialization, or a decline in variability cannot be considered as positive compensatory effects.

Graphical Abstract

1. Introduction

The Doñana National Park holds a system of ponds on aeolian sands that is considered one of the best examples of Mediterranean temporary ponds in Europe [1,2]. The aeolian sands extend out of the protected area and with it also the system of temporary ponds. Ponds start flooding with autumn and winter rains and dry-up during spring and summer [3]. Groundwater seepage feeds the ponds, and the oscillation of the water table level is also responsible for ponds flooding and drying up [1,3,4]. Ponds differ in size and also in the length of time they get flooded every year. They go from rain puddles that remain flooded for a few days -or weeks on rainy winters- to some ponds, up to 187 ha in size, that keep permanent water in part of their basin most of the years. On the other hand, there is a large variability among years on when each pond gets flooded and how long it keeps the water [1]. We denominate hydroperiod the duration of the flooding period of a pond in an annual cycle [5]. This is an important ecological parameter as it determines which plants and animals can use the pond for breeding and to complete their life cycles [2]. Some species have short breeding cycles and require ponds with short hydroperiods where they can out-compete other species with longer cycles. In general, temporary ponds, that dry up completely at some time in their annual cycle do not allow the colonization of species that require permanent water, like fish, and are comparatively free of predators [2].
The aeolian sands system rests on top of a single, partially confined, aquifer known as the “Almonte-Marismas” aquifer [6]. It is considered a single aquifer from the perspective of groundwater dynamic although it has been recently divided in 5 sections from a water management point of view [7]. The Doñana’s system of Mediterranean temporary ponds on the aeolian sands has received substantial attention from the perspective of its hydrology [8,9], geomorphology [10,11,12], relation to groundwater [4,13,14], water chemistry [15,16], biodiversity [17,18,19,20] but most studies have been conducted in a few ponds and only for a few years.
Remote sensing is a cost-effective technique to study the temporal dynamic of wetlands. It has the capacity to study simultaneously numerous water bodies over a large area, has the potential to scale-up in-situ data, and can provide retrospective information on wetland evolution using time series of medium resolution sensors [21,22,23,24]. The Landsat satellites have been acquiring images with comparable sensors from 1972 to the present [25]. In the case of Doñana’s aeolian sands more than 400 Landsat satellite images have been acquired during the last 35 years.
Optical images are an adequate tool to discriminate flooded areas and can delineate continental water bodies taking advantage of the low reflectance of water in the Near Infrared (NIR) and Short Wave Infrared (SWIR) regions [26,27]. They have been used to monitor inland water characteristics [28,29], but its use is less frequent for shallow and small waters bodies [30,31]. In particular, turbid waters or waters covered by aquatic vegetation may be more difficult to map due to the reflectance peaks between 700 and 900 nm related to the presence of phytoplankton and/or suspended sediments [29,32], and the reflectance of vegetation in the NIR region. So small, shallow temporary ponds sometimes turbid or covered by vegetation can be challenging to map accurately [33].
Optical sensors on board of Landsat satellites have a medium resolution, 60 m pixels in the MSS and 30 m pixels in TM, ETM+ and OLI multispectral bands. They overpass the area every 16 days and, as it is possible to have a long time series of comparable images, we think it is possible to study the hydrological behavior and temporal trends of relatively small temporary ponds on the Doñana’s aeolian sands. The medium resolution of Landsat sensors do not allow to distinguish by their shape small artificial ponds that have been built for irrigation or recreational purposes from natural temporary ponds. Also, some natural ponds in the area may have been transformed and be suffering an artificial flooding regime. For that reason we use the term “temporary water body” to refer to all areas identified by remote sensing that show some kind of periodic flooding regime on the aeolian sands, and we limit the use of the term “temporary pond” to those that have a natural flooding regime.
There are areas in the Doñana’s aeolian sands that have received different levels of legal protection. We expected that the areas that have had a higher level of protection and have been protected for a longer time will have a better preserved system of temporary ponds with a temporal dynamic that is less affected by water use, and more dependent on precipitation. As water use for irrigated agriculture and water consumption by residential areas are the main pressures for groundwater abstraction or surface water diversion in the area [9,13,14,34], we expected that ponds close to irrigated agriculture or residential areas would be the ones more affected in their hydroperiods and would show a sharper declining trend with time. As agriculture and people “compete” for water with the ponds, we expected that ponds that were close to irrigated agriculture and residential areas would have shorter hydroperiods. Also, as water demand by agriculture and people increase, the mean hydroperiod and total surface flooded would show a stronger declining trend with time in areas with a higher surface of irrigated agriculture or residential areas.

2. Study Area

The study area is the Doñana’s aeolian sands, also known as the coastal aeolian mantle (“Manto Eólico Litoral de Doñana”) in the province of Huelva, Spain (Figure 1). The area is delimited by the aeolian sand deposits as extracted from the REDIAM geomorphological map of Andalusia 1:400,000 [35]. These are coastal aeolian sand deposits of Quaternary age on top of Miocene marine marls. Part of the aeolian sands are included in Doñana National Park that was created in 1969 with an extension of 507.2 km2. The protected area was expanded in 1980 with the creation of a buffer zone the “Preparque Doñana” (542.5 km2), that was later transformed into the Doñana Natural Park by the Andalusian regional government in 1989. UNESCO recognized Doñana National Park as a Biosphere Reserve with an extent of 772.6 km2 in 1980. In 1982 Doñana wetlands were included in the list of the Ramsar Convention. In July 2012, UNESCO approved the extension of the biosphere reserve to the present limits covering an area of 2550 km2. The aeolian sands are affected by these different levels of protection. An area of 239 km2 has been included within the National Park since 1969. A total of 255 km2 were included in the buffer zone in 1980 and, currently, are part of the Doñana Natural Park. An extent of 316 km2 were included in the 2012 expansion of Doñana Biosphere Reserve. The rest of the aeolian sands (108 km2) have no legal protection, with the exception of some water-table lakes that are included in the “Paraje Natural Lagunas de Palos y Las Madres” (1989), and have an extent of 6.9 km2 (Figure 1). We consider in this paper four sections of the aeolian sands that have been receiving different levels of legal protection considering both the time since they were legally protected and the different legal protection forms. They go, from more to less protection: National Park > Natural Park > Biosphere reserve > Unprotected area.
The main transformations in the area that affect water use are residential areas and agriculture. In 2007 residential areas were occupying 8.8 km2, outside the National and Natural Park, but extended close to their limits. Irrigated agriculture had been increasing during the last three decades, mainly in the area that is unprotected but also within the Doñana Biosphere Reserve, and was covering an extent of 101 km2 [36]. The “Almonte-Marismas” aquifer extends over an area of 3290 km2 and includes 88.6 km2 of residential use and 617 km2 of irrigated agriculture (Figure 2).

3. Material and Methods

3.1. Satellite Imagery

We acquired every available Landsat TM, and ETM+ image for the Doñana’s aeolian sands (Path: 202, Row: 34) from 1985 to 2014. We decided to exclude the images from the Landsat MSS sensor because most temporary ponds in the aeolian sands are shallow and small in size, the MSS pixel size (79 × 57 m2) is too coarse to identify them, MSS has no SWIR band that is the most adequate for identification of shallow waters, and from 1975 to 1984 the number of satellite images per year was small (0–8 images) so hydroperiod estimation of temporary ponds had too much error. We used 325 valid satellite images from 1985 to 2014 providing 30 years of continuous observation using the TM and ETM+ sensor from Landsat 4 to 7. Some images were discarded due to acquisition problems such as missing lines and columns, radiometric incoherence or lines shifts, and cloud cover. Number of valid scenes per year ranged from 4 to 22 with a median of 11 scenes (Table S1).
In order to make the time series comparable, a semi-automatic robust and coherent pre-processing was applied. The complete procedure is detailed by Díaz-Delgado et al. [37] but in short it included metadata retrieval from raw 1 G data images, geometric and atmospheric corrections followed by time series radiometric normalization using ordinary least squares regression towards a reference cloud-free and clear atmosphere image (18 July 2002) using in pseudo-invariant areas [38,39].

3.2. Hydroperiod Estimation

We extracted the SWIR band, b5 (1.550–1.750 in TM μm) of the TM and ETM+ images, and obtained a flood mask by considering as flooded any pixel with a normalized reflectance <0.186 (Kappa = 0.66, global accuracy = 94% [37]). This procedure has been successfully tested to classify flooded areas in shallow waters affected by turbidity and aquatic vegetation [29,30,37], and was employed by Gómez-Rodríguez et al. [33] to map small temporary ponds in the same area.
We call pixel annual hydroperiod the number of days a 30 m pixel remains flooded in an annual hydrological cycle. As rains in the area start after the northern hemisphere summer our cycles go from 1 September to 31 August of the next year. By convention we name the cycles with the end year, so the cycle from 1 September 1984 to 31 August 1985 is the “1985 cycle”. To calculate the pixel annual hydroperiod we reclassified all flooded pixels in image i to the number of days spanned from the previous image (I − 1), provided the pixel was flooded in the previous image, and added up all the reclassified flood masks in the cycle. So pixel annual hydroperiod (H) can be computed according to Equation (1).
Hc = i = 1 n 1 ( DoC i + 1 DoC i ) = ( DoC 2 DoC 1 ) + + ( DoC n DoC n 1 )
where DoC stands for Day of Cycle. For every annual cycle, the i-th image corresponds to the first date the pixel was classified as flooded while the n image is the last date when the pixel was classified as flooded. The n value may eventually be equal the total number of flooding masks per cycle if the pixel is classified as flooded at every flooding mask for that annual cycle. The procedure assumes that a pixel classified as flooded between two consecutive dates has not dried-up in between. If a pixel is dry at an intermediate date the formula discounts the days lapsed from the previous flooded image.
This method to calculate the pixel annual hydroperiod tends to underestimate the true value, because the pixel should flood before the first image in the cycle we classify it as flooded, and should dry-up after the last image in the cycle we classify it as flooded. As this underestimation is influenced by the number and the dates of valid images in a cycle, this tends to create an underestimation in the initial cycles that have a smaller number of images. The existence of this bias is clear by calculating the temporal trend of hydroperiod (H) in pixels in the Atlantic Ocean that show a positive, although not significant, temporal trend (slope b = 0.87, Student-t = 0.875, p = 0.4). To correct this bias we performed a linear stretch to the hydroperiod estimate multiplying Hi by 365/Hmax. Being Hi the value of the pixel, Hmax the maximum value of H in the annual cycle that correspond to pixels extracted from the ocean, and 365 days the expected maximum value of hydroperiod in the cycle. Correction factors for each cycle (C) are provided in Table S1.

3.3. Study Area Limits

We created a mask of the Doñana’s aeolian sands by reclassifying all the geomorphological units on the Geological map of Andalusia [35] on top of the “Almonte-Marismas” aquifer that are classified as active dunes, vegetated dunes, aeolian sand mantles, and sand deposits, including the small geomorphological units in-between classified as creeks or lakes, but excluding the impermeable clay deposits that form the Doñana and the Tinto and Odiel marshes (Figure 1).
Protected area limits were taken from the official published limits of the Doñana Natural Space that include the National and the Natural Park and the Doñana Biosphere Reserve as of 2016. The “Rocina protection area” of Doñana National Park is considered with the same level of protection as the Doñana National Park.
The official limits of the “Almonte-Marismas” aquifer were extracted from the web map services from the Ministry of Agriculture and Environment from Spain [40] considering all the aquifer sections below the Doñana’s aeolian sands.

3.4. Potential Temporary Waters and Temporary Water Bodies Definitions

Any pixel on the aeolian sands that was classified as flooded once during the study period is considered as potentially holding temporary water ponds, and we denominate this as potential temporary waters. Temporary ponds can be small in size and the spatial resolution of TM and ETM+ (30 m pixels) suggests that water bodies bellow 900 m2 in area can be overlooked and will be omission errors [17,33]. The procedure is also prone to errors of commission caused by cloud, topographic, and tree shadows, specially during winter when the sun is low. All pixels that were classified as potential temporary waters in any annual cycle were overlaid to create a mask of total potential temporary water that was used to define the maximum extent of potential temporary waters on the aeolian sands (Figure 1) to estimate mean pixel hydroperiod values, temporal trends, and number of temporary water bodies.
We consider as a single temporary water body the one formed by all pixels from the maximum extent of temporary waters that are spatially connected and separated from other water bodies. Pixels that are connected are probably not independent water bodies and analyzing the data on a pixel by pixel basis would be pseudoreplication. This definition of water body we use implies that ponds that are considered as independent ponds in other studies may have been considered as a single water body in our analyses if they are connected at high water levels. This, together with the lower size limit of 900 m2 to identify a temporary water body, should be taken into account when using the number of water bodies estimated with our methods to compare them with those provided by other studies in the area [41]. Initial vectorization of the potential temporary waters mask without correcting for topographic shadows or pine woods identified a total of 3638 water bodies.

3.5. Correction of Potential Temporary Waters for Topographic Shadows and Pine Woods

Preliminary exploration of the data suggested that topographic shadows and dense pine woods could be the cause of commission errors in the maximum extent of temporary waters.
To correct for topographic shadows we created a mask of mean topographic shadows in December, when solar elevation is lowest. We calculated a topographic hill-shading image on 21 December, at the time of Landsat overpass, using the 5-m DEM from CNIG [42] and resampled it at 30 m (with ArcGis). We considered as shadows all pixels with relative illumination below 80 (less that 2% of the lowest values of the hill-shading image). Shadows did not seem to be a general a source of commission errors as could be observed overlaying a vector file of potential temporary waters pixels on top of the hill-shading image. To evaluate the extent of commission errors caused by shadows we extracted the pixels under shadows and only a few temporary water bodies had 100% of the pixels in shadows (33 out of 3638 temporary water bodies, less than 1%).
To correct for errors caused by dense pine woods we reclassified two Landsat images, one from December 2013 and one from December 2010, using a density slice of NDVI values. We selected a time of the year when pine woods show maximum NDVI values [43]. The cut-point in NDVI was chosen in each case to select dense pine woods, avoiding selecting shrubland or pastures, using as reference well known areas in the Doñana Biological Reserve. A potential temporary waters vector file was overlaid on the Landsat pine maps. Some large areas of pine woods in the south of Doñana tended to be classified as potential temporary waters. It has to be taken into account that if a pixel covered by pines is identified as water in a single mask it will appear in this maximum extent of temporary waters. Out of 3638 water bodies a total of 173 had 100% of pixels covered by pines in both dates, 69 were covered by pines in 2010, and another 176 were covered by pines in 2013. So a total of 418 (11%) temporary water bodies could be commission errors caused by pine woods.
To reduce commission errors caused by pine woods, and by topographic shadows, we excluded from the maximum extent of temporary waters all pixels covered by pine woods or affected by shadows using these masks.
After masking the coverage was vectorized to obtain a total 3939 temporary water bodies. After this we corrected the coverage by eliminating water bodies that corresponded to tidal zones (variations in tide level on the coast) and those that were connected to the marsh, to the river, or to the sea. The final number of temporary water bodies on the aeolian sands was reduced to 3678.

3.6. Environmental Predictors

3.6.1. Precipitation

We obtained the mean accumulated precipitation on top of the “Almonte-Marismas” aquifer for each annual cycle (1 September to 31 August) summing monthly precipitation surfaces obtained from REDIAM [44]. Surfaces have been calculated by REDIAM by interpolating from all available automatic meteorological station using the inverse distance algorithm in ArcGis. From this gis raster cover we extracted the mean precipitation for each hydroperiod cycle from September 1984 to August 2012 and calculated the mean for the aquifer polygon using SAGA 2.1.4 Statistics for Grids [45]. The data from 1984 to 2009 are from the “serie consolidada” [44] and from 2010 to 2012 from “serie provisional” [46]. To complete the series we obtained the monthly precipitation data from 1 September to 31 August from the I.F.A.P.A. automatic agroclimatic stations [47] inside the aquifer polygon for the cycles 2013 and 2014 that had no interpolated data by REDIAM.

3.6.2. Distance to Irrigated Agriculture and Residential Areas

We assume that distance to irrigated agriculture and distance to residential areas are the best predictors for the main sources for groundwater abstraction and surface water diversion for agriculture and human consumption respectively. Not all wells used for groundwater pumping are known, volumes pumped are not well documented, and many wells are illegal. The distribution of irrigated agriculture and residential areas has been done by photo-interpretation and so less prone to error. We used as source for irrigated agriculture and residential areas the 2007 1:25,000 Land use-land cover coverage from Andalusia from REDIAM [36]. The area was clipped with a polygon defined calculating a 5 km buffer on the “Almonte-Marismas” aquifer limits (Figure 2). We reclassified all polygons representative of irrigated agriculture and calculated distance in meters to the closest polygon with IDRISI Selva module DISTANCE [48]. In a similar way we reclassified all land uses representative of residential areas, including gardens and golf courses and calculated distance in m to the closest residential area pixel. Reclassification criteria are provided in Tables S2 and S3. Pixels included in irrigated agriculture polygons have a distance value of 0 m to irrigated agriculture and pixels included in residential polygons have a distance value of 0 m to residential areas.

3.6.3. Altitude above Sea Level

Altitude above sea level can be an important characteristics of a temporary pond as it can determine the elevation in relation to the groundwater table and in relation to the surface drainage network, and explain differences in flooding behavior. We used the 5 m DEM from GNIG [42] resampled to 30 m to calculate the mean altitude above sea level for each water body.

3.6.4. Spatial Location of Water Bodies

Trends in hydroperiod or in the correlation of hydroperiod with precipitation could be due to other spatial trends that are not explained by the environmental predictors we selected (e.g., distance to irrigation cultures or to residential areas). We calculated the X and Y coordinates of the centroid for each temporary water body and tested for a significant spatial trend in the data. Spatial trends can be tested, or corrected, in statistical models in several ways. For example, testing for a linear or polynomial trend in X or in Y. In our models we used smoothing splines surfaces, that is a method of fitting a smooth surface to noisy data using splines, that are piece-wise polynomial functions. The degrees of freedom (d.f.) of the spline surface indicate how much flexibility the model is allowed to follow the spatial pattern in the data. We tested in each model a smoothing spline surface of the spatial coordinates with up to 3 d.f. [49,50]. This is a way of testing if the trends observed in relation to distance to threats or protection level are still significant when correcting for more general spatial trends.

3.7. Photo-Interpretation of Randomly Selected Temporary Water Bodies

Not all water bodies identified with this procedure are temporary ponds. To estimate the proportion of different types of water bodies in our coverage we photo-interpreted a random selection of water bodies. We selected at random 100 water bodies from each of the four levels of protection (National Park, Natural Park, Biosphere reserve and Unprotected area). A random number was generated for each water body and the first 100 from each level were selected for photo-interpretation. One of us (DA) blindly photo-interpreted the water body vector limit using Bing and Google Earth images and classified them in the following categories: (1) natural pond; (2) stream; (3) irrigation pond; (4) artificial pond; (5) dense pine wood; (6) other dense vegetation; (7) disperse vegetation; (8) non vegetated, and also classified the water body neighborhood as: (1) natural vegetation; (2) cereal agriculture; (3) irrigated agriculture; (4) greenhouses; (5) industrial; (6) urban; (7) beach-tidal; and (8) other. The criteria used by the photo-interpreter are provided in Tables S4 and S5.

3.8. Statistical Analyses

We calculated the mean hydroperiod for each pixel included in the coverage of maximum extent of temporary waters on the aeolian sands (excluding topographic shadows, pine woods and tidal areas) for the whole time series (1985–2014) by adding the annual hydroperiod images and dividing by the number of years with IDRISI Selva Image Calculator [48].
With the module PROFILE of IDRISI Selva [48] we extracted a temporal profile of the maximum surface flooded, mean hydroperiod, and standard deviation of the hydroperiod for each annual cycle for pixels in each of the four protection levels, using only the pixels included in the mask of maximum extent of temporary waters.
The annual hydroperiod of temporary ponds in the aeolian sands is heavily dependent on accumulated precipitation [33]. To reduce the noise in pixel hydroperiod temporal trend caused by variable precipitation in each annual cycle we calculated for each pixel an ordinary least squares regression using as response the annual hydroperiod and as predictor the mean accumulated precipitation. From this regression we calculated R-squared at the pixel level, and obtained the residual images of the linear regression using IDRISI Selva module CORRELATE [48]. With the residual images we calculated the Mann-Kendall rank correlation (or Kendall’s tau) and Theil-Sen slope with IDRISI Selva module KENDALL. Positive Kendall’s tau indicate pixels that have positive trends in hydroperiod once corrected for precipitation (water bodies that did not exist previously or that are increasing in hydroperiod more than it could be expected from precipitation alone) while negative values indicate declining temporal trends in hydroperiod (that cannot be explained by a decline in precipitation). The Theil-Sen slope is a robust method for fitting a linear trend to a set of points that chooses the median slope among all lines through pairs of two points. It is much less sensitive to outliers and in our case indicates how fast is hydroperiod declining or increasing once corrected for precipitation.
All statistical models were fitted in R (version 3.2.5) [51]. We used linear models and generalized additive models (GAMs), fitted with stats and gam 1.12 R packages respectively, to test which environmental predictors could explain the relationship between precipitation and hydroperiod, and the temporal trends in hydroperiod using all temporary water bodies. For each water body and year we extracted the mean value of the R-squared of the regression with precipitation, mean Mann-Kendall rank correlation (Kendall’s tau), and mean Theil-Sen slope, using all pixels included in the water body. These were used as response variables in the models. The R-squared value (range 0–1) was arcsine transformed, calculating the arcsine of the square root, as this transformation normalized the distribution of errors. For each water body we extracted the protection level (National Park, Natural Park, Biosphere Reserve, Unprotected area), mean distance to irrigated agriculture, mean distance to residential areas, mean altitude above sea level, and X and Y coordinates of the centroid (in m), to be used as predictors. Distances were log transformed (distance + 1).
We started fitting a linear model that included as predictors: protection level as a factor, and distance to irrigated agriculture and distance to residential areas as a linear terms using normal errors. This is the model defined by our a priori hypothesis, that water abstraction pressures and protection level are the main factors influencing hydroperiod behavior in temporary ponds. The model was tested in stepwise fashion using step.gam (gam 1.12 package) for simplification or for increased complexity (backwards-forward procedure) by removing terms or by substitution of the linear term with 2 d.f. and 3 d.f. smoothing splines of the distance to irrigated agriculture or distance to residential areas, a linear term of the altitude above sea level, or a smoothing spline surface of the spatial coordinates (2 d.f. and 3 d.f). The best model obtained in this procedure, the one with minimum AIC, was selected. Δ AIC values for each predictor in each model are given in Tables S6–S9.
We expected that water bodies could show a different temporal trend depending on their natural or artificial origin. To test if natural or artificial water bodies could be showing different temporal trends in hydroperiod we fitted a full linear model with protection level, type of flooding (natural or artificial), distance to irrigated agriculture and distance to residential areas with all interaction terms, and performed a backward simplification until all remaining terms and interactions were significant according an F-test. Then we refitted the models using a random sample of water bodies that had been photo-interpreted as with natural flooding (natural ponds, streams, and open vegetation far from agricultural or residential areas) to test if the models found using all data were representative of natural temporary ponds or were mainly caused by artificial water bodies.

3.9. Structure of Statistical Analyses

First, we analyze, within each level of protection, the trends from 1985 to 2014 in: mean hydroperiod, coefficient of variation of hydroperiod, and maximum surface flooded, for the pixels that were included in the maximum extent of temporary waters (Figure 1), and we compare these trends with the trend in precipitation. Then we correlate precipitation with mean hydroperiod and surface flooded for each level of protection to test to what extent precipitation explains the differences in trends.
Second, we fit models at the level of water body to test if distances to threats (irrigation and residential areas) and protection level statistically explain the differences in correlation between hydroperiod an precipitation, and the differences in hydroperiod trend once corrected by precipitation.
Third, from all water bodies we select a random sample from each protection level that are photo-interpreted and classified into natural water bodies, artificial water bodies and errors. They allow us to calculate commission errors, and to estimate proportions of natural ponds and artificial ponds in each protection level.
Fourth, with the sample of photo-interpreted water bodies, using as response variable hydroperiod trend (Kendall’s tau) we test for an interaction between flooding regime (natural or artificial) and protection level, under the hypothesis that natural ponds may be showing different trends in hydroperiod depending on the protection level of the area.
Fifth, we refit models to the natural ponds identified in the random sample to test to what extent correlation models developed for all water bodies are representative of the hydroperiod behavior of natural temporary ponds.

4. Results

4.1. Density and Distribution of Temporary Water Bodies in Relation to Protection Level

We identified a total of 3678 temporary water bodies ranging in size from 900 m2 to 1,879,200 m2 and covering a total surface of 29.6 km2. In Table 1 we provide some statistics of the distribution and density of water bodies and maximum surface of temporary waters in relation to the protection level in the aeolian sands.

4.2. Temporal Trends in Hydroperiod in Relation to Protection Level

We analyzed the trends in the mean and coefficient of variation of hydroperiod of all pixels of temporary waters within a certain protection level. From 1985 to 2014 the mean pixel hydroperiod of temporary waters showed a non statistically significant declining temporal trend in the areas with higher protection levels, both in the National Park (high protection, slope b = −0.084, F = 0.06, p = 0.8, Figure 3A) and in the Natural Park (medium protection, b = −0.23, F = 0.43, p = 0.52, Figure 3A). Meanwhile, the trend in mean pixel hydroperiod showed a significant increase in the Biosphere Reserve (low protection, b = 0.31, F = 4.63, p = 0.04, Figure 3B) and a highly significant increase in the unprotected area (unprotected, b = 1.65, F = 40.65, p < 0.001, Figure 3C). During the same period, the mean accumulated annual precipitation on top of the aquifer showed a non significant increasing trend (b = 2.74, F = 0.41, p = 0.53, Figure 3D).
The coefficient of variation of hydroperiod of temporary water pixels within each protection level has maintained a stable temporal trend in the protected areas (high protection, b = 1.14, F = 0.18, p = 0.67, Figure 4A; medium protection, b = 2.42, F = 0.67, p = 0.46, Figure 4B; low protection, b = −3.88, F = 2.27, p = 0.14, Figure 4C); but a highly significant declining trend out of the protected areas (unprotected, b = −2.15, F = 28.33, p < 0.001, Figure 4D).
The maximum flooded surface per year at each protection level shows a similar behavior to that of the mean hydroperiod. The only significant trend is an increasing trend out of the protected areas (unprotected, b = 4.0, F = 13.55, p < 0.001, Figure S1).
Annual accumulated precipitation on top of the “Almonte-Marismas” aquifer explains a large amount of the variation in mean pixel hydroperiod in the areas with higher levels of protection, the National and the Natural Park (high protection: adjusted R-squared = 34.16%, b = 0.047, F = 16.04, p < 0.001, Figure 5A; medium protection: adjusted R-squared = 27.2%, b = 0.047, F = 11.84, p = 0.002, Figure 5B); but no significant amount of mean pixel hydroperiod variation in the Biosphere Reserve, or out of the protected area (low protection: adjusted R-squared = 0.88%, b = 0.007, F = 1.26, p = 0.27, Figure 5C; unprotected: adjusted R-squared = 0.87%, b = 0.015, F = 0.75, p = 0.39, Figure 5D).
Annual accumulated precipitation also explained to a great extent the annual maximum flooded surface in the areas with higher levels of protection, the National and the Natural Park (high protection: adjusted R-squared = 47.08%, b = 0.913, F = 26.8, p < 0.0001; medium protection: adjusted R-squared = 37.15%, b = 0.584, F = 18.14, p = 0.0002, respectively); but no significant amount in the Biosphere Reserve, or out of the protected area (Figure S2).

4.3. Correlation between Hydroperiod and Precipitation

At each 30 m pixel of the maximum extent of temporary waters we calculated an ordinary least squares regression between annual hydroperiod and mean annual accumulated precipitation on the aquifer. We calculated the R-squared value for each temporary water body as the mean R-squared of all the pixels that constitute the water body. The R-squared value indicates if a water body has an hydroperiod very dependent on annual precipitation (high values) or not correlated at all (low values). If we make a density plot to see how frequent are water bodies of each type we see multi modal distributions that differ in shape among the protection levels. Water bodies with low correlation with precipitation are more frequent in the areas with low levels of protection (Biosphere Reserve and unprotected, Figure 6).

4.4. Temporal Trends in Hydroperiod of Water Bodies

The Kendall’s tau indicates if a pixel has a positive or negative trend in hydroperiod with year, once, corrected for differences in annual precipitation. They are extracted for each water body as the mean value of all the pixels that constitute the water body. A map of the spatial distribution of Kendall’s tau values for water bodies suggest that water bodies with positive trends (blue circles) concentrate in the unprotected area while the more negative trends (red circles) show up in the Natural park (Figure 7). A density plot of mean Kendall’s tau of water bodies also shows multi modal distributions that differ among protection levels. Protected areas show a dominance of water bodies with negative trends in hydroperiod, while in the unprotected area water bodies with positive trends are more abundant (Figure S3).

4.5. Predictive Model of Hydroperiod Correlation with Precipitation for All Water Bodies

The best model for the correlation between hydroperiod and precipitation (R-squared) included distance to irrigated cultures, distance to residential areas, altitude, protection level, and the spatial coordinates (Figure 8). The shape of the partial responses to distance to irrigated cultures indicates that water bodies very close to irrigated cultures show almost no influence of precipitation in hydroperiod, influence increases steeply in the first 1–2 km from irrigation and then declines and stabilizes at an average value (Figure 8A). Close to residential areas water bodies show a mean low correlation between hydroperiod and precipitation, but with wide confidence intervals, that indicate great variability in responses (Figure 8B). The relationship with altitude indicates that water bodies higher above sea level, farther from the groundwater table, are more correlated with precipitation, while those at lower altitude, closer to the groundwater table, fluctuate less with precipitation (Figure 8C). The relationship with protection level indicates that areas with higher protection levels (National and Natural Park) have water bodies with hydroperiods that are more dependent on precipitation while those in the Biosphere Reserve reserve, or out of the protected area, fluctuate less with precipitation (Figure 8D). The effect of spatial coordinates indicates that there is some spatial variability that is not adequately explained by previous factors (Model significance table in Table S6).

4.6. Predictive Model of Hydroperiod Trend for All Water Bodies

We used the residuals of the ordinary least squares regression of annual hydroperiod with accumulated precipitation to calculate a non-parametric correlation (Kendall’s tau) with year. Positive values of the Kendall’s tau indicate that the hydroperiod of the water body, once corrected for annual fluctuations in precipitation, has increased with time, and negative correlations that it has decreased with time. The best model for the Kendall’s tau obtained by stepwise selection included as predictors: distance to irrigated cultures, protection level, altitude, and the spatial coordinates (Figure 9). The shape of the model partial response with distance to irrigated cultures indicates that water bodies very close to irrigated cultures have increased their hydroperiod, and that this effect steeply declines with distance to irrigated agriculture (Figure 9A). The effect with altitude indicates that water bodies at lower altitudes, those more dependent on the groundwater table, have declined more with time than those at higher altitudes that are farther from the water-table, and show a more variable trend (wider confidence interval, Figure 9B). The effect of the protection level shows a peculiar pattern. Hydroperiods in general have declined more with time with decreasing protection, but in the unprotected area the hydroperiods have increased (Figure 9C). The spatial coordinates indicates that there is spatial variation that is unexplained by environmental predictors (Model significance table in Table S7).
The same procedure was used for the Theil-Sen slope that is the median slope of hydroperiod with year, and indicates the sign and speed of change in hydroperiod. The Theil-Sen slope of each water body is correlated with Kendall’s tau, but values are not normally distributed having long tails of positive and negative slopes. The best statistical model for the Theil-Sen slope obtained by stepwise selection of predictors included: protection level, distance to irrigated cultures, distance to residential areas, and spatial coordinates. This model is not independent from that of Kendall’s tau, and results are similar. Details are presented in Supplementary Material.

4.7. Photo-Interpretation of Temporary Water Bodies in Relation to Protection Level

A random selection of 100 temporary water bodies for each protection level was photo-interpreted and classified into different water bodies types (Table 2). This classification showed that temporary water bodies were dominated by areas with bare soil and disperse or herbaceous vegetation in which a clear pond basin could not be identified (62%). This is coherent with a majority of the water bodies corresponding to rain puddles or ephemeral water bodies that last only for a few days or weeks after heavy rainfalls, and perhaps are not formed every year. They were followed by 20% natural temporary ponds (ponds with a clearly defined basin), 6% streams, 10% irrigation and artificial ponds, and 3.5% of dense pine woods (probably commission errors not solved by our masking process).
This classification was simplified into two types of flooding regimes, (1) natural flooding: (natural ponds and water-table lakes, streams, disperse and non vegetated areas surrounded by natural vegetation or dryland agriculture); (2) artificial flooding (irrigation ponds, artificial ponds, disperse and non vegetated areas surrounded by irrigated agriculture, greenhouses, industrial or urban); and (3) classification errors (pine woods, or bare soil from beach or tidal areas that cannot be considered temporary ponds) (Table 3). This new classification indicated a 4.75% of commission errors in our procedure to identify temporary pools. It also shows that water bodies of artificial origin represent a different proportion of water bodies depending on the protection level. They are 42% of water bodies in the unprotected area, 21% in the Biosphere Reserve, 2% in the Natural Park and are absent (or below 1%) in the National Park.

4.8. Interaction between Natural or Artificial Flooding Regime and Protection Level

Using the water bodies that were photo-interpreted we tested for a significant interaction between protection level and type of flooding (natural or artificial) in relation to hydroperiod trend (Kendall’s tau). The model (Table 4) indicated that there were no significant linear effects with distance to irrigated cultures or distance to residential areas, but a significant interaction between protection level and type of flooding (natural or artificial). Artificial water bodies have increased their hydroperiod in areas with low levels of protection (Biosphere reserve, or unprotected area). Natural water bodies show decreasing trends in hydroperiod with decreasing protection level (the decline is greater in the Biosphere reserve compared to the Natural and National Park) but an opposite increasing trend in the unprotected area (Figure 10).

4.9. Predictive Model of Hydroperiod Correlation with Precipitation for Natural Water Bodies

As the correlation of hydroperiod with precipitation and the trend of hydroperiod with time could be affected by the presence of artificial water bodies which could show a different behavior to that of natural water bodies, and differ in their proportion in the different protection levels, we repeated the stepwise model selection using only the random sample of water bodies that had been photo-interpreted as with natural flooding (n = 316).
The best model for the correlation between hydroperiod and precipitation (R-squared) included distance to irrigated cultures, distance to residential areas, protection level, and spatial coordinates (Figure 11). The shape of the partial effect of the model indicates that correlation with precipitation is low close to irrigated cultures (Figure 11A), or residential areas (0 distance), but has wide confidence intervals (large variability in behavior, Figure 11B). Then, correlation declines with distance to irrigated cultures but is not affected by distance to residential areas. The relationship with protection level indicates a higher correlation with precipitation in areas with higher levels of protection (National and Natural Park) than with lower levels of protection (Biosphere Reserve or unprotected) (Figure 11C).

4.10. Predictive Model of Hydroperiod Trend for Natural Water Bodies

The best model for hydroperiod temporal trend (Kendall’s tau) included as predictors: distance to residential areas, protection level, and spatial coordinates (Figure 12). The model with distance to residential areas indicates a sharp temporal decline in hydroperiod duration close to residential areas but with wide confidence intervals (Figure 12A). The relationship with protection level indicates a clear decline in hydroperiod in the Biosphere Reserve but an increase in the unprotected area, compared to the National and Natural Park (Figure 12B). The results for the Theil-Sen slope were similar to those of Kendall’s tau. This model, as well as the statistical significance tables for R-squared and Kendall’s tau, are provided in Supplementary Material.

5. Discussion

Mediterranean temporary ponds are a priority habitat in the European Union (Code 3170, Habitats Directive 92/43/ECC). They are of great conservation importance because they harbor a large part of the aquatic biodiversity at the landscape scale, and are home to rare and endangered species [1,19,52]. Climate change, together with impacts by urbanization and agriculture are threatening these habitats at the global scale [52]. The system of ponds on Doñana’s aeolian sands is a good example of Mediterranean temporary ponds. They have been declining with climate change by an increased aridification since the “Little Ice Age” [10,11,12,53,54] but have also been suffering an increased anthropogenic pressure [1,55]. It has been documented how groundwater abstraction for agriculture [13] or for tourism resorts has negatively influenced or dessicated certain ponds in the system [34].
In this paper we show that the time series of Landsat satellite images is a great source of information for surface water variability and temporal trends. It can be used to study temporal trends in the system of temporary ponds on Doñana’s aeolian sands, its hydrological evolution, and to monitor the landscape effect of different conservation actions. Even in a system where many temporary ponds can be smaller than a Landsat pixel [41] the high number of images allows to detects trends and evaluate changes in the flooding regime. Remote sensing has the advantage of being able to study the system as a whole providing a synoptic view, and the length of the Landsat time series allows to study the hydrological dynamic for a period of 30 years.
Of course, our study has some limitations. The 30 m pixel size of Landsat precludes a clear identification of many ponds that are below this pixel size, also the mean interval between images, 30 days, is longer of the annual hydroperiod of many ephemeral ponds, so it is difficult to obtain accurate estimates of annual hydroperiod duration for this kind of ponds. With our definition of water body, many water bodies are sets of smaller ponds that at our resolution appear as joined. In the Doñana National park another study by Gomez-Rodríguez et al. [41] with higher spatial resolution images (5 m pixels) located twice as many temporary ponds as we did. Also, although the methods used to identify flooded areas have been validated in other papers [27,37] and we have performed a validation of our method to identify temporary water bodies by photo-interpreting a random sample, this validation is only partial. We cannot retrospectively test if a single pixel that we classified as temporarily flooded, or dry, has had water, or not, coincident with a Landsat overpass in the 30 years period. Our mean pixel hydroperiod values (Table 1) are calculated including areas that flood only occasionally, and are values much lower than those provided by other works that estimate hydroperiod duration from initial flooding to complete dry-up of each pond [33,34,41].
Even though, we think our methods provide a new synoptic approach to the study of temporary ponds on the Doñana’s aeolian sands, point out that surface water availability is created not only by natural ponds but also by a diverse gradient from natural to artificial water bodies. All this suggest that the hydrodynamic behavior detected by remote sensing could be also used to identify and monitor the natural or artificial flooding regime of water bodies in the area.
The main results of our study show that there is a higher density, extension, and variability of temporary water bodies in the areas of the Doñana’s aeolian sands that have received higher levels of environmental protection (Table 1). As our analyses with Landsat images start in 1985, once the area was already protected, we cannot discriminate if it was that the area with a higher richness and diversity of temporary ponds was selected for protection, or if legal protection has prevented a stronger decline in the National park compared to the Natural park or the Biosphere reserve (that is the dryer one, with a lower extension and density of temporary water bodies).
The Doñana protected areas (National and Natural park) have temporary water bodies with wider fluctuations in hydroperiod and more dependent on annual precipitation. Protected areas show a slight decline in mean hydroperiod while precipitation has been slightly increasing during the last 30 years, although neither trends are statistically significant (Figure 3A,D).
Hydroperiod of temporary water bodies has increased in the areas with low or no protection (Biosphere reserve and unprotected area Figure 3C,D), that are the areas receiving the greater impact of irrigated agriculture and residential use (Figure 2), and there is no correlation between precipitation and mean hydroperiod (Figure 5C,D) or annual flooded surface in these areas.
The coefficient of variation of the hydroperiod has significantly declined in the unprotected areas (water bodies showing less variability in hydrological regime, Figure 4D) and the surface flooded has increased also (Figure S1).
A random selection and photo-interpretation of water bodies suggest low commission errors (only less than 4% water bodies could be errors) and validates partially our approach to study this system of small temporary water bodies on the aeolian sands. It also shows that we are detecting different types of water bodies both according to their origin (natural ponds, artificial irrigation ponds, ponds transformed into water reservoirs for agriculture) and to their hydrodynamics (rain puddles, epigenetic and hypogenetic temporary ponds, water-table lakes, etc.). The photo-interpretation indicates that the contribution of the different typologies of water bodies to the sample is variable depending on protection level. In the aeolian sands out of the protected areas 22% of the water bodies that hold surface water are artificial irrigation ponds or ponds that have been modified for irrigation (pond that have been artificialized and show irrigation pumps and pipes). In the Biosphere reserve artificial ponds are 14% of the water bodies (Table 2). These artificial water bodies have been increasing during the last 30 years, and, at least partially, explain the increase in surface flooded and in mean hydroperiod, and the decline in the coefficient of variation of the hydroperiod out of the protected area. We checked the 20 irrigation ponds identified in the random sample (Table 2) in historical aerial photos and no one was visible in 1984 (one was a natural pond at that time) while 13 years later, in 1997, 17 of the irrigation ponds (85%) had already been built. Artificial ponds are no substitute for natural temporary ponds, but provide alternative habitats for certain aquatic organisms and can act, in certain situations, as complementary habitats of conservation value [56]; but, in general, artificial water bodies have a reduced biodiversity value [57].
The models for correlation of hydroperiod with precipitation, and the temporal trends in hydroperiod once corrected for precipitation, indicate important effects of distance to water abstraction pressures (irrigated cultures and residential areas) and of protection level in the hydrological behavior of temporary water bodies. First of all, distance to irrigated cultures is an important driver but its overall effect is contrary to what we expected a priori. Water bodies close to irrigated cultures show an increasing trend in hydroperiod, mainly because of the construction of new irrigation ponds. But also because natural ponds outside the protected areas have a trend of increasing hydroperiod (Figure 10 and Figure 12B). This change has affected natural ponds in the unprotected area, but not in the Biosphere Reserve that show declining hydroperiods (Figure 10). One possible explanation is that the excess of water for irrigation finally flows into natural ponds close to the fields and these ponds increase their hydroperiods, and show lower fluctuations with annual precipitation. In fact, a transfer of 4.99 hm3 of surface water from other basins has been authorized for the irrigation cultures on the wester sector of the aeolian sands, in the unprotected area (Huelva province), to reduce the illegal use of groundwater [58,59,60,61]. Distance to residential areas has a smaller impact, but tends to affect negatively the trend in hydroperiod (Figure 12A). In any case, the wide confidence intervals close to residential areas and irrigated cultures, mainly in relation to the correlation between hydroperiod and precipitation point to the existence of different types of water bodies that can be influenced in different ways. New artificial water bodies are created in gardens and golf courses that are artificially flooded and do not depend on annual precipitation, while water-table lakes close to extraction wells suffer stronger fluctuations with the water table and show a stronger correlation with annual precipitation [34]. At the same time, other water-table lakes have hydroperiods that depend on the accumulated rainfall over several years [33], and have a smaller correlation with annual precipitation.

6. Conclusions

Our study shows that it is possible to study the long-term flooding dynamic of a system of small temporary ponds using the Landsat satellites time-series of optical images to evaluate the effect of protection level, and the threats caused by water abstraction by irrigation cultures and residential areas. The results indicate that the protected area holds a better preserved system of temporary ponds but that its flooding dynamic seems to be affected by water abstraction by irrigated cultures and, to a smaller extent, by residential areas. The effect of irrigation and residential areas is not simple at a landscape scale. Mean hydroperiod has increased by the creation of new irrigation ponds and recreation ponds in agricultural and residential areas, respectively. The flooding regime of some natural ponds is altered by groundwater table fluctuation generated by water abstraction, while other ponds are probably receiving excess irrigation water, or are directly transformed into irrigation reservoirs. Our study suggests that mean hydroperiod of water bodies over an extensive area is probably not a good indicator of the degree of conservation of temporary ponds, but low correlation with precipitation and a low coefficient of variation of hydroperiod could probably be indicators of the degree of artificialization of an area. At the same time, our study suggests that flooding dynamics and flooding anomalies (differences from mean values or mean trends), as determined by remote sensing data, could be used to classify and identify different types of water bodies and evaluate their degree of artificialization.

Supplementary Materials

The following are available online at www.mdpi.com/2072-4292/8/10/867/s1, Figure S1: Temporal trends in the maximum flooded surface of potential temporary waters within protection levels. (A) Red circles = National Park; (B) green triangles = Natural Park; (C) blue squares = Biosphere Reserve; (D) gray circles = unprotected. Only statistically significant trend lines are shown, Figure S2: Relationship between annual accumulated precipitation and the maximum flooded surface of potential temporary waters within protection levels. (A) Red circles = National Park; (B) green triangles = Natural Park; (C) blue squares = Biosphere Reserve; (C) gray circles = unprotected. Only statistically significant trend lines are shown, Figure S3: Density plot of the Kendall’s tau of water bodies. Positive values indicate increasing hydroperiods and negative values decreasing hydroperiods, once corrected for annual precipitation. Protection level indicated with colors: high protection = red, medium protection = green, low protection = blue, unprotected = gray. The rug plots below indicate the distribution of water bodies. High protection (red) and unprotected (gray) lines are drawn thicker for clarity, Figure S4: Density plot of the Theil-Sen slope of water bodies. Positive values indicate increasing hydroperiods and negative values decreasing hydroperiods, once corrected for annual precipitation. Protection level indicated with colors: high protection = red, medium protection = green, low protection = blue, unprotected = gray. The rug plots below indicate the distribution of water bodies. High protection (red) and unprotected (gray) lines are drawn thicker for clarity, Figure S5: Partial effect of predictors in the best gam model for annual hydroperiod trend corrected for precipitation (Theil-Sen slope) for all water bodies, Figure S6: Partial effects in the best gam model for temporal trend in annual hydroperiod corrected for precipitation (Theil-Sen slope) for natural water bodies, Figure S7: Distribution of water bodies according to their size. Size in m2. Sample size for each size interval indicated on the column, Figure S8: Example of hydroperiod distribution maps for a small area inside the National park for two cycles that have similar precipitation but a different number of Landsat images. (A) 1996; (B) 2010. Darker blue indicates longer hydroperiod for the pixel. The black line indicates the maximum extent of temporary waters (pixels that have been flooded at least once in 30 years). The distribution of images in the flooding cycle is shown, Table S1: Number of images used, distribution by months, and correction factor (C) for each cycle, Table S2: Reclassification of land use /land cover map as used to calculate distance to irrigated agriculture. First a clip of the area of the polygon defined by the “Almonte-Marismas” aquifer with a 5 km buffer was taken from the 2007 REDIAM MUCVA_25_07_EscDetalle (Land use and land cover vegetation map of Anadalusia 2007). Then in Qgis the following land uses were reclassified as irrigated agriculture. The reclassified shapefile is named W_regadio_ED50-PS, Table S3: Reclassification of land use /land cover map as used to calculate distance to residential areas. First a clip of the area of the polygon defined by the “Almonte-Marismas” aquifer with a 5 km buffer was taken from the 2007 REDIAM MUCVA_25_07_EscDetalle (Land use and land cover vegetation map of Anadalusia 2007). Then in Qgis the following land uses were reclassified as residential areas. The reclassified shapefile is named W_residencial_ED50-PS, Table S4: Criteria used in the photo-intrepretation of water bodies. The polygons selected were transformed into a klm file and overlayed in Google Earth and Bing using the mosrt recent images of the area. The interior part of the polygon was classified exclusively in one of the following categories, Table S5: Criteria used in the photo-intrepretation of water bodies. The surroundings of the polygon was classified in one of the following categories, Table S6: Relative importance of predictors in the best gam model for the correlation between annual hydroperiod and precipitation (R-squared) for all water bodies. Irrigation = smoothing spline of log distance to irrigated cultures, residential = smoothing spline of log distance to residential areas, protection = protection level (factor, four levels), altitude = altitude above sea level, latitude, longitude = surface smoothing spline of spatial coordinates in m. Predictors in decreasing order of relative importance based on Δ AIC estimated by single term simplification from the full model, Table S7: Relative importance of predictors in best gam model for annual hydroperiod trend corrected for precipitation (Kendall’s tau) for all water bodies. Irrigation = smoothing spline of log distance to irrigated cultures, protection = protection level (factor, four levels), altitude = altitude above sea level, latitude, longitude = surface smoothing spline of spatial coordinates in m. Predictors in decreasing order of relative importance based on Δ AIC estimated by single term simplification from the full model, Table S8: Relative importance of predictors in the best gam model for the correlation between annual hydroperiod and precipitation (R-squared) for a random selection of natural water bodies. Irrigation = smoothing spline of log distance to irrigated cultures, residential = smoothing spline of log distance to residential areas, protection = protection level (factor, four levels), latitude, longitude = surface smoothing spline of spatial coordinates in m. Predictors in decreasing order of relative importance based on Δ AIC estimated by single term simplification from the full model, Table S9: Relative importance of predictors in the best gam model for temporal trend in annual hydroperiod corrected for precipitation (Kendall’s tau) for natural water bodies. Residential = smoothing spline of log distance to residential areas, protection = protection level (factor, four levels), latitude, longitude = surface smoothing spline of spatial coordinates in m. Predictors in decreasing order of relative importance based on Δ AIC estimated by single term simplification from the full model, Table S10: Relative importance of predictors in best gam model for annual hydroperiod trend corrected for precipitation (Theil-Sen slope) for all water bodies. Irrigation = smooth spline of log distance to irrigated cultures, residential = smooth spline of log distance to residential areas, protection = protection level (factor, four levels), latitude, longitude = surface smooth spline of spatial coordinates in m. Predictors in decreasing order of relative importance based on Δ AIC estimated by single term simplification from the full model, Table S11: Relative importance of predictors in the best gam model for temporal trend in annual hydroperiod corrected for precipitation (Theil-Sen slope) for natural water bodies. protection = protection level (factor, four levels), latitude, longitude = surface smooth spline of spatial coordinates in m. Predictors in decreasing order of relative importance based on Δ AIC estimated by single term simplification from the full model.

Acknowledgments

This study has been funded by the Spanish Ministry of Science and Innovation through the research projects HYDRA (#CGL2006-02247/BOS) and HYDRA2 (CGL2009-09801/BOS), and by funding from the European Union’s Horizon 2020 research and innovation program under grant agreement No. 641762 to ECOPOTENTIAL project. We acknowledge comments by C. Díaz-Paniagua and R. Díaz-Delgado on previous drafts of the paper. We are grateful to many students who helped in image processing.

Author Contributions

All the authors are intellectually responsible for the conducted research. Javier Bustamante (J.B.) designed the work, run statistical analyses, and prepared the manuscript. David Aragonés (D.A.), Isabel Afán and J.B. processed and analyzed the time series of images and participated in ground-truth data collection. D.A. photo-interpreted a random sample of temporary water bodies. J.B. wrote a first draft of the manuscript and all authors participated in the final writing.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Díaz-Paniagua, C.; Fernández-Zamudio, R.; Serrano, L.; Florencio, M.; Gómez-Rodríguez, C.; Sousa, A.; Sńachez Castillo, P.; García Murillo, P.; Siljestrom, P. El Sistema de Lagunas Temporales de Doñana, una Red de Hábitats Acuáticos Singulares; Naturaleza y Parques Nacionales, Serie Técnica; Organismo Autónomo Parques Nacionales: Madrid, Spain, 2015. [Google Scholar]
  2. Díaz-Paniagua, C.; Zamudio, R.F.; Florencio, M.; Murillo, P.G.; Rodríguez, C.G.; Portheault, A.; Martín, L.S.; Siljestrom, P. Temporay ponds from Doñana National Park: A system of natural habitats for the preservation of aquatic flora and fauna. Limnetica 2010, 29, 41–58. [Google Scholar]
  3. Serrano, L.; Reina, M.; Martín, G.; Reyes, I.; Arechederra, A.; León, D.; Toja, J. The aquatic systems of Doñana (SW Spain): Watersheds and frontiers. Limnetica 2006, 25, 11–32. [Google Scholar]
  4. Sacks, L.A.; Herman, J.S.; Konikow, L.F.; Vela, A.L. Seasonal dynamics of groundwater-lake interactions at Doñana National Park, Spain. J. Hydrol. 1992, 136, 123–154. [Google Scholar] [CrossRef]
  5. Grillas, P.; Gauthier, P.; Yavercovski, N.; Perennou, C. Mediterranean Temporary Pools; Station biologique de la Tour du Valat: Arles, France, 2004. [Google Scholar]
  6. Llamas, M.R. Geohydrology of the eolian sands of the Doñana National Park (Spain). Catena Suppl. 1990, 18, 145–154. [Google Scholar]
  7. Plan Hidrológico de la DH del GUADALQUIVIR (2015–2021). Real Decreto 1/2016, de 8 de Enero, por el que se Aprueba la Revisión de los Planes Hidrológicos de las Demarcaciones Hidrográficas del Cantábrico Occidental, Guadalquivir, Ceuta, Melilla, Segura y Júcar, y de la parte Española de las Demarcaciones Hidrográficas del Cantábrico Oriental, Miño-Sil, Duero, Tajo, Guadiana y Ebro. Available online: https://www.boe.es/boe/dias/2016/01/19/pdfs/BOE-A-2016-439.pdf (accessed on 27 May 2016).
  8. Lozano, E.; Delgado, F.; Manzano, M.; Custodio, E.; Coleto, C. Hydrochemical characterisation of ground and surface waters in “the Cotos” area, Doñana National Park, southwestern Spain. Groundw. Hum. Dev. 2005, 6, 217–231. [Google Scholar]
  9. Manzano, M.; Custodio, E.; Lozano, E.; Higueras, H. Relationships between wetlands and the Doñana coastal aquifer (SW Spain). Groundw. Ecosyst. 2013, 169, 169–182. [Google Scholar]
  10. Sousa, A.; García-Barrón, L.; Morales, J.; García-Murillo, P. Post-Little Ice Age warming and desiccation of the continental wetlands of the aeolian sheet in the Huelva region (SW Spain). Limnetica 2006, 25, 57–70. [Google Scholar]
  11. Sousa, A.; García-Murillo, P.; Sahin, S.; Morales, J.; García-Barrón, L. Wetland place names as indicators of manifestations of recent climate change in SW Spain (Doñana Natural Park). Clim. Chang. 2010, 100, 525–557. [Google Scholar] [CrossRef]
  12. Sousa, A.; García-Murillo, P. Changes in the wetlands of Andalusia (Doñana Natural Park, SW Spain) at the end of the Little Ice Age. Clim. Chang. 2003, 58, 193–217. [Google Scholar] [CrossRef]
  13. Suso, J.; Llamas, M. Influence of groundwater development on the Doñana National Park ecosystems (Spain). J. Hydrol. 1993, 141, 239–269. [Google Scholar] [CrossRef]
  14. Guardiola-Albert, C.; Jackson, C.R. Potential impacts of climate change on groundwater supplies to the doñana Wetland, Spain. Wetlands 2011, 31, 907–920. [Google Scholar] [CrossRef] [Green Version]
  15. Serrano, L.; Toja, J. Limnological description of four temporary ponds in the Doñana National Park (SW, Spain). Arch. Hydrobiol. 1995, 133, 497–516. [Google Scholar]
  16. Serrano, L. Leaching from vegetation of soluble polyphenolic compounds, and their abundance in temporary ponds in the Doñana National Park (SW Spain). Hydrobiologia 1992, 229, 43–50. [Google Scholar] [CrossRef]
  17. Gómez-Rodríguez, C.; Bustamante, J.; Koponen, S.; Díaz-Paniagua, C. High-resolution remote-sensing data in amphibian studies: Identification of breeding sites and contribution to habitat models. Herpetol. J. 2008, 18, 103–113. [Google Scholar]
  18. Gómez-Rodríguez, C.; Díaz-Paniagua, C.; Bustamante, J.; Portheault, A.; Florencio, M. Inter-annual variability in amphibian assemblages: Implications for diversity assessment and conservation. Aquat. Conserv. Mar. Freshw. Ecosyst. 2010, 20, 668–677. [Google Scholar] [CrossRef]
  19. Florencio, M.; Díaz-Paniagua, C.; Gómez-Rodríguez, C.; Serrano, L. Biodiversity patterns in a macroinvertebrate community of a temporary pond network. Insect Conserv. Divers. 2014, 7, 4–21. [Google Scholar] [CrossRef]
  20. Florencio, M.; Díaz-Paniagua, C.; Serrano, L.; Bilton, D.T. Spatio-temporal nested patterns in macroinvertebrate assemblages across a pond network with a wide hydroperiod range. Oecologia 2011, 166, 469–483. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  21. Almeida, T.I.R.; Penatti, N.C.; Ferreira, L.G.; Arantes, A.E.; Amaral, C.H. Principal component analysis applied to a time series of MODIS images: The spatio-temporal variability of the Pantanal wetland, Brazil. Wetl. Ecol. Manag. 2015, 2015, 1–12. [Google Scholar] [CrossRef]
  22. Domínguez Gómez, J.A.; Chuvieco Salinero, E.; Sastre Merlín, A. Monitoring transparency in inland water bodies using multispectral images. Int. J. Remote Sens. 2009, 30, 1567. [Google Scholar] [CrossRef]
  23. Sokoletsky, L.G.; Lunetta, R.S.; Wetz, M.S.; Paerl, H.W. MERIS retrieval of water quality components in the turbid albemarle-pamlico sound estuary, USA. Remote Sens. 2011, 3, 684–707. [Google Scholar] [CrossRef]
  24. Yan, Y.E.; Ouyang, Z.T.; Guo, H.Q.; Jin, S.S.; Zhao, B. Detecting the spatiotemporal changes of tidal flood in the estuarine wetland by using MODIS time series data. J. Hydrol. 2010, 384, 156–163. [Google Scholar] [CrossRef]
  25. Röder, A.; Kuemmerle, T.; Hill, J. Extension of retrospective datasets using multiple sensors. An approach to radiometric intercalibration of Landsat TM and MSS data. Remote Sens. Environ. 2005, 95, 195–210. [Google Scholar] [CrossRef]
  26. Johnston, R.; Barson, M. Remote sensing of Australian wetlands: An evaluation of Landsat TM data for inventory and classification. Mar. Freshw. Res. 1993, 44, 235–252. [Google Scholar] [CrossRef]
  27. Bustamante, J.; Díaz-Delgado, R.; Aragonés, D. Determinación de las características de masas de aguas someras en las marismas de Doñana mediante teledetección. Rev. Teledetec. 2005, 24, 107–111. [Google Scholar]
  28. Belluco, E.; Camuffo, M.; Ferrari, S.; Modenese, L.; Silvestri, S.; Marani, A.; Marani, M. Mapping salt-marsh vegetation by multispectral and hyperspectral remote sensing. Remote Sens. Environ. 2006, 105, 54–67. [Google Scholar] [CrossRef]
  29. Bustamante, J.; Pacios, F.; Díaz-Delgado, R.; Aragonés, D. Predictive models of turbidity and water depth in the Doñana marshes using Landsat TM and ETM+ images. J. Environ. Manag. 2009, 90, 2219–2225. [Google Scholar] [CrossRef] [PubMed]
  30. Gardiner, N.; Díaz-Delgado, R. Trends in Selected Biomes, Habitats and Ecosystems: Inland Waters. In Sourcebook on Remote Sensing and Biodiversity Indicators; Secretariat of the Convention on Biological Diversity: Montreal, QC, Canada, 2007; pp. 83–102. [Google Scholar]
  31. Davranche, A.; Poulin, B.; Lefebvre, G. Mapping flooding regimes in Camargue wetlands using seasonal multispectral data. Remote Sens. Environ. 2013, 138, 165–171. [Google Scholar] [CrossRef] [Green Version]
  32. Binding, C.E.; Bowers, D.G.; Mitchelson-Jacob, E.G. Estimating suspended sediment concentrations from ocean colour measurements in moderately turbid waters; the impact of variable particle scattering properties. Remote Sens. Environ. 2005, 94, 373–383. [Google Scholar] [CrossRef]
  33. Gómez-Rodríguez, C.; Bustamante, J.; Díaz-Paniagua, C. Evidence of Hydroperiod Shortening in a Preserved System of Temporary Ponds. Remote Sens. 2010, 2, 1439–1462. [Google Scholar] [CrossRef] [Green Version]
  34. Díaz-Paniagua, C.; Aragonés, D. Permanent and temporary ponds in Doñana National Park (SW Spain) are threatened by desiccation. Limnetica 2015, 34, 407–424. [Google Scholar]
  35. Mapa Geomorfológico de Andalucía: Fisiografía, Unidades Geomorfológicas, Procesos Geomorfológicos y Cuencas Marinas. (Consejería de Medio Ambiente, Junta de Andalucía). Scale 1:400.000. Available online: http://www.juntadeandalucia.es/medioambiente/mapwms/REDIAM_Geologico_Andalucia? (accessed on 13 November 2015).
  36. Mapa de usos y coberturas vegetales del suelo de Andalucía. (Consejería de Medio Ambiente, Junta de Andalucía). Scale 1:25.000. Available online: http://www.juntadeandalucia.es/medioambiente/servtorrent/MUCVA25_2007.torrent (accessed on 13 November 2015).
  37. Díaz-Delgado, R.; Aragonés, D.; Afán, I.; Bustamante, J. Long-term monitoring of the flooding regime and hydroperiod of Doñana marshes with Landsat time series (1974–2014). Remote Sens. 2016. [Google Scholar] [CrossRef]
  38. Hall, F.G.; Strebel, D.E.; Nickeson, J.E.; Goetz, S.J. Radiometric rectification: Toward a common radiometric response among multidate, multisensor images. Remote Sens. Environ. 1991, 35, 11–27. [Google Scholar] [CrossRef]
  39. Schott, J.R.; Salvaggio, C.; Volchok, W.J. Radiometric scene normalization using pseudoinvariant features. Remote Sens. Environ. 1988, 26, 1–16. [Google Scholar] [CrossRef]
  40. Web Map Service: Agua—Planificación Hidrográfica—Masas de Agua Subterránea Artículo 5 Directiva Marco del Agua. Available online: http://wms.magrama.es/sig/agua/MasasAguaSubterraneaDMA2005/wms.aspx (accessed on 13 November 2015).
  41. Gómez-Rodríguez, C.; Díaz-Paniagua, C.; Bustamante, J. Cartografía de lagunas temporales del Parque Nacional de Doñana; Agencia Andaluza del Agua, Consejería de Medio Ambiente, Junta de Andalucía: Sevilla, Spain, 2011. [Google Scholar]
  42. DTM5: Digital elevation model, 5 m pixel, from the Plan Nacional de Ortofotografía Aérea (PNOA). Reference system ETRS89. Proyecction UTM. (Instituto Geográfico Nacional). Scale 1:5.000. Available online: https://www.ign.es/ign/layoutIn/modeloDigitalTerreno.do (accessed on 27 May 2016).
  43. Aragonés, D.; Díaz-Delgado, R.; Afán, I.; Bustamante, J.; Moreno, E. Estudio de la fenología de los pinares de pino piñonero en Andalucía occidental mediante el análisis de imágenes MODIS. In Biogeografía de Sistemas Litorales. Dinámica y Conservación; Cámara, R., Rodríguez, B., Muriel, J.L., Eds.; Universidad de Sevilla: Sevilla, Spain, 2014; pp. 417–421. [Google Scholar]
  44. Precipitación mensual en Andalucía (serie consolidada). Raster format (GeoTIFF). (Consejería de Medio Ambiente, Junta de Andalucía). Available online: http://www.juntadeandalucia.es/medioambiente/servtorrent/Precip_mensual.torrent (accessed on 13 November 2015).
  45. Conrad, O.; Bechtel, B.; Bock, M.; Dietrich, H.; Fischer, E.; Gerlitz, L.; Wehberg, J.; Wichmann, V.; Böhner, J. System for Automated Geoscientific Analyses (SAGA) v. 2.1.4. Geosci. Model. Dev. 2015, 8, 1991–2007. [Google Scholar] [CrossRef]
  46. Precipitación mensual en Andalucía. Raster format (GeoTIFF). (Consejería de Medio Ambiente, Junta de Andalucía). Available online: http://www.juntadeandalucia.es/medioambiente/servtorrent/Precip_mensual_prov.torrent (accessed on 13 November 2015).
  47. Estaciones Agroclimáticas Junta de Andalucía. (Consejería de Agricultura, Pesca y Desarrollo Rural, Junta de Andalucía). Available online: http://www.juntadeandalucia.es/agriculturaypesca/ifapa/ria/servlet/FrontController?action=Init (accessed on 10 May 2009).
  48. Eastman, J.R. IDRISI Selva Tutorial. Idrisi. Prod. Clark Labs-Clark Univ. 2012, 45, 354–356. [Google Scholar]
  49. Chambers, J.M.; Hastie, T. Statistical Models in S; Wadsworth & Brooks/Cole Advanced Books & Software: Pacific Grove, CA, USA, 1992. [Google Scholar]
  50. Hastie, T.J.; Tibshirani, R.J. Generalized Additive Models; Chapman & Hall: London, UK, 1990. [Google Scholar]
  51. R Core Team. R: A Language and Environment for Statistical Computing; R Foundation for Statistical Computing: Vienna, Austria, 2016. [Google Scholar]
  52. Zacharias, I.; Zamparas, M. Mediterranean temporary ponds. A disappearing ecosystem. Biodivers. Conserv. 2010, 19, 3827–3834. [Google Scholar] [CrossRef]
  53. Sousa, A.; Morales, J.; García-Barrón, L.; García-Murillo, P. Changes in the Erica ciliaris Loefl. ex L. peat bogs of southwestern Europe from the 17th to the 20th centuries AD. Holocene 2013, 23, 255–269. [Google Scholar] [CrossRef]
  54. Sousa, A.; Morales, J.; García-Murillo, P.; García-Barrón, L. Los arroyos costeros del SW de España (Manto Eólico Litoral onubense) como indicadores de cambios en las tendencias de precipitación. In Clima, Ciudad y Ecosistemas. Asociación Española de Climatología; Publicaciones Asociación Española de Climatología; Asociación Española de Climatología: Barcelona, Spain, 2010. [Google Scholar]
  55. Sousa, A.; García-Murillo, P.; Morales, J.; García-Barrón, L. Anthropogenic and natural effects on the coastal lagoons in the southwest of Spain (Doñana National Park). ICES J. Mar. Sci. J. Cons. 2009, 66, 1508–1514. [Google Scholar] [CrossRef]
  56. Kloskowsky, J.; Green, A.; Polak, M.; Bustamante, J.; Krogulec, J. Complementary use of natural and artificial wetlands by waterbirds wintering in Doñana, south-west Spain. Aquat. Conserv. 2009, 19, 815–826. [Google Scholar] [CrossRef]
  57. Sebastián–González, E.; Green, A.J. Reduction of avian diversity in created versus natural and restored wetlands. Ecography 2015. [Google Scholar] [CrossRef]
  58. Europa Press. Autorizado Riego Superficial Para 724 Hectáreas en El Condado (Huelva). Available online: http://www.europapress.es/andalucia/huelva-00354/noticia-gobierno-autoriza-riego-agua-superficial-724-hectareas-condado-huelva-20160315151103.html (accessed on 26 July 2016).
  59. Compromiso cumplido—Detalle de noticia—CHG. Available online: http://www.chguadalquivir.es/-/compromiso-cumpli-1 (accessed on 26 July 2016).
  60. La Vanguardia. Gobierno autoriza riego agua superficial entorno Doñana e inicia cierre pozos. Available online: http://www.lavanguardia.com/local/sevilla/20160315/40446242387/gobierno-autoriza-riego-agua-superficial-entorno-donana-e-inicia-cierre-pozos.html (accessed on 26 July 2016).
  61. Junta de Andalucía. La Junta invierte 16,2 millones en el trasvase que garantizará el abastecimiento de agua al Condado onubense (Consejería de Medio Ambiente y Ordenación del Territorio, Junta de Andalucía). Available online: http://goo.gl/kNAdnr (accessed on 26 July 2016).
Figure 1. (A) Study area. The colored areas represent the Doñana’s aeolian sands and the different colors represent different levels of protection. The maximum extent of temporary waters is overlaid on the aeolian sands. Lines indicate the limits of the “Almonte-Marismas” aquifer, Doñana National Park, Natural Park, and Biosphere Reserve; (B) Location of the study area in southwestern Spain. Numbers locate some areas mentioned in the text: (1) Paraje Natural Lagunas de Palos y Las Madres; (2) Tinto and Odiel marshes; (3) Doñana marshes; (4) Rocina protection area.
Figure 1. (A) Study area. The colored areas represent the Doñana’s aeolian sands and the different colors represent different levels of protection. The maximum extent of temporary waters is overlaid on the aeolian sands. Lines indicate the limits of the “Almonte-Marismas” aquifer, Doñana National Park, Natural Park, and Biosphere Reserve; (B) Location of the study area in southwestern Spain. Numbers locate some areas mentioned in the text: (1) Paraje Natural Lagunas de Palos y Las Madres; (2) Tinto and Odiel marshes; (3) Doñana marshes; (4) Rocina protection area.
Remotesensing 08 00867 g001
Figure 2. Distribution of water abstraction pressures in 2007 using a 5 km buffer on top of the “Almonte-Marismas” aquifer. Irrigated agriculture in yellow, residential areas in black. The other colored areas identify the Doñana’s aeolian sands with each color representing a level of protection. Lines indicate the limits of the “Almonte-Marismas” aquifer, Doñana National Park, Natural Park and Biosphere Reserve.
Figure 2. Distribution of water abstraction pressures in 2007 using a 5 km buffer on top of the “Almonte-Marismas” aquifer. Irrigated agriculture in yellow, residential areas in black. The other colored areas identify the Doñana’s aeolian sands with each color representing a level of protection. Lines indicate the limits of the “Almonte-Marismas” aquifer, Doñana National Park, Natural Park and Biosphere Reserve.
Remotesensing 08 00867 g002
Figure 3. Temporal trends of mean pixel hydroperiod of temporary waters in each protection level and annual precipitation on the “Almonte-Marismas” aquifer. (A) Red circles = National Park, green triangles = Natural Park; (B) blue squares = Biosphere Reserve; (C) gray circles = unprotected; (D) precipitation trend. Only statistically significant trend lines are shown.
Figure 3. Temporal trends of mean pixel hydroperiod of temporary waters in each protection level and annual precipitation on the “Almonte-Marismas” aquifer. (A) Red circles = National Park, green triangles = Natural Park; (B) blue squares = Biosphere Reserve; (C) gray circles = unprotected; (D) precipitation trend. Only statistically significant trend lines are shown.
Remotesensing 08 00867 g003
Figure 4. Temporal trends in the coefficient of variation of the hydroperiod of temporary water pixels in each protection levels. (A) Red circles = National Park; (B) green triangles = Natural Park; (C) blue squares = Biosphere Reserve; (D) gray circles = unprotected. Only significant trend lines are shown.
Figure 4. Temporal trends in the coefficient of variation of the hydroperiod of temporary water pixels in each protection levels. (A) Red circles = National Park; (B) green triangles = Natural Park; (C) blue squares = Biosphere Reserve; (D) gray circles = unprotected. Only significant trend lines are shown.
Remotesensing 08 00867 g004
Figure 5. Relationship between annual accumulated precipitation and mean pixel hydroperiod for each level of protection. (A) Red circle = National Park; (B) green triangles = Natural Park, (C) blue squares = Biosphere Reserve, (D) gray circles = unprotected. Only statistically significant trend lines are shown.
Figure 5. Relationship between annual accumulated precipitation and mean pixel hydroperiod for each level of protection. (A) Red circle = National Park; (B) green triangles = Natural Park, (C) blue squares = Biosphere Reserve, (D) gray circles = unprotected. Only statistically significant trend lines are shown.
Remotesensing 08 00867 g005
Figure 6. Density plot of the mean R-squared value of water bodies. The R-squared comes from an ordinary least squares regression between mean annual accumulated precipitation on the aquifer and mean hydroperiod of the pixels classified as temporary waters. Higher values indicate greater dependence of hydroperiod from annual precipitation. Arcsine transformed mean R-squared values on the X axis (radians). Protection level indicated with colors: high protection = red, medium protection = green, low protection = blue, unprotected = gray. The rug plots below indicate the distribution of water bodies. High protection (red) and unprotected (gray) lines are drawn thicker for clarity.
Figure 6. Density plot of the mean R-squared value of water bodies. The R-squared comes from an ordinary least squares regression between mean annual accumulated precipitation on the aquifer and mean hydroperiod of the pixels classified as temporary waters. Higher values indicate greater dependence of hydroperiod from annual precipitation. Arcsine transformed mean R-squared values on the X axis (radians). Protection level indicated with colors: high protection = red, medium protection = green, low protection = blue, unprotected = gray. The rug plots below indicate the distribution of water bodies. High protection (red) and unprotected (gray) lines are drawn thicker for clarity.
Remotesensing 08 00867 g006
Figure 7. Hydroperiod trend (Kendall’s tau) of water bodies on the Doñana’s aeolian sands. Points represent the centroid of the water body. Blue colors represent increasing hydroperiods and red colors decreasing hydroperiods. The reference image indicates with color the levels of protection.
Figure 7. Hydroperiod trend (Kendall’s tau) of water bodies on the Doñana’s aeolian sands. Points represent the centroid of the water body. Blue colors represent increasing hydroperiods and red colors decreasing hydroperiods. The reference image indicates with color the levels of protection.
Remotesensing 08 00867 g007
Figure 8. Partial effects of predictors in the best gam model for the correlation between annual hydroperiod and precipitation (R-squared) for all water bodies. (A) Irrigation = smoothing spline of log distance to irrigated cultures; (B) residential = smoothing spline of log distance to residential areas, (C) altitude = altitude above sea level; (D) protection = protection level (factor, four levels).
Figure 8. Partial effects of predictors in the best gam model for the correlation between annual hydroperiod and precipitation (R-squared) for all water bodies. (A) Irrigation = smoothing spline of log distance to irrigated cultures; (B) residential = smoothing spline of log distance to residential areas, (C) altitude = altitude above sea level; (D) protection = protection level (factor, four levels).
Remotesensing 08 00867 g008
Figure 9. Partial effect of predictors in the best gam model for annual hydroperiod trend corrected for precipitation (Kendall’s tau) for all water bodies; (A) Irrigation = smoothing spline of log distance to irrigated cultures; (B) altitude = altitude above sea level; (C) protection = protection level (factor, four levels).
Figure 9. Partial effect of predictors in the best gam model for annual hydroperiod trend corrected for precipitation (Kendall’s tau) for all water bodies; (A) Irrigation = smoothing spline of log distance to irrigated cultures; (B) altitude = altitude above sea level; (C) protection = protection level (factor, four levels).
Remotesensing 08 00867 g009
Figure 10. Interaction plot of hydroperiod trend (Kendall’s tau) between protection level (high = National Park, medium = Natural Park, low = Biosphere Reserve) and type of flooding: natural = blue, or artificial = orange. n = 381 photo-interpreted water bodies.
Figure 10. Interaction plot of hydroperiod trend (Kendall’s tau) between protection level (high = National Park, medium = Natural Park, low = Biosphere Reserve) and type of flooding: natural = blue, or artificial = orange. n = 381 photo-interpreted water bodies.
Remotesensing 08 00867 g010
Figure 11. Partial effects of predictors in the best gam model for correlation between annual hydroperiod and precipitation (R-squared) for a random selection of natural water bodies. (A) Irrigation = smoothing spline of log distance to irrigated cultures; (B) residential = smoothing spline of log distance to residential areas; (C) protection = protection level (factor, four levels).
Figure 11. Partial effects of predictors in the best gam model for correlation between annual hydroperiod and precipitation (R-squared) for a random selection of natural water bodies. (A) Irrigation = smoothing spline of log distance to irrigated cultures; (B) residential = smoothing spline of log distance to residential areas; (C) protection = protection level (factor, four levels).
Remotesensing 08 00867 g011
Figure 12. Partial effects in the best gam model for temporal trend in annual hydroperiod corrected for precipitation (Kendall’s tau) for natural water bodies. (A) Residential = smoothing spline of log distance to residential areas; (B) protection = protection level (factor, four levels).
Figure 12. Partial effects in the best gam model for temporal trend in annual hydroperiod corrected for precipitation (Kendall’s tau) for natural water bodies. (A) Residential = smoothing spline of log distance to residential areas; (B) protection = protection level (factor, four levels).
Remotesensing 08 00867 g012
Table 1. Distribution and density of temporary water bodies in the Doñana’s aeolian sands in relation to protection level.
Table 1. Distribution and density of temporary water bodies in the Doñana’s aeolian sands in relation to protection level.
Protection LevelProtection FormTotal Area (km2)No. Temporary Water BodiesMaximum Extent Temporary Waters (km2)Mean Pixel Hydroperiod 1985–2014 (Days)Density (Ponds/km2)% Area Flooded × 102
highNational Park239157513.611.66.595.69
mediumNatural Park25510739.111.44.213.57
lowBiosphere Reserve3166122.715.21.940.84
unprotected 1084184.223.93.873.93
Total 918367829.613.54.013.2
Table 2. Distribution of a random selection of 100 temporary water bodies from each protection level photo-interpreted into different types of water bodies.
Table 2. Distribution of a random selection of 100 temporary water bodies from each protection level photo-interpreted into different types of water bodies.
Protectio LevelNatural PondStreamNo VegetationDisperse VegetationArtificial PondIrrigation PondPine WoodTotal
High244065007100
Medium381254005100
Low12147534100100
Unprotected74224312102100
Total812331215162014400
%20.255.757.7553.75453.5
Table 3. Distribution of a random selection of 100 temporary water bodies from each protection level classified according to type of inundation.
Table 3. Distribution of a random selection of 100 temporary water bodies from each protection level classified according to type of inundation.
Protection LevelArtificialErrorNaturalTotal
high0793100
medium2593100
low21079100
unprotected42751100
Total6519316400
%16.254.7579
Table 4. Linear model for temporal trend in hydroperiod corrected for precipitation (Kendall’s tau) obtained by simplification from a full model with interactions considering the factors: protection level (four levels) and flood type (two levels), distance to irrigated cultures, and distance to residential areas (log transformed). n = 381 water bodies.
Table 4. Linear model for temporal trend in hydroperiod corrected for precipitation (Kendall’s tau) obtained by simplification from a full model with interactions considering the factors: protection level (four levels) and flood type (two levels), distance to irrigated cultures, and distance to residential areas (log transformed). n = 381 water bodies.
EstimateStd. Errort Valuep
Intercept−0.0410.032−1.2380.217
protection(medium)−0.0660.105−0.6320.527
protection(low)0.1110.0452.4590.014
protection(unprotected)0.1360.0255.521<0.0001
flood type(natural)0.0170.0290.5650.573
protection(medium): flood type(natural)0.1090.1051.0320.303
protection(low): flood type(natural)−0.1410.046−3.1080.002
protection(unprotected): flood type(natural)NANANANA
Residual standard error0.141374 df
Adjusted R-squared0.141
F-statistic (6374 df)11.42p < 0.0001

Share and Cite

MDPI and ACS Style

Bustamante, J.; Aragonés, D.; Afán, I. Effect of Protection Level in the Hydroperiod of Water Bodies on Doñana’s Aeolian Sands. Remote Sens. 2016, 8, 867. https://0-doi-org.brum.beds.ac.uk/10.3390/rs8100867

AMA Style

Bustamante J, Aragonés D, Afán I. Effect of Protection Level in the Hydroperiod of Water Bodies on Doñana’s Aeolian Sands. Remote Sensing. 2016; 8(10):867. https://0-doi-org.brum.beds.ac.uk/10.3390/rs8100867

Chicago/Turabian Style

Bustamante, Javier, David Aragonés, and Isabel Afán. 2016. "Effect of Protection Level in the Hydroperiod of Water Bodies on Doñana’s Aeolian Sands" Remote Sensing 8, no. 10: 867. https://0-doi-org.brum.beds.ac.uk/10.3390/rs8100867

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

Article Metrics

Back to TopTop