Next Article in Journal
A GIS-Based Hydrological Modeling Approach for Rapid Urban Flood Hazard Assessment
Next Article in Special Issue
Numerical Representation of Groundwater-Surface Water Exchange and the Effect on Streamflow Contribution Estimates
Previous Article in Journal
Microalgal-Bacterial Granular Sludge Process in Non-Aerated Municipal Wastewater Treatment under Natural Day-Night Conditions: Performance and Microbial Community
Previous Article in Special Issue
Methods in Capturing the Spatiotemporal Dynamics of Flow and Biogeochemical Reactivity in Sandy Beach Aquifers: A Review
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Comparative Analysis of Runoff and Evaporation Assessment Methods to Evaluate Wetland–Groundwater Interaction in Mediterranean Evaporitic-Karst Aquatic Ecosystem

Center of Hydrogeology, Department of Geology, Facutly of Science, University of Málaga, E-29071 Málaga, Spain
*
Author to whom correspondence should be addressed.
Submission received: 14 April 2021 / Revised: 17 May 2021 / Accepted: 18 May 2021 / Published: 25 May 2021

Abstract

:
This work compares the applicability of several free-surface evaporation and runoff equations in simulating water level variations of small Mediterranean wetlands. The Amarga and Jarales wetland are two pilot sites with an evaporite-karst genesis located in southern Spain. The water level was continuously recorded in both wetlands, and exhaustive weather monitoring was performed. The combined datasets have permitted quantification of the surficial elements of their water budget (precipitation, runoff, and evaporation). Several campaigns of groundwater level measurements were also done to characterize the direction of groundwater flows. The morphometrical analysis of the Jarales wetland was accurately performed based on a LiDAR dataset. A total of 225 limnimetric simulations of the Jarales (90) and Amarga (135) wetlands were performed, combining different evaporation and runoff equations. During the study period, the curve number method, coupled with the Penman equation, reached the Jarales wetland’s best calibrations. The Vardavas–Fountoulakis modification of the Penman model fit better with the Amarga wetland record. The obtained results permit specification of the water budget of both wetlands during several years and confirm that the groundwater–surface water relationship affects the wetland hydric dynamic to different degrees. Nonetheless, the limnimetric models were calibrated for a short period, including dry years, making it necessary to extend the control period longer and validate the models under different hydroclimatic conditions. Finally, the differences between wetland functioning are explained in a conceptual hydrological model that can be useful for wetland conservation and management of related aquatic ecosystems. The understanding of the origin and fate of water in wetlands permits assessment of how future scenarios would affect hydric functioning and suggests adequate conservation measurements.

1. Introduction

Wetlands are valuable aquatic ecosystems with a permanent or temporal presence of water at the surface or in the soil [1] that provides a broad list of environmental and human services (e.g., [2,3,4]). Wetlands are frequently related to groundwater, particularly in semi-arid areas, where rainfall presents a seasonal pattern and its annual amount is lower than the evaporation rate [5]. Groundwater–surface water interactions in wetlands are complex, as they are temporarily and spatially variable and sensitive to topographical, geological, and climatic conditions [6,7,8]. During the last decades, some methods have been developed to quantify the groundwater–wetland exchange. However, this remains challenging due to the temporal variation of the water storage and fluxes in wetlands and their sediments’ heterogeneity and low permeability [9,10].
Much of the research indirectly determines the net groundwater–wetland exchange (groundwater input–wetland infiltration) as the wetland water budget residual [11,12,13,14,15,16,17,18,19]. Such an approach is common in areas where the groundwater fluxes are difficult to determine with direct measurements. However, its accuracy dramatically depends on the uncertainties in estimating the rest of the budget components: direct precipitation, overland runoff, evaporation, surface water inflow and outflow, and the changes in water storage [10]. As the estimation of volumes (storage, inputs, and outputs) is dependent on the wetland area, the error can be more extensive in wetlands located in low-gradient basins, where the flooded surface can broadly vary (e.g., [9,20]). The uncertainties of the wetland water budget can be partially reduced using precise digital elevation models based on modern tools for elevation data acquisition, such as LiDAR (Laser Imaging Detection and Ranging), which results in more reliable wetland bathymetry [21].
In semi-arid areas, such as the Mediterranean region, wetlands are particularly vulnerable to climate change due to the expected increase in the evaporation rates and the precipitation pattern changes [22]. Thus, defining the hydrogeological functioning of such aquatic ecosystems is vital for their management and conservation. For that reason, it is necessary to develop models that accurately estimate the weight of each of the components in the budget to compute future hydrological changes in this sort of wetlands. In that sense, several comparisons of estimation methods for free-surface water evaporation are found in the literature [23,24,25,26], but equations for overland runoff are more rarely discussed. However, this is a crucial aspect to consider in small wetlands in semi-arid regions [10,27].
The Amarga and Jarales are two wetlands located in the province of Cordoba (southern Spain, Figure 1). The regional government protects them as Natural Reserves, they are part of the Ramsar List of Wetlands, and they are a Special Protection Area (SPA) under the European Union Directive on the Conservation of Wild Birds. The genesis of both wetlands is related to dissolution/karstification processes affecting the underlying evaporite rocks (gypsum, anhydrite, and halite) embedded in a Triassic clayey matrix that constitutes the best part of the so-called Chaotic Subbetic Complexes (CSC) [28]. Despite having a similar origin, both wetlands have different hydrological functioning. Previous research carried out in the area [12,29,30,31,32] described their hydrological functioning and quantified their water budget components. Most of the research reached similar conclusions regarding the groundwater–wetland interaction, even though they used diverse methods for evaporation and runoff estimation. They performed limnimetric simulations based on utterly different soil parameters, partly because neither the observation period nor the frequency of control/modeling was the same.
This study aimed to reach more profound insights into the hydrologic functioning of evaporite karst wetlands in semi-arid regions, emphasizing the surface water–groundwater relationship. The applicability of several methods used in estimating the water budget component is also discussed, particularly those for runoff and evaporation computation. The final goal was to create a hydrological conceptual model that helps to understand groundwater’s contribution to the wetlands and its temporal variations as a tool for environmental management of the dependent aquatic ecosystems.

2. Geological and Hydrological Settings

The Amarga and Jarales wetlands are located on a plateau of 80 km2 between the rivers Genil and Anzur (Figure 1a). The outcropping rocks mainly belong to the CSC, formed by a clayey-evaporitic megabreccia in which olistoliths (blocks) of diverse nature and size are embedded, including limestones, dolostones, sandstones, and massive gypsum enclaves, among others [28]. Discordantly over the CSC rocks, post-orogenic deposits are constituted by sandstones, conglomerates, marls, calcarenites, and unconsolidated sediments from Miocene to Quaternary ages (Figure 1b,c).
Figure 1. (a) Location and general settings of the study area, including groundwater table contours. Geology of the wetland watersheds, after Martín-Serrano [33]: (b) Amarga and (c) Jarales.
Figure 1. (a) Location and general settings of the study area, including groundwater table contours. Geology of the wetland watersheds, after Martín-Serrano [33]: (b) Amarga and (c) Jarales.
Water 13 01482 g001
The landscape is characterized by gentle slopes and plenty of endorheic areas, some of which are occupied by wetlands (Figure 1a) with different hydroperiods. Despite the abundance of low permeability materials in the CSC, the presence of small outflows, wells, and boreholes evidence groundwater flows in the media [30,34]. Moreover, the plateau is drained through three brine springs emplaced by the Genil and Anzur Rivers (Figure 1a), between 247 and 300 m above sea level (asl) [35]. Dissolution and karstification processes promoted the subsidence and collapse of the terrain, giving rise to closed-depression landforms. Furthermore, based on the alignment of fractures and slopes, previous research suggests that tectonics could also play an essential role in the development of karst depression [29,36].
The Jarales wetland is located in the central part of the area (Figure 1c) at 406.6 m asl and has a flat-bottom basin. Its average flooding surface is 5 ha, while the maximum water level recorded is about 3 m. Its flooding pattern is irregular because it commonly gets dry during the summer of average years (from the standpoint of rain amount) but not when the hydrological year is considered wet. On the contrary, no flooding is observed during dry years, not even between autumn and springtime, when most of the annual rainfall occurs. The Amarga wetland is located 3 km NW of the Jarales wetland (Figure 1a), at an altitude of 366 m asl. It has a deep basin and a permanent hydroperiod. Its water stage largely varies from 2.5 m to nearly 6 m high, and its average flooding area extends 4 ha.
The watershed of the Jarales wetland extends 90.4 ha over post-orogenic deposits (Miocene-Quaternary) outcrops (mainly marls, sandstones, and conglomerates) but also on clayey-evaporitic rocks of the CSC (Figure 1b,c). The same materials are on the surface of the Amarga wetland watershed (256.9 ha), although it presents a higher proportion of clays of the CSC, as well as dolostone olistoliths. The morphological characteristics of both wetlands’ watersheds (Table 1) reveal similar shape characteristics, although Amarga‘s is relatively less elongated and has a deeper bottom. The protected area around the wetlands is defined by a 500 m buffer around the maximum flooding surface’s perimeter, although only the first 50 m have the highest level of protection given by the Natural Reserve (Figure 1).
The climate in the area is Mediterranean, with arid and warm summers. The historical mean values of annual precipitation (1963/64–2016/17) and air temperature (1980/81–2016/17) are 508 mm and 17.2 °C, respectively [35]. The monitoring period used in the present study comprises the hydrological years 2014/15, 2015/16, and 2016/17, and it is considered as dry since the annual rainfall registered in a weather station located in the central part of the area (Figure 1a) was 324 mm, 403 mm, and 366 mm, in each case. This dry period took place after some wet years: 2009/10 (913 mm), 2010/11 (675 mm), and 2012/13 (821 mm), according to the historical rainfall distribution pattern (Figure 2; [35]).

