Next Article in Journal
Response of Landscape and Ecological Characteristics to the Optimal Rainwater Harvesting Dual-Element Mulch Covered Soil Model in Beijing
Next Article in Special Issue
Behavior of TiO2 and CeO2 Nanoparticles and Polystyrene Nanoplastics in Bottled Mineral, Drinking and Lake Geneva Waters. Impact of Water Hardness and Natural Organic Matter on Nanoparticle Surface Properties and Aggregation
Previous Article in Journal
Climate Change Impacts on Drought-Flood Abrupt Alternation and Water Quality in the Hetao Area, China
Previous Article in Special Issue
An Improved Method Constructing 3D River Channel for Flood Modeling
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Dominant Influencing Factors of Groundwater Recharge Spatial Patterns in Ergene River Catchment, Turkey

1
Civil Engineering department, Hydraulics division at Yildiz Technical University, Istanbul 34349, Turkey
2
Civil Engineering department, Yildiz Technical University, Istanbul 34349, Turkey
*
Author to whom correspondence should be addressed.
Water 2019, 11(4), 653; https://doi.org/10.3390/w11040653
Submission received: 12 February 2019 / Revised: 15 March 2019 / Accepted: 26 March 2019 / Published: 29 March 2019
(This article belongs to the Special Issue Advances in Hydrogeology: Trend, Model, Methodology and Concepts)

Abstract

:
Groundwater is of great significance in sustaining life on planet earth. The reliable estimation of groundwater recharge is the key understanding the groundwater reservoir and forecasting its potential accessibility. The main objective of this study was to assess the groundwater recharge and its controlling factors at the Ergene river catchment. A grid-based water balance model was adopted to determine the spatially distributed long-term groundwater recharge and other water budget components, relying upon the hydro-climatic variables, land-use, soil, geology, and relief of the investigated area. The model calculations were performed for the hydrological reference horizon of 20 years at a spatial resolution of 100 × 100 m. The base flow index (BFI) separation concept was applied to split up the simulated total runoff into groundwater recharge and direct runoff. Subsequently, the statistical methods of Pearson product–moment correlation and principal component analysis (PCA) were combined for identifying the dominant catchment and meteorological factors influencing the recharge. The average groundwater recharge over the investigated area amounts to 95 mm/year. The model validation and statistical analysis indicate that the difference between simulated and observed total runoff and recharge values is generally under 20% and no significant inconsistency was observed. PCA indicated that recharge is controlled, in order of significance, by land-use, soil, and climate variables. The findings of this research highlight the key role of spatial variables in recharge determination. In addition, the generated outputs may contribute to groundwater resource management in the Ergene river catchment.

1. Introduction

Groundwater is the capital source of freshwater for nearly half of the earth’s population for irrigation and domestic water needs [1]. It is identified as a renewable water resource for supporting agricultural, industrial, environmental, and municipal domestic water demands [2]. The estimation of groundwater recharge is the key to understanding the groundwater reservoir and forecasting its potential accessibility and sustainability, even though other elements have to be taken into account; e.g., social, economic, and hydrogeological considerations [3]. Nevertheless, groundwater recharge is the least understood component, mainly due to its extensive spatial and temporal variation and that it requires much effort for direct measurement [4]. The studies about groundwater recharge in a catchment or on a regional scale are crucial for determining the quantitative status of groundwater resources.
Existing studies about groundwater in the Ergene catchment focused only on the groundwater quality [5,6], not on the quantity. The spatially distributed water budgets, including the recharge, have not been broadly ascertained and properly understood. Despite this, the groundwater is the central and vital source of municipal domestic water for 1,150,000 people and industrial water demand for 2500 leather and textile businesses; furthermore, agriculture covers 73% of the total land of the Ergene catchment [7]. Sustainable water resource management, where the influence of agricultural activity over a catchment area is high, aims to guarantee a durable, convenient, and long-standing water supply system to achieve the plants’ water needs, along with an increase in a number of industries and municipal drinking water needs, through mitigating harmful environmental effects. Fundamental water resource management is directly linked with the estimation of accessibility of water within a drainage area and implies the understanding of catchment groundwater hydrologic processes. In this regard, the spatial-temporal pattern of the groundwater recharge generated in this study will have a significant role in supporting the water resource management scheme of the Ergene catchment. Groundwater recharge refers to deep percolation of water from the land surface through the vadose zone and replenishing of the groundwater aquifer. Its quantitative assessment plays a central role in the awareness of large-scale hydrological processes and groundwater supply management and sustainability. To expand the approach of water resources management in the Ergene River catchment, the water budget components, e.g., groundwater recharge and direct runoff, have to be separated and calculated. Estimation of groundwater recharge rates are indispensable for sustainable groundwater management and proper planning in the Ergene catchment.
Various methods may be used to estimate the recharge [8], including direct measurement techniques, e.g., seepage meters and tracer methods. However, the field measurement techniques seem to be extremely hard [9], and are influenced by assessment errors, spatial irregularity patterns, and are operated at high cost [10]. The recorded baseflow from the river gauging stations and its analysis can be used to estimate the level of groundwater recharge [11,12] as a comparison method. Base flow is taken into account as a groundwater recession flow, which is the portion of groundwater flow to streams during rainless seasons. The main consideration in this approach, in order to determine the recharge levels, is that base flow is taken as equivalent to the total groundwater recharge of a basin and the groundwater discharge is roughly considered as the recharge [13]. The base flow analysis method does not consider the origin of the water; however, it considers that great rates of base flow are generated by flow coming from the groundwater reservoir. Thus, for a long period of time, base flow rates were generally taken as watershed groundwater recharge [14]. This method is broadly useful to validate groundwater models which are based on empirical water balance techniques [15]. However, this approach involves the complexity of directly comparing the recharge and base flow, while the latter approach estimates the substitution of groundwater discharge and thus of real recharge [13]. Computer-based models, which rely on a soil–water balance approach, may also be applied for determining the aquifer recharge when other quantities of water balance are well-known with high accuracy [16]. However, the outputs of the model are limited in the absence of model calibration and model validation due to the considerable uncertainty in input parameters. In view of this, for simulating spatial and temporal long-term groundwater recharge, empirical models have proved their authenticity in macroscale catchments [11,17]. Geohydrological models have also been used in regional estimation of groundwater recharge [18], where their applications include zoning via meteorological landscape and soil attributes. Recharge has been considered to be deeply spatially and temporally variable and not a linear function of climate, soil group types, land-use types, and other parameters of a drainage area [19]. Against this background, techniques have been developed in order to combine the spatial and temporal patterns of recharge [20,21]. Therefore, the climate and land-use effects on water budget quantities [22,23] can be spatially and temporally assessed [17] at the catchment or regional scale.
For the determination of groundwater recharge based on an empirical water balance model, a large variety of climate and catchment data (land cover, soils groups, geologic data, and topography) is required. The variety of the input data, e.g., land use, soils, and topography with various spatial and temporal ranges, together with various hydrological data requirements, complicates the interpretation of the outputs in terms of dominating factors. To generate common tendencies, the complexity and the dimensional space in data may be reduced in order to make interpretation easier [24]. Against this background, principal component analysis (PCA) can be applied, defined as a multivariate statistical method which applies orthogonal alteration to turn a set of the most correlated parameters into a set of linearly noncorrelated parameters, known as principal components or eigenvectors, while keeping the maximum amount of total variance in the dataset [24]. Every component is taken from a set of interrelated variables, which are affected by identical process, e.g., recharge and catchment data.
The main objective of this study is to quantify the groundwater recharge and to assess the recharge-controlling factors. In this study, a newly developed grid-based water balance model, mGROWA [11], was applied, due to its accuracy in simulating the spatially distributed water budget components, including groundwater recharge at a large scale, to determine the actual evapotranspiration, total runoff, direct runoff, and recharge in the Ergene catchment. The results of the model were interpreted by developing a methodology for assessing the most dominant factors influencing the groundwater recharge. For this purpose, principal component analysis (PCA) was applied to filter the data; accordingly, important independent basin factors presenting high variance in groundwater recharge were detected. The study was split into five sequential tasks: (i) to create a regular and reliable GIS input database made up of climate data (e.g., precipitation), land cover data, type of the soil, and topographical and geological data at the catchment scale; (ii) to perform, at a macroscale catchment, water budget analysis with the purpose of determining the renewable water resources, e.g., total runoff; (iii) to quantify predominating runoff components (direct runoff and groundwater recharge) by applying the base flow index (BFI) technique; (iv) to calibrate and validate the mGROWA model outputs by referring to a study area database of recorded runoff from different gauging stations; and (v) finally, to apply statistical methods to determine the dominant catchment and climate factors influencing groundwater recharge. To our knowledge, this is the first research to determine the temporally and spatially distributed groundwater recharge for sustainable groundwater reservoir management in the Ergene catchment.

