Next Article in Journal
Public Participation in Water Planning in the Ebro River Basin (Spain) and Tucson Basin (U.S., Arizona): Impact on Water Policy and Adaptive Capacity Building
Previous Article in Journal
Water Sensitive Urban Design: An Investigation of Current Systems, Implementation Drivers, Community Perceptions and Potential to Supplement Urban Water Services
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Multi-Model Approach Using Statistical Index and Information Criteria to Evaluate the Adequacy of the Model Geometry in a Fissured Carbonate Aquifer (Italy)

1
School of Science and Technology, University of Camerino, Geology Division, Via Gentile III da Varano, 62032 Camerino (MC), Italy
2
DISTeVA, Dpartment of Earth, Life and Environmental Sciences, Via Cà le Suore, 2/4, 61029 Urbino (PU), Italy
3
Faculty of Mining, Geology and Petroleum Engineering, Department of Geology and Geological Engineering, University of Zagreb, Pierottijeva 6, Zagreb 10000, Croatia
*
Author to whom correspondence should be addressed.
Water 2016, 8(7), 271; https://doi.org/10.3390/w8070271
Submission received: 30 March 2016 / Revised: 21 June 2016 / Accepted: 23 June 2016 / Published: 28 June 2016

Abstract

:
A conceptual model related to a mountain aquifer that is characterized by a lack of data of hydrogeological parameters and boundary conditions, which were based on a single available observational dataset used for calibration, was studied using numerical models. For the first time, a preliminary spatial-temporal analysis has been applied to the study area in order to evaluate the real extension of the aquifer studied. The analysis was based on four models that were characterized by an increasing degree of complexity using a minimum of two zones and a maximum of five zones, which consequently increased the number of adjustable parameters from a minimum of 10 to a maximum of 22, calibrated using the parameter estimation code PEST. Statistical index and information criteria were calculated for each model, which showed comparable results; the information criteria indicated that the model with the low number of adjustable parameters was the optimal model. A comparison of the simulated and observed spring hydrographs showed a good shape correspondence but a general overestimation of the discharge, which indicated a good fit with the rainfall time series and a probably incorrect extension of the aquifer structure: the recharge contributes more than half of the total outflow at the springs but is not able to completely feed the springs.

1. Introduction

The application of numerical modeling to the study of fractured systems is an interesting field of scientific study [1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16] in relation to the complexity of their conceptualizations and to the quantity and quality of available data that are required [4,17].
The spatial distributions of hydraulic conductivities, fault and joint system geometries and their interconnections are still important factors to explore, and, moreover, fractured aquifers are usually located in mountainous areas characterized by complex geological settings and often by the scarcity or total absence of data regarding their physical-mechanical and/or hydrogeological parameters.
There are several methodologies [18,19] useful for modeling fractured systems, but each method is necessarily conditioned by available data and consequently produces uncertainties on the numerical analysis.
Uncertainty is the key issue in hydrogeological modeling, and it is related to the quantity and quality of data, scale of choice, boundary conditions and initial conditions. In addition, the models are often subject to errors associated with the geometry and calibration process [20,21,22].
In some cases, the hydrogeological parameters involved in the simulation are assessed using a laboratory analysis, and, of course, those analyses cannot be immediately transferred to numerical models but require up-scaling. The up-scaling process alone is possibly insufficient to transfer the small-scale parameters to valid distribution parameters at a basin scale, and errors consequently will then be incurred. Consider also that models with multiple parameters increase uncertainty because the observational information is distributed to all of the parameters.
Uncertainty can be reduced using parsimonious models [20], which are characterized by a limited number of adjustable parameters. In the literature, there are several approaches to this problem based on statistical methods [20,23,24,25,26,27] in order to understand the validity of model geometries and the spatial distributions of parameters.
In groundwater modeling analyses, parallel to the calculation of the main statistical index (e.g., variance and standard error), the evaluation of the information criteria is the most useful [20,28,29,30,31,32,33,34]. One of these criteria is the Akaike information criterion (AIC) [35,36] that gives a relative distance value between the fitted model and the unknown true model [37]; this criterion is based on the Kullback-Leibler information to calculate the relative distance. The other methods are the Bayesian information criterion (BIC) [38] and Kashyap’s information criterion (KIC) [39]; the first method was developed because the AIC sometimes requires more parameters than needed, and the second method was developed on the basis of Bayesian theory but instead assumes an asymptotic approximation to the model likelihood [40]. In general, these indexes become large as more parameters are added, and smaller values indicate a more accurate model [20]. Carrera and Neumann [23] first successfully applied those parameters to discriminate between different parameter sets of a groundwater model.
The present study considers a mountain aquifer that is characterized by a limited number of available data and a consequent lack of information on the boundary and initial conditions. The observation data are related to the spring discharges and were collected for approximately 10 years; these are the only available data for model calibration.
The analysis was conducted by increasing the number of zones in the model domain and consequently the number of adjustable parameters, and statistical index and information criteria were calculated for each model under consideration.
For the first time, this methodology was applied to the study area to carried out a preliminary spatial-temporal analysis in order to evaluate the real extension of the aquifer studied.

2. Materials and Methods

2.1. Study Area

The study area is located in the northeast sector of the Sibillini Mountains, which is a calcareous massif located in the outermost part of the Umbria-Marche Ridge, and the study area corresponds to the mountainous portion of the Tennacola stream catchment between Mount Valvasseto (1526 m a.s.l.) and Mount Castel Manardo (1977 m a.s.l.) (Figure 1a).
The outcropping geological formations are mainly composed of calcareous, marly and siliciclastic lithotypes (from Trias to Miocene in age) that belong to the Umbria-Marche Succession [41,42]. From bottom to top, the following formations can be observed (Figure 2):
-
Calcare massiccio formation (Hettangian-Pliensbachian p.p.): stratified limestone in tough blocks without noticeable vertical heterogeneities, approximately 700 m thick.
-
Bugarone group (Pliensbachian p.p.-Titonian p.p.): comprised of limestones and marly limestones with frankly greenish intercalations of marls in the marly calcareous facies; it outcrops only locally in substitution of the three previous formations, with an overall thickness of 30–40 meters.
-
Corniola formation (Sinemurian p.p.-Toarcian p.p.): finely stratified micritic limestones, with brown or light grey chert and with grey-green pelitic intercalations that are rather abundant both in the upper part and at the base; the thickness varies from 150 to 400 meters.
-
Bosso formation (Toarcian p.p.-Bajocian p.p.): comprised of a marly and marly calcareous unit with medium-thick strata (Marne del Serrone and Rosso Ammonitico); a third unit is made of limestones and marly limestones (Calcari a Posidonia). The overall thickness varies from 40 to 120 meters.
-
Calcari Diasprini formation (Bajocian p.p.-Titonian p.p.): limestones and cherty limestones with thick stratification; thickness can reach 100 meters.
-
Maiolica formation (Titonian p.p.-Aptian p.p.): whitish micritic limestones, in medium strata, with thin pelitic intercalations that increase towards the top of the formation. The thickness ranges from 150 and 400 meters.
-
Marne a Fucoidi formation (Aptian p.p.-Albian p.p.): the lower portion is made up of multicolored, closely layered marls and, subordinately, limestones that towards the top become cherty with frequent marly levels; 80–100 meters thick.
-
Scaglia rosata formation (Albian p.p.-Bartonian p.p.): limestones and marly limestones with polychrome chert, in medium and thin layers that vary in color from pink to red, (200–400 m);
-
Scaglia cinerea formation (Bartonian p.p.-Aquitanian p.p.): marly limestones, calcareous marls and clayey marls in thin layers (150–200 m).
-
Bisciaro formation (Aquitanian p.p.-Burdigalian p.p.): limestones, cherty limestones, marly limestones and calcareous and clayey marls, in medium and, sometimes, thick layers; the thickness ranges from 20 and 60 meters.
The current structural setting is the result of alternating compressive and extensive tectonic phases. The last configuration is a direct consequence of a former compressive phase that developed between the Miocene and early Pliocene that created east-verging folds and thrusts within the bedrock. Subsequently, a second phase, which had been active since the early Pleistocene, was characterized by an extensive regime associated with an intense tectonic uplift; this last phase generated vast active fault systems towards the Apennines, and the systems sometimes exploited the same old thrust planes [43,44]. The geological context described above is greatly complicated by the presence of the Sibillini thrust, which is characterized by a curved geometry and an approximately NW-SE direction (Figure 1a).
The presence of lithotypes with different permeabilities (aquifers/aquicludes/aquitards) allowed the formation of three main hydrogeological complexes that overlapped and are characterized by springs with different discharges and regimes (Figure 1a,b and Figure 2). The springs investigated herein emerge at different heights above the present thalweg, on the hydrographic right of the Tennacola stream.
Detailed surveys performed during the last three years have allowed the springs to be differentiated based on their proximities and hydrogeological characteristics (Figure 1b) into two main groups:
-
the Tennacola high group of springs; and
-
the Tennacola low group of springs.
The springs were captured and measured for groups of springs by the Tennacola Water Consortium (C.I.I.T.); The University of Camerino began detailed monitoring in 2014 by instrumenting the Tennacola low group of springs for the purpose of evaluating the single contribution.
In this study, only the Tennacola low group of springs, which is representative of the system under study, were considered. The Tennacola low group of springs emerges between approximately 700 m and 765 m a.s.l., at the contact between the Scaglia aquifer and the Scaglia Cinerea impermeable formation, and it represents the most important resource in the area, with discharges that range from 0.050 to 0.250 m3/s. In this case, the Scaglia aquifer likely forms a unique hydrogeological complex with the Maiolica aquifer, which is stratigraphically reversed here (Figure 1b). The Marne a Fucoidi aquiclude, which is usually interposed between the two aquifers, is partially elided or absent due to the compressive tectonics and the presence of the thrust plane. Moreover, preliminary water budgets led us to additionally hypothesize a minor contribution (approximately 10%) from the Basal aquifer on the western side of the area (Figure 1b).