3. Methods

Four field campaigns of water table measurements in 39 wells and boreholes were carried out under different hydroclimatic situations (May and September 2015, January and March 2016). Those measurements permitted the drawing of schematic water-table contour maps that characterize the direction of groundwater flows in the study area.
Two automatic water-depth data loggers were installed next to respective staff gauges in the Amarga and Jarales wetlands to obtain daily records of water surface fluctuations. The register extended from October 2014 to September 2017 for the Amarga wetland, while it lasted from March 2014 to July 2015 in the Jarales wetland, when it got dry. Additionally, single staff lectures were done in the Jarales wetland until September 2017, when occasional flooding occurred.
The changes in the volume of water stored in a wetland (V) can be approximated as a function of the variations in the wetland stage (h) and the flooding area (A) [10]:
Δ V = ˜ Δ h A + Δ A 2
Thus, the change in the wetland storage between two heights, h1 and h2, can be calculated given the size of their corresponding flooding areas (A1 and A2, respectively):
V h 1 h 2 = ˜ A 1 A 2 2 · h 1 h 2
The flooding surface related to different wetland stages was calculated, considering height intervals of 0.5 m. To that end, a Digital Terrain Model (DTM) of the area, with a resolution of 0.5 m × 0.5 m, was created, based on LiDAR data [40] obtained from the National Geographical Institute of Spain (IGN). This dataset has a mean square error in vertical measures below 10 cm and a nominal data density of 0.5 points/m2. This density could be low for the DTM resolution; however, minor deviations are observed when interpolating the DTM, and therefore this issue has low implications on the volume calculation. The point cloud cover was generated from a photogrammetric flight performed in summer 2015 when the Jarales wetland was dry, allowing for the whole wetland basin’s morphological characterization. In Figure 3a,b, the considered h values are shown versus their corresponding stored volumes (V) and flooding areas (A). The equations defining those relationships (the filling curve and the hypsometric curve, respectively) were obtained using the software CurveExpert. During the recording period, the Amarga wetland’s stage descended below the water table’s height registered by the photogrammetric flight. Consequently, LiDAR data do not provide helpful elevation information for assessing the Amarga wetland storage changes during the control period. Nevertheless, García-Ferrer et al. [41] published bathymetric data of the Amarga wetland used by Aljibe Consultores [31] to obtain the filling curve and the hypsometric curve defining the wetland stage relationship with the storage and flooding surface (Figure 3c,d). Additionally, the delimitation of the watershed of both wetlands was created based on the DTM obtained from the LiDAR dataset. To that end, flow direction and flow accumulation covers were calculated, and all of the micro-watersheds flowing toward the wetlands were merged later.
As the drainage network in the study area is not well defined and there are no permanent surface water inflows to any of the wetlands, the only two surficial water inputs considered were direct precipitation and diffuse overland runoff. Direct rainfall inputs were obtained by multiplying the wetland surface at each time and the precipitation was measured at a weather station located in the central part of the plateau (Figure 1). Surficial runoff was estimated using simple rainfall-runoff models that are commonly used in wetland water budget studies: the curve number (CN) method [16,42,43] and the soil water balance (SWB) [12,22,32,44,45].
Figure 3. Wetland stage (h) versus stored volume (V) and flooding area (A) of the Jarales (a,b) and Amarga wetlands (c,d). The equations of the Amarga wetland were obtained from Aljibe Consultores [31].
Figure 3. Wetland stage (h) versus stored volume (V) and flooding area (A) of the Jarales (a,b) and Amarga wetlands (c,d). The equations of the Amarga wetland were obtained from Aljibe Consultores [31].
Water 13 01482 g003
The CN method was developed by the US Soil Conservation Service [46]. It considers that runoff starts when a rainfall event exceeds a particular threshold value called initial abstraction (Ia), independent of whether the soil is saturated with water. Then, runoff (Q) depends on precipitation (P) and Ia:
Q = P I a 2 P 4 I a     for     P > I a
Q = 0       for     P I a
The Ia value depends on land use, slope percentage, and soil texture. In general terms, a low value of Ia (nearly 0 mm) corresponds to low permeability soils (clayey), scarce vegetation, and steep slopes. On the contrary, Ia has high values (>100 mm) if the conditions are less propitious for runoff generation: permeable soils (sandy), dense vegetation, and smooth relief. The Ia values were selected according to the criteria of the Spanish Ministry of Public Works [47], and they are relatively low (10, 15, and 19 mm), since the wetland watersheds have gentle slopes and low permeability soils, occupied mainly by sparse natural vegetation or olive orchards. The agricultural practices related to olive growing include plowing and eliminating other arboreal or herbaceous plants.
The SWB [48] assumes that the soil has a specific water storage capacity in the root zone, called “available water-holding capacity” (AWC). The water held within the soil initially varies depending on the input due to precipitation and the output provoked by the actual evapotranspiration (ETA), which withdraws water from the soil until running out or until its maximum potentiality, called potential evapotranspiration (ETP). When the volume of water stored in the soil and the precipitation exceeds ETP, there is an excess of water that produces either surficial runoff, infiltration, or both. Such excess is known as effective precipitation (Pe). To compare different methods, the values of daily ETP were estimated by applying different equations: Thornthwaite [49], Hargreaves [50], and Blaney–Criddle [51]. The selected AWC values were 50, 60, 75, 100, and 125 mm. As the permeability of the bedrock in the study area is generally low, the runoff coefficients (RC) applied to Pe were high: 50%, 75%, and 100%.
The computation of free-surface water evaporation was made using the climatic data recorded in the weather station placed in the center of the study area (Figure 1a) by applying the Penman equation [52]. This evaporation model is one of the best for estimating evaporation from shallow wetlands (<2 m of water column) according to a meta-study by McMahon et al. [26], in which many previous studies were reviewed ([23,25,53,54,55,56], among others). However, the Penman equation does not work so well when applied to water bodies more than 2 m depth, as it does not take into account the heat transfer toward the water mass, which buffers the evaporation effect [57,58,59,60,61]. Vardavas and Fountoulakis [62] proposed an easy-to-apply modified version of the Penman model that considers the net energy exchange between the atmosphere and the water body, based on the monthly temperature variations measured in the wetland water and the monthly mean wetland stage. Since the water column in the Amarga wetland is several meters high, the Vardavas and Fountoulakis modification was also considered in the present work. Additionally, the Penman–Monteith equation [63] was included in this study to compare the results obtained from different methods.
The simulation of the daily changes on wetland storage during the observation period was made using their respective hypsometric and filling curves (Figure 3), rainfall recorded at the weather station (Figure 1a), and overland runoff and evaporation obtained from the different methods mentioned. The underground components of the wetland water budgets were inferred from the differences between the recorded and the computed limnimetric evolutions.

4. Results

