An Integrated Approach for Modeling Wetland Water Level: Application to a Headwater Wetland in Coastal Alabama, USA

Headwater wetlands provide many benefits such as water quality improvement, water storage, and providing habitat. These wetlands are characterized by water levels near the surface and respond rapidly to rainfall events. Driven by both groundwater and surface water inputs, water levels (WLs) can be above or below the ground at any given time depending on the season and climatic conditions. Therefore, WL predictions in headwater wetlands is a complex problem. In this study a hybrid modeling approach was developed for improved WL predictions in wetlands, by coupling a watershed model with artificial neural networks (ANNs). In this approach, baseflow and stormflow estimates from the watershed draining to a wetland are first estimated using an uncalibrated Soil and Water Assessment Tool (SWAT). These estimates are then combined with meteorological variables and are utilized as inputs to an ANN model for predicting daily WLs in wetlands. The hybrid model was used to successfully predict WLs in a headwater wetland in coastal Alabama, USA. The model was then used to predict the WLs at the study wetland from 1951 to 2005 to explore the possible teleconnections between the El Nino Southern Oscillation (ENSO) and WLs. Results show that both precipitation and the variations in WLs are partially affected by ENSO in the study area. A correlation analysis between seasonal precipitation and the Nino 3.4 Index suggests that winters are wetter during El Nino in Coastal Alabama. Analysis also revealed a significant negative correlation between WLs and the Nino 3.4 Index during the El Nino phase for spring. The findings of this study and the developed methodology/tools are useful to predict long-term WLs in wetlands and construct more accurate restoration plans under a variable climate.


Introduction
Headwater wetlands are among the most important wetland types in the coastal plains of the southeastern U.S. Since they are often the source of first and second order perennial streams, they act as natural filters for water quality improvements and therefore are crucial in protecting downstream aquatic systems [1]. Headwater wetlands are characterized by water tables at or near the surface that respond rapidly to rainfall events [2]. Driven by both groundwater and surface water inputs, water levels (WLs) can be above or below the ground at any given time, depending on the season and climatic conditions [3]; this dynamic behavior results in a highly complex ecosystem. The seasonal patterns of WLs in a wetland (i.e., hydroperiod) exhibit annual variation in response to climatic and developed a hybrid modeling structure (by utilizing Soil and Water Assessment Tool SWAT and ANNs) to predict runoff generation from prairie landscapes.
We develop a coupled SWAT-ANN model in this paper to predict WLs in headwater wetlands. We then apply the model to a headwater wetland in coastal Alabama, USA to predict the WL at the deepest point of the wetland. The study wetland is a typical headwater wetland in Coastal Alabama and is representative of many in the region since this wetland is one of those 15 selected by Barksdale et al. [19], to properly represent the surrounding LULC of the region. Also, groundwater discharge (from surrounding uplands) is a major hydrologic input to these headwater wetlands [2]. Baseflow and stormflow estimates from the SWAT model were fed into the ANN model as forcing data to capture slow and quick flow components from the contributing watershed to the study wetland. Observed WL data were then used to train and test the ANN model. The coupled model of this paper overcomes some of the limitations of the two previously developed WL prediction models by Rezaeianzadeh et al. [3], which either required WL data from nearby wetlands (continuous model) or antecedent WL data from the target wetland itself (event based model) as input. None of those models are suited for long term predictions, which is needed to explore the impacts of climate variability on wetland hydrology. The utility of the developed model is demonstrated with an application, where we explored the potential teleconnections between El Niño Southern Oscillation (ENSO) and WL fluctuations in the study wetland. To the best of our knowledge, this study is the first attempt to look at the hydrologic function of a headwater wetland in Coastal Alabama under varying ENSO conditions. Although the WL data used to train and test the model was limited, the example application still provides some useful insights into the impact of climatic variability on WLs in the study wetland.

