Abstract: A hydrochemical study employing modelling techniques, was carried out using samples taken at 65 points (springs and wells) in Triassic carbonate aquifers (Lújar-Gádor Unit, Alpujárride Complex, Betic Cordillera). These aquifers are made up of limestones and dolomites with some gypsum scattered or interbedded. Though the area is semi-arid, recharge is relatively high because of their mountainous nature. The carbonate rocks contain dense microfissuration; the groundwater flow regime is predominantly diffuse. The karstic forms are in general poorly developed. Two main hydrochemical processes have been identified in these aquifers. One is incongruent dissolution of dolomite that determines the chemical composition of the less mineralised water. The other is dedolomitisation (dolomite dissolution together with calcite precipitation caused by dissolution of gypsum), which becomes predominant when the flow encounters interbedded gypsum. This reaction is also frequently associated with low temperature thermalism, and can play a part in more intense local karstification (cavities, sinkholes, high transmisivity in wells) observed in the sectors of these aquifers where gypsum is more abundant. A reaction path model has been used to simulate the geochemical processes through a hypothetical aquifer (with similar lithology to the Alpujárride carbonate aquifers). Successive stages of evolution through the carbonate sequence, represented by different saturation states with respect to calcite, dolomite gypsum and CO2, have been modelled and then compared with the field data.