Abstract

In order to understand the geological and hydrogeological processes influencing the hydrogeochemical behavior of the Cerro Prieto Geothermal Field (CP) aquifer, Mexico, a characterization of the water samples collected from geothermal wells was carried out. Different hydrochemical diagrams were used to evaluate brine evolution of the aquifer. To determine pressure conditions at depth, a calculation was performed using hydrostatic and lithostatic properties from CP, considering geological characteristics of the study area, and theoretical information about some basin environments. Groundwater shows hydrogeochemical and geological evidences of the diagenetic evolution that indicate overpressure conditions. Some physical, chemical, textural, and mineralogical properties reported in the lithological column from CP are explained understanding the evolutionary process of the sedimentary material that composes the aquifer.

Dedicated in memory of Dr. Ramiro Rodriguez C.

1. Introduction

Different authors have observed that the origin of some brines around the world can be caused by diagenesis evolution and porewater trapped during burial [13]. Commonly, porewater shows characteristics according to the depositional environment [4]. The different stages of diagenesis are composed of burial and compaction of material in sedimentary basins; the sandstones and shales undergo changes in their physical, chemical, textural, and mineralogical properties which are reflected in the sediment density, the compaction of the granular package, and loss of porosity [5]; additionally, a chemical alteration occurs in sandstones generating cementation and lithification as a product of chemical precipitation affecting detrital grains, dissolution, recrystallization, or mineralogical alteration [6]. Silica is the most abundant cementing agent in sandstones, compared to calcium carbonate, either calcite or aragonite, since the latter dissolves more easily in contact with groundwater, mainly the second [7]. In the diagenetic stages, interstitial fluids are constantly lost in most shales [8]; this phenomenon can occur at low depth and can last for thousands or hundreds of thousands of years or have a longer duration and greater depth; but when it suddenly occurs and the permeability of the medium is so low that it does not allow the interstitial fluid to leak, an increasing degree of stress can be generated [9]; this geological event is known as overpressure [10]. Overpressure is an abnormally high pressure in the subsoil that exceeds the hydrostatic pressure at a certain depth; this pressure is carried out in the pores where the pressure of the interstitial fluids increases as the overcoat increases [11, 12]. Overpressure indicates that the high pressures developed during compaction do not dissipate efficiently [10] and can generate hydraulic fracturing in the system (stress applied to compressible rock and fluid expansion). These processes generate overpressured fluids related to porosity reduction, changes on porewater flow, and diagenetic reactions due to compaction and disequilibrium in a sedimentary basin [1214]. Also overpressure can be related to chemical compaction due to changes in mineralogy (e.g., ion exchange, dissolution/precipitation) or to diagenetic processes and fluid expansion by thermal water in pore space.

Cerro Prieto Geothermal Field (CP), located in northwestern Mexico (32°24′43′′N, 115°14′41′′W), is a brine with high-temperature geothermal system characteristics. Several studies about the origin and behavior of CP groundwater have been reported [1518]. According to geological evidences, a large accumulation of sedimentary material from a continental and marine origin, overlying the depositional basin, is related to the origin of brine [17]. The sedimentary material shows diagenetic evolution evidences and porewater trapped between sediment grains during burial processes. The porewater is saline with high Cl, Na+, Ca2+, and K+ concentration; in geothermal brines this characteristic is common [1, 2, 10]; the composition depends mainly on the primary origin, mineralogical composition of the sediments and their modification due to diagenetic processes (e.g., facies distribution), and hydrothermal characteristics [15, 19]. Among the most evident diagenetic processes in CP are cementation, mineral replacement, recrystallization, authigenesis, and growth of concretions and nodules [16]. On the other hand, in deep sedimentary basins as CP, mechanical processes of deformation related to burial mechanisms are common; likewise, hydrostatic and lithostatic conditions increase with depth due to an increase of the superposed fluids hydraulically connected, through the pore and the pressure exerted by sediments overload [10]. If pore pressure in deep aquifers like CP is higher than expected from hydrostatic conditions, anomalous pressure (overpressure) can be generated; overpressure is common mainly within 2–4.5 km depth [1012].