Study Area and Data Sets
The study wetland is located in Baldwin County, Alabama, USA ( Figure 1). Mild winters and hot and humid summers are the main climatic characteristics of the region. Mean annual temperature and precipitation are 19.5 • C and 1650 mm, respectively. Headwater wetlands in coastal Alabama are groundwater driven and are usually located in the headwater reaches of first order streams [2]. Hummock/hollow microtopography and gradual slopes are the terrain characteristics. The soils of these wetlands are classified as "wet loamy alluvial lands" [26]. The study wetland is a representative of the headwater wetlands throughout southern Alabama. The soil consists of two main groups: Lv (i.e., local alluvial land; C) and Hb (i.e., Hyde, Bayboro, and Muck soils; A/D) [27]. Figure 1a shows the 15 headwater wetlands that have been originally studied by . The minimum and maximum elevation of the study wetland is 7.12 m and 10.24 m, respectively. The study wetland is located within the Fish River watershed in the southwest portion of the county, which drains into Weeks Bay, a sub-estuary of Mobile Bay. The delineated watershed area (based on the outlet of the wetland) draining to the wetland is 83.57 ha with a mixed land use/cover (Figure 1b). Pasture/hay and evergreen forests are the dominant land uses in the watershed (Figure 1c). The dominant hydrologic soil groups are C (i.e., slow infiltration rate covering 41.3% of the area) and A (with high infiltration rate and 40% of the area), respectively. According to the information derived from National Wetlands Inventory (NWI), the study wetland is comprised of three sections, which includes the PFO1C (0.737 ha), PFO4C (0.485 ha), and PFO1C (0.314 ha) classification codes. These three sections are called Freshwater Forested/Shrub Wetland types by NWI and are classified as inland wetlands. Precipitation, minimum and maximum temperature data were obtained from the NOAA Robertsdale (USC00016988) climatic station. Wetland WLs (relative to the ground surface) were monitored at an hourly time scale at the deepest point of the wetland in shallow wells using In-Situ Mini-Troll 500 pressure transducers from February 2011 to March 2012 [19]. During the summer and fall (June through November) of 2011, drought conditions persisted across the study area with the driest period occurring from May-June [28]. However, there was no WL data available for the study wetland from 25 April 2011 to 17 July 2011 ( Figure 2). Based on the long-term average of precipitation, five months including March, July, August, and September 2011 as well as February 2012 were above the average and the rest were below the average precipitation.

Soil and Water Assessment Tool (SWAT)
SWAT is a semi-distributed, process-based watershed model that operates on a daily basis [29]. Division of a watershed into a number of subbasins is needed and done when different parts of the watershed are dominated by land uses or soils dissimilar enough to impact hydrology [30]. Hydrologic response units (HRUs) are lumped land areas within the subbasin that consist of unique land cover, soil, and management combinations [30]. In this study, the thresholds for land use, soil, and slope overlay, were all set to 10%, and the number of subbasins and HRUs were 5 and 74, respectively. The Soil Survey Geographic Database (SSURGO, [31]) was used to extract soil parameters. Land use/cover related data came from the 2011 National Land Cover Database (NLCD; [32]). Lu et al. [33] recommended Hamon method [34] for regional applications of potential evapotranspiration (PET) estimations in the southeastern United States. Thus, the Hamon method was used for the PET calculation, which utilizes daily mean air temperature. Daily average temperatures were obtained by simply averaging daily minimum and maximum temperatures. A three-year warm-up period was considered for SWAT run, and the daily baseflow and stormflow estimates generated from SWAT (no calibration) were regarded as two hydrologic inputs to the ANN model. Observed flow data (outflow from the watershed or inflow to the wetland) were not available for model calibration, which is commonly the case for many natural wetlands.

Artificial Neural Networks (ANNs)
An ANN consists of a number of neurons organized into three basic layers (input, hidden, and output). The multilayer feed forward is one of the most common ANNs and it is therefore utilized in this study. Input combinations, the number of neurons in the hidden layer, transfer function, and training algorithms are the most important components of an ANN structure. The input combinations and transfer function are discussed in the following sections. For this study, average daily WLs were calculated from the hourly WL data, available from February 2011 to March 2012; we Water level recorded at the study wetland; there was no water level (WL) data available for the study wetland from 25 April 2011 to 17 July 2011. Note that there were two other periods of missing data for the study site. Altogether, water level data were available for 279 days.