2. Materials and Methods

2.1. Study Area

The Ergene catchment is situated in the middle of the Thrace zone, at the north-west part of Turkey. The catchment is bordered by the north Marmara catchment, Meriç catchment, and Bulgaria, and is located approximately between 40° 48’ 0” and 42° 09’ 0” N (latitude) and 26° 19’ 0” and 28° 11’ 0” E (longitude) (Figure 1). It is spread over 1.4% of Turkey and covers an area of about 11,021.5 km2. Its topography reveals that the minimum and maximum elevations are 2.3 m and 1020.6 m above sea level, respectively. The main river of the catchment flows from the West part in the Istranca Mountains to the East part of the basin with 265 km of the length and the mean annual flow rate of 28.7 m3/s. The Ergene River and its affluent streams are significant sources of water in the north-western part of Turkey (Figure 1). The Ergene river catchment is dominated by terrestrial climate, summers are hot and arid in the north, whereas the winters are cold. In the south part, it is characterized by a Mediterranean climate where the summers are hot and dry and winters are warm and rainy. During the year, January and February are the coldest months, with an average minimum temperature of 4 degrees Celsius, whereas June and July are the hottest months, with an average maximum temperature of 20 degrees Celsius. The average annual temperature of the Ergene river basin is around 13 degrees Celsius. The mean annual long-term precipitation is about 590 mm, the highest precipitation amounts occur in November and the minimum in August. In the south, annual precipitation distribution increases from 650 mm/year to 730 mm/year, whereas in the northern part, the precipitation ranges from 600 mm/year to 900 mm/year [25]. In the Ergene river basin, since 1991, the precipitation information has been collected and no remarkable trend was found in precipitation [26,27]. The annual average evapotranspiration ranges from 450 to 950 mm. Hydrogeological units in the Ergene river basin area mainly composed by metamorphites (gneiss, schist, and marble) and clastic rocks [5,25]. The southeastern part is dominated by basalt fields, whereas alluvium exists in the valley floors, and tectonic lines cover the northern and northeast parts of the investigated area [5,25].
Figure 1 shows the Ergene river catchment geographic context, subcatchments, the location of available river gauging stations and monitoring wells. In Figure 1, the investigated sub-catchments, including Babaeski (488.08 km2), Inanlı (0.99 km2), Lulerburgaz (2697.16 km2), soğucak (71.45 km2), Çayırderekoy (54.48 km2), Ayvacık (1.3 km2), Payralı (90.38 km2), K.yoncak (3.37 km2), Inecık (100 km2), and Üzunköprü (10,026 km2), are presented. These subcatchments were delineated according to the available stream gauging stations with the ArcGIS watershed tool. The 32 monitoring wells existing in the investigated area are also presented in Figure 1. In this study, the groundwater depth data were collected from these wells for the period 1991–2010, where at least one record per month was available for each well.

2.2. Water Balance Modelling Based on the mGROWA Model

mGROWA is a grid-based hydrological water balance model developed for high resolution, long-term simulation of direct runoff, evapotranspiration, and groundwater recharge, and is applicable in macroscale catchments and on a regional scale [11]. mGROWA, in the first run, physically determines the total runoff and actual evapotranspiration in daily steps and the dynamics of soil water in the root zone (Equation (1)).
Δ S = p + q c r e t a q t
where S is the storage or soil moisture content in the root zone (mm) of areas covered with vegetation or the water stored on impervious surfaces in the case of urban areas. p stands for the precipitation in a grid cell (mm/day), qcr is the capillary rise from the groundwater, eta is the actual evapotranspiration (mm/day) determined by using Equation (2), and qt is the total runoff (mm/day).
e t a = e t 0 K L N f ( β , γ ) f ( s )
where et0 is the grass reference evapotranspiration (mm/T), which, in this study, was estimated using a modified Hargreaves and Samani approach [28], KLN stands for crop coefficiency and corrects et0 respective of land-use. From this perspective, it considers various species of vegetation with properties different from the reference grass area. f(β,γ) is a terrain characteristic used to modify eta according to the slope γ and exposition β, and f(s) stands for a storage function.
f ( β , γ ) = γ [ 1.605 10 2 sin ( β 90 ) 2.5 10 4 ] + 1
f(β,γ) was obtained by applying Equation (3) [29].
Figure 2 presents the aspect and slope of the Ergene river catchment. It can be seen that the minimum slope in the study area is 0%, whereas the maximum slope is 66.5% in the northern part of the catchment. The arithmetic mean slope was observed to be 5% with a standard deviation of 4.3%. The aspect reflects the structure of the land surface in the investigated area. The spatial model’s resolution is arbitrary and can be freely chosen. For the Ergene catchment, a 100-m grid is considered suitable, resulting in a total number of 1,036,861 grid cells. The water balance quantities were calculated individually for each grid cell and only the vertical motion is considered as the grid cells are not interconnected; therefore, the lateral flow is not counted in mGROWA.
In the second step, mGROWA empirically divides the total runoff into direct runoff and groundwater recharge based on monthly time steps as the groundwater management does not require a higher temporal resolution [11]. The separation of total runoff into dominant components (direct runoff and groundwater recharge) is based on the theory of base flow indices (BFI) as described in Kunkel and Wendland, (2002) (Equation (4)) [16].
q t = B F I q t + ( 1 B F I ) q t = q r + q d
where qd represents direct runoff and qr is the groundwater recharge.
The attribution of area-differentiated baseflow indices by mGROWA is based on BFI from literature depending on the type of rocks (unconsolidated and hard rock types).
Table 1 shows the baseflow indices values used in mGROWA for unconsolidated rock areas. As can be seen in Table 1, the most important attributes for estimating baseflow indices are groundwater table, water logging tendency of the soil, the degree of sealing, and the area-related slope gradient. In the solid rock areas, the baseflow index values presented in Table 1, may not be applied due to the influence of geological subsoil conditions that predominate the soil properties in quantifying the runoff fractions [30].
Table 2 presents the categories of geological characteristics of hard rocks, such as the permeability, as stated by Bogena et al. (2003) [34], and aquifer typologies initiated by Wendland et al. (2007) [35]. The baseflow indices approach is established by assuming that aquifers possessing similar petrographic characteristics exhibit similar or comparable hydrodynamic and hydrologic conditions regardless of their regional occurrence worldwide. Then, baseflow indices are given for single model grid cells of a groundwater recharge model depending the geological information. This practice has had successful application in Germany [11], Turkey [36], Slovenia [37], Greece [38], and Cameroon [39], yielding realistic groundwater recharge values.
Figure 3 shows the base flow index (BFI) methodology applied in mGROWA for simulating the groundwater recharge in the Ergene catchment. In the determination of groundwater recharge, the hierarchical method is split into five tasks. Firstly, in every grid cell, it is supposed that recharge is negligible for an impervious surface (1). Subsequently, it checks the remaining grid cells for the existence of significant artificial drainage, such as channel or ditch drainage, and its corresponding base flow index is determined (2). On the condition that there is no impervious surface or artificial drainage in the investigated area, the analysis of the type of the rocks (consolidated or unconsolidated rock) is done (3). For consolidated rock, attributes related to the hydrogeology of the case study are applied in BFI estimation (4). For unconsolidated rock regions, the BFI is estimated by taking into account the groundwater table levels, the slope, and water logging tendency (5). Accordingly, the impervious surface, e.g., urban areas, display a total runoff which is equal to the direct runoff. Under that surface, the recharge is equivalent to zero as on the impervious layer; no penetration of the water occurs. In the pervious surface, e.g., vegetation cover, the groundwater table level plays an important role in separation. In a shallow zone, where the influence of capillary rise is observed, the net groundwater recharge will be determined. Where the influence of capillary rise is not considered, the rock types control the recharge. Depending on the types of rocks, the BFI can be calculated [16].

2.3. Baseflow Determination

