The nested superblock approach for regional-scale analytic element models

A new approach is presented for improving the computational efficiency of regional-scale ground water models based on the analytic element method (AEM). The algorithm is an extension of the existing "superblock" algorithm, which combines the effects of multiple analytic elements into Laurent series and Taylor series (superblock expansions). With the new "nested superblock" formulation, Laurent series are nested in a hierarchical (quad-tree) data structure with direct mathematical relationships between parent and child superblock coefficients. Nested superblocks significantly accelerate the evaluation of the complex potential and discharge function in models that contain a large number of analytic elements at multiple scales. This evaluation process, the primary computational cost of AEM models, is required to determine the element coefficients, generate contour plots, and trace pathlines. The performance of the nested superblocks is demonstrated with a simplified model based on the Lake Ontario watershed geometry comprising thousands of hydrogeologic features at multiple geographic scales.