A deterministic approach to the coupled analysis of karst springs hydrographs and chemographs

During the chemically based recession flow phase of karstic springs the carbonate (dissolved limestone) concentration can be expressed as negative power of the flow rate. The empirically determined Conc/Q relationship allows two parameters (alpha and A) to be defined, of which one (alpha) depends on the geometric dimensions of the saturated (submerged) karstic network. In this paper we present a deterministic model which simulates the concentration of carbonate at the outlet of a network of circular rectilinear conduits as a function of flow rate. This model, based on hydraulic principles and the calcite dissolution kinetics, allows the sensitivity of the alpha and A parameters to be studied under different chemical, physical and geometric scenarios.