Soil and Water Assessment Tool (SWAT)
SWAT is a semi-distributed, process-based watershed model that operates on a daily basis [29]. Division of a watershed into a number of subbasins is needed and done when different parts of the watershed are dominated by land uses or soils dissimilar enough to impact hydrology [30]. Hydrologic response units (HRUs) are lumped land areas within the subbasin that consist of unique land cover, soil, and management combinations [30]. In this study, the thresholds for land use, soil, and slope overlay, were all set to 10%, and the number of subbasins and HRUs were 5 and 74, respectively. The Soil Survey Geographic Database (SSURGO, [31]) was used to extract soil parameters. Land use/cover related data came from the 2011 National Land Cover Database (NLCD; [32]). Lu et al. [33] recommended Hamon method [34] for regional applications of potential evapotranspiration (PET) estimations in the southeastern United States. Thus, the Hamon method was used for the PET calculation, which utilizes daily mean air temperature. Daily average temperatures were obtained by simply averaging daily minimum and maximum temperatures. A three-year warm-up period was considered for SWAT run, and the daily baseflow and stormflow estimates generated from SWAT (no calibration) were regarded as two hydrologic inputs to the ANN model. Observed flow data (outflow from the watershed or inflow to the wetland) were not available for model calibration, which is commonly the case for many natural wetlands.

Artificial Neural Networks (ANNs)
An ANN consists of a number of neurons organized into three basic layers (input, hidden, and output). The multilayer feed forward is one of the most common ANNs and it is therefore utilized in this study. Input combinations, the number of neurons in the hidden layer, transfer function, and training algorithms are the most important components of an ANN structure. The input combinations and transfer function are discussed in the following sections. For this study, average daily WLs were calculated from the hourly WL data, available from February 2011 to March 2012; we obtained the best architecture to be (8, 6, 1). A simple trial and error procedure was carried out to find the best training algorithm and the Levenberg-Marquardt training algorithm was chosen for this study. A momentum constant (0 to 1) was used to avoid getting stuck in a local minimum [35]. Before applying the ANN models, the observed data were normalized to [0.05, 0.95] using a linear transformation [3,36]: where, X n and X are the normalized and original inputs and X min and X max are the minimum and maximum values of the original data, respectively. The normalized data were employed to train the ANN model. The model outputs were then transformed back to the original scale and then the performance indices were computed for the training and test data sets. In this study, roughly 70% of the data was used for training and the remaining 30% was used for testing. To distribute the high, mean, and low WL values into both training and testing datasets, the data were randomly selected [3]. This randomization eliminates potential biases that may arise when there are distinct dry and wet periods. The randomization is done by randomly distributing the data into training and testing datasets until the p-value in the t-test is above 5%; i.e., there is no statistically significant difference between the means from the training and testing datasets. Note that the t-test is performed based on the equal variances assumption by conducting a Levene's test. The methodology is described in detail in Rezaeianzadeh et al. (2015) [3], therefore it is not repeated here. Previous exercises by Rezaeian-Zadeh et al. [36] and Rezaeianzadeh et al. [3] showed the efficiency of this method to optimally train the ANNs. The codes were written in MATLAB for the ANN simulation [37].
Although it is generally accepted that ANNs cannot extrapolate beyond the range of the training data [38], utilization of sigmoidal-type and linear transfer functions in hidden and output layers, respectively, has been recommended for extrapolation purposes [39,40]. For instance, Cigizoglu [41], and Rezaeian-Zadeh et al. [42], respectively showed the extrapolation ability of ANNs in predicting daily streamflows and hourly air temperatures. Hence, in the current study, logistic sigmoid and linear transfer functions were utilized in hidden and output layers, respectively. The logistic sigmoid transfer function is defined for a variable s as: We utilized sigmoid and linear transfer functions for the hidden and output layers to build a more generalized model capable of capturing variations on unseen datasets including extrapolations. However, we need a longer dataset to increase the reliability of the developed model. Note that data limitation should not result in overlooking the capability of the introduced methodology in this study.