There are only ten streamflow observation stations available with a long hydrograph record for the period of 1970–2014. The related gauging stations and subcatchments were delineated by means of the hydrology tool in Arc-GIS on the basis of the digital elevation model and river grid map of the Ergene catchment (Figure 1). For the validation assessment of mGROWA outputs, the daily mean of total runoff at selected gauging stations was evaluated. The ten streamflow stations constitute different areas within the basin, from 1 to 10,026 km2, and for a great variety of climatological, geological, and site soil conditions. The observed mean long-term total runoffs were applied to check the validity of the simulated mean long-term total runoff. Groundwater recharge can only indirectly be determined from the recorded streamflow by applying relevant separation methods. The base-flow for selected subcatchments was separated from streamflow based on monthly low water flow rates (MoLR) [40]. Wundt (1958) [40] demonstrated that in unconsolidated aquifers, for prolonged series of years (20 years), the mean MoLR values are good estimates of groundwater recharge. This method is referred to as the MoMLR method, where MoMLR stands for Monthly Mean Low water discharge for prolonged series of years. In hard rock aquifer areas, however, the MoLR-values do not fit to the average groundwater recharge, as MoLR includes considerable components of direct runoff; e.g., direct runoff and interflow, considered also by the MoMLR method. Therefore, Kille (1970) [41] and Demuth (1993) [42] developed a suitable hydrograph separation method applicable to consolidated aquifers, to permit a reduction of MoMLR values by the interflow quantity.
After determining the regional averages of the base flow index values for the investigated basin areas based on river discharge, it was followed by the disaggregation and the transfer base flow indexes to areas where there is no available suitable gauge data. Thus, there was a consideration of parameters and runoff-effective regional features where the hierarchical approach in Figure 3 was taken into account, in which the value of a site condition is taken into account to be exclusively estimating the base flow index value. In the procedure of iterative calibration [36], BFI was calculated by adopting Kille’s (1970) techniques [41], and various characteristics (degree of sealing, groundwater depth, water logging tendency, slope, and permeability) were identified. Among the total gauging stations with possessing data, half of the drainage points were used for calibration and the other half were used for validation of groundwater recharge [43]. Accordingly, a set of base flow index values is established and guides the most favorable fit for the investigated area.

2.4. Statistical Analysis

The comparisons of simulated and recorded water balance quantities in the Ergene catchment were assessed using the Nash–Sutcliffe index [44] to spatially evaluate divergences in simulated groundwater recharge and total runoff with recorded hydrographs in subcatchments possessing sufficient stream flow observations. The Nash–Sutcliffe index, also referred to as the efficiency index (Ef), has been widely applied and is a valid statistic indicator for evaluating the goodness of fit of a model [45]. Ef may vary in the range of −∞ and 1. The acceptable range of agreement varies from 0 to 1, a negative value reveals disagreement of the performance of the hydrological models, and a value of 1 shows excellent performance. Herrmann et al. (2015) [11] presented a different concept for simultaneously applying Ef to the all considered stream flow stations and their corresponding drainage areas within a catchment by applying Equation (5).
E f = 1 i = 1 n ( A i ( Q o b s e r v e d , i Q s i m u l a t e d ,     i ) 2 ) i = 1 n ( A i ( Q o b s e r v e d , i Q o b s e r v e d , A V ) 2 )
where Qobserved stands for the recorded stream flow value per unit of area, Qsimulated the related modeled stream flow, Q o b s e r v e d , A V the recorded overall flow rate per unit of all sub-basins taken into account. A is the specific basin area and i the basin indicator.
In this study, a multivariate technique called principal component analysis (PCA) was adopted through the XL-STAT software. PCA alters an original set of dependent variables into an equivalent number of noncorrelated factors, referred to as principal components (PCs). In this study, PCA was adopted to reduce the amount of redundant data and to determine the dominant variables that show most of the variance of groundwater recharge, thus the factors most importantly influencing the spatial variation of recharge. Its purpose is to withdraw the important information from the dataset. This approach squeezes the data by maintaining only the most significant information, facilitating the explanation, and investigating the form and constituents of the data [46]. A varimax rotation method was conducted in order to increase the variation between the parameters for every component [47].
In the first step, the strength of existing correlations between groundwater recharge and catchment attributes were determined by the Pearson product–moment correlation coefficient [48]. For the matrix of interrelated variables, the linear relationship among the catchment attribute was tested by measuring the correlation between parameters in the range of [−1,1]. Greater coefficients show the importance of the correlation and high interdependence among the variables. A positive coefficient represents similarity and harmony among the interrelated variables, and a negative sign reflects the opposite [49].
In the second step, the variables that displayed significant correlation with the recharge were analyzed by applying PCA. The first PC reflects a great portion of the variance in the original information. The PCs possessing an eigenvalue greater than one and containing a high factor loading (>0.7) were considered as the variables that control the spatial variation of recharge [12,17]. Following PCA, the highly weighted variable in a principal component was determined based on the level of correlation coefficients and their sums [48]. This approach minimized the number of variables to produce a small data set that highly influenced recharge.

2.5. Input Data

All the climatic data were recorded from 41 stations located in the Ergene catchment and were provided by the Turkish General Directorate of State Hydraulic Works (DSI). The database consists of daily precipitation and daily extreme temperature data. The latter data was used to determine the reference grass evapotranspiration level at corresponding stations by applying Hargreaves’ method [50]. This method was established by only referring to the mean air temperature of a region, temperature extremes, and extraterrestrial radiation. It is approved by the Food and Agriculture Organization of the United Nations (FAO) for its reliability and accuracy in estimating reference evapotranspiration where there are gaps in climatic data; e.g., in relative humidity, net solar radiation, and wind speed [50]. However, in humid regions, the Hargreaves equation requires calibration [51]. Cobaner et al. (2017) [28] proposed a modified Hargreaves equation to improve its estimation throughout Turkey. The dataset was produced as digital datasets for the hydrologic period of 1991 to 2010 in a gridded pattern for the study area.
The data was used to produce catchment-wide 100 × 100 m data by applying geographic information system (GIS) interpolation techniques saved in an American standard code for information interchange (ASCII) text grid file applying the unique value. Precipitation and grass reference evapotranspiration values were interpolated using inverse distance weighted techniques (IDW) [52] and their values were allocated to every single grid (100 × 100 m) of the catchment in ArcGIS. This method was applied to produce the area distribution of the precipitation and grass reference evapotranspiration in the Ergene catchment. The climate input data was developed spatially and temporally, such that long-term mean precipitation and grass reference evapotranspiration for the period of 1991–2010 over the study area was approximated to 591.29 ± 66.9 mm/ year and 1035.26 ± 34.46 mm/year, respectively. The more standard deviations signify that the spatial variability of climatic data is significant over the catchment. Then, the spatially distributed climatic data were used in mGROWA for the entire study catchment. The spatially distributed daily mean climatic data was aggregated into monthly values and the grids of monthly sums of precipitation and grass reference evapotranspiration in mm were prepared as model input. Table 3 shows the prepared dataset of mGROWA for the study area.
The land use information was required in order to estimate the real vegetation evapotranspiration and was collected from the digital landscape model. The concept of CORINE (coordination of information on the environment) land cover contains 44 land cover classes, of which 26 occur in the study area (Figure 4).
Figure 4 presents the land use types of the Ergene river catchment for the year 2012. In Figure 4, one can see that the majority of the Ergene catchment is covered by vegetation, i.e., nonirrigated arable land of code 211 (approx. 54.9% of the total study area), followed by permanently irrigated land of code 212 (approx. 8%), land principally occupied by agriculture of code 243 (approx. 7.1%), transitional woodland-shrub of code 324 (6.8%), broard-leaved forest of code 311 (6%), pastures of code 231 (3.9%), natural grasslands of code 321 (3.2%), rice fields of code 213 (2.6%), and discontinuous urban land (112) makes up 1.7%. The share of other categories exhibits low percentages (below 1%). For Ergene water balance modeling, the information of imperviousness, rooting depth of plants, and land-use-specific actual evapotranspiration factors were allocated to the Ergene digital landscape model correspondingly.
Figure 5 shows the soil groups in the Ergene river catchment. As can been observed from Figure 5, in the Ergene catchment, eight soil groups exist. In the north of the study area, brown forest soil without lime mostly dominates. The central part of the Ergene catchment is covered by brown soils without lime and vertisol groups. The river beds and its surrounds are characterized by the existence of alluvial soils. The southern part of the study area is covered by brown forest soils with and/or without lime. The urban soils are distributed in small quantities over the catchment. The soil groups mainly occurring in the study area are brown soils without lime, brown forest soils without lime, and vertisols, respectively contributing 33.8, 25.2, and 23.6%, followed by alluvial soils (9.2%), brown forest soils (5.8%), urban soils (2.1%), and rendzinas soil (0.02%). The parameterization of Ergene soil profiles has been done by applying pedo-transfer functions. The rooting depth and effective field capacity were derived to present the storing and releasing capacity of the Ergene soils.

