Saltwater intrusion of the Costa de Hermosillo aquifer, Sonora, Mexico: A numerical simulation
Main Article Content
Abstract
Inversion of gravity data enabled us to obtain the basement topography of the Costa de Hermosillo aquifer. The depth to the basement ranges from 300 to 3500 m, and has a complex character. Between Hermosillo and Kino Bay, the basement constitutes a subcircular structural high. Around this structural high basement is characterized by a series of horsts and grabens. Several grabens are interconnected, allowing the flow of seawater inland. The horsts act as barriers to the seawater flow. The pattern of total dissolved solids in groundwater supports this interpretation. Gravity data also suggest the existence of tectonic lineaments correlating with the southward prolongation of the Cerro Prieto fault, and with the Craton margin lineament. From electrical resistivity distribution the seawater front extends inland more than 25 km and is more than 200 m deep. Based on borehole information and on the electrical resistivity distribution we inferred the stratigraphy of the aquifer. Geophysical, geological and geohydrological information were used to derive a geological model for the Costa de Hermosillo aquifer, which was subsequently used in a 2-D numerical simulation of salt water intrusion. A non-dimensional formulation based in the equations of Darcy, of continuity (for the water and for the salt, respectively) and in a constitutive relation (fluid density and salt concentration), was used to simulate the seawater intrusion. The resulting equations were solved, through a finite difference scheme, for the parametrized geologic model along selected profiles.
We find an advective pattern where no barrier exists to the seawater flow. However, those sections separated by intrusives or horsts show convective cells. The aquifer flow is complicated by the changing nature of the sedimentary infill and the overexploitation of the aquifer in the zone. Assuming the present rate of withdrawal the steady-state simulated reproduces salt concentration agrees with the observed values.
Article Details
This work is licensed under a Creative Commons Attribution-NonCommercial-ShareAlike 4.0 International License.
References
ABDOH, A., D. COWAN and M. PILKINGTON, 1990. 3D gravity inversion of the Cheshire basin. Geophys. Prospect. 38, 399-1011. DOI: https://doi.org/10.1111/j.1365-2478.1990.tb01887.x
ALENCASTER, G., 1961. Estratigrafía del Triásico superior de la parte central del estado de Sonora. Paleontol. Mex., 1-18.
ANDERSON, T. H. and L. T. SILVER, 1977. U-Pb isotope ages of granitic plutons near Cananea, Sonora. Econ. Geol., 72, 827-836. DOI: https://doi.org/10.2113/gsecongeo.72.5.827
ANDERSON, T. H. and L. T. SILVER, 1981. An overview of Precambriam rocks in Sonora. Rev. Inst. Geología, UNAM, 5, 131-139.
ANDREWS, R. W., 1981. Salt-water intrusion in the Costa de Hermosillo, Mexico: A numerical analysis of water management proposals. Ground Water, 19, 635-647 . DOI: https://doi.org/10.1111/j.1745-6584.1981.tb03518.x
ARREGUIN, J., G. FIGUEROA and S. PEÑA, 1968. Estudio hidrogeológico completo de los acuíferos de la Costa de Hermosillo, Sonora. Ariel Consultores y Dirección de Aguas Subterráneas. Secretaría de Recursos Hidráulicos, México. vol. 1, 1-18.
BEAR, J. and G. DAGAN, 1964. Moving interface in coastal aquifers. Hydr. Div. ASCE., 106, 193-200. DOI: https://doi.org/10.1061/JYCEAJ.0001065
BEAUVAIS, L. and T. E. STUMP, 1976. Corals, molluses and paleogeography of late Jurassic strata of the Poso Serna region, Sonora, Mexico. Palaeogeogr., Palaeo-clim., Palaeoecol., 19, 275-301. DOI: https://doi.org/10.1016/0031-0182(76)90030-4
CRUICKSHANK, C. and R. CHAVEZ-GUILLEN, 1969. Reporte. Ingeniería Hidráulica en México, XXXIII: 31-42.
FERRASI, M. and A. MARINELLI, 1996. An extended formulation of the integrated finite difference method for groundwater flow and transport. J. Hydrol., 453-471. DOI: https://doi.org/10.1016/S0022-1694(96)80020-5
FLORES-MARQUEZ, E. L. 1992. Transfert de chaleur et de masse en milieu sédimentaire et fracturé. Modélisation numérique de la convection naturelle autour du site géothermique de Soultz (Graben du Rhin). Ph. D: Institut National Polytechnique de Lorraine, 229 pp.
FLORES-MARQUEZ, E. L., 1998. Two Dimensional Solute transport model in an Anisotropic and Heterogeneous Porous Medium. In preparation.
DE LA FUENTE, D. M., C. L. AIKEN, J. M. MENA and R. W. SIMPSON, 1994. Cartas gravimétricas de la República Mexicana. Instituto de Geofísica, UNAM, México D.F.
GASTIL, R. G. and D. KRUMMENACHER, 1977. Reconnaissance geology of coastal Sonora between Puerto Lobos and Bahía Kino. Geol. Soc. Amer. Bull., 88, 189-198. DOI: https://doi.org/10.1130/0016-7606(1977)88<189:RGOCSB>2.0.CO;2
GHASSEMI, F., A. J. JAKEMAN, G. JACOBSON and K. F. HOWARD, 1996. Simulation of seawater intrusion with 2D and 3D models: Nauru Island case study. Hydrogeol. J., 4, 4-22. DOI: https://doi.org/10.1007/s100400050251
GIPOULOUX, O., J.-P. CALTAGIRONE and P. MOREL, 1991. Approche numérique du comportement de la loi de Darcy en fonction du nombre de Reynolds. Comptes Rendus de l’Académie des Sciences (Paris), Série II, 313: 1493-1498.
GOMEZ, J. 1971. Sobre la presencia de estratos marinos del Mioceno en el estado de Sonora, México. Rev. Inst. Mexicano del Petróleo, 3, 77-78.
GRANT, F. S. and G. F. WEST, 1965. Interpretation theory in applied geophysics. New York, 584 pp.
HENRY, H. R., 1959. Salt Intrusion into Freshwater Aquifers. J. Geophys. Res., 64, 1911-1959. DOI: https://doi.org/10.1029/JZ064i011p01911
HENRY, H. R., 1964. Effect of dispersion on salt encroachment in coastal aquifers. In: H. H. C. et al. (Editors), Sea Water in Coastal Aquifers. US Geological Survey Water, Supply Paper, pp. 1613-C.
KING, P. B., 1939. Geologic Reconnaissance in northern Sierra Madre Occidental of Mexico. Geol. Soc. Amer. Bull., 50, 1624-1722. DOI: https://doi.org/10.1130/GSAB-50-1625
LEE, C.-H. and R. T.-S. CHENG, 1974. On Seawater Encroachment in Coastal Aquifers. Water Resour. Res., 10, 1039-1043. DOI: https://doi.org/10.1029/WR010i005p01039
MERRIAM, R., 1972. Reconnaisance geologic map of the Sonoyta quadrangle, northwest Sonora, Mexico. Geol. Soc. Amer. Bull., 83, 3523-3536. DOI: https://doi.org/10.1130/0016-7606(1972)83[3533:RGMOTS]2.0.CO;2
MUALEM, Y. and J. BEAR, 1974. The shape of the interface in steady flow in a stratified aquifer. Water Resour. Res., 10, 1207. DOI: https://doi.org/10.1029/WR010i006p01207
O’NEILL, D. J., 1975. Improved linear filter coefficients for application, in apparent resistivity computation. Bull. Aust. Soc. Explor. Geophys., 6, 104-109. DOI: https://doi.org/10.1071/EG975104
PARKER, R. L., 1973. The rapid calculation of potential anomalies. Geophys. J. Royal Astron. Soc., 31, 447-455. DOI: https://doi.org/10.1111/j.1365-246X.1973.tb06513.x
PILKINGTON, M. and D. J. CROSSLEY, 1986. Determination of crustal interface topography from potential fields. Geophysics, 51, 1277-1284. DOI: https://doi.org/10.1190/1.1442180
PINDER, G. F. and H. H. A. COOPER, 1970. A numerical technique for calculating the transient position of the salt water front. Water Resour. Res., 6, 875. DOI: https://doi.org/10.1029/WR006i003p00875
RANGIN, C. 1978. Consideraciones sobre la evolución geológica de la parte septentrional del estado de Sonora. Libreto Guía del Primer Simposio Sobre la Geologíay Potencial Minero del estado de Sonora. Inst. de Geología, UNAM, 35-56.
RANGIN, C., 1981. Late Triassic to Late Cretaceous geological evolution of northwestern Mexico: 77th annual international meeting of the Geological Society of America, Geol. Sur. of America.
SEGOL, G. and A. G. F. PINDER, 1976. Transient simulation of saltwater intrusion in south-eastern Florida. Water Resour. Res., 12, 65. DOI: https://doi.org/10.1029/WR012i001p00065
SEGOL, G., G. F. PINDER and W. G. GRAY, 1975. A Galerkin finite element technique for calculating the transient position of the saltwater front. Water Resour. Res., 11, 343. DOI: https://doi.org/10.1029/WR011i002p00343
SHAMIR, U. and G. DAGAN, 1971. Motion of sea water interface in coastal aquifers: A numerical solution. Water Resour. Res., 7, 644. DOI: https://doi.org/10.1029/WR007i003p00644
STEINICH, B., I. SIMON, J. A. CHAVARRIA and L. E. MARIN, 1997. Geophysical investigations of the vadose zone in the Valley of Hermosillo aquifer, Sonora, Mexico. Geofísica Internacional, 36, 191-200. DOI: https://doi.org/10.22201/igeof.00167169p.1997.36.3.650
WOOLLARD, G. P., 1979. The new gravity system changes international gravity base values and anomalies. Geophysics, 44, 1352-1366. DOI: https://doi.org/10.1190/1.1441012