Modeling transient stream/aquifer interaction with the non-linear Boussinesq equation and its analytical solution

• Published in: Journal of hydrology (Amsterdam) Vol. 206; no. 3; pp. 245 - 255
• Main Authors: Serrano, Sergio E., Workman, S.R.
• Format: Journal Article
• Language: English
• Published: Amsterdam Elsevier B.V 01.05.1998; Elsevier Science
• Subjects: Analytical solution; aquifers; Earth sciences; Earth, ocean, space; equations; Exact sciences and technology; Fluctuating rivers; groundwater; groundwater flow; hydraulic conductivity; Hydrogeology; Hydrology. Hydrogeology; Mathematical models; Non-linear equation; rivers; stream flow; Stream/aquifer interaction; streams; Unconfined aquifers; water table; United States; Unconfined aquifers; Mathematical models; Non-linear equation; Stream/aquifer interaction; Analytical solution; Fluctuating rivers; transient flow; models; canals; streams; mathematical models; ground water; water table; ground-water recharge; North America; transmissivity; hydraulic conductivity; propagation; aquifers; daily average; solution; unconfined aquifers
• ISSN: 0022-1694; 1879-2707
• DOI: 10.1016/S0022-1694(98)00111-5

A new analytical solution of the non-linear transient groundwater flow equation subject to time variable river boundaries was used to simulate stream/aquifer interactions in an alluvial valley aquifer. The differential equations were solved using the method of decomposition. The mathematical model required relatively few parameters to simulate groundwater elevations: hydraulic conductivity, specific yield, and recharge. The model was physically based and could simulate the process of a flood wave propagation into an unconfined aquifer. The model was tested using observed water table elevations at three locations across a 2 km wide alluvial valley aquifer. The average daily deviation between observed and simulated water table elevations was approximately 0.09 m. The transient redistribution of water in the aquifer was simulated well with the model. The non-linear form of the Boussinesq equation was shown to better simulate cases when the transmissivity of the aquifer could not be assumed to be a constant.