During the four groundwater-level measurement campaigns, minor variations of the groundwater table were observed, with an average drop of 0.6 m from May 2015 to March 2016. In most wells, a general decreasing trend was observed (Figure 4), although in some cases, a groundwater level rise was registered between September 2015 and March 2016. That was mainly seen in wells near the Jarales wetland, particularly those emplaced at higher altitudes than the wetland itself (Figure 4).
The water table contour outlined in January 2016 is presented as an example in Figure 1a since the results obtained from the rest of the campaigns and the data retrieved from previous research [30] are relatively similar. According to the outline, the water table’s morphology is somehow adapted to the land surface, typical of low permeability media. Thus, a potentiometric dome exists in the plateau’s central area, creating a hydrogeological divide from which groundwater flows in a radial and centrifugal way toward the edges of the system.
The fit of each limnimetric simulation to the wetland stage record series was evaluated using the root-mean-square error (RMSE) and the correlation coefficient (R2), whose values are shown in Figure 5. The RMSE measures the magnitude of the predictive error and evaluates the fit accuracy based on the differences between the computed values and the observed values. The R2 tests the relationship between the two data series, indicating the covariance of both variables independent of the scale.
Regarding the runoff models, the worst results were obtained using the Thornthwaite equation into the SWB, which provides high RMSE values and low R2 (Figure 5). Hargreaves and Blaney–Criddle methods give more relevant results, although they are slightly better in the second case. The simulations derived from the CN equation are generally well correlated with the limnimetric record, and their RMSE values are low (Figure 5). As for the evaporation models, the Penman equation generates better results in the Jarales wetland than those obtained from the Penman–Monteith method, while in the Amarga wetland, the best fittings were obtained with the Vardavas–Fountoulakis correction to the Penman model.
Figure 5. The root-mean-square error (RMSE) and the correlation coefficient (R2) of each limnimetric simulation computed in the Jarales and Amarga wetlands.
Figure 5. The root-mean-square error (RMSE) and the correlation coefficient (R2) of each limnimetric simulation computed in the Jarales and Amarga wetlands.
Water 13 01482 g005
The four best limnimetric simulations in the Jarales wetland for the control period, according to the results of Figure 5, are represented in Figure 6. For the limnimetric simulation of the Jarales wetland, the CN runoff with an Ia of 19 mm was included twice: once coupled with the Penman evaporation model (JW1) and the other using the Penman–Monteith equation (JW2). The other two computed series considered the Penman model with the runoff estimated by mean of the SWB, using the Blaney–Criddle ETP: one with an AWC of 100 mm (JW3) and the other with an AWC of 75 mm and a runoff coefficient of 50% (JW4). Although these simulations were created for the whole control period, the RMSE and the R2 values of those simulations were calculated only from April 2014 to July 2015, as no continuous record of the wetland stage is available after. None of the represented series predict flooding during the hydrological year 2016/17, agreeing with the observed data.
The JW3 simulation has the lowest RMSE value (0.19 m) and one of the best R2 (0.957). However, in that simulation, the water stored within the soil never exceeds the AWC (nor in the rest of the simulations with identical RMSE and R2 values in Figure 5). Consequently, no runoff generates during the heaviest rain events, contrary to in situ observations. If the field capacity is lowered to 75 mm and a RC of 50% is applied (JW4), the model simulates a certain runoff amount, but it does not replicate the flooding starting in November 2015 (Figure 3). The two simulations that include the CN runoff (JW1 and JW2) better replicate the water level rises recorded, particularly in November 2015 (Figure 6). The RMSE and R2 values of JW1 are slightly better than those of JW2, partly because the last model does not predict the wetland’s desiccation toward the end of the hydrological year 2014/15.
Figure 7 shows the four selected simulations of the evolution of the Amarga wetland stage. Three of them were obtained applying the CN method with an Ia of 15 mm and the evaporation according to Penman (AW1), Penman–Monteith (AW2), and Vardavas–Fountoulakis (AW3). In this case, the best fitting corresponds to an Ia lower than the one in the Jarales basin (19 mm). That is in agreement with the lithological characteristics of the catchment area (Figure 1b,c), which is less permeable in the Amarga wetland case due to a more significant outcropping proportion of CSC rocks that confer the soils a clayey texture. The other one (AW4) was computed based on the SWB, with the ETP after Blaney–Criddle (AWC: 50 mm, RC: 75%) and the evaporation of Vardavas and Fountoulakis. AW4 is one of the simulations with the lowest error (Figure 5), as it only has slight differences with the actual values for most of the recording period, except for part of 2016 (Figure 7). However, the computed water table rises do not resemble the observed ones, and some of the actual ascents were not replicated. The latter occurs because the runoff generating such water level changes is caused by rainfall occurring when the soil was not saturated. Contrary to the SWB, the CN method (used in AW1, AW2, and AW3) predicts the runoff generation in those circumstances, as it does not depend on the previous state of the AWC but on the total amount of rain.
Based on the previous analysis, JW1 and AW3 simulations were closer to the actual records. Thus, each wetland’s water budget (Table 2 and Table 3) was calculated on a monthly and annual scale, using direct rainfall precipitation inputs and the evaporation and runoff simulation used in JW1 and AW3. The difference between the actual and simulated water storage volume was attributed to the net groundwater–wetland exchange (ΔG), meaning the difference between groundwater inputs and outputs due to infiltration in the wetland bed. Thus, positive values of ΔG stand for a net underground input. During 2014/15, evaporation was the most significant component in the Jarales wetland budget (34.9 dam3; Table 2). Direct precipitation inputs were 13.5 dam3, while the runoff contribution was 1 order of magnitude lower (1.6 dam3), and it only provided inflows in October and November 2014. The monthly values of ΔG vary between −4.9 dam3 (October 2014) and 2.7 dam3 (November 2014). The annual water budget of the Jarales wetland in 2014/15 indicates that there would be a net groundwater input of 0.6 dam3.
In the Amarga wetland, evaporation accounts for 74.2 dam3 in 2015/16 and 66.7 dam3 in 2016/17 (Table 3), while precipitation inputs are 20.3 dam3 and 15.0 dam3, respectively. Runoff contributions vary notably, from 13.1 dam3 in 2015/16 to 2.3 dam3 in 2016/17. This last value must be an underestimate, as the simulated level rise in November 2016 is much lower than the actual one (Figure 7). That underestimation could relate to the fact that the intensity of rainfall was lower in this period than in the previous years; thus, the CN method did not predict much runoff. The total simulated variation of stored water (ΔV) between 2015/16 and 2016/17 is −90.2 dam3 (−40.8 and −49.4 dam3, respectively), which would reduce the stored volume at the end of the simulation period (September 2017) to 37.7 dam3 (Table 2), in contrast to the actual storage (91.7 dam3).

5. Discussion

5.1. Choosing the Runoff and Evaporation Equations for the Limnimetric Simulation

The correlation degree observed between the actual limnimetric record and the simulations made with the CN method is generally good, and the magnitude of the errors is low (Figure 5). That could be due to the low permeability of most of the rocks outcropping in the wetland watersheds. Thus, water does not infiltrate rapidly into the soil, resulting in runoff generation during heavy rain events even if it is not completely saturated. Such a process fits with the theoretical basis of the CN. The dry hydroclimatic conditions registered during the study period would have favored that the water held in the soil rarely exceeded the AWC so that the SWB was not able to predict the runoff generation in certain moments. The best simulations obtained applying the SWB considered low values of AWC (50–75 mm, Figure 5) that were coherent not with the characteristics of the outcropping rocks but with the limited thickness of the soil and its limited radical depth. Both factors combined with the semi-arid regime of the study area, and the climatic conditions (dry years) would have enhanced the weight of the ETP, which impedes the SWB to be suitable to the studied reality. Based on all that, the CN method is the most adequate for estimating runoff in the explained context.
Regarding the estimation of free-water surface evaporation, the simulations integrating the Penman model fit slightly better in the Jarales wetland (Figure 5). Additionally, the JW1 simulation predicts the wetland’s desiccation in the summer of 2015 using this equation, while JW2 does not replicate this issue using the Penman–Monteith evaporation (Figure 6). The good applicability of Penman equation is in line with other research comparing methods for the computation of free-water surface evaporation [24,26]. In shallow lakes, like the Jarales wetland, the Penman equation can be used to estimate lake evaporation, as the advected energy and changes in seasonal stored energy can be ignored. However, the Penman–Monteith method was designed to calculate reference-crop evaporation, and it includes simplifications of the atmospheric pressure and the latent heat of vaporization [63] that affect open-water evaporation.
On the contrary, the RMSE values derived from the simulations using Penman–Monteith in the Amarga wetland are generally lower than those obtained with the Penman model, although the best adjustment was reached using Penman plus Vardavas and Fountoulakis correction (Figure 5). The Penman equation’s assumptions are not entirely applicable to deep water bodies [26,60], such as the Amarga wetland. On one hand, the runoff inputs to the Amarga wetland can create rises of tens of centimeters in the water column. Such inflows may have a different temperature than the average water temperature in the water body, so advected energy cannot be neglected. On the other hand, the water body can store and transport heat, so heating and cooling must also be considered as energy fluxes [64]. As the depth of the mixing varies in time and space, it is hard to estimate changes at a short time step, although a number of methods offer satisfactory approaches [57,59,61]. The Penman modification by Vardavas and Fountoulakis is one of them [62] and is based on monthly surface water temperatures, which are easy to acquire. The method reached satisfactory results in Australia, including case studies with Mediterranean and semi-arid climates. In the Amarga wetland, with a similar climatic setting, the limnimetric simulation is also improved when applying the modification. Thus, it is highly recommended to check its applicability in other wetlands in similar contexts.
Previous work carried out in the study area [12,29,32] performed limnimetric simulations and wetland water budgets using the SWB. They reached satisfactory results using AWC values between 103 and 260 mm, which are high considering the limited soil thickness. Said works considered average monthly data and/or wet years from the standpoint of precipitation amount. The use of those values for the period 2014/15–2016/17, integrated with dry and medium years, would not have resulted in daily-stepped runoff. In wet years, the rain widely exceeds the AWC, even if it has relatively high values. Additionally, the monthly-stepped simulations use accumulated values of precipitation and ETP, so they are less precise in predicting the time distribution of the soil saturation status. Therefore, the hydroclimatic conditions and the control and simulation steps are essential for selecting both the equations and the fitting parameters. Thus, generating long daily data series of wetland stage variation would be advisable for obtaining a model able to predict the limnimetric evolution of the wetlands under any climatic conditions.

