Statistical and analytical study of the propagation of flood-induced groundwater rise in an alluvial aquifer

Journal of

Hydrology E L S EV I E R

Journal of Hydrology 205 (1998) I 12- 125

Statistical and analytical study of the propagation of flood-induced groundwater rise in an alluvial aquifer Z. Vekerdy*, A.M.J. Meijerink International lnstitute j b r Aerospace Survey and Earth Sciences (1TC), P.O. Box 6, 7500 AA Enschede, The Netherlands Received 22 January 1997; revised 22 October 1997; accepted 20 November 1997

Abstract

The propagation of stage rises of the river Danube in an adjoining alluvial aquifer (Kisalf'old, NW Hungary) has been studied by preparing a map which shows the lag times corresponding to maximum correlation values between the hydrographs of groundwater observation wells and the river stages. As expected, the lag times for the unconfined aquifer generally exceed those for the semi-confined part. The complex pattern suggests that other rivers in the area also play a role and that transmissivities and resistances of a cover layer are not sufficiently well known, despite the presence of a relatively dense network of bore holes. A systematic increase of the lag times with the distance from the river was noted at two sections. Nine flood events have been selected to study how well observed head rises could be predicted using equations for the calculation of the propagation of the head for phreatic and for confined conditions. The results suggest that the water levels in the complex flood plain and the river bed conductance of silted up branches during the rise of the flood play an important role. © 1998 Elsevier Science B.V.

Keywords: Alluvial aquifer; Groundwater; Bank storage

I. I n t r o d u c t i o n

A rise in the water level o f a river traversing an alluvial aquifer with permeable layers induces flow from the river into the aquifer. When the river is bordered by alluvial terraces, most o f the water is returned to the river during low stages. This is known as bank storage and was quantitatively modelled, e.g. by Moench et al. (1974) and Hunt (1990). When the river flows on an alluvial fan, groundwater flow is generally away from the river, with the result that only little water returns back to the river when the stage is falling. Hence there is a net flux from the river to the groundwater, as is the case in the studied area. * Corresponding author. 0022-1694/98/$19.00 © 1998 Elsevier Science B.V. All rights reserved PH S 0 0 2 2 - 1 6 9 4 ( 9 7 ) 0 0 1 4 8 - 0

Early works concentrated on the setup and analytical solution o f formulae describing the river groundwater interaction, e.g. the draft on a river by pumping a well adjacent to it (e.g. Glover and Balmer, 1954). More recent developments include the use o f a three-dimensional coupled surface water and groundwater flow model, based on the numerical solution o f open channel flow and groundwater flow equations and allowing the modelling o f the surface w a t e r groundwater interaction. An overview o f the development o f such models is given e.g. by Vasiliev (1987). One o f the examples o f the recent developments in this field is shown by Swain (1994), who describes the coupling o f US Geological Survey computer codes B R A N C H and M O D F L O W to calculate the surface w a t e r - g r o u n d w a t e r interactions.

Z. Vekerdy, A.M.J. Mei/erink/dournal o/ IJvdrolog~ 205 (1998) 112-125

An analysis of the effects of the river stage rises could precede the complex modelling to get information about the spatial inhomogeneity of hydrologic parameters for grid design or for the calibration of a groundwater model. The model calibration is usually difficult in alluvial areas in the zones along the river due to the complicated interactions between river stages and groundwater levels. Simplification and generalisation are usually required. It is useful to investigate whether simple analytical approaches can be of assistance. A general method to understand the response of groundwater heads on stage changes in rivers may follow two steps: first a statistical description of the phenomenon is made, then a hydraulic analysis. This approach was used by Sophocleous (1991), who investigated the pressure propagation along buried paleochannels in Kansas, USA. He demonstrated by conceptual modelling that a pressure wave generated by a river may propagate over tens of kilometres in the confined paleochannels. A simplification of a complex regional problem is possible by reducing the number of dimensions. This is feasible in the case of rivers since the floodinduced head wave propagates in a well-defined direction, namely perpendicularly away from the river, superposed on the general groundwater flow pattern. In case of a single aquifer with uniform vertical distribution of the heads, a one-dimensional analysis is feasible. A further simplification in the presence of much spatial variation of hydrogeologic parameters, often not known from available data sets, is to use equations with generic parameters (in the following referred to as simple equations). Complex river

river

ones would require much guesswork in the absence of a proper observation network. The results obtained by the simple equations could help in finding out what aspects should be included for further improvement. Analytical solutions of the flow equations are known for such cases (e.g. Edelman, 1947; Heij, 1989; Leenen, 1994; Zhang, 1992). The first objective of this study is to analyse the pattern of the propagation of high water levels in the aquifer away from the river after a linear rise of the stage by using lag correlation. After that, the results of equations for the propagation of the heads for phreatic (Edelman, 1947) and for confined (Leenen, 1994) conditions are compared with observed heads at selected regions.

2. River-aquifer boundary on alluvial regions and description of the study area In case the aquifer is predominantly composed of coarse grained materials, essentially four types of boundary conditions can be recognised along a river, as sketched in Fig. 1. Most alluvial areas are seldom uniform in plan and profile because of changing climatic conditions in the provenance areas during the time of deposition. Lithologic logs of shallow boreholes often show important variations at relatively short distances. Therefore, the type of boundary condition along the river can be difficult to define and the same is true for the variations in transmissivity. The presence of a cover layer [Fig. l(c)] with fine grained deposits resting on thick, coarser grained, braided river deposits is characteristic for many river

river

I\.oo'.ol ~"

12. a

a.) b.) c.) d.)