3. Results

3.1. Simulation of Water Balance Quantities

The long-term annual mean of spatially distributed water budget components calculated with the mGROWA model over the study area for the period 1991–2010 are presented in this section.
Figure 6 presents the water budget components: the actual evapotranspiration (ETa), runoffs, and recharge. The runoff maps reveal the individual runoff quantities produced in 100 × 100 m grid cells; therefore, the Ergene River flow rates passing through the investigated area are not included. The actual evapotranspiration is predominantly affected by vegetation cover, soil, and climate. It can be noticed that the area-distributed ETa presents considerable variability in land cover across the study area. The lowest ETa (less than 300 mm/year) was observed where the impervious areas dominate, especially in urbanized areas. A high ETa in the range of 400–450 mm/year appears in the highest altitudes of the northern part where the land is covered by forest. In the low altitude land, ETa is also influenced by the land-use types. In the nonirrigated arable land (dominating in the catchment (55%)), the ETa varies in the range of 350–400 mm/year. The arithmetic mean value of the real evapotranspiration is 450 mm/year over the study area.
Regarding the runoffs (Figure 6), because of the high level of imperviousness, the urbanized areas display the highest total runoff values, with direct runoff levels in the ranges of 350–400 mm/year and 200–300 mm/year, respectively. The low permeability of urbanized areas promotes runoff. The total runoff levels steadily decrease in the high-altitude land located in the northern part due to the forest cover which generates high actual evapotranspiration (less than 150 mm/year). Accordingly, the highland areas in the north-eastern part display total runoff values less than 100 mm/year, whereas in the north-western part of study area, values in the interval of 100–150 mm/year prevail. The generated direct runoff shows similar spatial patterns to the total runoff and both are highly influenced by the land-use type. Across the Ergene catchment (approx. 1,091,645 grid cells), the arithmetic averages for total runoff and direct runoff are 200 mm/year and 85 mm/year, respectively.
Regarding simulated long-term groundwater recharge in Figure 6, the Ergene catchment can be mainly characterized into three zones: high, moderate, and low groundwater recharge rates depending on the zone variety parameters, e.g., soil, geology, the ability of the aquifer to store water, hydraulic conductivities. The influence of hydrogeology is obvious, i.e., the alluvial plains in the lowland of the catchment prevail a high groundwater recharge levels greater than 350 mm/year, because of high infiltration rates in such unconsolidated sand aquifers. Besides, the southern part receives high precipitation, thus, some grids located in alluvial soils display the maximum groundwater recharge values. In the moderate zone, groundwater recharge displays mainly values in the interval of 100–200 mm/year. In the northern zone of the catchment, the groundwater recharge rates are generally diminishing, down to the lowest values (less than 50 mm/year). The calculated long-term average groundwater recharge level over the study area is 95 mm/year.
The results for the average monthly groundwater recharge for the period of 1991–2010 are presented as a map in Figure 7. It can be observed that the groundwater recharge rates start to occur in November and end in April. The highest rates were observed in January. In April, the groundwater recharge starts to decrease as the vegetation time begins.
In Figure 8 and Figure 9, the average annual water budget quantities—actual evapotranspiration, direct runoff, and groundwater recharge—are displayed in terms of existing CORINE land cover classes and soils group types in the Ergene catchment. Their corresponding coverage areas in the catchment are also shown.
In Figure 8, the groundwater recharge shows great variation depending on the land-use. The artificial surface, represented by CORINE_2012 codes 111, 112, 121, and 122 (continuous-discontinuous urban fabric and industrial areas) carry a low groundwater recharge value due to the high level of imperviousness characterizing these types of land-use. Consequently, these land-use types generate high runoff. Descriptions of the codes are presented in Figure 4. Artificial land covers 3% of the area of the Ergene catchment. The higher levels of groundwater recharge can be observed in land use types such as mineral extraction sites, dump sites, and constructions sites (represented by Corine codes 131, 132, and 133, respectively) which only comprise 0.3% of the study area. These land-use types (131, 132, and 133) affect the natural hydrogeology quite profoundly, where the permeability may increase up to 100 times [53,54]. For agricultural areas (211, 212, 213, 221, 222, 231, 242, and 243), the actual evapotranspiration dominates other water balance quantities, which is also the case in the forest and seminatural cover areas (311, 312, 313, 321, 323, 324, and 334). These types of land-use cover displayed a moderate average level of recharge. The exception was observed in bunt cover areas (334), where groundwater recharge (240 mm/year) is higher than actual evapotranspiration.
In Figure 9, the groundwater recharge is shown to be dependent on soil group types in the study area, which is dominated by alluvial soils. This type of soil presents a high ability of water infiltration to replenish the groundwater aquifer. A low average groundwater recharge of 50.3 mm/year is observed for the urban soil type. Urban soil is a type of anthropic soils, widely formed by human activities, in urban areas.
Figure 10 illustrates the descriptive statistics of base flow obtained from the Wundt and Kille separation techniques. The highest median was noticed in Soğucak, Çayirderekoy, Ayvacık, and Poyralı subcatchments of 90.4, 77, 73.5, and 49 mm/year, respectively. The mentioned subcatchments are found in the northern part characterized by high topographic elevation and slope (Figure 2). The relief significantly affects the baseflow and recharge. The increase in high relief, e.g., slope angle, amplifies the level of groundwater transmission [55]. Other subcatchments, especially Luleburgaz, Babaeski, and Inanlı, displayed the lowest minimum median in estimated base flow with 33, 37, and 38 mm/year respectively. These subcatchments are located in areas of medium slope and medium elevation. On the other hand, the low value of base flow in these subcatchments may be justified by their location on land covered by agricultural land (73%) and are frequently watered from groundwater. A portion of urban areas existing in these subcatchments may also negatively influence base flow due to the change in the natural system’s surface and subsurface capacity of to drain water. As a result, the interaction between groundwater recharge and rivers is deteriorated [56]. In the Ergene catchment, it can be concluded that catchment characteristics (e.g., topography, slope land-use types) and the interaction between surface and subsurface water manage the recharge-base flow interconnection.

3.2. Model Evaluation

mGROWA cannot be characterized as a runoff concentration model; therefore, its validation does not rely on comparing daily simulated and observed water balance quantities [11]. However, validation of the mGROWA model outputs is done by comparing the long-term observed and simulated runoffs and recharges. Accordingly, Figure 11 shows the scatter plots with a 1:1 line, which presents the comparisons between modeled total runoff and groundwater recharge with their corresponding observations for the hydrological period 1991–2010. It can be observed from the left of Figure 11 that the comparison displays a quantitative Nash–Sutcliffe statistic (Ef) of 0.94 for observed-simulated total runoff, which indicates a very good performance of the hydrological model. On the right side of Figure 11, the good performance of the mGROWA model (Ef = 0.76) for calculated groundwater recharge and base flow obtained from the Wundt–Kille separation techniques was achieved and taken as observed data. As stated by Moriasi et al. (2007) [57], a Nash–Sutcliffe efficiency value greater than 0.75 and 0.65 shows a very good and a good performance of the hydrological model in an individual basin, respectively. As can be observed, the deviations between observed and simulated total runoff and groundwater recharge for nearly all considered subcatchments are below ±20%. For only three of the sub-catchments considered are the deviations between the recorded and modeled groundwater recharge levels higher than 20%. The inaccuracy in simulated values and observations may be caused by scale discrepancy of the empirical model [20]. In this comparison between the observed and the model outputs, there are obvious differences. These inaccuracies can be caused by unreliability in the complexity of the model-data environment, such as the modelling approach, errors occurring in observation, uncertainties occurring in soils groups, land-use types, topography, and climatic data. In addition, the estimated groundwater divides in some considered subcatchments might not be representative for real corresponding river drainage area due to artificial channels, irrigation practices, dams, and pumped groundwater wells that can perturb the natural water budget of a catchment. Besides, the fact of assuming that recharge is equivalent to base flow implies the equivalence between the groundwater shed and the surface watershed, where the groundwater cannot discharge across boundaries. Although groundwater can run into an aquifer, most of the time it overlaps many basins, and within basins, groundwater and surface water boundaries frequently differ from each other, and this may change in time. Therefore, identifying the groundwater shed area that contributes to a river is a hard work. Usually, a dissimilarity between the observed base flow in the river and the estimated recharge may be caused by a lack of consideration of the lateral fluxes in such models that consider only vertical fluxes, and the effect of scale.