According to geochemical evidences, the origin of geothermal brine at CP could be governed by mixing processes related to a hydrothermal environment and the sedimentary material located at depth which shows burial diagenesis evolution with hydrogeochemical evidences of an overpressurized environment. The aim of the present study was to evaluate hydrogeochemical behavior of geothermal groundwater and its relation with diagenetic processes, including overpressure.

2. Localization

Cerro Prieto Geothermal Field (CP), located in Mexicali Valley, SE of Mexicali City, in Baja California State, Mexico (Figure 1), is a Basin of Salton Sea [17]. Climate is arid with temperatures up to 40°C in July and to 4°C in winter. The average annual precipitation is 55 mm/year and the average annual evaporation is 2200 mm/year [23]. Groundwater at CP is extracted from geothermal wells that are in constant exploitation to generate electricity. The Comisión Federal de Electricidad of Mexico (CFE) operates and manages the Geothermal Field. CP power production is up to 720 MW and is composed of five individual units: CP1, CP2, CP3, CP4, and CP5; each unit has a total capacity of production with a specific number of production wells. Sampled wells are indicated in Figure 1. Their localization in the five individual CP units is included in Table 1. All the wells are located in zone “beta,” 1500 to 3100 m depth [17].

2.1. CP Geology

The lithology around CP is composed of gneiss (quartz-feldspars), shale (quartz-mica), marble, amphibolite, and quartzite from Permic to Jurassic [24] and metamorphic, granitic, and granodiorite rocks which are intruded by batholitic rocks [16], together with dacite and andesite from Miocene and rhyodacite from Quaternary [25]. The tectonic basin was filled by sedimentary material that, due to burial, compaction, and diagenesis processes, evolved to gray shales from Late Miocene (shales and silt shales that vary from light gray to black); this unit overlies the granitic basement and the mafic intrusive and is interlayered by permeable sandstones (composed of quartz and feldspars; arkoses type). The thickness is near 3000 m [21, 26]. Immediately above a layer of brown shale (shale and silt shale), the gray shale covers interlayered permeable sandstones and sands cemented by carbonates, about 500 m thickness [21]. In these zones, a rapid distribution of geothermal fluids enhances the recharge. Erratic mudstone layer and unconsolidated clastic sediments (clays, silts, sand, and little gravel) overlie all the previous units. The thickness of these units is between 400 and 2500 m [26]. Sedimentary material shows evidence of the diagenetic evolution and recrystallization processes due to an incipient low grade metamorphism.

The geological evolution of CP is a complex blend of rifting, rapid deltaic sedimentation, and large scale strike-slip faulting located within the Salton Basin [16, 19, 21, 27]. The Geothermal Field is placed in a shear zone where NW-SE and NE-SW fault systems intersect. The more important faults are Cucapá, Imperial, Cerro Prieto, and Michoacán [21]. This fault system is part of a major regional lineament that penetrates deep into the crustal and basement rocks and serves as conduit for geothermal flow. The system originates in a tectonic basin of 5200 m depth, filled by alluvial and deltaic sediments from Tertiary to Quaternary [21].

Vonder Haar and Howard [27] observed that, in sandstone and shale units, a mineral dissolution/precipitation took place along microfractures, originating secondary porosity and newly precipitated hydrothermal minerals, causing a reduction of permeability. Likewise, Elders et al. [20] observed cementation, mineral replacement, recrystallization, authigenesis, and growth of concretions and nodules; these processes are related to diagenesis.

2.2. CP Hydrogeology

Some authors consider that CP brine may have been formed from marine evaporite dissolved and partly by evaporated Colorado River water [16, 19, 21, 2831]. However, according to geological evidences, a large accumulation of sedimentary material, overlying the depositional basin, from a continental and marine origin and mixing with meteoric water [3033], was related to the origin of brine [17]. The sedimentary material shows diagenetic evolution evidences and porewater trapped between grains during burial processes. Isotopic evaluations (18O, 2H) and chemical analysis (Cl and Br) elaborated by Coplen [34] and Birkle et al. [35] suggest that Salton Sea was the probable predecessor of high chlorinated groundwater of CP.

3. Methodology