b

t 13

c

•

D"

d

Highly resistant riverbed

River in highly resistant cover layer (Semi-)confined groundwater, permeable riverbed Phreatic groundwater, permeable riverbed

Fig. 1. Possible contacts between a river and the groundwater table in alluvial regions.

114

Z. Vekerdy, A.M.J. Meijerink/Journal ~?]'Hydrolog, y 205 (1998) 112-125

alluvial areas. The cover layer varies usually in thickness and has variations in textures because it is composed of overbank suspended load deposits and channel fill deposits. This situation is found on the alluvial fan of the Danube river in the KisaW61d region (Small Plain, in the north-western part of Hungary). In the studied area the cover layer is coarse in the upstream part of the fan and the groundwater is phreatic. In the downstream part of the fan the cover layer is fine grained, resulting in semi-confined groundwater. The available data indicate that mainly river-aquifer boundary types (c) and (d) of Fig. 1 occur, i.e. the main branch of the Danube is not sealed from the groundwater by an impervious river bed. It should be mentioned that, despite the availability of a fairly dense network of shallow (approximately l0 m deep) boreholes with lithologic descriptions, there remains much uncertainty concerning the local spatial variation of the hydrogeologic variables (Vekerdy, 1996).

heavier textures of the cover layer. As expected, the propagation under confined or semi-confined conditions was faster (i.e. at the downstream part of the alluvial fan) than in the phreatic (upstream) part. There was a tendency for lag times to increase with distance away from the river. However, some erratic values can be noted on the map. These can be attributed mainly to the influence of the Moson-Danube River (see Fig. 2). Other factors can also have effects, such as the spatially variable transmissivity of the lower gravels and sands; the spatially variable resistivity of the cover layer; non-steady flow at the beginning of the head rises; and the transformation of phreatic to (semi-)confined conditions at places, where rising groundwater from the lower gravels enters into the cover layer during the rise. Along two sections, line #6, located in the (semi-)confined part, and line #2, located in the phreatic part of the aquifer, the lag time increases essentially regularly with distance. These two lines have been selected for further scrutiny.

3. Statistical analysis of the groundwater head propagation

3.2. Time lags along selected lines