3.3. Determination of Dominant Factors

Prior to principal component analysis, the Pearson correlation test was applied to indicate the degree of linear association between catchment characteristics (soils, land use, area, topography, and climate) and the groundwater recharge. Twenty catchment attributes were assembled for processing by means of Arc-GIS tools. Table 4 shows the catchment characteristics, symbols, units, and their correlation with groundwater recharge. The catchment characteristics are mostly presented as a percentage of their corresponding areas relative to the total area of the Ergene catchment. In Table 4, the correlation between the catchments’ characteristics and groundwater recharge are presented. The correlation assessment in Table 5 indicates that among the variables, significant correlations were observed. Consequently, some variables are absolutely redundant. In Table 4, only six catchment characteristics were found to be significantly correlated with groundwater recharge. The climate parameters appeared to influence the area distribution of recharge. The highest positive correlation was obtained by precipitation (0.817), and this was found to be a natural principle in many hydrogeological studies [12].
The actual evapotranspiration showed a negative correlation with recharge. This means that high actual evapotranspiration will deny the aquifer maximum groundwater recharge and low actual evapotranspiration will cause an increase in groundwater recharge. The same finding has been previously presented [58]. It was followed by Brown forest soils without lime, which was detected to be inversely correlated with groundwater recharge. Two types of soil appear to be strong in a relationship with recharge. Alluvial soil presents a strong positive correlation with recharge. The Ergene basin’s alluvial aquifers are formed principally by sandy soil which possesses intermediate to high hydraulic conductivity, thus has high infiltration capacity, leading to an increase of groundwater recharge. Brown forest soils without lime appear to correlate negatively with recharge. The most important constituent of this type of soil is the clay which has low permeability. Two of the five land-use types indicated a significant correlation with recharge. Agriculture land and forests were found to have significant interdependency with recharge. The agglomeration of the forest biomass is likely to increase the evapotranspiration and decrease the total runoff, thus negatively affecting the amount of water that replenishes the aquifer [59]. The area of the catchment and its topography did not show any significant correlation with the groundwater recharge. Principal component analysis (PCA) was applied to the variables most closely correlated with groundwater recharge. AS, BFSWL, AA, F, P, and AE displayed significantly interdependency with the recharge. Table 6 shows the principal component analysis of catchment factors displaying a high correlation with recharge. Eigenvalues greater than one were considered as significant in this study (Table 6). The spatial variance of the recharge for the two considered principal components was 84%. Percentage of forest land cover yielded the largest loading factor, with agricultural land cover coming second, followed by brown forest soils without lime and actual evapotranspiration, respectively, in PC1. The precipitation impacts the recharge more in PC2.
Catchment factors displaying high positive weighting on PC1 include alluvial soil (AS), agriculture area (AA), and precipitation (P) (Table 6). Precipitation showed a strong positive correlation with recharge since precipitation is the main source of groundwater recharge, as confirmed by many studies [17]. Alluvial soils increase the amount of recharge, as this soil type is associated with high rates of sand which allows infiltration into alluvial aquifers [60]. Therefore, it is plausible that alluvial soils are related to an increase in recharge. The positive correlation between recharge and agricultural land use may increase the recharge depending on its patterns of infiltration and evapotranspiration; nevertheless, without knowing the details of these land covers in the Ergene catchment, the direct mechanism cannot be ascertained. Although these interrelations identify the catchment factors which can control the spatial patterns in recharge, further study is recommended to determine the exact processes by which these factors contribute to the increase of recharge across the catchment.
The PCA outputs show that recharge is lower in catchments where brown forest soils without lime, forest, and actual evapotranspiration dominate in Table 6. Brown forest soils without lime were found to have a considerable rate of clay soils, which are highly related with evapotranspiration and runoff, which leads to decreased groundwater recharge. Actual evapotranspiration negatively affects the groundwater recharge by diminishing the available quantity of water for recharge. Forest cover is known to assist in infiltration and possesses a relatively high amount of evapotranspiration, thus low recharge. In accordance with PCA, the recharge variation is controlled, in order of significance, by vegetation land-use, soil group types, and climate. The great influence of vegetation on groundwater recharge was previously confirmed in a global assessment of groundwater recharge [61]. Globally, vegetated land provides the greatest amount of water input, such as precipitation and irrigation water, that becomes recharge [61]. The recharge in coarse soil is much greater than in fine textured soil, thus the movement of the water in the soil is highly controlled by soil texture and structure. This result is in agreement with previous studies [12]. The great influence of climate on recharge was also confirmed by Jan et al. (2007), in their study [62] about rainfall intensity and distribution impact on groundwater level changes, in which they found that the increase or decrease in level of groundwater is mainly a result of the presence or absence of precipitation, respectively. Thus, groundwater directly relies on accumulated precipitation.

4. Conclusions

The Ergene River catchment is one of the important basins in Turkey that are exposed to a high level of groundwater water exploitation due to industrial and agricultural activities. Thus, quantifying the groundwater recharge is important for its management. The spatial groundwater recharge is improperly considered in the modeling of groundwater recharge. However, hydrograph separation methods show spatial recharge distribution. With the development of GIS and the accessibility of spatial inputs, the recharge can be determined as the outcome of water budget modeling. The mGROWA model can, therefore, accurately determine the spatially distributed long-term groundwater recharge of the Ergene River catchment. The study outcomes show that the groundwater recharge at the scale of the Ergene catchment has significant variation, while in small-scale catchments, the variation is less noticeable. Mineral extraction sites, dumps, and construction sites display the highest recharge rates, followed by agriculture land. Regarding soil structure and texture, alluvial soils produce the highest groundwater recharge levels. Urban soil displays the lowest recharge rate. Accuracy is challenging in recharge determination. In this respect, a base flow assessment of ten available river gauging stations was used as a validation approach where the average base flow fluxes for a long period can be considered as a proxy for mean groundwater recharge rates. The assessment of groundwater recharge together with the catchment attributes, such as land-use, soil types, and climate variables, displayed negative and positive correlations. In accordance with PCA, the spatial groundwater recharge variation is controlled in order of significance by vegetation land-use, soil group types, and climate. To aid water resource management in Turkey, the approach developed in this study is suggested to be applied in different catchments across the country. The results may contribute to the establishment of a quantitative status of groundwater water reservoir exploitation and enhancement of the water resource management of the region.

Author Contributions

Conceptualization, E.R. and A.D.; methodology, E.R. and A.D.; software, E.R. and A.D.; validation, E.R. and A.D.; formal analysis, E.R. and A.D.; investigation, E.R. and A.D.; writing—original draft preparation, E.R. and A.D.

Funding

This research was funded through Scientific Research Projects of Yildiz Technical University with project #2015-05-01-KAP04 as well as partially supported by Scientific and Technological Research Council of Turkey (TUBITAK) with project number 115Y008.

Acknowledgments