2.2. Water Budget

The water budget was evaluated for a period of 10 hydrologic years from 2004 to 2014 using meteorological data recorded at the Pintura di Bolognola weather station and fluvial and spring discharge data (Figure 3, Figure 4 and Figure 5).
The evapotranspiration was evaluated using different methods [45,46,47,48,49] to evaluate the differences in the methods.
The results (Table 1) are comparable to those of other studies carried out at a larger scale [50]; in particular the evapotranspiration ranged from 390 mm/year to 624 mm/year (Figure 6) with an effective precipitation ranging from 890 mm/year to 1124 mm/year.
The methods of Turc and Thornthwaite provided similar average values of approximately 500 mm/year, while the daily Hargreaves-Samani and Penman-Monteith methods yielded average values of approximately 400 mm/year and 620 mm/year, respectively; the Penman-Monteith method, recommended by the Food and Agriculture Organization of the United Nations (F.A.O.), can be, however, a double-edged weapon because the greater amount of input measured data could lead to larger uncertainty. The evapotranspiration value was not considered during negligible conditions that avoid the evapotranspiration process [47].
The runoff was evaluated through the use of bibliographic and monitoring data; the first fluvial discharges were measured by Perrone [51], who considered river sections located downstream of the study area. Over the years, another analysis, made by the University of Rome La Sapienza [50,52,53], showed a negligible runoff.
Three measurement campaigns were conducted by the University of Camerino along the Tennacola stream between August 2013 and January 2014 in order to monitor the river discharge during the dry and wet seasons (Figure 4). For the monitoring, 16 measuring stations were used that were evenly distributed from 1100 and 600 m a.s.l., which corresponded to the highest and lowest group of analyzed springs. The data collected during the summer periods (August 2013) showed that the river became nearly dry along the investigated reach. In fact, only limited sectors of the river (a few tens of meters) had discharges of 0.002–0.003 m3/s that, however, were suddenly absorbed by the large thicknesses of the coarse materials (pebbles and blocks) that form the riverbed. Slightly more significant discharges (approximately 0.006 m3/s) were recorded at the lowest station, which was located at an altitude of 600 m a.s.l. These conditions certainly persisted until the first half of October when, during the periodical surveys, no significant changes were observed.
The situation changed significantly after the heavy rainfall that was recorded in the area on 11, 12 and 13 November 2013. The second campaign of measurements (December 2013) showed significant increases in fluvial discharge. Although it was not possible to reach the first three stations (due to high snow avalanche and landslide hazards), the monitoring detected significant increases in the fluvial reach between 875 and 860 m a.s.l. However, these quantities were rapidly absorbed from the riverbed that, up to an altitude of 747 m a.s.l., showed considerable discharge variations, with alternating reaches completely dry or with reduced flows. The major increases were recorded starting from an altitude of 723 m a.s.l. (in correspondence with the Tennacola low group of springs), where the values reached 0.320 m3/s and then decreased significantly (0.287 m3/s) at the river station at 675 m a.s.l. (due to the large amounts of coarse materials in the river bed). After that, the flow increased again at the last station (altitude 600 m a.s.l.), where the fluvial discharge exceeded 0.450 m3/s. A similar trend was observed during the January 2014 measurement campaign.
The bibliographic and monitoring data showed negligible or no runoff, and any runoff accompanied snow melting and significant weather events (as confirmed by measurements made during December 2013); the river is therefore probably connected to the piezometric field, as observed at river station 713 (Figure 4), and flow under the debris is associated with the formation of a riverbed aquifer.
To elaborate the water budget, it was therefore decided to assume negligible runoff and apply a coefficient of potential infiltration (C.I.P.) of 95% [54]. The effective infiltration (Table 2) ranged on average from 843.87 mm/year to 1068.19 mm/year. Considering the bibliographic [50,52] and monitoring data, a maximum effective infiltration of approximately 1068 mm/year was assumed.
A recession analysis of the springs was conducted examining their depletion curves for the period 2004–2014. The analysis showed that the system’s depletion constant ranged from 0.003 to 0.008 with an average value of 0.006. Based on those results, the groundwater recharge was estimated from the spring discharges for the same period [55].
The Tennacola low group of springs from 2004 to 2014 was characterized by a dynamic change in reserves of approximately 827,471 m3, with a recharge of approximately 55,515,801 m3 and an annual average flow of approximately 5,046,891 m3/year (Table 3).

2.3. Conceptual Model

Based on the geological and hydrogeological characteristics described so far, a conceptual model for the study area was defined (Figure 1a,b):
-
a western side, where the Maiolica and Calcare Massiccio aquifers are separated by an aquiclude comprised of the Bugarone, Calcari a Posidonia and Calcari Diasprigni formations: based on the water budget discussed in the previous section, a partial connection between the two aquifers (with a very limited contribution from the basal aquifer) was considered;
-
an eastern side, where the thrust plane and the underlying Scaglia cinerea low permeability formation, constitute an impermeable barrier;
-
a northern side, where the Tennacola low group of springs is located, and with no contribution from the Tennacola stream; and
-
a southern side, where inflow was excluded because of the presence of a structural underground water divide.
The main aquifer is comprised of the Maiolica and Scaglia formations; the Marne a Fucoidi formation, which is usually interposed between them, was not detected during the monitoring, probably because it has been partly erased by the thrust. For this reason, here it must be considered as an aquitard that produces a general delay in the flow of groundwater to the springs.
The spring hydrograph analysis confirmed the complexity of the whole system. The low value of the depletion constant, α = 0.006, and the curve homogeneity describe a well-defined catchment (also confirmed by the strict correlation between precipitation spring discharge) and a groundwater circuit that mainly occurs as a gradual flow through a widespread network of fractures and a low permeability matrix [56,57]. In relation to the features of the geological formations, the dimensions of the study area and the results from the spring hydrograph analysis, the equivalent continuum approach was applied in the present study, with the same value for the fractures and matrix system.
The geological and hydrogeological monitoring allowed a recharge area of 3.8 km2 to be considered (Figure 1a); therefore, in relation to the water budget and annual discharge at the springs, the effective infiltration was determined to range from 70% to 80% of the mean yearly recharge.
The role of the river is uncertain. The fluvial discharge monitoring carried out during 2012–2013 evidenced that the river is completely dry for most of the year, except for a few days following snow melt and extreme rainfall events; this fact can be related to the high amount of coarse debris present along the river bed. Therefore, the river flow under the debris is probably connected to the groundwater flow.