SWAT-ANN Coupling
To enhance WL predictions by the ANN model, a coupled SWAT-ANN model was developed. To that end, baseflow and stormflow components of streamflow are first simulated at the watershed outlet by the SWAT model without calibration (i.e., default SWAT parameter values are used). Note that even if they have a well-defined inlet, these wetlands typically receive surface and subsurface water from various points between the wetland inlet and outlet. Those dominantly un-regulated inflows to the wetlands impact the hydrology of these ecosystems in that area substantially. Hence, we decided to consider the outlet point of the wetland instead of the inlet point for SWAT analysis. This is achieved by picking the outlet of the wetland and delineating the watershed draining to that point as if there is no wetland there. Then, we generated streamflow at that point with SWAT using default model parameters. No calibration was carried out because streamflow data was not available. In addition to streamflow, SWAT also provides model outputs to help partition the streamflow into stormflow and baseflow components [30]. We did not separate baseflow from streamflow. SWAT calculates baseflow, subsurface flow and surface runoff components. We considered surface runoff and subsurface flow as stormflow. We note that in many instances the wetland WLs are below the ground surface whereby groundwater bypasses the wetland and discharges directly to the main channel as baseflow. The simulated values for baseflow and stormflow at day t are then fed into ANN to predict WL on day t. A trial and error procedure was utilized to examine various input combinations, and their accuracies were evaluated with the Nash-Sutcliffe (E NS ) and root mean square error (RMSE) performance metrics. We evaluated the performance indices from both training and testing phases; however, the final input combination was decided based on the best validation matrices. To predict WL on day t, this combination was identified: where, WL(t) is water depth at the deepest point in the wetland at time t; and Q b , Q s , P and PET denote baseflow (i.e., groundwater contribution to streamflow), stormflow (or quick flow, i.e., surface runoff+interflow), precipitation, and potential evapotranspiration, respectively. These variables were found to serve as the best input combinations and were imported into ANN to predict WLs at their deepest point in the wetland for baseline and future periods. Figure 3 shows a flowchart describing the developed methodology of this study.

El Niño Southern Oscillation (ENSO) Effects on WL Variations
Since rainfall and streamflow are influenced by ENSO in the southeastern US [13], one can also expect ENSO affecting WLs in wetlands. Because long-term observations are not available to study the impacts of climate variability, the developed SWAT-ANN model can be used to generate WLs in the study wetland. To predict WL from 1951 to 2005, observed data from the Fairhope 2 station (USC00012813), which is the nearest climatic station (<10 Km) with long-term historical data, from 1948 to 2005 were used to run the SWAT model. The SWAT model was run with a three-year warm-up period (1948)(1949)(1950) from 1951 to 2005 to prepare the input data required for the trained ANN model. By importing baseflow and stormflow from the SWAT model, as well as precipitation and potential evapotranspiration to the trained ANN, WLs were predicted for this period (i.e., 1951-2005). The WLs predicted by the SWAT-ANN model over the period 1951-2005 were associated with ENSO indices to explore the existence of such a teleconnection at the study wetland. Several correlation methods are available to assess the strength of the relationship between large-scale atmospheric circulation patterns (such as ENSO) and hydro-climatic variables, but the two most commonly applied ones are Pearson's correlation and Spearman's rank correlation. Pearson's correlation needs both variables to be normally distributed while no such assumption is necessary for the application of Spearman's rank correlation [43]. For this study, Spearman's correlation coefficient between the Niño 3.4 index (as an indicator of ENSO strength) and the hydro-climatic variables precipitation and WL were examined. The Niño 3.4 index represents the sea-surface temperature (SST) anomalies in the Niño 3.4 region (5 • N-5 • S, 120-170 • W) and is based on a 3-month running average [13]. The Niño 3.4 indices for the period 1951 to 2005 (to capture approximately 10 ENSO cycles) were obtained from the National Weather Service Climate Prediction Center.
When the Niño 3.4 index is between −0.5 • C and 0.5 • C, it is considered to be a Neutral phase, and indices above 0.5 • C and below −0.5 • C are listed as El Niño and La Niña phases of ENSO, respectively.   Figure 3. Model development flowchart. Q b , Q s , P and PET denote baseflow, stormflow, precipitation, and potential evapotranspiration, respectively. Data were randomly selected in which Levene's test and t-test (Rezaeianzadeh et al. [3]) were used to obtain the optimal datasets for training and testing phases. Figure 4 displays the simulated and observed daily average WLs for the training and testing phases in scatter plots and exceedance curves (i.e., probability of exceedance of "depth below ground surface") formats. The E NS and RMSE values were equal to 0.73, 14.5 cm, and 0.52, 17.8 cm, respectively, for the training and testing phases. Performance indices showed better results in the training phase compared to the testing phase, as expected. Better results regarding performance metrics were achieved too, but taking generalization and overfitting problems into consideration, the reported performances and the resulting model was the best and the most reliable. Figure 4 displays the simulated and observed daily average WLs for the training and testing phases in scatter plots and exceedance curves (i.e., probability of exceedance of "depth below ground surface") formats. The ENS and RMSE values were equal to 0.73, 14.5 cm, and 0.52, 17.8 cm, respectively, for the training and testing phases. Performance indices showed better results in the training phase compared to the testing phase, as expected. Better results regarding performance metrics were achieved too, but taking generalization and overfitting problems into consideration, the reported performances and the resulting model was the best and the most reliable.