The authors are grateful to Frank Wendland and Frank Hermann at Forschungszentrum Jülich for providing the mGROWA model and their assistant during the mGROWA model development in the Ergene River catchment. The authors thank also the editor and four anonymous reviewers for constructive and valuable comments.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Shah, T.; Burke, J.; Villholth, K.G.; Angelica, M.; Custodio, E.; Daibes, F.; Hoogesteger, J.; Giordano, M.; Girman, J.; Van Der Gun, J.; et al. Groundwater: A Global Assessment of Scale and Significance; International Water Management Institute (IWMI): Supriyo Das, India, 2007. [Google Scholar]
  2. Yenehun, A.; Walraevens, K.; Batelaan, O. Spatial and temporal variability of groundwater recharge in Geba bain. Northern Ethiopia. J. Afr. Earth Sci. 2017, 134, 198–212. [Google Scholar] [CrossRef]
  3. Seiler, K.P.; Gat, J.R. Groundwater Recharge from Run-Off, Infiltration and Percolation; Water Science and Technology Library: New York, NY, USA, 2007; Volume 55. [Google Scholar]
  4. Simmers, I. Estimation of Natural Groundwater Recharge; D. Reidel Publishing Company: Dordrecht, The Netherlands, 1988; Volume 70, p. 510. [Google Scholar]
  5. Arkoc, O. Heavy metal concentrations of groundwater in the east of Ergene Basin, Turkey. Bull. Environ. Contam. Toxicol. 2014, 93, 429–433. [Google Scholar] [CrossRef]
  6. Füsun, E.; Zehra, B. Linear Alkylbenzene Sulfonates in the Groundwater and Surface Waters: Ergene Basin Case Study. J. Agric. Sci. Technol. A 2017, 7, 369–376. [Google Scholar] [CrossRef]
  7. Arkoç, O. Assessment of scaling properties of groundwater with elevated sulfate concentration: A case study from Ergene Basin, Turkey. Arab. J. Geosci. 2012, 6, 4377–4385. [Google Scholar] [CrossRef]
  8. Scanlon, B.R.; Healy, R.W.; Cook, P.G. Choosing appropriate techniques for quantifying groundwater recharge. Hydrogeol. J. 2002, 10, 18–39. [Google Scholar] [CrossRef]
  9. Sharma, M.L. Measurement and prediction of natural groundwater recharge—An overview. J. Hydrol. 1989, 25, 49–56. [Google Scholar]
  10. Assefa, K.A.; Woodbury, A.D. Transient, spatially varied groundwater recharge modeling. Water Resour. Res. 2013, 49, 4593–4606. [Google Scholar] [CrossRef] [Green Version]
  11. Herrmann, F.; Keller, L.; Kunkel, R.; Vereecken, H.; Wendland, F. Determination of spatially differentiated water balance components including groundwater recharge on the Federal State level—A case study using the mGROWA model in North Rhine-Westphalia (Germany). J. Hydrol. Reg. Stud. 2015, 4, 294–312. [Google Scholar] [CrossRef]
  12. Zomlot, Z.; Verbeiren, B.; Huysmans, M.; Batelaan, O. Spatial distribution of groundwater recharge and base flow: Assessment of controlling factors. J. Hydrol. Reg. Stud. 2015, 4, 349–368. [Google Scholar] [CrossRef] [Green Version]
  13. Risser, D.W.; Gburek, W.J.; Folmar, G.J. Comparison of recharge estimates at a small watershed in east-central Pennsylvania, USA. Hydrogeol. J. 2008, 17, 287–298. [Google Scholar] [CrossRef]
  14. Meyer, S.C. Analysis of base flow trends in urban streams, northeastern illinois, USA. Hydrogeol. J. 2004, 13, 871–885. [Google Scholar] [CrossRef]
  15. Delin, G.N.; Healy, R.W.; Lorenz, D.L.; Nimmo, J.R. Comparison of local- to regional-scale estimates of ground-water recharge in Minnesota, USA. J. Hydrol. 2007, 334, 231–249. [Google Scholar] [CrossRef] [Green Version]
  16. Kunkel, R.; Wendland, F. The GROWA98 model for water balance analysis in large river basins-the river Elbe case study. J. Hydrol. 2002, 259, 152–162. [Google Scholar] [CrossRef]
  17. Rumsey, C.A.; Miller, M.P.; Susong, D.D.; Tillman, F.D.; Anning, D.W. Regional scale estimates of baseflow and factors influencing baseflow in the Upper Colorado River Basin. J. Hydrol. Reg. Stud. 2015, 4, 91–107. [Google Scholar] [CrossRef] [Green Version]
  18. Grinevskii, S.O.; Pozdnyakov, S.P. Principles of regional estimation of infiltration groundwater recharge based on geohydrological models. Water Resour. 2010, 37, 638–652. [Google Scholar] [CrossRef]
  19. Healy, R.W.; Cook, P.G. Using groundwater levels to estimate recharge. Hydrogeol. J. 2002, 10, 91–109. [Google Scholar] [CrossRef]
  20. Bogena, H.; Kunkel, R.; Schöbel, T.; Schrey, H.P.; Wendland, F. Distributed modeling of groundwater recharge at the macroscale. Ecol. Model. 2005, 187, 15–26. [Google Scholar] [CrossRef]
  21. Cooper, D.J.; Wolf, E.C.; Ronayne, M.J.; Roche, J.W. Effects of groundwater pumping on the sustainability of a mountain wetland complex, Yosemite National Park, California. J. Hydrol. Reg. Stud. 2015, 3, 87–105. [Google Scholar] [CrossRef] [Green Version]
  22. Rukundo, E.; Dogan, A. Assessment of Climate and Land Use Change Projections and their Impacts on Flooding. Pol. J. Environ. Stud. 2016, 25, 2541–2551. [Google Scholar] [CrossRef] [Green Version]
  23. Ehlers, L.; Herrmann, F.; Blaschek, M.; Duttmann, R.; Wendland, F. Sensitivity of mGROWA-simulated groundwater recharge to changes in soil and land use parameters in a Mediterranean environment and conclusions in view of ensemble-based climate impact simulations. Sci. Total Environ. 2016, 543, 937–951. [Google Scholar] [CrossRef]
  24. Bücker, A.; Crespo, P.; Frede, H.-G.; Vaché, K.; Cisneros, F.; Breuer, L. Identifying Controls on Water Chemistry of Tropical Cloud Forest Catchments: Combining Descriptive Approaches and Multivariate Analysis. Aquat. Geochem. 2009, 16, 127–149. [Google Scholar] [CrossRef]
  25. Karataş, A. Identifying Surface Runoff Distribution and Amount in Stream Basins: Ergene River Basin. Turk. J. Water Sci. Manag. 2018, 2, 40–81. [Google Scholar] [CrossRef]
  26. Bostan, P.A. Analysis and Modeling of Spatially and Temporally Varying Meteorological Parameter: Precipitation Over Turkey. Ph.D. Thesis, Middle East Technical University, Ankara, Turkey, 2013. [Google Scholar]
  27. Bostan, P.A.; Heuvelink, G.B.M.; Akyurek, S.Z. Comparison of regression and kriging techniques for mapping the average annual precipitation of Turkey. Int. J. Appl. Earth Obs. Geoinf. 2012, 19, 115–126. [Google Scholar] [CrossRef]
  28. Cobaner, M.; Citakoğlu, H.; Haktanir, T.; Kisi, O. Modifying Hargreaves–Samani equation with meteorological variables for estimation of reference evapotranspiration in Turkey. Hydrol. Res. 2017, 48, 480–497. [Google Scholar] [CrossRef]
  29. Golf, W. Ermittlung der Wasserressourcen im Mittelgebirge. Wasserwirtsch. Wassertech. 1981, 30, 90–95. [Google Scholar]
  30. Kunkel, R.; Wendland, F. Der Landschaftswasserhaushalt im Flubeinzugsgebiet der Elbe—Verfahren. Datengrundlagen und Bilanzgröben; Schriften des Forschungszentrum Jülich. Reine Umwelt.: Jülich, Germany, 1998; Volume 12. [Google Scholar]
  31. Dörhöfer, G.; Josopait, V. Eine Methode zur flächendifferenzierten Ermittlung der Grundwasserneubildungsrate. Geol. Clays 1980, 27, 45–65. [Google Scholar] [CrossRef]
  32. Hennings, V. Methodendokumentation Bodenkunde—Auswertemethoden zur Beurteilung der Empfindlichkeit und Belastbarkeit von Böden; Geologisches Jahrbuch, Reihe F: Hannover, Germany, 2000. [Google Scholar]
  33. Wessolek, G.; Facklam, M. Standorteigenschaften und Wasserhaushalt von versiegelten Flächen. Pflanzenernähr. Bodenk. 1997, 160, 41–46. [Google Scholar] [CrossRef]
  34. Bogena, H.; Kunkel, R.; Schöbel, T.; Schrey, H.P.; Wendland, F. Die Grundwasserneubildung in Nordrhein-Westfalen; FZ Jülich GmbH: Jülich, Germany, 2003; Volume 37, p. 148. [Google Scholar]
  35. Wendland, F.; Blum, A.; Coetsiers, M.; Gorova, R.; Griffioen, J.; Grima, J.; Hinsby, K.; Kunkel, R.; Marandi, A.; Melo, T.; et al. European aquifer typology: A practical framework for an overview of major groundwater composition at European scale. Environ. Geol. 2007, 55, 77–85. [Google Scholar] [CrossRef]
  36. Karpuzcu, M.A.; Alpaslan, N.; Engin, G.; Görmann, H.; Gunduz, M.; Kocal, M.; Oncel, S.; Pekdeger, A.; Voigt, H.; Kuhr, P.; et al. Integrated Modelling of Nutrients in Selected Rşver Basins of Turkey; Forschungszentrum Jülich GmbH: Jülich, Germany, 2008; Volume 17. [Google Scholar]
  37. Tetzlaff, B.; Andjelov, M.; Kuhr, P.; Uhan, J.; Wendland, F. Model-based assessment of groundwater recharge in Slovenia. Environ. Earth Sci. 2015, 74, 6177–6192. [Google Scholar] [CrossRef]
  38. Panagopoulos, A.; Arampatzis, G.; Kuhr, P.; Kunkel, R.; Tziritis, E.; Wendland, F. Area-differentiated modeling of water balance in Pinios river basin, central Greece. Glob. Nest 2015, 17, 221–235. [Google Scholar]
  39. Cheo, A.E.; Voigt, H.J.; Wendland, F. Modeling groundwater recharge through rainfall in the Far-North region of Cameroon. Groundw. Sustain. Dev. 2017, 5, 118–130. [Google Scholar] [CrossRef]
  40. Wundt, W. Die Kleinstwasserf¨ uhrung der Fl¨ usse als Maß f¨ ur die verf¨ ugbaren Grundwassermengen; Forsch Deutsch Landeskunde: Jülich, Germany, 1958; Volume 104. [Google Scholar]
  41. Kille, K. Das Verfahren MoMNQ, ein Beitrag zur Berechnung der mittleren langjährigen Grundwasserneubildung mit Hilfe der monatlichen Niedrigwasserabflüsse. Zeitschrift der Deutschen Geologischen Gesellschaft 1970, 122, 89–95. [Google Scholar]
  42. Demuth, S. Untersuchung zum Niedrigwasser in West-Europa. Freiburger Schriften xur Hydrologie; Universitat Freiburg: Freiburg, Germany, 1993. [Google Scholar]
  43. Andjelov, M.; Mikulic, Z.; Tetzlaff, B.; Uhan, J.; Wendland, F. Groundwater recharge in Slovenia: Results of a Bilateral German-Slovenia Research Project; Forschungszentrums Jülich GmbH: Jülich, Germany, 2016; Volume 339. [Google Scholar]
  44. Nash, J.E.; Sutcliffe, J.V. River flow forecasting through conceptiual models part I—A discussion of principles. J. Hydrol. 1970, 10, 282–290. [Google Scholar] [CrossRef]
  45. McCuen, R.H.; Knight, Z.; Cutter, A.G. Evaluation of the Nash–Sutcliffe Efficiency Index. J. Hydrol. Eng. 2006, 11, 597–602. [Google Scholar] [CrossRef]
  46. Abdi, H.; Williams, L.J. Principal component analysis. Wiley Interdiscip. Rev. Comput. Stat. 2010, 2, 433–459. [Google Scholar] [CrossRef] [Green Version]
  47. Knudson, E.J.; Duewer, D.L.; Christian, G.D.; Larson, T.V. Application of Factor Analysis to the Study of Rain Chemistry in the Puget Sound Region; ACS Symposium Series; Amercan Chemestry Society: Washington, DC, USA, 1977; pp. 80–116. [Google Scholar]
  48. Wang, Y.; Shao, M.A.; Liu, Z.; Horton, R. Regional-scale variation and distribution patterns of soil saturated hydraulic conductivities in surface and subsurface layers in the loessial soils of China. J. Hydrol. 2013, 487, 13–23. [Google Scholar] [CrossRef]
  49. Abou Zakhem, B.; Al-Charideh, A.; Kattaa, B. Using principal component analysis in the investigation of groundwater hydrochemistry of Upper Jezireh Basin, Syria. Hydrol. Sci. J. 2017, 62, 2266–2279. [Google Scholar] [CrossRef]
  50. George, H.H.; Zohrab, A.S. Reference Crop Evapotranspiration from Temperature. Appl. Eng. Agric. 1985, 1, 96–99. [Google Scholar] [CrossRef]
  51. Berti, A.; Tardivo, G.; Chiaudani, A.; Rech, F.; Borin, M. Assessing reference evapotranspiration by the Hargreaves method in north-eastern Italy. Agric. Water Manag. 2014, 140, 20–25. [Google Scholar] [CrossRef]
  52. Lu, G.Y.; Wong, D.W. An adaptive inverse-distance weighting spatial interpolation technique. Comput. Geosci. 2008, 34, 1044–1055. [Google Scholar] [CrossRef]
  53. Dudgeon, C.R. Effects of non-Darcy flow and partial penetration on water levels near open-pit excavations. In Proceedings of the 18th Congress of the International Association of Hydrogeologists, Cambridge, UK, 8–13 September 1985; pp. 122–132. [Google Scholar]
  54. Younger, P.L.; Banwart, S.A.; Hedin, R.S. Mine Water. Hdrology, Pollution, Remediation; Springer Science & Business Media: London, UK, 2002. [Google Scholar]
  55. Buttle, J.M.; Hazlett, P.W.; Murray, C.D.; Creed, I.F.; Jeffries, D.S.; Semkin, R. Prediction of groundwater characteristics in forested and harvested basins during spring snowmelt using a topographic index. Hydrol. Process. 2001, 15, 3389–3407. [Google Scholar] [CrossRef]
  56. Price, K. Effects of watershed topography, soils, land use, and climate on baseflow hydrology in humid regions: A review. Prog. Phys. Geogr. Earth Environ. 2011, 35, 465–492. [Google Scholar] [CrossRef]
  57. Moriasi, D.N.; Arnold, J.G.; Van Liew, M.W.; Bingner, R.L.; Harmel, R.D.; Veith, T.L. Model evaluation guidelines for systematic quantification of accuracy in watershed simulations. Am. Soc. Agric. Biol. Eng. 2007, 50, 885–900. [Google Scholar]
  58. Richard, K. Development of Groundwater Recharge Model for the Sumanpa Catchment at Ashanti-Mampong-Ashanti Area in Ghana. Sci. Res. 2015, 3, 289. [Google Scholar] [CrossRef]
  59. Krishnaswamy, J.; Bonell, M.; Venkatesh, B.; Purandara, B.K.; Rakesh, K.N.; Lele, S.; Kiran, M.C.; Reddy, V.; Badiger, S. The groundwater recharge response and hydrologic services of tropical humid forest ecosystems to use and reforestation: Support for the “infiltration-evapotranspiration trade-off hypothesis”. J. Hydrol. 2013, 498, 191–209. [Google Scholar] [CrossRef]
  60. Wolock, D.M.; Winter, T.C.; McMahon, G. Delineation and Evaluation of Hydrologic-Landscape Regions in the United States Using Geographic Information System Tools and Multivariate Statistical Analyses. Environ. Manag. 2004, 34, S71–S88. [Google Scholar] [CrossRef]
  61. Kim, J.H.; Jackson, R.B. A Global Analysis of Groundwater Recharge for Vegetation, Climate, and Soils. Vadose Zone J. 2012, 11. [Google Scholar] [CrossRef] [Green Version]
  62. Jan, C.-D.; Chen, T.-H.; Lo, W.-C. Effect of rainfall intensity and distribution on groundwater level fluctuations. J. Hydrol. 2007, 332, 348–360. [Google Scholar] [CrossRef]