A hydrogeological and hydrogeochemical study was carried out in geothermal groundwater samples from CP following standard methods of APHA-AWWA [38]. Water samples were collected from geothermal wells that are in constant exploitation. Temperature, pH, and conductivity were measured in the field during the summer of 2010 and calibrated to the water temperature at each site. The chemical analyses for the major elements, B and SiO2, were performed at the Analytical Chemistry Laboratory of the Geophysics Institute, UNAM, Mexico (the laboratory participates in international calibration exercises of chemical analyzes of geothermal waters). Boron was determined by colorimetry through reactions with carminic acid (Method 4500-B C) APHA-AWWA [38]. SiO2 was determined by atomic absorption spectrophotometry with flame and UV-visible spectroscopy (molybdosilicic acid method). Major ions were analyzed following standard methods [38]. and were determined by volumetry (titrating with HCl), Ca2+ and Mg2+ were determined by volumetry (titrating with EDTA), Cl was determined by potentiometry with selective electrodes (4500-Cl) [38], Na+ and K+ were determined by atomic emission spectrophotometry (3500-Na+ and K+), and was determined by turbidimetry (4500-). Analytical quality was assessed through ionic balance (less than 10%) and the use of certified (NIST) reference solutions.

In order to evaluate brine evolution of CP, different techniques were used: (a) Carpenter [39] evaluated the behavior of major elements using a plot with concentrations as a function of dissolved chloride concentration, considering the composition of seawater during evaporation and diagenesis; using chemical results of CP groundwater a similar evaluation was elaborated; (b) Davisson and Criss [40] devised a diagram to determine the geochemical evolution of mineralogy in brines applying an evaluation of and in water samples; hydrogeochemical results of samples from CP were evaluated with this diagram; (c) Boschetti [3] considers B-Cl concentrations to explain the evolutionary process in the groundwater; B-Cl diagram was used to determine the dominant geological environment in CP. A geochemical simulation with the measured water concentrations was carried out using the Phreeqc© program to determine saturation indexes.

In order to determine pressure conditions, a calculation was performed using hydrostatic and lithostatic properties from CP, considering geological characteristics of the study area and theoretical information about some basin environments. To estimate pressure conditions () exerted by a geological column at a depth () offshore in a geologic formation, (1) was used [12]:It is necessary to take into account pure water column weight at sea level , depth of seawater column , seawater density  kg/m3, gravity constant (m/s2), and material density (rock or sediments) . The integral about overburden weight of sediments can be replaced by a sum of the individual weights of layers [12]:where it is important to consider the number of layers () with thickness ( is the layer number measured in m), rocks density (km/m3), porosity , and water density (which can change with salinity variation, while temperature and pressure dependence is relatively small or negligible).

To estimate pressure in a reservoir unit from CP it was considered that the study area is located onshore, few meters above seawater level; therefore the integral including the weight of seawater column is zero; the distinctive strata units, their thickness, and depth of each sampled well were considered; the porosity values reported by Olson [41] and Hiriart Le Bert [42] in the geologic material from CP were used (between 0.15 and 0.25), for lithological units without reported porosity; theoretical values proposed by distinct authors studied unlike geological material and depth were used (e.g., mudstone and shales, slate, quartz, feldspars, uncemented sandstones, or sandstones as reservoirs) [4346].

4. Results and Discussion

4.1. Major Elements in CP Groundwater

Table 1 shows the major solutes concentration in CP groundwater.

In Table 1 chemical results from the studied wells are shown. Total dissolved solids were calculated from conductivity measures showing values between 12945.5 and 49251.2 mg/L. These values are comprised between saline water (>1,000 mg/L) and brine (>35,000). All the sampled wells are located at beta reservoir [17, 47]. To explain geochemical variations in CP brine it is necessary to evaluate the composition and removal of solutes by salt precipitation according to diagrams proposed by Carpenter [39] where the circle represents the solute-chloride composition of normal seawater; the line represents the limit between evaporation-dilution curves of seawater and freshwater (Figure 2).