SWAT-ANN Model Performance
Although model performance metrics are not very high, as can be seen in Figure 4 the model has mostly captured the observed WL fluctuations well. Despite the model's poor estimations of some extremely high WLs, the testing phase of Figure 4 also confirms the ability of the model to predict WLs accurately. The WL exceedance curve in Figure 4 regarding the training phase illustrates that in general, the model slightly overestimated WLs below −70 cm. On the contrary, it slightly underestimated WLs above −60 cm. Despite the results of the model in the training phase, the model in the testing phase overestimated the extreme high WLs. Similar to the training phase, WLs below −70 cm were overestimated in the testing phase, too. To utilize the developed model for long-term WL predictions, we elaborate on the generalization ability of the model. In making long-term predictions, some extreme values may be observed beyond the range of available time series used for model training, which would require paying attention to the generalization ability of the ANNs. The generalization ability including both interpolations and extrapolations is defined as the model's skill to perform well on a dataset that is not utilized for its calibration [44]. Although it is generally accepted that the ANNs cannot extrapolate beyond the range of the training data [38], utilization of sigmoidal-type, and linear transfer functions in hidden and Although model performance metrics are not very high, as can be seen in Figure 4 the model has mostly captured the observed WL fluctuations well. Despite the model's poor estimations of some extremely high WLs, the testing phase of Figure 4 also confirms the ability of the model to predict WLs accurately. The WL exceedance curve in Figure 4 regarding the training phase illustrates that in general, the model slightly overestimated WLs below −70 cm. On the contrary, it slightly underestimated WLs above −60 cm. Despite the results of the model in the training phase, the model in the testing phase overestimated the extreme high WLs. Similar to the training phase, WLs below −70 cm were overestimated in the testing phase, too.
To utilize the developed model for long-term WL predictions, we elaborate on the generalization ability of the model. In making long-term predictions, some extreme values may be observed beyond the range of available time series used for model training, which would require paying attention to the generalization ability of the ANNs. The generalization ability including both interpolations and extrapolations is defined as the model's skill to perform well on a dataset that is not utilized for its calibration [44]. Although it is generally accepted that the ANNs cannot extrapolate beyond the range of the training data [38], utilization of sigmoidal-type, and linear transfer functions in hidden and output layers, respectively, has been recommended for extrapolation purposes [39,40]. For instance, Rezaeian-Zadeh et al. [42], successfully tested the generalization ability of ANNs to predict hourly air temperatures by using the same type of transfer functions. The other successful example of the extrapolation ability of ANNs (compared to multi non-linear regression) was reported by Cigizoglu [41], on daily river flow data. Hence, in the current study, logistic sigmoid and linear transfer functions were utilized in hidden and output layers, respectively.
Generalization ability is also defined as a function of the ratio of the number of training samples to the number of connection weights. If this ratio is too small, continued training can result in overfitting of the training data [36]. Overfitting (i.e., having high variance) can happen when there are too many features (i.e., inputs) but insufficient amounts of observed data; in which case the learned hypothesis may fit the training set very well, but fail to generalize to new datasets. For this study, we obtained the best architecture to be (8,6,1), where those values refer to the number of inputs, and the number of neurons in hidden and output layers, respectively. Hence, we had 8 × 195 (=1560; 195 corresponds to 70% of daily WL data used in the training phase) training patterns (samples), 54 weights and (6 + 1 = 7) "bias" weights, which constitute a total of 61 connection weights (i.e., free model parameters). Hence, the ratio for the developed model in this study was higher than 25, which confirms a reliable ANN with regard to the generalization ability and overfitting problem. To ensure a good generalization ability, various ratios have been suggested by researchers such as 2 to 1 [45], or 10 to 1 (e.g., [46]). Amari et al. [47], suggested that overfitting does not occur if the above ratio exceeds 30.