5.2. Integrating the Results into a Hydrogeological Conceptual Model

During most of the observation period, the Jarales wetland’s hydrological dynamic was mainly influenced by the surficial components of its water budget, particularly evaporation and direct precipitation to a lesser extent (Table 2). The groundwater–wetland interaction only represents a small volume of the whole budget but explains the difference between the actual limnimetric record and the computed series (Figure 6). These results agree with the functioning proposals based on average data generated by previous research [29,30]. According to the monthly water budget of the Jarales wetland during the first six months of control (May to October 2014), the loss of water storage is greater than the evaporation outputs (Figure 6), which suggests that there is net wetland water infiltration (ΔG < 0, Table 2). On the contrary, throughout most of the year 2014/15, net groundwater inputs were estimated. This input would explain the nearly steady state of the wetland water level from December 2014 to May 2015, which does not fit the downward trends simulated, particularly since March (Figure 6).
The Jarales wetland is placed in an intermediate position with respect to the water table contour outline (Figure 1a), between the potentiometric dome in the central part of the study site and the discharge areas situated at lower altitudes. Thus, Jarales would be a flow-through wetland concerning groundwater [5,7,65], meaning it would receive groundwater inputs in some locations (SW), while infiltration losses would occur in others (NE). The water-table changes nearby the wetland can modify the hydraulic gradient around it, making it possible to explain variations in the net groundwater–wetland exchange direction [6,66]. The high water level in the Jarales wetland at the beginning of the control period (Figure 6) resulted from essential inputs during the previous year. This would have caused a hydraulic gradient rise, favoring water infiltration toward areas with a lower table. That explains the negative values of ΔG from May to October 2014 and the difference between the actual and the simulated limnimetric evolutions. In 2014/15, the evaporation rate extensively exceeded precipitation and runoff inputs, causing a pronounced fall of the wetland stage. The low permeability of the terrain would have smoothed the water table’s descent, and therefore a higher gradient would have been established between recharge areas (potentiometric dome, Figure 1) and the wetland. The result is a positive volume of net groundwater–wetland exchange (groundwater input) from November 2014 to May 2015 (Table 2). Eventually, the water table would drop downflow the wetland during summer 2015 (Figure 4), favoring water infiltration and leading to wetland desiccation. The subsequent water table rise upflow the wetland in the first half of the hydrological year 2015/16 may have helped the intermittent swamp of the area, although the general descending trend observed in the downflow wells would have favored infiltration and impeded long-lasting flooding. In any case, this whole interpretation is made based on a short control period, and it will be necessary to confirm it with a more extended observation period that includes both dry and wet years and more detailed groundwater table monitoring.
The permanent hydroperiod in the Amarga wetland is only possible if a considerable net groundwater input exists, as deduced from its ΔG values (Table 3). The estimated volume of that contribution ranges between 20.9 dam3/year (2015/16) and 33.0 dam3/year (2016/17), in agreement with previous research. Aljibe Consultores [31] calculated 28.94 dam3 of average groundwater contribution, and Moya [67] estimated 33.6 dam3 for the hydrological year 1983/84 (average from the rainfall amount standpoint). Nonetheless, an overestimation of the groundwater input in 2016/17 is possible due to the undervaluation of runoff mentioned before. The Amarga wetland’s location toward the potentiometric dome border, together with its deep basin, favors the water table to be always above the wetland bed, as observed in the nearby wells (Figure 4). As a result, groundwater input is permanent and exceeds the water infiltration outputs. Though the gross groundwater inflow to the Amarga wetland was not computed, it could be higher than the net value, as the presence of discharge pointing toward lower altitudes would enable water infiltration downflow. Additionally, locals declare to have observed water infiltration through karst swallow holes formed in the wetland side, and the hydraulic connection between nearby endorheic areas and one of the springs was proven by dye tracer tests [36].
Thus, the differences in the Jarales and Amarga wetlands’ hydric functioning are due to hydrological aspects related to their watershed characteristics and hydrogeological ones linked to wetland–groundwater interactions. A hydrogeological conceptual model (Figure 8) can explain the wetlands’ functioning based on a previous model proposed by Andreo et al. [65] and proved after by Gil-Márquez et al. [68], which describes the general hydrogeological functioning of the CSC. The higher emplacement of the Jarales wetland and its flatter bed impede the permanent wetland–groundwater interaction. When the water table reaches the wetland bed, either net groundwater input or infiltration outputs may occur, depending on the wetland’s water level. While groundwater inflows would be related to local flows, the water infiltration occurring in the wetland bed would be incorporated into the regional groundwater flow net. On the other hand, the Amarga wetland position toward the potentiometric dome border would permit permanent groundwater inputs from nearby areas but would also recharge the region at the highest part of the plateau. The infiltration that may occur in the Amarga wetland would incorporate into the regional flow system that would eventually be drained through the brine springs existing at the plateau’s border [65,68].
Previous research has collected hydrobiogeochemical evidence that agrees with the proposed conceptual model (Figure 8). The Cl concentration and the enrichment in δ18O of both wetlands’ waters are coherent with the hydrochemical and isotopic signature of groundwater from nearby wells after undergoing evaporation [35]. Thus, wetlands are likely receiving groundwater inflows. On the other hand, the analysis and joint interpretation of 3H, 3He, 4He, and 14C in groundwater samples indicate that several flow components with different ages converge at the brine springs located in the system’s border [68]. Finally, the determination of δ15N and δ18O of dissolved NO3− and δ13CDIC in groundwater samples revealed that nitrate removal processes occur at regional scale, linked to organic carbon oxidation taking place in the wetlands [69]. Therefore, water that is infiltrated in the wetlands is incorporated into flowpaths of different lengths until reaching the discharge areas. Not only does the result derived from the present work fit with the above statement, but it also permits an explanation of the hydric variability of the wetlands.
Local diversity of wetland hydroperiods is common in many regions worldwide [18,22,70,71,72,73,74,75]. Differences in surficial inflows and climatic factors (evaporation and precipitation) often explain the variability of flooding periods and lengths [22,76,77,78,79]. However, groundwater exchange is commonly crucial for explaining the different hydroperiods. Several factors may condition the wetland–groundwater connection, including the basin geometry [75], the lithology [80], the wetland position with respect to the groundwater table [5,18,74], and its location within the regional groundwater flow networks [73,81]. In this last sense, Smith and Townley [82] simulate 39 flow regime geometries and found that flow-through wetlands are dominant within regional flow systems, particularly when shallow lakes intersect high-gradient regional flow systems, typical of low-medium permeability media. Neff et al. [83] proved that the presence of depressional wetlands within a broader landscape modifies the regional groundwater flow regime, causing the appearance of flow-through or discharge wetlands far from the regional discharge area. The presence of numerous depressional wetlands within the CSC and the relatively high groundwater gradient can favor the wetland–groundwater interaction, as is suggested by this research.
Previous research proved the wetland–groundwater interaction in the Amarga wetland [12,31], but the present work also reveals the groundwater connection with the Jarales wetland. In both cases, the protected areas are buffers that do not fit with the actual watersheds (Figure 1), so the management measures are not able to act on processes that can affect the wetland through surficial runoff, particularly the carryover of agricultural pollutants and sediments that can potentially lead to the wetlands’ clogging. Moreover, as both wetlands receive regional groundwater inflows in different proportions, it is crucial to advance on the definition of the wetland catchment area, as human activities can affect their hydrological functioning beyond their surficial watershed. Future management, conservation, and restoration policies must take that into account.