Figure 1. Ergene catchment geographic context and available gauging stations.
Figure 1. Ergene catchment geographic context and available gauging stations.
Water 11 00653 g001
Figure 2. Aspect and slope maps of the Ergene River catchment.
Figure 2. Aspect and slope maps of the Ergene River catchment.
Water 11 00653 g002
Figure 3. Base flow index method for total runoff-separation in mGROWA [11].
Figure 3. Base flow index method for total runoff-separation in mGROWA [11].
Water 11 00653 g003
Figure 4. Land use types of the Ergene River catchment for the year 2012.
Figure 4. Land use types of the Ergene River catchment for the year 2012.
Water 11 00653 g004
Figure 5. Soil groups in the Ergene River catchment.
Figure 5. Soil groups in the Ergene River catchment.
Water 11 00653 g005
Figure 6. mGROWA simulated water budget quantities in the Ergene catchment.
Figure 6. mGROWA simulated water budget quantities in the Ergene catchment.
Water 11 00653 g006
Figure 7. Long-term monthly groundwater recharge for the period 1991–2010.
Figure 7. Long-term monthly groundwater recharge for the period 1991–2010.
Water 11 00653 g007
Figure 8. Average annual water budget quantities as in terms of CORINE_2012 land-use classifications.
Figure 8. Average annual water budget quantities as in terms of CORINE_2012 land-use classifications.
Water 11 00653 g008
Figure 9. Average annual water budget quantities in terms of soil groups in the Ergene catchment.
Figure 9. Average annual water budget quantities in terms of soil groups in the Ergene catchment.
Water 11 00653 g009
Figure 10. Descriptive statistics of base flow by subcatchments.
Figure 10. Descriptive statistics of base flow by subcatchments.
Water 11 00653 g010
Figure 11. Comparison of the yearly long-term (observed-simulated) total runoff and groundwater recharge for the period 1991–2010 in ten subcatchments.
Figure 11. Comparison of the yearly long-term (observed-simulated) total runoff and groundwater recharge for the period 1991–2010 in ten subcatchments.
Water 11 00653 g011
Table 1. Base flow index (BFI) values of the unconsolidated rock areas according to Dörhöfer and Josopait (1980) [31], Hennings (2000) [32], and Wessolek and Facklam (1997) [33].
Table 1. Base flow index (BFI) values of the unconsolidated rock areas according to Dörhöfer and Josopait (1980) [31], Hennings (2000) [32], and Wessolek and Facklam (1997) [33].
Degree of SealingGroundwater DepthWaterlogging TendencySlopeBaseflow Indices
<2 mNo water logging<1%1
1.3–2 m1 (very low)1–3.5%0.9
I (10–45%) 0.82
3.5–7%0.67
7–10%0.59
0.8–1.3 m2 (low)10–13%0.5
0.4–0.8 m3–4 (medium – high)13–15%0.44
<0.4 m5 (very high)>15%0.4
II (45–75%) 0.33
III (75–90%) 0.28
IV (>90%) 0.2
Table 2. The classification of the hydrogeological properties of the hard rocks in North Rhine-Westphalia and the associated BFI values obtained by calibration [20].
Table 2. The classification of the hydrogeological properties of the hard rocks in North Rhine-Westphalia and the associated BFI values obtained by calibration [20].
Hydrogeological ClassPermeabilityHydraulic ConductivityBaseflow Indices
IVery high>10−2 m/sec0.9
IIHigh>10−3–10−2 m/sec0.6
IIIMedium>10−4–10−3 m/sec0.57
IVModerate>10−5–10−4 m/sec0.3
VLow>10−7–10−5 m/sec0.29
VIVery low>10−9–10−7 m/sec0.18
VIIExtremely low<10−9 m/sec0.12
Table 3. Database of the mGROWA water balance model for the Ergene catchment.
Table 3. Database of the mGROWA water balance model for the Ergene catchment.
DataData BaseScale/Spatial ResolutionsData Source
Climate data Precipitation (1991 to 2010)
Temperature (1991 to 2010)
100 × 100 mGeneral Directorate of State Hydraulic Works (DSI)
Soil dataMain soil groups1/25,000National soil database by General Directorate of State Hydraulic Works (DSI)
The depth of the groundwater table,
Perching water influence
Root depth, effective field capacity
Derived based on pedo-transfer functionsNational soil database by General Directorate of State Hydraulic Works (DSI) and fields studies
Capillary rise from the groundwater table
Land coverLand use types
Percentage imperviousness
1/25,000National soil database by General Directorate of State Hydraulic Works (DSI)
HydrogeologyGeological map1/500,000General Directorate of State Hydraulic Works (DSI) and fields studies
ReliefDigital elevation model SRTM
Ground surface slope
Ground surface exposition
30 × 30 mSRTM/X-SAR by National imaging and mapping agency
Runoff Station data Monthly and daily resolutionGeneral Directorate of State Hydraulic Works (DSI) and fields studies
Table 4. Catchment attributes, units, and their Pearson correlation with groundwater recharge.
Table 4. Catchment attributes, units, and their Pearson correlation with groundwater recharge.
Catchment AttributesAbbreviationUnitsCorrelation with the Groundwater Recharge
Soil types
Alluvial soilAS%0.671
Brown forest soilsBFS%−0.293
Brown forest soils without limeBFSWL%0.752
RendzinasR%−0.277
Brown soils without limeBSWL%−0.020
VertisolV%0.302
Urban soilsUS%−0.062
Land cover characteristics
Artificial surface LandASL%0.130
Agriculture areaAA%0.676
ForestF%0.659
Semi natural areaANA%−0.300
Wetland and water bodiesWWB%−0.072
Area and topography
Basin drainage areaBDAKm20.541
Basin slopeBS%−0.058
ElevationEM−0.469
Depth to water tableDWTM0.185
Climate
PrecipitationPmm/year0.817
Actual evapotranspirationAEmm/year0.686
Note: values in bold represent correlation with a significance level alpha = 0.05.
Table 5. Correlation analysis of catchment characteristics.
Table 5. Correlation analysis of catchment characteristics.
VariablesGRBSWLASBFSVBFSWLUSRASAAFANAWWBBDABSPAEDWTE
GR1.00
BSWL−0.021.00
AS0.670.211.00
BFS−0.29−0.610.111.00
V0.30−0.04−0.16−0.431.00
BFSWL0.75−0.23−0.30−0.12−0.531.00
US−0.060.200.47−0.21−0.250.231.00
R−0.28−0.30−0.420.41−0.200.07−0.161.00
AS0.130.430.55−0.470.15−0.130.62−0.231.00
AA0.680.410.33−0.340.730.85−0.08−0.460.341.00
F0.66−0.44−0.370.380.720.830.020.47−0.421.001.00
ANA−0.300.200.18−0.33−0.160.36−0.03−0.160.18−0.240.201.00
WWB−0.070.29−0.15−0.570.45−0.08−0.27−0.32−0.160.30−0.290.571.00
BDA0.540.240.53−0.280.14−0.150.28−0.150.370.34−0.36−0.030.061.00
BS−0.06−0.390.070.49−0.560.43−0.25−0.17−0.59−0.540.580.350.08−0.241.00
P0.82−0.220.510.59−0.01−0.53−0.35−0.25−0.330.34−0.29−0.17−0.070.190.441.00
AE0.69−0.42−0.600.25−0.430.69−0.300.300.720.800.830.250.11−0.570.64−0.191.00
DWT0.190.310.12−0.520.70−0.500.06−0.130.680.640.690.150.280.120.69−0.260.681.00
E0.47−0.58−0.390.45−0.570.720.020.40−0.450.850.87−0.24−0.55−0.230.43−0.170.680.751.00
Note: values in bold represent correlation with a significance level alpha = 0.05.
Table 6. Principal component analysis (PCA) of catchment attributes that are most correlated with groundwater recharge.
Table 6. Principal component analysis (PCA) of catchment attributes that are most correlated with groundwater recharge.
Statistical ProcedureIndicatorsPC1PC2
PC AnalysisEigenvalue3.9461.090
Variability (%)65.77418.174
Cumulative (%)65.77483.948
Variables
AS 0.5820.621
BFSWL 0.8940.087
AA 0.939−0.277
F 0.9410.289
P 0.5160.719
AE 0.8800.139
Note: bold type indicates variables with a loading factor greater than 0.70.

Share and Cite

MDPI and ACS Style

Rukundo, E.; Doğan, A. Dominant Influencing Factors of Groundwater Recharge Spatial Patterns in Ergene River Catchment, Turkey. Water 2019, 11, 653. https://doi.org/10.3390/w11040653

AMA Style

Rukundo E, Doğan A. Dominant Influencing Factors of Groundwater Recharge Spatial Patterns in Ergene River Catchment, Turkey. Water. 2019; 11(4):653. https://doi.org/10.3390/w11040653

Chicago/Turabian Style

Rukundo, Emmanuel, and Ahmet Doğan. 2019. "Dominant Influencing Factors of Groundwater Recharge Spatial Patterns in Ergene River Catchment, Turkey" Water 11, no. 4: 653. https://doi.org/10.3390/w11040653

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