Teleconnection between ENSO and WLs
We had only one year of data to train and test the ANN model. Although model performance during both phases are good, the shortness of the data period used in training the model is still a limitation of the model and potentially raises a question on the reliability of extrapolation. Therefore, the results presented below should be interpreted with caution. Note that there are many successful ANN applications in hydrology in the literature with similarly short data (e.g., [48,49]). This fact was acknowledged by the ASCE Task Committee on Application of the ANNs in Hydrology (ASCE 2000) as "Very often, we may have no alternative but to proceed with limited data." Figure 5 displays the Spearman's correlation coefficient between the Niño 3.4 index (as an indicator of ENSO strength) and monthly precipitation and WLs from 1951 to 2005. In August, September, and October, monthly precipitations are negatively correlated with the Niño 3.4 index in all three phases, but not at a statistically significant level (α = 5%). From November to July, the correlations between monthly precipitations and the Niño 3.4 index are positive for almost all months and phases. However, only January in the neutral phase (p-value = 0.036) and February in the La Niña phase (p-value = 0.030) had statistically significant correlations (approximately r = 0.28 for both cases). Average monthly WLs in March and April have a statistically significant negative correlation with the Niño index in both El Niño and La Niña phases, with a stronger correlation with the former. WLs in May also showed significant negative correlation with the Niño index in an El Niño phase. There are no statistically significant correlations between monthly WLs and the Niño index during the neutral phase.  Figure 6 displays the Spearman correlation coefficients between the Niño index and precipitation/WL at the seasonal level. Correlations between the Niño index representing various ENSO phases and seasonal precipitations are non-significant except for a positive correlation in winter during the El Niño phase. This also aligns with what others reported (e.g., [50,51]). Correlations between the Niño index and seasonal WLs are non-significant except for spring during the El Niño phase, which has a negative correlation. Hence, lower WLs can be expected during the El Niño phase in spring. In late spring, Alabama is anomalously dry during the resurgence phase of El Niño. In the late spring of a developing El Niño, low-level winds that bring moist air from the Gulf of Mexico to the U.S. shift westward [52], which again confirms having less precipitation (than normal). Mearns et al. [53], found that precipitation over the southeastern U.S. shows different seasonality over the region.  Figure 6 displays the Spearman correlation coefficients between the Niño index and precipitation/WL at the seasonal level. Correlations between the Niño index representing various ENSO phases and seasonal precipitations are non-significant except for a positive correlation in winter during the El Niño phase. This also aligns with what others reported (e.g., [50,51]). Correlations between the Niño index and seasonal WLs are non-significant except for spring during the El Niño phase, which has a negative correlation. Hence, lower WLs can be expected during the El Niño phase in spring. In late spring, Alabama is anomalously dry during the resurgence phase of El Niño. In the late spring of a developing El Niño, low-level winds that bring moist air from the Gulf of Mexico to the U.S. shift westward [52], which again confirms having less precipitation (than normal). Mearns et al. [53], found that precipitation over the southeastern U.S. shows different seasonality over the region.