6. Conclusions

This work analyzed the relation between groundwater and two Mediterranean wetlands with evaporite karst origin through the interpretation of water table contour outline and the simulation of wetland level fluctuations. The latter has allowed quantification of the wetland water budget components and estimation of the existence and weight of the wetland–groundwater interaction. A hydrological functioning model was proposed in which the relation of the wetlands (mainly Jarales) with the regional groundwater flow network much condition its hydric behaviors. Thus, it is necessary to control the hydrological factors affecting its water budget, including the groundwater surface, so that the hydrogeological context can be fully understood.
The use of different methods for the estimation of runoff and open-water evaporation permitted a contrast of their applicability to the simulation of water level variations of wetland in a semi-arid context. Results evidence that the Penman model is better for evaporation computation, although the correction proposed by Vardavas and Fountoulakis improves it in deep wetlands. The curve number method provided better results than the soil water balance method in dry years, in which the water stored in the soil rarely exceeds the field capacity, but runoff generates nonetheless when rainfall overcomes the soil’s capacity (predominately clayey) to infiltrate water.
The difference between the calibration parameters included in the simulation presented in this work and previous studies evidence the need for performing a continuous limnimetric record in wetlands for long and uninterrupted periods. The obtained series can be the base for developing and calibrating hydric functioning models that apply to any hydroclimatic conditions and that help reach a more accurate knowledge of the wetland–groundwater relationship. That kind of model will be essential for the management, protection, and restoration of wetlands, particularly in areas that are particularly sensitive to climate change, such as the Mediterranean area. Only through an understanding of the origin and fate of water in wetlands is it possible to assess how future scenarios would affect hydric functioning and to propose adequate conservation measurements.

Author Contributions

Conceptualization, J.M.G.-M. and B.A.; methodology, J.M.G.-M. and M.M.; formal analysis, J.M.G.-M. and B.A.; investigation, J.M.G.-M. and M.M.; resources, B.A. and M.M.; data curation, J.M.G.-M.; writing—original draft preparation, J.M.G.-M.; writing—review and editing, J.M.G.-M., M.M., and B.A.; visualization, J.M.G.-M. and M.M.; supervision, B.A.; project administration, B.A.; funding acquisition, B.A. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by JUNTA DE ANDALUCÍA, project numbers P10-RNM-6895R and P11-RNM-8087, and Universidad de Málaga covered the publication fees.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available on request from the corresponding author. The data are not publicly available due to ongoing research using it.

Acknowledgments

