Analytic Element Modeling (AEM) Software Compared


  • Comprehensive transient flow: Most things that can be modeled in steady state can also be modeled in transient simulations. Important parameters such as specified heads, discharges, and recharge rates can vary with time during the simulation.
  • Multi-level models: Simulation of multiple levels with vertical leakage and resistance between levels, as opposed to single-layer (2D) simulation.
  • Layer scheme transitions: Different layering schemes in different regions of the model. For example, a single layer in one area and four layers in the area of interest.
  • Horizontal anisotropy: K anisotropy in the plane of the layer, and the angle and ratio can be heterogeneous within a layer.
  • Fresh/salt interface aquifers: Aquifers with a sharp fresh-salt interface using the Ghyben-Herzberg equation.
  • Multi-level wells: Wells that draw from multiple layers in a multi-layer model, but have a consistent head.
  • Uniform area sinks – circular: Uniform rate of recharge applied within a circular area.
  • Uniform area sinks – irregular polygon: Uniform rate of recharge applied within a polygon area.
  • Spatially-variable leakage/storage: Smoothly-varying recharge, leakage, (or storage in transient models).
  • Top/bottom head-dependent leakage: Recharge or leakage in/out the top or bottom of the model based on head difference from a specified head.
  • Heterogeneity (interdomain) line elements: Line elements across which the aquifer definition changes.
  • Circular/elliptical heterogeneities: Same as above but for circle and/or ellipse-shaped boundaries.
  • Head-specified line elements: Elements that extract at a rate to meet specified head condition(s) on the line.
  • Normal flux-specified line elements: Elements that create a specified discharge normal to the line (e.g. impermeable).
  • Head-dependent discharge (river) line elements: Internal elements where the discharge per length is proportional to the difference between a specified stage and the underlying aquifer head, and proportional to a specified resistance.
  • Head-dependent normal discharge line elements: Like the river elements described above, but for an external lateral boundary. Like the GHB boundary in MODFLOW.
  • Discharge-specified line elements: Internal elements where the discharge/length or total discharge is specified.
  • Conductive drain/fracture line elements: Internal elements where the flux carrierd by the line is proportional to a specified conductance and the head gradient along the line.
  • Resistant barrier line elements: Internal elements where the flux normal to the line is proportional to a specified conductance and the head difference across the line.
  • High order line elements: Each line segment may have at least 3 degrees of freedom (parameters), the number of which can be adjusted to improve accuracy.
  • Conjunctive surface water – groundwater: Solutions that also solve for river or lake stage, based on leakage rates between the surface water and aquifer.
  • Graphic user interface: User interface with integrated menus, data entry forms, graphics, etc. Programs without a GUI are run from text commands.
  • Detailed tools to analyze solution: Solution may be examined in detail to check accuracy of solution along boundaries, discharges along elements, hydrographs, head profiles, and more.

Software Notes:

  • Visual AEM is a user interface for several AEM computational engines, including Bluebird, Cardinal, Split, and dated versions of TimML. The current version of TimML is not supported in Visual AEM, so TimML capabilities are listed separately.
  • Winflow is the AEM portion of Aquifer Win32, which is a user interface for several programs.

Sources other than AnAqSim (accessed July 24-26, 2016):

  • GFLOW:
  • TimML:
  • TTIM:
  • Visual AEM:
  • ModAEM:,
  • WhAEM:
  • WinFlow: Demo download with documentation available at