Next Article in Journal
Impact of Hybrid Water Supply on the Centralised Water System
Previous Article in Journal
Removal of Heavy Metals from Urban Stormwater Runoff Using Bioretention Media Mix
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Statistical and Fractal Approaches on Long Time-Series to Surface-Water/Groundwater Relationship Assessment: A Central Italy Alluvial Plain Case Study

1
Department of Engineering and Geology “InGeo”, “G. D’Annunzio” University of Chieti-Pescara, 66013 Chieti, Italy
2
National Research Council (CNR), Research Institute for the Hydrogeological Protection (IRPI), 70126 Bari, Italy
*
Author to whom correspondence should be addressed.
Water 2017, 9(11), 850; https://doi.org/10.3390/w9110850
Submission received: 6 October 2017 / Revised: 25 October 2017 / Accepted: 26 October 2017 / Published: 3 November 2017

Abstract

:
In this research, univariate and bivariate statistical methods were applied to rainfall, river and piezometric level datasets belonging to 24-year time series (1986–2009). These methods, which often are used to understand the effects of precipitation on rivers and karstic springs discharge, have been used to assess piezometric level response to rainfall and river level fluctuations in a porous aquifer. A rain gauge, a river level gauge and three wells, located in Central Italy along the lower Pescara River valley in correspondence of its important alluvial aquifer, provided the data. Statistical analysis has been used within a known hydrogeological framework, which has been refined by mean of a photo-interpretation and a GPS survey. Water–groundwater relationships were identified following the autocorrelation and cross-correlation analyses. Spectral analysis and mono-fractal features of time series were assessed to provide information on multi-year variability, data distributions, their fractal dimension and the distribution return time within the historical time series. The statistical–mathematical results were interpreted through fieldwork that identified distinct groundwater flowpaths within the aquifer and enabled the implementation of a conceptual model, improving the knowledge on water resources management tools.

1. Introduction

The careful and correct management of water resources in densely populated alluvial areas is becoming necessary because of the groundwater depletion caused by increased exploitation associated with demographic growth. Previous scientific research on long time series was focused mainly on hydrogeological processes concerning relationships between precipitation and spring discharge or river level, quite often considering fractured aquifer systems [1,2,3,4,5,6]; these analyses provided complementary information to the hydrogeological properties of a natural system, to the transit time through the subsurface’s unsaturated zone, and to the recharge of the aquifer, supporting the correct implementation of hydrogeological conceptual and numerical models [7,8,9,10]. On the other hand, the behavior analysis of long piezometric level time series was almost neglected and limited to short periods that rarely exceeded a year (e.g., contaminated site studies). Recently, Cai and Ofterdinger [11] evaluated the statistical relationships between two-year rainfall and piezometric time series logged consecutively every 15-min in fractured aquifers, while Misstear et al. [12] deepened the relationships between well hydrographs and rainfall recharge in five-year periods in sandy-gravel aquifers. Furthermore, the definition of relationships between piezometric level fluctuation and its regulating factors is useful to characterize alluvial aquifers and better quantify and manage groundwater resources, especially in dry periods when piezometric heads decrease and, similarly, in wet periods when their increase can contribute to groundwater flooding.
Rainfall should be considered an input or independent regulating variable for the aquifer and for the hydrographic network, as it is a key factor affecting respectively infiltration recharge and runoff, while the river level and the groundwater piezometric level should be considered as outputs or dependent variables. At the same time, the piezometric level cannot be considered to be independent from the river stage, because of surface-water/groundwater interactions [3,6,13,14,15,16,17]. In this study, 24-year time series of rainfall, river level and groundwater piezometric level measurements, acquired every three days, were analyzed with various univariate and bivariate statistical methods, including autocorrelation, cross-correlation, and spectral analyses. In addition, the residual component of the data has been elaborated to deepen the random relationships between the variables. The analyses were undertaken in order to evaluate memory effects, the delay of piezometric level response to rainfall [18] and river head changes, and the periodical components of the time series. The methods developed by Mangin [1] can describe the capability of the unsaturated and saturated zones of the aquifer in filtering, amplifying or attenuating and, accordingly, transforming the input functions (rainfall ajd river level) into the output functions (piezometric level).
In addition to the standard statistical analyses, the feasibility of using fractal tools [19,20] to characterize the variability of rainfall, hydrometric data and piezometric level of selected wells was explored. These methods for the analysis of hydrogeological data are unusual in the literature; Yu et al. [21] used the fractal analysis to relate piezometry data only to the geological profile and not to other hydrological time series. In this research, a mono-fractal analysis, based on the Box-Counting and Cantor-Dust methods, was carried out [22,23,24], and then some indications on events return time (DRT) were given.
The methodology is applied to an alluvial aquifer located in the Italian Central Adriatic area (Figure 1). The paper is subdivided in four sections: in the “Study Area” section, the territory is described from the geological, morphological and hydrogeological point of view. In the sections “Datasets” and “Methods” the time series and the statistical–mathematical analysis methods are illustrated. In Section 4, “Results”, the most significant results are shown. The final sections “Discussion” and “Conclusions” offer a wide result comparative analysis in which relationships with the hydrological and hydrogeological contests are described and criticisms and future research development are outlined.

2. Study Area

The studied alluvial aquifers are located between the Apennine Mountains and the Adriatic Sea. The more elevated areas in this sector are constituted of Mio-Pliocenic terrigenous pelitic-arenaceous and pelitic-marly deposits [25,26,27]. The alluvial aquifers are sedimentary structures that include various orders of terraces (I to IV order) and the current alluvial plain. These deposits are mostly silty and sandy in nature with gravelly-sandy lenses of the period between medium-upper Pleistocene and Holocene age. The spatio-temporal evolution of the fluvial depositional environments, driven by interactions between climatic and glacio-eustatic changes and the plio-quartenary tectonic phases, has conditioned the spatial distribution of the characteristic lithotypes [28]. The presence of a mostly pelitic substrate, acting as an aquiclude (k values ranging from 10−9 to 10−7 m/s), below the more permeable alluvial deposits (k values ranging from 10−5 to 10−3 m/s) suggests the presence of shallow, single and multi-layered aquifers [29,30].
These aquifers often interact with the streams and rivers through the hyphorheic zone. Hydraulic connection between the stream and the groundwater can occur with a transfer of mass and/or pressure head, based on the local hydrogeological conditions [31]. The alluvial aquifers are fed by the limestone aquifers of the Apennine Mountains in the internal areas [32,33,34,35,36], while recharge from rainfall has a major role in the total recharge to the shallow alluvial aquifers moving downstream towards the sea, particularly in correspondence of the terraced deposits that are hydraulically connected to the lower valleys [37,38,39,40]. The wells considered in this study fall in the terminal areas of the Pescara River alluvial valley, near the Adriactic Sea.
The Pescara River has a basin of 3190 km2 (Figure 1), which can be subdivided in two main parts: the first one is the western mountainous sector, and the second one is the eastern sector, where the alluvial aquifer under study is located. The Pescara River mean discharge measured in the Pescara Santa Teresa gauge (PST) is about 50 m3/s (31 m3/s comes from the western part of its basin).
The alluvial deposits making up the aquifer present four orders of terraces, which have different elevations with respect to the base of the valley. The first and second order terraces are at higher elevations and are mainly constituted of gravelly deposits with silty-sandy lenses; the third and four order terraces are at lower elevations and are constituted of gravelly deposits with extensive silty-sandy and silty-clayey lenses. Only the fourth order terrace runs on both sides of the river, with thin gravelly-sandy deposits, considered to be the current alluvial deposits, being directly connected to this terrace (Figure 2). The overall thickness of the quaternary alluvial deposits increases moving downstream. Thickness increases from 12 to 16 m in the middle section of the Pescara River valley to 40–50 m near the coast. As the thickness increases, the nature of the deposits varies from mainly gravelly to silty-sandy.
The groundwater flow direction follows that of the river in the most upstream areas of the valley, although groundwater flow will be dependent on the local spatial distribution of the lithotypes (presence of paleostreams, etc.) and on the morphology of the substrate, which will be connected to erosional history and tectonic deformations. The groundwater flow is toward the river moving downstream where recharge to the aquifer is clearly originating from the terraces (Figure 1). The morphology of the substrate has a strong role in controlling groundwater flow, since groundwater divides give rise to secondary flow circuits within the aquifer. The average hydraulic gradient of the shallow groundwater below the river bed varies from 4‰ to 6‰ [29,41]. A meandering hydrographic pattern is observed in the area, which, together with the changes in local river characteristics (i.e., river head slope and river width variations), controls the relationship between groundwater and surface water (Figure 1 and Figure 3).
The geometric characteristics of the river section under investigation, assessed by means of aerial photo-interpretation and in-situ GPS surveys, indicate that the width of the river is quite variable. However, in general (Figure 3), the width increases moving downstream. There are two discontinuities along the river profile, which correspond to river dams for hydroelectric purposes, and the river head slope is very low upward of these artificial impoundments. These hydroelectric plants were constructed after the period in which the data was collected and therefore had no influence on the collected datasets.
The selected wells are located on the right respect to the river bank, at variable distances from the river bed (Figure 4, Table 1). These wells have been drilled to total depth through silty-sandy deposits and their exploitation is sporadic and limited to backyard irrigation during summer season.

3. Datasets

Datasets used in this study, made available by the Hydrographic Service of Abruzzo Region, refer to “Spoltore” automatic rain gauge (Sp), “Pescara S. Teresa” hydrometer (PST) and three wells (Figure 1 and Figure 2, and Table 1), “Surricchio” (Su), “Sanità” (Sa) and “De Nicola” (DN). While rainfall data were available on a daily basis, the piezometric and hydrometric levels instead were measured every three days, which means that, for every month, 10 measures were made available. For this reason, to compare the time series, daily rainfall datasets where cumulated on a three-day basis. Since the data available are relative to the 24-year period 1986–2009, each time series consisted of 2880 records. Well characteristics and stratigraphy were identified from previous geognostic and field surveys (Table 1, and Figure 4). The sea level time series has not been considered because of the distance of each monitoring point from the sea (Table 1).
The three-day time series are plotted against time in Figure 5. River level and rainfall data monitored at the Pescara S. Teresa (PST) and Spoltore (Sp) sites, and piezometric heads measured from the Sa and DN wells show that the variability appears to be mainly seasonal, while at first sight a plurennial cyclic behavior is evident in the piezometric level of Su well.