Figure 2 shows that groundwater in CP is dominated by high concentrations of Na+, K+, Ca2+, and Cl. When Cl concentration increases Na+, K+, and Ca2+ concentrations also increase in ratios 1 : 1, 1 : 1, and 2 : 1, respectively (Figures 2(a), 2(b), and 2(c)). Bicarbonate ion, Mg2+, and have a low concentration (Figures 2(d), 2(e), and 2(f)), Na+ values are parallel and near coincident with their respective evaporation-dilution curve (Figure 2(a)), Ca2+ values cross the evaporation-dilution curve (Figure 2(b)), and K+ also increases with Cl but the values are enlarged relative to the seawater evaporation-dilution curve (Figure 2(c)). In CP brine, Mg2+ concentrations lie well below the seawater evaporation trajectory indicating significant depletion of the element. According to Hanor [48] and Kharaka and Hanor [1] Mg2+ concentrations in brines decrease when temperature increases in the subsurface, and when alkalinity decreases [49]. Evaporation of continental waters has a more variable concentration range in water samples from CP due to evaporating mixtures of continental and marine water with meteoric water [50].

In CP brine depleted and enriched concentrations of some major elements are a consequence of reactions linked with the hydrothermal processes and water-rock interactions. Major elements concentrations are controlled by the alteration and formation of minerals like feldspar-K, plagioclases, quartz, biotite, amphibole, chlorite, pyrite, wairakite, prehnite, muscovite, epidote, and talc, as reported in CP by Elders et al. [20] and Izquierdo et al. [16] and shown by calculated saturation index values (Figures 3 and 4).

The K+ origin is restrained by alteration of feldspar-K, illite, and biotite and by muscovite formation. Very low and concentrations in CP could be inhibited by the water interactions with anhydrite, dolomite, talc, pyrite, or calcite (Table 1; Figures 2 and 3); Elders et al. [20] report pyrite formation at depth. Low concentration of Mg2+ and high concentration of Ca2+ could be related to dolomitization of limestone as major source of Ca2+; and low Mg2+ contents are associated with the evolution of chlorites and micas when temperature and depth increase according to mineralogy reported in CP (Figure 3); similar behavior has been reported previously in hydrothermal brines with diagenetic evolution evidences [51] (Figure 3). Albite reactions at depth at high temperature can be linked with the slight Na+ decrease [16, 30, 31].

4.2. Saturation Index

Results of saturation index calculations are shown in Figure 4. From these results amorphous silica (SiO2 am), albite, k-feldspar, and in some cases k-mica show a behavior close to equilibrium with the fluid. Besides, quartz, chalcedony, talc, and crysocole are oversaturated. Dolomite, calcite, and aragonite are undersaturated at some sites and oversaturated at others, in agreement with Figures 2 and 3.

4.3. Na(deficit)-Ca(excess) Plot

Figure 5 was used to explain the initial composition of brines and the nature of fluid-rock interactions. In the diagram, Basinal Fluid Line (BFL) is a straight line with a unit slope that indicates a 2Na-1Ca exchange relationship [40]; BFL represents the effect of plagioclase albitization on water composition. Seawater evaporation trajectory is a representation of the natural trends for seawater evaporation which is formed by large positive and small negative ; reactions involving seawater evaporation follow a vertical descent and afterwards produce large deficits along a horizontal trend. Halite dissolution in seawater or freshwater can produce negative values along a slope of 1 : 4.

To determine the origin and geochemical evolution of CP brine, an evaluation of and was applied to explain the initial composition and nature of fluid-rock interaction (Figure 5). All the analyzed water samples from CP were located right and over the seawater evaporation trajectory (SET) in the - diagram (Figure 5) indicating that the fluids are a product of brine that passed the point of halite precipitation evaporated (along of axis). The horizontal line of CP has a large positive and a small negative , but the fluid is more enriched in than expected from seawater evaporation.

Dolomitization produces elevated Ca contents, increasing without changing the , and can explain the phenomenon observed in CP. Other dissolved minerals could interact in the evolutionary processes to brine aside from the above-mentioned ones, such as calcite, anhydrite, quartz, halite, and illitization to smectite transformation. Often, these geochemical processes are treated separately but in some situations the origin is linked to a mixing of processes and is not mutually exclusive [1]. Illite formation involves reactions relevant for diagenesis: (a) release of water during transformation of feldspar to kaolinite and smectite; (b) potassium and silica budget [2, 52]. Some of these processes can occur in CP.