Discussion
Headwater wetlands of coastal Alabama have two distinct aspects: (i) they respond quickly to rainfall events, which represent the stormflow component, and (ii) they are groundwater driven in which this component apparently plays a crucial rule in these ecosystems. The idea behind the coupled hybrid model was to capture the hydrologic inputs from the contributing watershed to the study wetland. In this way and by utilizing the SWAT model, which partitions the simulated flows into stormflow and baseflow, we gained better insight into the hydrology in the study wetland. To elaborate more on this, spring precipitation and WL can be a typical example of the situation. There is a (low) positive correlation between spring precipitation and the Niño index in an El Niño phase and while one can expect to see an increase (although insignificant) in the spring WLs, this is not the case here. Spring WLs have been probably affected more by the lagged response of subsurface hydrology and groundwater recharge to decreased precipitation at the end of summer during the El Niño phase and early autumn than spring precipitation. The developed hybrid model in this study requires neither WL data from nearby wetlands nor antecedent WL data from the target wetland itself as input, which makes it practical for long term predictions. To the best of the authors' knowledge, the developed methodology in this study is the first of its own in wetland hydrology studies, which feeds the simulated surface and groundwater components from a semi-distributed watershed hydrology model (SWAT) with a neural network model.

Discussion
Headwater wetlands of coastal Alabama have two distinct aspects: (i) they respond quickly to rainfall events, which represent the stormflow component, and (ii) they are groundwater driven in which this component apparently plays a crucial rule in these ecosystems. The idea behind the coupled hybrid model was to capture the hydrologic inputs from the contributing watershed to the study wetland. In this way and by utilizing the SWAT model, which partitions the simulated flows into stormflow and baseflow, we gained better insight into the hydrology in the study wetland. To elaborate more on this, spring precipitation and WL can be a typical example of the situation. There is a (low) positive correlation between spring precipitation and the Niño index in an El Niño phase and while one can expect to see an increase (although insignificant) in the spring WLs, this is not the case here. Spring WLs have been probably affected more by the lagged response of subsurface hydrology and groundwater recharge to decreased precipitation at the end of summer during the El Niño phase and early autumn than spring precipitation. The developed hybrid model in this study requires neither WL data from nearby wetlands nor antecedent WL data from the target wetland itself as input, which makes it practical for long term predictions. To the best of the authors' knowledge, the developed methodology in this study is the first of its own in wetland hydrology studies, which feeds the simulated surface and groundwater components from a semi-distributed watershed hydrology model (SWAT) with a neural network model. Slight changes in wetland hydrologic conditions can affect soil biogeochemistry and consequently nutrient cycles; biota may respond with considerable changes in species composition and richness and in ecosystem productivity [4]. In forested wetlands, hydrologic changes can also have substantial effects on forest productivity and carbon cycling [18], directly through duration, intensity, and frequency of flooding, and through changes in vegetation [54]. Barksdale et al. [18], found a strong correlation between WL and the total carbon content at the soil surface using data from 15 headwater wetlands (which included the wetland studied in this paper). According to their findings, wetlands with lower median groundwater levels showed surface soils with less organic matter and reduced carbon storage. WLs are expected to drop in spring during El Niño in the study wetland and this will potentially lead to reduced organic material and carbon stock and the same impacts could be expected in the other headwater wetlands of Baldwin County, Alabama. Reduction in WLs could also result in changes in the vegetation communities since the depth and percentage of time a location remains inundated has been shown to have a direct influence on the vegetative communities [55]. This could cause the establishment of invasive plant species, which are major threats to local diversity and other ecosystem functions [20]. With lowered WLs and the consequent increase in oxygen availability, a reduction in CH 4 emission and an increase in CO 2 production will be the most probable scenarios expected. With respect to the CH 4 emission, Zhu et al. [56], concluded that an El Niño event can trigger water table reduction and increased CH 4 soil consumption, leading to a reduction in CH 4 emission from tropical wetlands. Hydrologic alteration of headwater wetlands has the potential to impact some wildlife species, but the most severe impacts would be to amphibians since they are highly vulnerable to wetland drainage [2].
In this study, the potential impacts of climate variability on wetlands through alteration in one of the hydrologic indices (i.e., wetland water level) were discussed only. Other climate-related variables could have noticeable impacts such as increased temperature and altered evapotranspiration, altered biogeochemistry, altered amount and patterns of suspended sediment loadings and fire [57]. Also, human-induced alterations (which were beyond the scope of this study) including wetland drainage, filling/dredging, water diversions, introducing pollutants, and unsustainable developments need to be taken into account in the evaluation of potential impacts on wetlands and to avoid the detrimental impacts of these significant stressors.

