We study the large-scale evolution and flow in a fractured karst aquifer by means of a newly developed numerical method. A karst aquifer is discretized into a set of irregularly spaced nodal points, which are connected to their set of natural neighbors to simulate a network of interconnected conduits in two dimensions. The conduits are allowed to enlarge by solutional widening. The geometric flexibility of this method, along with a simplified model for the dissolution kinetics within the system water-carbon dioxide-calcite, enables us to study both laminar and turbulent flow in a karst aquifer during its early phase of evolution. A sensitivity analysis is conducted for parameters such as conduit diameter, hydraulic pressure differences, and recharge conditions along the surface of the aquifer and shows that passage evolution depends strongly on the recharge condition and the amount of water available. Under fixed hydraulic head boundary conditions an early single-passage system develops under laminar conditions and is transformed into a maze-like passage system after the onset of turbulence. Fixed recharge boundary conditions are more likely to result in a branchwork-like passage system, although the addition of distributed recharge may lead to a maze-like system of secondary passages. ¿ 1999 American Geophysical Union