2.4. Model Design

The numerical model related to the Tennacola low group of springs was developed using the FEFLOW [58], which is widely used to describe groundwater flow and transport processes in porous and fractured media.
As discussed in the previous sections, the model domain is very complex because of the presence of numerous lithological classes and structural elements. Therefore, the local contexts in the numerical model involve simplifications.
The study area can be divided into five different sectors (Figure 7):
-
Sector 1, which is characterized by two overlapping aquifer systems: the Maiolica aquifer (sMai) on the top and the Basal aquifer (sBas) below, which are separated by low permeability formations;
-
Sector 2, which corresponds to the Maiolica aquifer (sMai);
-
Sector 3, which is characterized by the Marne a Fucoidi aquiclude (sFuc);
-
Sector 4, which is characterized by the Scaglia aquifer (sScr); and
-
Sector 5, which corresponds to the valley floor of the Tennacola stream thick fluvial deposits (sTen).
The modeling approach used four different simulations of increasing complexity (Figure 7):
-
two zone model, M2Z (sBas, sMai-sFuc-sScr-sTen);
-
three zone model, M3Z (sBas, sMai-sFuc-sScr, sTen);
-
four zone model, M4Z (sBas, sMai-sScr, sFuc, sTen); and
-
five zone model, M5Z (sBas, sMai, sFuc, sScr, sTen).
The model domain, which was comprised of one layer, extended over an area of approximately 3.8 km2 and was discretized into a mesh of 2085 elements and 2290 nodes. Because stratigraphic and water table data were missing in the study area, the initial piezometric surface was reconstructed using interpretative hydrogeological cross-sections, bibliographic data and the results of fluvial discharge monitoring carried out in the mountainous portion of the Tennacola stream.

2.5. Hydrogeological Properties

The hydrogeological parameters (Table 4) were estimated from bibliographic data [59,60,61,62], in relation to the features of the study area. The calibration was focused on conductivity, specific yield and transfer rate parameters.
The range of hydrogeological parameters related to the five zones is shown in Table 4. The parameters have been made vary during the calibration process within a specific range (indicated with min and max values in Table 4) starting from a reference value (indicated with ref in the same table). The number of zones was progressively increased from two to five and the range was refined in relation to the results of previous simulations.

2.6. Boundary Conditions

The choice of boundary conditions (BC) followed the conceptual model described above (Figure 8). No-flow boundaries were used for the western, eastern and southern sides of the study area, but the Tennacola stream was represented using a constant head boundary (Cauchy BC). The infiltration was assumed to be time varying and was applied homogeneously to the entire study sector, which was assumed to be the recharge area. The numerous springs were grouped together and simulated using a Seepage Face (Dirichlet BC with a constraint), and the outflow from the model occurred only when the water table exceeded the elevation of the springs.

2.7. Model Calibration