Geology of CP shows equilibrium with clays (Figure 4); groundwater shows high concentrations of K+. Also reactions between calcite, illite, and K+ to form K-feldspar could be involved in brine evolution. According to the results (Figures 2 and 5) the origin of CP groundwater is a consequence of an evolution of dissolved evaporative products (e.g., halite), residual water remaining during the dissolution precipitation of seawater evaporites, and different water-rock interactions (e.g., clays, siltstones, and shales) but could have a slight contribution to exchange reactions 1Ca for 2Na in the aquifer by albitization of plagioclases, which could change the ionic composition of fluids (Figure 5), according to the following reaction proposed by Carpenter [39], Hanor [48], and Demir and Seyler [53]:

4.4. Overpressure
4.4.1. B-Cl Plot

Some hydrogeochemical evidences (Figures 2 and 5) shown in a B-Cl plot (Figure 6) indicate that overpressurized fluids participate in the evolutionary process of the CP geothermal brine. Boron behavior helps to define these geologic environments because B is adsorbed by clay minerals and is released into the fluid in a deep environment, mainly when tectonic stresses by vertical and/or lateral compaction are high and temperature increases with depth, and stronger geologic deformations are generated [10, 12, 54].

4.4.2. Pressure Conditions Estimated due to Overburden

Some authors [55, 56] estimated pressure at depth at CP considering the host rock density at a specific depth. The obtained values were in the range between 0.5 and 42 MPa. Those calculations considered only geothermal water conditions from the extraction zone and no lithologic information was used from the stratigraphic columns of each site. To confirm hydrogeochemical evidence of overpressure, a calculation to determine pressure conditions was performed using hydrostatic and lithostatic properties from CP and theoretical information about some basin environments (Figure 7). It is necessary to consider that normal pressure increases with depth according to the hydrostatic pressure gradient (10 MPa/Km), higher or lower values of this gradient, and their associated depths are considered abnormal pressures (e.g., overpressure). Lithostatic pressure is equivalent to the total charge of the overlaying sediments in a geological formation and increases according to the lithostatic pressure gradient (23 MPa/Km) [1012, 57].

In general, overpressured systems can take place when porewater is not expelled from rock at a proper rate, remaining under hydrostatic pressure. Overpressured volume rock must be trapped by low permeability layers where fluids move slowly, even when there is high pressure in the environment. The overpressure affects the effective stress that acts between the grains within the rock and generates a change in the compaction. In many areas of active sedimentation rate around the world, porewater pressure in deep groundwater (>1 km) is higher than would be expected from hydrostatic circumstances [10, 5759]. It is necessary to consider thermal expansion in the pore space (increasing volume) and increment of the system temperature by thermal conditions and by fluid movement and mineral phases transition [11]. According to Swarbrick et al. [11] and Kauerauf and Hantschel [12] secondary overpressure by chemical cementation may occur at large depths when porosity is reduced by dissolution-diffusive transport-precipitation of silica cement (temperature affects diffusion-precipitation rate) or by fluid expansion processes when gas or thermal solutions are originated in highly permeable facies interconnected locally at certain depth levels generating compaction, rearrangement of grains, and reduction in the pore space. Some of these conditions could occur in CP aquifer [22].

The results obtained confirm that in CP overpressure is present (Figures 7 and 8). Positive anomalies with an increase of depth are shown between 0.6 and 3.1 km; this phenomenon can be caused by a hydraulic seal. Overpressure coincides with lithology variation in CP [22], when sedimentary material composed of sand and clay changes to shale and sandstone due to clay and sand diagenesis and shale compaction (Figure 7). In CP, sandstone is hosted by a low permeability layer of shale and siltstones creating adequate conditions for this process. Overpressure in CP brine could occur as a consequence of a rapid sedimentation, deposition, and accumulation rates of fine-grained material along time and due to the stress increase in sediments, compressibility, or expansion of fluids by hydrothermalism.

