Premium
Low‐dimensional modelling of numerical groundwater flow
Author(s) -
Vermeulen P. T. M.,
Heemink A. W.,
te Stroet C. B. M.
Publication year - 2004
Publication title -
hydrological processes
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 1.222
H-Index - 161
eISSN - 1099-1085
pISSN - 0885-6087
DOI - 10.1002/hyp.1424
Subject(s) - covariance , dimension (graph theory) , mathematics , computer science , flow (mathematics) , eigenvalues and eigenvectors , state space , scale (ratio) , reduction (mathematics) , groundwater flow , covariance matrix , grid , matrix (chemical analysis) , mathematical optimization , algorithm , groundwater , geology , geometry , statistics , physics , quantum mechanics , aquifer , pure mathematics , composite material , geotechnical engineering , materials science
Numerical models are often used for simulating groundwater flow. Written in state‐space form, the dimension of these models is of the order of the number of grid cells used and can be very high (more than a million). As a result, these models are computationally very demanding, especially if many different scenarios have to be simulated. In this paper we introduce a model reduction approach to develop an approximate model with a significantly reduced dimension. The reduction method is based upon several simulations of the large‐scale numerical model. By computing the covariance matrix of the model results, we obtain insight into the variability of the model behaviour. Moreover, by selecting the leading eigenvectors of this covariance matrix, we obtain the spatial patterns that represent the directions in state space where the model variability is dominant. These patterns are also called empirical orthogonal functions. We can project the original numerical model onto those dominant spatial patterns. The result is a low‐dimensional model that is still able to reproduce the dominant model behaviour. Copyright © 2004 John Wiley & Sons, Ltd.