4. Methods

4.1. Analyses in Time Domain

With reference to the statistical theory [1,2,3,4,42,43,44,45,46], univariate autocorrelation and bivariate cross-correlation functions were applied to the time series as well as their residual component, obtained with the additive Seasonal-Trend decomposition procedure Census I.

4.1.1. Seasonal-Trend Decomposition Procedure Census I

This statistical method was applied to raw data to remove the periodical calendar related fluctuations (seasonal and cyclical) and the trend component, representing the non-casual long term changes in the level of the time series. In this way, it was possible to extract from the original data the casual residual component describing the random short term fluctuations, neither systematic nor predictable [47,48,49].
The selected statistical approach is based on the assumption that the observed value x t at time t of the time series under consideration (or, in the case of a discrete ensemble of elements, the t-th data measured at the time t, thus 1 t N ) can be thought as consisting of four different components, combined to each other in an essentially additive model:
x t = T t × C t + S t + ε t
where T t is the trend component; S t is the seasonal component that occurs at regular (seasonal) intervals; C t is the cyclical component that varies from cycle to cycle and that have usually a longer duration than the seasonal component; and ε t is the random, error, or irregular component. In the Census I method, the trend and cyclical components are customarily combined into a trend-cycle component T t × C t [47,48].

4.1.2. The Autocorrelation Function

The Autocorrelation Function (ACF) provides the opportunity to evaluate the linear dependency of successive values of a single parameter belonging to a time series. For this reason, the method is univariate and quantifies the “memory effect” that corresponds to the temporal reciprocal influence on subsequent data. Therefore, the ACF correlogram describes the time necessary to “forget” its own initial conditions [1,3,4,44,50,51].
In general, a correlogram with a gentle slope can be an indicator of data series persistence, while, if it drops off rapidly, it describes the random nature of the values [45]; furthermore, the slower decline of ACF function for hydrological data describes an aquifer with low draining properties, low permeability or major groundwater storage, conversely a faster decline of the autocorrelation function indicates a more rapid flow of water through the aquifer and/or its limited storage capacity [1,2,5,11,18,52,53].
To compare the “inertia” of signals, a delay (lag) is identified in correspondence of the decorrelation coefficient 0.2, since events with a lower ACF coefficient can be considered quasi-independent [1,6,54].

4.1.3. The Cross-Correlation Function

The Cross-Correlation Function (CCF), bivariate analysis performed between rainfall and river level input signals and between rainfall and piezometric output signals, describes the time lag necessary for the latter to reach the maximum CCF value in response of the input signals. Furthermore, the function’s slope characterizes an aquifer in terms of infiltration rate, in terms of draining capacity and storage and, by doing so, it describes the system’s modulating power and consequently its inertial nature [1,2,3,13,42].
The cross-correlation function, to be acceptable, has to be characterized by a significant correlation at the 95% confidence level, which means a correlation coefficient superior to the standard error ~ 2 N = 0.037, where N is 2880 and corresponds to the number of values in the selected datasets [55,56].
Applying this function to the raw data, we obtain information on both the casual and non-casual components relationships of input and output variables, as well as their importance [3]. On the other hand, cross-correlation applied only to residuals describes the relationship between the random components of the data, giving an idea about the response of the system against a unitary impulse, hence the cross-correlogram can be used as the “unitary hydrograph” in surface hydrology [2,3]. Furthermore, the delay of the response gives an estimation of the pressure pulse transfer times (the piezometric level increase due to hydrostatic pressure increase in the aquifer) and of the particle travel times through the aquifer [5,6,9,44,48,49,50,54,57]. In our study, because of the small surface area considered, the input rainfall signal, referred to SP rain gauge, is common to all wells and consequently the variability of the transfer function can be related only to site-specific structural and functional variability of the hydrological processes [57].

4.2. Analyses in Frequency Domain