Summary and Conclusions
Predicting WLs in wetlands where the hydrology is controlled by both surface runoff and groundwater is a challenging task. Acquiring knowledge about the probable changes in long-term WL fluctuations in wetlands and identifying their trends are even more complicated because of the complexity of hydrological processes in these ecosystems. In this study, we built upon the study of Rezaeianzadeh et al. [3], where the need for WL data from nearby wetlands as inputs to the model has been eliminated. Using simulated baseflow and stormflow components from the SWAT model run by default model parameters as hydrologic inputs to an ANN, this semi-distributed watershed hydrology model was externally coupled with an ANN. Although coupling data-driven models with semi-distributed hydrologic models (such as SWAT) are not new, to the best of the authors' knowledge, the developed model for the prediction of daily WLs is a novel approach in wetland hydrology studies. The ANN component of the hybrid model requires only the WL measurements over a certain period of time for training of the ANN.
The utility of the proposed model was demonstrated in a headwater wetland in coastal Alabama. The model performance statistics E NS and RMSE were 0.73 and 14.5 cm, respectively, for the training phase, and 0.52 and 17.8 cm, respectively, for the testing phase. We had only one year of data to train and test the ANN model, and with this limited amount of data the model had good skill in predicting WLs. We acknowledge that although the availability of a longer wetland WL data series is desirable to better demonstrate the reliability of the developed model, the length of the WL data record was limited in this study and thus was the only option to gain some insights into the variations in the WLs in the study wetland. To illustrate the potential use of wetland WL prediction, the developed hybrid model was used to generate WLs for the period 1951-2005 to explore the potential teleconnections between ENSO and WLs in the study wetland.
In summary, the following conclusions are drawn from this study: • Coupled SWAT-ANN model was proved to be a viable tool to simulate WLs for wetlands at a daily time scale. By examining the ratio of the number of training samples to the number of connection weights, the model development was conducted cautiously to avoid the overfitting problem that is commonly ignored in the application of machine learning techniques in the hydrology and water resources field. Note that there are more sophisticated techniques to avoid overfitting, such as adding regularization to the cost function, and dropout regularization [58], but those have not been considered in this study.

•
Correlations between the Niño index representing various ENSO phases and seasonal precipitations were non-significant at the study wetland, except for a positive correlation in winter during the El Niño phase. Correlations between the Niño index and seasonal WLs were non-significant, except for spring during the El Niño phase, which had a negative correlation. Hence, the findings suggest that winter gets wetter with regards to precipitation and spring gets drier in terms of WL over the El Niño phase in the study area.

•
The teleconnection between WLs and ENSO phases shown in this study can also have important implications for the wetland vegetation and the functioning of wetlands. Specifically, in our study wetland, a reduction in WLs are expected in spring during El Niño and this would potentially lead to reduced organic material and carbon stock and the same impacts could be expected in the other headwater wetlands of Baldwin County, Alabama.
The methodology proposed in this paper should be transferable to other wetlands and their draining watersheds within different geographic and climatic regions to evaluate probable WL fluctuations in response to climatic variations. Studies such as the one presented in this paper are important in order to understand the potential hydrological changes associated with climate variability. These findings are of interest to wetland scientists, eco-hydrologists, wetlands managers, as well as domestic/international policy makers.