This work is a contribution of the Research Group RNM-308 of Junta de Andalucía. The authors want to thank the collaboration provided by the directors and staff of the Natural Reserves of the Jarales and Amarga wetlands during the development of this work.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Mitsch, W.J.; Gosselink, J.G. Wetlands, 5th ed.; Wiley & Sons: Hoboken, NJ, USA, 2015; ISBN 9788578110796. [Google Scholar]
  2. National Research Council. Wetlands: Characteristics and Boundaries; National Academy Press: Washington, DC, USA, 1995; ISBN 0309051347.
  3. Millenium Ecosystem Assessment. Ecosystems and Humanwell-Being: Wetlands and Water Syntesys; World Resources Institute: Washington, DC, USA, 2005. [Google Scholar]
  4. LePage, B.A. Wetlands; LePage, B.A., Ed.; Springer: Dordrecht, The Netherlands, 2011; ISBN 9789400705500. [Google Scholar]
  5. Jolly, I.D.; McEwan, K.L.; Holland, K.L. A review of groundwater–surface water interactions in arid/semi-arid wetlands and the consequences of salinity for wetland ecology. Ecohydrology 2008, 1, 43–52. [Google Scholar] [CrossRef]
  6. Rosenberry, D.O.; Winter, T.C. Dynamics of water-table fluctuations in an upland between two prairie-pothole wetlands in North Dakota. J. Hydrol. 1997, 191, 266–289. [Google Scholar] [CrossRef]
  7. Winter, T.C.; Rosenberry, D.O.; Labaugh, J.W. Where Does the Ground Water in Small Watersheds Come From? GroundWater 2003, 41, 989–1000. [Google Scholar] [CrossRef]
  8. House, A.R.; Thompson, J.R.; Acreman, M.C. Projecting impacts of climate change on hydrological conditions and biotic responses in a chalk valley riparian wetland. J. Hydrol. 2016, 534, 178–192. [Google Scholar] [CrossRef] [Green Version]
  9. Mitchell, J.; Jawitz, J.W. Wetland Water Budgets. In Methods in Biogeochemistry of Wetlands; DeLaune, R.D., Reddy, K.R., Richardson, C.J., Megonigal, J.P., Eds.; SSSA: Madison, WI, USA, 2013; pp. 919–936. [Google Scholar]
  10. Rosenberry, D.O.; Hayashi, M. Assesing and measuring wetland hydrology. In Wetland Techniques. Vol. 1: Foundations; Anderson, J.T., Davis, C.A., Eds.; Springer Sciences Business Media: Dordrecht, The Netherlands, 2013; pp. 87–225. [Google Scholar]
  11. Walton, R.; Chapman, R.S.; Davis, J.E. Development and application of the wetlands dynamic water budget model. Wetlands 1996, 16, 347–357. [Google Scholar] [CrossRef]
  12. Rodríguez-Rodríguez, M.; Fernández, A.; Moral, F. Hydrological regime and modeling of three ponds of the Mediterranean area (south of Córdoba, Spain). Hydrobiologia 2016, 782, 155–168. [Google Scholar] [CrossRef]
  13. Rimmer, A.; Gal, G. Estimating the saline springs component in the solute and water balance of Lake Kinneret, Israel. J. Hydrol. 2003, 284, 228–243. [Google Scholar] [CrossRef]
  14. O′Driscoll, M.A.; Parizek, R.R. The hydrologic catchment area of a chain of karst wetlands in central Pennsylavania, USA. Wetlands 2003, 23, 171–179. [Google Scholar] [CrossRef]
  15. Min, J.H.; Perkins, D.B.; Jawitz, J.W. Wetland-groundwater interactions in subtropical depressional wetlands. Wetlands 2010, 30, 997–1006. [Google Scholar] [CrossRef]
  16. McKillop, R.; Kouwen, N.; Soulis, E.D. Modeling the rainfall-runoff response of a headwater wetland. Water Resour. Res. 1999, 35, 1165–1177. [Google Scholar] [CrossRef]
  17. Hood, J.L.; Roy, J.W.; Hayashi, M. Importance of groundwater in the water balance of an alpine headwater lake. Geophys. Res. Lett. 2006, 33, 1–5. [Google Scholar] [CrossRef] [Green Version]
  18. Castañeda, C.; García-Vera, M.Á. Water balance in the playa-lakes of an arid environment, Monegros, NE Spain. Hydrogeol. J. 2008, 16, 87–102. [Google Scholar] [CrossRef] [Green Version]
  19. Ludwig, A.L.; Hession, W.C. Groundwater influence on water budget of a small constructed floodplain wetland in the Ridge and Valley of Virginia, USA. J. Hydrol. Reg. Stud. 2015, 4, 699–712. [Google Scholar] [CrossRef] [Green Version]
  20. Hayashi, M.; Van Der Kamp, G. Simple equations to represent the volume-area-depth relations of shallow wetlands in small topographic depressions. J. Hydrol. 2000, 237, 74–85. [Google Scholar] [CrossRef]
  21. Los Huertos, M.; Smith, D. Wetland bathymetry and mapping. In Wetland Techniques. Vol. 1: Foundations; Anderson, J.T., Davis, C.A., Eds.; Springer: Dordrecht, The Netherlands, 2013; pp. 49–86. [Google Scholar]
  22. Ray, A.M.; Sepulveda, A.J.; Irvine, K.M.; Wilmoth, S.K.C.; Thoma, D.P.; Patla, D.A. Wetland drying linked to variations in snowmelt runoff across Grand Teton and Yellowstone national parks. Sci. Total Environ. 2019, 666, 1188–1197. [Google Scholar] [CrossRef]
  23. Keijman, J.Q.; Koopmans, R.W.R. A comparison of several methods of estimating the evaporation of Lake Flevo. Int. Assoc. Hydrol. Sci. 1973, 109, 225–232. [Google Scholar]
  24. Rosenberry, D.O.; Winter, T.C.; Buso, D.C.; Likens, G.E. Comparison of 15 evaporation methods applied to a small mountain lake in the northeastern USA. J. Hydrol. 2007, 340, 149–166. [Google Scholar] [CrossRef]
  25. Yao, H. Long-Term Study of Lake Evaporation and Evaluation of Seven Estimation Methods: Results from Dickie Lake, South-Central Ontario, Canada. J. Water Resour. Prot. 2009, 1, 59–77. [Google Scholar] [CrossRef] [Green Version]
  26. McMahon, T.A.; Peel, M.C.; Lowe, L.; Srikanthan, R.; McVicar, T.R. Estimating actual, potential, reference crop and pan evaporation using standard meteorological data: A pragmatic synthesis. Hydrol. Earth Syst. Sci. 2013, 17, 1331–1363. [Google Scholar] [CrossRef] [Green Version]
  27. Rodríguez-Rodríguez, M.; Moral, F.; Benavente, J.; Beltrán, M. Developing hydrological indices in semi-arid playa-lakes by analyzing their main morphometric, climatic and hydrochemical characteristics. J. Arid Environ. 2010, 74, 1478–1486. [Google Scholar] [CrossRef]
  28. Vera, J.A.; Martín-Algarra, A. Cordillera Bética. In Geología de España; Vera, J.A., Ed.; IGME-SGE: Madrid, Spain, 2004; pp. 345–464. [Google Scholar]
  29. CHG. Definición del Contexto Hidrogeológico de Humedales de la Campiña Andaluza Central; Confederación Hidrográfica del Guadalquivir: Seville, Spain, 2008. [Google Scholar]
  30. CMA. Definición del Contexto Hidrogeológico de Humedales Andaluces; Conserjería de Medio Ambiente de la Junta de Andalucía: Málaga, Spain, 2005. [Google Scholar]
  31. Aljibe-Consultores. Estudio Hidrogeológico de las Lagunas del sur de Córdoba: Lagunas de Rincón y Santiago, Lagunas Amarga y Dulce; Junta de Andalucía: Córdoba, Spain, 2005. [Google Scholar]
  32. Moral, F.; Rodríguez-Rodríguez, M.; Beltrán, M.; Benavente, J.; Cifuentes, V.J. Water Regime of Playa Lakes from Southern Spain: Conditioning Factors and Hydrological Modeling. Water Environ. Res. 2013, 85, 632–642. [Google Scholar] [CrossRef]
  33. Martín-Serrano, A. Mapa Geológico de la Hoja no 1006 (Benamejí). Mapa Geológico de España, E. 1:50.000. Segunda Serie (MAGNA), Primera Edición; IGME: Madrid, Spain, 1986. [Google Scholar]
  34. Gil-Márquez, J.M.; Andreo, B.; Mudarra, M. Combining hydrodynamics, hydrochemistry, and environmental isotopes to understand the hydrogeological functioning of evaporite-karst springs. An example from southern Spain. J. Hydrol. 2019, 576, 299–314. [Google Scholar] [CrossRef]
  35. Gil-Márquez, J.M. Caracterización Hidrogeológica de Humedales y Manantiales Salinos Asociados a Acuíferos Kársticos Evaporíticos del Sector Central del Subbético; Universidad de Granada: Granada, Spain, 2018. [Google Scholar]
  36. Gil-Márquez, J.M.; Barberá, J.A.; Mudarra, M.; Andreo, B.; Prieto-Mera, J.; Sánchez, D.; Rizo-Decelis, L.D.; Argamasilla, M.; Nieto, J.M.; De la Torre, B. Karst development of an evaporitic system and its hydrogeological implications inferred from GIS-based analysis and tracing techniques. Int. J. Speleol. 2017, 46, 219–235. [Google Scholar] [CrossRef] [Green Version]
  37. Miller, V.C. A Quantitative Geomorphic Study of Drainage Basin Characteristics in the Clinch Mountain Area Virginia and Tennessee (Proj. NR 389-402, Technical Report 3); Columbia University: New York, NY, USA, 1953. [Google Scholar]
  38. Schumm, S.A. Geological Society of America Bulletin Evolution of Drainage Systems and Slopes in Badlands at Perth Amboy, New Jersey. Geol. Soc. Am. Bull. 1956, 67, 597–646. [Google Scholar] [CrossRef]
  39. Horton, R.E. Drainage-basin characteristics. Trans. Am. Geophys. Union 1932, 13, 350. [Google Scholar] [CrossRef]
  40. IGN Plan Nacional de Ortofotografía Aerea. Available online: https://pnoa.ign.es/el-proyecto-pnoa-lidar (accessed on 1 March 2021).
  41. García-Ferrer, A.; Recio Espejo, J.M.; Sánchez de la Orden, M. Medidas morfométricas y batimétricas de las lagunas Amarga y del Rincón (Córdoba). Boletín la Estac. Cent. Ecol. 1983, 12, 51–55. [Google Scholar]
  42. Voldseth, R.A.; Johnson, W.C.; Gilmanov, T.; Guntenspergen, G.R.; Millett, B.V. Model estimation of land-use effects on water levels of northern prairie wetlands. Ecol. Appl. 2007, 17, 527–540. [Google Scholar] [CrossRef]
  43. Gurrieri, J.T.; Furniss, G. Estimation of groundwater exchange in alpine lakes using non-steady mass-balance methods. J. Hydrol. 2004, 297, 187–208. [Google Scholar] [CrossRef]
  44. Krasnostein, A.L.; Oldham, C.E. Predicting wetland water storage. Water Resour. Res. 2004, 40, 1–12. [Google Scholar] [CrossRef]
  45. Dansgaard, W. Stable isotopes in precipitation. Tellus 1964, 16, 436–468. [Google Scholar] [CrossRef]
  46. Mockus, V. Estimation of direct runoff from storm rainfall. In SCS, National Engineering Handbook; USDA Soil Conservation Services: Washington DC, USA, 1964; p. 30. [Google Scholar]
  47. Públicas, M.d.O. Cálculo Hidrometeorológico de Caudales Máximos en Pequeñas Cuencas Naturales; Ministerio de Obras Públicas: Madrid, Spain, 1987. [Google Scholar]
  48. Thornthwaite, C.W.; Mather, J.R. The water balance. Publ. Climatol. 1955, 8, 86. [Google Scholar]
  49. Thornthwaite, C.W. An approach toward rational classification of climate. Geogr. Rev. 1948, 38, 55–94. [Google Scholar] [CrossRef]
  50. Hargreaves, G.H.; Samani, Z.A. Reference crop evapotranspiration from temperature. Appl. Eng. Agric. 1985, 1, 96–99. [Google Scholar] [CrossRef]
  51. Blaney, H.F.; Criddle, W.D. Determining Water Requirements in Irrigated Areas from Climatological and Irrigation Data. Technical Report-96; USDA: Washington, DC, USA, 1950. [Google Scholar]
  52. Penman, H.L. Natural evaporation from open water, bare soil and grass. Proc. R. Soc. Lond. A 1948, 193, 120–145. [Google Scholar]
  53. Souch, C.; Grimmond, C.S.B.; Wolfe, C.P. Evapotranspiration rates from wetlands with different disturbance histories: Indiana Dunes National Lakeshore. Wetlands 1998, 18, 216–229. [Google Scholar] [CrossRef]
  54. Finch, J.W. A comparison between measured and modelled open water evaporation from a reservoir in south-east England. Hydrol. Process. 2001, 15, 2771–2778. [Google Scholar] [CrossRef]
  55. Wossenu, A. Evaporation Estimation for Lake Okeechobee in South Florida. J. Irrig. Drain. Eng. 2001, 127, 140–147. [Google Scholar] [CrossRef]
  56. Rosenberry, D.O.; Morin, R.H. Use of an Electromagnetic Seepage Meter to Investigate Temporal Variability in Lake Seepage. Ground Water 2004, 42, 68–77. [Google Scholar] [CrossRef]
  57. Kohler, M.A.; Parmele, L.H. Generalized estimates of free-water evaporation. Water Resour. Res. 1967, 3, 997–1005. [Google Scholar] [CrossRef]
  58. Monteith, J.L. Evaporation and surface temperature. Q. J. R. Meteorol. Soc. 1981, 107, 1–27. [Google Scholar] [CrossRef]
  59. Morton, F.I. Practical estimates of lake evaporation. J. Clim. Appl. Meteorol. 1986, 25, 371–387. [Google Scholar] [CrossRef]
  60. Fennessey, N.M. Estimating average monthly lake evaporation in tile northeast united states. J. Am. Water Resour. Assoc. 2000, 36, 759–769. [Google Scholar] [CrossRef]
  61. McJannet, D.L.; Webster, I.T.; Stenson, M.P.; Sherman, B.S. Estimating Open Water Evaporation for the Murray-Darling Basin: A Report to the Australian Government from the CSIRO Murray-Darling Basin Sustainable Yields Project; CSIRO: Canbera, Australia, 2008; Volume 50. [Google Scholar]
  62. Vardavas, I.M.; Fountoulakis, A. Estimation of lake evaporation from standard meteorological measurements: Application to four Australian lakes in different climatic regions. Ecol. Modell. 1996, 84, 139–150. [Google Scholar] [CrossRef]
  63. Allen, R.G.; Pereira, L.S.; Raes, D.; Smith, M. Crop Evapotranspiration—Guidelines for Computing Crop Water Requirements; FAO Irrigation and Drainage Paper 56; FAO: Rome, Italy, 1998. [Google Scholar]
  64. Finch, J.; Calver, A. Methods for the Quantification of Evaporation from Lakes. Available online: http://nora.nerc.ac.uk/id/eprint/14359/1/wmoevap_271008.pdf (accessed on 1 March 2021).
  65. Andreo, B.; Gil-Márquez, J.M.; Mudarra, M.; Linares, L.; Carrasco, F. Hypothesis on the hydrogeological context of wetland areas and springs related to evaporitic karst aquifers (Málaga, Córdoba and Jaén provinces, Southern Spain). Environ. Earth Sci. 2016, 75, 1–19. [Google Scholar] [CrossRef]
  66. Doss, P.K. The nature of a dynamic water table in a system of non-tidal, freshwater coastal wetlands. J. Hydrol. 1993, 141, 107–126. [Google Scholar] [CrossRef]
  67. Moya, J.L. Hidrología de las Zonas Húmedas del sur de la Provincia de Córdoba; Universidad de Córdoba: Córdoba, Spain, 1988. [Google Scholar]
  68. Gil-Márquez, J.M.; Sültenfuß, J.; Andreo, B.; Mudarra, M. Groundwater dating tools (3H, 3He, 4He, CFC-12, SF6) coupled with hydrochemistry to evaluate the hydrogeological functioning of complex evaporite-karst settings. J. Hydrol. 2020, 580, 124263. [Google Scholar] [CrossRef]
  69. Valiente, N.; Gil-Márquez, J.M.; Gómez-Alday, J.J.; Andreo, B. Unraveling groundwater functioning and nitrate attenuation in evaporitic karst systems from southern Spain: An isotopic approach. Appl. Geochem. 2020, 123, 104820. [Google Scholar] [CrossRef]
  70. Cameron, C.R.; Hancock, M.C.; Carr, D.W.; Hurst, M.K.; Campbell, D.E.; Venning, T.J.; Tara, P.D.; Holzwart, K.R. Hydroperiods of Cypress Domes in West-Central Florida, USA. Wetlands 2020, 40, 2225–2234. [Google Scholar] [CrossRef]
  71. Montgomery, J.S.; Hopkinson, C.; Brisco, B.; Patterson, S.; Rood, S.B. Wetland hydroperiod classification in the western prairies using multitemporal synthetic aperture radar. Hydrol. Process. 2018, 32, 1476–1490. [Google Scholar] [CrossRef]
  72. Hayashi, M.; van der Kamp, G.; Rosenberry, D.O. Hydrology of Prairie Wetlands: Understanding the Integrated Surface-Water and Groundwater Processes. Wetlands 2016, 36, 237–254. [Google Scholar] [CrossRef]
  73. Tóth, Á.; Havril, T.; Simon, S.; Galsa, A.; Monteiro Santos, F.A.; Müller, I.; Mádl-Szőnyi, J. Groundwater flow pattern and related environmental phenomena in complex geologic setting based on integrated model construction. J. Hydrol. 2016, 539, 330–344. [Google Scholar] [CrossRef]
  74. Manzano, M.; Custodio, E. Relationships between weltands and the Doñana coastal aquifer (SW Spain). In Groundwater and Ecosystems; IAH, Ed.; CRC Press: London, UK, 2013; pp. 169–182. ISBN 978-1-138-00033-9. [Google Scholar]
  75. Bryant, R.G.; Rainey, M.P. Investigation of flood inundation on playas within the Zone of Chotts, using a time-series of AVHRR. Remote Sens. Environ. 2002, 82, 360–375. [Google Scholar] [CrossRef]
  76. Brooks, R.T. Potential impacts of global climate change on the hydrology and ecology of ephemeral freshwater systems of the forests of the northeastern United States. Clim. Chang. 2009, 95, 469–483. [Google Scholar] [CrossRef]
  77. Ma, J.Z.; Ding, Z.; Gates, J.B.; Su, Y. Chloride and the environmental isotopes as the indicators of the groundwater recharge in the Gobi Desert, northwest China. Environ. Geol. 2008, 55, 1407–1419. [Google Scholar] [CrossRef] [Green Version]
  78. Su, M.; Stolte, W.J.; Van Der Kamp, G. Modelling Canadian prairie wetland hydrology using a semi-distributed streamflow model. Hydrol. Process. 2000, 14, 2405–2422. [Google Scholar] [CrossRef]
  79. Hayashi, M.; Rosenberry, D.O. Effects of Ground Water Exchange on the Hydrology and Ecology of Surface Water. Ground Water 2002, 40, 309–316. [Google Scholar] [CrossRef]
  80. Acreman, M. Wetlands and Hydrology; Crivelli, J.A., Albert, J., Eds.; Station Biologuique de la Tour du Valat: Arles, France, 2000; ISBN 2-910368-32-7. [Google Scholar]
  81. Euliss, N.H.; Labaugh, J.W.; Fredrickson, L.H.; Mushet, D.M.; Laubhan, M.K.; Swanson, G.A.; Winter, T.C.; Rosenberry, D.O.; Nelson, R.D. the Wetland Continuum: A Conceptual Framework for Interpreting Biological Studies. Wetlands 2004, 24, 448–458. [Google Scholar] [CrossRef]
  82. Smith, A.J.; Townley, L.R. Influence of regional setting on the interaction between shallow lakes and aquifers. Water Resour. Res. 2002, 38. [Google Scholar] [CrossRef]
  83. Neff, B.P.; Rosenberry, D.O.; Leibowitz, S.G.; Mushet, D.M.; Golden, H.E.; Rains, M.C.; Brooks, J.R.; Lane, C.R. A Hydrologic Landscapes Perspective on Groundwater Connectivity of Depressional Wetlands. Water 2019, 12, 50. [Google Scholar] [CrossRef] [Green Version]