Geological evidence in CP indicates that mineralogical changes occur at depth, mainly by diagenetic processes (Figure 5). According to the geochemical behavior observed in CP, it is possible that mineral dissolution precipitation processes generated this phenomenon; osmosis, buoyancy, and tectonic or magmatic process can generate changes in mineralogy (e.g., feldspar to illite or smectite to illite conversion), diagenesis, and carbonate or silicate cementation. In CP, sandstones and shales units reveal clogging, mineral dissolution, and mineral precipitation along microfractures as indicated by Vonder Haar and Howard [27]. In the study area the calcite dissolution and/or cementation more likely develops where cold water interacts with hotter rocks or precipitation of quartz and k-feldspar occurs when hot waters interact with colder rocks [16, 20]. Exchange of Na+ by K+ with a higher ion radius is carried out in mineral transformations and increases the volume of the solid matrix [12]; this process is controlled by temperature and K+ availability in minerals.

The permeability of sandstones in CP facilitates reactions between rocks and hydrothermal fluids (e.g., dissolution of some minerals); these reactions can reduce or increase the porosity and also generate secondary fracturing or microfracturing and modify physical properties [16, 21], which could cause overpressure (Figures 7 and 8). Sandy shale and siltstone facies in CP are most amenable to increased microfracturing. In sandstones (where high temperature dominates) and shales from CP, mineral dissolution precipitation takes place along microfractures, which originates secondary porosity and newly precipitated hydrothermal minerals, causing a reduction of permeability.

5. Geothermometry

Alkali feldspar geothermometers are the most used tool to determine chemical equilibrium in fluids at depth in a geothermal system [60]. Na/K and Na-K-Ca geothermometers were developed to evaluate the temperature in high enthalpy geothermal systems [61]; these geothermometers are less affected by chemical reequilibration in mixing zones, but the calculated temperature may be affected by mixing with cold water or by deposition of aluminum-rich minerals or clays [62]. Na/K geothermometers are adapted for temperatures between 180 and 350°C. Na/K geothermometers have been generally used to estimate temperatures at CP [36, 37, 62]. We applied (4) and (5) [61, 62] to evaluate the temperature of CP reservoir:The temperature range of CP reservoir, calculated with Na/K geothermometer, varies between 236 and 306°C for (4) and between 251 and 318°C for (5) (Table 2).

6. Conclusions

CP shows intermediate brine characteristics (Na-Ca-Cl) with high K+, Ca2+, and salinity contents; the relation Na+/Cl is less than 1.

The porewater composition in CP evolved from its primary origin and was modified by the interaction with minerals of the sedimentary material. Brine characteristics were acquired by deep-burial diagenesis processes and low grade metamorphism at high temperatures. Results show geochemical evidence of overpressured fluids due to compaction.

Groundwater samples from CP show a mixing of marine and continental water; this situation is partially related to a continental and evaporative precursor. The hydrogeochemical evidence indicates that the sedimentary material has porewater between grains, which was trapped during burial processes. The diagenetic processes could have generated high concentrations of Cl, Na+, K+, and Ca2+. Calcium enrichment, Na+ depletion, and K+ release could have a relation with a contribution of exchange reactions 1Ca for 2Na in the aquifer by albitization of plagioclases or by illitization processes, respectively, and precipitation of secondary minerals. High K+ and low Mg2+ contents are related to alteration of feldspar-K, illite, biotite, and muscovite formation. Ca-Na exchange with plagioclases could be a geochemical control on the fluids of CP and may directly explain slight and in the brine.

Overpressure in CP is related to burial mechanisms; secondary overpressure is related to chemical pressure by mineralogical changes and by fluid expansion, which increases with depth. The magnitude of overpressure may be produced by some characteristics of deposit formation (burial), permeability evolution of sedimentary material, and the compressibility of rock and fluid. Secondary overpressure in the system related to chemical pressure and porosity changes due to mineral dissolution can be generated at large depths. Fluid expansion takes place in the reservoir which generates compaction, rearrangement of grains, and reduction of pore space. Boron release at overpressure conditions can be related to high contents of K+ in water.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

The authors thank Aguayo A., Ceniceros N., and Cruz O. for chemical determinations. The authors acknowledge Comision Federal de Electricidad for support on sampling within the Cerro Prieto Geothermal Field.