Records of daily water levels in observation wells, of the entire period from 1979 to 1989, have been used for lag correlation. The correlation coefficients (Pearson's r) have been calculated for the stage heights of the river and the heads in the wells with different time lags (1,2,...,31,33,35 ..... 85 days) for each individual well.

3.1. The regional pattern of the propagation as expressed by lag correlation The characteristic time lag was determined as the highest value of the correlation coefficient for each observation well. A choice was made to consider only those wells that have a correlation coefficient of more than 0.2, because with lower values of the coefficient it was difficult to define the characteristic time lag. The characteristic time lags have been plotted on a map (Fig. 2) together with the identification numbers of the wells and well-lines. When the result of the time lag analysis is compared with the lithological description of bore logs, it is apparent that lower lag times are associated with

The lag correlation coefficients have been plotted against the lag times for wells at various distances from the Danube, along line #6 (Fig. 3) and line #2 (Fig. 4). In both sections, the curves of the correlation coefficients show a clear maximum or culmination only at short distances from the river, particularly for the confined condition. Considerable dampening of the effect of the pressure propagation takes place with increasing distances, which is expressed by a decrease of the lag correlation coefficient and by a broadening of the culmination of the curves. For the unconfined part of the aquifer, culminations of the curves can be easily noted for the three wells closest to the river. For wells further away the culminations are broad and the maximum value was determined by eye. The time lag corresponding to the culmination (maximum value) was considered as the characteristic lag time and is plotted on the map of Fig. 2. The generally flat culmination of these curves is due to the dampening effect, and possibly to the effect of surface recharge. When the characteristic lag times, as deduced from the graphs of Fig. 3 and Fig. 4, are plotted against the distance of the wells from the river along line #6

Z. Vekerdy, A.M.J. Me(jerink/Journalqf'Hydrology 205 (1998) 112-125

115

N

#1 #2 Well-line #3

~,,

,®

Well-line#4 /

0

~

lO km

Groundwater well identifier Characteristic time lag

/

d

/

I

Well-line #5

I

/

Well-line #6

Well-line #7

Fig. 2. Characteristic time lags between the river stages in the Danube and the groundwaterlevels measured in the observationwells. and line #2, a consistent pattern emerges (Fig. 5). The lag time for the phreatic conditions increases with distance. The overall relationship is not fully linear. For the semi-confined conditions sufficient data are not available to express the relationship for larger distances. However, the graph shows that a rise in the stage of the river is propagated faster than under the phreatic condition.

4. Analysis of flood-induced groundwater head rise based on hydraulic equations The selection of wells along the two lines, where the heads showed some consistency in the response to rises in the river stage, has been done to allow comparison of

observed heads with heads predicted by equations based on theory. They will now be discussed. 4.1. Semi-confined condition

The analytical solution of the equation of the onedimensional groundwater flow given by Leenen (1994) was used for the analysis of the groundwater response to the stage rise in the Danube. The derivation of the equation is described in Appendix A. The solution gives the groundwater head rise caused by the rise of the flood stage. Analysis showed that the heads are sensitive to the distance from the river, the transmissivity of the aquifer and the resistance of the cover layer (Fig. 6), but less sensitive to the elastic storage of the aquifer.

[ 16

Z Vekerdy, A.M.J. Meijerink/Journal q/Hvdrology 205 (1998) 112 125

0.8

0.6

0.4

0.2

0

5

10

15

20

25

30

35

40 45 50 lag [days]

55

60

65

70

75

80

Well 63 (1630 m from the Danube) ~

Well 64 (1930 m from the Danube)

Well 65 (2170 m from the Danube) ~

Well 66 (2050 m from the Danube)

Well 62 (2450 m from the Danube) ~

Well 68 (3640 m from the Danube)

- - I - Well 69 (4680 m from the Danube) ~

Well 71 (5700 m from the Danube)

85

Well 72 (5800 m from the Danube)

Fig. 3. Lag correlation coefficients (r) calculated with different time lags along well-line #6 (semi-confined part of the aquifer).

0.8 i

0,6

"0.4

0.2

,==:I:;::,,,:=I::::I:::==;=;=I====I=:;;1===~ . . . . . . . . . . . =I= :: :I : : : : I ; = : ; I:: : :I: : ::I :==:==;=: 0

10

15

20

25

30

35

40 45 50 Lag [day]

55

60

65

70

75

80

Well #129 ( 250 m from the Danube)

Well #128 ( 480 m from the Danube)

Well #126 (1150 m from the Danube)

Well #137 (1600 m from the Danube)

- = - Well #135 (3700 m from the Danube)

- + - Well #134 (6100 m from the Danube)

Well #133 (7600 m from the Danube)

- , * - Well #132 (8800 m from the Danube)

85

--v-- Well #112 (16 500 m from the Danube)

Fig. 4. Lag correlation coefficients (r) calculated with different time lags along well-line #2 (phreatic part of the aquifer).

Z. Vekerdy. A M.J. Mei/erink/Journal ()['Hvdrologv 205 (1998) 112 125 70 A

60 .•.'...

50 40 --I ~

30

....................................

20

......................................

10

....................................

i,•

o ~ • : n.~z 0

,

i

,

i

,

i

2000 4000 6000 8000 Distance from the Danube [m] O Well-line # 6 (confined aquifer)

I

117

determination of Kz•cVThe hydraulic properties along the well-line were determined using a map of hydraulic terrain units (HU), which was compiled by crossing of the map with classes of the cover layer resistance (Fig. 8) and the map with classes of aquifer transmissivity (Fig. 9) in a geographical information system (GIS). A schematic cross-section (Fig. 10) illustrates the derivation of HUs. Since the pressure head propagation in a semi-confined aquifer is fast (exceeding several kilometres per day), it has been assumed that in the neighbourhood of the river the characteristic times are equal in each consecutive HU: t , l ~ t , 2 --~ . . . ~ l,j 7

• Well-line # 2 (phreatic aquifer) I

Fig. 5. D e p e n d e n c e o f the characteristic time lag on the distance from the Danube.

Some observation wells on the semi-confined, downstream part of the Danube fan [along and close to well-line #6, Fig. 7(b)] have been selected for comparing heads calculated by Eq. (6), with the observed ones. The confining cover layer in this part of the aquifer has a high resistance (80-400 day), as determined by Eq. (2), based on lithologic logs of the cover layer• Data from literature describing infiltration measurements and other field data of the Kisalf61d ( Erd61yi, 1971; Maginecz, 1993; Papp, 1988; Sz~kely, 1977; (lj faludy, 1986) were used for the

®

0.8 .~0.6 0.4 0.2

=®

1. continuous (quasi-)linear rise of the river stage; 2. relatively stable low water period preceding the flood; 3. the selected floods represent a range of different maximum stages. The results (Fig. 11) show a reasonable fit for smaller floods, but there is a systematic deviation of

g

-~ 1

o

where n is the number of HUs. Elastic storage was assumed to be 0.001, and the river-aquifer boundary as sketched in Fig. l(c). Nine flood waves were selected for the analysis of the head propagation (Table 1), using the following criteria:

0

~

.~0.6 ~ 0.4

P "o 0.2 = 0

z

I

~ 0.8

z

0

t)

J Fig. 6. Sensitivity o f the normalised g r o u n d w a t e r head to the different parameters in the calculation method o f Leenen (1994) (S~ = elastic storage o f the aquifer, T = transmissivity o f the aquifer, t* = normalised duration o f the flood, col = resistance o f the c o v e r layer).

Z. Vekerdy, A.M.J. Meijerink/Journal o["Hydrology 205 (1998) 112-125

118

b

groundwater

~25

~

observation well

dike

river ~"

0

10 km

small stream

Fig. 7. Location of the wells used in the hydraulic equation-based analysis of flood-induced groundwater head rise: (a) phreatic part of the aquifer; (b) semi-confined part of the aquifer.

calculated heads compared to the observed ones for higher floods. The main explanation of the deviation lies in the flood plain geometry (Fig. 12). When the flood peak is not higher than the river bank, the distance between the infiltration front of river water and the observation wells is rather constant. At higher

stages the flood waters start spreading over the flood plain, first in the flood plain channels and later over the full width of the flood plain, up to the dikes. Because of the complex geometry this process is gradual, but on average the distance from the infiltration front to the wells is shortened. It is difficult to

N 10km

0

•

Groundwater observation well Dike

Confining layer resistance: ?iii~i~ a

1 - 20 day 20-

40 day

80 day 80- 160 day 160- 320 day 320-640 day 40-

BIB

Fig. 8. Confined parts of the aquifer with the resistance of the confining cover layer. The map was compiled from data provided by the Geological Survey of Hungary and the Water Resources Research Centre (VITUKI), Budapest.

119

Z. 1,2,kera~v,A.M.J. Me(/erink/Journal o/"Hydrolog3' 205 (1998) 112 125

0

e

i

10 km

Groundwater observation well Dike

Transmissivity of the aquifer: 1-

200 m2/day

iii!i

200 -

400 m2/day

'!~,~i,

400 -

800 m2/day

800 -

1600 m2/day

' I

1600 -

3200 m2/day

3200-

6400m2/day

I

6400 - 12 800 m2/day

I

12 800 - 25 600 m2/day

I

25 600 - 51 200 m2/day

I

> 51 200 m2/day

Fig. 9. Transmissivity of the aquifer. The map was compiled from data provided by the Geological Survey of Hungary and the Water Resources Research Centre (VITUKI), Budapest. consider this complex geometry in the calculations, but it could explain why the calculated heads are lower than the measured ones. 4.2. P h r e a t i c c o n d i t i o n

Edelman (1947) gives an analytical solution o f the equation o f one-dimensional groundwater flow in an

aquifer bordered by a water body. The most important steps o f the solution are given in Appendix A. The solution gives the rise o f the groundwater table during the period o f the river stage rise (when c~ is positive, see Appendix A), or river stage decline [when c~ is negative in Eq. (12)]. The phreatic model is sensitive to transmissivity and distance from the river, as shown in Fig. 13. It

HU1 Xl

~

°

a.)

T1

HU2 "

"~'~

°i.

x2

HU3 :

~~

x3

"

~_"

"T2 "i. = "T3°

o ~

b.)

Fig. 10. Scheme of the derivation of the hydraulic terrain units (HU): (a) plan; (b) section.

120

Z. Vekerdy, A.M.J. Meijerink/Journal ()/'Hydrologo, 205 (1998) 112 125

Table 1 Main characteristics of the floods selected for the analysis Flood duration From

To

01-Feb-85 19-Jul-81 18-Dec-87 10-Mar-81 12-Jan-86 02-Feb-80 04-Feb-81 28-Feb-83 28-Nov-83

05-Feb-85 25-Jul-81 2 l-Dec-87 16-Mar-81 16-Jan-86 08-Feb-80 1O-Feb-81 03-Mar-83 0 l-Dec-83

Duration of rise, tp (day)

Full rise,

4 6 3 6 4 6 6 3 3

4.5 4.45 4.09 3.88 2.98 2.74 2.63 1.58 1.35

Hp (m)

station (at Rajka) were projected to the side branches. Distances from the corresponding river branch to the wells were taken perpendicularly to the river. Transmissivity values encountered along the flow lines have been calculated based on data provided by the Water Resources Research Centre (VITUKI), Budapest. The specific yield was assumed to be 0.25. The calculated head rises in the wells have been plotted against the measured ones, for the nine events mentioned above, see Fig. 14. The positions of the curves reflect the lower heads at increasing distance from the river. The increase of calculated head tallies almost with that of the

is not particularly sensitive to the variations in the specific yield, for which a general value can be adopted. In the upstream part of the Danube alluvial fan the groundwater is mostly phreatic. The main river reach, where direct river stage measurements are available, is surrounded by a wide flood plain up to the dikes. Large side branches of the river are located on the flood plain just at the foot of the dikes. The selected wells #125, #128, #129 and #140 are in the vicinity of the side branches of the river on the protected side of the dikes, in the phreatic part of the aquifer [Fig. 7(a)]. The river stage data of a nearby gauging

J i

"E 2

,,"/

':/"

"o Q

/, .~

.,"

X

./

:

........ 0"~

0

,

L

~

X o o

i

,

,

J

,

I 2 Measured rise [m]

Wel164 (1930 m from the Danube)

~

,

,

t

3

~ Wel165 (2170 m from the

Well 66 (2050 m from the Danube)

Fig. l l. Comparison of" measured and calculated groundwater rise in wells as a result of nine selected floods of the Danube (semi-confined aquifer). Every symbol represents one flood event.

Z. k2~kerdv, A.M.J. Me!/erink/Journal qf 14vdrolo~, 205 (I 998) 112-125 X2

:< ;

X1

< Hp2;

l Xl

121

~.

~2

H

p

> x2

groundwater head Fig. 12. Scheme of the different distances (x) between the infiltration front of the river water and the observation well at different flood peaks

(Hp). measured head, because the slopes of the curves show an approximate one-to-one relation, except for the farthest well at 820 m distance, where the calculated rise was nil. However, the graph shows that small floods of the river, represented by data points on the left side of Fig. 14, do not affect the observed water levels in the wells• This can be explained by the presence of low dams across the upstream end of the side branches, so small floods cannot enter these branches from their upstream end. The river stage in the branches rises as a result of backwater from their

downstream end. As Fig. 7(a) shows, the upstream ends of the branches have no contact with the main channel. For larger floods, the calculated values are higher than the observed ones, particularly for the wells at distances of 150 m (well #140) and 480 m (well #128) from the Danube. This can be partly explained by the delay of the actual water rise in the branches caused by the low dams at the upstream end of the branches, and partly by the fact that the river bed conductance was neglected in the calculations. The branches are

O) rCa

t-

5

~5

4

~" ~(.1 4

3

"~3

2 tO

0

"

t9

+,,t

otl~

the

river[~]

Fig. 13. Sensitivity of the groundwater head to the changes of different parameters in the calculation method of Edelman (1947).

Z Vekerdv, A.M.J. Me!ierink/Journal of Hydrology 205 (1998) 112 125

122

3

0

×

1

:

,

<×

0.~== 0

',

',

1

i 2

:

i

,

3

Measured [m]

i

Wel1140 (150 Well 129 (250 Well 128 (480 Well 125 (820

m m m m

far from the Danube) far from the Danube) far from the Danube) far from the Danube)

Fig. 14. Comparison of measured and calculated groundwater rise in wells as a result of nine selected floods of the Danube (phreatic aquifer). Every symbol represents one flood event.

partly silted up because they are protected from the main current. In the main river bed current velocities are high and the gravels are dredged for the maintenance of the navigational route. Lower actual river bed conductance of the branches may result in overestimation of the groundwater rise in the wells using Eq. (12). This is in agreement with the conclusion of Hantush (1965) who developed a procedure for the calculation of drawdown distribution and the rate and total volume of river depletion resulting from pumping near a stream with semi-pervious bed. He points out that if the effect of semi-perviousness of the river bed is ignored, the calculated head will be too high.

5. S u m m a r y

and c o n c l u s i o n s

A data set of river level rises after periods with low stages of the Danube River in north-western Hungary was used to study the propagation of high river stages in the aquifer. Relatively great numbers of geologic bore hole logs and observation well records were available. The spatial pattern of the groundwater head changes caused by the Danube was analysed statistically using lag correlation. A complex pattern was obtained when the time lag for the highest lag correlation between river stages and the selected well rises

was plotted for the groundwater wells on a map of the study area. It was possible to separate the response for the phreatic part of the aquifer and that for the part where confined or semi-confined conditions occur. The latter could be judged from the textures and the groundwater table with respect to the thickness of the cover layer. Deviations of observed characteristic time lags from expected spatial pattern could be attributed to the presence of smaller rivers, drainage channels and - - unknown - - local variables of hydrogeological conditions. Study of the lag times, as observed along lines of wells in the phreatic and in the semi-confined part, showed clear responses only at relatively short distances from the river, especially under semi-confined conditions. For the phreatic part of the aquifer, the response was difficult to detect beyond a distance of about 9 km because of the strong dampening of the river-induced pressure wave, and probably the superimposed effect of the surface recharge. However, when the lag times corresponding to the highest correlation coefficient were plotted against the distance from the river Danube, it could be demonstrated that a defined functional relationship existed. In other words, physical laws which govern propagation of a high river stage into the aquifer had an effect on the heads. In the phreatic part the effect could be detected up to at least a distance of 8800 m. In an attempt to predict the propagation of river stage rises, two equations have been selected because of their simplicity, one for confined conditions (Leenen, 1994) and one for phreatic conditions (Edelman, 1947). For the part of the aquifer where the equation for confined conditions has been used, the calculated results matched the observed ones only for river stage rises not higher than the river bank. For overbank flows the calculated values were lower than the measured ones. This can, at least partly, be attributed to the complex geometry of the flood plain. Uncertainties in the determination of the hydraulic parameters of the cover layer and the aquifer may play a role too. In the phreatic part of the aquifer the calculated heads were generally higher than what was observed. Relatively small floods measured at the closest gauging station had no observed effects on the groundwater, whilst the calculations predicted some groundwater head rise. It can be assumed that the

Z Vekerdy, A.M.J. Me!/erink/Journal q/Hvdrologv 205 (1998) 112 125

small dams across the upstream end of the river branches prohibited the smaller floods from entering these branches. Larger floods had smaller observed effects than predicted by the calculations that are attributed to the delay of river stage rise in the branches near the wells. The possibly low river bed conductance can explain why calculated rises are too high. The one-dimensional hydraulic analysis has shown that, apart from local variations of transmissivity, accurate data about the surface hydraulic configuration are needed. The analysis indicated what type of additional field data are needed, and showed that there can be little confidence in calibration results of a three-dimensional groundwater model if those data are not available.

Acknowledgements We thank Mr J.C. van Dam for remembering and retrieving the work of Edelman used in this study and for his critical attitude in our application of the method. Mr D. Kovfics attended for the application of Leenen's work. The support of Mr O. Starosolsky, Director General of the Water Resources Research Centre (VITUKI, Budapest) and Mr Gy. Jakus, Director of the North-Transdanubian Water Authority (I~DUVIZIG, Gy6r) is acknowledged, in providing access to their data sets. We are grateful for the constructive comments of the two reviewers of this paper: Mr M.A. Sophocleous and an anonymous referee.

Appendix A This appendix contains the most important steps of the derivations of the one-dimensional groundwater flow equations. Appendix A. 1 Semi-confined condition

For an aquifer covered by a semi-pervious layer and bordered by a river with no river bed resistance [Fig. l(c)], the heads can be calculated from the one-dimensional flow equation, as a function of the river stage and the distance from the river in a

123

section perpendicular to the river (Leenen, 1994): Se 04,

4, --

T Ot

024,=

Tecl

Ox2

0

(A1)

and the boundary conditions, assuming linear rise of the river stage:

x=0

4,=Hpt

(t-
tp

04,

X -----~oe :

----*0

0x

where Se is the elastic storage coefficient or storativity of the aquifer, T is the aquifer transmissivity, 4, is the head in the aquifer, t is the time, tp is the duration of the rising phase of the river stage, c'cl is the resistance of the cover layer, x is the distance from the river, Hp is the full rise of the river stage. The resistance of the cover layer can be calculated from the thickness of the cover layer (dcl) and its hydraulic conductivity (Kz,d) by: dcl

Cc~=

(12)

Kz, cl

It can be proved by a dimensional analysis of Eq. (1) that the process is determined by the following characteristic time (tu) and length (l.) units: tu = Secd

(A3)

lu = ~

(14)

It is possible to rewrite Eq. (1) in a dimensionless form taking: X , = - -X; lu

t * = --" t tp'

4,*=

4, Hp

and defining: SeCcl

tu

tp

tp

Eq. (1) in dimensionless form is:

~ °Ot**+ 4,. 024,*

- Ox~=O

The boundary conditions: x.=0: 4,.=t. (t. ~< 1) X, --* ac :

04,. 0x,

-----~

0

(15)

124

z. Vekerdv, A.M.J. Meijerink/Journal o]"Hydrolog?. ' 205 (1998) 112-125

L e e n e n (1994) gives an approximate analytical solution o f Eq. (5): 1

0 , = < t , - ~/3x,) e

where b is a constant and x u= - -

(A10)

(16) a n d J ( u ) has to satisfy the following equation:

Appendix A.2 Phreatic conditions

f " ( u ) + 2 u f ' ( u ) - 2 n f (u) : 0

A s s u m i n g a vertical river bank, an initial horizontal g r o u n d w a t e r table and no surface recharge, the o n e - d i m e n s i o n a l g r o u n d w a t e r flow in an u n c o n f i n e d aquifer with a water body at one o f its boundaries can be described by:

E d e l m a n (1947) gives solutions for the cases n = 0,1,2 and 3. Case n = 2 describes a linearly rising head in the river (H=c~t, so o~ = H/t):

00 at

T 0z4' = 0 ix Ox2

(A7)

4,=0

(A12)

w h e r e f z ( u ) is a function o f u. E d e l m a n (1947) gives the analytical form off2(u): .f2 (u) = ~ u e

and the b o u n d a r y conditions are: t:0;

4~= o~t~(u)

(0<-x<-~)

q~=0

- u2 - (2u 2 + 1)( 1 - E(u))

(A 13)

e - " du

(A14)

where E(u)=

x:o<

(A 1 1)

0

(0<-t-<~)

E d e l m a n (1947) provides the corresponding n u m e r i cal values o f u and fz(u) in tabular form too. x=0;

0=Hp

(0-t-
where T is the aquifer transmissivity, 4' is the head in the aquifer (the datum is the head at t = 0), t is the time, x is the distance from the river b a n k and IX is the specific yield o f the aquifer. The Dupuit assumptions must be met: 1. the equipotential lines are vertical; 2. the slope o f the g r o u n d w a t e r table is equivalent to the hydraulic gradient and is invariant with depth. E d e l m a n (1947) gave analytical solutions for Eq. (7) to calculate the head rise in a phreatic aquifer. The most important steps are as follows. He introduced: T ~-= --t IX

(A8)

Substituting ~- into Eq. (7): a2~

a~

0x 2

OT

- 0

Solutions o f the following form satisfy Eq. (9): C~: bTn/2/'(U)

(A9)

References

Edelman, J.H., 1947. Over de berekening van grondwaterstroomingen (About the calculation of groundwater flow). Ph.D. Thesis, Delft University of Technology, The Netherlands. Erd61yi, M., I971. A Nyugat-Dun~.ntfil~s a Kisalf61d vizf61dtana (Hydrogeology of Western Transdanubia and the KisalF61d), Hidrol6giai K6zliSny, 11,485 490. Glover, R.E. and Balmer, G.G., 1954. River depletion resulting from pumping a well near a river, Trans. Am. Geophys. Union, 35 (3), 468 470. Hantush, M.S., 1965. Wells near streams with semipervious beds, J. Geophys. Res., 70 (12), 2829 2838. Heij, G.J., 1989. River groundwater relationships in the lower parts of The Netherlands, J. Hydrol., 108, 35-62. Hunt, B., 1990. An approximation of the bank storage effect, Water Resource Res., 26, 2769 2775. Leenen, J.D., 1994. Voortplanting van grondwaterstandsfluctuaties, met het effect van een hoogwatergolf als bi]zondere toepassing (Propagation of groundwater fluctuations specially applied to the effect of a high water wave), H20, 27 (5), 128 131. Maginecz, J., 1993. A Szigetk6z talajvizviszonyainak alakulfisa (The groundwater regime of the Szigetk6z). VITUKI research report, Budapest, Hungary. Moench, A.F., Sauer, V.B. and Jennings, E.M., 1974. Modification of routed streamflow by channel loss and base flow, Water Resource Res., 10, 963 968.

Z Vekerdy, A.M.J. Me(jerink/Journal (~fHydrology 205 (1998) 112-125 Papp, F., 1988. Szivfirgfisi terepkis~rlet eredm~nyeinek 8rt~kel6se (Evaluation of the results of an infiltration measurement). Manuscript, VIZITERV, Budapest, Hungary. Sophocleous, M.A., 1991. Stream-floodwave propagation through the Great Bent alluvial aquifer, Kansas: Field measurements and numerical simulations, J. Hydrol., 124, 207 228. Swain, E.D., 1994. Implementation and use of direct-flow connections in a coupled ground-water and surface-water model, Ground Water, 32 (1), 139 144. Sz6kely, F., 1977. A Dunakiliti vizl+pcs6 hatfisa a Szigetk~3z talajv izviszonyaira (The effect of the Dunakiliti dam on the groundwater regime of the Szigetk6z). V1TUKI research report, Budapest, Hungary.

125

tJjfaludy, L., 1986. Terepkis6rletek a Szigetk6z szivfirgfisi viszonyainak feltfir~.sa c~ljab61 (In situ experiments for the determination of the infiltration characteristics of the Szigetk6z). VITUKI research report, Budapest, Hungary. Vasi[iev, O.F., 1987. System modelling of the interaction between surface and ground waters in problems of hydrology, Hydrol. Sci. J., 32 (3), 297-311. Vekerdy, Z., 1996. Geographical information system based hydrological modelling of alluvial regions using the example of the Kisalf'old, Hungary. ITC publ. 42, ITC, Enschede, The Netherlands. Zhang, W.Z., 1992. Transient groundwater flow in an aquifer aquitard system in response to water level changes in rivers or canals, J. Hydrol., 133,233 257.