The frequency domain is useful for data manipulation to extract frequency signals and provide a picture of the frequency contents that are unclear in the time domain data [58]. Furthermore, this analysis allows detecting average periodicities in long term series, to see if there are any dominants oscillation modes and to identify the frequency bands of interest [42,59,60].
Assuming that the time series under study { x t } is a linear combination of sinusoidal functions, each representing the multiple periodical nature of the phenomena, the discrete Fourier’s analysis, also called harmonic analysis, can be employed. This univariate analysis decomposes the time series with cyclical components, into few underlying co-sinusoidal and sinusoidal functions of all possible (discrete) frequencies, each characterized by an amplitude related with the variance of the time series. Accordingly, when N is odd:
{ x t = x ¯ + n = 1 q 2 N { a n cos ( n 2 π N t ) + b n sin ( n 2 π N t ) } + ε t a 0 = x ¯ = ( t = 1 N x t ) / N a n = t = 1 N x t cos ( n 2 π N t ) b n = t = 1 N x t sin ( n 2 π N t )
where x ¯ is the mean of the time series, ν n = n / N is the n-th frequency (harmonics) multiple of the fundamental frequency 1/N; and ε t is the random component associated to the t-th data and q = ( N 1 ) / 2 . If N is even, q = N / 2 , a q = ( t = 1 N ( 1 ) t x t ) / N and b q = 0 .
It worth noting that the cosine and the sine parameters, respectively, a n and b n in Equation (2), are obtained applying the least squares estimation method and indicate the degree to which the respective cosine and sine functions are correlated with data. To explore how much each cyclic phenomenon affects the overall time series trend, the amplitude square of each component, function of the related frequency (spectral Analysis), was identified as an “importance measure” parameter.
The Fourier Transform, usually employed through its fast computational version (FFT; a particular gain is obtained if the number of data is equal to N = 2 n 1 × 3 n 2 × 5 n 3 × , in this paper N = 2880 = 2 6 × 3 2 × 5 , is a common approach aimed at performing univariate spectral analysis, providing the feasibility to switch from the temporal to the frequency domain. Accordingly, the following expression: I ( ν n ) = ( a n 2 + b n 2 ) 2 / N , 0 n N / 2 ( 0 υ n 1 / 2 ), called “periodogram” or the “sample spectrum”, was introduced. After some straightforward manipulations, a link between the sample spectrum and the estimate of the autocovariance, called also “power spectrum” [42,44] emerges. Moreover, in order to compare time series obtained at different scales of measurements, it is useful to introduce the estimated (for this reason characterized by a circumflex accent) “spectral density function” S D ^ F ( ν ) , normalizing the power spectrum respect to the variance σ ^ x 2 . Finally, the following expression was exploited:
S D ^ F ( ν ) = 2 { c ^ x x ( 0 ) σ ^ x 2 + 2 k = 1 N 1 c ^ x x ( k ) σ ^ x 2 cos ( 2 π ν k ) }   0 ν 1 2
where c ^ x x ( k ) is the estimate of the k-th lag auto-covariance and c ^ x x ( 0 ) σ ^ x 2 is the estimate of the variance (for more details see [42,44]). That is, S D ^ F ( ν ) could be interpreted as the Fourier transform of the autocorrelation A C ^ F ( k ) of the time series. The physical meaning of the S D ^ F ( ν ) is that the area under a portion of the total plot, related to a selected frequency range, represents the contribution to variance of components with frequencies within the range. When the entire spectrum is drawn, the equation indicates that the area underneath the curve is equal to the variance of the process. However, the simple Fourier analysis is based on the assumption of fixed amplitudes, frequencies and phases, while time series are affected by random changes of these kind of parameters.
A consequence is that Expression (3) fluctuates uncontrollably and, for this reason, is not useful for any meaningful interpretations. Accordingly, the S D ^ F ( ν ) expressed by Equation (3) is not a consistent estimator [43] for an ensemble of variables taken from a discrete purely random process (de-trended time series). A way to solve this kind of issue is to smooth the periodogram, exploiting a sort of truncation of the autocorrelation function, giving less weight to the values of the A C ^ F ( k ) as k increases. As discussed in [42,43] , the following truncated S D ^ F ( ν ) T appears to be a consistent estimator to consider properly the randomness nature of the time series under study:
S D ^ F ( ν ) T = 2 { λ ( 0 ) c ^ x x ( 0 ) σ ^ x 2 + 2 k = 1 M λ ( k ) c ^ x x ( k ) σ ^ x 2 cos ( 2 π ν k ) }   0 ν 1 2
where { λ ( k ) } is an ensemble of weights called the “lag window” and M < N is called the “truncation point” [43]. In this paper, we adopted, among many others, the Blackman–Tukey window expression:
λ ( k ) = 1 2 { 1 + cos ( k π M ) } 0 k M
The value of M, following the detailed discussion reported in [16,43], was selected equal to 5. Hereafter, to simplify the symbology, ACF, CCF and SDF stand for A C ^ F ( k ) , C C ^ F ( k ) , and S D ^ F ( ν ) T .
Summarizing, the power spectrum of a time series is the simple translation of its autocovariance function (i.e., the ratio of autocorrelation and variance) in the frequency; complementary to correlation analysis, it may reveal certain characteristics present in a time series that are buried in noise and not discernable otherwise [61,62].

4.3. Fractal Analysis

The data clustering of time series, expressed by only the number of events occurring within a period, cannot fully describe the possible large time irregularities of events [24]. For this reason, a preliminary fractal approach application was proposed in this study, in addition to the time series statistical approach described above. The purpose of the fractal analysis was to gain more information in relation to the variability of the data, and therefore provide insights into this issue.
One of the most important properties of a fractal object, time series included, is that it can be subdivided into parts (or intervals of time), each of which is (usually only approximately) a reduced-size copy of the whole. This approximate property is termed “self-affine” and implies that properties such as frequencies of selected events are almost invariant across spatial and temporal scales. Their irregular or complex occurrences, in general, can be “measured” to some extent, by the distance of the Hausdorff–Besicovitch dimension D, called hereafter the “fractal dimension”, from an integer number [19,20,63,64,65]. However, not all time series data are truly fractal [66]. Therefore, a preliminary analysis of the assessed data is recommended in order to verify the occurrence of fractal behavior. A self-affine performance of a time series could be detected, for example, by long-term correlations resulting from the analysis of the ACF and by the evidence of a power-law relation between the power spectrum S(f) and the time frequency of the occurrence of a selected event: S ( f ) f β , resulting from the Spectral Analysis [66]. The exponent β is linked to the Hurst’s coefficient H and to the fractal dimension D [67,68,69].
The fractal features of different hydrological time series, including rainfall and river dynamics, were highlighted in previous papers [19,20,22,23,63,65,66,67,68,69,70,71,72]. To have a fast response, we carried out a fast fractal analyses as a first step. Nevertheless, as we demonstrated in the following sections, some interesting results were obtained. The feasibility of the fractal analysis application to the time series studied in this paper emerge from the inspection of the results obtained by both Autocorrelation Function and Spectral Analysis. Different types of fractal dimension and the methods of estimating these have been proposed in literature [65,73].
In this paper, an algorithm based on the Box-Counting and Cantor-Dust methods, successfully applied in [22,23,70,71], was developed. As a first attempt, only raw data were considered. Subsequently, the values of the data belonging to each time series, i.e. rainfall intensities, river hydrometry and three wells piezometric levels, were normalized with respect to the full range of their variability within their own time series. The resulting percentages (accordingly from 0 to 100%, to be considered as fluctuations above the minimum values) were considered as the “events” that can be compared to each other, although originating from different sources. In the following the related analytical expression is given:
f i , j ( % ) = x i , j x min , j x max , j x min , j × 100
where f i , j ( % ) (“fluctuation”) is the “event” related to the i-th data x i , j belonging to the j-th time series (j = 1,2,3,4,5, since the time series under consideration are 5); and x min , j and x max , j are, respectively, the minimum and the maximum numerical values of the j-th time series data. According to the algorithm procedures, different values of f i , j ( % ) , hereafter “F_threshold”, should be introduced. No particular constraints are required, thus, in this paper 7–9 values between 0% up to 90% for both the river hydrometry and wells, while just only up to 40% for Spoltore rainfall were selected. Then the time range was divided into increasing integer numbers (1, 2, 3,…) of equal-time interval groups, Ng. Each interval group contains the same integer number ε of measured data. Ng values that did not satisfy the N g × ε = 2880 requirement was excluded. Therefore, the resulting Ng numbers are all the integer divisors of 2880. For each accepted number Ng of intervals, the number of interval time N ( ε ) , in which at least one event above the selected F_threshold occurred, was counted. We introduced the parameter “n_events” as the number of the fluctuations f i , j ( % ) that, in the selected interval time, satisfy a selected comparison relationship (≤ or ≥) with a selected F_threshold. To be more clear, if over the interval time under consideration (i.e., from 9-th May 1990 to 10-th June 1990) the rainfall fluctuations f i , j ( % ) , measured at Spoltore station, were equal to or above the selected F_threshold (i.e., 60%, thus F_threshold = 0.6) for 5 times, accordingly it results that n_events = 5 and the value of N ( ε ) is increased by just only 1. The following power law relation [22,23,65] was assumed:
N ( ε ) = m ( 1 ε ) D
where m is a constant (equal to 1 in this paper) and D is a parameter that is identified as the fractal dimension if:
D = lim ε 0 [ ln N ( ε ) ln m ] / ln ( 1 / ε )
ln N ( ε ) was plotted versus ln ( 1 / ε ) and D was estimated from the slope of the best fit straight line, obtained by means of the Minimum Square Root Method. Generally, the D parameter is not constant and depends on the selected F_threshold values. Accordingly, the value of the D parameter calculated in this paper should be intended, in general, as an average fractal D dimension. Moreover, the comparison of the fractal analyses to the results obtained in the previous statistical approaches should be considered as qualitative, rather than quantitative, although some quantitative considerations are also made.
To better clarify the physical meaning of the variability of the fractal dimension D by the mean of the Cantor Dust approach (DF in Figures 13 and 15), the following example can be useful. Two rainfall time series are analyzed. The data were obtained by daily measurements over 30 years (10,920 values), at two gauges, A and B, representative of two different sites (A and B). We are interested to study how the occurrences of daily rainfall above or equal to 50 mm (F_threshold) are distributed over the selected interval of time. If we consider a time scale of one year, the entire interval is subdivided into 30 subintervals, each lasting one year. Accordingly, Ng = 30 and ε = 365, as resulting from the application of the Cantor Dust method. Suppose that for 40 times (n-events = 40) the intensity of the rainfall, measured at A and B stations exceeded the value of the selected F_threshold over, respectively, 2 and 22 different years, even non-continuously. Consequently, N(ε) = 2 and N(ε) = 22, respectively, for time series A and B.
It could be reasonably inferred that, within the site where Station A is located, rainfall values exceeding the considered F_threshold are rather uncommon events over the selected time scale of 30 years. Their occurrences are concentrated during just two years. On the other hand, within the site where Station B is located, this kind of events is quite common and the expected return time would be reasonable equal to about one year. Then, consider a further partition of the whole time range into 60 subintervals (Ng = 60 and ε = 182). Moreover, assume that, as the time scale decreases, N(ε) = 2 remains the same for Station A, while N(ε) = 37, increases (thus almost 1 for each semester) for Station B. Repeating this operation by the means of the Cantor Dust approach, the resulting value of the D fractal dimension related to the time series A would be low (in this case close to 0), while for the other time series, D would be close to 1.
As highlighted by the previous example and according to Mazzarella [23,24], we can infer that, for the time series considered in this paper, the D Cantor Dust fractal dimension can range from 0, corresponding to randomly distributed events, just as isolated points to the lim ε 0 D = 0 , to 1 ( lim ε 0 D = 1 : a line by a geometrical point of view) corresponding to a distribution of the selected events occurrence that remains almost the same within all the considered time scales.

5. Results

5.1. Descriptive Statistics

The variability of the data, described by the standard deviation analysis (Table 2), points out that the piezometric level fluctuations within the Su and Sa wells are higher than in the DN well. Their average water depth (unsaturated zone thickness) is calculated from the difference between the altitude and the average piezometric level for the 1986–2009 period.
Averaging for each month the 10 available piezometry, hydrometry and cumulated rainfall data, over the 24-year period of 1986–2009, it is possible to represent the monthly average of each variable (Figure 6). River level at the PST site presents the maximum value in April and the minimum value in July, while the Sp rain gauge presents a maximum in December and a minimum in July. Furthermore, the standardization of the computed piezometric levels through the normalization, based on the mean and standard deviation, gives us the chance to compare at the same scale the entire 1986–2009 group of data, showing that all wells present the lowest piezometric level in August and the highest in January, February and March, respectively, for the DN, Su and Sa wells.

5.2. The Residual Analysis

The comparison at a three-day resolution between the 2880 piezometric residuals data for the whole 1986–2009 period identifies in Su site the records with the greatest piezometric variability; furthermore, a similar residual distribution especially for the Su well, the Sp rain gauge and the PST river level can be identified, since the peaks fall together on several dates (Figure 7, left side).
The periods during which the casual components are concurrently more influent can be identified by calculating the average annual series on the basis of a monthly arrangement of the data, which consisted in averaging the 10 monthly three-day resolution records, over the 24 years of 1986–2009 period (Figure 7, right side); it is evident that the months of March–April, October–November and June–July, in correspondence of the spring, autumn and summer periods, respectively, present the greatest random variability.
Some exceptions can be highlighted: for example, a greater residual concentration during the 2007 summer is self-evident at the DN location and is characterized by a strong variability of the raw data. In December 1994, an intense rainfall event was observed at the Spoltore rain gauge and a correspondent peak was identified in the rainfall residual. However, it was not identified in the datasets of river stage and piezometric head residuals, probably because of the local character of this event.

5.3. Autocorrelation

The correlogram relative to rainfall observed at the Sp rain gauge (Figure 8A) shows how, after a three-day time lag, the ACF function immediately decreases to uncorrelated values. This confirms its random behavior because events with duration longer than three days are rare and, therefore, the correlation between the data for intervals longer than the reference lag is nil.
The river level reaches the decorrelation value threshold ACF = 0.2 after 81 days, highlighting a strong memory effect (Figure 8B). However, the two factors influencing the ACF function shape are the runoff velocity and the catchment area; indeed, the curve shape in the first three days, when the ACF value decreases to 0.7, is characterized by a greater slope and a more rapid decorrelation rate down, meaning a stronger independence of the data. In fact, in this short period, the river level ACF curve has the same shape of the rainfall one because it describes the flooding action caused by direct runoff, while after three days the more autocorrelated shape of the river level ACF is controlled by the baseflow due to the large recharge area and groundwater contribution [74].
The hydrogeological time series as the piezometric levels of our study are strongly autocorrelated in all sites because of inertia and carryover processes in the physical system, that impart a correlation to successive piezometry levels. In detail, Su reaches 0.2 threshold value after 880 days, while the ACF curves for Sa and DN reach that level after respectively 105 and 90 days (Figure 8). Su well stronger autocorrelation might be related to site-specific physical and structural conditions of the aquifer, such as the hydraulic conductivity influenced by changes in the grain size of the alluvial aquifer or its degree of compaction [18,75], because, of course, as suggested by Imagawa et al. [53], persistence can be an indicator of the local storage capacity of alluvial aquifers, provided that the assumption of homogeneity in the hydraulic conductivity throughout the aquifer is respected.

5.4. Raw Data Cross-Correlation

The piezometric level fluctuation in answer to rainfall (Figure 9) highlights low correlations value for all the wells. The Su and DN wells behave the same, since the shape of the curves are similar, with a maximum CCF value around 0.16, reached at both sites after, respectively, a time lag of nine and six days.
The cross-correlation analysis between groundwater levels and river level highlights a major dependency of groundwater level fluctuations to the river stage, rather than to rainfall. In fact, the maximum CCF values vary from 0.20 in Su, 0.25 in DN and 0.37 in Sa with corresponding time lags of three days (Su, DN) and six days in Sa (Figure 9). The DN well, which is 1960 m from Pescara River, presents the intermediate maximum cross-correlation value.

5.5. Residual Cross-Correlation

The CCF between rainfall and river level residual components shows an initial high positive correlation value (0.34), indicating that an immediate positive fluctuation of river stage corresponds to the three-day precipitation impulse. After three days, this correlation becomes negative, highlighting that river level continues to rise positively for three days after the end of the rainfall event (Figure 10). It is important to remember that the river level is affected by both the local precipitation events (monitored at the Sp rain gauge) and by those events happened upward in the river drainage basin, far from the selected rain gauge.
The rainfall–piezometric and river–piezometric levels residual cross-correlation are similarly shaped, even though the river impulse affects the well fluctuations more than rainfall, especially for the Su and Sa wells, which are closer to the river. In fact, in Su well, the residuals CCF precipitation–piezometric level and river level–piezometric level are, respectively, 0.28 and 0.36 (Figure 11); these maximum values are reached immediately as a consequence of the first impulse event, while after six days the correlation becomes negative (Figure 11) pointing out that the piezometric level continues to rise for six days after the rainfall–river stage pulse.
At Sa well, the maximum CCF value is lower (0.05–0.1), while DN well level presents slight correlations to rainfall and river stage (0.06–0.05).

5.6. Spectral Analysis

The power spectrum of precipitation, containing homogeneously unimodal distributed amplitudes, is dominated only by an annual cycle representing the wet/dry seasonal cycle and is void of other prevailing periodic component. Hydrometric and piezometric levels spectrum, instead, are not homogeneously distributed, and present a bimodal behavior, because in addition to the seasonal cycle present high amplitude oscillations at lower-frequencies (higher periods) concentrated particularly between five and eight years, with broader peaks more evident in correspondence of the interannual 72-month cycle (≈6 years) and of the interdecadal 144-month cycle (≈12 years), especially in Su and Sa wells located closer to the river (Figure 12).
Hydrometric and piezometric level response to the rainfall main driver results from the passage of infiltrated water through several subsurface filters that, depending on the scale range, have the effect to either amplify or attenuate the precipitation variability modes [76,77]. Accordingly, the hydrometric and piezometric level output signals can have significantly different spectral characteristics than the input precipitation signal; in our study comparing their power spectrum in a log–log scale graph it is evident how strongly is related the frequency of piezometry and hydrometry, while this relationship is much weaker with precipitation frequency. Consequently, the unimodal spectrum of the latter is not affected by a power law type, whilst the previous are characterized by a bimodal power spectra shifted towards low frequency scales (Figure 12, in red).
Spectral density and amplitude in Figure 12 indicate that piezometric trends seem to follow approximately a decaying exponential function S ( f ) f β = T β (T is the period and f is the frequency). They are localized in two intervals for Su and Sa, 12 months < T < (72~144) months and 0 months < T < 12 months, and only in the latter interval for DN. These occurrences further support the applicability of fractal analysis, as discussed above. Moreover, according to many authors [78,79,80] due to the variation of the trends observed in Figure 12 (related to the piezometry Spectral densities, before and after one year), it could be argued that the time series of the piezometric levels and just weakly for the hydrometry, show a multifractal behaviors.

5.7. Cantor-Dust Fractal Analysis

The characteristic of the selected F_threshold events and their occurrences are detailed in Table 3 for each time series. More details regarding the meaning of parameters are discussed in the text, in particular in the Discussion section (Section 6.3. Fractal analysis).
In Figure 13, for Sp rainfall and PST hydrometry respectively, the highest values of the average fractal dimension, ranging from D = 1.000 to D = 0.851 and D = 0.906 were reached for level values up to 5% and 10%, corresponding to ΔL = 100 mm and ΔL = 0.52 m fluctuations (Table 3).
The outcomes of piezometric level variations are plotted in Figure 14. Irregular distributions of points occur for all the three wells at high values of F_threshold events with normalized level variations higher than 80–90%. This kind of drawback, due to the sensitivity of the method to small differences in the box size, can be solved using the Shifted Box Counting Method [22]. For the DN time series (Table 3), normalized events, related to values ranging from 1% to 20% (Figure 14) of the total level variability, show a similar variation (average fractal dimensions from D = 1.000 to D = 0.9986) of the distribution level values, corresponding to ΔL = 0.27 m. For the Su and Sa wells, the highest values of the average fractal dimension, ranging from D = 1.000 to D = 0.963 and D = 0.988, respectively, were reached for level values up to 20%, corresponding to ΔL = 1.4 m and ΔL = 0.58 m level fluctuations (Table 3).
It is roughly possible to identify several “breaking points” in Figure 13 and Figure 14 and the underlying data on which the figures are based. Such breaking time scales are specific for each F_threshold value and for each time series source, for which the relation between the number of covering boxes and the interval time (box size), shows a departure from a linear trend. The latter has a fractal dimension that is very close to unity. As discussed above, if the fractal dimension of a dataset is close to unity for a specific time scale, the data show similar distributions within time scales higher than the identified one. Accordingly, this time-scale could be considered for the introduction of a “distribution return time” (DRT) of the selected event, which should not be confused with the commonly used “event return time”, because the last parameter is calculated by means of different approaches. Therefore, DRT values (estimated through inspection of Figure 13 and Figure 14 and the underlying data) of different f i , j ( % ) values are reported in Figure 15.
It is interesting to note that the Sp Rainfall and the DN estimated distribution curve of DRT show a similar sigmoidal regression trend of the following type:
T = a / { 1 + exp [ ( N v p b ) / c ] }
where T is the DRT (days); Nvp is the Normalized variability percentage (that is f i , j ( % ) ); and “a”, “b”, and “c” are three parameters, the values of which are different rainfall and DN data. The regression of DRT related to the other sources follows exponential trends:
T = a + b exp ( c N v p )
where “a”, “b”, and “c” are the three regression parameters.

6. Discussion

The relationship between rainfall, river level and piezometric levels was assessed by means of statistical–mathematical analyses and a conceptual model of the local hydrogeology.

6.1. Analyses in Time Domain

The DN well is the furthest from the main stream, therefore recharge effects from the Pescara River are considered to be negligible, while rainfall infiltration represents the main local recharge; this hypothesis is confirmed by the CCF analysis of the rainfall and piezometry data (Figure 9A), which shows that DN reaches the highest CCF value. An additional confirmation is provided by the CCF of the residuals, which indicates the absence of a relationship between the DN well and the river (Figure 11B) meaning that extreme rainfall events appear to have less influence on the piezometry, as suggested by its low residual component. The higher CCF value between the well and the river level (Figure 9B) can be explained by the simultaneous changes of piezometric and river levels, due to the seasonal cycle, even if the two oscillations are not physically correlated. This is confirmed by the hydrogeology of the area. In fact, the water table depth in the DN well is quite low (Figure 4; Table 2). Similarly, the difference between the head in the DN well (mean 10.46 m a.s.l.) and the elevation of the Pescara River (mean 3.33 m a.s.l.) is significant. In addition, the groundwater hydrodynamic indicates the presence of a groundwater divide, which has a direction that is approximately parallel to the direction of the river (Figure 1). This excludes that the water level changes in the well can be induced by the river.
Conversely, groundwater and surface water interaction is likely in the area near the Su and Sa wells, since they are located near the Pescara River (Figure 16).
The Su well is located near the river (Figure 4, Table 1) and has piezometric heads that are comparable to the elevations of the river section that falls just upstream of the well itself. The piezometric level spatial distribution indicates that groundwater flows from the river to the aquifer, thus recharging the latter (Figure 16). The analysis of the raw data CCF (Figure 9A) indicates that there is a low correlation between rainfall and piezometric level fluctuations. This is in agreement with both the high average depth of the water table (6.70 m; Table 2) and with the idea of aquifer recharge from the river. The value of the river/well raw data CCF (Figure 9B) is generally lower compared to the corresponding residual component CCF (Figure 11B). A possible explanation is the position of the gauging station compared to the well location; the analysis of the hydraulic conditions (Figure 3) indicates that the river section in correspondence of the gauging station is wider compared to the section just upstream of the well. Therefore, considering the same discharge, a wider section should imply smaller river level changes. The higher cross-correlation of the residuals confirms this. In fact, when high discharges occur during flood events, the higher river level recorded downstream of the well is indicative of even higher changes upstream where the river section is narrower. However, the strong impulsive changes in the piezometric level within the well, clearly represented by the strong oscillation of the residual component (Figure 11B), is due to a pressure transfer rather than a mass transfer, given the low hydraulic conductivity of the aquifer (about k = 10–5 m/s).
The Sa well is also located near the Pescara River (Figure 4) and has piezometric head that is constantly higher than the river level upstream of the well itself (Figure 16). This indicates a constant discharge from the aquifer to the river. The piezometric level supports this idea and the flow lines are constantly towards the river. The cross-correlation with the rainfall shows a very low value (Figure 9A) that, together with the high water depth (mean 3.89 m; Table 2), allows assessing the low contribution of the direct recharge from local rainfall. The CCF with the river level (Figure 9B) shows a relatively high value. This apparent anomaly can be explained with the presence of a significant annual cycle in both datasets, based on the hydrogeological understanding of the area. A further confirmation is provided by the CCF function of the residuals, which shows very low values for both the rainfall and the river level fluctuation (Figure 11).

6.2. Analyses in Frequency Domain

The spectral analysis has highlighted a clear seasonal cycle (one year) in piezometry, hydrometry and rainfall spectral densities distribution. Furthermore, evident multi-year periodicities have been detected in the piezometric and river levels fluctuations (pointed out in Figure 12). The rainfall spectrum is characterized by a unimodal and homogeneous power spectrum distribution and, even though multi-year cycles are present, they are not dominant compared to the annual cycle. On the other hand, piezometric and river levels present a clearly bimodal and not homogeneous spectral density distribution, characterized by a sharp spectral power decrease with frequency increase (Figure 12).
This effect is well known in the literature and is due to low-frequency pass filter effect of the aquifer and of the river basin, on the piezometric and river levels, respectively [53,59,76,78,81]. Thus, these two measures modulate the rainfall input signal, lowering the short period fluctuations and highlighting the long period ones (seasonal and multi-year) that are in association to large scale climatic oscillations, as many researchers have reported [59,79,82,83,84,85,86,87].
The water volumes that are at stake proportionally influence the weight of this filtering effect [16,75,76,77], therefore extended aquifers and river basins can be considered as excellent climate signal proxies, at regional scale [18,88,89]. Similar results are described by Luque-Espinar [16] and Fendeková et al. [90], who found that piezometers, located close to stream with leakage to the aquifer, or along the main river in those areas where it recharges the aquifer, present long period cycles. In addition, Chiaudani et al. [91] found out the presence of multi-year and decadal cycles in Abruzzo Region in 11 wells belonging to alluvial aquifers located near to their main stream.

6.3. Fractal Analysis

It was possible to estimate a sort of “pseudo” fractal dimension through the fractal analysis discussed in this paper. This fractal dimension can be regarded as average values of trends, within the appropriate log–log domain for this kind of investigation. As expected, the rainfall plot (left of Figure 13) shows that these types of events are more “isolated” and random than “events” associated to river level (right of Figure 13) and piezometric levels (Figure 14). This can be explained by the fact that the river level and piezometric variations are more continuous in time. They are the result of a much more complex phenomena related to not only rainfall at one site, but also to the hydrological dynamics of the entire basin.
The river level plot shows that almost the same event distribution is related to variations of levels lower than 10% (<0.52 m). This distribution of events, which corresponds to F_threshold = 3.26 m level (see Table 3) occurs within all time scales. It is interesting to note that the qualitative trends, shown in Figure 14, are in agreement with the trends obtained for the Warta River flows in [22], even if the phenomena dynamics are different. The occurrence of the same fractal dimension equal to 1 means that each distribution of such level fluctuations remains similar in all the selected time interval scales; however, distributions related to different sources of data are not necessary equal. Consequently, it would be hard to identify any correlation between different time series, where only the value of the fractal dimension obtained by the Cantor Dust approach is observed.
Figure 15 and Table 3 show that the DRT distributions of both rainfall and the DN piezometric data are well interpolated by sigmoidal curves, while Su, Sa and river stage are better interpolated by exponential ones. This result is in agreement with the statistical analysis that identifies two distinct groups of data with respect to well defined hydrological relationships. The shift of the sigmoidal curve of DN well respect to the rainfall one indicates that the piezometric levels shows a random behavior for higher L_thresholds, as the breaking points in Figure 13 and Figure 15 and the previous discussion highlight. This could be related to the modulator effect of the aquifer which gives higher stationarity to DN time series. Moreover, the interpolating sigmoidal and exponential expressions could be compared to the classical return-time laws proposed in literature ([92], and references therein).

6.4. Conceptual Model

All the applied statistical–mathematical analyses, considering hydrogeological knowledge, allowed defining a final conceptual model of the system hydrodynamic (Figure 17).
In particular, the DN well is located in a Pescara aquifer part where a secondary flow path has been recognized. Here, the recharge is completely related with local rainfall, described by the Sp rain gauge, and with the runoff infiltration at the clayey slope base.
On the other hand, the main flow path has been identified near the Pescara River. In this area, there are clear surface-water/groundwater relationships, considering either an aquifer recharge from the main stream or the opposite. In detail, upstream, where the Su well is located, the main river recharges the aquifer. Furthermore, when huge flooding events occur the Su well detects outright the increase of the hydrometric level, probably due to a pressure transfer. Contrariwise, downstream near the Sa well, the piezometric head is constantly higher than river level, even during huge flooding events. Thus, it is clear that the aquifer discharges groundwater into the Pescara River.
Finally, the local rainfall monitored in Sp rain gauge, near the Pescara River, does not seem to influence substantially the piezometric fluctuation, due to the importance of aquifer baseflow.

7. Conclusions

The assessment of the relationships between groundwater and both rainfall and river recharge influencing the alluvial aquifer hydrodynamic needs a detailed hydrogeological framework definition. Long time series are useful for this purpose, and, in this study, hydrological data series, monitored for the 24-year period 1986–2009, have been analyzed.
The statistical–mathematical analysis of rainfall, river level and groundwater level time series data has produced useful information for the understanding of their relationships existing within the studied hydrogeological context. To better understand the hydraulic conditions and to make statistical results coherent, aerial photo-interpretation analyses and in-situ GPS surveys have been realized.
The “memory effect”, representing the self-coherency indicator of each time series, was assessed with the Auto-Correlation Function and with the fractal analysis pointing out how stationary behaviors were higher for piezometric and river levels than for rainfall. Interdependency between different hydrological parameters was found by means of the Cross-Correlation Function, highlighting a strong groundwater/surface water interaction monitored between the Pescara River and the nearest wells. On the other hand, the farthest well showed a clear correlation with local rainfall. This occurrence could also be observed by fractal analysis from the regression trends of the “distribution return times”, for high time scales.
Moreover, the Cross-Correlation Function applied to residual data, after seasonal cycle and trend removal, points out a strong pressure transfer from the river to groundwater during flooding events, when the river leakage recharges the aquifer.
The spectral analysis identified in all data series a predominant annual cycle, linked to seasonal fluctuations, while multi-year climatic related regimen are detected in the hydrometry and in piezometry affected respectively from the low-pass filter action of the river basin and the aquifer. These results are supported also from autocorrelation stationary analysis and the fractal approach that shows a more systematic behavior (high D values) both in hydrometry and in piezometry than in rainfall (low D values).
These research results highlighted that statistical–mathematical analyses performed on long time series, when coupled with detailed assessments of the local hydrogeological conditions, can provide a quantitative understanding of the spatial-temporal relationships that exist between input and output parameters. In this respect, the application of such tools can be also useful for numerical modeling and in all those cases (e.g., water resource management, groundwater vulnerability, etc.), where it is important to understand the type and the entity of the relationships between surface-water and groundwater.

Acknowledgments

The authors would like to thank the “Servizio Idrografico e Mareografico of Abruzzo Region” for having provided the piezometric, rainfall and river level data. We are also thankful to the learned referees and the editorial board for the important suggestions, which helped in improving the manuscript.

Author Contributions

All the authors were involved in all phases of the article development, setup and revision.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Mangin, A. Pour une meilleure connaissance des systemes hydrologiques à partir des analyses corrélatoire et spectrale. J. Hydrol. 1984, 67, 25–43. [Google Scholar] [CrossRef]
  2. Padilla, A.; Pulido-Bosch, A. Study of hydrographs of karstic aquifers by means of correlation and cross-spectral analysis. J. Hydrol. 1995, 168, 73–89. [Google Scholar] [CrossRef]
  3. Larocque, M.; Mangin, A.; Razack, M.; Banton, O. Contribution of correlation and spectral analyses to the regional study of a karst aquifer (Charente, France). J. Hydrol. 1998, 205, 217–231. [Google Scholar] [CrossRef]
  4. Thanh Tam, V.; De Smedt, F.; Batelaan, O.; Dassargues, A. Characterization of a cavern conduit system in Vietnam by time series correlation, cross-spectrum and wavelet analyses. Hydrol. Sci. J. 2004, 49, 879–900. [Google Scholar] [CrossRef]
  5. Panagopoulos, G.; Lambrakis, N. The contribution of time series analysis to the study of the hydrodynamic characteristics of the karst systems: Application on two typical karst aquifers of Greece (Trifilia, Almyros Crete). J. Hydrol. 2006, 329, 368–376. [Google Scholar] [CrossRef]
  6. Delbart, C.; Valdes, D.; Barbecot, F.; Tognelli, A.; Richon, P.; Couchoux, L. Temporal variability of karst aquifer response time established by the sliding-windows cross-correlation method. J. Hydrol. 2013, 511, 580–588. [Google Scholar] [CrossRef]
  7. Pasculli, A.; Audisio, C.; Lollino, G. Application of cellular automaton model for river morphological studies: CAESAR and the Pellice River (Piedmont, Italy). Rend. Online Soc. Geol. Ital. 2010, 11, 118–119. [Google Scholar]
  8. Pasculli, A.; Audisio, C.; Sciarra, N. Application of CAESAR for catchment and river evolution. Rend. Online Soc. Geol. Ital. 2015, 35, 224–227. [Google Scholar] [CrossRef]
  9. Audisio, C.; Pasculli, A.; Sciarra, N. Conceptual and numerical models applied on the river Pellice (North western Italy). In River Basins, Reservoir Sedimentation and Water Resources Engineering. Proceedings of the IAEG XII Congress; Springer International Publishing: Torino, Italy, 2015; Volume 3, pp. 327–330. [Google Scholar] [CrossRef]
  10. Pasculli, A.; Audisio, C. Cellular Automata modelling of fluvial evolution: Real and parametric numerical results comparison along river pellice (NW Italy). Environ. Model. Assess. 2015, 20, 425–441. [Google Scholar] [CrossRef]
  11. Cai, Z.; Ofterdinger, U. Analysis of groundwater-level response to rainfall and recharge estimates in fractured hard rock aquifers, NW Ireland. J. Hydrol. 2016, 535, 71–84. [Google Scholar] [CrossRef] [Green Version]
  12. Misstear, B.D.R.; Brown, L.; Johnston, P.M. Estimation of groundwater recharge in a major sand and gravel aquifer in Ireland using multiple approaches. Hydrogeol. J. 2009, 17, 693–706. [Google Scholar] [CrossRef]
  13. Polemio, M. Il Regime Della Falda Costiera Ionica di Metaponto. In Proceedings of the Third Conference of Geotechnical Engineering Researchers, “Il Ruolo dei Fluidi nei Problemi di Ingegneria Geotecnica”, Mondovì (Cuneo-Italy), Italy, 6–7 September 1994; pp. 135–149. Available online: www.earth-prints.org/bitstream/2122/6025/1/falda-metaponto.pdf (accessed on 15 October 2017).
  14. Lee, J.Y.; Lee, K.K. Use of hydrologic time series data for identification of recharge mechanism in a fractured bedrock aquifer system. J. Hydrol. 2000, 229, 190–201. [Google Scholar] [CrossRef]
  15. Polemio, M.; Dragone, V. Gli effetti antropici e climatici sulla disponibilità delle acque sotterranee della piana ionico-lucana. In Proceedings of the First National Congress—Italian Association of Applied and Environmental Geology, Chieti, Italy, 19–20 February 2003. [Google Scholar]
  16. Luque-Espinar, J.A.; Chica-Olmo, M.; Pardo-Igùzquiza, E.; Garcìa-Soldado, M.J. Influence of climatological cycles on hydraulic heads across a Spanish aquifer. J. Hydrol. 2008, 354, 35–52. [Google Scholar] [CrossRef]
  17. Cozzolino, M.; Piscopo, V.; Scarelli, A.; Summa, G. Analisi delle serie temporali dei dati di precipitazione e portata per la determinazione della portata sostenibile di un pozzo funzionante a carico costante. Ital. J. Eng. Geol. Environ. 2010, 2, 41–50. [Google Scholar] [CrossRef]
  18. Duvert, C.; Jourde, H.; Raiber, M.; Cox, M.E. Correlation and spectral analyses to assess the response of a shallow aquifer to low and high frequency rainfall fluctuations. J. Hydrol. 2010, 527, 894–907. [Google Scholar] [CrossRef]
  19. Mandelbrot, B.B. The Fractal Geometry of Nature; Macmillan: New York, NY, USA, 1983. [Google Scholar]
  20. Turcotte, D.L. Fractal and Chaos in Geology and Geophysics; Cambridge University Press: New York, NY, USA, 1997. [Google Scholar]
  21. Yu, X.; Ghasemizadeh, R.; Padilla, I.Y.; Kaeli, D.; Alshawabkeh, A. Patterns of temporal scaling of groundwater level fluctuation. J. Hydrol. 2016, 536, 485–495. [Google Scholar] [CrossRef]
  22. Radziejewiski, M.; Kundzewicz, Z. Fractal analysis of flow of the river Warta. J. Hydrol. 1997, 200, 280–294. [Google Scholar] [CrossRef]
  23. Mazzarella, A. The time clustering of floodings in Venice and the Cantor dust method. Theor. Appl. Climatol. 1998, 59, 147–150. [Google Scholar] [CrossRef]
  24. Mazzarella, A. Multifractal dynamic rainfall processes in Italy. Theor. Appl. Climatol. 1999, 63, 73–78. [Google Scholar] [CrossRef]
  25. Casnedi, R. L’avanfossa abruzzese fra i fiumi Vomano e Pescara nel Pliocene inf.: Rapporti fra sedimentazione e tettonica. Studi Geologici Camerti. 1991, 357–379. [Google Scholar]
  26. Centamore, E.; Bigi, S.; Berti, D.; Micarelli, A.; Morelli, C. Nuovi dati sui depositi neogenici di avanfossa del Pescarese. Boll. Soc. Geol. Ital. 1992, 111, 437–447. [Google Scholar]
  27. D’Alessandro, L.; Miccadei, E.; Piacentini, T. Morphostructural elements of central–eastern Abruzzi: Contributions to the study of the role of tectonics on the morphogenesis of the Apennine chain. Quat. Int. 2003, 101–102, 115–124. [Google Scholar] [CrossRef]
  28. Parlagreco, L.; Miccadei, E.; Antonioli, F.; Gianolla, D.; Devoti, D.; Leoni, G.; Silenzi, S. New data on Holocene relative sea level along the Abruzzo coast (central Adriatic, Italy). Quat. Int. 2011, 232, 179–186. [Google Scholar] [CrossRef]
  29. Desiderio, G.; Nanni, T.; Rusi, S. La pianura alluvionale del fiume Pescara (Abruzzo): Idrogeologia e vulnerabilità dell’acquifero. Mem. Soc. Geol. Ital. 2001, 56, 197–211. [Google Scholar]
  30. Desiderio, G.; Nanni, T.; Rusi, S. La pianura del fiume Vomano (Abruzzo): Idrogeologia, antropizzazione e suoi effetti sul depauperamento della falda. Boll. Soc. Geol. Ital. 2003, 122, 421–434. [Google Scholar]
  31. Anibas, C.; Verbeiren, B.; Buis, K.; Chormanski, J.; De Doncker, L.; Okruszko, T.; Meire, P.; Batelaan, O. A hierarchical approach on groundwater-surface water interaction in wetlands along the upper Biebrza River, Poland. Hydrol. Earth Syst. Sci. 2012, 16, 2329–2346. [Google Scholar] [CrossRef]
  32. Boni, C.; Bono, P.; Capelli, G. Schema idrogeologico dell’Italia centrale. Mem. Soc. Geol. Ital. 1986, 35, 991–1012. [Google Scholar]
  33. Conese, M.; Nanni, T.; Peila, C.; Rusi, S.; Salvati, R. Idrogeologia della Montagna del Morrone (Appennino Abruzzese): Dati preliminari. Mem. Soc. Geol. Ital. 2001, 56, 181–196. [Google Scholar]
  34. Massoli-Novelli, R.; Petitta, M.; Salvati, R. La situazione idrogeologica e ambientale della Gole di Popoli (Abruzzo): Primi risultati e prospettive di ricerca. Mem. Soc. Geol. Ital. 1998, 53, 563–584. [Google Scholar]
  35. Nanni, T.; Rusi, S. Idrogeologia del massiccio carbonatico della montagna della Majella (Appennino centrale). Boll. Soc. Geol. Ital. 2003, 122, 173–202. [Google Scholar]
  36. Petitta, M.; Tallini, M. Idrodinamica sotterranea del massiccio del Gran Sasso (Abruzzo): Nuove indagini idrologiche, idrogeologiche e idrochimiche (1994–2001). Boll. Soc. Geol. Ital. 2001, 121, 343–363. [Google Scholar]
  37. Palmucci, W.; Rusi, S.; Di Curzio, D. Mobilisation processes responsible for iron and manganese contamination of groundwater in Central Adriatic Italy. Environ. Sci. Pollut. Res. 2016, 23, 11790–11805. [Google Scholar] [CrossRef] [PubMed]
  38. Palmucci, W.; Rusi, S. Boron-rich groundwater in Central Eastern Italy: A hydrogeochemical and statistical approach to define origin and distribution. Environ. Earth Sci. 2015, 72, 5139–5157. [Google Scholar] [CrossRef]
  39. Desiderio, G.; Folchi Vici D’Arcevia, C.; Nanni, T.; Rusi, S. Hydrogeological mapping of the highly anthropogenically influenced Peligna Valley intramontane basin (Central Italy). J. Maps. 2012, 8, 165–168. [Google Scholar] [CrossRef]
  40. Caschetto, M.; Barbieri, M.; Galassi, D.M.P.; Mastrorillo, L.; Rusi, S.; Stoch, F.; Di Cioccio, A.; Petitta, M. Human alteration of groundwater–surface water interactions (Sagittario River, Central Italy): Implication for flow regime, contaminant fate and invertebrate response. Environ. Earth Sci. 2014, 71, 1791–1807. [Google Scholar] [CrossRef]
  41. Desiderio, G.; Rusi, S.; Tatangelo, F. Caratterizzazione idrogeochimica delle acque sotterranee abruzzesi e relative anomalie. Ital. J. Geosci. 2010, 129, 207–222. [Google Scholar] [CrossRef]
  42. Jenkins, G.M.; Watts, D.G. Spectral Analysis and Its Applications; Holden Day: San Francisco, CA, USA, 1968. [Google Scholar]
  43. Chatfield, C. The Analysis of Time Series, 4th ed.; Chapman and Hall: London, UK, 1991. [Google Scholar]
  44. Box, G.E.P.; Jenkins, G.M.; Reinsel, G.C. Time Series Analysis: Forecasting and Control, 3rd ed.; Prentice Hall: Englewood Cliffs, NJ, USA, 1994. [Google Scholar]
  45. Massei, N.; Dupont, J.P.; Mahler, B.J.; Laignel, B.; Fournier, M.; Valdes, D.; Ogier, S. Investigating transport properties and turbidity dynamics of a karst aquifer using correlation, spectral, and wavelet analysis. J. Hydrol. 2006, 329, 244–257. [Google Scholar] [CrossRef]
  46. Fiorillo, F.; Doglioni, A. The relation between karst spring discharge and rainfall by cross-correlation analysis (Campania, southern Italy). Hydrogeol. J. 2010, 18, 1881. [Google Scholar] [CrossRef]
  47. Makridakis, S.G.; Wheelwright, S.C.; McGee, V.E. Forecasting: Methods and Applications, 2nd ed.; Wiley: New York, NY, USA, 1983. [Google Scholar]
  48. Makridakis, S.G.; Wheelwright, S.C. Forecasting Methods for Management, 5th ed.; Wiley: New York, NY, USA, 1989. [Google Scholar]
  49. Polemio, M.; Dragone, V. Serie storiche piezometriche delle unità idrogeologiche pugliesi: Regime piezometrico, effetti climatici ed antropici. Quad. Geol. Appl. 1999, 4, 153–162. [Google Scholar]
  50. Benavente, J.; Pulido-Bosch, A.; Mangin, A. Application of correlation and spectral procedures to the study of the discharge in a karstic system (Eastern Spain). In Proceedings of the Karst Water Resources Symposium, Ankara-Antalya, Turkey, 7–19 July 1985; IAHS: Wallingford, UK, 1985; Volume 161, pp. 67–75. [Google Scholar]
  51. Polemio, M.; Dragone, V.; Di Cagno, M. Effetti antropici e naturali sul degrado quantitativo delle acque sotterranee del Tavoliere. Quad. Geol. Appl. 1999, 4, 143–152. [Google Scholar]
  52. Bouchaou, L.; Mangin, A.; Chauve, P. Turbidity mechanism of water from a karstic spring: Example of the Ain Asserdoune spring (Beni Mellal Atlas, Morocco). J. Hydrol. 2002, 265, 34–42. [Google Scholar] [CrossRef]
  53. Imagawa, C.; Takeuchi, J.; Kawachi, T.; Chono, S.; Ishida, K. Statistical analyses and modeling approaches to hydrodynamic characteristics in alluvial aquifer. Hydrol. Process. 2013, 27, 4017–4027. [Google Scholar] [CrossRef]
  54. Amraoui, F.; Razack, M.; Bouchaou, L. Turbidity dynamics in karstic systems. Example of Ribaa and Bittit springs in the Middle Atlas (Morocco). Hydrol. Sci. J. 2003, 48, 971–984. [Google Scholar] [CrossRef]
  55. Diggle, P.D. Time Series: A Biostatistical Introduction; Oxford Statistical Science Series 5; Oxford Science Publications-Oxford University Press: Oxford, UK, 1990. [Google Scholar]
  56. Lee, L.J.E.; Lawrence, D.S.L.; Price, M. Analysis of water-level response to rainfall and implications for recharge pathways in the Chalk aquifer, SE England. J. Hydrol. 2006, 330, 604–620. [Google Scholar] [CrossRef]
  57. Delbart, C.; Valdes, D.; Barbecot, F.; Tognelli, A.; Couchoux, L. Spatial organization of the impulse response in a karst aquifer. J. Hydrol. 2016, 537, 18–26. [Google Scholar] [CrossRef]
  58. Acworth, R.I.; Rau, G.C.; McCallum, M.; Anderson, M.S.; Cuthbert, M.O. Understanding connected surface-water/groundwater systems using Fourier analysis of daily and sub-daily head fluctuations. Hydrogeol. J. 2015, 23, 143–159. [Google Scholar] [CrossRef]
  59. Andreo, B.; Jiménez, P.; Durán, J.J.; Carrasco, F.; Vadillo, I.; Mangin, A. Climatic and hydrological variations during the last 117–166 years in the south of the Iberian Peninsula, from spectral and correlation analyses and continuous wavelet analyses. J. Hydrol. 2005, 324, 24–39. [Google Scholar] [CrossRef]
  60. Gill, L.W.; Naughton, O.; Johnston, P.M.; Basu, B.; Ghosh, B. Characterisation of hydrogeological connections in a lowland karst network using time series analysis of water levels in ephemeral groundwater-fed lakes (turloughs). J. Hydrol. 2013, 499, 289–302. [Google Scholar] [CrossRef]
  61. Box, G.E.P.; Jenkins, G.M.; Reinsel, G.C. Time Series Analysis: Forecasting and Control, 4th ed.; Wiley: Hoboken, NJ, USA, 2008. [Google Scholar]
  62. Percival, D.B.; Walden, A.T. Spectral Analysis for Physical Applications: Multitaper and Conventional Univariate Techniques; Cambridge University Press: Cambridge, NY, USA, 1993. [Google Scholar]
  63. Rodriguez-Iturbe, I.; Rinaldo, A. Fractal River Basins: Chance and Self-Organization; Cambridge University Press: New York, NY, USA, 2001. [Google Scholar]
  64. Falconer, K.J. Fractal Geometry: Mathematical Foundations and Applications; John Wiley & Sons: New York, NY, USA, 1990. [Google Scholar]
  65. Zhou, X. Fractal and Multifractal Analysis of Runoff Time Series and Stream Networks in Agricultural Watersheds. Ph.D. Dissertation, Faculty of the Virginia Polytechnic Institute and State University, Blacksburg, VA, USA, October 2004. [Google Scholar]
  66. Kantelhardt, J.W. Fractal and multifractal time series. In Encyclopedia of Complexity and Systems Science; Meyers, R.A., Ed.; Springer: New York, NY, USA, 2009; pp. 3754–3779. [Google Scholar]
  67. Svensson, C.; Olsson, J.; Berntsson, R. Multifractal properties of daily rainfall in two different climates. Water Resour. 1996, 32, 2463–2472. [Google Scholar] [CrossRef]
  68. Menabde, M.; Harris, D.; Seed, A.; Austin, G.; Stow, D. Multiscaling properties of rainfall and bounded random cascade. Water Resour. 1997, 33, 2823–2830. [Google Scholar] [CrossRef]
  69. Pandey, G.; Lovejoy, S.; Schertzer, D. Multifractal analysis of daily river flows including exstremes for basins of five to two million square kilometres, one day to 75 years. J. Hydrol. 1998, 208, 62–81. [Google Scholar] [CrossRef]
  70. Mazzarella, A.; Diodato, N. The alluvial events in the last two centuries at Sarno, Southern Italy: Their classification and power-law time-occurrence. Theor. Appl. Climatol. 2002, 72, 55–64. [Google Scholar] [CrossRef]
  71. Mazzarella, A.; Rapetti, F. Scale invariance laws in the recurrence interval of extreme floods: An application to the Upper Po River Valley (Northern Italy). J. Hydrol. 2004, 288, 264–271. [Google Scholar] [CrossRef]
  72. Hurst, H.E. The long term storage capacity of reservoirs. Trans. Am. Soc. Civ. Eng. 1951, 116, 770–808. [Google Scholar]
  73. Pasculli, A.; Sciarra, N. A fragile fracture study by fractal analysis. In Proceedings of the Annual Conference of the International Association for Mathematical Geology IAMG 2005, Toronto, ON, Canada, 21–26 August 2005. [Google Scholar]
  74. Liu, L.; Chen, X.; Xu, G.; Shu, L. Use of hydrologic time-series data for identification of hydrodynamic function and behavior in a karstic water system in China. Hydrogeol J. 2011, 19, 1577–1585. [Google Scholar] [CrossRef]
  75. Zhang, Y.; Schilling, K. Temporal scaling of hydraulic head and river base flow and its implication for groundwater recharge. Water Resour. Res. 2004, 40, W03504. [Google Scholar] [CrossRef]
  76. Shun, T.; Duffy, C.J. Low-frequency oscillations in precipitation, temperature, and runoff on a west facing mountain front: A hydrogeologic interpretation. Water Resour. Res. 1999, 35, 191–201. [Google Scholar] [CrossRef]
  77. Markovic, D.; Koch, M. Stream response to precipitation variability: A spectral view based on analysis and modelling of hydrological cycle components. Hydrol Process. 2015, 29, 1806–1816. [Google Scholar] [CrossRef]
  78. Zhang, Y.-K.; Li, Z. Temporal scaling of hydraulic head fluctuations: Nonstationary spectral analyses and numerical simulations. Water Resour. Res. 2005, 41, W07031. [Google Scholar] [CrossRef]
  79. Neves, M.C.; Costa, L.; Monteiro, J.P. Climatic and geologic controls on the piezometry of the Querença-Silves karst aquifer, Algarve (Portugal). Hydrogeol J. 2016, 24, 1015–1028. [Google Scholar] [CrossRef]
  80. Molénat, J.; Davy, P.; Gascuel-Odoux, C.; Durand, P. Study of three subsurface hydrologic systems based on spectral and cross-spectral analysis of time series. J. Hydrol. 1999, 222, 152–164. [Google Scholar] [CrossRef]
  81. Slimani, S.; Massei, N.; Mesquita, J. Combined climatic and geological forcings on the spatio-temporal variability of piezometric levels in the chalk aquifer of Upper Normandy (France) at pluridecennal scale. Hydrogeol. J. 2009, 17, 1823–1832. [Google Scholar] [CrossRef]
  82. Gámiz-Fortis, S.; Pozo-Vázquez, D.; Trigo, R.M.; Castro-Díez, Y. Quantifying the Predictability of Winter River Flow in Iberia. Part I: Interannual Predictability. J. Clim. 2008, 21, 2484–2502. [Google Scholar] [CrossRef]
  83. Labat, D. Cross wavelet analyses of annual continental freshwater discharge and selected climate. J. Hydrol. 2010, 385, 269–278. [Google Scholar] [CrossRef]
  84. Holman, I.P.; Rivas-Casado, M.; Bloomfield, J.P. Identifying non-stationary groundwater level response to North Atlantic ocean-atmosphere teleconnection patterns using wavelet coherence. Hydrogeol. J. 2011, 19, 1269–1279. [Google Scholar] [CrossRef] [Green Version]
  85. Tremblay, L.; Larocque, M.; Anctil, F.; Rivard, C. Teleconnections and interannual variability in Canadian groundwater levels. J. Hydrol. 2011, 410, 178–188. [Google Scholar] [CrossRef]
  86. El Janyani, S.; Massei, N.; Dupont, J.P.; Fournier, M.; Dörfliger, N. Hydrological responses of the chalk aquifer to the regional climatic signal. J. Hydrol. 2012, 464, 485–493. [Google Scholar] [CrossRef]
  87. Briciu, A.E.; Mihăilă, D. Wavelet analysis of some rivers in SE Europe and selected climate indices. Environ. Monit. Assess. 2014, 186, 6263–6286. [Google Scholar] [CrossRef] [PubMed]
  88. Probst, J.L.; Tardy, Y. Long range streamflow and world continental runoff fluctuations since the beginning of this century. J. Hydrol. 1987, 94, 289–311. [Google Scholar] [CrossRef]
  89. Labat, D. Oscillations in land surface hydrological cycle. Earth Planet. Sci. Lett. 2006, 242, 143–154. [Google Scholar] [CrossRef]
  90. Fendeková, M.; Pekárová, P.; Fendek, M.; Pekár, J.; Škoda, P. Global drivers effect in multi-annual variability of runoff. J. Hydrol. Hydromech. 2014, 62, 169–176. [Google Scholar] [CrossRef]
  91. Chiaudani, A.; Di Curzio, D.; Di Lena, B.; Palmucci, W.; Rusi, S. Preliminary analysis regarding long period rainfall-piezometric levels relationships in Abruzzo Region Peri-Adriatic area. In Proceedings of the Acts XIX Agrometeorology National Congress: New Adversities and New Services for Agroecosystems, Bologna, Italy, 14–16 June 2016. [Google Scholar]
  92. ISPRA-Istituto Superiore per la Protezione e la Ricerca Ambientale. Linee Guida per L’analisi e L’elaborazione Statistica di Base Delle Serie Storiche di Dati Idrologici; Manuali e Linee Guida 84/13; ISPRA: Roma, Italy, 2013.
Figure 1. (Top) The study area; and (Bottom) the hydrogeological framework and the monitoring network are shown (see Table 1 for initials and abbreviations).
Figure 1. (Top) The study area; and (Bottom) the hydrogeological framework and the monitoring network are shown (see Table 1 for initials and abbreviations).
Water 09 00850 g001
Figure 2. Schematic geolithological section with a 10× vertical exaggeration (Cross section 1–1′ in Figure 1).
Figure 2. Schematic geolithological section with a 10× vertical exaggeration (Cross section 1–1′ in Figure 1).
Water 09 00850 g002
Figure 3. Morphologic and hydraulic characteristics of the river profile assessed in this study: (A) distribution of the river width; (B) plan view of the elevation of the river water profile (GPS survey, April 2015); and (C) elevation profile of the river water head.
Figure 3. Morphologic and hydraulic characteristics of the river profile assessed in this study: (A) distribution of the river width; (B) plan view of the elevation of the river water profile (GPS survey, April 2015); and (C) elevation profile of the river water head.
Water 09 00850 g003
Figure 4. Locations, stratigraphy and piezometric levels or the assessed wells.
Figure 4. Locations, stratigraphy and piezometric levels or the assessed wells.
Water 09 00850 g004
Figure 5. Rainfall (A); river level (B); and piezometric level (C) for 1986–2009 three-day time series.
Figure 5. Rainfall (A); river level (B); and piezometric level (C) for 1986–2009 three-day time series.
Water 09 00850 g005
Figure 6. Monthly average: for rainfall and hydrometric level (A); and for normalized piezometric level (B), 1986–2009 period.
Figure 6. Monthly average: for rainfall and hydrometric level (A); and for normalized piezometric level (B), 1986–2009 period.
Water 09 00850 g006
Figure 7. Residual time series (left); and year type (right).
Figure 7. Residual time series (left); and year type (right).
Water 09 00850 g007
Figure 8. Autocorrelation functions: rainfall-river level (A); and piezometric levels (B). The dash-dotted line indicates the decorrelation threshold. The initials and the abbreviations refer to Table 1.
Figure 8. Autocorrelation functions: rainfall-river level (A); and piezometric levels (B). The dash-dotted line indicates the decorrelation threshold. The initials and the abbreviations refer to Table 1.
Water 09 00850 g008
Figure 9. Cross-Correlations Function: rainfall–piezometric level (A); and river level–piezometric level (B). The dash-dotted lines indicate the significant thresholds.
Figure 9. Cross-Correlations Function: rainfall–piezometric level (A); and river level–piezometric level (B). The dash-dotted lines indicate the significant thresholds.
Water 09 00850 g009
Figure 10. Residual Cross-Correlation Function: rainfall–hydrometric level. The dash-dotted lines indicate the significant thresholds.
Figure 10. Residual Cross-Correlation Function: rainfall–hydrometric level. The dash-dotted lines indicate the significant thresholds.
Water 09 00850 g010
Figure 11. Residual Cross-Correlation Functions: rainfall–piezometric level (A); and river level–piezometric level (B). The dash-dotted lines indicate the significant thresholds.
Figure 11. Residual Cross-Correlation Functions: rainfall–piezometric level (A); and river level–piezometric level (B). The dash-dotted lines indicate the significant thresholds.
Water 09 00850 g011
Figure 12. Spectral density functions: rainfall and hydrometry levels (A); and piezometric levels (B).
Figure 12. Spectral density functions: rainfall and hydrometry levels (A); and piezometric levels (B).
Water 09 00850 g012
Figure 13. Rainfall and river level outcomes; fractal dimensions related to different F_threshold events; rainfall F_threshold > 40%: DF 0.00 ; river level F_threshold < 5%: DF 0.9999 .
Figure 13. Rainfall and river level outcomes; fractal dimensions related to different F_threshold events; rainfall F_threshold > 40%: DF 0.00 ; river level F_threshold < 5%: DF 0.9999 .
Water 09 00850 g013
Figure 14. Su, Sa, DN outcomes related to piezometric levels; Fractal dimensions related to different F_threshold events.
Figure 14. Su, Sa, DN outcomes related to piezometric levels; Fractal dimensions related to different F_threshold events.
Water 09 00850 g014
Figure 15. Estimated DRT vs. f i , j ( % ) values.
Figure 15. Estimated DRT vs. f i , j ( % ) values.
Water 09 00850 g015
Figure 16. Comparison of borehole hydraulic head changes and the distribution of river level elevations. The borehole piezometric heads, which vary most frequently between the 25th and 75th percentile, correspond to the river bed elevations along the blue sections of the river profile.
Figure 16. Comparison of borehole hydraulic head changes and the distribution of river level elevations. The borehole piezometric heads, which vary most frequently between the 25th and 75th percentile, correspond to the river bed elevations along the blue sections of the river profile.
Water 09 00850 g016
Figure 17. Refined conceptual model of groundwater flow paths and groundwater/surface water relationships.
Figure 17. Refined conceptual model of groundwater flow paths and groundwater/surface water relationships.
Water 09 00850 g017
Table 1. Topographical data and type of monitored sites.
Table 1. Topographical data and type of monitored sites.
NameIDTypeElevation (m a.s.l.)Distance from River (m)Distance from Sea (km)
SpoltoreSpRain gauge165.0026504.4
SurricchioSuWell18.102508.3
SanitàSaWell8.681704.1
De NicolaDNWell11.6819603.7
Pescara S. TeresaPSTHydrometer4.5107.0
Table 2. Average and standard deviation of the piezometric level, rainfall and hydrometric level time series data.
Table 2. Average and standard deviation of the piezometric level, rainfall and hydrometric level time series data.
TypeNamePiezometric Level (m a.s.l.)Water Depth (m)
AverageSt. dev.MinMaxAverage
WellSu12.061.429.2916.286.70
Sa5.490.434.417.283.89
DN10.460.199.8511.182.02
TypeNameRainfall (mm)
AverageSt. dev.MinMax
Rain GaugeSp5.3511.90199.2
TypeNameHydrometric Level (m a.s.l.)
AverageSt. dev.MinMax
Hydrometer (m)PST3.330.372.747.97
Table 3. Details of the Fractal analysis; d: days, w: week, m: months, y: years, F_threshold: the actual level (m) corresponding to the selected F_threshold f i , j ( % ) level (percentage of the whole variability range).
Table 3. Details of the Fractal analysis; d: days, w: week, m: months, y: years, F_threshold: the actual level (m) corresponding to the selected F_threshold f i , j ( % ) level (percentage of the whole variability range).
Selected F_threshold LevelSp RainfallPST HydrometrySu WellSa WellDN Well
Range =(0–199.2) (mm)(2.74–7.97) (m)(9.29–16.28) (m)(4.41–7.28) (m)(9.85–11.18) (m)
L_thresh = min + 0.5 × range/100 (m) =0.0012.779.334.429.86
Fluctuation (m) =0.0010.030.040.010.01
n_events: fi,j ≥ F_threshold = 0.5% =11702876287728792879
percentage (%) =1170/2880 = 40.699.8699.8999.9699.97
distribution return time (DRT) (d) =17 (2 w)3333
L_thresh = min + 1.0 × range/100 (m) =0.0022.799.364.449.86
Fluctuation (m) =0.0020.050.070.030.01
n_events::fi,j ≥ F_threshold = 1.0% =9982860287428792879
percentage (%) =34.799.3199.7999.9699.97
DRT (d) =30 (1 m)3333
L_thresh = min + 2.5 × range/100 (m) =0.0052.879.474.489.88
Fluctuation (m) =0.0050.130.180.070.03
n_events: fi,j ≥ F_threshold = 2.5% =7302792286228722878
percentage (%) =25.396.9499.3899.7299.93
DRT (d) =47 (1.7 m)3333
L_thresh = min + 5.0 × range/100 (m) =0.013.009.644.559.92
Fluctuation (m) =0.010.260.310.140.07
n_events::fi,j ≥ F_threshold = 5% =4742491279428662876
percentage (%) =6.586.4997.0199.5199.86
DRT (d) =71 (2.4 m)3333
L_thresh = min + 10 × range/100 (m) =0.023.269.994.709.98
Fluctuation (m) =0.020.520.70.290.13
n_events: fi,j ≥ F_threshold = 10% =2521472268427982861
percentage (%) =8.851.1193.1997.1599.34
DRT (d) =107 (3.6 m)95 (3.2 m)333
L_thresh = min + 20 × range/100 (m) =0.043.7910.694.9910.12
Fluctuation (m) =0.041.051.40.580.27
n_events: fi,j ≥ F_threshold = 20% =74218231326712837
percentage (%) =2.67.56980.3192.7498.51
DRT (d) =361 (12 m)428 (14.3 m)69 (2.3 m)33
L_thresh =min + 30 × range/100 (m) =0.064.3211.405.2710.25
Fluctuation (m) =0.061.582.110.860.40
n_events: fi,j ≥ F_threshold = 30% =2647182419472618
percentage (%) =0.91.63263.3367.6090.9
DRT (d) =727 (24 m)572 (19 m)176 (5.9 m)130 (4.3 m)30
L_thresh = min + 40 × range/100 (m) =0.084.8412.105.5610.38
Fluctuation (m) =0.082.12.811.150.53
n_events: fi,j ≥ F_threshold = 40% =817138211391859
percentage (%) =0.280.59047.9839.5564.55
DRT (d) =1995 (67 m)861 (29 m)589 (19.6 m)360 (12 m)282
L_thresh = min + 60 × range/100 (m) =0.125.8813.486.1310.65
Fluctuation (m) =0.123.144.191.720.8
n_events: fi,j ≥ F_threshold = 60% =17518216457
percentage (%) =0.0350.24317.987.515.87
DRT (d) =8640 (24 y)1392 (46 m)764 (25.5 m)1272 (42 m)727 (24 m)
L_thresh = min + 80 × range/100 (m) =0.1606.9214.886.7110.91
Fluctuation (m) =0.1604.185.592.31.06
n_events: fi,j ≥ F_threshold = 80% =13583563
percentage (%) =0.0350.1042.0141.2152.19
DRT (d) =8640 (24 y)2694 (7.5 y)2694 (7.5 y)2860 (7.9 y)8640 (24 y)
L_thresh = min + 90 × range/100 (m) =0.1807.4515.586.9911.05
Fluctuation (m) =0.1804.716.292.581.20
n_events: fi,j ≥ F_threshold = 90% =1291321
percentage (%) =0.0350.0690.3130.4510.73
DRT (d) =8640 (24 y)4310 (12 y)3636 (10 y)4186 (11.5 y)8640 (24 y)

Share and Cite

MDPI and ACS Style

Chiaudani, A.; Di Curzio, D.; Palmucci, W.; Pasculli, A.; Polemio, M.; Rusi, S. Statistical and Fractal Approaches on Long Time-Series to Surface-Water/Groundwater Relationship Assessment: A Central Italy Alluvial Plain Case Study. Water 2017, 9, 850. https://doi.org/10.3390/w9110850

AMA Style

Chiaudani A, Di Curzio D, Palmucci W, Pasculli A, Polemio M, Rusi S. Statistical and Fractal Approaches on Long Time-Series to Surface-Water/Groundwater Relationship Assessment: A Central Italy Alluvial Plain Case Study. Water. 2017; 9(11):850. https://doi.org/10.3390/w9110850

Chicago/Turabian Style

Chiaudani, Alessandro, Diego Di Curzio, William Palmucci, Antonio Pasculli, Maurizio Polemio, and Sergio Rusi. 2017. "Statistical and Fractal Approaches on Long Time-Series to Surface-Water/Groundwater Relationship Assessment: A Central Italy Alluvial Plain Case Study" Water 9, no. 11: 850. https://doi.org/10.3390/w9110850

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