The calibration was performed using the parameter estimation code PEST [63], which is an automated model calibration that uses an inverse parameter estimation process based on the Gauss-Marquardt-Levenberg method. This methodology allowed us to find an optimal parameter distribution to reduce the differences between the observed and simulated data by minimizing the sum of the squared weighted residuals.
The residual analysis allowed the quality of the model to be evaluated from statistical indexes such as the sum of squared weighted residuals, Φ, and the standard error, σ [20,64].
The sum of the squared weighted residuals between the observed and simulated values is expressed by the following relationship:
Φ = ( c c 0 J ( b b 0 ) ) t Q ( c c 0 J ( b b 0 ) ) = i = 1 m [ w i ( c i c i ' ) ] = i = 1 m ( w i r i ) 2
where wi is the weight of the ith-observation; c i is the ith-observation; c i ' is the simulated value corresponding to c i ; and r i is the residual value.
The standard error can be calculated starting from the variance σr2 [64]:
σ r 2 = Φ min ( n m )
where Φ min is the minimum objective function; n is the observations number; and m is the number of adjustable parameters.
The square root of the previous equation gives the standard error, σ:
σ = σ r 2 w i
The aforementioned indexes cannot adequately indicate the convergence of a model, especially as the number of adjustable parameters increases. To ensure a greater precision, four information criteria were employed [20,35,36,37,64,65]:
AIC = n ln ( σ r 2 ) + 2 k
AICc = AIC + 2 k ( k + 1 ) n k 1
BIC = n ln ( σ r 2 ) + k ln ( n )   and
KIC = ( n ( k 1 ) ) ln ( σ r 2 ) ( k 1 ) ln ( 2 π ) + ln | J t QJ |
where k is the number of estimable parameters; J is the Jacobian matrix; Q is the weight matrix; and JtQJ is the normal matrix.
If the statistics for a model with few parameters are better than a model with more parameters, it is preferable to select the model with few parameters, unless there is no information on the validity of the more complex model [20].
These parameters are usually evaluated through the difference Δj of the AIC, which is also generally used to indicate differences in the other parameters, BIC and KIC [37]. This value is then used together with the Akaike weight:
AIC w j = exp ( 1 2 · Δ j ) j = 1 R exp ( 1 2 · Δ j )
where Δj is the difference of the AIC of a specific model to the smallest AIC from all of the considered models.
The value wj indicates the best model among all of the evaluated models [37]. The difference value is inversely proportional to the weight wj, and for low values of the weight, the probability that the considered model is the best is null.
The good fit between the simulated and observed data can be graphically evaluated through the correlation coefficient, R [20,63]:
R = ( w i c i m ) ( w i c i m ) [ ( w i c i m ) ( w i c i m ) ( w i c i m ) ( w i c i m ) ] 1 / 2
where c i is the ith-observation value; c i ' is the model-generated counterpart to the ith-observation value; m is the mean of the weighted observed values; m is the mean of the weighted model-generated values; and wi is the weight associated with the ith-observation.
The correlation coefficient is independent of the number of observations and the uncertainty level associated with the observations.
The PEST software allows the sensitivity of the selected parameters against the observed data to be evaluated. The PEST sensitivity for a generic parameter i is calculated as follows [63]:
s i = ( J t QJ ) ii m
where m is the number of observations with a nonzero weight.
The composite sensitivity sj (related to observation j) is instead calculated as follows [63]:
s j = { Q ( J t J ) } jj n
where JtJ is the Hessian matrix; j is the number of observations; and n is the number of adjustable parameters.
The model calibration was performed using the spring discharge values from 2004 and 2014. In particular, automatic inverse modeling was carried out starting from the simple model (M2Z) to obtain the more complex model (M5Z) that had a larger number of adjustable parameters. The best model was finally chosen based on the derived statistical indexes.

3. Results

3.1. Uncertainty Analysis

The uncertainty analysis was performed using four models with increasing complexity, with a minimum of 10 to a maximum of 22 adjustable parameters.
Figure 9 shows the simulated vs. observed discharges for all of the models. Table 5 shows the results from the application of the main statistical indexes to the weighted residuals of each model. The weight used was evaluated considering the inverse of the standard deviation from the PEST analysis [20,66]. The reference variance should therefore be equal to 1 [20,63]; lower values represent a good fit between simulated and observed data, and higher values indicate a bad fit.
A detailed comparison between the simulated and observed hydrographs was made for the hydrologic years 2006, 2008, 2010 and 2012 (Figure 10); the figure shows that all of the models exhibited a similar trend and lack in the discharge (Table 6).
The regression analysis showed a bad fit between the simulated and observed data based on the slope of the best fit line (Figure 11).
The analysis of the AIC, BIC and KIC coefficients allowed a determination of the best of the four proposed configurations: the coefficient AICc, in this case equal to AIC, was not considered. The information criteria show relative values, and the absolute values are nonsense; the scale of the y-axis in Figure 12 is therefore omitted. Finally, the Akaike weight (Equation (8)) was calculated to evaluate the likelihood of the models; values of 100% indicate the best model, and a value close to zero indicates a negligible model (Table 7).

3.2. Sensitivity Analysis

The dimensional sensitivity was evaluated in relation to the observation group for each of the considered models (Figure 13).
The sensitivity analysis was very useful for evaluating the capacity of the models to represent the observed data.
For the proposed models, the sensitivity analysis showed comparable trends, but in particular the value was higher for the model with few parameters. In general, increasing the number of adjustable parameters decreased the sensitivity to the observation data within the model calibration [34].
All the proposed models showed a higher sensitivity for the observations matching the rainfall peaks (i.e., those of 5 August 2004 and 5 October 2013). In more detail, the models with three and four zones showed a higher sensitivity for the peak of rainfall on 5 October 2013, but the models with two and five zones were more sensitive for the peak on 05 August 2004.

4. Discussion

Based on the aforementioned analyses, it was possible to evaluate the goodness of the proposed models (Table 5). The variance, which was less than 1 in all of the considered models, indicates a good correspondence; the lowest value, in particular, corresponds to the model with two zones and a limited number of adjustable parameters. On the contrary, R, which was constant for all of the models, is apparently lower than the recommended value [66]. This result was also confirmed by the scatterline graph, in which the slope of the best-fit line was different than 1 (Figure 11). Even the simulated hydrographs (Figure 9) for all of the models had average errors from 0.038 and 0.040 m3/s, which were too high with respect to the observed discharges at the springs. The error is comparable among the models but slightly increases with the number of zones (and consequently with the number of adjustable parameters); this is due to the problem of equifinality, in which increasing the parameter space can increase the number of optimum solutions available and reducing the performance of all “optimums” [67]. As a consequence, the model with two zones and a low number of adjustable parameters remained the most suitable among the models investigated in relation to the available data. The analysis of the AIC, BIC and KIC parameters confirmed these results.
In addition, the AIC wj, which was 100%, confirmed the goodness of the model with two zones among the models tested; on the contrary, the other models had values close to zero (Table 7). Engelhardt et al. [34] carried out an analysis of the significance of the BIC, KIC and AIC parameters. The first two parameters indicated that the best model had the lowest number of adjustable parameters, whereas AIC indicated that the best model had a higher number of parameters; the Akaike weight definitively confirmed the result expressed by the AIC parameter. However, because AIC ignores the potential error related to the geometry of the model, the authors stated that the choosing the KIC parameter could be the best solution, which confirmed the result achieved by [40]. Foglia et al. [30] conducted a multi-model analysis using AIC and BIC supported by a cross-validation, and they concluded that both information criteria selected the same optimal model with only three adjusted parameters.
The solution obtained in the present study could be related to the use of a single type of observational data; in fact the information criteria considers different model as optimal if a large number of observational data are available [34].
The detailed analyses carried out for single hydrologic years (Figure 10) instead showed good approximations of the hydrograph shape (discrete spatial distribution of the hydrogeological parameters); nevertheless, the simulated discharge showed a general underestimation in the observed discharge of approximately 1,500,000 m3/year on average (Table 6).
The selected model simulated a mean annual discharge of approximately 44,021,000 m3 over the entire period of investigation, but the mean annual discharge calculated by means of the spring hydrograph analysis and water budget, on the other hand, indicated a value of 55,515,801 m3. The results confirmed that infiltration in the study area plays a very important role because it feeds approximately 80% of the total discharges to the springs, and, as a consequence, approximately 20% should come from a neighboring aquifer.

5. Conclusions

The analysis was conducted by increasing complexity through increasing the number of adjustable parameters using the FEFLOW computer program and PEST automatic calibration; for each of the models under evaluation, the differences between the simulated and observed data were evaluated using the main statistical index and information criteria to evaluate the associated errors and uncertainties. The optimal model was evaluated using the differences between the simulated and observed water budgets to understand the goodness of the hypothesized conceptual model.
All of the models under consideration had small values of the main statistical index as well as those for the reference variance and the correlation index, which highlights the complexity of the geometry that controls the flow at the springs. The differences between observed and simulated data ranges, on average, between 0.038 and 0.040 m3/s; therefore, it seems probably high compared to the discharges at the springs. In relation to the four models evaluated, the information criteria in parallel with the Akaike weight identified the model with two zones, i.e., the model with the smallest number of adjustable parameters, to be the most suitable. Using that model, a total discharge for the whole period 2004–2014 of approximately 44,021,000 m3 was determined, compared with the average total observed at the springs of approximately 54,688,330 m3, which implied a contribution from infiltration of approximately 80% and a minor contribution of approximately 20% that was probably derived from the nearby basal aquifer on the western side of the area.
The comparison between the simulated and observed spring hydrographs showed a good shape correspondence but a general overestimation of the discharge, which indicates a good correspondence with the rainfall time series and consequently with the estimated water budget and a probably incorrect extension of the aquifer structure. The statistical analysis showed similar results with comparable discharge differences for each model, which indicated the same lack in available data independent of the geometrical model. In conclusion, the study confirmed that recharge contributes more than half of the total outflow at the springs but is not able to completely feed the springs, and it is necessary to invoke an external contributor, which is probably flow from the close Basal aquifer on the western side.
In conclusion, the aim of the paper would be not to provide a definitive solution, rather to foster debate and restrict uncertainty about assessment of different alternative models. Therefore the study presented is preliminary and designed to assess, in relation of available data, the real extension of the model. The analysis allowed an understanding of the uncertainties and errors associated with the hypothesized conceptual model in relation to the available data and its consequent limitations, showed the necessity of acquiring new data (e.g., new fluvial discharge monitoring, at least once per year, refining the measured discharges at the springs using the instruments installed by the University of Camerino in 2014 and geochemical surveys) and modifying the conceptual model (e.g., extending the study area) in order to improve and verify the achieved results.

Acknowledgments

This study was supported by “2014 Academic funds for Research” of the University of Camerino. We thank the anonymous reviewers for the useful suggestion in the former phase of revision.

Author Contributions

The research presented here was carried out in collaboration between all authors. Marco Giacopetti and Ezio Crestaz conceived and designed the study; Marco Giacopetti performed the numerical and statistical analysis; Marco Giacopetti, Ezio Crestaz, Marco Materazzi, Gilberto Pambianchi and Kristijan Posavec analyzed the results; and Marco Giacopetti, Ezio Crestaz, Marco Materazzi and Kristijan Posavec wrote and revised the paper. All authors discussed the structure and commented on the manuscript at all stages.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

The following abbreviations are used in this manuscript:
AICAkaike information criterion.
AICcAkaike information criterion.
AIC wjAkaike weight.
BCBoundary condition.
BICBayesian information criterion.
C.I.I.T.Tennacola Water Consortium.
C.I.P.Potential Infiltration Coefficient.
FAOFood and Agriculture Organization of the United Nations.
FEFLOWFinite element groundwater flow and transport modeling tool.
KICKashyap’s information criterion.
RCorrelation coefficient.
PESTModel-Independent Parameter Estimation and Uncertainty Analysis.

References

  1. Kiraly, L. La régularisation de l’Areuse (Jura Suisse) simulée par modèle mathématique. Hydrogeol. Karstic Terrains 1984, 1, 94–99. [Google Scholar]
  2. Kiraly, L. Modelling karst aquifers by the combined discrete channel and continuum approach. Bull. d’Hydrogéol. 1998, 16, 77–98. [Google Scholar]
  3. Teutsch, G. An extended double-porosity concept as a practical modeling approach for a karstified terraine. In Hydrogeological Processes in Karst Terraines, Proceedings of the Antalya Symposium and Field Seminar, Antalya, Turkey, 7–17 October 1990.
  4. Sauter, M. Double porosity models in karstified limestone aquifers: Field validation and data provision. Hydrogeol. Process. Karst Terranes 1993, 207, 261–279. [Google Scholar]
  5. Boni, C.F.; Mastrorillo, L.; Preziosi, E. Simulazione numerica di acquiferi carbonatici: L’esempio della struttura, M. Maggio-M. Penna (Nocera Umbra). Geol. Romana 1994, 30, 27–36. [Google Scholar]
  6. Leonardi, V.; Arthaud, F.; Grillot, J.C.; Avetissian, V.; Bochnaghian, P. Modélisation d’un aquifère basaltique fracturé tenant compte de données géologiques, climatiques et hydrauliques: cas des basaltes perchés de Garni (Armenie). J. Hydrol. 1996, 179, 87–109. [Google Scholar] [CrossRef]
  7. Angelini, P.; Dragoni, W. The problem of modeling limestone springs: the case of Bagnara (north Appennines, Italy). Ground Water 1997, 35, 612–618. [Google Scholar] [CrossRef]
  8. Barrett, M.E.; Charbeneau, R.J. A parsimonious model for simulating of flow in a karst aquifer. J. Hydrol. 1997, 196, 47–65. [Google Scholar] [CrossRef]
  9. Eisenlohr, L.; Kiraly, L.; Bouzelboudjen, M.; Rossier, Y. Numerical simulation as a tool for checking the interpretation of karst spring hydrographs. J. Hydrol. 1997, 193, 306–315. [Google Scholar] [CrossRef]
  10. González-Herrera, R.; Sánchez-y-Pinto, I.; Gamboa-Vargas, J. Groundwater-flow modeling in the Yucatan karstic aquifer, Mexico. Hydrogeol. J. 2002, 10, 539–552. [Google Scholar] [CrossRef]
  11. Scanlon, B.R.; Mace, R.E.; Barrett, M.E.; Smith, B. Can we simulate regional groundwater flow in a karst system using equivalent porous media models? Case study, Barton Springs Edwards aquifer, USA. J. Hydrol. 2003, 276, 137–158. [Google Scholar] [CrossRef]
  12. Jaquet, O.; Siegel, P.; Klubertanz, G.; Benabderrhamane, H. Stochastic discrete model of karstic networks. Adv. Water Resour. 2004, 27, 751–760. [Google Scholar] [CrossRef]
  13. Cheng, J.; Chen, C.X. An integrated linear/non-linear flow model for the conduit-fissure-pore media in the karst triple void aquifer system. Environ. Geol. 2005, 47, 163–174. [Google Scholar] [CrossRef]
  14. Gurwin, J.; Lubczynski, M. Modeling of complex multi-aquifer systems for groundwater resources evaluation—Swidnica study case (Poland). Hydrogeol. J. 2005, 13, 627–639. [Google Scholar] [CrossRef]
  15. Nastev, M.; Lefebvre, R.; Rivera, A.; Martel, R. Numerical simulation of groundwater flow in regional rock aquifers, southwestern Quebec, Canada. Hydrogeol. J. 2005, 13, 835–848. [Google Scholar] [CrossRef]
  16. Rani, F.M.; Chen, Z.H. Numerical modeling of groundwater flow in karst aquifer, Makeng mining area. Am. J. Environ. Sci. 2010, 6, 78–82. [Google Scholar] [CrossRef]
  17. Ford, D.; Williams, P. Karst Hydrogeology and Geomorphology, 2nd ed.; Wiley: Chichester, UK, 2007; p. 562. [Google Scholar]
  18. Dietrich, P.; Helmig, R.; Sauter, M.; Hötzl, H.; Köngeter, J.; Teutsch, G. Flow and Transport in Fractured Porous Media, 1st ed.; Springer-Verlag Berlin Heidelberg: Berlin, Germany, 2005; p. 447. [Google Scholar]
  19. Stevanović, Z. Karst Aquifers—Characterization and Engineering, 1st ed.; Springer International Publishing: Cham, Switzerland, 2015; p. 698. [Google Scholar]
  20. Hill, M.C.; Tiedeman, C.R. Effective Calibration of Groundwater Models: With Analysis of Data, Sensitivities, Predictions, and Uncertainty, 1st ed.; Wiley: New York, NY, USA, 2007; p. 480. [Google Scholar]
  21. Doherty, J.; Welter, D. A short exploration of structural noise. Water Resour. Res. 2010, 46. [Google Scholar] [CrossRef]
  22. Sepúlveda, N.; Doherty, J. Uncertainty analysis of a groundwater flow model in East-central Florida. Ground Water 2015, 53, 464–474. [Google Scholar] [CrossRef] [PubMed]
  23. Carrera, J.; Neumann, S.P. Estimation of aquifer parameters under transient and steady state conditions: 1. Maximum likelihood method incorporating prior information. Water Resour. Res. 1986, 22, 199–242,. [Google Scholar] [CrossRef]
  24. Hernandez, A.F.; Neuman, S.P.; Guadagnini, A.; Carrera, J. Inverse stochastic moment analysis of steady state flow in randomly heterogeneous media. Water Resour. Res. 2006, 42. [Google Scholar] [CrossRef]
  25. Tsai, F.T.C.; Li, X. Inverse groundwater modeling for hydraulic conductivity estimation using Bayesian model averaging and variance window. Water Resour. Res. 2008, 44. [Google Scholar] [CrossRef]
  26. Kavetski, D.; Fenicia, F.; Clark, M.P. Impact of temporal data resolution on parameter inference and model identification in conceptual hydrological modeling: Insights from an experimental catchment. Water Resour. Res. 2011, 47. [Google Scholar] [CrossRef]
  27. Gupta, H.V.; Clark, M.P.; Vrugt, J.A.; Abramowitz, G.; Ye, M. Towards a comprehensive assessment of model structural adequacy. Water Resour. Res. 2012, 48. [Google Scholar] [CrossRef]
  28. Poeter, E.P.; Anderson, D.R. Multimodel ranking and inference in ground water modeling. Ground Water 2005, 43, 597–605. [Google Scholar] [CrossRef] [PubMed]
  29. Hill, M.C. The practical use of simplicity in developing ground water models. Ground Water 2006, 44, 775–781. [Google Scholar] [CrossRef] [PubMed]
  30. Foglia, L.; Mehl, S.W.; Hill, M.C.; Burlando, P. Evaluation of information criteria, linear statistics, cross validation, and predictions using alternative groundwater models. Ground Water 2007, 45, 627–641. [Google Scholar] [CrossRef] [PubMed]
  31. Parker, A.H.; West, L.J.; Odling, N.E.; Bown, R.T. A forward modeling approach for interpreting impeller flow logs. Ground Water 2010, 48, 79–91. [Google Scholar] [CrossRef] [PubMed]
  32. Singh, A.; Mishra, S.; Ruskauff, R. Model averaging techniques for quantifying conceptual model uncertainty. Ground Water 2010, 48, 701–715. [Google Scholar] [CrossRef] [PubMed]
  33. Ye, M.; Pohlmann, K.F.; Chapman, J.B.; Pohll, G.M.; Reeves, D.M. A model-averaging method for assessing groundwater conceptual model uncertainty. Ground Water 2010, 48, 716–728. [Google Scholar] [CrossRef] [PubMed]
  34. Engelhardt, I.; De Aguinaga, J.G.; Mikat, H.; Schüth, C.; Lenz, O.; Liedl, R. Complexity vs. Simplicity: Groundwater Model Ranking Using Information Criteria. Ground Water 2013, 52, 573–583. [Google Scholar] [CrossRef] [PubMed]
  35. Akaike, H. Information Theory and an Extension of the Maximum Likelihood Principle; Breakthroughs in Statistics, 1; Springer-Verlag: New York, NY, USA, 1992; pp. 610–624. [Google Scholar]
  36. Akaike, H. A new look at statistical model identification. IEEE Trans. Autom. Control 1974, 19, 716–723. [Google Scholar] [CrossRef]
  37. Burnham, K.P.; Anderson, D.R. Model Selection and Multimodel Inference: A Practical Information-Theoretic Approach, 2nd ed.; Springer-Verlag: New York, NY, USA, 2002; p. 488. [Google Scholar]
  38. Schwarz, G. Estimating the dimension of a model. Ann. Stat. 1978, 6, 461–464. [Google Scholar] [CrossRef]
  39. Kashyap, R.L. Optimal choice of AR and MA parts in autoregressive moving average models. IEEE Trans. Pattern Anal. Mach. Intell. 1982, 4, 99–104. [Google Scholar] [CrossRef] [PubMed]
  40. Ye, M.; Meyer, P.D.; Neuman, S.P. On model selection criteria in multimodel analysis. Water Resour. Res. 2008, 44. [Google Scholar] [CrossRef]
  41. Centamore, E.; Deiana, G. La Geologia delle Marche. Studi Geol. Camerti 1986, spec. vol., 145. [Google Scholar]
  42. Pierantoni, P.P.; Deiana, G.; Galdenzi, S. Stratigraphic and structural features of the Sibillini Mountains (Umbria-Marche Apennines, Italy). Ital. J. Geosci. 2013, 132, 497–520. [Google Scholar] [CrossRef]
  43. Calamita, F.; Deiana, G. The arcuate shape of the Umbria-Marche Apennines (Central Italy). Tectonophysics 1988, 146, 139–147. [Google Scholar] [CrossRef]
  44. Cantalamessa, G.; Di Celma, C. Sequence response to syndepositional regional uplift: Insights from high-resolution sequence stratigraphy of late early Pleistocene strata, Periadriatic Basin, central Italy. Sediment. Geol. 2004, 164, 283–309. [Google Scholar] [CrossRef]
  45. Thornthwaite, C.W. An Approach toward a Rational Classification of Climate. Am. Geoph. Soc. 1948, 38, 55–94. [Google Scholar] [CrossRef]
  46. Turc, L. The water balance of soils. Relation between precipitation, evaporation and flow. Ann. Agron. 1954, 5, 491–569. [Google Scholar]
  47. Allen, R.G.; Pereira, L.S.; Raes, D.; Smith, M. Crop evapotranspiration: Guide-lines for computing crop water requirements. In FAO Irrigation and Drainage Paper; FAO: Rome, Italy, 1998; p. 56. [Google Scholar]
  48. Hargreaves, G.H.; Samani, Z. Reference Crop Evapotranspiration from Temperature. J. Appl. Eng. Agric. 1985, 1, 96–99. [Google Scholar] [CrossRef]
  49. Samani, Z. Estimating Solar Radiation and Evapotranspiration Using Minimum Climatological Data. J. Irrig. Drain Eng. 2000, 126, 265–267. [Google Scholar] [CrossRef]
  50. Boni, C.F.; Baldoni, T.; Banzato, F.; Cascone, D.; Petitta, M. Studio idrogeologico per l’identificazione, la caratterizzazione e la gestione degli acquiferi del parco nazionale dei monti sibillini. Ital. J. Eng. Geol. Environ. 2010, 2, 21–39. [Google Scholar] [CrossRef]
  51. Perrone, E. Reno, Lamone, Fiumi Uniti, Savio, Marecchia, Foglia, Metauro, Esino, Potenza, Chienti, Tenna, Aso ed Altri Minori: Pozzi Artesiani Della Pianura Emiliana, dal Panaro Alla Marecchia e del Litorale Adriatico, da Porto Corsini a Pesaro; Ministero di Agricoltura, Industria e Commercio: Roma, Italy, 1910; Volume 35, p. 439. [Google Scholar]
  52. Boni, C.F.; Bono, P.; Capelli, G. Schema Idrogeologico dell’Italia Centrale. Mem. Soc. Geol. Ital. 1986, 35, 991–1012. [Google Scholar]
  53. Baldoni, T. Idrogeologia dei bacini Ambro e Tennacola. Bachelor’s Thesis, Università La Sapienza, Roma, Italy, 2002. [Google Scholar]
  54. Celico, P. Prospezioni Idrogeologiche Vol. II, 1st ed.; Liguori Editore: Napoli, Italy, 1988; p. 536. [Google Scholar]
  55. Korkmaz, N. The estimation of groundwater recharge from spring hydrographs. Hydrol. Sci. J. 1990, 35, 209–217. [Google Scholar] [CrossRef]
  56. Schôeller, H. Hydrodynamique dans le karst (écoulementet emmagasinement). Proc. Dubrov. Symp. Hydrol. Fract. Rocks 1967, 1, 3–20. [Google Scholar]
  57. Bonacci, O. Karst springs hydrographs as indicators of karst aquifers. Hydrol. Sci. J. 1993, 38, 51–62. [Google Scholar] [CrossRef]
  58. Diersch, H.J. FEFLOW, Finite Element Modeling of Flow, Mass and Heat Transport in Porous and Fractured Media; Springer-Verlag: Berlin/Heidelberg, Germany, 2014; p. 384. [Google Scholar]
  59. Morris, D.A.; Johnson, A.I. Summary of Hydraulic and Physical Properties of Rock and Soil Materials as Anlayzed by the Hydraulic Laboratory of the U.S. Geological Survey; US Geological Survey: Denver, CO, USA, 1967; pp. D1–D42. [Google Scholar]
  60. Domenico, P.A.; Schwartz, F.W. Physical and Chemical Hydrogeology, 1st ed.; Wiley: New York, NY, USA, 1990; p. 528. [Google Scholar]
  61. Anderson, M.P.; Woessner, W.W. Applied Groundwater Modeling: Simulation of Flow and Advective Transport, 1st ed.; Accademic Press Inc.: San Diego, CA, USA, 1992; p. 381. [Google Scholar]
  62. Caprari, M.; Nanni, T.; Scesi, L. Contributo all’analisi della circolazione idrica, mediante l’applicazione del metodo di Kiraly, in idrostrutture carbonatiche dell’Appennino adriatico. Boll. Soc. Geol. Ital. 2002, 121, 99–120. [Google Scholar]
  63. Doherty, J. PEST—Model-Independent Parameter Estimation. User’s Manual, 5th ed.; Watermark Numerical Computing: Brisbane, Australia, 2010; p. 279. [Google Scholar]
  64. Doherty, J. Addendum to the PEST Manual; Watermark Numerical Computing: Brisbane, Australia, 2014. [Google Scholar]
  65. Doherty, J.; Christensen, S. Use of paired simple and complex models to reduce predictive bias and quantify uncertainty. Water Resour. Res. 2012, 47. [Google Scholar] [CrossRef]
  66. Hill, M.C. Methods and Guidelines for Effective Model Calibration, 1st ed.; US Geological Survey: Denver, CO, USA, 1998. [Google Scholar]
  67. Beven, K.J. A Manifesto for the Equifinality Thesis. J. Hydrol. 2006, 320, 18–36. [Google Scholar] [CrossRef] [Green Version]
Figure 1. (a) Hydrogeological sketch of the study area; and (b) hydrogeological conceptual model for the two springs on study.
Figure 1. (a) Hydrogeological sketch of the study area; and (b) hydrogeological conceptual model for the two springs on study.
Water 08 00271 g001
Figure 2. Hydrogeological complexes of the Umbria-Marche succession.
Figure 2. Hydrogeological complexes of the Umbria-Marche succession.
Water 08 00271 g002
Figure 3. The pluviometric and termometric regime for the studied period at the Pintura di Bolognola station.
Figure 3. The pluviometric and termometric regime for the studied period at the Pintura di Bolognola station.
Water 08 00271 g003
Figure 4. River discharge monitoring (m3/s).
Figure 4. River discharge monitoring (m3/s).
Water 08 00271 g004
Figure 5. Tennacola low group of springs hydrograph compared to the rainfall measured to Pintura di Bolognola station.
Figure 5. Tennacola low group of springs hydrograph compared to the rainfall measured to Pintura di Bolognola station.
Water 08 00271 g005
Figure 6. Evapotranspiration compared with the rainfall data for the period 2004–2014.
Figure 6. Evapotranspiration compared with the rainfall data for the period 2004–2014.
Water 08 00271 g006
Figure 7. Location of the five zones considered.
Figure 7. Location of the five zones considered.
Water 08 00271 g007
Figure 8. Boundary conditions.
Figure 8. Boundary conditions.
Water 08 00271 g008
Figure 9. Simulated versus observed spring hydrograph: (a) model with two zones; (b) model with three zones; (c) model with four zones; and (d) model with five zones.
Figure 9. Simulated versus observed spring hydrograph: (a) model with two zones; (b) model with three zones; (c) model with four zones; and (d) model with five zones.
Water 08 00271 g009
Figure 10. Simulated versus observed spring hydrograph for the years 2006, 2008, 2010 and 2012: (a) model with two zones; (b) model with three zones; (c) model with four zones; and (d) model with five zones.
Figure 10. Simulated versus observed spring hydrograph for the years 2006, 2008, 2010 and 2012: (a) model with two zones; (b) model with three zones; (c) model with four zones; and (d) model with five zones.
Water 08 00271 g010
Figure 11. Regression analysis: (a) model with two zones; (b) model with three zones; (c) model with four zones; and (d) model with five zones.
Figure 11. Regression analysis: (a) model with two zones; (b) model with three zones; (c) model with four zones; and (d) model with five zones.
Water 08 00271 g011
Figure 12. Information criterion and sum of squared weighted residuals assessment of the calibrated models.
Figure 12. Information criterion and sum of squared weighted residuals assessment of the calibrated models.
Water 08 00271 g012
Figure 13. Sensitivity analysis.
Figure 13. Sensitivity analysis.
Water 08 00271 g013
Table 1. Evapotranspiration assessment using different methods and next effective precipitation evaluation. P, rainfall; T, temperature; Etr, real evapotranspiration; Etp, potential evapotranspiration; Et0, reference evapotranspiration; Peff, effective rainfall.
Table 1. Evapotranspiration assessment using different methods and next effective precipitation evaluation. P, rainfall; T, temperature; Etr, real evapotranspiration; Etp, potential evapotranspiration; Et0, reference evapotranspiration; Peff, effective rainfall.
Hydrologic YearP (mm)T (°C)TurcThornthwaiteHargreaves-SamaniHargreaves-Samani CorrectedPenman-Monteith
Etr (mm)Peff (mm)Etp (mm)Peff (mm)ET0 (mm)Peff (mm)ET0 (mm)Peff (mm)ET0,PM (mm)Peff (mm)
20042027.807.91507.551520.25552.551475.25411.531616.27379.661648.14556.481471.32
20051839.007.82501.711337.29519.131319.87379.691459.31348.501490.50522.461316.54
20061401.407.85490.86910.54521.31880.09400.341001.06368.901032.50588.04813.36
2007823.809.80489.18334.62610.29213.51516.76307.04476.75347.05769.3454.46
20081411.208.30504.34906.86527.35883.85424.12987.08392.981018.22633.10778.1
20091814.408.26514.931299.47540.551273.45369.991444.41341.481472.92567.631246.77
20101314.007.77484.96829.44496.12818.28407.39907.01375.90938.50595.47718.93
20111266.608.18494.10772.50522.15743.65464.16802.44428.42838.18718.42548.18
20121342.408.62510.63831.77573.85770.75523.60818.80483.35859.05752.39590.01
20131331.608.13495.93835.67519.31812.09353.48978.12326.271005.33567.77763.83
20142058.008.74534.841523.16533.801523.40368.471689.53339.841718.16588.361469.64
Mean value1511.878.31502.641009.23537.86974.02419.961091.92387.461124.41623.59888.29
Table 2. Effective infiltration assessment using different methods for the period between 2004 and 2014.
Table 2. Effective infiltration assessment using different methods for the period between 2004 and 2014.
Hydrologic YearParameterTurcThornthwaiteHargreaves-SamaniHargreaves-Samani CorrectedPenman-Monteith
2004P (mm)2027.82027.82027.82027.82027.8
Et (mm)507.6552.6411.5379.7556.5
Peff (mm)1520.31475.31616.31648.11471.3
R (mm)76.073.880.882.473.6
Ieff (mm)1444.21401.51535.51565.71397.8
2005P (mm)1839.01839.01839.01839.01839.0
Et (mm)501.7519.1379.7348.5522.5
Peff (mm)1337.31319.91459.31490.51316.5
R (mm)66.966.073.074.565.8
Ieff (mm)1270.41253.91386.41416.01250.7
2006P (mm)1401.41401.41401.41401.41401.4
Et (mm)490.9521.3400.3368.9588.0
Peff (mm)910.5880.11001.11032.5813.4
R (mm)45.544.050.051.640.7
Ieff (mm)865.0836.1951.0980.9772.7
2007P (mm)823.8823.8823.8823.8823.8
Et (mm)489.2610.3516.8476.8769.3
Peff (mm)334.6213.5307.0347.154.5
R (mm)16.710.715.417.42.7
Ieff (mm)317.9202.8291.7329.751.7
2008P (mm)1411.21411.21411.21411.21411.2
Et (mm)504.3527.4424.1393.0633.1
Peff (mm)906.9883.9987.11018.2778.1
R (mm)45.344.249.450.938.9
Ieff (mm)861.5839.7937.7967.3739.2
2009P (mm)1814.41814.41814.41814.41814.4
Et (mm)514.9540.6370.0341.5567.6
Peff (mm)1299.51273.91444.41472.91246.8
R (mm)65.063.772.273.762.3
Ieff (mm)1234.51210.21372.21399.31184.4
2010P (mm)1314.41314.41314.41314.41314.4
Et (mm)485.0496.1407.4375.9595.5
Peff (mm)829.4818.3907.0938.5718.9
R (mm)41.540.945.447.036.0
Ieff (mm)788.0777.4861.7892.0683.0
2011P (mm)1266.61266.61266.61266.61266.6
Et (mm)494.1522.2464.2428.4718.4
Peff (mm)772.5744.5802.4838.2548.2
R (mm)38.637.240.141.927.4
Ieff (mm)733.9707.2762.3796.3520.8
2012P (mm)1342.41342.41342.41342.41342.4
Et (mm)510.6573.9523.6483.4752.4
Peff (mm)831.8768.6818.8859.1590.0
R (mm)41.638.440.943.029.5
Ieff (mm)790.2730.1777.9816.1560.5
2013P (mm)1331.61331.61331.61331.61331.6
Et (mm)495.9519.3353.5326.3567.8
Peff (mm)835.7812.3978.11005.3763.8
R (mm)41.840.648.950.338.2
Ieff (mm)793.9771.7929.2955.1725.6
2014P (mm)2058.02058.02058.02058.02058.0
Et (mm)534.9533.8368.5339.9588.4
Peff (mm)1523.21524.21689.51718.21469.6
R (mm)76.276.284.585.973.5
Ieff (mm)1447.01448.01605.11632.31396.2
AverageP (mm)1511.91511.91511.91511.91511.9
Et (mm)502.6537.9420.0387.5623.6
Peff (mm)1009.2974.01091.91124.4888.3
R (mm)50.548.754.656.244.4
Ieff (mm)958.8925.31037.31068.2843.9
Table 3. Tennacola low group of springs budget analysis. Q0, discharge at the beginning of the period measurement; Qt, discharge at the time t; V0, dynamic reserve at the beginning of the period; Vt, dynamic reserve at the end of the period; ΔV, dynamic reserve change during the water year; Qya, groundwater discharge volume during the water year; R, groundwater recharge.
Table 3. Tennacola low group of springs budget analysis. Q0, discharge at the beginning of the period measurement; Qt, discharge at the time t; V0, dynamic reserve at the beginning of the period; Vt, dynamic reserve at the end of the period; ΔV, dynamic reserve change during the water year; Qya, groundwater discharge volume during the water year; R, groundwater recharge.
Hydrologic YearInitial DischargesDynamic ReserveDynamic Reserve ChangeDischargeRecharge
Q0 (m3/s)Qt (m3/s)V0 (m3)Vt (m3)ΔV (m3)Qya (m3)R (m3)
20040.0814050.1672311,004,770.532,064,108.341,059,337.85,224,3516,283,689
20050.1636780.0999152,020,254.571,233,234.84−787,019.725,832,7175,045,698
20060.0946870.0957031,168,703.871,181,246.2412,542.3675,484,0575,496,599
20070.0785920.063809970,045.879787,591−182,454.882,407,6262,225,171
20080.0735670.087314908,021.5411,077,701.27169,679.734,746,1554,915,835
20090.0738310.115963911,288.6951,431,319.19520,030.495,599,0066,119,037
20100.1171080.0977391,445,448.31,206,378.76−239,069.545,530,2125,291,143
20110.0975660.0962141,204,238.861,187,551.71−16,687.1485,501,5485,484,861
20120.0797990.144573984,948.7431,784,445.76799,497.023,794,0604,593,557
20130.1538520.0952751,898,976.221,175,969.27−723,006.955,600,4134,877,406
20140.0862820.103671,064,962.681,279,584.59214,621.924,968,1845,182,806
Total value1.1003661.16740613,581,659.914,409,131827,471.0954,688,33055,515,801
Average value0.0921890.0967031,137,8821,193,595.1255,713.564,971,6665,046,891
Table 4. The range of hydrogeological parameters for the five zones used in the analysis. Kxx,yy,zz = hydraulic conductivity for the different space orientations; Sy = specific yield; TR IN/OUT = transfer rate input/output; minimum (min), reference and maximum (max) values of hydrogeological parameters.
Table 4. The range of hydrogeological parameters for the five zones used in the analysis. Kxx,yy,zz = hydraulic conductivity for the different space orientations; Sy = specific yield; TR IN/OUT = transfer rate input/output; minimum (min), reference and maximum (max) values of hydrogeological parameters.
ModelZoneSector IDKxx (m/s)Kyy (m/s)Kzz (m/s)Sy (-)TR IN/OUT (1/day)
MinRefMaxMinRefMaxMinRefMaxMinRefMaxMinRefMax
M2Z1sBas1e−081e−061e−051e−081e−061e−051e−091e−071e−060.0080.010.051e−018.64e−041e+10
2sSca,sMai,sFuc,sTen1e−045e−041e−031e−045e−041e−031e−055e−051e−040.0020.0050.01
M3Z1sBas1e−081e−051e−031e−081e−051e−031e−081e−051e−030.0010.0050.1
2sSca, sMai,sFuc1e−091e−071e−051e−091e−071e−051e−091e−071e−050.0050.050.5
3sTen1e−031e−021e−011e−031e−021e−011e−031e−021e−010.050.10.5
M4Z1sBas1e−093.5e−071e−051e−095.8e−081e−051e−091.6e−071e−050.0050.0090.5
2sSca-SMai1e−084e−051e−031e−081.6e−051e−031e−082.5e−051e−030.0010.0050.5
3sTen1e−034.5e−031e−011e−031.7e−021e−011e−031.2e−021e−010.050.30.5
4sFuc1e−081e−071e−051e−081e−071e−051e−081e−071e−050.0010.0020.1
M5Z1sBas1e−093.5e−071e−051e−095.8e−081e−051e−091.6e−071e−050.0050.0090.5
2sSca1e−061e−051e−031e−061e−051e−031e−061e−051e−030.010.050.5
3sTen1e−034.5e−031e−011e−031.7e−021e−011e−031.2e−021e−010.050.30.5
4sFuc1e−081e−071e−051e−081e−071e−051e−081e−071e−050.0010.0050.1
5sMai1e−061e−041e−031e−061e−041e−031e−061e−041e−030.010.050.5
Table 5. Statistical index derived from PEST analysis. MEw, Average value of the weighted residuals; MAEw, Absolute average value of the weighted residuals; σr2, variance of weight residuals; σw, standard deviation; Φ (SSWR), sum of squared of the weighted residuals; R, correlation index.
Table 5. Statistical index derived from PEST analysis. MEw, Average value of the weighted residuals; MAEw, Absolute average value of the weighted residuals; σr2, variance of weight residuals; σw, standard deviation; Φ (SSWR), sum of squared of the weighted residuals; R, correlation index.
ModelMEwMAEwσr2σΦR
M2Z−0.3000.3900.140.371330.73
M3Z−0.3710.5090.400.602250.72
M4Z−0.5560.7640.590.805170.73
M5Z−0.5500.7400.550.744910.72
Table 6. Difference between simulated and observed total discharge for the single year considered (ΔQ).
Table 6. Difference between simulated and observed total discharge for the single year considered (ΔQ).
Model2006200820102012
ΔQ (m3)ΔQ (m3)ΔQ (m3)ΔQ (m3)
M2Z1,119,3291,428,2741,489,0131,343,655
M3Z1,035,0461,455,9531,534,6581,403,962
M4Z1,054,7141,284,5871,538,1681,258,872
M5Z978,5681,353,9701,313,5711,239,845
Average1,046,9141,380,6961,468,8531,311,584
Table 7. Difference of the Akaike information criterion (AIC), Bayesian information criterion (BIC), and Kashyap’s information criterion (KIC) values and likelihood of the flow models from the Akaike Weights (AIC wi).
Table 7. Difference of the Akaike information criterion (AIC), Bayesian information criterion (BIC), and Kashyap’s information criterion (KIC) values and likelihood of the flow models from the Akaike Weights (AIC wi).
ModelAICKICBICAIC wi
M2Z0001.00
M3Z310.85328.26357.710.00
M4Z797.66710.75817.640.00
M5Z775.75827.97706.830.00

Share and Cite

MDPI and ACS Style

Giacopetti, M.; Crestaz, E.; Materazzi, M.; Pambianchi, G.; Posavec, K. A Multi-Model Approach Using Statistical Index and Information Criteria to Evaluate the Adequacy of the Model Geometry in a Fissured Carbonate Aquifer (Italy). Water 2016, 8, 271. https://doi.org/10.3390/w8070271

AMA Style

Giacopetti M, Crestaz E, Materazzi M, Pambianchi G, Posavec K. A Multi-Model Approach Using Statistical Index and Information Criteria to Evaluate the Adequacy of the Model Geometry in a Fissured Carbonate Aquifer (Italy). Water. 2016; 8(7):271. https://doi.org/10.3390/w8070271

Chicago/Turabian Style

Giacopetti, Marco, Ezio Crestaz, Marco Materazzi, Gilberto Pambianchi, and Kristijan Posavec. 2016. "A Multi-Model Approach Using Statistical Index and Information Criteria to Evaluate the Adequacy of the Model Geometry in a Fissured Carbonate Aquifer (Italy)" Water 8, no. 7: 271. https://doi.org/10.3390/w8070271

APA Style

Giacopetti, M., Crestaz, E., Materazzi, M., Pambianchi, G., & Posavec, K. (2016). A Multi-Model Approach Using Statistical Index and Information Criteria to Evaluate the Adequacy of the Model Geometry in a Fissured Carbonate Aquifer (Italy). Water, 8(7), 271. https://doi.org/10.3390/w8070271

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