STEADY GROUNDWATER-FLOW AS A DYNAMICAL SYSTEM

Dynamical systems describe the time evolution of moving spatial points under the influence of a smooth, bounded vector field. The theory of these systems is focused on the global properties of their flow paths, not on their integration, and so gives general, qualitative informati on that does not depend on the details of the spatial variability of t he vector field. This approach was applied to describe the global cons equences of the Darcy law for steady groundwater flows in isotropic, h eterogeneous aquifers. No particular model of the spatial variability of the hydraulic conductivity (K) was assumed. Vorticity in the flow p aths was shown to exist wherever isoconductivity and equipotential sur faces intersect transversely, and the importance of the Lamb vector (t he vector product of vorticity and specific discharge) for the geometr y of flow paths was established. Because steady groundwater flows gove rned by the Darcy law have zero helicity, they cannot exhibit tangled vorticity lines or become chaotic. The absence of chaos is related clo sely to the impossibility of closed flow paths, the asymptotic stabili ty of isolated minima of the hydraulic head (H), and the existence of a function H(K, H) on whose level surfaces all flow paths are confined . This last function also permits groundwater flows to be represented by moving points in the K, H plane, with motion there generated by a f orm of Hamilton's equations. The results obtained are not related to a ny stochastic approach to aquifer spatial variability, but instead may be applied to constrain stochastic models on purely dynamical grounds .