Figure 2. Evolution of annual rainfall and its accumulated deviation in the Jarales area (recording period 1963/64 to 2016/2017). Modified from Gil-Márquez [35]. References to wet, intermediate, and dry hydrological years are indicated.
Figure 2. Evolution of annual rainfall and its accumulated deviation in the Jarales area (recording period 1963/64 to 2016/2017). Modified from Gil-Márquez [35]. References to wet, intermediate, and dry hydrological years are indicated.
Water 13 01482 g002
Figure 4. Groundwater level evolution in wells near the Jarales and Amarga wetlands (May 2015–March 2016). Figure 1 shows the location of the wells (in red). The evolution of the Amarga and Jarales stages is also shown.
Figure 4. Groundwater level evolution in wells near the Jarales and Amarga wetlands (May 2015–March 2016). Figure 1 shows the location of the wells (in red). The evolution of the Amarga and Jarales stages is also shown.
Water 13 01482 g004
Figure 6. The best four limnimetric simulations of the Jarales wetland (March 2014–September 2017) and the daily wetland stage and rainfall records. JW1, CN runoff with an Ia of 19 mm and Penman evaporation; JW2, CN runoff with an Ia of 19 mm and Penman–Monteith evaporation; JW3, SWB runoff with Blaney–Criddle ETP, AWC of 100 mm, and Penman evaporation; J4, SWB runoff with Blaney–Criddle ETP, AWC of 75 mm, 50% RC, and Penman evaporation.
Figure 6. The best four limnimetric simulations of the Jarales wetland (March 2014–September 2017) and the daily wetland stage and rainfall records. JW1, CN runoff with an Ia of 19 mm and Penman evaporation; JW2, CN runoff with an Ia of 19 mm and Penman–Monteith evaporation; JW3, SWB runoff with Blaney–Criddle ETP, AWC of 100 mm, and Penman evaporation; J4, SWB runoff with Blaney–Criddle ETP, AWC of 75 mm, 50% RC, and Penman evaporation.
Water 13 01482 g006
Figure 7. The best four limnimetric simulations of the Amarga wetland (October 2014–September 2017) and the daily wetland stage and rainfall records. AW1, CN runoff with an Ia of 15 mm and Penman evaporation; AW2, CN runoff with an Ia of 15 mm and Penman–Monteith evaporation; AW3, CN runoff with an Ia of 15 mm and Penman evaporation corrected by Vardavas and Fountoulakis; J4, SWB runoff with Blaney–Criddle ETP, AWC of 50 mm, 75% RC, and Penman evaporation corrected by Vardavas and Fountoulakis.
Figure 7. The best four limnimetric simulations of the Amarga wetland (October 2014–September 2017) and the daily wetland stage and rainfall records. AW1, CN runoff with an Ia of 15 mm and Penman evaporation; AW2, CN runoff with an Ia of 15 mm and Penman–Monteith evaporation; AW3, CN runoff with an Ia of 15 mm and Penman evaporation corrected by Vardavas and Fountoulakis; J4, SWB runoff with Blaney–Criddle ETP, AWC of 50 mm, 75% RC, and Penman evaporation corrected by Vardavas and Fountoulakis.
Water 13 01482 g007
Figure 8. Hydrogeological conceptual model of the Jarales and Amarga wetlands. AW, Amarga wetland; JW, Jarales wetland.
Figure 8. Hydrogeological conceptual model of the Jarales and Amarga wetlands. AW, Amarga wetland; JW, Jarales wetland.
Water 13 01482 g008
Table 1. Morphological information of the wetland watersheds obtained from the DTM.
Table 1. Morphological information of the wetland watersheds obtained from the DTM.
JARALES WETLAND
ShapeSlopeHeight
Surface90.4 haCircularity Ratio (Rc) a0.33Maximum11.4%Maximum464.5 m asl
Perimeter5857 mElongation Ratio (Re) b0.62Average3.6%Minimum406.6 m asl
Maximum Length1736 mForm Factor (Ff) c0.30 Difference50.6 m
AMARGA WETLAND
ShapeSlopeHeight
Surface256.9 haCircularity Ratio (Rc) a0.34Maximum30.2%Maximum440.4 m asl
Perimeter9730 mElongation Ratio (Re) b0.79Average3.8%Minimum370.1 m asl
Maximum Length2280 mForm Factor (Ff) c0.49 Difference70.3 m
a Miller’s circularity ratio [37]: Rc = 4πA/P2, where A = area of the basin (m2), P = perimeter (m). b Schumm’s elongation ratio [38]: Re = D/L = 1.128√A, where D = diameter of a circle of the same area (A) as the basin, A = area of the basin (m2), L = basin length (m). c Horton’s form factor [39]: Ff = A/L2, where A = area of the basin (m2), L = basin length (km).
Table 2. Monthly and annual water budget of the Jarales wetland in the study period. Units in dam3.
Table 2. Monthly and annual water budget of the Jarales wetland in the study period. Units in dam3.
Month/Hydrol. YearV0aPbRcEdSimulated Vf eSimulated ΔV fActual V fActual ΔVΔG g
Monthly water budget
May 201494.80.40.014.680.6−14.277.6−17.2−3.0
June 201477.60.70.014.563.8−13.861.2−16.3−2.5
July 201461.20.00.014.546.7−14.542.4−18.8−4.4
August 201442.40.00.011.231.2−11.227.0−15.4−4.2
September 201427.00.30.05.422.0−5.119.2−7.8−2.7
October 201419.22.51.13.519.30.114.4−4.8−4.9
November 201414.44.80.51.618.13.720.86.42.7
December 201420.80.80.01.220.4−0.519.8−1.0−0.6
January 201519.81.70.01.520.00.222.02.22.0
February 201522.01.20.02.121.1−0.922.40.41.4
March 201522.41.40.04.119.7−2.822.0−0.42.3
April 201522.01.20.05.517.6−4.418.2−3.80.6
May 201518.20.00.07.410.8−7.412.2−5.91.5
June 201512.20.00.06.85.5−6.81.0−11.2−4.4
July 20151.00.00.01.00.0−1.00.0−1.00.0
August 20150.00.00.00.00.00.00.00.00.0
September 20150.00.00.00.00.00.00.00.00.0
Annual water budget
2014/1519.213.51.634.90−19.80−19.20.6
a Initial volume of water stored in the wetland. b Total volume of precipitation inputs. c Total volume of runoff inputs. d Total volume of evaporation outputs. e Final volume of water stored in the wetland. f Change of the volume of water stored in the wetland. g Volume of the net groundwater–wetland exchange.
Table 3. Monthly and annual water budgets of the Amarga wetland in the study period. Units in dam3.
Table 3. Monthly and annual water budgets of the Amarga wetland in the study period. Units in dam3.
Month/Hydrol. YearV0 aPbRcEdSimulated Vf eSimulated ΔV fActual V fActual ΔVΔG g
Monthly water budget
December 2014187.981.100.002.40186.69−1.30187.13−0.850.45
January 2015187.132.460.002.94186.65−0.48188.581.441.93
February 2015188.581.610.002.60187.59−0.99188.620.041.03
March 2015188.621.960.004.63185.95−2.67188.25−0.362.31
April 2015188.251.690.006.69183.26−5.00187.57−0.684.31
May 2015187.570.000.0010.52177.05−10.52177.12−10.450.07
June 2015177.120.020.0012.01165.14−11.99163.27−13.85−1.86
July 2015163.270.000.0013.02150.25−13.02149.39−13.89−0.87
August 2015149.390.050.0010.42139.02−10.36138.19−11.20−0.83
September 2015138.190.950.007.84131.30−6.89127.93−10.26−3.37
October 2015127.933.330.004.67126.59−1.34127.83−0.101.24
November 2015127.832.384.662.75132.124.29135.727.883.59
December 2015135.720.130.002.75133.10−2.62133.78−1.940.68
January 2016133.784.015.162.03140.937.15145.2411.474.32
February 2016145.242.100.002.83144.52−0.72146.040.801.52
March 2016146.041.220.004.42142.85−3.19143.49−2.550.65
April 2016143.492.852.674.82144.190.69145.131.640.94
May 2016145.133.270.587.61141.37−3.76143.39−1.752.01
June 2016143.390.000.0010.97132.42−10.97131.76−11.63−0.67
July 2016131.760.810.0012.05120.51−11.24123.72−8.043.21
August 2016123.720.020.0011.67112.07−11.65112.98−10.740.91
September 2016112.980.180.007.62105.55−7.43108.07−4.912.52
October 2016108.071.590.005.15104.52−3.55109.591.525.08
November 2016109.593.301.012.78111.121.53121.9212.3310.80
December 2016121.922.390.132.24122.190.27128.826.906.63
January 2017128.820.300.002.25126.86−1.96124.57−4.24−2.29
February 2017124.572.040.002.09124.52−0.06127.042.462.52
March 2017127.042.190.003.68125.54−1.49127.640.612.10
April 2017127.641.070.006.19122.52−5.12122.34−5.31−0.18
May 2017122.340.700.007.77115.27−7.06116.58−5.761.31
June 2017116.580.000.0010.22106.36−10.22105.72−10.86−0.64
July 2017105.720.010.009.5996.14−9.5895.48−10.24−0.66
August 201795.481.421.148.1189.93−5.5596.200.716.27
September 201796.200.000.006.6089.60−6.6091.66−4.542.06
Annual water budget
2015/16127.920.313.174.287.1−40.8108.1−19.920.9
2016/17108.115.02.366.758.7−49.491.7−16.433.0
a Initial volume of water stored in the wetland. b Total volume of precipitation inputs. c Total volume of runoff inputs. d Total volume of evaporation outputs. e Final volume of water stored in the wetland. f Change of the volume of water stored in the wetland. g Volume of the net groundwater–wetland exchange.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Gil-Márquez, J.M.; Andreo, B.; Mudarra, M. Comparative Analysis of Runoff and Evaporation Assessment Methods to Evaluate Wetland–Groundwater Interaction in Mediterranean Evaporitic-Karst Aquatic Ecosystem. Water 2021, 13, 1482. https://0-doi-org.brum.beds.ac.uk/10.3390/w13111482

AMA Style

Gil-Márquez JM, Andreo B, Mudarra M. Comparative Analysis of Runoff and Evaporation Assessment Methods to Evaluate Wetland–Groundwater Interaction in Mediterranean Evaporitic-Karst Aquatic Ecosystem. Water. 2021; 13(11):1482. https://0-doi-org.brum.beds.ac.uk/10.3390/w13111482

Chicago/Turabian Style

Gil-Márquez, José Manuel, Bartolomé Andreo, and Matías Mudarra. 2021. "Comparative Analysis of Runoff and Evaporation Assessment Methods to Evaluate Wetland–Groundwater Interaction in Mediterranean Evaporitic-Karst Aquatic Ecosystem" Water 13, no. 11: 1482. https://0-doi-org.brum.beds.ac.uk/10.3390/w13111482

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