Assessment of groundwater flow in volcanic faulted areas. A study case in Queretaro, Mexico

array:22 ["pii" => "S0016716915000197" "issn" => "00167169" "doi" => "10.1016/j.gi.2015.04.016" "estado" => "S300" "fechaPublicacion" => "2015-07-01" "aid" => "16" "copyrightAnyo" => "2015" "documento" => "article" "crossmark" => 0 "licencia" => "http://creativecommons.org/licenses/by-nc-nd/4.0/" "subdocumento" => "fla" "cita" => "Geofisica Internacional. 2015;54:199-220" "abierto" => array:3 [ "ES" => true "ES2" => true "LATM" => true ] "gratuito" => true "lecturas" => array:2 [ "total" => 2603 "formatos" => array:3 [ "EPUB" => 109 "HTML" => 1763 "PDF" => 731 ] ] "itemSiguiente" => array:18 [ "pii" => "S0016716915000215" "issn" => "00167169" "doi" => "10.1016/j.gi.2015.04.018" "estado" => "S300" "fechaPublicacion" => "2015-07-01" "aid" => "18" "documento" => "article" "crossmark" => 0 "licencia" => "http://creativecommons.org/licenses/by-nc-nd/4.0/" "subdocumento" => "fla" "cita" => "Geofisica Internacional. 2015;54:221-44" "abierto" => array:3 [ "ES" => true "ES2" => true "LATM" => true ] "gratuito" => true "lecturas" => array:2 [ "total" => 3046 "formatos" => array:3 [ "EPUB" => 97 "HTML" => 2236 "PDF" => 713 ] ] "en" => array:11 [ "idiomaDefecto" => true "titulo" => "Mesozoic gliding and Tertiary basin and range tectonics in eastern Sonora, Mexico" "tienePdf" => "en" "tieneTextoCompleto" => "en" "tieneResumen" => array:2 [ 0 => "es" 1 => "en" ] "paginas" => array:1 [ 0 => array:2 [ "paginaInicial" => "221" "paginaFinal" => "244" ] ] "contieneResumen" => array:2 [ "es" => true "en" => true ] "contieneTextoCompleto" => array:1 [ "en" => true ] "contienePdf" => array:1 [ "en" => true ] "resumenGrafico" => array:2 [ "original" => 0 "multimedia" => array:7 [ "identificador" => "fig0045" "etiqueta" => "Figure 9" "tipo" => "MULTIMEDIAFIGURA" "mostrarFloat" => true "mostrarDisplay" => false "figura" => array:1 [ 0 => array:4 [ "imagen" => "gr9.jpeg" "Alto" => 1641 "Ancho" => 1448 "Tamanyo" => 866681 ] ] "descripcion" => array:1 [ "en" => "Intense deformation (a and b) recorded in the base of the Cañada de Tarachi Unit. Note the scale of deformation. Reverse and normal faults plus folds affect the sequence. Z folds recorded in the Cañada de Tarachi Unit (c, e, d) suggest transport to the south-west. In a, the bed where the man stand is the base of the detachment. The S fold on f crops out along the Arivechi-Tarachi road, also suggesting transport to the west." ] ] ] "autores" => array:1 [ 0 => array:2 [ "autoresLista" => "José Luis Rodríguez Castañeda, Jaime Roldán Quintana, Amabel Ortega Rivera" "autores" => array:3 [ 0 => array:2 [ "nombre" => "José Luis Rodríguez" "apellidos" => "Castañeda" ] 1 => array:2 [ "nombre" => "Jaime Roldán" "apellidos" => "Quintana" ] 2 => array:2 [ "nombre" => "Amabel Ortega" "apellidos" => "Rivera" ] ] ] ] ] "idiomaDefecto" => "en" "EPUB" => "https://multimedia.elsevier.es/PublicationsMultimediaV1/item/epub/S0016716915000215?idApp=UINPBA00004N" "url" => "/00167169/0000005400000003/v2_201509120041/S0016716915000215/v2_201509120041/en/main.assets" ] "en" => array:19 [ "idiomaDefecto" => true "titulo" => "Assessment of groundwater flow in volcanic faulted areas. A study case in Queretaro, Mexico" "tieneTextoCompleto" => true "paginas" => array:1 [ 0 => array:2 [ "paginaInicial" => "199" "paginaFinal" => "220" ] ] "autores" => array:1 [ 0 => array:4 [ "autoresLista" => "Gil Humberto Ochoa-González, Dora Carreón-Freyre, Mariano Cerca, Margarita López-Martínez" "autores" => array:4 [ 0 => array:3 [ "nombre" => "Gil Humberto" "apellidos" => "Ochoa-González" "referencia" => array:1 [ 0 => array:2 [ "etiqueta" => "a" "identificador" => "aff0005" ] ] ] 1 => array:4 [ "nombre" => "Dora" "apellidos" => "Carreón-Freyre" "email" => array:1 [ 0 => "freyre@geociencias.unam.mx" ] "referencia" => array:2 [ 0 => array:2 [ "etiqueta" => "b" "identificador" => "aff0010" ] 1 => array:2 [ "etiqueta" => "*" "identificador" => "cor0005" ] ] ] 2 => array:3 [ "nombre" => "Mariano" "apellidos" => "Cerca" "referencia" => array:1 [ 0 => array:2 [ "etiqueta" => "c" "identificador" => "aff0015" ] ] ] 3 => array:3 [ "nombre" => "Margarita" "apellidos" => "López-Martínez" "referencia" => array:1 [ 0 => array:2 [ "etiqueta" => "d" "identificador" => "aff0020" ] ] ] ] "afiliaciones" => array:4 [ 0 => array:3 [ "entidad" => "Posgrado en Ciencias de la Tierra Centro de Geociencias Universidad Nacional Autónoma de México Campus Juriquilla Blvd. Juriquilla 3001 Col. Juriquilla, Querétaro, C.P. 76230 Querétaro, México" "etiqueta" => "a" "identificador" => "aff0005" ] 1 => array:3 [ "entidad" => "Laboratorio de Mecánica de Geosistemas Centro de Geociencias Universidad Nacional Autónoma de México Campus Juriquilla Blvd. Juriquilla 3001 Col. Juriquilla, Querétaro, C.P. 76230 Querétaro, México" "etiqueta" => "b" "identificador" => "aff0010" ] 2 => array:3 [ "entidad" => "Laboratorio de Mecánica de Geosistemas Centro de Geociencias Universidad Nacional Autónoma de México Campus Juriquilla Campus Juriquilla Blvd. Juriquilla 3001 Col. Juriquilla, Querétaro, C.P. 76230 Querétaro, México" "etiqueta" => "c" "identificador" => "aff0015" ] 3 => array:3 [ "entidad" => "Departamento de Geología Centro de Investigación Científica y Educación Superior de Ensenada Carretera Ensenada-Tijuana No. 3918 Zona Playitas C.P. 22860 Ensenada, B.C., México" "etiqueta" => "d" "identificador" => "aff0020" ] ] "correspondencia" => array:1 [ 0 => array:3 [ "identificador" => "cor0005" "etiqueta" => "⁎" "correspondencia" => "Corresponding author." ] ] ] ] "resumenGrafico" => array:2 [ "original" => 0 "multimedia" => array:7 [ "identificador" => "fig0045" "etiqueta" => "Figure 9" "tipo" => "MULTIMEDIAFIGURA" "mostrarFloat" => true "mostrarDisplay" => false "figura" => array:1 [ 0 => array:4 [ "imagen" => "gr9.jpeg" "Alto" => 1160 "Ancho" => 1808 "Tamanyo" => 208188 ] ] "descripcion" => array:1 [ "en" => "Measured drawdown during the pumping test at El Nabo (blue squares) and calculated curve using the implemented model. The derivative of the data obtained by applying the fixed-end-point method (<a class="elsevierStyleCrossRef" href="#bib0040">Bourdet et al., 1989</a>) is presented in a orange line. Derivative drawdown was analyzed by the semi-analytical solution of <a class="elsevierStyleCrossRef" href="#bib0050">Butler and Liu (1991)</a> obtaining 1:2 and 1:1 slopes." ] ] ] "textoCompleto" => "IntroductionUnderstanding groundwater flow systems is crucial for the management of groundwater resources in Mexico because many areas in arid zones depend on groundwater as the main source for urban, agricultural, and industrial development. In the Valley of Queretaro Aquifer (VQA), many wells are drilled near or along major faults and the high variability of piezometric levels in relatively short distances evidence changes on hydraulic properties (<a class="elsevierStyleCrossRef" href="#bib0060">Carreón-Freyre et al., 2005;</a><a class="elsevierStyleCrossRef" href="#fig0005">Figure 1</a>). Furthermore, the spatial distribution of extraction wells and drawdown zones in the VQA during the last 30 years suggests that faults compartmentalize the aquifer system (<a class="elsevierStyleCrossRef" href="#bib0060">Carreón-Freyre et al., 2005</a>). In the VQA, groundwater is extracted from depths greater than two hundred meters; at this depth, geological heterogeneity should be considered for evaluating groundwater flow patterns to allow an accurate estimation of groundwater resources. Thus, the VQA could be considered a faulted aquifer, where faults and fracture zones strongly influence the overall flow dynamics of the system (<a class="elsevierStyleCrossRefs" href="#bib0055">Caine et al., 1996; Allen and Michel, 1999</a>) and the assessment of the hydrogeologic behavior of major structures can be as important as obtaining the hydraulic properties of the stratigraphic sequences. This work examines the role of major structures (e.g., faults and volcanic or plutonic bodies) as hydrogeological independent units, rather than only as hydraulic boundaries, and presents a numerical modeling of the groundwater flow conditions in the faulted VQA. Integration of representative geological data in the model included: (a) vertical variations by assigning a hydraulic conductivity (K) to each geological unit; (b) faults; (c) lateral truncation of the hydro-stratigraphic units and; (d) the presence of intrusive rocks (<a class="elsevierStyleCrossRefs" href="#bib0175">Walton, 1970; Ingebritsen et al., 2006; Zoback, 2007</a>), in particular dikes (<a class="elsevierStyleCrossRef" href="#bib0115">Kulkarni et al., 2000</a>).<elsevierMultimedia ident="fig0005"></elsevierMultimedia>The study area is located in the northwestern part of the VQA system (<a class="elsevierStyleCrossRef" href="#fig0005">Figure 1</a>). The Valley of Queretaro is located in the central part of the Miocene to Recent Trans Mexican Volcanic Belt (TMVB; <a class="elsevierStyleCrossRefs" href="#bib0085">Ferrari, 2000; Alaniz-Alvarez et al., 2002</a>) and is delimited by groups of N-S and NE-SW trending faults active at least since the Miocene (<a class="elsevierStyleCrossRef" href="#bib0015">Alaniz-Alvarez et al., 2001</a>). We aim at having a better understanding on the role of major structures in groundwater flow patterns in the VQA by integrating them in a numerical model.Faulted aquifer systemsPrevious works have documented the variations in hydraulics properties induced by the complex architecture of regional faults (<a class="elsevierStyleCrossRef" href="#bib0125">MaClay and Small, 1983</a>; López and Smith, 1995; López and Smith, 1996; <a class="elsevierStyleCrossRefs" href="#bib0055">Caine et al., 1996; Allen and Michel, 1999; Carreón-Freyre et al., 2005; Anderson, 2006; Mayer et al., 2007; Delinom, 2009; Rafini and Larocque, 2012</a>). Anisotropy of hydraulic conductivity in faulted and stratified geological media determines groundwater flow patterns (<a class="elsevierStyleCrossRef" href="#bib0100">Hsieh and Neuman, 1985</a>). Regional faults represent an important heterogeneity influencing overall aquifer dynamics, not only determining preferential flow directions but also modifying their storage capacity (<a class="elsevierStyleCrossRefs" href="#bib0135">Pacheco, 2002; Burbey, 2008; Bense et al., 2013</a>). Lithological changes and fractures in volcanic rocks, on the other hand, modify hydrogeological properties of geological units. A fault can be a preferential conduit for vertical flow, i.e., a “leaky fault”, when the damaged zone is well developed and is laterally sealed (<a class="elsevierStyleCrossRef" href="#bib0185">Zeidouni, 2012</a>). Fault zones can also behave as leakage areas between aquifers by juxtaposing aquifers that would otherwise be separated by impermeable layers.Major faults can behave as flow barriers or preferential channels depending on the hydraulic properties developed during deformation history (<a class="elsevierStyleCrossRefs" href="#bib0125">MaClay and Small, 1983; Caine et al., 1996; Carreón-Freyre et al., 2005; Mayer et al., 2007; Delinom, 2009</a>, <a class="elsevierStyleCrossRef" href="#bib0140">Rafini and Larocque, 2012</a>), or they can only separate units of different hydraulic conductivity (<a class="elsevierStyleCrossRefs" href="#bib0025">Allen and Michel, 1999; Mayer et al., 2007</a>); for instance, the vertical displacement of a fault can juxtapose permeable and impermeable layers in the direction perpendicular to the fault plane. In early models, faults were commonly considered as border conditions when assessing hydraulic parameters, even in stratified leaky systems (<a class="elsevierStyleCrossRef" href="#bib0145">Rathod and Rushton, 1991</a>).The style of deformation and rupture of a fault zone strongly influences variations of its hydraulic properties. <a class="elsevierStyleCrossRef" href="#bib0055">Caine et al. (1996)</a> proposed four conceptual schemes for fault-related fluid-flow according to the distribution of conduits and barriers (localized and distributed) within fault plain. <a class="elsevierStyleCrossRef" href="#bib0030">Anderson (2006)</a> identified three essential hydraulic elements of a fault plane: (1) a low hydraulic conductivity core characterized by fault breccias, cataclastic rocks, and gouge; (2) an adjacent damaged zone characterized by an increase in the hydraulic conductivity caused by brecciation and fracturing and; (3) fresh unaltered rock with a characteristic hydraulic conductivity that can change in both sides of the fault because of the fault displacement.A decrease in the hydraulic conductivity within a fault plane (flow barrier) can be associated with (<a class="elsevierStyleCrossRefs" href="#bib0130">Mayer et al., 2007; MaClay and Small, 1983</a>): (1) cataclasis or grain size reduction; (2) juxtaposition of permeable and impermeable layers, displacements along the fault plane can decrease the possibility of flow through the fault and/or deviating flow along the fault; (3) rotation of plane and elongate detritus that became parallel to fault plane reducing the hydraulic conductivity in the orientation perpendicular to the plane; (4) presence of sedimentary materials with low hydraulic conductivity (e. g., clays) in the fault zone; and (5) chemical precipitation of mineral in the fault zone reducing permeability.Furthermore, major faults often present associate or secondary planes and fractures that may cause contrasting hydraulic behaviors. In anisotropic media, each geological unit is represented by a hydraulic conductivity tensor considering three dimensional variations of hydraulic conductivity. In fractured rocks the tensor can align following the preferential orientation of fractures (<a class="elsevierStyleCrossRef" href="#bib0110">Jourde et al., 2007</a>).Geology of the Study Area (Northwest of the Valley of Queretaro Aquifer, VQA)The regional geology of the VQA has been extensively studied in the last decades (<a class="elsevierStyleCrossRefs" href="#bib0165">Trejo-Moedano, 1989; Valdez-Moreno et al., 1998; Aguirre-Díaz and Lopez-Martinez, 2001</a>; Alaníz-Alvarez et al., 2001, 2002; <a class="elsevierStyleCrossRefs" href="#bib0060">Carreón-Freyre et al., 2005; Aguirre-Díaz et al., 2005; Xu et al., 2011</a>). The Valley of Queretaro is a basin delimited by groups of N-S and NE-SW trending faults and volcanoes, located in the central part of the Trans Mexican Volcanic Belt (TMVB; <a class="elsevierStyleCrossRefs" href="#bib0085">Ferrari, 2000; Alaniz-Alvarez et al., 2002</a>; <a class="elsevierStyleCrossRef" href="#fig0005">Figure 1</a>). The N-S trending west-dipping Central fault is the eastern limit of the regional Queretaro graben (<a class="elsevierStyleCrossRef" href="#fig0005">Figure 1</a>). The other major faulting family (trending NE-SW) produced a nearly orthogonal pattern with the N-S system that has resulted in a mosaic formed by horsts, grabens, and half-grabens (Alaníz-Alvarez et al., 2001, 2002; <a class="elsevierStyleCrossRefs" href="#bib0060">Carreón-Freyre et al., 2005; Aguirre-Díaz et al., 2005; Xu et al., 2011</a>) with varying vertical displacements, in some areas up to 400 m. Both fault systems have been active at least since the Miocene (e. g., <a class="elsevierStyleCrossRefs" href="#bib0015">Alaniz-Alvarez et al., 2001; Zuñiga et al., 2003; Aguirre-Díaz et al., 2005</a>) and both controlled the deposition of sediments and volcanic rocks in the VQA.This study is focused on the northwestern part of the VQA (dashed square in <a class="elsevierStyleCrossRefs" href="#fig0005">Figures 1 and 2</a>), a zone of groundwater recharge because of faulting (<a class="elsevierStyleCrossRef" href="#bib0070">Conagua, 2003</a>). The available information was compiled and integrated in a geological map with the hydrogeologic information of the VQA and surroundings. The local stratigraphy was established from the reinterpretation of well logs, 40Ar*/39Ar dating of magmatic rocks, and geologic mapping. The structural and stratigraphic setting of this part of the basin can be defined with sufficient detail and indicates that the geological evolution is more complex than hitherto recognized (<a class="elsevierStyleCrossRef" href="#fig0010">Figure 2</a>a). In this area, the Cenozoic volcanic and colluvial sequences cover partially the Cretaceous units that were uplifted by a shortening episode and thus are well exposed in the hanging wall. A schematic W-E geological cross section AA’ was constructed using the lithological logs of water extraction wells available in the zone (Mompani II, El Nabo, UNAM, and El Salitre) (<a class="elsevierStyleCrossRef" href="#fig0010">Figure 2</a>).<elsevierMultimedia ident="fig0010"></elsevierMultimedia>The oldest sequence (Cretaceous) consists of volcano-sedimentary rocks, documented here for the first time in the Valley of Queretaro (Ki MsV, see below), a thick Albian limestone and shale sequence (Ki Cz, of 800 m thick), and a late Cretaceous shale and sandstone interbedded with argillaceous limestone (Ki Cz-Ar). The Cretaceous sequence is pervasively folded by shortening deformation.The Cenozoic sequence consists of a continental sedimentary unit with predominance of conglomerates and sandstone (Te Cg-Ar), overlain by a volcanic sequence with pyroclastic and sedimentary interbedded deposits. The Cenozoic sequence is fractured and affected by extensional faults. The Te Cg-Ar layer of probable Eocene age is present at the Mompani II well (at a depth of 200 m), and at the bottom of UNAM-Juriquilla and El Salitre wells (<a class="elsevierStyleCrossRef" href="#fig0015">Figure 3</a>a). The volcanic rock sequence includes lavas such as the Salitre Andesite (ToA) and silicic domes (Tom Da; Juriquilla dacite; <a class="elsevierStyleCrossRef" href="#bib0020">Alaniz-Alvarez et al., 2002</a>) of the Oligocene and early Miocene. An intrusive body of Oligocene age (ca. 30 Ma) is exposed in the zone of Juriquilla within the study area (<a class="elsevierStyleCrossRef" href="#fig0010">Figure 2</a>a). The pyroclastic and lacustrine sequences in the Queretaro area consist of tuffs and other pyroclastic deposits of 80 m in average of thickness (Tom PyLac); and this unit includes the Ezequiel Montes Pumice, an important stratigraphic marker dated between 7.5 and 5.6 Ma (<a class="elsevierStyleCrossRef" href="#bib0005">Aguirre-Díaz and López-Martínez, 2001</a>). The coarse facies of this sequence (Tom PyLac) was the main hydrogeological unit exploited in the regional VQA, but is scarcely exposed in the study area, it is found in the upper part of the Mompani II log (<a class="elsevierStyleCrossRef" href="#fig0015">Figure 3</a>). Fissure basalts and lava flows (Tmt AB), ranging in age from 12 to 5 Ma (<a class="elsevierStyleCrossRef" href="#bib0005">Aguirre-Díaz and López-Martínez, 2001</a>) erupted from nearby volcanoes. The largest volcanic edifice, located 13km to the NW from the studied area, is the La Joya Volcano ranging in age from 12 to 10 Ma (<a class="elsevierStyleCrossRef" href="#bib0170">Valdez-Moreno et al., 1998</a>).<elsevierMultimedia ident="fig0015"></elsevierMultimedia>The truncation and throw of the N-S trending extensional faults was inferred by the elevation differences among stratigraphic levels from the correlation of lithological logs presented in <a class="elsevierStyleCrossRef" href="#fig0015">Figure 3</a>. The vertical displacements are smaller than 200 m. From <a class="elsevierStyleCrossRef" href="#fig0010">Figure 2</a>b note that two additional buried faults nearby, and to west of UNAM well, were inferred in this work. Structurally, the study area is a block delimited by the major N-S trending Balvanera and El Nabo fault zones dipping to the northeast (<a class="elsevierStyleCrossRef" href="#fig0020">Figure 4</a>a). The South fault and other NE-SW trending fault segments delineate the limits of the block toward the south. Both faults systems have nearly vertical fault planes and their kinematics indicate dip-slip towards the center of the Queretaro Valley since Miocene time. These major faults include several parallel fault planes affecting areas with more than 50 m in width (<a class="elsevierStyleCrossRef" href="#fig0020">Figure 4</a>). Faults are partially exposed and their trace length was inferred by topography, field observations, and lithological correlations of well logs. The relative age of faulting in the region was analyzed in detail by <a class="elsevierStyleCrossRef" href="#bib0015">Alaniz-Alvarez et al (2001)</a>; the older system corresponds to the South Fault that favors the emplacement of the Oligocene intrusive body. The Balvanera and El Nabo fault systems are younger and favor the Miocene volcanism through dykes emplaced along the faults planes (<a class="elsevierStyleCrossRef" href="#fig0020">Figure 4</a>). The El Nabo Fault is of particular relevance to this work, it is a nearly vertical fault partially exposed to the north of the study area where plateau lavas fill the hanging wall of the fault (solid line with thick marks in <a class="elsevierStyleCrossRef" href="#fig0010">Figure 2</a>a). Toward its southern part the fault is not clearly observed on surface but it can be inferred from the topographic relief and the interpretation of the stratigraphic records (dotted line in <a class="elsevierStyleCrossRef" href="#fig0010">Figure 2</a>a). A dyke emplaced along a fault plane parallel to El Nabo fault was observed at the surface northeastward of the El Nabo well (<a class="elsevierStyleCrossRef" href="#fig0020">Figure 4</a>b). All this information was integrated to constrain the numerical model.<elsevierMultimedia ident="fig0020"></elsevierMultimedia>MethodologyStratigraphic correlationThe conceptual model of the study area is based on geological units with contrasting hydraulic properties, granular and fractured rock layers, truncated by major faults and volcanic edifices/bodies that may interrupt or favor lateral flows. <a class="elsevierStyleCrossRef" href="#bib0060">Carreón-Freyre et al. (2005)</a> correlated the hydrogeological units with lithostratigraphic units to describe the VQA and proposed a conceptual model of multilayer aquifer with groundwater flow varying from local to regional. The stratigraphic correlation of lithological logs for the Cenozoic was based on the reconnaissance of an andesitic-basaltic unit of late Miocene age, named Tmt AB, as an index layer. For instance, as shown in the record of the Ejido Jurica well (<a class="elsevierStyleCrossRef" href="#fig0015">Figure 3</a>a) the Tmt AB layer commonly is overlaid by a Miocene-Pliocene sedimentary unit, named Tmp Cg-Ar-Lm, and underlied by an Oligocene-Miocene pyroclastic and lacustrine unit, named Tom PyLac.The Cretaceous units are more relevant to groundwater flow in this area and include sandstone and argillaceous sequences, limestone rocks, and a volcano-sedimentary sequence. These units were found during drilling of El Nabo well, which is the deepest in the study area and is located within a block where Tertiary sequences are partially eroded (<a class="elsevierStyleCrossRef" href="#fig0015">Figure 3</a>b). The contrasting stratigraphy and the vertical difference in elevation of the Ki Cz unit from 1680 to 1100 masl among the Mompani II and El Nabo records separated by approximately one kilometer suggested the presence of the El Nabo Fault at this location.Drilling of the first 600 m depth of El Nabo well initiated in 2000 at this time the production of water was approximately 6 l/s. Extended drilling in El Nabo well from 600 to 1000 m depth was performed in 2007, the production of the well increased dramatically to 47 l/s after drilling at 850 m depth, that coincided with a lithological change from limestone to a volcano-sedimentary sequence (volcanic rocks and limestone) bearing biotite minerals). Further evidence of the fault and its relevance to groundwater flow was obtained by the geochronological data of rock samples from the El Nabo well presented below.Biotite age determinationIn order to define the age of the basal volcanic unit, rock samples collected each 3 m during drilling from 800 to 1050 m depth were classified by the use of a 20x-40x stereoscopic magnifier, and some samples of interest were selected for thin section description. Biotite concentrates were selected for dating at depths between 850 and 950 m. These biotite concentrates were analyzed by the 40Ar*/39Ar technique at the Geochronology Lab in CICESE using a VG5400 mass spectrometer and a laser extraction system. The analytical details of the method can be consulted elsewhere (e.g., <a class="elsevierStyleCrossRef" href="#bib0065">Cerca et al., 2007</a>). Mineral concentrates were previously irradiated at the nuclear reactor of the University of McMaster in Hamilton, Ontario, Canada, using as monitor of irradiation two known-age samples of sanidine, FCT-2 of 27.84 ± 0.04 Ma and TCR-2 of 27.87 ± 0.04 Ma. All the experiments were made by step-heating by increasing the laser power from 0.2 to 8 watts and measuring the isotopic composition of argon released during each step. The reported age is calculated for each individual fraction and the integrated results adding individual steps in accordance to <a class="elsevierStyleCrossRef" href="#bib0090">Hall (1981)</a>, and considering the corrections by discrimination, radioactive decay (<a class="elsevierStyleCrossRef" href="#bib0155">Steiger & Jäger, 1977</a>), the line background, and interference reactions from isotopes derived from Ca, Cl, and K. Results were grouped according to depth and size of the samples and are presented in a graph named age spectra. Correlation diagrams in which the relation (40Ar*/39ArK) is given by the inverse of the intercept of the line defined by the data and is equivalent to the age in Ma are presented in a supplementary Figure (S1).Integration of conceptual and numerical modelsThe conceptual model resulting from the interpretation of lithological logs, geological mapping, analysis of samples, and isotopic data is shown in <a class="elsevierStyleCrossRef" href="#fig0025">Figure 5</a> (a, b). The wells integrated in the numerical model are El Nabo and Mompani II, 1100 and 250 m depth, respectively (<a class="elsevierStyleCrossRef" href="#fig0025">Figure 5</a>a). For the implementation of the numerical model, the volcanic sequences, sedimentary materials, and geological faults were considered as independent hydrogeological units; each one with a hydraulic conductivity tensor.<elsevierMultimedia ident="fig0025"></elsevierMultimedia>Groundwater simulation was accomplished with the Visual-Modflow software (Schlumberger Water Services, Inc) and the Modflow 2000 code (<a class="elsevierStyleCrossRef" href="#bib0095">Harbaugh et al., 2000</a>). In order to represent the hydraulic behavior of the mapped geological structures the numerical model was built with 100,000 cells, distributed in a horizontal grid of 100 x 100 cells in the N-S and E-W direction, each one represents a square of 2500 m2. The model box is a square of 25km2 with a top surface with an average elevation of 1945 masl reproducing the actual topography and a flat surface at the bottom of the model (800 masl). The model box was divided in 10 horizontal layers of cells with variable thickness from the surface, according to the topography, to reach a bottom flat layer. Cells with similar hydraulic behavior were grouped to represent the geological units included in the model considering stratigraphic sequences and major structures. According to their physical characteristics, spatial distribution, and thickness, four geological units were defined in the conceptual model and represented in the numerical model shown in <a class="elsevierStyleCrossRef" href="#fig0025">Figure 5</a>a:<ul class="elsevierStyleList" id="lis0005"><li class="elsevierStyleListItem" id="lsti0005">(1)TomDa-Rhyolite, this is the upper layer of the model, and groups the Cenozoic volcanic rocks including the volcanic and pyroclastic fractured sequences located above the piezometric level in the study area. These units, observed in the lithological records of the Ejido Jurica, UNAM-Juriquilla, and Ejido El Nabo wells, were not considered in the conceptual model because either they are located above the regional groundwater level (vadose zone) or are absent nearby the El Nabo well area (<a class="elsevierStyleCrossRef" href="#fig0015">Figure 3</a>). These units include: Late Miocene Andesites and Basalts (Tmt AB), Early Miocene Conglomerates and Sandstones Sandstones (Tmp Ar-Cg), and Quaternary Alluvium (Qal; <a class="elsevierStyleCrossRef" href="#fig0010">Figure 2</a>);</li><li class="elsevierStyleListItem" id="lsti0010">(2)Te Cg Ar Conglomerate and sandstone, a sedimentary unit of conglomerate and sandstone found at the west of the El Nabo fault, was represented by a dark layer shown in the model of <a class="elsevierStyleCrossRef" href="#fig0025">Figure 5</a>a;</li><li class="elsevierStyleListItem" id="lsti0015">(3)Ki Cz Ar Limestone and shale, this unit includes Ki Cz and also the Ki Cz-Ar of sandstone and argillaceous limestone for modeling purposes, and;</li><li class="elsevierStyleListItem" id="lsti0020">(4)Ki MsV Metasedimentary and volcanic (<a class="elsevierStyleCrossRef" href="#fig0025">Figure 5</a>), located at the bottom of the model, corresponds in nature to the Cretaceous volcano-sedimentary unit (Ki MsV).</li></ul>Faults are represented by groups of 50 m width vertically aligned cubic cells with specific hydraulic conductivity (gray cells in <a class="elsevierStyleCrossRef" href="#fig0025">Figure 5</a>b). The massive intrusive/subvolcanic body of 30 Ma is represented by a group of cells with a dome shape in the located central part of the model (black cells in <a class="elsevierStyleCrossRef" href="#fig0025">Figure 5</a>b). Major geological structures were considered as four individual units in the model: (5) the El Nabo Fault, (6) the Balvanera Fault, (7) the South Fault and, (8) the intrusive body placed in the lower part of the model (<a class="elsevierStyleCrossRef" href="#fig0010">Figure 2</a> and <a class="elsevierStyleCrossRef" href="#fig0025">Figure 5</a>b).Faulting and the intrusive body are thought to greatly influence groundwater by channelizing or arresting flow, or by interrupting the lateral continuity of the hydro stratigraphic units, such as the conglomerate and sandstone unit (Te Cg-Ar). Truncation by faulting was suitable represented because each geological unit was composed by independent cells.Once designed the geometry of the model, the initial and border conditions were established considering regional variations of the piezometric surface interpolated from 2004 and 2006 measurements (<a class="elsevierStyleCrossRef" href="#fig0030">Figure 6</a>). Previous work by Carreon-Freyre et al. (2005) indicated a regional gradient in the groundwater flow with a nearly north to south direction in the VQA. Regional piezometric were interpolated by krigging of the data from 40 water extraction wells located in the VQA, systematically monitored by the local agency of water (Comisión Estatal de Aguas de Querétaro, CEAQ) since 1996 (<a class="elsevierStyleCrossRef" href="#bib0060">Carreón-Freyre et al., 2005</a>). The interpolations show a persistent regional gradient with an approximate North-South trending. The wells located in the study area are El Nabo and Mompani II (<a class="elsevierStyleCrossRef" href="#fig0005">Figure 1</a> and <a class="elsevierStyleCrossRef" href="#fig0025">Figure 5</a>a). A steady state numerical simulation was performed to establish the initial and boundary conditions (<a class="elsevierStyleCrossRef" href="#fig0035">Figure 7</a>b). Boundary effects were minimized by leaving a minimum radius of 2000 m from each well to the model borders.<elsevierMultimedia ident="fig0030"></elsevierMultimedia><elsevierMultimedia ident="fig0035"></elsevierMultimedia>Simulation of hydrodynamic conditions using the numerical modelA numerical simulation was performed to assess the hydrodynamic behavior of the system by the use of partial results of a pumping test. In the computational experiment the hydraulic properties were adjusted to correlate calculated and measured drawdown in the pumping and in one observation well.A 260,000 s (72 hour) pumping test performed in the El Nabo well was simulated using the implemented numerical model. The location of the pumping test area is shown in <a class="elsevierStyleCrossRef" href="#fig0010">Figure 2</a>; the data record (water level and temperature) was made by a CTD Diver (Schlumberger) datalogger. This well is the single one in the region with about 1000 m in depth and therefore it was not possible to locate any observation well. The Mompani II well is located 1000 m to the SW, with a depth of 250 m. Because of the presence of the El Nabo Fault between the two wells no change in the hydraulic head during the pumping test was recorded at the Mompani II and thus, we assumed that groundwater flow in both wells were independent (possibly local and regional systems).To accomplish the simulation, the hydraulic conductivities of the model and measured drawdown were closely fitted. We used a segment of the pumping test to simulate transient conditions generated by the pumping test. Continuous measurements of drawdown and temperature every 15 s are presented in <a class="elsevierStyleCrossRef" href="#fig0040">Figure 8</a> (blue and green line respectively); blue squares indicate samples from the continuous measurements used for simulation. Because of the well depth only a part of the pumping test could be considered for adjustment. The first 10,800 s data (3hours) were not considered because the pumping rate could not be maintained stable during this period and/or a well storage effect should be considered (<a class="elsevierStyleCrossRef" href="#bib0025">Allen and Michel, 1999</a>), the same issues were considered for the interpretation of the recovery data.<elsevierMultimedia ident="fig0040"></elsevierMultimedia>Simulated and observed drawdown in the segment between 10,800 and 250,000 s (3 and 69hours) of the El Nabo test were fitted by adjusting the hydrogeological properties of the main units (e.g., <a class="elsevierStyleCrossRef" href="#bib0025">Allen and Michel, 1999</a>). The estimation of these properties was achieved using the PEST (Parameter Estimation Software) algorithm (<a class="elsevierStyleCrossRefs" href="#bib0095">Harbaugh et al., 2000; Doherty, 2005</a>) trough Visual-Modflow. The finite differences Modflow algorithm solved groundwater flow by adjusting hydraulic conductivities in three perpendicular directions, corresponding with the orientation of the grid to North-South, East-West, and vertically for this study case. The initial value of hydraulic conductivity (K) was established at 1e-04 m/s, corresponding to reported values for massive and fractured rocks (<a class="elsevierStyleCrossRefs" href="#bib0175">Walton, 1970; Ingebritsen et al., 2006; Zoback, 2007</a>).The best fit of the observed and calculated drawdowns was achieved by using the hydraulic conductivities (K) values presented in <a class="elsevierStyleCrossRef" href="#tbl0005">Table 1</a>, along two directions: horizontal E-W and vertically. The North-South values of K remained fixed at the initial value. The sensitivity values for the obtained hydraulic conductivities (<a class="elsevierStyleCrossRef" href="#tbl0005">Table 1</a>) of the model vary with the distance between the geological units/faults and the pumping test zone. High sensitivity means a high influence of the geological unit in drawdown. The “limestone and shale” and “El Nabo fault” units show the highest sensitivity values (<a class="elsevierStyleCrossRef" href="#tbl0005">Table 1</a>b). The high values in the El Nabo Fault are due to the fact that the pumping well is located along this fault. The lower obtained sensitivities correspond to the Balvanera and South faults, meaning that the conductivity values of these units did not affect the calculated drawdown.<elsevierMultimedia ident="tbl0005"></elsevierMultimedia>The results of the pumping test were analyzed by the derivative of the data obtained by applying the fixed-end-point method (orange line) proposed by <a class="elsevierStyleCrossRef" href="#bib0040">Bourdet et al. (1989)</a> (orange line in <a class="elsevierStyleCrossRef" href="#fig0045">Figure 9</a>). The use of the derivative allows analyzing in detail pressure changes in drawdown plots (<a class="elsevierStyleCrossRefs" href="#bib0040">Bourdet et al., 1989; Spane and Wurstner, 1993; Tiab, 2005</a>). The derivative was smoothed by the use of the third point before and after the point of interest to reduce the sample rate noise. The intersection between dyke/fault and the El Nabo well was analyzed by the semi-analytical solution of <a class="elsevierStyleCrossRef" href="#bib0050">Butler and Liu (1991)</a>.<elsevierMultimedia ident="fig0045"></elsevierMultimedia>Discussion of resultsThe implemented model allowed representing the faulted system of the study area by the hydraulic conductivities assigned to specified geological units, including faults. In order to constrain the overall flow conditions regional measurements of hydraulic heads and geological observations were integrated.Geological constraints for the imple-mentation of the modelThe Cretaceous sequence drilled in El Nabo well, composed by a limestone and shale unit (Ki Cz-Ar) that overlies a volcano-sedimentary sequence (Ki Mvs), was of the most importance to establish the geological constraints of the model. The volcanic character of the basal volcano-sedimentary sequence was interpreted from the drilling samples and was confirmed by the isotopic ages. The overlying limestone and shale unit is more than 800 m thick and is inter-bedded both with sandstone or clay. In nearby superficial outcrops this unit is deformed by folding and contains calcite veins and fractures. Rock samples obtained from depths below 850 m at the El Nabo well of the volcanic sequence contained a large quantity of biotite. Differences in the size and aspect of biotite defined two size populations: (1) between 425 and 710μm and (2) between 250 and 425μm. The age spectra obtained for both sizes confirms that they correspond also to two different age populations.<a class="elsevierStyleCrossRef" href="#fig0050">Figure 10</a>a shows the age spectra for the biotite population 1 indicating a Miocene age. The weighted average for argon fractions was made by considering the consecutive fractions with greater radiogenic argon content and more than 10% of the total 39Ar released. A representative age of t = 9.5 ± 0.1 Ma was selected from the weighted average of 13 fractions from different experiments in the age spectra diagram. This biotite population is characterized by a ratio 37ArCa/39ArK ≤ 0.5. The correlation diagrams obtained for these samples (<a class="elsevierStyleCrossRef" href="#fig0065">Figure S1</a>) show two isochron lines also converging at an age around 9.5 Ma. The Andesitic lava flows (<a class="elsevierStyleCrossRef" href="#bib0170">Valdez-Moreno et al., 1998</a>) cover most of the study area and have the same Miocene age obtained for biotite 1 population, the coincidence in age of the lavas and the biotite sampled at depth strongly support the interpretation that the lavas were feed by dikes emplaced along faults.<elsevierMultimedia ident="fig0050"></elsevierMultimedia><elsevierMultimedia ident="fig0065"></elsevierMultimedia>The biotite population 2 (<a class="elsevierStyleCrossRef" href="#fig0050">Figure 10</a>b) yielded older ages of ∼100 Ma with age spectra diagram indicating a severe argon loss. A reliable age cannot be obtained but the sequence is older than the minimum age obtained. This behavior was reproduced in the 855 m depth Biotite in duplicate experiments. A higher ratio of 37ArCa/39ArK > 0.5 compared with biotite population 1 clearly indicates a different composition. Because of the different isotopic behavior of these samples the biotite populations were plotted separately in the correlation diagram (<a class="elsevierStyleCrossRef" href="#fig0065">Figure S1</a>), in which the equivalent age is shown in the abscissa axis. The correlation using the fractions from 3 to 6 of the 1st experiment and the fractions 3 to 8 of the duplicate experiment in the biotite 2 population also indicates an age of ∼100 Ma, for the Ki Ms V rock sequence.The documented difference in age of the biotite populations strongly supports our interpretation of the presence of a Miocene magmatic dyke intruding the Cretaceous sequence as shown in the schema of <a class="elsevierStyleCrossRef" href="#fig0050">Figure 10</a>c. The presence of this fault has important implications in groundwater flow by the increasing of the hydraulic conductivity within the fault plane that was observed during drilling.It has been documented that dikes open weakened fractures to conduct magma toward the surface and during their emplacement they can cause further brecciation of the fractured zone and the increase of hydraulic conductivity (<a class="elsevierStyleCrossRefs" href="#bib0190">Zoback, 2007; Anderson, 2006</a>). Beside the presence of the dike, further evidences supporting the circulation of hydrothermal fluids along the fault include the presence of hydrothermally altered calcite veins in exposures and in drilling samples of the El Nabo well (<a class="elsevierStyleCrossRef" href="#fig0020">Figure 4</a>b). Present day hydrothermal activity along the fault was confirmed by the increase from 25 to 63°C in the temperature records during pumping test (green line, <a class="elsevierStyleCrossRef" href="#fig0040">Figure 8</a>).In order to consider the role of major faults in groundwater flow we integrated the El Nabo, Balvanera and South faults in the model. El Nabo and Balvanera are sub-parallel normal faults with NW-SE orientation (panoramic view, <a class="elsevierStyleCrossRef" href="#fig0020">Figure 4</a>a). The Balvanera Fault has a throw of more than 200 m after the Miocene lavas emplacement. The El Nabo Fault is located in a block delimited by the regional graben faults as shown in the cross-section of <a class="elsevierStyleCrossRef" href="#fig0010">Figure 2</a>b. Several faults and fractures with similar orientation were observed around the zone of the El Nabo well. The Cretaceous sedimentary sequences (Ki Cz-Ar) have a low hydraulic conductivity as shown by the low flow rate measured in the El Nabo well above 800 m depth (6 l/s). Nevertheless, flow rate increased to 47 l/s when drilling intersected the El Nabo Fault and dike at 850 m depth. We infer that the dike caused brecciation, fracturing, and hydrothermal alteration of the host rock and consequently an increase in the hydraulic conductivity. Evidence of other dikes emplaced on the El Nabo Fault area affecting other stratigraphic levels was observed in exposed parallel fault planes. For modeling purposes, the Balvanera Fault features are assumed to be similar to the El Nabo Fault and was calibrated as an independent unit.The exposed plane of the older age South Fault (<a class="elsevierStyleCrossRef" href="#fig0020">Figure 4</a>b) contains open fractures filled with sediments and calcite that precipitated by the circulation of hydrothermal fluids. The South Fault favored the location of the Oligocene sub-volcanic body (30 Ma, TomDa) of silicic composition (Alanis-Alvarez et al., 2002) observed as a dome structure at the surface (<a class="elsevierStyleCrossRef" href="#fig0010">Figure 2</a>). Given its large spatial dimensions this structure was considered as an individual unit in the numerical model (<a class="elsevierStyleCrossRef" href="#fig0025">Figure 5</a>b).The geological evidences suggest that groundwater flow through faults is significant in the study area: thus, the model integrated an increase of conductivity along the El Nabo fault.Assessment of groundwater flow in faulted and high heterogeneous areasThe numerical model was constrained by the results of the pumping test in El Nabo well. The record of the test was sub-sampled and the obtained drawdown plot is presented in <a class="elsevierStyleCrossRef" href="#fig0045">Figure 9</a> (blue squares) along its derivative (<a class="elsevierStyleCrossRef" href="#bib0040">Bourdet et al., 1989</a>). Considering an infinite linear strip embedded in a host material (matrix) with differing hydraulic properties <a class="elsevierStyleCrossRef" href="#bib0050">Butler and Liu (1991)</a> presented derivative curves with slopes varying from 1:2 and 1:4 for a case in which hydraulic conductivity is much larger than the host material (matrix). In their analysis these slopes corresponded to an initial linear flow within the strip that changed to bilinear flow when the flow from host material contributes to discharge. The slopes obtained from the drawdown derivatives of the El Nabo well pumping test are shown in <a class="elsevierStyleCrossRef" href="#fig0045">Figure 9</a>; a slope 1:1 (green dashed line) for the 20,000-70,000 s (5-19 hour) interval and a slope 1:2 (black dashed line) for the last part of the data. The first slope (1:1) corresponds to the primary flow drawdown that may be explained by two conditions: (a) groundwater flow contributed by the limestone unit presenting high K values in the E-W direction and (b) is enclosed within the fault walls (<a class="elsevierStyleCrossRef" href="#tbl0005">Table 1</a>) or, (2) non lineal flow inside the fault walls and constrained by the regional gradient, which magnitude is large enough to oppose flow toward the well.The 1:2 slope may be related either to discharge of a lineal regional flow within the El Nabo Dyke/Fault (i.e. strip unit, <a class="elsevierStyleCrossRef" href="#bib0050">Butler and Liu, 1991</a>) or a new source of regional flow. Support evidence of the contribution of a hydrothermal deep groundwater flow system through the fault was provided by the increase in temperature while pumping (25 to 63°C, green line in <a class="elsevierStyleCrossRef" href="#fig0040">Figure 8</a>). The initial increase of temperature in a time interval of 83minutes suggests mixing of local fresh groundwater and hydrothermal flow, with the presence of a preferential conduct through the magmatic dyke/fault.The calibrated hydraulic conductivities (obtained using PEST software) were used to simulate the pumping test in the numerical model (red line; <a class="elsevierStyleCrossRef" href="#fig0045">Figure 9</a>). The wellbore storage may have caused the rapid initial drawdown of the El Nabo pumping test (the well is more than 1000 m depth and 6 inches in diameter). Time required for pumping stabilization (10,000 s, or approximate 3hours, <a class="elsevierStyleCrossRef" href="#fig0040">Figure 8</a>) was not considered in the model, drawdown initiated from a relative zero value. A pseudo-steady flow regime observed in the 10,800-35,000 s (3-10 hour) intervals (<a class="elsevierStyleCrossRefs" href="#fig0040">Figures 8 and 9</a>) was not completely adjusted by the calculated drawdown. From 20 hour (75,000 s) the numerical solution converges with the data and the model results are better fitted (<a class="elsevierStyleCrossRef" href="#fig0045">Figure 9</a>). In spite of the initial differences, we consider that model curve captures the overall behavior of the pumping test suggesting a similarity of hydraulic behavior between model and test.The model allows discussing the contribution of each geological element in the overall hydraulic behavior of the system. Faults were considered as vertical planar structures with varying hydraulic conductivities according to the E-W or vertical groundwater flow direction. The block model of <a class="elsevierStyleCrossRef" href="#fig0055">Figure 11</a> summarize the main groundwater flow results using the K values presented in <a class="elsevierStyleCrossRef" href="#tbl0005">Table 1</a>, and shows the piezometric levels reached at the end of the pumping test simulation. Colors indicate piezometric levels, which range from 1925 in the north area to 1820 masl toward the south. The higher gradient is obtained toward the east of the El Nabo Fault. Arrows, indicating flow directions, show clearly that major structures (El Nabo and South faults) favor ground-water flow patterns along planes and also delimit compartments with differences in hydraulic behavior. Field observations show that the hydraulic interaction of the fault with the El Nabo well starts at 850 m of depth and the presence of hydrothermal groundwater can be interpreted in terms of ongoing deep lateral flow facilitated by the heterogeneous permeability within the fault-zone probably due to the dyke intrusion (López and Smith, 1995; <a class="elsevierStyleCrossRef" href="#bib0120">Lopez and Smith, 1996</a>). Horizontal flow predominates in the model while the flow patterns deviate along the geometric boundaries of the intrusive body, interrupting the connection between El Nabo and South faults, and following the regional N-S gradient.<elsevierMultimedia ident="fig0055"></elsevierMultimedia><a class="elsevierStyleCrossRef" href="#fig0060">Figure 12</a> shows one horizontal and two vertical cross sections (in the E-W and N-S directions) of the model (<a class="elsevierStyleCrossRef" href="#fig0055">Figure 11</a>) that allow to analyze the influence of major structures (El Nabo and South faults and the intrusive body) in groundwater flow patterns. In the horizontal cross section at 1400 masl of <a class="elsevierStyleCrossRef" href="#fig0060">Figure 12</a>a the northern segment of the El Nabo Fault act as a barrier by interrupting groundwater flow from north to south within the limestone and shale unit, also shown in the vertical BB’ profile (<a class="elsevierStyleCrossRef" href="#fig0060">Figure 12</a>c). The hydraulic head observed at 1900 m north from the El Nabo fault decreases to 1840 m in the southern block. In the East-West cross section (AA¿ in <a class="elsevierStyleCrossRef" href="#fig0060">Fig. 12</a>b) the color scale indicates a difference in piezometric elevations of 40 m at both sides of the fault. In the model the regional gradient is interrupted by the El Nabo Fault that acts as a lateral flow barrier in the east-west direction, the vertical equipotential lines indicate vertical flow within the fault. The NE-SW orientation of the South Fault seems to favor channelized lateral flow towards the SW along the main fault plane. The South Fault yielded the highest K value in the horizontal E-W direction (<a class="elsevierStyleCrossRef" href="#tbl0005">Table 1</a>). The model is not sensitive to this structure and this high value may be an artifact produced in the model by the equalization of drawdown and regional gradient. However, field data suggest also a higher K for the South Fault than for El Nabo Fault.<elsevierMultimedia ident="fig0060"></elsevierMultimedia>The 1:2 slope of the derivative curve (orange line, <a class="elsevierStyleCrossRef" href="#fig0045">Figure 9</a>) was interpreted as lineal flow within a high permeability channel with semi-permeable walls, interpretation consistent with the presence of El Nabo Fault-Dyke (<a class="elsevierStyleCrossRef" href="#fig0060">Figure 12</a> b, c). The semi-permeable walls condition in the model is represented by the equipotential lines at both sides of the faults. Smooth gradients within the stratigraphic units increase when crossing the fault indicating a high hydraulic gradient but with flow continuity. The limestone and shale unit also show hydraulic continuity with the fault to the west; whereas in the eastern part have a higher elevation suggesting a greater water contribution from the fault (<a class="elsevierStyleCrossRef" href="#fig0060">Figure 12</a>c).Given that only one observation point could be considered, the El Nabo well is the deepest well drilled in the Queretaro area with most of the well depths about 300 m, analysis of representativeness of the obtained hydraulic conductivity was not possible to do. Nevertheless, the information yielded by the pumping test was useful for the calibration of the numerical model allowing simulating groundwater flow in the study area.In summary, this work attempts to highlight the importance of integrating an accurate geological setting to understand groundwater flow dynamics in models of faulted aquifers in particular in cases were piezometric levels decrease constantly.Fault zones in the study area have more than 50 m in thickness and several fault planes were observed where high and low K segments may coexist. Because the style and relative age of faulting determines variations in its hydraulic properties we consider that in faulted and deep aquifers (more than 200 m depth) the assessment of the overall flow dynamics is as, or more, important as obtaining specific hydraulic properties of the hydro-stratigraphic units.Field evidences of present day hydrothermal water circulation within fault planes in the study area and the model results allow proposing that faults collect groundwater flow and explain the flow rate increase to 47 l/s at El Nabo Fault. Moreover, the Cretaceous sedimentary sequences show an overall low hydraulic conductivity. The magmatic dikes caused brecciation, open fractures, and hydrothermal alteration of the Cretaceous host rock and consequently enhancing its hydraulic conductivity (e.g., <a class="elsevierStyleCrossRef" href="#bib0185">Zeidouni, 2012;</a><a class="elsevierStyleCrossRef" href="#tbl0005">Table 1</a>). The South fault system favored the emplacement of the Oligocene intrusive body. Field observations suggest that this fault also favored fluid-flow over time along open fractures partially filled with sediments and secondary fault planes in a zone more than 50 m in width. We have observed precipitation of calcite in fractures (<a class="elsevierStyleCrossRef" href="#fig0020">Figure 4</a>b), and probably infiltration of superficial water. Correspondingly, after the simulation the higher K values obtained in the East-West direction of the South Fault allowed the hydraulic continuity of the system, according with the regional gradient (<a class="elsevierStyleCrossRef" href="#fig0055">Figure 11</a>).ConclusionsImplementation of the geological structure in a numerical groundwater flow model, along with the estimation of hydraulic conductivities by a pumping test, produced relevant information about flow dynamics in the faulted aquifer system at the north of the Valley of Queretaro. This methodology can be applied to study groundwater flow in other areas of the Queretaro Valley and elsewhere in central Mexico, where water is extracted from depths greater than 200 m and high uncertainty exist about the influence of major structures in groundwater flow conditions.The study area is a block delimited by faults where the older geological units in the Queretaro Valley are exposed. Considering faults and truncation of hydrostratigraphic units in the sequence leaded to the reconnaissance of compartments in the aquifer system also reproduced in the model. The identification of groundwater flow patterns through the aquifer when modeling drawdown was consistent with measured piezometric levels. The derivative analysis of the pumping test suggests the presence of two different groundwater flow systems: (1) a flow, with slope 1:1, related to the limestone and shale sequence and, (2) a regional linear flow, with slope 1:2, transporting water through faults.The numerical model along with the drawdown field observations allowed estimating relative values of hydraulic conductivity for eight geological units (<a class="elsevierStyleCrossRef" href="#tbl0005">Table 1</a>). As expected, the model is sensitive to the fault located near the El Nabo well (<a class="elsevierStyleCrossRef" href="#fig0005">Figure 1</a>). Of importance to the study is the presence of volcano-sedimentary sequences below thick limestone units. Isotopic dating confirmed the presence of a Miocene magmatic dyke emplaced along the El Nabo Fault. Both the Balvanera and El Nabo fault systems favored the emplacement of Miocene volcanism through dykes along the faults planes (<a class="elsevierStyleCrossRef" href="#fig0020">Figure 4</a>). The presence of the fault-dyke was corroborated with stratigraphic correlation of a well log, and the results confirm its relevance to modify groundwater flow. The increased high flow rate in the El Nabo well pumping test (47 l/s) is an evidence of the high influence of the dike that allows groundwater to flow through the fault.The obtained values agree with the conceptual model but connectivity between major faults cannot be evaluated only by the analysis of pumping tests. The numerical model allowed exploring the relevance of geological structures in the simulation of groundwater flow (<a class="elsevierStyleCrossRefs" href="#fig0055">Figures 11 and 12</a>). In particular, the model results highlight how faults and intrusive bodies emplaced along faults modify groundwater flow by changing hydraulic properties and consequently flow direction. The model allowed improving the interpretation of the pumping tests. Numerical modeling is useful for a better understanding of complex aquifer systems that cannot be evaluated by analytical methods." "textoCompletoSecciones" => array:1 [ "secciones" => array:14 [ 0 => array:3 [ "identificador" => "xres551697" "titulo" => "Resumen" "secciones" => array:1 [ 0 => array:1 [ "identificador" => "abst0005" ] ] ] 1 => array:2 [ "identificador" => "xpalclavsec569110" "titulo" => "Palabras clave" ] 2 => array:3 [ "identificador" => "xres551696" "titulo" => "Abstract" "secciones" => array:1 [ 0 => array:1 [ "identificador" => "abst0010" ] ] ] 3 => array:2 [ "identificador" => "xpalclavsec569111" "titulo" => "Keywords" ] 4 => array:2 [ "identificador" => "sec0005" "titulo" => "Introduction" ] 5 => array:2 [ "identificador" => "sec0010" "titulo" => "Faulted aquifer systems" ] 6 => array:2 [ "identificador" => "sec0015" "titulo" => "Geology of the Study Area (Northwest of the Valley of Queretaro Aquifer, VQA)" ] 7 => array:3 [ "identificador" => "sec0020" "titulo" => "Methodology" "secciones" => array:4 [ 0 => array:2 [ "identificador" => "sec0025" "titulo" => "Stratigraphic correlation" ] 1 => array:2 [ "identificador" => "sec0030" "titulo" => "Biotite age determination" ] 2 => array:2 [ "identificador" => "sec0035" "titulo" => "Integration of conceptual and numerical models" ] 3 => array:2 [ "identificador" => "sec0040" "titulo" => "Simulation of hydrodynamic conditions using the numerical model" ] ] ] 8 => array:2 [ "identificador" => "sec0045" "titulo" => "Discussion of results" ] 9 => array:3 [ "identificador" => "sec0050" "titulo" => "Geological constraints for the imple-mentation of the model" "secciones" => array:1 [ 0 => array:2 [ "identificador" => "sec0055" "titulo" => "Assessment of groundwater flow in faulted and high heterogeneous areas" ] ] ] 10 => array:2 [ "identificador" => "sec0060" "titulo" => "Conclusions" ] 11 => array:2 [ "identificador" => "xack186041" "titulo" => "Acknowledgements" ] 12 => array:1 [ "titulo" => "Further reading" ] 13 => array:1 [ "titulo" => "References" ] ] ] "pdfFichero" => "main.pdf" "tienePdf" => true "fechaRecibido" => "2013-01-15" "fechaAceptado" => "2014-09-24" "PalabrasClave" => array:2 [ "es" => array:1 [ 0 => array:4 [ "clase" => "keyword" "titulo" => "Palabras clave" "identificador" => "xpalclavsec569110" "palabras" => array:6 [ 0 => "Flujo de aguas subterráneas" 1 => "sistema acuífero fallado" 2 => "modelación numérica" 3 => "estructuras volcánicas" 4 => "diques" 5 => "México" ] ] ] "en" => array:1 [ 0 => array:4 [ "clase" => "keyword" "titulo" => "Keywords" "identificador" => "xpalclavsec569111" "palabras" => array:6 [ 0 => "Groundwater flow" 1 => "faulted aquifers" 2 => "numerical modeling" 3 => "volcanic structures" 4 => "dykes" 5 => "Mexicos" ] ] ] ] "tieneResumen" => true "resumen" => array:2 [ "es" => array:2 [ "titulo" => "Resumen" "resumen" => "En este trabajo se integra la estratigrafía local y fallas como restricciones para un modelo numérico de flujo del agua subterránea en la parte norte del acuífero fallado del Valle de Querétaro. La información geológica de base fue establecida a partir de la reinterpretación de registros litológicos de pozos, datación isotópica 40Ar*/39Ar de rocas y cartografía de campo. En particular, la datación isotópica permitió obtener el tiempo de emplazamiento de un dique magmático a lo largo de la falla de El Nabo. El aumento de la tasa de flujo y la temperatura en el pozo El Nabo sugieren que la falla y el dique transportan un flujo regional. Un modelo numérico fue implementado para evaluar las contribuciones de fallas y estructuras volcánicas en el flujo general de agua subterránea. El modelo realizado con el software Visual-Modflow y código Modflow 2000 fue parcialmente restringido por los resultados de una prueba de bombeo en pozo “El Nabo”, el pozo más profundo (Lat. 20° 42 ‘14 “N, Long. 100° 28’45” O, aproximadamente de 1000 m de profundidad) ubicado en la zona de Querétaro. La interpretación de los resultados del modelo y las observaciones de campo sugieren que las fallas reciben un flujo regional, compartimentalizan el sistema, y son responsables de un importante aumento del gasto, de 6 a 47 l/s, cuando se la perforación del pozo El Nabo bajo de 850 m profundidad. El estudio que aquí se presenta permite la evaluación de fallas normales y intrusivos que determinan el flujo del agua subterránea y podría ser relevante para el conocimiento de la dinámica de flujo valles volcánicos cercanos. El modelo numérico implementado basado en información geológica, permitió una mejor comprensión de los sistemas acuíferos complejos, lo que no se puede evaluar por métodos analíticos. Por otra parte puede ser una aplicación útil para una mejor interpretación de pruebas de bombeo y para añadir datos complementarios de otras pruebas y/o de análisis de hidrogeoquímica." ] "en" => array:2 [ "titulo" => "Abstract" "resumen" => "This work integrates local stratigraphy and faults to constraint a numerical model of groundwater flow in the North area of the Valley of Queretaro aquifer. Basic geological information was established from the reinterpretation of well logs, 40Ar*/39Ar dating of rocks, and field mapping. In particular, dating allowed to obtain the timing of emplacement of a magmatic dyke inferred to be emplaced along the El Nabo fault. The increase of flow rate and temperature in the El Nabo well suggest that the fault and dike transport a regional flow. A numerical model was implemented to evaluate the contributions of faulting and volcanic structures in the overall groundwater flow. The model, accomplished with the Visual-Modflow software and Modflow 2000 code, was partially constrained by the results of a pumping test in the El Nabo well, the deepest well (Lat. 20°42’14”N, Long. 100°28’45”W, ca. 1000 m depth) located in the Queretaro area. The interpretation of the model results and field observations suggest that faults collect the regional flow, compartmentalize the reservoir, and are responsible of the dramatic increase of flow rate from 6 to 47 l/s when drilling the El Nabo well below a depth of 850 m. The study presented here allowed the assessment of normal faults and intrusives determining groundwater flow and might be relevant to the knowledge of flow dynamics in nearby volcanic valleys. The analysis of the pumping test suggest the presence of two different groundwater flow systems: (1) a flow related with the limestone and shale unit and (2) a regional linear flow through faults. The implemented numerical model, based on geological information, allowed a better understanding of complex aquifer systems that cannot be evaluated by analytical methods. Moreover it can be a useful application for a better interpretation of pumping tests and to add complementary data of other pumping tests and/or hydrogeochemical analysis." ] ] "multimedia" => array:14 [ 0 => array:7 [ "identificador" => "fig0005" "etiqueta" => "Figure 1" "tipo" => "MULTIMEDIAFIGURA" "mostrarFloat" => true "mostrarDisplay" => false "figura" => array:1 [ 0 => array:4 [ "imagen" => "gr1.jpeg" "Alto" => 1939 "Ancho" => 1579 "Tamanyo" => 560678 ] ] "descripcion" => array:1 [ "en" => "(a) The Valley of Queretaro is located in the central part of Mexico; (b) Geological map of the Valley of Queretaro. Dashed red square line on the geological map indicates the area considered in modeling, black lines are major faults. Green polygon represents the Cretaceous rocks exposed in the zone. Medium and light gray units represent Oligocene and Miocene volcanic rocks, respectively. (Modified form Carreon-Freyre et al., 2005)." ] ] 1 => array:7 [ "identificador" => "fig0010" "etiqueta" => "Figure 2" "tipo" => "MULTIMEDIAFIGURA" "mostrarFloat" => true "mostrarDisplay" => false "figura" => array:1 [ 0 => array:4 [ "imagen" => "gr2.jpeg" "Alto" => 2554 "Ancho" => 1572 "Tamanyo" => 555033 ] ] "descripcion" => array:1 [ "en" => "(a) Simplified geological map of the study area, the main structures considered in the model are the normal faults El Nabo, Balvanera and South. Logs of water extracting wells were projected to the A-A¿ section shown in (b). The main stratigraphic units are presented in the cross section." ] ] 2 => array:7 [ "identificador" => "fig0015" "etiqueta" => "Figure 3" "tipo" => "MULTIMEDIAFIGURA" "mostrarFloat" => true "mostrarDisplay" => false "figura" => array:1 [ 0 => array:4 [ "imagen" => "gr3.jpeg" "Alto" => 2011 "Ancho" => 1428 "Tamanyo" => 260191 ] ] "descripcion" => array:1 [ "en" => "Stratigraphic correlation of the lithological logs of extraction wells used for the interpretation of the main structures and construction of the conceptual model. (a) Correlation of wells in the NW-SE direction including Mompani II, UNAM-Juriquilla, Ejido Jurica, and (CEA) El Salitre; (b) correlation in the N-S direction Ejido El Nabo, Mompani II, and El Nabo. The initial static water level in the El Nabo well is also shown." ] ] 3 => array:7 [ "identificador" => "fig0020" "etiqueta" => "Figure 4" "tipo" => "MULTIMEDIAFIGURA" "mostrarFloat" => true "mostrarDisplay" => false "figura" => array:1 [ 0 => array:4 [ "imagen" => "gr4.jpeg" "Alto" => 1137 "Ancho" => 1771 "Tamanyo" => 357479 ] ] "descripcion" => array:1 [ "en" => "(a) General view to the NW of the study area, the escarp of the Balvanera fault can be observed, the trace of the El Nabo fault is inferred by a dotted line because it delimits NE side of the uplifted bock (horst). As reference the La Joya volcano is shown toward the NW of the study area; (b) close view in outcrop of the South Fault (pen for scale) showing the contact between the intrusive body at the left and evidence of hydrothermal activity; (c) outcrop of the intrusive body; (d) outcrop of a plane associated to the El Nabo fault zone, the contact with the intrusive body and the structure of a dyke with brecciated material are shown." ] ] 4 => array:7 [ "identificador" => "fig0025" "etiqueta" => "Figure 5" "tipo" => "MULTIMEDIAFIGURA" "mostrarFloat" => true "mostrarDisplay" => false "figura" => array:1 [ 0 => array:4 [ "imagen" => "gr5.jpeg" "Alto" => 1528 "Ancho" => 1657 "Tamanyo" => 308203 ] ] "descripcion" => array:1 [ "en" => "(a) Model box formed by an horizontal square grid of 25km2 and variable thickness layers. El Nabo and Mompani II wells were integrated in the model. Cells were grouped to represent hydro-stratigraphic layers with specific hydraulic properties. (b) Major geological structures in the model: (1) El Nabo fault-dyke, (2) Balvanera fault, (3) South fault and, (4) an intrusive body. The topographic surface is considered in the model." ] ] 5 => array:7 [ "identificador" => "fig0030" "etiqueta" => "Figure 6" "tipo" => "MULTIMEDIAFIGURA" "mostrarFloat" => true "mostrarDisplay" => false "figura" => array:1 [ 0 => array:4 [ "imagen" => "gr6.jpeg" "Alto" => 942 "Ancho" => 1378 "Tamanyo" => 236148 ] ] "descripcion" => array:1 [ "en" => "Piezometric levels measured in the Queretaro Valley for (a) 2004 and (b) 2006. The 2006 piezometric levels were used as initial condition for the model. Because the North-South hydraulic gradient, shown in the dashed square (modeled area), is persistent during 2004 and 2006 we considered steady state conditions." ] ] 6 => array:7 [ "identificador" => "fig0035" "etiqueta" => "Figure 7" "tipo" => "MULTIMEDIAFIGURA" "mostrarFloat" => true "mostrarDisplay" => false "figura" => array:1 [ 0 => array:4 [ "imagen" => "gr7.jpeg" "Alto" => 1528 "Ancho" => 1344 "Tamanyo" => 257995 ] ] "descripcion" => array:1 [ "en" => "Piezometric conditions shown in the model box with (a) border and (b) initial conditions after the steady state measured in 2006 (see <a class="elsevierStyleCrossRef" href="#fig0030">figure 6</a>b)." ] ] 7 => array:7 [ "identificador" => "fig0040" "etiqueta" => "Figure 8" "tipo" => "MULTIMEDIAFIGURA" "mostrarFloat" => true "mostrarDisplay" => false "figura" => array:1 [ 0 => array:4 [ "imagen" => "gr8.jpeg" "Alto" => 1198 "Ancho" => 1787 "Tamanyo" => 223018 ] ] "descripcion" => array:1 [ "en" => "Drawdown measurements collected during the El Nabo pumping test. Blue line indicates continuous datalogger measurements, blue squares the sampled points for calibration. Green line indicates increasing temperature in the extracted water." ] ] 8 => array:7 [ "identificador" => "fig0045" "etiqueta" => "Figure 9" "tipo" => "MULTIMEDIAFIGURA" "mostrarFloat" => true "mostrarDisplay" => false "figura" => array:1 [ 0 => array:4 [ "imagen" => "gr9.jpeg" "Alto" => 1160 "Ancho" => 1808 "Tamanyo" => 208188 ] ] "descripcion" => array:1 [ "en" => "Measured drawdown during the pumping test at El Nabo (blue squares) and calculated curve using the implemented model. The derivative of the data obtained by applying the fixed-end-point method (<a class="elsevierStyleCrossRef" href="#bib0040">Bourdet et al., 1989</a>) is presented in a orange line. Derivative drawdown was analyzed by the semi-analytical solution of <a class="elsevierStyleCrossRef" href="#bib0050">Butler and Liu (1991)</a> obtaining 1:2 and 1:1 slopes." ] ] 9 => array:7 [ "identificador" => "fig0050" "etiqueta" => "Figure 10" "tipo" => "MULTIMEDIAFIGURA" "mostrarFloat" => true "mostrarDisplay" => false "figura" => array:1 [ 0 => array:4 [ "imagen" => "gr10.jpeg" "Alto" => 1901 "Ancho" => 1028 "Tamanyo" => 319713 ] ] "descripcion" => array:1 [ "en" => "(a) 40Ar/39Ar* age spectra obtained for the biotite population 1 (small size) collected at different depths between 846 and 979 m in El Nabo well. The age obtained for this population is interpreted as evidence for a Miocene dyke intrusion. (b) 40Ar/39Ar* age spectra obtained for the biotite population 2 collected at the same depth interval of population 1. Although a reliable age cannot be assigned to the samples, a minimum age of 100 Ma is interpreted for the volcano-sedimentary sequence hosting the Miocene intrusion, (c) schematic section showing the intrusion of the volcanic dyke (Tmt AB, 10 Ma) within the Cretaceous sedimentary and volcanic units (host rocks: Ki Cz and Ki Ms V) that changes dramatically the hydraulic conductivity of the system by brecciation." ] ] 10 => array:7 [ "identificador" => "fig0055" "etiqueta" => "Figure 11" "tipo" => "MULTIMEDIAFIGURA" "mostrarFloat" => true "mostrarDisplay" => false "figura" => array:1 [ 0 => array:4 [ "imagen" => "gr11.jpeg" "Alto" => 1105 "Ancho" => 1877 "Tamanyo" => 328498 ] ] "descripcion" => array:1 [ "en" => "Model box of the study area presenting drawdown elevations and groundwater flow vectors after simulation of El Nabo pumping test using calibrated of hydraulic conductivity values shown in <a class="elsevierStyleCrossRef" href="#tbl0005">Table 1</a>. Faults act as flow channel interfering with the regional flow gradient and compartmentalizing the aquifer system. Arrows show the direction and sense of flow. The large intrusive dome modifies flow in the fault. The El Nabo fault-dyke exerts a major influence in groundwater flow." ] ] 11 => array:7 [ "identificador" => "fig0060" "etiqueta" => "Figure 12" "tipo" => "MULTIMEDIAFIGURA" "mostrarFloat" => true "mostrarDisplay" => false "figura" => array:1 [ 0 => array:4 [ "imagen" => "gr12.jpeg" "Alto" => 2376 "Ancho" => 1882 "Tamanyo" => 535765 ] ] "descripcion" => array:1 [ "en" => "Horizontal and vertical sections of the model box are presented in <a class="elsevierStyleCrossRef" href="#fig0055">Figure 11</a>. Colors indicate variation of piezometric heads from 1940 to 1820 masl. (a) Horizontal cross section at 1400 masl (location shown in <a class="elsevierStyleCrossRef" href="#fig0060">Figure 12</a>b), the black lines represent the analyzed faults walls and intrusive body. (b) and (c) are two cross sections of the model box in the NS and EW directions respectively." ] ] 12 => array:7 [ "identificador" => "fig0065" "etiqueta" => "Figure S1" "tipo" => "MULTIMEDIAFIGURA" "mostrarFloat" => true "mostrarDisplay" => false "figura" => array:1 [ 0 => array:4 [ "imagen" => "gr13.jpeg" "Alto" => 1272 "Ancho" => 1648 "Tamanyo" => 178229 ] ] "descripcion" => array:1 [ "en" => "Correlation diagrams obtained for samples for Biotite Population 1 and 2. Isochron lines (A and B) yield an age around 9.5 Ma for samples of Biotite Population 1; whereas Biotite Population 2 (C and D) converge at an age around 100 Ma." ] ] 13 => array:7 [ "identificador" => "tbl0005" "etiqueta" => "Table 1" "tipo" => "MULTIMEDIATABLA" "mostrarFloat" => true "mostrarDisplay" => false "tabla" => array:1 [ "tablatextoimagen" => array:1 [ 0 => array:2 [ "tabla" => array:1 [ 0 => """ <table border="0" frame="\n \t\t\t\t\tvoid\n \t\t\t\t" class=""><thead title="thead"><tr title="table-row"><th class="td" title="table-head " align="left" valign="top" scope="col" style="border-bottom: 2px solid black">Hydraulic conductivity (m/s) \t\t\t\t\t\t\n \t\t\t\t</th><th class="td" title="table-head " align="left" valign="top" scope="col" style="border-bottom: 2px solid black">E-W \t\t\t\t\t\t\n \t\t\t\t</th><th class="td" title="table-head " align="left" valign="top" scope="col" style="border-bottom: 2px solid black">Vertical \t\t\t\t\t\t\n \t\t\t\t</th></tr></thead><tbody title="tbody"><tr title="table-row"><td class="td" title="table-entry " align="left" valign="top">Rhyolite \t\t\t\t\t\t\n \t\t\t\t</td><td class="td" title="table-entry " align="left" valign="top">2.8E-05 \t\t\t\t\t\t\n \t\t\t\t</td><td class="td" title="table-entry " align="left" valign="top">6.1E-06 \t\t\t\t\t\t\n \t\t\t\t</td></tr><tr title="table-row"><td class="td" title="table-entry " align="left" valign="top">Conglomerate and sandstone \t\t\t\t\t\t\n \t\t\t\t</td><td class="td" title="table-entry " align="left" valign="top">3.4E-05 \t\t\t\t\t\t\n \t\t\t\t</td><td class="td" title="table-entry " align="left" valign="top">4.0E-05 \t\t\t\t\t\t\n \t\t\t\t</td></tr><tr title="table-row"><td class="td" title="table-entry " align="left" valign="top">Limestone and shale \t\t\t\t\t\t\n \t\t\t\t</td><td class="td" title="table-entry " align="left" valign="top">5.6E-02 \t\t\t\t\t\t\n \t\t\t\t</td><td class="td" title="table-entry " align="left" valign="top">1.3E-05 \t\t\t\t\t\t\n \t\t\t\t</td></tr><tr title="table-row"><td class="td" title="table-entry " align="left" valign="top">Basement Limestone \t\t\t\t\t\t\n \t\t\t\t</td><td class="td" title="table-entry " align="left" valign="top">9.9E-05 \t\t\t\t\t\t\n \t\t\t\t</td><td class="td" title="table-entry " align="left" valign="top">8.8E-06 \t\t\t\t\t\t\n \t\t\t\t</td></tr><tr title="table-row"><td class="td" title="table-entry " align="left" valign="top">Nabo-Fault NW-SE \t\t\t\t\t\t\n \t\t\t\t</td><td class="td" title="table-entry " align="left" valign="top">4.4E-07 \t\t\t\t\t\t\n \t\t\t\t</td><td class="td" title="table-entry " align="left" valign="top">4.9E-04 \t\t\t\t\t\t\n \t\t\t\t</td></tr><tr title="table-row"><td class="td" title="table-entry " align="left" valign="top">South fault NE-SW \t\t\t\t\t\t\n \t\t\t\t</td><td class="td" title="table-entry " align="left" valign="top">1.2E+02 \t\t\t\t\t\t\n \t\t\t\t</td><td class="td" title="table-entry " align="left" valign="top">6.1E-06 \t\t\t\t\t\t\n \t\t\t\t</td></tr><tr title="table-row"><td class="td" title="table-entry " align="left" valign="top">Balvanera-fault NW-SE \t\t\t\t\t\t\n \t\t\t\t</td><td class="td" title="table-entry " align="left" valign="top">5.1E-04 \t\t\t\t\t\t\n \t\t\t\t</td><td class="td" title="table-entry " align="left" valign="top">9.7E-06 \t\t\t\t\t\t\n \t\t\t\t</td></tr><tr title="table-row"><td class="td" title="table-entry " align="left" valign="top">Intrusive \t\t\t\t\t\t\n \t\t\t\t</td><td class="td" title="table-entry " align="left" valign="top">6.15E-05 \t\t\t\t\t\t\n \t\t\t\t</td><td class="td" title="table-entry " align="left" valign="top">1.39E-05 \t\t\t\t\t\t\n \t\t\t\t</td></tr><tr title="table-row"><td class="td" title="table-entry " align="left" valign="top">a) \t\t\t\t\t\t\n \t\t\t\t</td><td class="td" title="table-entry " align="" valign="top"> \t\t\t\t\t\t\n \t\t\t\t</td><td class="td" title="table-entry " align="" valign="top"> \t\t\t\t\t\t\n \t\t\t\t</td></tr><tr title="table-row"><td class="td" title="table-entry " align="left" valign="top">Sensitivity \t\t\t\t\t\t\n \t\t\t\t</td><td class="td" title="table-entry " align="left" valign="top">E-W \t\t\t\t\t\t\n \t\t\t\t</td><td class="td" title="table-entry " align="left" valign="top">Vertical \t\t\t\t\t\t\n \t\t\t\t</td></tr><tr title="table-row"><td class="td" title="table-entry " align="left" valign="top">Rhyolite \t\t\t\t\t\t\n \t\t\t\t</td><td class="td" title="table-entry " align="left" valign="top">4.0E-03 \t\t\t\t\t\t\n \t\t\t\t</td><td class="td" title="table-entry " align="left" valign="top">0.0E+00 \t\t\t\t\t\t\n \t\t\t\t</td></tr><tr title="table-row"><td class="td" title="table-entry " align="left" valign="top">Conglomerate and sandstone \t\t\t\t\t\t\n \t\t\t\t</td><td class="td" title="table-entry " align="left" valign="top">4.0E-03 \t\t\t\t\t\t\n \t\t\t\t</td><td class="td" title="table-entry " align="left" valign="top">4.0E-03 \t\t\t\t\t\t\n \t\t\t\t</td></tr><tr title="table-row"><td class="td" title="table-entry " align="left" valign="top">Limestone and shale \t\t\t\t\t\t\n \t\t\t\t</td><td class="td" title="table-entry " align="left" valign="top">1.2E-01 \t\t\t\t\t\t\n \t\t\t\t</td><td class="td" title="table-entry " align="left" valign="top">1.6E-02 \t\t\t\t\t\t\n \t\t\t\t</td></tr><tr title="table-row"><td class="td" title="table-entry " align="left" valign="top">Basement Limestone \t\t\t\t\t\t\n \t\t\t\t</td><td class="td" title="table-entry " align="left" valign="top">8.8E-03 \t\t\t\t\t\t\n \t\t\t\t</td><td class="td" title="table-entry " align="left" valign="top">1.3E-02 \t\t\t\t\t\t\n \t\t\t\t</td></tr><tr title="table-row"><td class="td" title="table-entry " align="left" valign="top">Nabo-Fault NW-SE \t\t\t\t\t\t\n \t\t\t\t</td><td class="td" title="table-entry " align="left" valign="top">1.0E+00 \t\t\t\t\t\t\n \t\t\t\t</td><td class="td" title="table-entry " align="left" valign="top">2.9E-01 \t\t\t\t\t\t\n \t\t\t\t</td></tr><tr title="table-row"><td class="td" title="table-entry " align="left" valign="top">South fault E-W \t\t\t\t\t\t\n \t\t\t\t</td><td class="td" title="table-entry " align="left" valign="top">1.2E-01 \t\t\t\t\t\t\n \t\t\t\t</td><td class="td" title="table-entry " align="left" valign="top">0.0E+00 \t\t\t\t\t\t\n \t\t\t\t</td></tr><tr title="table-row"><td class="td" title="table-entry " align="left" valign="top">Balvanera-fault NW-SE \t\t\t\t\t\t\n \t\t\t\t</td><td class="td" title="table-entry " align="left" valign="top">4.0E-03 \t\t\t\t\t\t\n \t\t\t\t</td><td class="td" title="table-entry " align="left" valign="top">0.0E+00 \t\t\t\t\t\t\n \t\t\t\t</td></tr><tr title="table-row"><td class="td" title="table-entry " align="left" valign="top">Intrusive \t\t\t\t\t\t\n \t\t\t\t</td><td class="td" title="table-entry " align="left" valign="top">0.0E+00 \t\t\t\t\t\t\n \t\t\t\t</td><td class="td" title="table-entry " align="left" valign="top">0.0E+00 \t\t\t\t\t\t\n \t\t\t\t</td></tr><tr title="table-row"><td class="td" title="table-entry " align="left" valign="top">b) \t\t\t\t\t\t\n \t\t\t\t</td><td class="td" title="table-entry " align="" valign="top"> \t\t\t\t\t\t\n \t\t\t\t</td><td class="td" title="table-entry " align="" valign="top"> \t\t\t\t\t\t\n \t\t\t\t</td></tr></tbody></table> """ ] "imagenFichero" => array:1 [ 0 => "xTab890448.png" ] ] ] ] "descripcion" => array:1 [ "en" => "(a) East-West and vertical hydraulic conductivities obtained from simulation. (b) Table of sensitivity variability for PEST analysis, higher sensitivity values of the model are related to the El Nabo fault. See text for details." ] ] ] "bibliografia" => array:2 [ "titulo" => "References" "seccion" => array:1 [ 0 => array:2 [ "identificador" => "bibs0005" "bibliografiaReferencia" => array:39 [ 0 => array:3 [ "identificador" => "bib0005" "etiqueta" => "Aguirre-Di¿az and Lo¿pez-Marti¿nez, 2001" "referencia" => array:1 [ 0 => array:2 [ "contribucion" => array:1 [ 0 => array:2 [ "titulo" => "The Amazcala caldera, Querétaro, Mexico. Geology and geochronology" "autores" => array:1 [ 0 => array:2 [ "etal" => false "autores" => array:2 [ 0 => "G.J. Aguirre-Di¿az" 1 => "M. Lo¿pez-Marti¿nez" ] ] ] ] ] "host" => array:1 [ 0 => array:1 [ "Revista" => array:5 [ "tituloSerie" => "Journal of Volcanology and Geothermal Research" "fecha" => "2001" "volumen" => "11" "paginaInicial" => "203" "paginaFinal" => "218" ] ] ] ] ] ] 1 => array:3 [ "identificador" => "bib0010" "etiqueta" => "Aguirre-Di¿az et al., 2005" "referencia" => array:1 [ 0 => array:2 [ "contribucion" => array:1 [ 0 => array:2 [ "titulo" => "Seismogenic Basin and Range and intra-arc normal faulting in the central Mexican Volcanic Belt, Querétaro México" "autores" => array:1 [ 0 => array:2 [ "etal" => false "autores" => array:3 [ 0 => "G.J. Aguirre-Di¿az" 1 => "J. Nieto-Obrego¿n" 2 => "R. Zu¿n¿iga" ] ] ] ] ] "host" => array:1 [ 0 => array:1 [ "Revista" => array:5 [ "tituloSerie" => "Geological Journal" "fecha" => "2005" "volumen" => "40" "paginaInicial" => "215" "paginaFinal" => "243" ] ] ] ] ] ] 2 => array:3 [ "identificador" => "bib0015" "etiqueta" => "Alaniz-Álvarez et al., 2001" "referencia" => array:1 [ 0 => array:1 [ "referenciaCompleta" => "Alaniz-Álvarez S.A., Nieto-Samaniego A.F., Reyes-Zaragoza M.A., Orozco-Esquivel M.T., Ojeda-García A.C., Vassallo F.L., 2001, Estratigrafía y deformación extensional en la región San Miguel de Allende-Querétaro, México (Stratigraphy and extensional deformation in the San Miguel de Allende-Queretaro region, Mexico). Revista Mexicana de Ciencias Geologicas, 18, 129-148." ] ] ] 3 => array:3 [ "identificador" => "bib0020" "etiqueta" => "Alaniz-Álvarez et al., 2002" "referencia" => array:1 [ 0 => array:1 [ "referenciaCompleta" => "Alaniz-Álvarez S.A., Nieto-Samaniego A.F., Orozco-Esquivel M.T., Vassallo F.L., Xu S., 2002, El sistema de fallas Taxco-San Miguel de Allende: implicaciones en la deformacio¿n post-eocénica del centro de México (The fault system Taxco-San Miguel de Allende: implications for post-Eocene deformation of central Mexico) Boletin de la Sociedad Geologica Mexicana, 55, 12-29." ] ] ] 4 => array:3 [ "identificador" => "bib0025" "etiqueta" => "Allen and Michel, 1999" "referencia" => array:1 [ 0 => array:2 [ "contribucion" => array:1 [ 0 => array:2 [ "titulo" => "Characterizing a Faulted Aquifer by Field Testing and Numerical Simulation" "autores" => array:1 [ 0 => array:2 [ "etal" => false "autores" => array:2 [ 0 => "D.M. Allen" 1 => "F.A. Michel" ] ] ] ] ] "host" => array:1 [ 0 => array:2 [ "doi" => "10.1111/j.1745-6584.1999.tb01164.x" "Revista" => array:6 [ "tituloSerie" => "Ground Water" "fecha" => "1999" "volumen" => "37" "paginaInicial" => "718" "paginaFinal" => "728" "link" => array:1 [ 0 => array:2 [ "url" => "https://www.ncbi.nlm.nih.gov/pubmed/19125925" "web" => "Medline" ] ] ] ] ] ] ] ] 5 => array:3 [ "identificador" => "bib0030" "etiqueta" => "Anderson, 2006" "referencia" => array:1 [ 0 => array:2 [ "contribucion" => array:1 [ 0 => array:2 [ "titulo" => "Analytical solutions for flow to a well through a fault" "autores" => array:1 [ 0 => array:2 [ "etal" => false "autores" => array:1 [ 0 => "E.I. Anderson" ] ] ] ] ] "host" => array:1 [ 0 => array:1 [ "Revista" => array:5 [ "tituloSerie" => "Advances in Water Resources" "fecha" => "2006" "volumen" => "29" "paginaInicial" => "1790" "paginaFinal" => "1803" ] ] ] ] ] ] 6 => array:3 [ "identificador" => "bib0035" "etiqueta" => "Bense et al., 2013" "referencia" => array:1 [ 0 => array:2 [ "contribucion" => array:1 [ 0 => array:2 [ "titulo" => "Fault zone hydrogeology" "autores" => array:1 [ 0 => array:2 [ "etal" => false "autores" => array:5 [ 0 => "V.F. Bense" 1 => "T. Gleeson" 2 => "S.E. Loveless" 3 => "O. Bour" 4 => "J. Scibek" ] ] ] ] ] "host" => array:1 [ 0 => array:1 [ "Revista" => array:5 [ "tituloSerie" => "Earth-Science Reviews" "fecha" => "2013" "volumen" => "127" "paginaInicial" => "171" "paginaFinal" => "192" ] ] ] ] ] ] 7 => array:3 [ "identificador" => "bib0040" "etiqueta" => "Bourdet et al., 1989" "referencia" => array:1 [ 0 => array:2 [ "contribucion" => array:1 [ 0 => array:2 [ "titulo" => "Use of pressure derivative in well-test interpretation" "autores" => array:1 [ 0 => array:2 [ "etal" => false "autores" => array:3 [ 0 => "D. Bourdet" 1 => "J.A. Ayoub" 2 => "Y.M. Pirard" ] ] ] ] ] "host" => array:1 [ 0 => array:1 [ "Revista" => array:4 [ "tituloSerie" => "SPE Formation Evaluation June" "fecha" => "1989" "paginaInicial" => "293" "paginaFinal" => "302" ] ] ] ] ] ] 8 => array:3 [ "identificador" => "bib0045" "etiqueta" => "Burbey, 2008" "referencia" => array:1 [ 0 => array:2 [ "contribucion" => array:1 [ 0 => array:2 [ "titulo" => "The Influence of Geologic Structures on Deformation due to Ground Water Withdrawal" "autores" => array:1 [ 0 => array:2 [ "etal" => false "autores" => array:1 [ 0 => "T.J. Burbey" ] ] ] ] ] "host" => array:1 [ 0 => array:2 [ "doi" => "10.1111/j.1745-6584.2007.00395.x" "Revista" => array:6 [ "tituloSerie" => "Ground Water" "fecha" => "2008" "volumen" => "46" "paginaInicial" => "202" "paginaFinal" => "211" "link" => array:1 [ 0 => array:2 [ "url" => "https://www.ncbi.nlm.nih.gov/pubmed/18194320" "web" => "Medline" ] ] ] ] ] ] ] ] 9 => array:3 [ "identificador" => "bib0050" "etiqueta" => "Butler and Liu, 1991" "referencia" => array:1 [ 0 => array:2 [ "contribucion" => array:1 [ 0 => array:2 [ "titulo" => "Pumping test in non-uniform aquifers - the linear strip case" "autores" => array:1 [ 0 => array:2 [ "etal" => false "autores" => array:2 [ 0 => "J.J. Butler" 1 => "W.Z. Liu" ] ] ] ] ] "host" => array:1 [ 0 => array:1 [ "Revista" => array:5 [ "tituloSerie" => "Journal of Hydrology" "fecha" => "1991" "volumen" => "128" "paginaInicial" => "69" "paginaFinal" => "99" ] ] ] ] ] ] 10 => array:3 [ "identificador" => "bib0055" "etiqueta" => "Caine et al., 1996" "referencia" => array:1 [ 0 => array:2 [ "contribucion" => array:1 [ 0 => array:2 [ "titulo" => "Fault zone architecture and permeability structure" "autores" => array:1 [ 0 => array:2 [ "etal" => false "autores" => array:3 [ 0 => "J.S. Caine" 1 => "J.P. Evans" 2 => "C.B. Forster" ] ] ] ] ] "host" => array:1 [ 0 => array:1 [ "Revista" => array:5 [ "tituloSerie" => "Geology" "fecha" => "1996" "volumen" => "24" "paginaInicial" => "1025" "paginaFinal" => "1028" ] ] ] ] ] ] 11 => array:3 [ "identificador" => "bib0060" "etiqueta" => "Carreón-Freyre et al., 2005" "referencia" => array:1 [ 0 => array:2 [ "contribucion" => array:1 [ 0 => array:2 [ "titulo" => "Influencia de la estratigrafía y estructura geológica en el flujo de agua subterránea del Valle de Querétaro (Influence of the stratigraphy and geologic structures on groundwater flow Querétaro Valley)" "autores" => array:1 [ 0 => array:2 [ "etal" => false "autores" => array:4 [ 0 => "D. Carreón-Freyre" 1 => "M. Cerca" 2 => "L. Luna-González" 3 => "F.J. Gámez-González" ] ] ] ] ] "host" => array:1 [ 0 => array:1 [ "Revista" => array:5 [ "tituloSerie" => "Revista Mexicana de Ciencias Geologicas" "fecha" => "2005" "volumen" => "22" "paginaInicial" => "1" "paginaFinal" => "18" ] ] ] ] ] ] 12 => array:3 [ "identificador" => "bib0065" "etiqueta" => "Cerca et al., 2007" "referencia" => array:1 [ 0 => array:1 [ "referenciaCompleta" => "Cerca M., Ferrari L., López-Martínez M., Martiny B., Iriondo A., 2007, Late Cretaceous shortening and early Tertiary shearing in the central Sierra Madre del Sur, southern Mexico: Insights into the evolution of the Caribbean–North American plate interaction: Tectonics, 26, TC3007, doi:10.1029/2006TC001981." ] ] ] 13 => array:3 [ "identificador" => "bib0070" "etiqueta" => "Conagua, 2003" "referencia" => array:1 [ 0 => array:2 [ "contribucion" => array:1 [ 0 => array:2 [ "titulo" => "Plan de manejo del Acuífero Valle de Querétaro" "autores" => array:1 [ 0 => array:2 [ "etal" => false "autores" => array:1 [ 0 => "Conagua" ] ] ] ] ] "host" => array:1 [ 0 => array:1 [ "Libro" => array:3 [ "fecha" => "2003" "editorial" => "Comisión Nacional del Agua" "editorialLocalizacion" => "México" ] ] ] ] ] ] 14 => array:3 [ "identificador" => "bib0075" "etiqueta" => "Delinom, 2009" "referencia" => array:1 [ 0 => array:2 [ "contribucion" => array:1 [ 0 => array:2 [ "titulo" => "Structural geology controls on groundwater flow: Lembang Fault case study, West Java" "autores" => array:1 [ 0 => array:2 [ "etal" => false "autores" => array:1 [ 0 => "R.M. Delinom" ] ] ] ] ] "host" => array:1 [ 0 => array:1 [ "Revista" => array:5 [ "tituloSerie" => "Indonesia. Hydrogeology Journal" "fecha" => "2009" "volumen" => "17" "paginaInicial" => "1011" "paginaFinal" => "1023" ] ] ] ] ] ] 15 => array:3 [ "identificador" => "bib0080" "etiqueta" => "Doherty, 2005" "referencia" => array:1 [ 0 => array:1 [ "referenciaCompleta" => "Doherty J., 2005, PEST Model-Independent Parameter Estimation User Manual: 5th Edition Watermark Numerical Computing, 336 pp." ] ] ] 16 => array:3 [ "identificador" => "bib0085" "etiqueta" => "Ferrari, 2000" "referencia" => array:1 [ 0 => array:2 [ "contribucion" => array:1 [ 0 => array:2 [ "titulo" => "Avances en el conocimiento de la Faja Volcánica Transmexicana durante la última década" "autores" => array:1 [ 0 => array:2 [ "etal" => false "autores" => array:1 [ 0 => "L. Ferrari" ] ] ] ] ] "host" => array:1 [ 0 => array:1 [ "Revista" => array:5 [ "tituloSerie" => "Boletin de la Sociedad Geologica Mexicana" "fecha" => "2000" "volumen" => "53" "paginaInicial" => "84" "paginaFinal" => "92" ] ] ] ] ] ] 17 => array:3 [ "identificador" => "bib0090" "etiqueta" => "Hall, 1981" "referencia" => array:1 [ 0 => array:2 [ "contribucion" => array:1 [ 0 => array:2 [ "titulo" => "The application of K-Ar and 40Ar/39Ar methods to the dating of recent volcanics and the Laschamp event. PhD Thesis" "autores" => array:1 [ 0 => array:2 [ "etal" => false "autores" => array:1 [ 0 => "C.M. Hall" ] ] ] ] ] "host" => array:1 [ 0 => array:1 [ "Libro" => array:3 [ "fecha" => "1981" "editorial" => "University of Toronto" "editorialLocalizacion" => "Canada" ] ] ] ] ] ] 18 => array:3 [ "identificador" => "bib0095" "etiqueta" => "Harbaugh et al., 2000" "referencia" => array:1 [ 0 => array:1 [ "referenciaCompleta" => "Harbaugh A.W., Banta E.R., Hill M.C., Mcdonald M.G., 2000, MODFLOW-2000, the U.S. Geological Survey modular ground-water model - User guide to modularization concepts and the Ground-Water Flow Process. U.S. Geological Survey Open-File Report 00-92, 121 pp." ] ] ] 19 => array:3 [ "identificador" => "bib0100" "etiqueta" => "Hsieh and Neuman, 1985" "referencia" => array:1 [ 0 => array:2 [ "contribucion" => array:1 [ 0 => array:2 [ "titulo" => "Field determi-nation of the three-dimensional hydraulic conductivity tensor of anisotropic media, 1 Theory" "autores" => array:1 [ 0 => array:2 [ "etal" => false "autores" => array:2 [ 0 => "P.A. Hsieh" 1 => "S.P. Neuman" ] ] ] ] ] "host" => array:1 [ 0 => array:1 [ "Revista" => array:5 [ "tituloSerie" => "Water Resources Research" "fecha" => "1985" "volumen" => "21" "paginaInicial" => "1655" "paginaFinal" => "1665" ] ] ] ] ] ] 20 => array:3 [ "identificador" => "bib0105" "etiqueta" => "Ingebritsen et al., 2006" "referencia" => array:1 [ 0 => array:1 [ "referenciaCompleta" => "Ingebritsen S.E., Sanford W.E., Neuzil C.E., 2006, Groundwater in Geologic Processes. 2nd Editio.: Cambridge University Press, 536 pp." ] ] ] 21 => array:3 [ "identificador" => "bib0110" "etiqueta" => "Jourde et al., 2007" "referencia" => array:1 [ 0 => array:2 [ "contribucion" => array:1 [ 0 => array:2 [ "titulo" => "Relationship between the geometrical and structural properties of layered fractured rocks and their effective permeability tensor. A simulation study" "autores" => array:1 [ 0 => array:2 [ "etal" => false "autores" => array:5 [ 0 => "H. Jourde" 1 => "P. Fenart" 2 => "M. Vinches" 3 => "S. Pistre" 4 => "B. Vayssade" ] ] ] ] ] "host" => array:1 [ 0 => array:1 [ "Revista" => array:5 [ "tituloSerie" => "Journal of Hydrology" "fecha" => "2007" "volumen" => "337" "paginaInicial" => "117" "paginaFinal" => "132" ] ] ] ] ] ] 22 => array:3 [ "identificador" => "bib0115" "etiqueta" => "Kulkarni et al., 2000" "referencia" => array:1 [ 0 => array:2 [ "contribucion" => array:1 [ 0 => array:2 [ "titulo" => "Hydrogeological framework for the Deccan basalt groundwater systems, west-central India" "autores" => array:1 [ 0 => array:2 [ "etal" => false "autores" => array:5 [ 0 => "H. Kulkarni" 1 => "S.B. Deolankar" 2 => "A. Lalwani" 3 => "B. Joseph" 4 => "S. Pawar" ] ] ] ] ] "host" => array:1 [ 0 => array:1 [ "Revista" => array:5 [ "tituloSerie" => "Hydrogeology Journal" "fecha" => "2000" "volumen" => "8" "paginaInicial" => "368" "paginaFinal" => "378" ] ] ] ] ] ] 23 => array:3 [ "identificador" => "bib0120" "etiqueta" => "Lopez and Smith, 1996" "referencia" => array:1 [ 0 => array:2 [ "contribucion" => array:1 [ 0 => array:2 [ "titulo" => "Fluid flow in fault zones: influence of hydraulic anysotropy and heterogeneity on the fluid flow and heat transfer regime" "autores" => array:1 [ 0 => array:2 [ "etal" => false "autores" => array:2 [ 0 => "D.A. Lopez" 1 => "L. Smith" ] ] ] ] ] "host" => array:1 [ 0 => array:1 [ "Revista" => array:5 [ "tituloSerie" => "Water Resources Research" "fecha" => "1996" "volumen" => "32" "paginaInicial" => "3227" "paginaFinal" => "3235" ] ] ] ] ] ] 24 => array:3 [ "identificador" => "bib0125" "etiqueta" => "Maclay and Small, 1983" "referencia" => array:1 [ 0 => array:2 [ "contribucion" => array:1 [ 0 => array:2 [ "titulo" => "Hydrostratigraphic subdivisions and fault barriers of the Edwards Aquifer, south-central Texas U.S.A" "autores" => array:1 [ 0 => array:2 [ "etal" => false "autores" => array:2 [ 0 => "R.W. Maclay" 1 => "T.A. Small" ] ] ] ] ] "host" => array:1 [ 0 => array:1 [ "Revista" => array:5 [ "tituloSerie" => "Journal of Hydrology" "fecha" => "1983" "volumen" => "61" "paginaInicial" => "127" "paginaFinal" => "146" ] ] ] ] ] ] 25 => array:3 [ "identificador" => "bib0130" "etiqueta" => "Mayer et al., 2007" "referencia" => array:1 [ 0 => array:2 [ "contribucion" => array:1 [ 0 => array:2 [ "titulo" => "Estimation of fault-zone conductance by calibration of a regional groundwater flow model: Desert Hot Springs California." "autores" => array:1 [ 0 => array:2 [ "etal" => false "autores" => array:4 [ 0 => "A. Mayer" 1 => "W. May" 2 => "C. Lukkarila" 3 => "J. Diehl" ] ] ] ] ] "host" => array:1 [ 0 => array:1 [ "Revista" => array:5 [ "tituloSerie" => "Hydrogeology Journal" "fecha" => "2007" "volumen" => "15" "paginaInicial" => "1093" "paginaFinal" => "1106" ] ] ] ] ] ] 26 => array:3 [ "identificador" => "bib0135" "etiqueta" => "Pacheco, 2002" "referencia" => array:1 [ 0 => array:2 [ "contribucion" => array:1 [ 0 => array:2 [ "titulo" => "Response to pumping of wells in sloping fault zone aquifers" "autores" => array:1 [ 0 => array:2 [ "etal" => false "autores" => array:1 [ 0 => "F.A.L. Pacheco" ] ] ] ] ] "host" => array:1 [ 0 => array:1 [ "Revista" => array:5 [ "tituloSerie" => "Journal of Hydrology" "fecha" => "2002" "volumen" => "259" "paginaInicial" => "116" "paginaFinal" => "135" ] ] ] ] ] ] 27 => array:3 [ "identificador" => "bib0140" "etiqueta" => "Rafini and Larocque, 2012" "referencia" => array:1 [ 0 => array:1 [ "referenciaCompleta" => "Rafini S., Larocque M., 2012, Numerical modeling of the hydraulic signatures of horizontal and inclined faults. Hydrogeology Journal. doi: 10.1007/s10040-011-0812-4." ] ] ] 28 => array:3 [ "identificador" => "bib0145" "etiqueta" => "Rathod and Rushton, 1991" "referencia" => array:1 [ 0 => array:2 [ "contribucion" => array:1 [ 0 => array:2 [ "titulo" => "Interpretation of Pumping from two-zone layered aquifers using a numerical model" "autores" => array:1 [ 0 => array:2 [ "etal" => false "autores" => array:2 [ 0 => "K.S. Rathod" 1 => "K.R. Rushton" ] ] ] ] ] "host" => array:1 [ 0 => array:1 [ "Revista" => array:5 [ "tituloSerie" => "Ground Water" "fecha" => "1991" "volumen" => "29" "paginaInicial" => "499" "paginaFinal" => "509" ] ] ] ] ] ] 29 => array:3 [ "identificador" => "bib0150" "etiqueta" => "Spane and Wurstner, 1993" "referencia" => array:1 [ 0 => array:2 [ "contribucion" => array:1 [ 0 => array:2 [ "titulo" => "DERIV: A computer program for calculating pressure derivatives for use in hydraulic test analysis" "autores" => array:1 [ 0 => array:2 [ "etal" => false "autores" => array:2 [ 0 => "F.A. Spane" 1 => "S.K. Wurstner" ] ] ] ] ] "host" => array:1 [ 0 => array:1 [ "Revista" => array:5 [ "tituloSerie" => "Ground Water" "fecha" => "1993" "volumen" => "31" "paginaInicial" => "814" "paginaFinal" => "822" ] ] ] ] ] ] 30 => array:3 [ "identificador" => "bib0155" "etiqueta" => "Steiger and Jäger, 1977" "referencia" => array:1 [ 0 => array:2 [ "contribucion" => array:1 [ 0 => array:2 [ "titulo" => "Sub commission on Geochronology: Convention on the use of decay constants in Geo and Cosmochronology" "autores" => array:1 [ 0 => array:2 [ "etal" => false "autores" => array:2 [ 0 => "R. Steiger" 1 => "E. Jäger" ] ] ] ] ] "host" => array:1 [ 0 => array:1 [ "Revista" => array:5 [ "tituloSerie" => "Earth and Planetary Science Letters" "fecha" => "1977" "volumen" => "36" "paginaInicial" => "359" "paginaFinal" => "362" ] ] ] ] ] ] 31 => array:3 [ "identificador" => "bib0160" "etiqueta" => "Tiab, 2005" "referencia" => array:1 [ 0 => array:2 [ "contribucion" => array:1 [ 0 => array:2 [ "titulo" => "Analysis of pressure derivative data of hydraulically fractured wells by the Tiab's Direct Synthesis technique" "autores" => array:1 [ 0 => array:2 [ "etal" => false "autores" => array:1 [ 0 => "D. Tiab" ] ] ] ] ] "host" => array:1 [ 0 => array:1 [ "Revista" => array:5 [ "tituloSerie" => "Journal of Petroleum Science and Engineering" "fecha" => "2005" "volumen" => "49" "paginaInicial" => "1" "paginaFinal" => "21" ] ] ] ] ] ] 32 => array:3 [ "identificador" => "bib0165" "etiqueta" => "Trejo-Moedano, 1989" "referencia" => array:1 [ 0 => array:1 [ "referenciaCompleta" => "Trejo-Moedano A., 1989, Estratigrafía y propiedades mecánicas del subsuelo del valle de la zona urbana de Querétaro (Stratigraphy and mechanical properties of the subsoil in the urban zone of Queretaro). Universidad Autónoma de Querétaro, México, 115 pp." ] ] ] 33 => array:3 [ "identificador" => "bib0170" "etiqueta" => "Valdez-Moreno et al., 1998" "referencia" => array:1 [ 0 => array:2 [ "contribucion" => array:1 [ 0 => array:2 [ "titulo" => "El volcán La Joya, estados de Querétaro y Guanajuato; un estratovolcán miocénico del Cinturón Volcánico Mexicano (The volcano La Joya, states of Queretaro and Guanajuato, a Miocene stratovolcano from Mexican Volcanic Belt)" "autores" => array:1 [ 0 => array:2 [ "etal" => false "autores" => array:3 [ 0 => "G. Valdez-Moreno" 1 => "G.J. Aguirre-Díaz" 2 => "M. López-Martínez" ] ] ] ] ] "host" => array:1 [ 0 => array:1 [ "Revista" => array:5 [ "tituloSerie" => "Revista Mexicana de Ciencias Geologicas" "fecha" => "1998" "volumen" => "15" "paginaInicial" => "181" "paginaFinal" => "197" ] ] ] ] ] ] 34 => array:3 [ "identificador" => "bib0175" "etiqueta" => "Walton, 1970" "referencia" => array:1 [ 0 => array:1 [ "referenciaCompleta" => "Walton W.C., 1970, Groundwater Resources Evaluation. Mc. Graw-Hill, New York. 185 pp." ] ] ] 35 => array:3 [ "identificador" => "bib0180" "etiqueta" => "Xu et al., 2011" "referencia" => array:1 [ 0 => array:2 [ "contribucion" => array:1 [ 0 => array:2 [ "titulo" => "Structural analysis of a relay ramp in the Queretaro Graben, central Mexico: Implications for relay ramp development" "autores" => array:1 [ 0 => array:2 [ "etal" => false "autores" => array:4 [ 0 => "S. Xu" 1 => "A.F. Nieto-Samaniego" 2 => "S.A. Alaniz-Alvarez" 3 => "L.M. Cerca-Martinez" ] ] ] ] ] "host" => array:1 [ 0 => array:1 [ "Revista" => array:5 [ "tituloSerie" => "Revista Mexicana de Ciencias Geologicas" "fecha" => "2011" "volumen" => "28" "paginaInicial" => "275" "paginaFinal" => "289" ] ] ] ] ] ] 36 => array:3 [ "identificador" => "bib0185" "etiqueta" => "Zeidouni, 2012" "referencia" => array:1 [ 0 => array:1 [ "referenciaCompleta" => "Zeidouni M., 2012, Analytical model of leakage through fault to overlying formations. Water Resources Research 48, W00N02, doi:10.1029/2012WR012582." ] ] ] 37 => array:3 [ "identificador" => "bib0190" "etiqueta" => "Zoback, 2007" "referencia" => array:1 [ 0 => array:1 [ "referenciaCompleta" => "Zoback M.D., 2007, Reservoir Geomechanics. Cambridge Press, 449 p." ] ] ] 38 => array:3 [ "identificador" => "bib0195" "etiqueta" => "Zuñiga et al., 2003" "referencia" => array:1 [ 0 => array:1 [ "referenciaCompleta" => "Zuñiga F.R, Pacheco J.F., Guzmán-Speziale M., Aguirre-Diaz G.J., Espindola V.H., Nava E., 2003, The Sanfandila earthquake sequence of 1998, Queretaro, México: activation of an undocumented fault in the northern edge of central Trans-Mexican Volcanic Belt. Tectonophysics, 361, 229-238. doi:10.1016/S0040-1951(02)00606-6." ] ] ] ] ] ] ] "lecturaRecomendada" => array:1 [ 0 => array:3 [ "vista" => "all" "titulo" => "Further reading" "seccion" => array:1 [ 0 => array:2 [ "vista" => "all" "bibliografiaReferencia" => array:2 [ 0 => array:3 [ "identificador" => "bib0200" "etiqueta" => "Lopez and Smith, 1995" "referencia" => array:1 [ 0 => array:2 [ "contribucion" => array:1 [ 0 => array:2 [ "titulo" => "Fluid flow in fault zones: analysis of theinterplay of convective circulation and topographically driven groundwater flow" "autores" => array:1 [ 0 => array:2 [ "etal" => false "autores" => array:2 [ 0 => "D.A. Lopez" 1 => "L. Smith" ] ] ] ] ] "host" => array:1 [ 0 => array:1 [ "Revista" => array:5 [ "tituloSerie" => "Water Resources Research" "fecha" => "1995" "volumen" => "31" "paginaInicial" => "1489" "paginaFinal" => "1503" ] ] ] ] ] ] 1 => array:3 [ "identificador" => "bib0205" "etiqueta" => "York et al., 2004" "referencia" => array:1 [ 0 => array:2 [ "contribucion" => array:1 [ 0 => array:2 [ "titulo" => "Unified equations for the slope, intercept, and standard errors of the best straight line" "autores" => array:1 [ 0 => array:2 [ "etal" => false "autores" => array:4 [ 0 => "D. York" 1 => "N. Evensen" 2 => "M. Lopez-Martinez" 3 => "D.J. De Basabe" ] ] ] ] ] "host" => array:1 [ 0 => array:1 [ "Revista" => array:5 [ "tituloSerie" => "American Journal of Physics" "fecha" => "2004" "volumen" => "73" "paginaInicial" => "367" "paginaFinal" => "375" ] ] ] ] ] ] ] ] ] ] ] "agradecimientos" => array:1 [ 0 => array:4 [ "identificador" => "xack186041" "titulo" => "Acknowledgements" "texto" => "The project FOMIX CONACYT-Querétaro, Ref. QRO-2004-C01-39, supported this research in the Valley of Queretaro. Comision Estatal de Aguas de Queretaro (CEAQ), the local water management agency, is gratefully acknowledged for funding the isotopic dating and for the technical help provided during pumping tests. F.J. Gamez-González (Conagua), I. Barrón-Medellin (CEAQ), and J. Díaz Escarcega (CEAQ) provided useful discussions for the groundwater flow interpretation. G. Aguirre-Díaz and an anonymous reviewer helped to clarify the manuscript. R. Carrizosa-González (UNAM), I. Ortiz-Villaseñor (CEAQ), and J. Perez-Villarreal helped in the field campaign and data compilation." "vista" => "all" ] ] ] "idiomaDefecto" => "en" "url" => "/00167169/0000005400000003/v2_201509120041/S0016716915000197/v2_201509120041/en/main.assets" "Apartado" => array:4 [ "identificador" => "40021" "tipo" => "SECCION" "en" => array:2 [ "titulo" => "Articles" "idiomaDefecto" => true ] "idiomaDefecto" => "en" ] "PDF" => "https://static.elsevier.es/multimedia/00167169/0000005400000003/v2_201509120041/S0016716915000197/v2_201509120041/en/main.pdf?idApp=UINPBA00004N&text.app=https://www.elsevier.es/" "EPUB" => "https://multimedia.elsevier.es/PublicationsMultimediaV1/item/epub/S0016716915000197?idApp=UINPBA00004N"]

Statistics

Gil Humberto Ochoa-Gonzáleza, Dora Carreón-Freyreb,

Corresponding author

freyre@geociencias.unam.mx

Corresponding author.

, Mariano Cercac, Margarita López-Martínezd

a Posgrado en Ciencias de la Tierra Centro de Geociencias Universidad Nacional Autónoma de México Campus Juriquilla Blvd. Juriquilla 3001 Col. Juriquilla, Querétaro, C.P. 76230 Querétaro, México

b Laboratorio de Mecánica de Geosistemas Centro de Geociencias Universidad Nacional Autónoma de México Campus Juriquilla Blvd. Juriquilla 3001 Col. Juriquilla, Querétaro, C.P. 76230 Querétaro, México

c Laboratorio de Mecánica de Geosistemas Centro de Geociencias Universidad Nacional Autónoma de México Campus Juriquilla Campus Juriquilla Blvd. Juriquilla 3001 Col. Juriquilla, Querétaro, C.P. 76230 Querétaro, México

d Departamento de Geología Centro de Investigación Científica y Educación Superior de Ensenada Carretera Ensenada-Tijuana No. 3918 Zona Playitas C.P. 22860 Ensenada, B.C., México

Read

6700

Times
was read the article

1566

Total PDF

5134

Total HTML

Share statistics

 array:22 [
  "pii" => "S0016716915000197"
  "issn" => "00167169"
  "doi" => "10.1016/j.gi.2015.04.016"
  "estado" => "S300"
  "fechaPublicacion" => "2015-07-01"
  "aid" => "16"
  "copyrightAnyo" => "2015"
  "documento" => "article"
  "crossmark" => 0
  "licencia" => "http://creativecommons.org/licenses/by-nc-nd/4.0/"
  "subdocumento" => "fla"
  "cita" => "Geofisica Internacional. 2015;54:199-220"
  "abierto" => array:3 [
    "ES" => true
    "ES2" => true
    "LATM" => true
  ]
  "gratuito" => true
  "lecturas" => array:2 [
    "total" => 2603
    "formatos" => array:3 [
      "EPUB" => 109
      "HTML" => 1763
      "PDF" => 731
    ]
  ]
  "itemSiguiente" => array:18 [
    "pii" => "S0016716915000215"
    "issn" => "00167169"
    "doi" => "10.1016/j.gi.2015.04.018"
    "estado" => "S300"
    "fechaPublicacion" => "2015-07-01"
    "aid" => "18"
    "documento" => "article"
    "crossmark" => 0
    "licencia" => "http://creativecommons.org/licenses/by-nc-nd/4.0/"
    "subdocumento" => "fla"
    "cita" => "Geofisica Internacional. 2015;54:221-44"
    "abierto" => array:3 [
      "ES" => true
      "ES2" => true
      "LATM" => true
    ]
    "gratuito" => true
    "lecturas" => array:2 [
      "total" => 3046
      "formatos" => array:3 [
        "EPUB" => 97
        "HTML" => 2236
        "PDF" => 713
      ]
    ]
    "en" => array:11 [
      "idiomaDefecto" => true
      "titulo" => "Mesozoic gliding and Tertiary basin and range tectonics in eastern Sonora&#44; Mexico"
      "tienePdf" => "en"
      "tieneTextoCompleto" => "en"
      "tieneResumen" => array:2 [
        0 => "es"
        1 => "en"
      ]
      "paginas" => array:1 [
        0 => array:2 [
          "paginaInicial" => "221"
          "paginaFinal" => "244"
        ]
      ]
      "contieneResumen" => array:2 [
        "es" => true
        "en" => true
      ]
      "contieneTextoCompleto" => array:1 [
        "en" => true
      ]
      "contienePdf" => array:1 [
        "en" => true
      ]
      "resumenGrafico" => array:2 [
        "original" => 0
        "multimedia" => array:7 [
          "identificador" => "fig0045"
          "etiqueta" => "Figure 9"
          "tipo" => "MULTIMEDIAFIGURA"
          "mostrarFloat" => true
          "mostrarDisplay" => false
          "figura" => array:1 [
            0 => array:4 [
              "imagen" => "gr9.jpeg"
              "Alto" => 1641
              "Ancho" => 1448
              "Tamanyo" => 866681
            ]
          ]
          "descripcion" => array:1 [
            "en" => "<p id="spar0055" class="elsevierStyleSimplePara elsevierViewall">Intense deformation &#40;a and b&#41; recorded in the base of the Ca&#241;ada de Tarachi Unit&#46; Note the scale of deformation&#46; Reverse and normal faults plus folds affect the sequence&#46; Z folds recorded in the Ca&#241;ada de Tarachi Unit &#40;c&#44; e&#44; d&#41; suggest transport to the south-west&#46; In a&#44; the bed where the man stand is the base of the detachment&#46; The S fold on f crops out along the Arivechi-Tarachi road&#44; also suggesting transport to the west&#46;</p>"
          ]
        ]
      ]
      "autores" => array:1 [
        0 => array:2 [
          "autoresLista" => "Jos&#233; Luis Rodr&#237;guez Casta&#241;eda, Jaime Rold&#225;n Quintana, Amabel Ortega Rivera"
          "autores" => array:3 [
            0 => array:2 [
              "nombre" => "Jos&#233; Luis Rodr&#237;guez"
              "apellidos" => "Casta&#241;eda"
            ]
            1 => array:2 [
              "nombre" => "Jaime Rold&#225;n"
              "apellidos" => "Quintana"
            ]
            2 => array:2 [
              "nombre" => "Amabel Ortega"
              "apellidos" => "Rivera"
            ]
          ]
        ]
      ]
    ]
    "idiomaDefecto" => "en"
    "EPUB" => "https://multimedia.elsevier.es/PublicationsMultimediaV1/item/epub/S0016716915000215?idApp=UINPBA00004N"
    "url" => "/00167169/0000005400000003/v2_201509120041/S0016716915000215/v2_201509120041/en/main.assets"
  ]
  "en" => array:19 [
    "idiomaDefecto" => true
    "titulo" => "Assessment of groundwater flow in volcanic faulted areas&#46; A study case in Queretaro&#44; Mexico"
    "tieneTextoCompleto" => true
    "paginas" => array:1 [
      0 => array:2 [
        "paginaInicial" => "199"
        "paginaFinal" => "220"
      ]
    ]
    "autores" => array:1 [
      0 => array:4 [
        "autoresLista" => "Gil Humberto Ochoa-Gonz&#225;lez, Dora Carre&#243;n-Freyre, Mariano Cerca, Margarita L&#243;pez-Mart&#237;nez"
        "autores" => array:4 [
          0 => array:3 [
            "nombre" => "Gil Humberto"
            "apellidos" => "Ochoa-Gonz&#225;lez"
            "referencia" => array:1 [
              0 => array:2 [
                "etiqueta" => "<span class="elsevierStyleSup">a</span>"
                "identificador" => "aff0005"
              ]
            ]
          ]
          1 => array:4 [
            "nombre" => "Dora"
            "apellidos" => "Carre&#243;n-Freyre"
            "email" => array:1 [
              0 => "freyre&#64;geociencias&#46;unam&#46;mx"
            ]
            "referencia" => array:2 [
              0 => array:2 [
                "etiqueta" => "<span class="elsevierStyleSup">b</span>"
                "identificador" => "aff0010"
              ]
              1 => array:2 [
                "etiqueta" => "<span class="elsevierStyleSup">&#42;</span>"
                "identificador" => "cor0005"
              ]
            ]
          ]
          2 => array:3 [
            "nombre" => "Mariano"
            "apellidos" => "Cerca"
            "referencia" => array:1 [
              0 => array:2 [
                "etiqueta" => "<span class="elsevierStyleSup">c</span>"
                "identificador" => "aff0015"
              ]
            ]
          ]
          3 => array:3 [
            "nombre" => "Margarita"
            "apellidos" => "L&#243;pez-Mart&#237;nez"
            "referencia" => array:1 [
              0 => array:2 [
                "etiqueta" => "<span class="elsevierStyleSup">d</span>"
                "identificador" => "aff0020"
              ]
            ]
          ]
        ]
        "afiliaciones" => array:4 [
          0 => array:3 [
            "entidad" => "Posgrado en Ciencias de la Tierra Centro de Geociencias Universidad Nacional Aut&#243;noma de M&#233;xico Campus Juriquilla Blvd&#46; Juriquilla 3001 Col&#46; Juriquilla&#44; Quer&#233;taro&#44; C&#46;P&#46; 76230 Quer&#233;taro&#44; M&#233;xico"
            "etiqueta" => "a"
            "identificador" => "aff0005"
          ]
          1 => array:3 [
            "entidad" => "Laboratorio de Mec&#225;nica de Geosistemas Centro de Geociencias Universidad Nacional Aut&#243;noma de M&#233;xico Campus Juriquilla Blvd&#46; Juriquilla 3001 Col&#46; Juriquilla&#44; Quer&#233;taro&#44; C&#46;P&#46; 76230 Quer&#233;taro&#44; M&#233;xico"
            "etiqueta" => "b"
            "identificador" => "aff0010"
          ]
          2 => array:3 [
            "entidad" => "Laboratorio de Mec&#225;nica de Geosistemas Centro de Geociencias Universidad Nacional Aut&#243;noma de M&#233;xico Campus Juriquilla Campus Juriquilla Blvd&#46; Juriquilla 3001 Col&#46; Juriquilla&#44; Quer&#233;taro&#44; C&#46;P&#46; 76230 Quer&#233;taro&#44; M&#233;xico"
            "etiqueta" => "c"
            "identificador" => "aff0015"
          ]
          3 => array:3 [
            "entidad" => "Departamento de Geolog&#237;a Centro de Investigaci&#243;n Cient&#237;fica y Educaci&#243;n Superior de Ensenada Carretera Ensenada-Tijuana No&#46; 3918 Zona Playitas C&#46;P&#46; 22860 Ensenada&#44; B&#46;C&#46;&#44; M&#233;xico"
            "etiqueta" => "d"
            "identificador" => "aff0020"
          ]
        ]
        "correspondencia" => array:1 [
          0 => array:3 [
            "identificador" => "cor0005"
            "etiqueta" => "&#8270;"
            "correspondencia" => "Corresponding author&#46;"
          ]
        ]
      ]
    ]
    "resumenGrafico" => array:2 [
      "original" => 0
      "multimedia" => array:7 [
        "identificador" => "fig0045"
        "etiqueta" => "Figure 9"
        "tipo" => "MULTIMEDIAFIGURA"
        "mostrarFloat" => true
        "mostrarDisplay" => false
        "figura" => array:1 [
          0 => array:4 [
            "imagen" => "gr9.jpeg"
            "Alto" => 1160
            "Ancho" => 1808
            "Tamanyo" => 208188
          ]
        ]
        "descripcion" => array:1 [
          "en" => "<p id="spar0055" class="elsevierStyleSimplePara elsevierViewall">Measured drawdown during the pumping test at El Nabo &#40;blue squares&#41; and calculated curve using the implemented model&#46; The derivative of the data obtained by applying the fixed-end-point method &#40;<a class="elsevierStyleCrossRef" href="#bib0040">Bourdet et al&#46;&#44; 1989</a>&#41; is presented in a orange line&#46; Derivative drawdown was analyzed by the semi-analytical solution of <a class="elsevierStyleCrossRef" href="#bib0050">Butler and Liu &#40;1991&#41;</a> obtaining 1&#58;2 and 1&#58;1 slopes&#46;</p>"
        ]
      ]
    ]
    "textoCompleto" => "<span class="elsevierStyleSections"><span id="sec0005" class="elsevierStyleSection elsevierViewall"><span class="elsevierStyleSectionTitle" id="sect0025">Introduction</span><p id="par0005" class="elsevierStylePara elsevierViewall">Understanding groundwater flow systems is crucial for the management of groundwater resources in Mexico because many areas in arid zones depend on groundwater as the main source for urban&#44; agricultural&#44; and industrial development&#46; In the Valley of Queretaro Aquifer &#40;VQA&#41;&#44; many wells are drilled near or along major faults and the high variability of piezometric levels in relatively short distances evidence changes on hydraulic properties &#40;<a class="elsevierStyleCrossRef" href="#bib0060">Carre&#243;n-Freyre <span class="elsevierStyleItalic">et al</span>&#46;&#44; 2005&#59;</a><a class="elsevierStyleCrossRef" href="#fig0005">Figure 1</a>&#41;&#46; Furthermore&#44; the spatial distribution of extraction wells and drawdown zones in the VQA during the last 30 years suggests that faults compartmentalize the aquifer system &#40;<a class="elsevierStyleCrossRef" href="#bib0060">Carre&#243;n-Freyre <span class="elsevierStyleItalic">et al</span>&#46;&#44; 2005</a>&#41;&#46; In the VQA&#44; groundwater is extracted from depths greater than two hundred meters&#59; at this depth&#44; geological heterogeneity should be considered for evaluating groundwater flow patterns to allow an accurate estimation of groundwater resources&#46; Thus&#44; the VQA could be considered a faulted aquifer&#44; where faults and fracture zones strongly influence the overall flow dynamics of the system &#40;<a class="elsevierStyleCrossRefs" href="#bib0055">Caine <span class="elsevierStyleItalic">et al</span>&#46;&#44; 1996&#59; Allen and Michel&#44; 1999</a>&#41; and the assessment of the hydrogeologic behavior of major structures can be as important as obtaining the hydraulic properties of the stratigraphic sequences&#46; This work examines the role of major structures &#40;e&#46;g&#46;&#44; faults and volcanic or plutonic bodies&#41; as hydrogeological independent units&#44; rather than only as hydraulic boundaries&#44; and presents a numerical modeling of the groundwater flow conditions in the faulted VQA&#46; Integration of representative geological data in the model included&#58; &#40;a&#41; vertical variations by assigning a hydraulic conductivity &#40;K&#41; to each geological unit&#59; &#40;b&#41; faults&#59; &#40;c&#41; lateral truncation of the hydro-stratigraphic units and&#59; &#40;d&#41; the presence of intrusive rocks &#40;<a class="elsevierStyleCrossRefs" href="#bib0175">Walton&#44; 1970&#59; Ingebritsen <span class="elsevierStyleItalic">et al</span>&#46;&#44; 2006&#59; Zoback&#44; 2007</a>&#41;&#44; in particular dikes &#40;<a class="elsevierStyleCrossRef" href="#bib0115">Kulkarni <span class="elsevierStyleItalic">et al</span>&#46;&#44; 2000</a>&#41;&#46;</p><elsevierMultimedia ident="fig0005"></elsevierMultimedia><p id="par0010" class="elsevierStylePara elsevierViewall">The study area is located in the northwestern part of the VQA system &#40;<a class="elsevierStyleCrossRef" href="#fig0005">Figure 1</a>&#41;&#46; The Valley of Queretaro is located in the central part of the Miocene to Recent Trans Mexican Volcanic Belt &#40;TMVB&#59; <a class="elsevierStyleCrossRefs" href="#bib0085">Ferrari&#44; 2000&#59; Alaniz-Alvarez <span class="elsevierStyleItalic">et al</span>&#46;&#44; 2002</a>&#41; and is delimited by groups of N-S and NE-SW trending faults active at least since the Miocene &#40;<a class="elsevierStyleCrossRef" href="#bib0015">Alaniz-Alvarez <span class="elsevierStyleItalic">et al</span>&#46;&#44; 2001</a>&#41;&#46; We aim at having a better understanding on the role of major structures in groundwater flow patterns in the VQA by integrating them in a numerical model&#46;</p></span><span id="sec0010" class="elsevierStyleSection elsevierViewall"><span class="elsevierStyleSectionTitle" id="sect0030">Faulted aquifer systems</span><p id="par0015" class="elsevierStylePara elsevierViewall">Previous works have documented the variations in hydraulics properties induced by the complex architecture of regional faults &#40;<a class="elsevierStyleCrossRef" href="#bib0125">MaClay and Small&#44; 1983</a>&#59; L&#243;pez and Smith&#44; 1995&#59; L&#243;pez and Smith&#44; 1996&#59; <a class="elsevierStyleCrossRefs" href="#bib0055">Caine <span class="elsevierStyleItalic">et al</span>&#46;&#44; 1996&#59; Allen and Michel&#44; 1999&#59; Carre&#243;n-Freyre <span class="elsevierStyleItalic">et al</span>&#46;&#44; 2005&#59; Anderson&#44; 2006&#59; Mayer <span class="elsevierStyleItalic">et al</span>&#46;&#44; 2007&#59; Delinom&#44; 2009&#59; Rafini and Larocque&#44; 2012</a>&#41;&#46; Anisotropy of hydraulic conductivity in faulted and stratified geological media determines groundwater flow patterns &#40;<a class="elsevierStyleCrossRef" href="#bib0100">Hsieh and Neuman&#44; 1985</a>&#41;&#46; Regional faults represent an important heterogeneity influencing overall aquifer dynamics&#44; not only determining preferential flow directions but also modifying their storage capacity &#40;<a class="elsevierStyleCrossRefs" href="#bib0135">Pacheco&#44; 2002&#59; Burbey&#44; 2008&#59; Bense <span class="elsevierStyleItalic">et al</span>&#46;&#44; 2013</a>&#41;&#46; Lithological changes and fractures in volcanic rocks&#44; on the other hand&#44; modify hydrogeological properties of geological units&#46; A fault can be a preferential conduit for vertical flow&#44; i&#46;e&#46;&#44; a &#8220;leaky fault&#8221;&#44; when the damaged zone is well developed and is laterally sealed &#40;<a class="elsevierStyleCrossRef" href="#bib0185">Zeidouni&#44; 2012</a>&#41;&#46; Fault zones can also behave as leakage areas between aquifers by juxtaposing aquifers that would otherwise be separated by impermeable layers&#46;</p><p id="par0020" class="elsevierStylePara elsevierViewall">Major faults can behave as flow barriers or preferential channels depending on the hydraulic properties developed during deformation history &#40;<a class="elsevierStyleCrossRefs" href="#bib0125">MaClay and Small&#44; 1983&#59; Caine <span class="elsevierStyleItalic">et al</span>&#46;&#44; 1996&#59; Carre&#243;n-Freyre <span class="elsevierStyleItalic">et al</span>&#46;&#44; 2005&#59; Mayer <span class="elsevierStyleItalic">et al</span>&#46;&#44; 2007&#59; Delinom&#44; 2009</a>&#44; <a class="elsevierStyleCrossRef" href="#bib0140">Rafini and Larocque&#44; 2012</a>&#41;&#44; or they can only separate units of different hydraulic conductivity &#40;<a class="elsevierStyleCrossRefs" href="#bib0025">Allen and Michel&#44; 1999&#59; Mayer <span class="elsevierStyleItalic">et al</span>&#46;&#44; 2007</a>&#41;&#59; for instance&#44; the vertical displacement of a fault can juxtapose permeable and impermeable layers in the direction perpendicular to the fault plane&#46; In early models&#44; faults were commonly considered as border conditions when assessing hydraulic parameters&#44; even in stratified leaky systems &#40;<a class="elsevierStyleCrossRef" href="#bib0145">Rathod and Rushton&#44; 1991</a>&#41;&#46;</p><p id="par0025" class="elsevierStylePara elsevierViewall">The style of deformation and rupture of a fault zone strongly influences variations of its hydraulic properties&#46; <a class="elsevierStyleCrossRef" href="#bib0055">Caine <span class="elsevierStyleItalic">et al</span>&#46; &#40;1996&#41;</a> proposed four conceptual schemes for fault-related fluid-flow according to the distribution of conduits and barriers &#40;localized and distributed&#41; within fault plain&#46; <a class="elsevierStyleCrossRef" href="#bib0030">Anderson &#40;2006&#41;</a> identified three essential hydraulic elements of a fault plane&#58; &#40;1&#41; a low hydraulic conductivity core characterized by fault breccias&#44; cataclastic rocks&#44; and gouge&#59; &#40;2&#41; an adjacent damaged zone characterized by an increase in the hydraulic conductivity caused by brecciation and fracturing and&#59; &#40;3&#41; fresh unaltered rock with a characteristic hydraulic conductivity that can change in both sides of the fault because of the fault displacement&#46;</p><p id="par0030" class="elsevierStylePara elsevierViewall">A decrease in the hydraulic conductivity within a fault plane &#40;flow barrier&#41; can be associated with &#40;<a class="elsevierStyleCrossRefs" href="#bib0130">Mayer <span class="elsevierStyleItalic">et al</span>&#46;&#44; 2007&#59; MaClay and Small&#44; 1983</a>&#41;&#58; &#40;1&#41; cataclasis or grain size reduction&#59; &#40;2&#41; juxtaposition of permeable and impermeable layers&#44; displacements along the fault plane can decrease the possibility of flow through the fault and&#47;or deviating flow along the fault&#59; &#40;3&#41; rotation of plane and elongate detritus that became parallel to fault plane reducing the hydraulic conductivity in the orientation perpendicular to the plane&#59; &#40;4&#41; presence of sedimentary materials with low hydraulic conductivity &#40;e&#46; g&#46;&#44; clays&#41; in the fault zone&#59; and &#40;5&#41; chemical precipitation of mineral in the fault zone reducing permeability&#46;</p><p id="par0035" class="elsevierStylePara elsevierViewall">Furthermore&#44; major faults often present associate or secondary planes and fractures that may cause contrasting hydraulic behaviors&#46; In anisotropic media&#44; each geological unit is represented by a hydraulic conductivity tensor considering three dimensional variations of hydraulic conductivity&#46; In fractured rocks the tensor can align following the preferential orientation of fractures &#40;<a class="elsevierStyleCrossRef" href="#bib0110">Jourde <span class="elsevierStyleItalic">et al</span>&#46;&#44; 2007</a>&#41;&#46;</p></span><span id="sec0015" class="elsevierStyleSection elsevierViewall"><span class="elsevierStyleSectionTitle" id="sect0035">Geology of the Study Area &#40;Northwest of the Valley of Queretaro Aquifer&#44; VQA&#41;</span><p id="par0040" class="elsevierStylePara elsevierViewall">The regional geology of the VQA has been extensively studied in the last decades &#40;<a class="elsevierStyleCrossRefs" href="#bib0165">Trejo-Moedano&#44; 1989&#59; Valdez-Moreno <span class="elsevierStyleItalic">et al</span>&#46;&#44; 1998&#59; Aguirre-D&#237;az and Lopez-Martinez&#44; 2001</a>&#59; Alan&#237;z-Alvarez <span class="elsevierStyleItalic">et al</span>&#46;&#44; 2001&#44; 2002&#59; <a class="elsevierStyleCrossRefs" href="#bib0060">Carre&#243;n-Freyre <span class="elsevierStyleItalic">et al</span>&#46;&#44; 2005&#59; Aguirre-D&#237;az <span class="elsevierStyleItalic">et al</span>&#46;&#44; 2005&#59; Xu <span class="elsevierStyleItalic">et al</span>&#46;&#44; 2011</a>&#41;&#46; The Valley of Queretaro is a basin delimited by groups of N-S and NE-SW trending faults and volcanoes&#44; located in the central part of the Trans Mexican Volcanic Belt &#40;TMVB&#59; <a class="elsevierStyleCrossRefs" href="#bib0085">Ferrari&#44; 2000&#59; Alaniz-Alvarez <span class="elsevierStyleItalic">et al</span>&#46;&#44; 2002</a>&#59; <a class="elsevierStyleCrossRef" href="#fig0005">Figure 1</a>&#41;&#46; The N-S trending west-dipping Central fault is the eastern limit of the regional Queretaro graben &#40;<a class="elsevierStyleCrossRef" href="#fig0005">Figure 1</a>&#41;&#46; The other major faulting family &#40;trending NE-SW&#41; produced a nearly orthogonal pattern with the N-S system that has resulted in a mosaic formed by horsts&#44; grabens&#44; and half-grabens &#40;Alan&#237;z-Alvarez <span class="elsevierStyleItalic">et al</span>&#46;&#44; 2001&#44; 2002&#59; <a class="elsevierStyleCrossRefs" href="#bib0060">Carre&#243;n-Freyre <span class="elsevierStyleItalic">et al</span>&#46;&#44; 2005&#59; Aguirre-D&#237;az <span class="elsevierStyleItalic">et al</span>&#46;&#44; 2005&#59; Xu <span class="elsevierStyleItalic">et al</span>&#46;&#44; 2011</a>&#41; with varying vertical displacements&#44; in some areas up to 400 m&#46; Both fault systems have been active at least since the Miocene &#40;e&#46; g&#46;&#44; <a class="elsevierStyleCrossRefs" href="#bib0015">Alaniz-Alvarez <span class="elsevierStyleItalic">et al</span>&#46;&#44; 2001&#59; Zu&#241;iga <span class="elsevierStyleItalic">et al</span>&#46;&#44; 2003&#59; Aguirre-D&#237;az <span class="elsevierStyleItalic">et al</span>&#46;&#44; 2005</a>&#41; and both controlled the deposition of sediments and volcanic rocks in the VQA&#46;</p><p id="par0045" class="elsevierStylePara elsevierViewall">This study is focused on the northwestern part of the VQA &#40;dashed square in <a class="elsevierStyleCrossRefs" href="#fig0005">Figures 1 and 2</a>&#41;&#44; a zone of groundwater recharge because of faulting &#40;<a class="elsevierStyleCrossRef" href="#bib0070">Conagua&#44; 2003</a>&#41;&#46; The available information was compiled and integrated in a geological map with the hydrogeologic information of the VQA and surroundings&#46; The local stratigraphy was established from the reinterpretation of well logs&#44; <span class="elsevierStyleSup">40</span>Ar&#42;&#47;<span class="elsevierStyleSup">39</span>Ar dating of magmatic rocks&#44; and geologic mapping&#46; The structural and stratigraphic setting of this part of the basin can be defined with sufficient detail and indicates that the geological evolution is more complex than hitherto recognized &#40;<a class="elsevierStyleCrossRef" href="#fig0010">Figure 2</a>a&#41;&#46; In this area&#44; the Cenozoic volcanic and colluvial sequences cover partially the Cretaceous units that were uplifted by a shortening episode and thus are well exposed in the hanging wall&#46; A schematic W-E geological cross section AA&#8217; was constructed using the lithological logs of water extraction wells available in the zone &#40;Mompani II&#44; El Nabo&#44; UNAM&#44; and El Salitre&#41; &#40;<a class="elsevierStyleCrossRef" href="#fig0010">Figure 2</a>&#41;&#46;</p><elsevierMultimedia ident="fig0010"></elsevierMultimedia><p id="par0050" class="elsevierStylePara elsevierViewall">The oldest sequence &#40;Cretaceous&#41; consists of volcano-sedimentary rocks&#44; documented here for the first time in the Valley of Queretaro &#40;Ki MsV&#44; see below&#41;&#44; a thick Albian limestone and shale sequence &#40;Ki Cz&#44; of 800 m thick&#41;&#44; and a late Cretaceous shale and sandstone interbedded with argillaceous limestone &#40;Ki Cz-Ar&#41;&#46; The Cretaceous sequence is pervasively folded by shortening deformation&#46;</p><p id="par0055" class="elsevierStylePara elsevierViewall">The Cenozoic sequence consists of a continental sedimentary unit with predominance of conglomerates and sandstone &#40;Te Cg-Ar&#41;&#44; overlain by a volcanic sequence with pyroclastic and sedimentary interbedded deposits&#46; The Cenozoic sequence is fractured and affected by extensional faults&#46; The Te Cg-Ar layer of probable Eocene age is present at the Mompani II well &#40;at a depth of 200 m&#41;&#44; and at the bottom of UNAM-Juriquilla and El Salitre wells &#40;<a class="elsevierStyleCrossRef" href="#fig0015">Figure 3</a>a&#41;&#46; The volcanic rock sequence includes lavas such as the Salitre Andesite &#40;ToA&#41; and silicic domes &#40;Tom Da&#59; Juriquilla dacite&#59; <a class="elsevierStyleCrossRef" href="#bib0020">Alaniz-Alvarez <span class="elsevierStyleItalic">et al</span>&#46;&#44; 2002</a>&#41; of the Oligocene and early Miocene&#46; An intrusive body of Oligocene age &#40;ca&#46; 30 Ma&#41; is exposed in the zone of Juriquilla within the study area &#40;<a class="elsevierStyleCrossRef" href="#fig0010">Figure 2</a>a&#41;&#46; The pyroclastic and lacustrine sequences in the Queretaro area consist of tuffs and other pyroclastic deposits of 80 m in average of thickness &#40;Tom PyLac&#41;&#59; and this unit includes the Ezequiel Montes Pumice&#44; an important stratigraphic marker dated between 7&#46;5 and 5&#46;6 Ma &#40;<a class="elsevierStyleCrossRef" href="#bib0005">Aguirre-D&#237;az and L&#243;pez-Mart&#237;nez&#44; 2001</a>&#41;&#46; The coarse facies of this sequence &#40;Tom PyLac&#41; was the main hydrogeological unit exploited in the regional VQA&#44; but is scarcely exposed in the study area&#44; it is found in the upper part of the Mompani II log &#40;<a class="elsevierStyleCrossRef" href="#fig0015">Figure 3</a>&#41;&#46; Fissure basalts and lava flows &#40;Tmt AB&#41;&#44; ranging in age from 12 to 5 Ma &#40;<a class="elsevierStyleCrossRef" href="#bib0005">Aguirre-D&#237;az and L&#243;pez-Mart&#237;nez&#44; 2001</a>&#41; erupted from nearby volcanoes&#46; The largest volcanic edifice&#44; located 13<span class="elsevierStyleHsp" style=""></span>km to the NW from the studied area&#44; is the La Joya Volcano ranging in age from 12 to 10 Ma &#40;<a class="elsevierStyleCrossRef" href="#bib0170">Valdez-Moreno <span class="elsevierStyleItalic">et al</span>&#46;&#44; 1998</a>&#41;&#46;</p><elsevierMultimedia ident="fig0015"></elsevierMultimedia><p id="par0060" class="elsevierStylePara elsevierViewall">The truncation and throw of the N-S trending extensional faults was inferred by the elevation differences among stratigraphic levels from the correlation of lithological logs presented in <a class="elsevierStyleCrossRef" href="#fig0015">Figure 3</a>&#46; The vertical displacements are smaller than 200 m&#46; From <a class="elsevierStyleCrossRef" href="#fig0010">Figure 2</a>b note that two additional buried faults nearby&#44; and to west of UNAM well&#44; were inferred in this work&#46; Structurally&#44; the study area is a block delimited by the major N-S trending <span class="elsevierStyleItalic">Balvanera</span> and <span class="elsevierStyleItalic">El Nabo</span> fault zones dipping to the northeast &#40;<a class="elsevierStyleCrossRef" href="#fig0020">Figure 4</a>a&#41;&#46; The South fault and other NE-SW trending fault segments delineate the limits of the block toward the south&#46; Both faults systems have nearly vertical fault planes and their kinematics indicate dip-slip towards the center of the Queretaro Valley since Miocene time&#46; These major faults include several parallel fault planes affecting areas with more than 50 m in width &#40;<a class="elsevierStyleCrossRef" href="#fig0020">Figure 4</a>&#41;&#46; Faults are partially exposed and their trace length was inferred by topography&#44; field observations&#44; and lithological correlations of well logs&#46; The relative age of faulting in the region was analyzed in detail by <a class="elsevierStyleCrossRef" href="#bib0015">Alaniz-Alvarez <span class="elsevierStyleItalic">et al</span> &#40;2001&#41;</a>&#59; the older system corresponds to the South Fault that favors the emplacement of the Oligocene intrusive body&#46; The Balvanera and El Nabo fault systems are younger and favor the Miocene volcanism through dykes emplaced along the faults planes &#40;<a class="elsevierStyleCrossRef" href="#fig0020">Figure 4</a>&#41;&#46; The <span class="elsevierStyleItalic">El Nabo Fault</span> is of particular relevance to this work&#44; it is a nearly vertical fault partially exposed to the north of the study area where plateau lavas fill the hanging wall of the fault &#40;solid line with thick marks in <a class="elsevierStyleCrossRef" href="#fig0010">Figure 2</a>a&#41;&#46; Toward its southern part the fault is not clearly observed on surface but it can be inferred from the topographic relief and the interpretation of the stratigraphic records &#40;dotted line in <a class="elsevierStyleCrossRef" href="#fig0010">Figure 2</a>a&#41;&#46; A dyke emplaced along a fault plane parallel to El Nabo fault was observed at the surface northeastward of the El Nabo well &#40;<a class="elsevierStyleCrossRef" href="#fig0020">Figure 4</a>b&#41;&#46; All this information was integrated to constrain the numerical model&#46;</p><elsevierMultimedia ident="fig0020"></elsevierMultimedia></span><span id="sec0020" class="elsevierStyleSection elsevierViewall"><span class="elsevierStyleSectionTitle" id="sect0040">Methodology</span><span id="sec0025" class="elsevierStyleSection elsevierViewall"><span class="elsevierStyleSectionTitle" id="sect0045">Stratigraphic correlation</span><p id="par0065" class="elsevierStylePara elsevierViewall">The conceptual model of the study area is based on geological units with contrasting hydraulic properties&#44; granular and fractured rock layers&#44; truncated by major faults and volcanic edifices&#47;bodies that may interrupt or favor lateral flows&#46; <a class="elsevierStyleCrossRef" href="#bib0060">Carre&#243;n-Freyre <span class="elsevierStyleItalic">et al</span>&#46; &#40;2005&#41;</a> correlated the hydrogeological units with lithostratigraphic units to describe the VQA and proposed a conceptual model of multilayer aquifer with groundwater flow varying from local to regional&#46; The stratigraphic correlation of lithological logs for the Cenozoic was based on the reconnaissance of an andesitic-basaltic unit of late Miocene age&#44; named Tmt AB&#44; as an index layer&#46; For instance&#44; as shown in the record of the Ejido Jurica well &#40;<a class="elsevierStyleCrossRef" href="#fig0015">Figure 3</a>a&#41; the Tmt AB layer commonly is overlaid by a Miocene-Pliocene sedimentary unit&#44; named Tmp Cg-Ar-Lm&#44; and underlied by an Oligocene-Miocene pyroclastic and lacustrine unit&#44; named Tom PyLac&#46;</p><p id="par0070" class="elsevierStylePara elsevierViewall">The Cretaceous units are more relevant to groundwater flow in this area and include sandstone and argillaceous sequences&#44; limestone rocks&#44; and a volcano-sedimentary sequence&#46; These units were found during drilling of <span class="elsevierStyleItalic">El Nabo</span> well&#44; which is the deepest in the study area and is located within a block where Tertiary sequences are partially eroded &#40;<a class="elsevierStyleCrossRef" href="#fig0015">Figure 3</a>b&#41;&#46; The contrasting stratigraphy and the vertical difference in elevation of the Ki Cz unit from 1680 to 1100 masl among the Mompani II and El Nabo records separated by approximately one kilometer suggested the presence of the El Nabo Fault at this location&#46;</p><p id="par0075" class="elsevierStylePara elsevierViewall">Drilling of the first 600 m depth of <span class="elsevierStyleItalic">El Nabo</span> well initiated in 2000 at this time the production of water was approximately 6 l&#47;s&#46; Extended drilling in <span class="elsevierStyleItalic">El Nabo</span> well from 600 to 1000 m depth was performed in 2007&#44; the production of the well increased dramatically to 47 l&#47;s after drilling at 850 m depth&#44; that coincided with a lithological change from limestone to a volcano-sedimentary sequence &#40;volcanic rocks and limestone&#41; bearing biotite minerals&#41;&#46; Further evidence of the fault and its relevance to groundwater flow was obtained by the geochronological data of rock samples from the El Nabo well presented below&#46;</p></span><span id="sec0030" class="elsevierStyleSection elsevierViewall"><span class="elsevierStyleSectionTitle" id="sect0050">Biotite age determination</span><p id="par0080" class="elsevierStylePara elsevierViewall">In order to define the age of the basal volcanic unit&#44; rock samples collected each 3 m during drilling from 800 to 1050 m depth were classified by the use of a 20x-40x stereoscopic magnifier&#44; and some samples of interest were selected for thin section description&#46; Biotite concentrates were selected for dating at depths between 850 and 950 m&#46; These biotite concentrates were analyzed by the <span class="elsevierStyleSup">40</span>Ar&#42;&#47;<span class="elsevierStyleSup">39</span>Ar technique at the Geochronology Lab in CICESE using a VG5400 mass spectrometer and a laser extraction system&#46; The analytical details of the method can be consulted elsewhere &#40;e&#46;g&#46;&#44; <a class="elsevierStyleCrossRef" href="#bib0065">Cerca <span class="elsevierStyleItalic">et al</span>&#46;&#44; 2007</a>&#41;&#46; Mineral concentrates were previously irradiated at the nuclear reactor of the University of McMaster in Hamilton&#44; Ontario&#44; Canada&#44; using as monitor of irradiation two known-age samples of sanidine&#44; FCT-2 of 27&#46;84 &#177; 0&#46;04 Ma and TCR-2 of 27&#46;87 &#177; 0&#46;04 Ma&#46; All the experiments were made by step-heating by increasing the laser power from 0&#46;2 to 8 watts and measuring the isotopic composition of argon released during each step&#46; The reported age is calculated for each individual fraction and the integrated results adding individual steps in accordance to <a class="elsevierStyleCrossRef" href="#bib0090">Hall &#40;1981&#41;</a>&#44; and considering the corrections by discrimination&#44; radioactive decay &#40;<a class="elsevierStyleCrossRef" href="#bib0155">Steiger &#38; J&#228;ger&#44; 1977</a>&#41;&#44; the line background&#44; and interference reactions from isotopes derived from Ca&#44; Cl&#44; and K&#46; Results were grouped according to depth and size of the samples and are presented in a graph named age spectra&#46; Correlation diagrams in which the relation &#40;<span class="elsevierStyleSup">40</span>Ar&#42;&#47;<span class="elsevierStyleSup">39</span>Ar<span class="elsevierStyleInf">K</span>&#41; is given by the inverse of the intercept of the line defined by the data and is equivalent to the age in Ma are presented in a supplementary Figure &#40;S1&#41;&#46;</p></span><span id="sec0035" class="elsevierStyleSection elsevierViewall"><span class="elsevierStyleSectionTitle" id="sect0055">Integration of conceptual and numerical models</span><p id="par0085" class="elsevierStylePara elsevierViewall">The conceptual model resulting from the interpretation of lithological logs&#44; geological mapping&#44; analysis of samples&#44; and isotopic data is shown in <a class="elsevierStyleCrossRef" href="#fig0025">Figure 5</a> &#40;a&#44; b&#41;&#46; The wells integrated in the numerical model are El Nabo and Mompani II&#44; 1100 and 250 m depth&#44; respectively &#40;<a class="elsevierStyleCrossRef" href="#fig0025">Figure 5</a>a&#41;&#46; For the implementation of the numerical model&#44; the volcanic sequences&#44; sedimentary materials&#44; and geological faults were considered as independent hydrogeological units&#59; each one with a hydraulic conductivity tensor&#46;</p><elsevierMultimedia ident="fig0025"></elsevierMultimedia><p id="par0090" class="elsevierStylePara elsevierViewall">Groundwater simulation was accomplished with the Visual-Modflow software &#40;Schlumberger Water Services&#44; Inc&#41; and the Modflow 2000 code &#40;<a class="elsevierStyleCrossRef" href="#bib0095">Harbaugh <span class="elsevierStyleItalic">et al</span>&#46;&#44; 2000</a>&#41;&#46; In order to represent the hydraulic behavior of the mapped geological structures the numerical model was built with 100&#44;000 cells&#44; distributed in a horizontal grid of 100 x 100 cells in the N-S and E-W direction&#44; each one represents a square of 2500 m<span class="elsevierStyleSup">2</span>&#46; The model box is a square of 25<span class="elsevierStyleHsp" style=""></span>km<span class="elsevierStyleSup">2</span> with a top surface with an average elevation of 1945 masl reproducing the actual topography and a flat surface at the bottom of the model &#40;800 masl&#41;&#46; The model box was divided in 10 horizontal layers of cells with variable thickness from the surface&#44; according to the topography&#44; to reach a bottom flat layer&#46; Cells with similar hydraulic behavior were grouped to represent the geological units included in the model considering stratigraphic sequences and major structures&#46; According to their physical characteristics&#44; spatial distribution&#44; and thickness&#44; four geological units were defined in the conceptual model and represented in the numerical model shown in <a class="elsevierStyleCrossRef" href="#fig0025">Figure 5</a>a&#58;<ul class="elsevierStyleList" id="lis0005"><li class="elsevierStyleListItem" id="lsti0005"><span class="elsevierStyleLabel">&#40;1&#41;</span><p id="par0095" class="elsevierStylePara elsevierViewall"><span class="elsevierStyleItalic">TomDa-Rhyolite</span>&#44; this is the upper layer of the model&#44; and groups the Cenozoic volcanic rocks including the volcanic and pyroclastic fractured sequences located above the piezometric level in the study area&#46; These units&#44; observed in the lithological records of the Ejido Jurica&#44; UNAM-Juriquilla&#44; and Ejido El Nabo wells&#44; were not considered in the conceptual model because either they are located above the regional groundwater level &#40;vadose zone&#41; or are absent nearby the El Nabo well area &#40;<a class="elsevierStyleCrossRef" href="#fig0015">Figure 3</a>&#41;&#46; These units include&#58; Late Miocene Andesites and Basalts &#40;Tmt AB&#41;&#44; Early Miocene Conglomerates and Sandstones Sandstones &#40;Tmp Ar-Cg&#41;&#44; and Quaternary Alluvium &#40;Qal&#59; <a class="elsevierStyleCrossRef" href="#fig0010">Figure 2</a>&#41;&#59;</p></li><li class="elsevierStyleListItem" id="lsti0010"><span class="elsevierStyleLabel">&#40;2&#41;</span><p id="par0100" class="elsevierStylePara elsevierViewall"><span class="elsevierStyleItalic">Te Cg Ar Conglomerate and sandstone</span>&#44; a sedimentary unit of conglomerate and sandstone found at the west of the El Nabo fault&#44; was represented by a dark layer shown in the model of <a class="elsevierStyleCrossRef" href="#fig0025">Figure 5</a>a&#59;</p></li><li class="elsevierStyleListItem" id="lsti0015"><span class="elsevierStyleLabel">&#40;3&#41;</span><p id="par0105" class="elsevierStylePara elsevierViewall"><span class="elsevierStyleItalic">Ki Cz Ar Limestone and shale</span>&#44; this unit includes Ki Cz and also the Ki Cz-Ar of sandstone and argillaceous limestone for modeling purposes&#44; and&#59;</p></li><li class="elsevierStyleListItem" id="lsti0020"><span class="elsevierStyleLabel">&#40;4&#41;</span><p id="par0110" class="elsevierStylePara elsevierViewall"><span class="elsevierStyleItalic">Ki MsV Metasedimentary and volcanic</span> &#40;<a class="elsevierStyleCrossRef" href="#fig0025">Figure 5</a>&#41;&#44; located at the bottom of the model&#44; corresponds in nature to the Cretaceous volcano-sedimentary unit &#40;Ki MsV&#41;&#46;</p></li></ul></p><p id="par0115" class="elsevierStylePara elsevierViewall">Faults are represented by groups of 50 m width vertically aligned cubic cells with specific hydraulic conductivity &#40;gray cells in <a class="elsevierStyleCrossRef" href="#fig0025">Figure 5</a>b&#41;&#46; The massive intrusive&#47;subvolcanic body of 30 Ma is represented by a group of cells with a dome shape in the located central part of the model &#40;black cells in <a class="elsevierStyleCrossRef" href="#fig0025">Figure 5</a>b&#41;&#46; Major geological structures were considered as four individual units in the model&#58; &#40;5&#41; the <span class="elsevierStyleItalic">El Nabo Fault</span>&#44; &#40;6&#41; the <span class="elsevierStyleItalic">Balvanera Fault</span>&#44; &#40;7&#41; the <span class="elsevierStyleItalic">South Fault</span> and&#44; &#40;8&#41; the intrusive body placed in the lower part of the model &#40;<a class="elsevierStyleCrossRef" href="#fig0010">Figure 2</a> and <a class="elsevierStyleCrossRef" href="#fig0025">Figure 5</a>b&#41;&#46;</p><p id="par0120" class="elsevierStylePara elsevierViewall">Faulting and the intrusive body are thought to greatly influence groundwater by channelizing or arresting flow&#44; or by interrupting the lateral continuity of the hydro stratigraphic units&#44; such as the conglomerate and sandstone unit &#40;<span class="elsevierStyleItalic">Te Cg-Ar&#41;</span>&#46; Truncation by faulting was suitable represented because each geological unit was composed by independent cells&#46;</p><p id="par0125" class="elsevierStylePara elsevierViewall">Once designed the geometry of the model&#44; the initial and border conditions were established considering regional variations of the piezometric surface interpolated from 2004 and 2006 measurements &#40;<a class="elsevierStyleCrossRef" href="#fig0030">Figure 6</a>&#41;&#46; Previous work by Carreon-Freyre <span class="elsevierStyleItalic">et al</span>&#46; &#40;2005&#41; indicated a regional gradient in the groundwater flow with a nearly north to south direction in the VQA&#46; Regional piezometric were interpolated by krigging of the data from 40 water extraction wells located in the VQA&#44; systematically monitored by the local agency of water &#40;<span class="elsevierStyleItalic">Comisi&#243;n Estatal de Aguas de Quer&#233;taro</span>&#44; CEAQ&#41; since 1996 &#40;<a class="elsevierStyleCrossRef" href="#bib0060">Carre&#243;n-Freyre <span class="elsevierStyleItalic">et al</span>&#46;&#44; 2005</a>&#41;&#46; The interpolations show a persistent regional gradient with an approximate North-South trending&#46; The wells located in the study area are El Nabo and Mompani II &#40;<a class="elsevierStyleCrossRef" href="#fig0005">Figure 1</a> and <a class="elsevierStyleCrossRef" href="#fig0025">Figure 5</a>a&#41;&#46; A steady state numerical simulation was performed to establish the initial and boundary conditions &#40;<a class="elsevierStyleCrossRef" href="#fig0035">Figure 7</a>b&#41;&#46; Boundary effects were minimized by leaving a minimum radius of 2000 m from each well to the model borders&#46;</p><elsevierMultimedia ident="fig0030"></elsevierMultimedia><elsevierMultimedia ident="fig0035"></elsevierMultimedia></span><span id="sec0040" class="elsevierStyleSection elsevierViewall"><span class="elsevierStyleSectionTitle" id="sect0060">Simulation of hydrodynamic conditions using the numerical model</span><p id="par0130" class="elsevierStylePara elsevierViewall">A numerical simulation was performed to assess the hydrodynamic behavior of the system by the use of partial results of a pumping test&#46; In the computational experiment the hydraulic properties were adjusted to correlate calculated and measured drawdown in the pumping and in one observation well&#46;</p><p id="par0135" class="elsevierStylePara elsevierViewall">A 260&#44;000 s &#40;72 hour&#41; pumping test performed in the <span class="elsevierStyleItalic">El Nabo</span> well was simulated using the implemented numerical model&#46; The location of the pumping test area is shown in <a class="elsevierStyleCrossRef" href="#fig0010">Figure 2</a>&#59; the data record &#40;water level and temperature&#41; was made by a CTD Diver &#40;<span class="elsevierStyleItalic">Schlumberger</span>&#41; datalogger&#46; This well is the single one in the region with about 1000 m in depth and therefore it was not possible to locate any observation well&#46; The <span class="elsevierStyleItalic">Mompani II</span> well is located 1000 m to the SW&#44; with a depth of 250 m&#46; Because of the presence of the El Nabo Fault between the two wells no change in the hydraulic head during the pumping test was recorded at the Mompani II and thus&#44; we assumed that groundwater flow in both wells were independent &#40;possibly local and regional systems&#41;&#46;</p><p id="par0140" class="elsevierStylePara elsevierViewall">To accomplish the simulation&#44; the hydraulic conductivities of the model and measured drawdown were closely fitted&#46; We used a segment of the pumping test to simulate transient conditions generated by the pumping test&#46; Continuous measurements of drawdown and temperature every 15 s are presented in <a class="elsevierStyleCrossRef" href="#fig0040">Figure 8</a> &#40;blue and green line respectively&#41;&#59; blue squares indicate samples from the continuous measurements used for simulation&#46; Because of the well depth only a part of the pumping test could be considered for adjustment&#46; The first 10&#44;800 s data &#40;3<span class="elsevierStyleHsp" style=""></span>hours&#41; were not considered because the pumping rate could not be maintained stable during this period and&#47;or a well storage effect should be considered &#40;<a class="elsevierStyleCrossRef" href="#bib0025">Allen and Michel&#44; 1999</a>&#41;&#44; the same issues were considered for the interpretation of the recovery data&#46;</p><elsevierMultimedia ident="fig0040"></elsevierMultimedia><p id="par0145" class="elsevierStylePara elsevierViewall">Simulated and observed drawdown in the segment between 10&#44;800 and 250&#44;000 s &#40;3 and 69<span class="elsevierStyleHsp" style=""></span>hours&#41; of the El Nabo test were fitted by adjusting the hydrogeological properties of the main units &#40;e&#46;g&#46;&#44; <a class="elsevierStyleCrossRef" href="#bib0025">Allen and Michel&#44; 1999</a>&#41;&#46; The estimation of these properties was achieved using the PEST &#40;<span class="elsevierStyleItalic">Parameter Estimation Software</span>&#41; algorithm &#40;<a class="elsevierStyleCrossRefs" href="#bib0095">Harbaugh <span class="elsevierStyleItalic">et al</span>&#46;&#44; 2000&#59; Doherty&#44; 2005</a>&#41; trough Visual-Modflow&#46; The finite differences Modflow algorithm solved groundwater flow by adjusting hydraulic conductivities in three perpendicular directions&#44; corresponding with the orientation of the grid to North-South&#44; East-West&#44; and vertically for this study case&#46; The initial value of hydraulic conductivity &#40;K&#41; was established at 1e-04 m&#47;s&#44; corresponding to reported values for massive and fractured rocks &#40;<a class="elsevierStyleCrossRefs" href="#bib0175">Walton&#44; 1970&#59; Ingebritsen <span class="elsevierStyleItalic">et al</span>&#46;&#44; 2006&#59; Zoback&#44; 2007</a>&#41;&#46;</p><p id="par0150" class="elsevierStylePara elsevierViewall">The best fit of the observed and calculated drawdowns was achieved by using the hydraulic conductivities &#40;K&#41; values presented in <a class="elsevierStyleCrossRef" href="#tbl0005">Table 1</a>&#44; along two directions&#58; horizontal E-W and vertically&#46; The North-South values of K remained fixed at the initial value&#46; The sensitivity values for the obtained hydraulic conductivities &#40;<a class="elsevierStyleCrossRef" href="#tbl0005">Table 1</a>&#41; of the model vary with the distance between the geological units&#47;faults and the pumping test zone&#46; High sensitivity means a high influence of the geological unit in drawdown&#46; The &#8220;limestone and shale&#8221; and &#8220;El Nabo fault&#8221; units show the highest sensitivity values &#40;<a class="elsevierStyleCrossRef" href="#tbl0005">Table 1</a>b&#41;&#46; The high values in the El Nabo Fault are due to the fact that the pumping well is located along this fault&#46; The lower obtained sensitivities correspond to the Balvanera and South faults&#44; meaning that the conductivity values of these units did not affect the calculated drawdown&#46;</p><elsevierMultimedia ident="tbl0005"></elsevierMultimedia><p id="par0155" class="elsevierStylePara elsevierViewall">The results of the pumping test were analyzed by the derivative of the data obtained by applying the fixed-end-point method &#40;orange line&#41; proposed by <a class="elsevierStyleCrossRef" href="#bib0040">Bourdet <span class="elsevierStyleItalic">et al</span>&#46; &#40;1989&#41;</a> &#40;orange line in <a class="elsevierStyleCrossRef" href="#fig0045">Figure 9</a>&#41;&#46; The use of the derivative allows analyzing in detail pressure changes in drawdown plots &#40;<a class="elsevierStyleCrossRefs" href="#bib0040">Bourdet <span class="elsevierStyleItalic">et al</span>&#46;&#44; 1989&#59; Spane and Wurstner&#44; 1993&#59; Tiab&#44; 2005</a>&#41;&#46; The derivative was smoothed by the use of the third point before and after the point of interest to reduce the sample rate noise&#46; The intersection between dyke&#47;fault and the El Nabo well was analyzed by the semi-analytical solution of <a class="elsevierStyleCrossRef" href="#bib0050">Butler and Liu &#40;1991&#41;</a>&#46;</p><elsevierMultimedia ident="fig0045"></elsevierMultimedia></span></span><span id="sec0045" class="elsevierStyleSection elsevierViewall"><span class="elsevierStyleSectionTitle" id="sect0065">Discussion of results</span><p id="par0160" class="elsevierStylePara elsevierViewall">The implemented model allowed representing the faulted system of the study area by the hydraulic conductivities assigned to specified geological units&#44; including faults&#46; In order to constrain the overall flow conditions regional measurements of hydraulic heads and geological observations were integrated&#46;</p></span><span id="sec0050" class="elsevierStyleSection elsevierViewall"><span class="elsevierStyleSectionTitle" id="sect0070">Geological constraints for the imple-mentation of the model</span><p id="par0165" class="elsevierStylePara elsevierViewall">The Cretaceous sequence drilled in El Nabo well&#44; composed by a limestone and shale unit &#40;Ki Cz-Ar&#41; that overlies a volcano-sedimentary sequence &#40;Ki Mvs&#41;&#44; was of the most importance to establish the geological constraints of the model&#46; The volcanic character of the basal volcano-sedimentary sequence was interpreted from the drilling samples and was confirmed by the isotopic ages&#46; The overlying limestone and shale unit is more than 800 m thick and is inter-bedded both with sandstone or clay&#46; In nearby superficial outcrops this unit is deformed by folding and contains calcite veins and fractures&#46; Rock samples obtained from depths below 850 m at the El Nabo well of the volcanic sequence contained a large quantity of biotite&#46; Differences in the size and aspect of biotite defined two size populations&#58; &#40;1&#41; between 425 and 710<span class="elsevierStyleHsp" style=""></span>&#956;m and &#40;2&#41; between 250 and 425<span class="elsevierStyleHsp" style=""></span>&#956;m&#46; The age spectra obtained for both sizes confirms that they correspond also to two different age populations&#46;</p><p id="par0170" class="elsevierStylePara elsevierViewall"><a class="elsevierStyleCrossRef" href="#fig0050">Figure 10</a>a shows the age spectra for the biotite population 1 indicating a Miocene age&#46; The weighted average for argon fractions was made by considering the consecutive fractions with greater radiogenic argon content and more than 10&#37; of the total <span class="elsevierStyleSup">39</span>Ar released&#46; A representative age of t &#61; 9&#46;5 &#177; 0&#46;1 Ma was selected from the weighted average of 13 fractions from different experiments in the age spectra diagram&#46; This biotite population is characterized by a ratio <span class="elsevierStyleSup">37</span>Ar<span class="elsevierStyleInf">Ca</span>&#47;<span class="elsevierStyleSup">39</span>Ar<span class="elsevierStyleInf">K</span> &#8804; 0&#46;5&#46; The correlation diagrams obtained for these samples &#40;<a class="elsevierStyleCrossRef" href="#fig0065">Figure S1</a>&#41; show two isochron lines also converging at an age around 9&#46;5 Ma&#46; The Andesitic lava flows &#40;<a class="elsevierStyleCrossRef" href="#bib0170">Valdez-Moreno <span class="elsevierStyleItalic">et al</span>&#46;&#44; 1998</a>&#41; cover most of the study area and have the same Miocene age obtained for biotite 1 population&#44; the coincidence in age of the lavas and the biotite sampled at depth strongly support the interpretation that the lavas were feed by dikes emplaced along faults&#46;</p><elsevierMultimedia ident="fig0050"></elsevierMultimedia><elsevierMultimedia ident="fig0065"></elsevierMultimedia><p id="par0175" class="elsevierStylePara elsevierViewall">The biotite population 2 &#40;<a class="elsevierStyleCrossRef" href="#fig0050">Figure 10</a>b&#41; yielded older ages of &#8764;100 Ma with age spectra diagram indicating a severe argon loss&#46; A reliable age cannot be obtained but the sequence is older than the minimum age obtained&#46; This behavior was reproduced in the 855 m depth Biotite in duplicate experiments&#46; A higher ratio of <span class="elsevierStyleSup">37</span>Ar<span class="elsevierStyleInf">Ca</span>&#47;<span class="elsevierStyleSup">39</span>Ar<span class="elsevierStyleInf">K</span> &#62; 0&#46;5 compared with biotite population 1 clearly indicates a different composition&#46; Because of the different isotopic behavior of these samples the biotite populations were plotted separately in the correlation diagram &#40;<a class="elsevierStyleCrossRef" href="#fig0065">Figure S1</a>&#41;&#44; in which the equivalent age is shown in the abscissa axis&#46; The correlation using the fractions from 3 to 6 of the 1st experiment and the fractions 3 to 8 of the duplicate experiment in the biotite 2 population also indicates an age of &#8764;100 Ma&#44; for the <span class="elsevierStyleItalic">Ki Ms V</span> rock sequence&#46;</p><p id="par0180" class="elsevierStylePara elsevierViewall">The documented difference in age of the biotite populations strongly supports our interpretation of the presence of a Miocene magmatic dyke intruding the Cretaceous sequence as shown in the schema of <a class="elsevierStyleCrossRef" href="#fig0050">Figure 10</a>c&#46; The presence of this fault has important implications in groundwater flow by the increasing of the hydraulic conductivity within the fault plane that was observed during drilling&#46;</p><p id="par0185" class="elsevierStylePara elsevierViewall">It has been documented that dikes open weakened fractures to conduct magma toward the surface and during their emplacement they can cause further brecciation of the fractured zone and the increase of hydraulic conductivity &#40;<a class="elsevierStyleCrossRefs" href="#bib0190">Zoback&#44; 2007&#59; Anderson&#44; 2006</a>&#41;&#46; Beside the presence of the dike&#44; further evidences supporting the circulation of hydrothermal fluids along the fault include the presence of hydrothermally altered calcite veins in exposures and in drilling samples of the El Nabo well &#40;<a class="elsevierStyleCrossRef" href="#fig0020">Figure 4</a>b&#41;&#46; Present day hydrothermal activity along the fault was confirmed by the increase from 25 to 63<span class="elsevierStyleHsp" style=""></span>&#176;C in the temperature records during pumping test &#40;green line&#44; <a class="elsevierStyleCrossRef" href="#fig0040">Figure 8</a>&#41;&#46;</p><p id="par0190" class="elsevierStylePara elsevierViewall">In order to consider the role of major faults in groundwater flow we integrated the <span class="elsevierStyleItalic">El Nabo&#44; Balvanera</span> and <span class="elsevierStyleItalic">South</span> faults in the model&#46; <span class="elsevierStyleItalic">El Nabo</span> and <span class="elsevierStyleItalic">Balvanera</span> are sub-parallel normal faults with NW-SE orientation &#40;panoramic view&#44; <a class="elsevierStyleCrossRef" href="#fig0020">Figure 4</a>a&#41;&#46; The <span class="elsevierStyleItalic">Balvanera Fault</span> has a throw of more than 200 m after the Miocene lavas emplacement&#46; The <span class="elsevierStyleItalic">El Nabo Fault</span> is located in a block delimited by the regional graben faults as shown in the cross-section of <a class="elsevierStyleCrossRef" href="#fig0010">Figure 2</a>b&#46; Several faults and fractures with similar orientation were observed around the zone of the El Nabo well&#46; The Cretaceous sedimentary sequences &#40;Ki Cz-Ar&#41; have a low hydraulic conductivity as shown by the low flow rate measured in the El Nabo well above 800 m depth &#40;6 l&#47;s&#41;&#46; Nevertheless&#44; flow rate increased to 47 l&#47;s when drilling intersected the El Nabo Fault and dike at 850 m depth&#46; We infer that the dike caused brecciation&#44; fracturing&#44; and hydrothermal alteration of the host rock and consequently an increase in the hydraulic conductivity&#46; Evidence of other dikes emplaced on the <span class="elsevierStyleItalic">El Nabo Fault</span> area affecting other stratigraphic levels was observed in exposed parallel fault planes&#46; For modeling purposes&#44; the <span class="elsevierStyleItalic">Balvanera Fault</span> features are assumed to be similar to the <span class="elsevierStyleItalic">El Nabo Fault</span> and was calibrated as an independent unit&#46;</p><p id="par0195" class="elsevierStylePara elsevierViewall">The exposed plane of the older age <span class="elsevierStyleItalic">South Fault</span> &#40;<a class="elsevierStyleCrossRef" href="#fig0020">Figure 4</a>b&#41; contains open fractures filled with sediments and calcite that precipitated by the circulation of hydrothermal fluids&#46; The <span class="elsevierStyleItalic">South Fault</span> favored the location of the Oligocene sub-volcanic body &#40;30 Ma&#44; TomDa&#41; of silicic composition &#40;Alanis-Alvarez <span class="elsevierStyleItalic">et al</span>&#46;&#44; 2002&#41; observed as a dome structure at the surface &#40;<a class="elsevierStyleCrossRef" href="#fig0010">Figure 2</a>&#41;&#46; Given its large spatial dimensions this structure was considered as an individual unit in the numerical model &#40;<a class="elsevierStyleCrossRef" href="#fig0025">Figure 5</a>b&#41;&#46;</p><p id="par0200" class="elsevierStylePara elsevierViewall">The geological evidences suggest that groundwater flow through faults is significant in the study area&#58; thus&#44; the model integrated an increase of conductivity along the El Nabo fault&#46;</p><span id="sec0055" class="elsevierStyleSection elsevierViewall"><span class="elsevierStyleSectionTitle" id="sect0075">Assessment of groundwater flow in faulted and high heterogeneous areas</span><p id="par0205" class="elsevierStylePara elsevierViewall">The numerical model was constrained by the results of the pumping test in El Nabo well&#46; The record of the test was sub-sampled and the obtained drawdown plot is presented in <a class="elsevierStyleCrossRef" href="#fig0045">Figure 9</a> &#40;blue squares&#41; along its derivative &#40;<a class="elsevierStyleCrossRef" href="#bib0040">Bourdet <span class="elsevierStyleItalic">et al</span>&#46;&#44; 1989</a>&#41;&#46; Considering an infinite linear strip embedded in a host material &#40;matrix&#41; with differing hydraulic properties <a class="elsevierStyleCrossRef" href="#bib0050">Butler and Liu &#40;1991&#41;</a> presented derivative curves with slopes varying from 1&#58;2 and 1&#58;4 for a case in which hydraulic conductivity is much larger than the host material &#40;matrix&#41;&#46; In their analysis these slopes corresponded to an initial linear flow within the strip that changed to bilinear flow when the flow from host material contributes to discharge&#46; The slopes obtained from the drawdown derivatives of the El Nabo well pumping test are shown in <a class="elsevierStyleCrossRef" href="#fig0045">Figure 9</a>&#59; a slope 1&#58;1 &#40;green dashed line&#41; for the 20&#44;000-70&#44;000 s &#40;5-19 hour&#41; interval and a slope 1&#58;2 &#40;black dashed line&#41; for the last part of the data&#46; The first slope &#40;1&#58;1&#41; corresponds to the primary flow drawdown that may be explained by two conditions&#58; &#40;a&#41; groundwater flow contributed by the limestone unit presenting high K values in the E-W direction and &#40;b&#41; is enclosed within the fault walls &#40;<a class="elsevierStyleCrossRef" href="#tbl0005">Table 1</a>&#41; or&#44; &#40;2&#41; non lineal flow inside the fault walls and constrained by the regional gradient&#44; which magnitude is large enough to oppose flow toward the well&#46;</p><p id="par0210" class="elsevierStylePara elsevierViewall">The 1&#58;2 slope may be related either to discharge of a lineal regional flow within the El Nabo Dyke&#47;Fault &#40;i&#46;e&#46; strip unit&#44; <a class="elsevierStyleCrossRef" href="#bib0050">Butler and Liu&#44; 1991</a>&#41; or a new source of regional flow&#46; Support evidence of the contribution of a hydrothermal deep groundwater flow system through the fault was provided by the increase in temperature while pumping &#40;25 to 63<span class="elsevierStyleHsp" style=""></span>&#176;C&#44; green line in <a class="elsevierStyleCrossRef" href="#fig0040">Figure 8</a>&#41;&#46; The initial increase of temperature in a time interval of 83<span class="elsevierStyleHsp" style=""></span>minutes suggests mixing of local fresh groundwater and hydrothermal flow&#44; with the presence of a preferential conduct through the magmatic dyke&#47;fault&#46;</p><p id="par0215" class="elsevierStylePara elsevierViewall">The calibrated hydraulic conductivities &#40;obtained using PEST software&#41; were used to simulate the pumping test in the numerical model &#40;red line&#59; <a class="elsevierStyleCrossRef" href="#fig0045">Figure 9</a>&#41;&#46; The wellbore storage may have caused the rapid initial drawdown of the El Nabo pumping test &#40;the well is more than 1000 m depth and 6 inches in diameter&#41;&#46; Time required for pumping stabilization &#40;10&#44;000 s&#44; or approximate 3<span class="elsevierStyleHsp" style=""></span>hours&#44; <a class="elsevierStyleCrossRef" href="#fig0040">Figure 8</a>&#41; was not considered in the model&#44; drawdown initiated from a relative zero value&#46; A pseudo-steady flow regime observed in the 10&#44;800-35&#44;000 s &#40;3-10 hour&#41; intervals &#40;<a class="elsevierStyleCrossRefs" href="#fig0040">Figures 8 and 9</a>&#41; was not completely adjusted by the calculated drawdown&#46; From 20 hour &#40;75&#44;000 s&#41; the numerical solution converges with the data and the model results are better fitted &#40;<a class="elsevierStyleCrossRef" href="#fig0045">Figure 9</a>&#41;&#46; In spite of the initial differences&#44; we consider that model curve captures the overall behavior of the pumping test suggesting a similarity of hydraulic behavior between model and test&#46;</p><p id="par0220" class="elsevierStylePara elsevierViewall">The model allows discussing the contribution of each geological element in the overall hydraulic behavior of the system&#46; Faults were considered as vertical planar structures with varying hydraulic conductivities according to the E-W or vertical groundwater flow direction&#46; The block model of <a class="elsevierStyleCrossRef" href="#fig0055">Figure 11</a> summarize the main groundwater flow results using the K values presented in <a class="elsevierStyleCrossRef" href="#tbl0005">Table 1</a>&#44; and shows the piezometric levels reached at the end of the pumping test simulation&#46; Colors indicate piezometric levels&#44; which range from 1925 in the north area to 1820 masl toward the south&#46; The higher gradient is obtained toward the east of the <span class="elsevierStyleItalic">El Nabo Fault</span>&#46; Arrows&#44; indicating flow directions&#44; show clearly that major structures &#40;El Nabo and South faults&#41; favor ground-water flow patterns along planes and also delimit compartments with differences in hydraulic behavior&#46; Field observations show that the hydraulic interaction of the fault with the El Nabo well starts at 850 m of depth and the presence of hydrothermal groundwater can be interpreted in terms of ongoing deep lateral flow facilitated by the heterogeneous permeability within the fault-zone probably due to the dyke intrusion &#40;L&#243;pez and Smith&#44; 1995&#59; <a class="elsevierStyleCrossRef" href="#bib0120">Lopez and Smith&#44; 1996</a>&#41;&#46; Horizontal flow predominates in the model while the flow patterns deviate along the geometric boundaries of the intrusive body&#44; interrupting the connection between El Nabo and South faults&#44; and following the regional N-S gradient&#46;</p><elsevierMultimedia ident="fig0055"></elsevierMultimedia><p id="par0225" class="elsevierStylePara elsevierViewall"><a class="elsevierStyleCrossRef" href="#fig0060">Figure 12</a> shows one horizontal and two vertical cross sections &#40;in the E-W and N-S directions&#41; of the model &#40;<a class="elsevierStyleCrossRef" href="#fig0055">Figure 11</a>&#41; that allow to analyze the influence of major structures &#40;El Nabo and South faults and the intrusive body&#41; in groundwater flow patterns&#46; In the horizontal cross section at 1400 masl of <a class="elsevierStyleCrossRef" href="#fig0060">Figure 12</a>a the northern segment of the <span class="elsevierStyleItalic">El Nabo Fault</span> act as a barrier by interrupting groundwater flow from north to south within the limestone and shale unit&#44; also shown in the vertical BB&#8217; profile &#40;<a class="elsevierStyleCrossRef" href="#fig0060">Figure 12</a>c&#41;&#46; The hydraulic head observed at 1900 m north from the El Nabo fault decreases to 1840 m in the southern block&#46; In the East-West cross section &#40;AA¿ in <a class="elsevierStyleCrossRef" href="#fig0060">Fig&#46; 12</a>b&#41; the color scale indicates a difference in piezometric elevations of 40 m at both sides of the fault&#46; In the model the regional gradient is interrupted by the <span class="elsevierStyleItalic">El Nabo Fault</span> that acts as a lateral flow barrier in the east-west direction&#44; the vertical equipotential lines indicate vertical flow within the fault&#46; The NE-SW orientation of the <span class="elsevierStyleItalic">South Fault</span> seems to favor channelized lateral flow towards the SW along the main fault plane&#46; The <span class="elsevierStyleItalic">South Fault</span> yielded the highest K value in the horizontal E-W direction &#40;<a class="elsevierStyleCrossRef" href="#tbl0005">Table 1</a>&#41;&#46; The model is not sensitive to this structure and this high value may be an artifact produced in the model by the equalization of drawdown and regional gradient&#46; However&#44; field data suggest also a higher K for the <span class="elsevierStyleItalic">South Fault</span> than for <span class="elsevierStyleItalic">El Nabo Fault</span>&#46;</p><elsevierMultimedia ident="fig0060"></elsevierMultimedia><p id="par0230" class="elsevierStylePara elsevierViewall">The 1&#58;2 slope of the derivative curve &#40;orange line&#44; <a class="elsevierStyleCrossRef" href="#fig0045">Figure 9</a>&#41; was interpreted as lineal flow within a high permeability channel with semi-permeable walls&#44; interpretation consistent with the presence of <span class="elsevierStyleItalic">El Nabo Fault-Dyke</span> &#40;<a class="elsevierStyleCrossRef" href="#fig0060">Figure 12</a> b&#44; c&#41;&#46; The semi-permeable walls condition in the model is represented by the equipotential lines at both sides of the faults&#46; Smooth gradients within the stratigraphic units increase when crossing the fault indicating a high hydraulic gradient but with flow continuity&#46; The limestone and shale unit also show hydraulic continuity with the fault to the west&#59; whereas in the eastern part have a higher elevation suggesting a greater water contribution from the fault &#40;<a class="elsevierStyleCrossRef" href="#fig0060">Figure 12</a>c&#41;&#46;</p><p id="par0235" class="elsevierStylePara elsevierViewall">Given that only one observation point could be considered&#44; the El Nabo well is the deepest well drilled in the Queretaro area with most of the well depths about 300 m&#44; analysis of representativeness of the obtained hydraulic conductivity was not possible to do&#46; Nevertheless&#44; the information yielded by the pumping test was useful for the calibration of the numerical model allowing simulating groundwater flow in the study area&#46;</p><p id="par0240" class="elsevierStylePara elsevierViewall">In summary&#44; this work attempts to highlight the importance of integrating an accurate geological setting to understand groundwater flow dynamics in models of faulted aquifers in particular in cases were piezometric levels decrease constantly&#46;</p><p id="par0245" class="elsevierStylePara elsevierViewall">Fault zones in the study area have more than 50 m in thickness and several fault planes were observed where high and low K segments may coexist&#46; Because the style and relative age of faulting determines variations in its hydraulic properties we consider that in faulted and deep aquifers &#40;more than 200 m depth&#41; the assessment of the overall flow dynamics is as&#44; or more&#44; important as obtaining specific hydraulic properties of the hydro-stratigraphic units&#46;</p><p id="par0250" class="elsevierStylePara elsevierViewall">Field evidences of present day hydrothermal water circulation within fault planes in the study area and the model results allow proposing that faults collect groundwater flow and explain the flow rate increase to 47 l&#47;s at <span class="elsevierStyleItalic">El Nabo Fault</span>&#46; Moreover&#44; the Cretaceous sedimentary sequences show an overall low hydraulic conductivity&#46; The magmatic dikes caused brecciation&#44; open fractures&#44; and hydrothermal alteration of the Cretaceous host rock and consequently enhancing its hydraulic conductivity &#40;e&#46;g&#46;&#44; <a class="elsevierStyleCrossRef" href="#bib0185">Zeidouni&#44; 2012&#59;</a><a class="elsevierStyleCrossRef" href="#tbl0005">Table 1</a>&#41;&#46; The South fault system favored the emplacement of the Oligocene intrusive body&#46; Field observations suggest that this fault also favored fluid-flow over time along open fractures partially filled with sediments and secondary fault planes in a zone more than 50 m in width&#46; We have observed precipitation of calcite in fractures &#40;<a class="elsevierStyleCrossRef" href="#fig0020">Figure 4</a>b&#41;&#44; and probably infiltration of superficial water&#46; Correspondingly&#44; after the simulation the higher K values obtained in the East-West direction of the South Fault allowed the hydraulic continuity of the system&#44; according with the regional gradient &#40;<a class="elsevierStyleCrossRef" href="#fig0055">Figure 11</a>&#41;&#46;</p></span></span><span id="sec0060" class="elsevierStyleSection elsevierViewall"><span class="elsevierStyleSectionTitle" id="sect0080">Conclusions</span><p id="par0255" class="elsevierStylePara elsevierViewall">Implementation of the geological structure in a numerical groundwater flow model&#44; along with the estimation of hydraulic conductivities by a pumping test&#44; produced relevant information about flow dynamics in the faulted aquifer system at the north of the Valley of Queretaro&#46; This methodology can be applied to study groundwater flow in other areas of the Queretaro Valley and elsewhere in central Mexico&#44; where water is extracted from depths greater than 200 m and high uncertainty exist about the influence of major structures in groundwater flow conditions&#46;</p><p id="par0260" class="elsevierStylePara elsevierViewall">The study area is a block delimited by faults where the older geological units in the Queretaro Valley are exposed&#46; Considering faults and truncation of hydrostratigraphic units in the sequence leaded to the reconnaissance of compartments in the aquifer system also reproduced in the model&#46; The identification of groundwater flow patterns through the aquifer when modeling drawdown was consistent with measured piezometric levels&#46; The derivative analysis of the pumping test suggests the presence of two different groundwater flow systems&#58; &#40;1&#41; a flow&#44; with slope 1&#58;1&#44; related to the limestone and shale sequence and&#44; &#40;2&#41; a regional linear flow&#44; with slope 1&#58;2&#44; transporting water through faults&#46;</p><p id="par0265" class="elsevierStylePara elsevierViewall">The numerical model along with the drawdown field observations allowed estimating relative values of hydraulic conductivity for eight geological units &#40;<a class="elsevierStyleCrossRef" href="#tbl0005">Table 1</a>&#41;&#46; As expected&#44; the model is sensitive to the fault located near the El Nabo well &#40;<a class="elsevierStyleCrossRef" href="#fig0005">Figure 1</a>&#41;&#46; Of importance to the study is the presence of volcano-sedimentary sequences below thick limestone units&#46; Isotopic dating confirmed the presence of a Miocene magmatic dyke emplaced along the El Nabo Fault&#46; Both the Balvanera and El Nabo fault systems favored the emplacement of Miocene volcanism through dykes along the faults planes &#40;<a class="elsevierStyleCrossRef" href="#fig0020">Figure 4</a>&#41;&#46; The presence of the fault-dyke was corroborated with stratigraphic correlation of a well log&#44; and the results confirm its relevance to modify groundwater flow&#46; The increased high flow rate in the El Nabo well pumping test &#40;47 l&#47;s&#41; is an evidence of the high influence of the dike that allows groundwater to flow through the fault&#46;</p><p id="par0270" class="elsevierStylePara elsevierViewall">The obtained values agree with the conceptual model but connectivity between major faults cannot be evaluated only by the analysis of pumping tests&#46; The numerical model allowed exploring the relevance of geological structures in the simulation of groundwater flow &#40;<a class="elsevierStyleCrossRefs" href="#fig0055">Figures 11 and 12</a>&#41;&#46; In particular&#44; the model results highlight how faults and intrusive bodies emplaced along faults modify groundwater flow by changing hydraulic properties and consequently flow direction&#46; The model allowed improving the interpretation of the pumping tests&#46; Numerical modeling is useful for a better understanding of complex aquifer systems that cannot be evaluated by analytical methods&#46;</p></span></span>"
    "textoCompletoSecciones" => array:1 [
      "secciones" => array:14 [
        0 => array:3 [
          "identificador" => "xres551697"
          "titulo" => "Resumen"
          "secciones" => array:1 [
            0 => array:1 [
              "identificador" => "abst0005"
            ]
          ]
        ]
        1 => array:2 [
          "identificador" => "xpalclavsec569110"
          "titulo" => "Palabras clave"
        ]
        2 => array:3 [
          "identificador" => "xres551696"
          "titulo" => "Abstract"
          "secciones" => array:1 [
            0 => array:1 [
              "identificador" => "abst0010"
            ]
          ]
        ]
        3 => array:2 [
          "identificador" => "xpalclavsec569111"
          "titulo" => "Keywords"
        ]
        4 => array:2 [
          "identificador" => "sec0005"
          "titulo" => "Introduction"
        ]
        5 => array:2 [
          "identificador" => "sec0010"
          "titulo" => "Faulted aquifer systems"
        ]
        6 => array:2 [
          "identificador" => "sec0015"
          "titulo" => "Geology of the Study Area &#40;Northwest of the Valley of Queretaro Aquifer&#44; VQA&#41;"
        ]
        7 => array:3 [
          "identificador" => "sec0020"
          "titulo" => "Methodology"
          "secciones" => array:4 [
            0 => array:2 [
              "identificador" => "sec0025"
              "titulo" => "Stratigraphic correlation"
            ]
            1 => array:2 [
              "identificador" => "sec0030"
              "titulo" => "Biotite age determination"
            ]
            2 => array:2 [
              "identificador" => "sec0035"
              "titulo" => "Integration of conceptual and numerical models"
            ]
            3 => array:2 [
              "identificador" => "sec0040"
              "titulo" => "Simulation of hydrodynamic conditions using the numerical model"
            ]
          ]
        ]
        8 => array:2 [
          "identificador" => "sec0045"
          "titulo" => "Discussion of results"
        ]
        9 => array:3 [
          "identificador" => "sec0050"
          "titulo" => "Geological constraints for the imple-mentation of the model"
          "secciones" => array:1 [
            0 => array:2 [
              "identificador" => "sec0055"
              "titulo" => "Assessment of groundwater flow in faulted and high heterogeneous areas"
            ]
          ]
        ]
        10 => array:2 [
          "identificador" => "sec0060"
          "titulo" => "Conclusions"
        ]
        11 => array:2 [
          "identificador" => "xack186041"
          "titulo" => "Acknowledgements"
        ]
        12 => array:1 [
          "titulo" => "<span class="elsevierStyleSectionTitle" id="sect0095">Further reading</span>"
        ]
        13 => array:1 [
          "titulo" => "References"
        ]
      ]
    ]
    "pdfFichero" => "main.pdf"
    "tienePdf" => true
    "fechaRecibido" => "2013-01-15"
    "fechaAceptado" => "2014-09-24"
    "PalabrasClave" => array:2 [
      "es" => array:1 [
        0 => array:4 [
          "clase" => "keyword"
          "titulo" => "Palabras clave"
          "identificador" => "xpalclavsec569110"
          "palabras" => array:6 [
            0 => "Flujo de aguas subterr&#225;neas"
            1 => "sistema acu&#237;fero fallado"
            2 => "modelaci&#243;n num&#233;rica"
            3 => "estructuras volc&#225;nicas"
            4 => "diques"
            5 => "M&#233;xico"
          ]
        ]
      ]
      "en" => array:1 [
        0 => array:4 [
          "clase" => "keyword"
          "titulo" => "Keywords"
          "identificador" => "xpalclavsec569111"
          "palabras" => array:6 [
            0 => "Groundwater flow"
            1 => "faulted aquifers"
            2 => "numerical modeling"
            3 => "volcanic structures"
            4 => "dykes"
            5 => "Mexicos"
          ]
        ]
      ]
    ]
    "tieneResumen" => true
    "resumen" => array:2 [
      "es" => array:2 [
        "titulo" => "Resumen"
        "resumen" => "<span id="abst0005" class="elsevierStyleSection elsevierViewall"><p id="spar0005" class="elsevierStyleSimplePara elsevierViewall">En este trabajo se integra la estratigraf&#237;a local y fallas como restricciones para un modelo num&#233;rico de flujo del agua subterr&#225;nea en la parte norte del acu&#237;fero fallado del Valle de Quer&#233;taro&#46; La informaci&#243;n geol&#243;gica de base fue establecida a partir de la reinterpretaci&#243;n de registros litol&#243;gicos de pozos&#44; dataci&#243;n isot&#243;pica 40Ar&#42;&#47;39Ar de rocas y cartograf&#237;a de campo&#46; En particular&#44; la dataci&#243;n isot&#243;pica permiti&#243; obtener el tiempo de emplazamiento de un dique magm&#225;tico a lo largo de la falla de El Nabo&#46; El aumento de la tasa de flujo y la temperatura en el pozo El Nabo sugieren que la falla y el dique transportan un flujo regional&#46; Un modelo num&#233;rico fue implementado para evaluar las contribuciones de fallas y estructuras volc&#225;nicas en el flujo general de agua subterr&#225;nea&#46; El modelo realizado con el software Visual-Modflow y c&#243;digo Modflow 2000 fue parcialmente restringido por los resultados de una prueba de bombeo en pozo &#8220;El Nabo&#8221;&#44; el pozo m&#225;s profundo &#40;Lat&#46; 20&#176; 42 &#8216;14 &#8220;N&#44; Long&#46; 100&#176; 28&#8217;45&#8221; O&#44; aproximadamente de 1000 m de profundidad&#41; ubicado en la zona de Quer&#233;taro&#46; La interpretaci&#243;n de los resultados del modelo y las observaciones de campo sugieren que las fallas reciben un flujo regional&#44; compartimentalizan el sistema&#44; y son responsables de un importante aumento del gasto&#44; de 6 a 47 l&#47;s&#44; cuando se la perforaci&#243;n del pozo El Nabo bajo de 850 m profundidad&#46; El estudio que aqu&#237; se presenta permite la evaluaci&#243;n de fallas normales y intrusivos que determinan el flujo del agua subterr&#225;nea y podr&#237;a ser relevante para el conocimiento de la din&#225;mica de flujo valles volc&#225;nicos cercanos&#46; El modelo num&#233;rico implementado basado en informaci&#243;n geol&#243;gica&#44; permiti&#243; una mejor comprensi&#243;n de los sistemas acu&#237;feros complejos&#44; lo que no se puede evaluar por m&#233;todos anal&#237;ticos&#46; Por otra parte puede ser una aplicaci&#243;n &#250;til para una mejor interpretaci&#243;n de pruebas de bombeo y para a&#241;adir datos complementarios de otras pruebas y&#47;o de an&#225;lisis de hidrogeoqu&#237;mica&#46;</p></span>"
      ]
      "en" => array:2 [
        "titulo" => "Abstract"
        "resumen" => "<span id="abst0010" class="elsevierStyleSection elsevierViewall"><p id="spar0010" class="elsevierStyleSimplePara elsevierViewall">This work integrates local stratigraphy and faults to constraint a numerical model of groundwater flow in the North area of the Valley of Queretaro aquifer&#46; Basic geological information was established from the reinterpretation of well logs&#44; <span class="elsevierStyleSup">40</span>Ar&#42;&#47;<span class="elsevierStyleSup">39</span>Ar dating of rocks&#44; and field mapping&#46; In particular&#44; dating allowed to obtain the timing of emplacement of a magmatic dyke inferred to be emplaced along the El Nabo fault&#46; The increase of flow rate and temperature in the El Nabo well suggest that the fault and dike transport a regional flow&#46; A numerical model was implemented to evaluate the contributions of faulting and volcanic structures in the overall groundwater flow&#46; The model&#44; accomplished with the <span class="elsevierStyleItalic">Visual-Modflow</span> software and Modflow 2000 code&#44; was partially constrained by the results of a pumping test in the El Nabo well&#44; the deepest well &#40;Lat&#46; 20&#176;42&#8217;14&#8221;N&#44; Long&#46; 100&#176;28&#8217;45&#8221;W&#44; ca&#46; 1000 m depth&#41; located in the Queretaro area&#46; The interpretation of the model results and field observations suggest that faults collect the regional flow&#44; compartmentalize the reservoir&#44; and are responsible of the dramatic increase of flow rate from 6 to 47 l&#47;s when drilling the El Nabo well below a depth of 850 m&#46; The study presented here allowed the assessment of normal faults and intrusives determining groundwater flow and might be relevant to the knowledge of flow dynamics in nearby volcanic valleys&#46; The analysis of the pumping test suggest the presence of two different groundwater flow systems&#58; &#40;1&#41; a flow related with the limestone and shale unit and &#40;2&#41; a regional linear flow through faults&#46; The implemented numerical model&#44; based on geological information&#44; allowed a better understanding of complex aquifer systems that cannot be evaluated by analytical methods&#46; Moreover it can be a useful application for a better interpretation of pumping tests and to add complementary data of other pumping tests and&#47;or hydrogeochemical analysis&#46;</p></span>"
      ]
    ]
    "multimedia" => array:14 [
      0 => array:7 [
        "identificador" => "fig0005"
        "etiqueta" => "Figure 1"
        "tipo" => "MULTIMEDIAFIGURA"
        "mostrarFloat" => true
        "mostrarDisplay" => false
        "figura" => array:1 [
          0 => array:4 [
            "imagen" => "gr1.jpeg"
            "Alto" => 1939
            "Ancho" => 1579
            "Tamanyo" => 560678
          ]
        ]
        "descripcion" => array:1 [
          "en" => "<p id="spar0015" class="elsevierStyleSimplePara elsevierViewall">&#40;a&#41; The Valley of Queretaro is located in the central part of Mexico&#59; &#40;b&#41; Geological map of the Valley of Queretaro&#46; Dashed red square line on the geological map indicates the area considered in modeling&#44; black lines are major faults&#46; Green polygon represents the Cretaceous rocks exposed in the zone&#46; Medium and light gray units represent Oligocene and Miocene volcanic rocks&#44; respectively&#46; &#40;Modified form Carreon-Freyre <span class="elsevierStyleItalic">et al</span>&#46;&#44; 2005&#41;&#46;</p>"
        ]
      ]
      1 => array:7 [
        "identificador" => "fig0010"
        "etiqueta" => "Figure 2"
        "tipo" => "MULTIMEDIAFIGURA"
        "mostrarFloat" => true
        "mostrarDisplay" => false
        "figura" => array:1 [
          0 => array:4 [
            "imagen" => "gr2.jpeg"
            "Alto" => 2554
            "Ancho" => 1572
            "Tamanyo" => 555033
          ]
        ]
        "descripcion" => array:1 [
          "en" => "<p id="spar0020" class="elsevierStyleSimplePara elsevierViewall">&#40;a&#41; Simplified geological map of the study area&#44; the main structures considered in the model are the normal faults El Nabo&#44; Balvanera and South&#46; Logs of water extracting wells were projected to the A-A¿ section shown in &#40;b&#41;&#46; The main stratigraphic units are presented in the cross section&#46;</p>"
        ]
      ]
      2 => array:7 [
        "identificador" => "fig0015"
        "etiqueta" => "Figure 3"
        "tipo" => "MULTIMEDIAFIGURA"
        "mostrarFloat" => true
        "mostrarDisplay" => false
        "figura" => array:1 [
          0 => array:4 [
            "imagen" => "gr3.jpeg"
            "Alto" => 2011
            "Ancho" => 1428
            "Tamanyo" => 260191
          ]
        ]
        "descripcion" => array:1 [
          "en" => "<p id="spar0025" class="elsevierStyleSimplePara elsevierViewall">Stratigraphic correlation of the lithological logs of extraction wells used for the interpretation of the main structures and construction of the conceptual model&#46; &#40;a&#41; Correlation of wells in the NW-SE direction including Mompani II&#44; UNAM-Juriquilla&#44; Ejido Jurica&#44; and &#40;CEA&#41; El Salitre&#59; &#40;b&#41; correlation in the N-S direction Ejido El Nabo&#44; Mompani II&#44; and El Nabo&#46; The initial static water level in the El Nabo well is also shown&#46;</p>"
        ]
      ]
      3 => array:7 [
        "identificador" => "fig0020"
        "etiqueta" => "Figure 4"
        "tipo" => "MULTIMEDIAFIGURA"
        "mostrarFloat" => true
        "mostrarDisplay" => false
        "figura" => array:1 [
          0 => array:4 [
            "imagen" => "gr4.jpeg"
            "Alto" => 1137
            "Ancho" => 1771
            "Tamanyo" => 357479
          ]
        ]
        "descripcion" => array:1 [
          "en" => "<p id="spar0030" class="elsevierStyleSimplePara elsevierViewall">&#40;a&#41; General view to the NW of the study area&#44; the escarp of the Balvanera fault can be observed&#44; the trace of the El Nabo fault is inferred by a dotted line because it delimits NE side of the uplifted bock &#40;horst&#41;&#46; As reference the La Joya volcano is shown toward the NW of the study area&#59; &#40;b&#41; close view in outcrop of the South Fault &#40;pen for scale&#41; showing the contact between the intrusive body at the left and evidence of hydrothermal activity&#59; &#40;c&#41; outcrop of the intrusive body&#59; &#40;d&#41; outcrop of a plane associated to the El Nabo fault zone&#44; the contact with the intrusive body and the structure of a dyke with brecciated material are shown&#46;</p>"
        ]
      ]
      4 => array:7 [
        "identificador" => "fig0025"
        "etiqueta" => "Figure 5"
        "tipo" => "MULTIMEDIAFIGURA"
        "mostrarFloat" => true
        "mostrarDisplay" => false
        "figura" => array:1 [
          0 => array:4 [
            "imagen" => "gr5.jpeg"
            "Alto" => 1528
            "Ancho" => 1657
            "Tamanyo" => 308203
          ]
        ]
        "descripcion" => array:1 [
          "en" => "<p id="spar0035" class="elsevierStyleSimplePara elsevierViewall">&#40;a&#41; Model box formed by an horizontal square grid of 25<span class="elsevierStyleHsp" style=""></span>km<span class="elsevierStyleSup">2</span> and variable thickness layers&#46; El Nabo and Mompani II wells were integrated in the model&#46; Cells were grouped to represent hydro-stratigraphic layers with specific hydraulic properties&#46; &#40;b&#41; Major geological structures in the model&#58; &#40;1&#41; El Nabo fault-dyke&#44; &#40;2&#41; Balvanera fault&#44; &#40;3&#41; South fault and&#44; &#40;4&#41; an intrusive body&#46; The topographic surface is considered in the model&#46;</p>"
        ]
      ]
      5 => array:7 [
        "identificador" => "fig0030"
        "etiqueta" => "Figure 6"
        "tipo" => "MULTIMEDIAFIGURA"
        "mostrarFloat" => true
        "mostrarDisplay" => false
        "figura" => array:1 [
          0 => array:4 [
            "imagen" => "gr6.jpeg"
            "Alto" => 942
            "Ancho" => 1378
            "Tamanyo" => 236148
          ]
        ]
        "descripcion" => array:1 [
          "en" => "<p id="spar0040" class="elsevierStyleSimplePara elsevierViewall">Piezometric levels measured in the Queretaro Valley for &#40;a&#41; 2004 and &#40;b&#41; 2006&#46; The 2006 piezometric levels were used as initial condition for the model&#46; Because the North-South hydraulic gradient&#44; shown in the dashed square &#40;modeled area&#41;&#44; is persistent during 2004 and 2006 we considered steady state conditions&#46;</p>"
        ]
      ]
      6 => array:7 [
        "identificador" => "fig0035"
        "etiqueta" => "Figure 7"
        "tipo" => "MULTIMEDIAFIGURA"
        "mostrarFloat" => true
        "mostrarDisplay" => false
        "figura" => array:1 [
          0 => array:4 [
            "imagen" => "gr7.jpeg"
            "Alto" => 1528
            "Ancho" => 1344
            "Tamanyo" => 257995
          ]
        ]
        "descripcion" => array:1 [
          "en" => "<p id="spar0045" class="elsevierStyleSimplePara elsevierViewall">Piezometric conditions shown in the model box with &#40;a&#41; border and &#40;b&#41; initial conditions after the steady state measured in 2006 &#40;see <a class="elsevierStyleCrossRef" href="#fig0030">figure 6</a>b&#41;&#46;</p>"
        ]
      ]
      7 => array:7 [
        "identificador" => "fig0040"
        "etiqueta" => "Figure 8"
        "tipo" => "MULTIMEDIAFIGURA"
        "mostrarFloat" => true
        "mostrarDisplay" => false
        "figura" => array:1 [
          0 => array:4 [
            "imagen" => "gr8.jpeg"
            "Alto" => 1198
            "Ancho" => 1787
            "Tamanyo" => 223018
          ]
        ]
        "descripcion" => array:1 [
          "en" => "<p id="spar0050" class="elsevierStyleSimplePara elsevierViewall">Drawdown measurements collected during the El Nabo pumping test&#46; Blue line indicates continuous datalogger measurements&#44; blue squares the sampled points for calibration&#46; Green line indicates increasing temperature in the extracted water&#46;</p>"
        ]
      ]
      8 => array:7 [
        "identificador" => "fig0045"
        "etiqueta" => "Figure 9"
        "tipo" => "MULTIMEDIAFIGURA"
        "mostrarFloat" => true
        "mostrarDisplay" => false
        "figura" => array:1 [
          0 => array:4 [
            "imagen" => "gr9.jpeg"
            "Alto" => 1160
            "Ancho" => 1808
            "Tamanyo" => 208188
          ]
        ]
        "descripcion" => array:1 [
          "en" => "<p id="spar0055" class="elsevierStyleSimplePara elsevierViewall">Measured drawdown during the pumping test at El Nabo &#40;blue squares&#41; and calculated curve using the implemented model&#46; The derivative of the data obtained by applying the fixed-end-point method &#40;<a class="elsevierStyleCrossRef" href="#bib0040">Bourdet et al&#46;&#44; 1989</a>&#41; is presented in a orange line&#46; Derivative drawdown was analyzed by the semi-analytical solution of <a class="elsevierStyleCrossRef" href="#bib0050">Butler and Liu &#40;1991&#41;</a> obtaining 1&#58;2 and 1&#58;1 slopes&#46;</p>"
        ]
      ]
      9 => array:7 [
        "identificador" => "fig0050"
        "etiqueta" => "Figure 10"
        "tipo" => "MULTIMEDIAFIGURA"
        "mostrarFloat" => true
        "mostrarDisplay" => false
        "figura" => array:1 [
          0 => array:4 [
            "imagen" => "gr10.jpeg"
            "Alto" => 1901
            "Ancho" => 1028
            "Tamanyo" => 319713
          ]
        ]
        "descripcion" => array:1 [
          "en" => "<p id="spar0060" class="elsevierStyleSimplePara elsevierViewall">&#40;a&#41; 40Ar&#47;39Ar&#42; age spectra obtained for the biotite population 1 &#40;small size&#41; collected at different depths between 846 and 979 m in El Nabo well&#46; The age obtained for this population is interpreted as evidence for a Miocene dyke intrusion&#46; &#40;b&#41; 40Ar&#47;39Ar&#42; age spectra obtained for the biotite population 2 collected at the same depth interval of population 1&#46; Although a reliable age cannot be assigned to the samples&#44; a minimum age of 100 Ma is interpreted for the volcano-sedimentary sequence hosting the Miocene intrusion&#44; &#40;c&#41; schematic section showing the intrusion of the volcanic dyke &#40;Tmt AB&#44; 10 Ma&#41; within the Cretaceous sedimentary and volcanic units &#40;host rocks&#58; Ki Cz and Ki Ms V&#41; that changes dramatically the hydraulic conductivity of the system by brecciation&#46;</p>"
        ]
      ]
      10 => array:7 [
        "identificador" => "fig0055"
        "etiqueta" => "Figure 11"
        "tipo" => "MULTIMEDIAFIGURA"
        "mostrarFloat" => true
        "mostrarDisplay" => false
        "figura" => array:1 [
          0 => array:4 [
            "imagen" => "gr11.jpeg"
            "Alto" => 1105
            "Ancho" => 1877
            "Tamanyo" => 328498
          ]
        ]
        "descripcion" => array:1 [
          "en" => "<p id="spar0065" class="elsevierStyleSimplePara elsevierViewall">Model box of the study area presenting drawdown elevations and groundwater flow vectors after simulation of El Nabo pumping test using calibrated of hydraulic conductivity values shown in <a class="elsevierStyleCrossRef" href="#tbl0005">Table 1</a>&#46; Faults act as flow channel interfering with the regional flow gradient and compartmentalizing the aquifer system&#46; Arrows show the direction and sense of flow&#46; The large intrusive dome modifies flow in the fault&#46; The El Nabo fault-dyke exerts a major influence in groundwater flow&#46;</p>"
        ]
      ]
      11 => array:7 [
        "identificador" => "fig0060"
        "etiqueta" => "Figure 12"
        "tipo" => "MULTIMEDIAFIGURA"
        "mostrarFloat" => true
        "mostrarDisplay" => false
        "figura" => array:1 [
          0 => array:4 [
            "imagen" => "gr12.jpeg"
            "Alto" => 2376
            "Ancho" => 1882
            "Tamanyo" => 535765
          ]
        ]
        "descripcion" => array:1 [
          "en" => "<p id="spar0070" class="elsevierStyleSimplePara elsevierViewall">Horizontal and vertical sections of the model box are presented in <a class="elsevierStyleCrossRef" href="#fig0055">Figure 11</a>&#46; Colors indicate variation of piezometric heads from 1940 to 1820 masl&#46; &#40;a&#41; Horizontal cross section at 1400 masl &#40;location shown in <a class="elsevierStyleCrossRef" href="#fig0060">Figure 12</a>b&#41;&#44; the black lines represent the analyzed faults walls and intrusive body&#46; &#40;b&#41; and &#40;c&#41; are two cross sections of the model box in the NS and EW directions respectively&#46;</p>"
        ]
      ]
      12 => array:7 [
        "identificador" => "fig0065"
        "etiqueta" => "Figure S1"
        "tipo" => "MULTIMEDIAFIGURA"
        "mostrarFloat" => true
        "mostrarDisplay" => false
        "figura" => array:1 [
          0 => array:4 [
            "imagen" => "gr13.jpeg"
            "Alto" => 1272
            "Ancho" => 1648
            "Tamanyo" => 178229
          ]
        ]
        "descripcion" => array:1 [
          "en" => "<p id="spar0075" class="elsevierStyleSimplePara elsevierViewall">Correlation diagrams obtained for samples for Biotite Population 1 and 2&#46; Isochron lines &#40;A and B&#41; yield an age around 9&#46;5 Ma for samples of Biotite Population 1&#59; whereas Biotite Population 2 &#40;C and D&#41; converge at an age around 100 Ma&#46;</p>"
        ]
      ]
      13 => array:7 [
        "identificador" => "tbl0005"
        "etiqueta" => "Table 1"
        "tipo" => "MULTIMEDIATABLA"
        "mostrarFloat" => true
        "mostrarDisplay" => false
        "tabla" => array:1 [
          "tablatextoimagen" => array:1 [
            0 => array:2 [
              "tabla" => array:1 [
                0 => """
                  <table border="0" frame="\n
                  \t\t\t\t\tvoid\n
                  \t\t\t\t" class=""><thead title="thead"><tr title="table-row"><th class="td" title="table-head  " align="left" valign="top" scope="col" style="border-bottom: 2px solid black"><span class="elsevierStyleBold">Hydraulic conductivity &#40;m&#47;s&#41;</span>&nbsp;\t\t\t\t\t\t\n
                  \t\t\t\t</th><th class="td" title="table-head  " align="left" valign="top" scope="col" style="border-bottom: 2px solid black"><span class="elsevierStyleBold">E-W</span>&nbsp;\t\t\t\t\t\t\n
                  \t\t\t\t</th><th class="td" title="table-head  " align="left" valign="top" scope="col" style="border-bottom: 2px solid black"><span class="elsevierStyleBold">Vertical</span>&nbsp;\t\t\t\t\t\t\n
                  \t\t\t\t</th></tr></thead><tbody title="tbody"><tr title="table-row"><td class="td" title="table-entry  " align="left" valign="top">Rhyolite&nbsp;\t\t\t\t\t\t\n
                  \t\t\t\t</td><td class="td" title="table-entry  " align="left" valign="top">2&#46;8E-05&nbsp;\t\t\t\t\t\t\n
                  \t\t\t\t</td><td class="td" title="table-entry  " align="left" valign="top">6&#46;1E-06&nbsp;\t\t\t\t\t\t\n
                  \t\t\t\t</td></tr><tr title="table-row"><td class="td" title="table-entry  " align="left" valign="top">Conglomerate and sandstone&nbsp;\t\t\t\t\t\t\n
                  \t\t\t\t</td><td class="td" title="table-entry  " align="left" valign="top">3&#46;4E-05&nbsp;\t\t\t\t\t\t\n
                  \t\t\t\t</td><td class="td" title="table-entry  " align="left" valign="top">4&#46;0E-05&nbsp;\t\t\t\t\t\t\n
                  \t\t\t\t</td></tr><tr title="table-row"><td class="td" title="table-entry  " align="left" valign="top">Limestone and shale&nbsp;\t\t\t\t\t\t\n
                  \t\t\t\t</td><td class="td" title="table-entry  " align="left" valign="top">5&#46;6E-02&nbsp;\t\t\t\t\t\t\n
                  \t\t\t\t</td><td class="td" title="table-entry  " align="left" valign="top">1&#46;3E-05&nbsp;\t\t\t\t\t\t\n
                  \t\t\t\t</td></tr><tr title="table-row"><td class="td" title="table-entry  " align="left" valign="top">Basement Limestone&nbsp;\t\t\t\t\t\t\n
                  \t\t\t\t</td><td class="td" title="table-entry  " align="left" valign="top">9&#46;9E-05&nbsp;\t\t\t\t\t\t\n
                  \t\t\t\t</td><td class="td" title="table-entry  " align="left" valign="top">8&#46;8E-06&nbsp;\t\t\t\t\t\t\n
                  \t\t\t\t</td></tr><tr title="table-row"><td class="td" title="table-entry  " align="left" valign="top">Nabo-Fault NW-SE&nbsp;\t\t\t\t\t\t\n
                  \t\t\t\t</td><td class="td" title="table-entry  " align="left" valign="top">4&#46;4E-07&nbsp;\t\t\t\t\t\t\n
                  \t\t\t\t</td><td class="td" title="table-entry  " align="left" valign="top">4&#46;9E-04&nbsp;\t\t\t\t\t\t\n
                  \t\t\t\t</td></tr><tr title="table-row"><td class="td" title="table-entry  " align="left" valign="top">South fault NE-SW&nbsp;\t\t\t\t\t\t\n
                  \t\t\t\t</td><td class="td" title="table-entry  " align="left" valign="top">1&#46;2E&#43;02&nbsp;\t\t\t\t\t\t\n
                  \t\t\t\t</td><td class="td" title="table-entry  " align="left" valign="top">6&#46;1E-06&nbsp;\t\t\t\t\t\t\n
                  \t\t\t\t</td></tr><tr title="table-row"><td class="td" title="table-entry  " align="left" valign="top">Balvanera-fault NW-SE&nbsp;\t\t\t\t\t\t\n
                  \t\t\t\t</td><td class="td" title="table-entry  " align="left" valign="top">5&#46;1E-04&nbsp;\t\t\t\t\t\t\n
                  \t\t\t\t</td><td class="td" title="table-entry  " align="left" valign="top">9&#46;7E-06&nbsp;\t\t\t\t\t\t\n
                  \t\t\t\t</td></tr><tr title="table-row"><td class="td" title="table-entry  " align="left" valign="top">Intrusive&nbsp;\t\t\t\t\t\t\n
                  \t\t\t\t</td><td class="td" title="table-entry  " align="left" valign="top">6&#46;15E-05&nbsp;\t\t\t\t\t\t\n
                  \t\t\t\t</td><td class="td" title="table-entry  " align="left" valign="top">1&#46;39E-05&nbsp;\t\t\t\t\t\t\n
                  \t\t\t\t</td></tr><tr title="table-row"><td class="td" title="table-entry  " align="left" valign="top">a&#41;&nbsp;\t\t\t\t\t\t\n
                  \t\t\t\t</td><td class="td" title="table-entry  " align="" valign="top">&nbsp;\t\t\t\t\t\t\n
                  \t\t\t\t</td><td class="td" title="table-entry  " align="" valign="top">&nbsp;\t\t\t\t\t\t\n
                  \t\t\t\t</td></tr><tr title="table-row"><td class="td" title="table-entry  " align="left" valign="top"><span class="elsevierStyleBold">Sensitivity</span>&nbsp;\t\t\t\t\t\t\n
                  \t\t\t\t</td><td class="td" title="table-entry  " align="left" valign="top"><span class="elsevierStyleBold">E-W</span>&nbsp;\t\t\t\t\t\t\n
                  \t\t\t\t</td><td class="td" title="table-entry  " align="left" valign="top"><span class="elsevierStyleBold">Vertical</span>&nbsp;\t\t\t\t\t\t\n
                  \t\t\t\t</td></tr><tr title="table-row"><td class="td" title="table-entry  " align="left" valign="top">Rhyolite&nbsp;\t\t\t\t\t\t\n
                  \t\t\t\t</td><td class="td" title="table-entry  " align="left" valign="top">4&#46;0E-03&nbsp;\t\t\t\t\t\t\n
                  \t\t\t\t</td><td class="td" title="table-entry  " align="left" valign="top">0&#46;0E&#43;00&nbsp;\t\t\t\t\t\t\n
                  \t\t\t\t</td></tr><tr title="table-row"><td class="td" title="table-entry  " align="left" valign="top">Conglomerate and sandstone&nbsp;\t\t\t\t\t\t\n
                  \t\t\t\t</td><td class="td" title="table-entry  " align="left" valign="top">4&#46;0E-03&nbsp;\t\t\t\t\t\t\n
                  \t\t\t\t</td><td class="td" title="table-entry  " align="left" valign="top">4&#46;0E-03&nbsp;\t\t\t\t\t\t\n
                  \t\t\t\t</td></tr><tr title="table-row"><td class="td" title="table-entry  " align="left" valign="top">Limestone and shale&nbsp;\t\t\t\t\t\t\n
                  \t\t\t\t</td><td class="td" title="table-entry  " align="left" valign="top">1&#46;2E-01&nbsp;\t\t\t\t\t\t\n
                  \t\t\t\t</td><td class="td" title="table-entry  " align="left" valign="top">1&#46;6E-02&nbsp;\t\t\t\t\t\t\n
                  \t\t\t\t</td></tr><tr title="table-row"><td class="td" title="table-entry  " align="left" valign="top">Basement Limestone&nbsp;\t\t\t\t\t\t\n
                  \t\t\t\t</td><td class="td" title="table-entry  " align="left" valign="top">8&#46;8E-03&nbsp;\t\t\t\t\t\t\n
                  \t\t\t\t</td><td class="td" title="table-entry  " align="left" valign="top">1&#46;3E-02&nbsp;\t\t\t\t\t\t\n
                  \t\t\t\t</td></tr><tr title="table-row"><td class="td" title="table-entry  " align="left" valign="top">Nabo-Fault NW-SE&nbsp;\t\t\t\t\t\t\n
                  \t\t\t\t</td><td class="td" title="table-entry  " align="left" valign="top">1&#46;0E&#43;00&nbsp;\t\t\t\t\t\t\n
                  \t\t\t\t</td><td class="td" title="table-entry  " align="left" valign="top">2&#46;9E-01&nbsp;\t\t\t\t\t\t\n
                  \t\t\t\t</td></tr><tr title="table-row"><td class="td" title="table-entry  " align="left" valign="top">South fault E-W&nbsp;\t\t\t\t\t\t\n
                  \t\t\t\t</td><td class="td" title="table-entry  " align="left" valign="top">1&#46;2E-01&nbsp;\t\t\t\t\t\t\n
                  \t\t\t\t</td><td class="td" title="table-entry  " align="left" valign="top">0&#46;0E&#43;00&nbsp;\t\t\t\t\t\t\n
                  \t\t\t\t</td></tr><tr title="table-row"><td class="td" title="table-entry  " align="left" valign="top">Balvanera-fault NW-SE&nbsp;\t\t\t\t\t\t\n
                  \t\t\t\t</td><td class="td" title="table-entry  " align="left" valign="top">4&#46;0E-03&nbsp;\t\t\t\t\t\t\n
                  \t\t\t\t</td><td class="td" title="table-entry  " align="left" valign="top">0&#46;0E&#43;00&nbsp;\t\t\t\t\t\t\n
                  \t\t\t\t</td></tr><tr title="table-row"><td class="td" title="table-entry  " align="left" valign="top">Intrusive&nbsp;\t\t\t\t\t\t\n
                  \t\t\t\t</td><td class="td" title="table-entry  " align="left" valign="top">0&#46;0E&#43;00&nbsp;\t\t\t\t\t\t\n
                  \t\t\t\t</td><td class="td" title="table-entry  " align="left" valign="top">0&#46;0E&#43;00&nbsp;\t\t\t\t\t\t\n
                  \t\t\t\t</td></tr><tr title="table-row"><td class="td" title="table-entry  " align="left" valign="top">b&#41;&nbsp;\t\t\t\t\t\t\n
                  \t\t\t\t</td><td class="td" title="table-entry  " align="" valign="top">&nbsp;\t\t\t\t\t\t\n
                  \t\t\t\t</td><td class="td" title="table-entry  " align="" valign="top">&nbsp;\t\t\t\t\t\t\n
                  \t\t\t\t</td></tr></tbody></table>
                  """
              ]
              "imagenFichero" => array:1 [
                0 => "xTab890448.png"
              ]
            ]
          ]
        ]
        "descripcion" => array:1 [
          "en" => "<p id="spar0080" class="elsevierStyleSimplePara elsevierViewall">&#40;a&#41; East-West and vertical hydraulic conductivities obtained from simulation&#46; &#40;b&#41; Table of sensitivity variability for PEST analysis&#44; higher sensitivity values of the model are related to the El Nabo fault&#46; See text for details&#46;</p>"
        ]
      ]
    ]
    "bibliografia" => array:2 [
      "titulo" => "References"
      "seccion" => array:1 [
        0 => array:2 [
          "identificador" => "bibs0005"
          "bibliografiaReferencia" => array:39 [
            0 => array:3 [
              "identificador" => "bib0005"
              "etiqueta" => "Aguirre-Di¿az and Lo¿pez-Marti¿nez&#44; 2001"
              "referencia" => array:1 [
                0 => array:2 [
                  "contribucion" => array:1 [
                    0 => array:2 [
                      "titulo" => "The Amazcala caldera&#44; Quer&#233;taro&#44; Mexico&#46; Geology and geochronology"
                      "autores" => array:1 [
                        0 => array:2 [
                          "etal" => false
                          "autores" => array:2 [
                            0 => "G&#46;J&#46; Aguirre-Di¿az"
                            1 => "M&#46; Lo¿pez-Marti¿nez"
                          ]
                        ]
                      ]
                    ]
                  ]
                  "host" => array:1 [
                    0 => array:1 [
                      "Revista" => array:5 [
                        "tituloSerie" => "Journal of Volcanology and Geothermal Research"
                        "fecha" => "2001"
                        "volumen" => "11"
                        "paginaInicial" => "203"
                        "paginaFinal" => "218"
                      ]
                    ]
                  ]
                ]
              ]
            ]
            1 => array:3 [
              "identificador" => "bib0010"
              "etiqueta" => "Aguirre-Di¿az et al&#46;&#44; 2005"
              "referencia" => array:1 [
                0 => array:2 [
                  "contribucion" => array:1 [
                    0 => array:2 [
                      "titulo" => "Seismogenic Basin and Range and intra-arc normal faulting in the central Mexican Volcanic Belt&#44; Quer&#233;taro M&#233;xico"
                      "autores" => array:1 [
                        0 => array:2 [
                          "etal" => false
                          "autores" => array:3 [
                            0 => "G&#46;J&#46; Aguirre-Di¿az"
                            1 => "J&#46; Nieto-Obrego¿n"
                            2 => "R&#46; Zu¿n¿iga"
                          ]
                        ]
                      ]
                    ]
                  ]
                  "host" => array:1 [
                    0 => array:1 [
                      "Revista" => array:5 [
                        "tituloSerie" => "Geological Journal"
                        "fecha" => "2005"
                        "volumen" => "40"
                        "paginaInicial" => "215"
                        "paginaFinal" => "243"
                      ]
                    ]
                  ]
                ]
              ]
            ]
            2 => array:3 [
              "identificador" => "bib0015"
              "etiqueta" => "Alaniz-&#193;lvarez et al&#46;&#44; 2001"
              "referencia" => array:1 [
                0 => array:1 [
                  "referenciaCompleta" => "Alaniz-&#193;lvarez S&#46;A&#46;&#44; Nieto-Samaniego A&#46;F&#46;&#44; Reyes-Zaragoza M&#46;A&#46;&#44; Orozco-Esquivel M&#46;T&#46;&#44; Ojeda-Garc&#237;a A&#46;C&#46;&#44; Vassallo F&#46;L&#46;&#44; 2001&#44; Estratigraf&#237;a y deformaci&#243;n extensional en la regi&#243;n San Miguel de Allende-Quer&#233;taro&#44; M&#233;xico &#40;Stratigraphy and extensional deformation in the San Miguel de Allende-Queretaro region&#44; Mexico&#41;&#46; <span class="elsevierStyleItalic">Revista Mexicana de Ciencias Geologicas</span>&#44; 18&#44; 129-148&#46;"
                ]
              ]
            ]
            3 => array:3 [
              "identificador" => "bib0020"
              "etiqueta" => "Alaniz-&#193;lvarez et al&#46;&#44; 2002"
              "referencia" => array:1 [
                0 => array:1 [
                  "referenciaCompleta" => "Alaniz-&#193;lvarez S&#46;A&#46;&#44; Nieto-Samaniego A&#46;F&#46;&#44; Orozco-Esquivel M&#46;T&#46;&#44; Vassallo F&#46;L&#46;&#44; Xu S&#46;&#44; 2002&#44; El sistema de fallas Taxco-San Miguel de Allende&#58; implicaciones en la deformacio¿n post-eoc&#233;nica del centro de M&#233;xico &#40;The fault system Taxco-San Miguel de Allende&#58; implications for post-Eocene deformation of central Mexico&#41; <span class="elsevierStyleItalic">Boletin de la Sociedad Geologica Mexicana</span>&#44; 55&#44; 12-29&#46;"
                ]
              ]
            ]
            4 => array:3 [
              "identificador" => "bib0025"
              "etiqueta" => "Allen and Michel&#44; 1999"
              "referencia" => array:1 [
                0 => array:2 [
                  "contribucion" => array:1 [
                    0 => array:2 [
                      "titulo" => "Characterizing a Faulted Aquifer by Field Testing and Numerical Simulation"
                      "autores" => array:1 [
                        0 => array:2 [
                          "etal" => false
                          "autores" => array:2 [
                            0 => "D&#46;M&#46; Allen"
                            1 => "F&#46;A&#46; Michel"
                          ]
                        ]
                      ]
                    ]
                  ]
                  "host" => array:1 [
                    0 => array:2 [
                      "doi" => "10.1111/j.1745-6584.1999.tb01164.x"
                      "Revista" => array:6 [
                        "tituloSerie" => "Ground Water"
                        "fecha" => "1999"
                        "volumen" => "37"
                        "paginaInicial" => "718"
                        "paginaFinal" => "728"
                        "link" => array:1 [
                          0 => array:2 [
                            "url" => "https://www.ncbi.nlm.nih.gov/pubmed/19125925"
                            "web" => "Medline"
                          ]
                        ]
                      ]
                    ]
                  ]
                ]
              ]
            ]
            5 => array:3 [
              "identificador" => "bib0030"
              "etiqueta" => "Anderson&#44; 2006"
              "referencia" => array:1 [
                0 => array:2 [
                  "contribucion" => array:1 [
                    0 => array:2 [
                      "titulo" => "Analytical solutions for flow to a well through a fault"
                      "autores" => array:1 [
                        0 => array:2 [
                          "etal" => false
                          "autores" => array:1 [
                            0 => "E&#46;I&#46; Anderson"
                          ]
                        ]
                      ]
                    ]
                  ]
                  "host" => array:1 [
                    0 => array:1 [
                      "Revista" => array:5 [
                        "tituloSerie" => "Advances in Water Resources"
                        "fecha" => "2006"
                        "volumen" => "29"
                        "paginaInicial" => "1790"
                        "paginaFinal" => "1803"
                      ]
                    ]
                  ]
                ]
              ]
            ]
            6 => array:3 [
              "identificador" => "bib0035"
              "etiqueta" => "Bense et al&#46;&#44; 2013"
              "referencia" => array:1 [
                0 => array:2 [
                  "contribucion" => array:1 [
                    0 => array:2 [
                      "titulo" => "Fault zone hydrogeology"
                      "autores" => array:1 [
                        0 => array:2 [
                          "etal" => false
                          "autores" => array:5 [
                            0 => "V&#46;F&#46; Bense"
                            1 => "T&#46; Gleeson"
                            2 => "S&#46;E&#46; Loveless"
                            3 => "O&#46; Bour"
                            4 => "J&#46; Scibek"
                          ]
                        ]
                      ]
                    ]
                  ]
                  "host" => array:1 [
                    0 => array:1 [
                      "Revista" => array:5 [
                        "tituloSerie" => "Earth-Science Reviews"
                        "fecha" => "2013"
                        "volumen" => "127"
                        "paginaInicial" => "171"
                        "paginaFinal" => "192"
                      ]
                    ]
                  ]
                ]
              ]
            ]
            7 => array:3 [
              "identificador" => "bib0040"
              "etiqueta" => "Bourdet et al&#46;&#44; 1989"
              "referencia" => array:1 [
                0 => array:2 [
                  "contribucion" => array:1 [
                    0 => array:2 [
                      "titulo" => "Use of pressure derivative in well-test interpretation"
                      "autores" => array:1 [
                        0 => array:2 [
                          "etal" => false
                          "autores" => array:3 [
                            0 => "D&#46; Bourdet"
                            1 => "J&#46;A&#46; Ayoub"
                            2 => "Y&#46;M&#46; Pirard"
                          ]
                        ]
                      ]
                    ]
                  ]
                  "host" => array:1 [
                    0 => array:1 [
                      "Revista" => array:4 [
                        "tituloSerie" => "SPE Formation Evaluation June"
                        "fecha" => "1989"
                        "paginaInicial" => "293"
                        "paginaFinal" => "302"
                      ]
                    ]
                  ]
                ]
              ]
            ]
            8 => array:3 [
              "identificador" => "bib0045"
              "etiqueta" => "Burbey&#44; 2008"
              "referencia" => array:1 [
                0 => array:2 [
                  "contribucion" => array:1 [
                    0 => array:2 [
                      "titulo" => "The Influence of Geologic Structures on Deformation due to Ground Water Withdrawal"
                      "autores" => array:1 [
                        0 => array:2 [
                          "etal" => false
                          "autores" => array:1 [
                            0 => "T&#46;J&#46; Burbey"
                          ]
                        ]
                      ]
                    ]
                  ]
                  "host" => array:1 [
                    0 => array:2 [
                      "doi" => "10.1111/j.1745-6584.2007.00395.x"
                      "Revista" => array:6 [
                        "tituloSerie" => "Ground Water"
                        "fecha" => "2008"
                        "volumen" => "46"
                        "paginaInicial" => "202"
                        "paginaFinal" => "211"
                        "link" => array:1 [
                          0 => array:2 [
                            "url" => "https://www.ncbi.nlm.nih.gov/pubmed/18194320"
                            "web" => "Medline"
                          ]
                        ]
                      ]
                    ]
                  ]
                ]
              ]
            ]
            9 => array:3 [
              "identificador" => "bib0050"
              "etiqueta" => "Butler and Liu&#44; 1991"
              "referencia" => array:1 [
                0 => array:2 [
                  "contribucion" => array:1 [
                    0 => array:2 [
                      "titulo" => "Pumping test in non-uniform aquifers - the linear strip case"
                      "autores" => array:1 [
                        0 => array:2 [
                          "etal" => false
                          "autores" => array:2 [
                            0 => "J&#46;J&#46; Butler"
                            1 => "W&#46;Z&#46; Liu"
                          ]
                        ]
                      ]
                    ]
                  ]
                  "host" => array:1 [
                    0 => array:1 [
                      "Revista" => array:5 [
                        "tituloSerie" => "Journal of Hydrology"
                        "fecha" => "1991"
                        "volumen" => "128"
                        "paginaInicial" => "69"
                        "paginaFinal" => "99"
                      ]
                    ]
                  ]
                ]
              ]
            ]
            10 => array:3 [
              "identificador" => "bib0055"
              "etiqueta" => "Caine et al&#46;&#44; 1996"
              "referencia" => array:1 [
                0 => array:2 [
                  "contribucion" => array:1 [
                    0 => array:2 [
                      "titulo" => "Fault zone architecture and permeability structure"
                      "autores" => array:1 [
                        0 => array:2 [
                          "etal" => false
                          "autores" => array:3 [
                            0 => "J&#46;S&#46; Caine"
                            1 => "J&#46;P&#46; Evans"
                            2 => "C&#46;B&#46; Forster"
                          ]
                        ]
                      ]
                    ]
                  ]
                  "host" => array:1 [
                    0 => array:1 [
                      "Revista" => array:5 [
                        "tituloSerie" => "Geology"
                        "fecha" => "1996"
                        "volumen" => "24"
                        "paginaInicial" => "1025"
                        "paginaFinal" => "1028"
                      ]
                    ]
                  ]
                ]
              ]
            ]
            11 => array:3 [
              "identificador" => "bib0060"
              "etiqueta" => "Carre&#243;n-Freyre et al&#46;&#44; 2005"
              "referencia" => array:1 [
                0 => array:2 [
                  "contribucion" => array:1 [
                    0 => array:2 [
                      "titulo" => "Influencia de la estratigraf&#237;a y estructura geol&#243;gica en el flujo de agua subterr&#225;nea del Valle de Quer&#233;taro &#40;Influence of the stratigraphy and geologic structures on groundwater flow Quer&#233;taro Valley&#41;"
                      "autores" => array:1 [
                        0 => array:2 [
                          "etal" => false
                          "autores" => array:4 [
                            0 => "D&#46; Carre&#243;n-Freyre"
                            1 => "M&#46; Cerca"
                            2 => "L&#46; Luna-Gonz&#225;lez"
                            3 => "F&#46;J&#46; G&#225;mez-Gonz&#225;lez"
                          ]
                        ]
                      ]
                    ]
                  ]
                  "host" => array:1 [
                    0 => array:1 [
                      "Revista" => array:5 [
                        "tituloSerie" => "Revista Mexicana de Ciencias Geologicas"
                        "fecha" => "2005"
                        "volumen" => "22"
                        "paginaInicial" => "1"
                        "paginaFinal" => "18"
                      ]
                    ]
                  ]
                ]
              ]
            ]
            12 => array:3 [
              "identificador" => "bib0065"
              "etiqueta" => "Cerca et al&#46;&#44; 2007"
              "referencia" => array:1 [
                0 => array:1 [
                  "referenciaCompleta" => "Cerca M&#46;&#44; Ferrari L&#46;&#44; L&#243;pez-Mart&#237;nez M&#46;&#44; Martiny B&#46;&#44; Iriondo A&#46;&#44; 2007&#44; Late Cretaceous shortening and early Tertiary shearing in the central Sierra Madre del Sur&#44; southern Mexico&#58; Insights into the evolution of the Caribbean&#8211;North American plate interaction&#58; <span class="elsevierStyleItalic">Tectonics</span>&#44; 26&#44; TC3007&#44; doi&#58;10&#46;1029&#47;2006TC001981&#46;"
                ]
              ]
            ]
            13 => array:3 [
              "identificador" => "bib0070"
              "etiqueta" => "Conagua&#44; 2003"
              "referencia" => array:1 [
                0 => array:2 [
                  "contribucion" => array:1 [
                    0 => array:2 [
                      "titulo" => "Plan de manejo del Acu&#237;fero Valle de Quer&#233;taro"
                      "autores" => array:1 [
                        0 => array:2 [
                          "etal" => false
                          "autores" => array:1 [
                            0 => "Conagua"
                          ]
                        ]
                      ]
                    ]
                  ]
                  "host" => array:1 [
                    0 => array:1 [
                      "Libro" => array:3 [
                        "fecha" => "2003"
                        "editorial" => "Comisi&#243;n Nacional del Agua"
                        "editorialLocalizacion" => "M&#233;xico"
                      ]
                    ]
                  ]
                ]
              ]
            ]
            14 => array:3 [
              "identificador" => "bib0075"
              "etiqueta" => "Delinom&#44; 2009"
              "referencia" => array:1 [
                0 => array:2 [
                  "contribucion" => array:1 [
                    0 => array:2 [
                      "titulo" => "Structural geology controls on groundwater flow&#58; Lembang Fault case study&#44; West Java"
                      "autores" => array:1 [
                        0 => array:2 [
                          "etal" => false
                          "autores" => array:1 [
                            0 => "R&#46;M&#46; Delinom"
                          ]
                        ]
                      ]
                    ]
                  ]
                  "host" => array:1 [
                    0 => array:1 [
                      "Revista" => array:5 [
                        "tituloSerie" => "Indonesia&#46; Hydrogeology Journal"
                        "fecha" => "2009"
                        "volumen" => "17"
                        "paginaInicial" => "1011"
                        "paginaFinal" => "1023"
                      ]
                    ]
                  ]
                ]
              ]
            ]
            15 => array:3 [
              "identificador" => "bib0080"
              "etiqueta" => "Doherty&#44; 2005"
              "referencia" => array:1 [
                0 => array:1 [
                  "referenciaCompleta" => "Doherty J&#46;&#44; 2005&#44; PEST Model-Independent Parameter Estimation User Manual&#58; 5th Edition Watermark Numerical Computing&#44; 336 pp&#46;"
                ]
              ]
            ]
            16 => array:3 [
              "identificador" => "bib0085"
              "etiqueta" => "Ferrari&#44; 2000"
              "referencia" => array:1 [
                0 => array:2 [
                  "contribucion" => array:1 [
                    0 => array:2 [
                      "titulo" => "Avances en el conocimiento de la Faja Volc&#225;nica Transmexicana durante la &#250;ltima d&#233;cada"
                      "autores" => array:1 [
                        0 => array:2 [
                          "etal" => false
                          "autores" => array:1 [
                            0 => "L&#46; Ferrari"
                          ]
                        ]
                      ]
                    ]
                  ]
                  "host" => array:1 [
                    0 => array:1 [
                      "Revista" => array:5 [
                        "tituloSerie" => "Boletin de la Sociedad Geologica Mexicana"
                        "fecha" => "2000"
                        "volumen" => "53"
                        "paginaInicial" => "84"
                        "paginaFinal" => "92"
                      ]
                    ]
                  ]
                ]
              ]
            ]
            17 => array:3 [
              "identificador" => "bib0090"
              "etiqueta" => "Hall&#44; 1981"
              "referencia" => array:1 [
                0 => array:2 [
                  "contribucion" => array:1 [
                    0 => array:2 [
                      "titulo" => "The application of K-Ar and 40Ar&#47;39Ar methods to the dating of recent volcanics and the Laschamp event&#46; PhD Thesis"
                      "autores" => array:1 [
                        0 => array:2 [
                          "etal" => false
                          "autores" => array:1 [
                            0 => "C&#46;M&#46; Hall"
                          ]
                        ]
                      ]
                    ]
                  ]
                  "host" => array:1 [
                    0 => array:1 [
                      "Libro" => array:3 [
                        "fecha" => "1981"
                        "editorial" => "University of Toronto"
                        "editorialLocalizacion" => "Canada"
                      ]
                    ]
                  ]
                ]
              ]
            ]
            18 => array:3 [
              "identificador" => "bib0095"
              "etiqueta" => "Harbaugh et al&#46;&#44; 2000"
              "referencia" => array:1 [
                0 => array:1 [
                  "referenciaCompleta" => "Harbaugh A&#46;W&#46;&#44; Banta E&#46;R&#46;&#44; Hill M&#46;C&#46;&#44; Mcdonald M&#46;G&#46;&#44; 2000&#44; MODFLOW-2000&#44; the U&#46;S&#46; Geological Survey modular ground-water model - User guide to modularization concepts and the Ground-Water Flow Process&#46; U&#46;S&#46; Geological Survey Open-File Report 00-92&#44; 121 pp&#46;"
                ]
              ]
            ]
            19 => array:3 [
              "identificador" => "bib0100"
              "etiqueta" => "Hsieh and Neuman&#44; 1985"
              "referencia" => array:1 [
                0 => array:2 [
                  "contribucion" => array:1 [
                    0 => array:2 [
                      "titulo" => "Field determi-nation of the three-dimensional hydraulic conductivity tensor of anisotropic media&#44; 1 Theory"
                      "autores" => array:1 [
                        0 => array:2 [
                          "etal" => false
                          "autores" => array:2 [
                            0 => "P&#46;A&#46; Hsieh"
                            1 => "S&#46;P&#46; Neuman"
                          ]
                        ]
                      ]
                    ]
                  ]
                  "host" => array:1 [
                    0 => array:1 [
                      "Revista" => array:5 [
                        "tituloSerie" => "Water Resources Research"
                        "fecha" => "1985"
                        "volumen" => "21"
                        "paginaInicial" => "1655"
                        "paginaFinal" => "1665"
                      ]
                    ]
                  ]
                ]
              ]
            ]
            20 => array:3 [
              "identificador" => "bib0105"
              "etiqueta" => "Ingebritsen et al&#46;&#44; 2006"
              "referencia" => array:1 [
                0 => array:1 [
                  "referenciaCompleta" => "Ingebritsen S&#46;E&#46;&#44; Sanford W&#46;E&#46;&#44; Neuzil C&#46;E&#46;&#44; 2006&#44; Groundwater in Geologic Processes&#46; 2nd Editio&#46;&#58; <span class="elsevierStyleItalic">Cambridge University Press</span>&#44; 536 pp&#46;"
                ]
              ]
            ]
            21 => array:3 [
              "identificador" => "bib0110"
              "etiqueta" => "Jourde et al&#46;&#44; 2007"
              "referencia" => array:1 [
                0 => array:2 [
                  "contribucion" => array:1 [
                    0 => array:2 [
                      "titulo" => "Relationship between the geometrical and structural properties of layered fractured rocks and their effective permeability tensor&#46; A simulation study"
                      "autores" => array:1 [
                        0 => array:2 [
                          "etal" => false
                          "autores" => array:5 [
                            0 => "H&#46; Jourde"
                            1 => "P&#46; Fenart"
                            2 => "M&#46; Vinches"
                            3 => "S&#46; Pistre"
                            4 => "B&#46; Vayssade"
                          ]
                        ]
                      ]
                    ]
                  ]
                  "host" => array:1 [
                    0 => array:1 [
                      "Revista" => array:5 [
                        "tituloSerie" => "Journal of Hydrology"
                        "fecha" => "2007"
                        "volumen" => "337"
                        "paginaInicial" => "117"
                        "paginaFinal" => "132"
                      ]
                    ]
                  ]
                ]
              ]
            ]
            22 => array:3 [
              "identificador" => "bib0115"
              "etiqueta" => "Kulkarni et al&#46;&#44; 2000"
              "referencia" => array:1 [
                0 => array:2 [
                  "contribucion" => array:1 [
                    0 => array:2 [
                      "titulo" => "Hydrogeological framework for the Deccan basalt groundwater systems&#44; west-central India"
                      "autores" => array:1 [
                        0 => array:2 [
                          "etal" => false
                          "autores" => array:5 [
                            0 => "H&#46; Kulkarni"
                            1 => "S&#46;B&#46; Deolankar"
                            2 => "A&#46; Lalwani"
                            3 => "B&#46; Joseph"
                            4 => "S&#46; Pawar"
                          ]
                        ]
                      ]
                    ]
                  ]
                  "host" => array:1 [
                    0 => array:1 [
                      "Revista" => array:5 [
                        "tituloSerie" => "Hydrogeology Journal"
                        "fecha" => "2000"
                        "volumen" => "8"
                        "paginaInicial" => "368"
                        "paginaFinal" => "378"
                      ]
                    ]
                  ]
                ]
              ]
            ]
            23 => array:3 [
              "identificador" => "bib0120"
              "etiqueta" => "Lopez and Smith&#44; 1996"
              "referencia" => array:1 [
                0 => array:2 [
                  "contribucion" => array:1 [
                    0 => array:2 [
                      "titulo" => "Fluid flow in fault zones&#58; influence of hydraulic anysotropy and heterogeneity on the fluid flow and heat transfer regime"
                      "autores" => array:1 [
                        0 => array:2 [
                          "etal" => false
                          "autores" => array:2 [
                            0 => "D&#46;A&#46; Lopez"
                            1 => "L&#46; Smith"
                          ]
                        ]
                      ]
                    ]
                  ]
                  "host" => array:1 [
                    0 => array:1 [
                      "Revista" => array:5 [
                        "tituloSerie" => "Water Resources Research"
                        "fecha" => "1996"
                        "volumen" => "32"
                        "paginaInicial" => "3227"
                        "paginaFinal" => "3235"
                      ]
                    ]
                  ]
                ]
              ]
            ]
            24 => array:3 [
              "identificador" => "bib0125"
              "etiqueta" => "Maclay and Small&#44; 1983"
              "referencia" => array:1 [
                0 => array:2 [
                  "contribucion" => array:1 [
                    0 => array:2 [
                      "titulo" => "Hydrostratigraphic subdivisions and fault barriers of the Edwards Aquifer&#44; south-central Texas U&#46;S&#46;A"
                      "autores" => array:1 [
                        0 => array:2 [
                          "etal" => false
                          "autores" => array:2 [
                            0 => "R&#46;W&#46; Maclay"
                            1 => "T&#46;A&#46; Small"
                          ]
                        ]
                      ]
                    ]
                  ]
                  "host" => array:1 [
                    0 => array:1 [
                      "Revista" => array:5 [
                        "tituloSerie" => "Journal of Hydrology"
                        "fecha" => "1983"
                        "volumen" => "61"
                        "paginaInicial" => "127"
                        "paginaFinal" => "146"
                      ]
                    ]
                  ]
                ]
              ]
            ]
            25 => array:3 [
              "identificador" => "bib0130"
              "etiqueta" => "Mayer et al&#46;&#44; 2007"
              "referencia" => array:1 [
                0 => array:2 [
                  "contribucion" => array:1 [
                    0 => array:2 [
                      "titulo" => "Estimation of fault-zone conductance by calibration of a regional groundwater flow model&#58; Desert Hot Springs California&#46;"
                      "autores" => array:1 [
                        0 => array:2 [
                          "etal" => false
                          "autores" => array:4 [
                            0 => "A&#46; Mayer"
                            1 => "W&#46; May"
                            2 => "C&#46; Lukkarila"
                            3 => "J&#46; Diehl"
                          ]
                        ]
                      ]
                    ]
                  ]
                  "host" => array:1 [
                    0 => array:1 [
                      "Revista" => array:5 [
                        "tituloSerie" => "Hydrogeology Journal"
                        "fecha" => "2007"
                        "volumen" => "15"
                        "paginaInicial" => "1093"
                        "paginaFinal" => "1106"
                      ]
                    ]
                  ]
                ]
              ]
            ]
            26 => array:3 [
              "identificador" => "bib0135"
              "etiqueta" => "Pacheco&#44; 2002"
              "referencia" => array:1 [
                0 => array:2 [
                  "contribucion" => array:1 [
                    0 => array:2 [
                      "titulo" => "Response to pumping of wells in sloping fault zone aquifers"
                      "autores" => array:1 [
                        0 => array:2 [
                          "etal" => false
                          "autores" => array:1 [
                            0 => "F&#46;A&#46;L&#46; Pacheco"
                          ]
                        ]
                      ]
                    ]
                  ]
                  "host" => array:1 [
                    0 => array:1 [
                      "Revista" => array:5 [
                        "tituloSerie" => "Journal of Hydrology"
                        "fecha" => "2002"
                        "volumen" => "259"
                        "paginaInicial" => "116"
                        "paginaFinal" => "135"
                      ]
                    ]
                  ]
                ]
              ]
            ]
            27 => array:3 [
              "identificador" => "bib0140"
              "etiqueta" => "Rafini and Larocque&#44; 2012"
              "referencia" => array:1 [
                0 => array:1 [
                  "referenciaCompleta" => "Rafini S&#46;&#44; Larocque M&#46;&#44; 2012&#44; Numerical modeling of the hydraulic signatures of horizontal and inclined faults&#46; Hydrogeology Journal&#46; doi&#58; 10&#46;1007&#47;s10040-011-0812-4&#46;"
                ]
              ]
            ]
            28 => array:3 [
              "identificador" => "bib0145"
              "etiqueta" => "Rathod and Rushton&#44; 1991"
              "referencia" => array:1 [
                0 => array:2 [
                  "contribucion" => array:1 [
                    0 => array:2 [
                      "titulo" => "Interpretation of Pumping from two-zone layered aquifers using a numerical model"
                      "autores" => array:1 [
                        0 => array:2 [
                          "etal" => false
                          "autores" => array:2 [
                            0 => "K&#46;S&#46; Rathod"
                            1 => "K&#46;R&#46; Rushton"
                          ]
                        ]
                      ]
                    ]
                  ]
                  "host" => array:1 [
                    0 => array:1 [
                      "Revista" => array:5 [
                        "tituloSerie" => "Ground Water"
                        "fecha" => "1991"
                        "volumen" => "29"
                        "paginaInicial" => "499"
                        "paginaFinal" => "509"
                      ]
                    ]
                  ]
                ]
              ]
            ]
            29 => array:3 [
              "identificador" => "bib0150"
              "etiqueta" => "Spane and Wurstner&#44; 1993"
              "referencia" => array:1 [
                0 => array:2 [
                  "contribucion" => array:1 [
                    0 => array:2 [
                      "titulo" => "DERIV&#58; A computer program for calculating pressure derivatives for use in hydraulic test analysis"
                      "autores" => array:1 [
                        0 => array:2 [
                          "etal" => false
                          "autores" => array:2 [
                            0 => "F&#46;A&#46; Spane"
                            1 => "S&#46;K&#46; Wurstner"
                          ]
                        ]
                      ]
                    ]
                  ]
                  "host" => array:1 [
                    0 => array:1 [
                      "Revista" => array:5 [
                        "tituloSerie" => "Ground Water"
                        "fecha" => "1993"
                        "volumen" => "31"
                        "paginaInicial" => "814"
                        "paginaFinal" => "822"
                      ]
                    ]
                  ]
                ]
              ]
            ]
            30 => array:3 [
              "identificador" => "bib0155"
              "etiqueta" => "Steiger and J&#228;ger&#44; 1977"
              "referencia" => array:1 [
                0 => array:2 [
                  "contribucion" => array:1 [
                    0 => array:2 [
                      "titulo" => "Sub commission on Geochronology&#58; Convention on the use of decay constants in Geo and Cosmochronology"
                      "autores" => array:1 [
                        0 => array:2 [
                          "etal" => false
                          "autores" => array:2 [
                            0 => "R&#46; Steiger"
                            1 => "E&#46; J&#228;ger"
                          ]
                        ]
                      ]
                    ]
                  ]
                  "host" => array:1 [
                    0 => array:1 [
                      "Revista" => array:5 [
                        "tituloSerie" => "Earth and Planetary Science Letters"
                        "fecha" => "1977"
                        "volumen" => "36"
                        "paginaInicial" => "359"
                        "paginaFinal" => "362"
                      ]
                    ]
                  ]
                ]
              ]
            ]
            31 => array:3 [
              "identificador" => "bib0160"
              "etiqueta" => "Tiab&#44; 2005"
              "referencia" => array:1 [
                0 => array:2 [
                  "contribucion" => array:1 [
                    0 => array:2 [
                      "titulo" => "Analysis of pressure derivative data of hydraulically fractured wells by the Tiab&#39;s Direct Synthesis technique"
                      "autores" => array:1 [
                        0 => array:2 [
                          "etal" => false
                          "autores" => array:1 [
                            0 => "D&#46; Tiab"
                          ]
                        ]
                      ]
                    ]
                  ]
                  "host" => array:1 [
                    0 => array:1 [
                      "Revista" => array:5 [
                        "tituloSerie" => "Journal of Petroleum Science and Engineering"
                        "fecha" => "2005"
                        "volumen" => "49"
                        "paginaInicial" => "1"
                        "paginaFinal" => "21"
                      ]
                    ]
                  ]
                ]
              ]
            ]
            32 => array:3 [
              "identificador" => "bib0165"
              "etiqueta" => "Trejo-Moedano&#44; 1989"
              "referencia" => array:1 [
                0 => array:1 [
                  "referenciaCompleta" => "Trejo-Moedano A&#46;&#44; 1989&#44; Estratigraf&#237;a y propiedades mec&#225;nicas del subsuelo del valle de la zona urbana de Quer&#233;taro &#40;Stratigraphy and mechanical properties of the subsoil in the urban zone of Queretaro&#41;&#46; Universidad Aut&#243;noma de Quer&#233;taro&#44; M&#233;xico&#44; 115 pp&#46;"
                ]
              ]
            ]
            33 => array:3 [
              "identificador" => "bib0170"
              "etiqueta" => "Valdez-Moreno et al&#46;&#44; 1998"
              "referencia" => array:1 [
                0 => array:2 [
                  "contribucion" => array:1 [
                    0 => array:2 [
                      "titulo" => "El volc&#225;n La Joya&#44; estados de Quer&#233;taro y Guanajuato&#59; un estratovolc&#225;n mioc&#233;nico del Cintur&#243;n Volc&#225;nico Mexicano &#40;The volcano La Joya&#44; states of Queretaro and Guanajuato&#44; a Miocene stratovolcano from Mexican Volcanic Belt&#41;"
                      "autores" => array:1 [
                        0 => array:2 [
                          "etal" => false
                          "autores" => array:3 [
                            0 => "G&#46; Valdez-Moreno"
                            1 => "G&#46;J&#46; Aguirre-D&#237;az"
                            2 => "M&#46; L&#243;pez-Mart&#237;nez"
                          ]
                        ]
                      ]
                    ]
                  ]
                  "host" => array:1 [
                    0 => array:1 [
                      "Revista" => array:5 [
                        "tituloSerie" => "Revista Mexicana de Ciencias Geologicas"
                        "fecha" => "1998"
                        "volumen" => "15"
                        "paginaInicial" => "181"
                        "paginaFinal" => "197"
                      ]
                    ]
                  ]
                ]
              ]
            ]
            34 => array:3 [
              "identificador" => "bib0175"
              "etiqueta" => "Walton&#44; 1970"
              "referencia" => array:1 [
                0 => array:1 [
                  "referenciaCompleta" => "Walton W&#46;C&#46;&#44; 1970&#44; Groundwater Resources Evaluation&#46; Mc&#46; Graw-Hill&#44; New York&#46; 185 pp&#46;"
                ]
              ]
            ]
            35 => array:3 [
              "identificador" => "bib0180"
              "etiqueta" => "Xu et al&#46;&#44; 2011"
              "referencia" => array:1 [
                0 => array:2 [
                  "contribucion" => array:1 [
                    0 => array:2 [
                      "titulo" => "Structural analysis of a relay ramp in the Queretaro Graben&#44; central Mexico&#58; Implications for relay ramp development"
                      "autores" => array:1 [
                        0 => array:2 [
                          "etal" => false
                          "autores" => array:4 [
                            0 => "S&#46; Xu"
                            1 => "A&#46;F&#46; Nieto-Samaniego"
                            2 => "S&#46;A&#46; Alaniz-Alvarez"
                            3 => "L&#46;M&#46; Cerca-Martinez"
                          ]
                        ]
                      ]
                    ]
                  ]
                  "host" => array:1 [
                    0 => array:1 [
                      "Revista" => array:5 [
                        "tituloSerie" => "Revista Mexicana de Ciencias Geologicas"
                        "fecha" => "2011"
                        "volumen" => "28"
                        "paginaInicial" => "275"
                        "paginaFinal" => "289"
                      ]
                    ]
                  ]
                ]
              ]
            ]
            36 => array:3 [
              "identificador" => "bib0185"
              "etiqueta" => "Zeidouni&#44; 2012"
              "referencia" => array:1 [
                0 => array:1 [
                  "referenciaCompleta" => "Zeidouni M&#46;&#44; 2012&#44; Analytical model of leakage through fault to overlying formations&#46; Water Resources Research 48&#44; W00N02&#44; doi&#58;10&#46;1029&#47;2012WR012582&#46;"
                ]
              ]
            ]
            37 => array:3 [
              "identificador" => "bib0190"
              "etiqueta" => "Zoback&#44; 2007"
              "referencia" => array:1 [
                0 => array:1 [
                  "referenciaCompleta" => "Zoback M&#46;D&#46;&#44; 2007&#44; Reservoir Geomechanics&#46; Cambridge Press&#44; 449 p&#46;"
                ]
              ]
            ]
            38 => array:3 [
              "identificador" => "bib0195"
              "etiqueta" => "Zu&#241;iga et al&#46;&#44; 2003"
              "referencia" => array:1 [
                0 => array:1 [
                  "referenciaCompleta" => "Zu&#241;iga F&#46;R&#44; Pacheco J&#46;F&#46;&#44; Guzm&#225;n-Speziale M&#46;&#44; Aguirre-Diaz G&#46;J&#46;&#44; Espindola V&#46;H&#46;&#44; Nava E&#46;&#44; 2003&#44; The Sanfandila earthquake sequence of 1998&#44; Queretaro&#44; M&#233;xico&#58; activation of an undocumented fault in the northern edge of central Trans-Mexican Volcanic Belt&#46; Tectonophysics&#44; 361&#44; 229-238&#46; doi&#58;10&#46;1016&#47;S0040-1951&#40;02&#41;00606-6&#46;"
                ]
              ]
            ]
          ]
        ]
      ]
    ]
    "lecturaRecomendada" => array:1 [
      0 => array:3 [
        "vista" => "all"
        "titulo" => "<span class="elsevierStyleSectionTitle" id="sect0095">Further reading</span>"
        "seccion" => array:1 [
          0 => array:2 [
            "vista" => "all"
            "bibliografiaReferencia" => array:2 [
              0 => array:3 [
                "identificador" => "bib0200"
                "etiqueta" => "Lopez and Smith&#44; 1995"
                "referencia" => array:1 [
                  0 => array:2 [
                    "contribucion" => array:1 [
                      0 => array:2 [
                        "titulo" => "Fluid flow in fault zones&#58; analysis of theinterplay of convective circulation and topographically driven groundwater flow"
                        "autores" => array:1 [
                          0 => array:2 [
                            "etal" => false
                            "autores" => array:2 [
                              0 => "D&#46;A&#46; Lopez"
                              1 => "L&#46; Smith"
                            ]
                          ]
                        ]
                      ]
                    ]
                    "host" => array:1 [
                      0 => array:1 [
                        "Revista" => array:5 [
                          "tituloSerie" => "Water Resources Research"
                          "fecha" => "1995"
                          "volumen" => "31"
                          "paginaInicial" => "1489"
                          "paginaFinal" => "1503"
                        ]
                      ]
                    ]
                  ]
                ]
              ]
              1 => array:3 [
                "identificador" => "bib0205"
                "etiqueta" => "York et al&#46;&#44; 2004"
                "referencia" => array:1 [
                  0 => array:2 [
                    "contribucion" => array:1 [
                      0 => array:2 [
                        "titulo" => "Unified equations for the slope&#44; intercept&#44; and standard errors of the best straight line"
                        "autores" => array:1 [
                          0 => array:2 [
                            "etal" => false
                            "autores" => array:4 [
                              0 => "D&#46; York"
                              1 => "N&#46; Evensen"
                              2 => "M&#46; Lopez-Martinez"
                              3 => "D&#46;J&#46; De Basabe"
                            ]
                          ]
                        ]
                      ]
                    ]
                    "host" => array:1 [
                      0 => array:1 [
                        "Revista" => array:5 [
                          "tituloSerie" => "American Journal of Physics"
                          "fecha" => "2004"
                          "volumen" => "73"
                          "paginaInicial" => "367"
                          "paginaFinal" => "375"
                        ]
                      ]
                    ]
                  ]
                ]
              ]
            ]
          ]
        ]
      ]
    ]
    "agradecimientos" => array:1 [
      0 => array:4 [
        "identificador" => "xack186041"
        "titulo" => "Acknowledgements"
        "texto" => "<p id="par0275" class="elsevierStylePara elsevierViewall">The project FOMIX CONACYT-Quer&#233;taro&#44; Ref&#46; QRO-2004-C01-39&#44; supported this research in the Valley of Queretaro&#46; Comision Estatal de Aguas de Queretaro &#40;CEAQ&#41;&#44; the local water management agency&#44; is gratefully acknowledged for funding the isotopic dating and for the technical help provided during pumping tests&#46; F&#46;J&#46; Gamez-Gonz&#225;lez &#40;Conagua&#41;&#44; I&#46; Barr&#243;n-Medellin &#40;CEAQ&#41;&#44; and J&#46; D&#237;az Escarcega &#40;CEAQ&#41; provided useful discussions for the groundwater flow interpretation&#46; G&#46; Aguirre-D&#237;az and an anonymous reviewer helped to clarify the manuscript&#46; R&#46; Carrizosa-Gonz&#225;lez &#40;UNAM&#41;&#44; I&#46; Ortiz-Villase&#241;or &#40;CEAQ&#41;&#44; and J&#46; Perez-Villarreal helped in the field campaign and data compilation&#46;</p>"
        "vista" => "all"
      ]
    ]
  ]
  "idiomaDefecto" => "en"
  "url" => "/00167169/0000005400000003/v2_201509120041/S0016716915000197/v2_201509120041/en/main.assets"
  "Apartado" => array:4 [
    "identificador" => "40021"
    "tipo" => "SECCION"
    "en" => array:2 [
      "titulo" => "Articles"
      "idiomaDefecto" => true
    ]
    "idiomaDefecto" => "en"
  ]
  "PDF" => "https://static.elsevier.es/multimedia/00167169/0000005400000003/v2_201509120041/S0016716915000197/v2_201509120041/en/main.pdf?idApp=UINPBA00004N&text.app=https://www.elsevier.es/"
  "EPUB" => "https://multimedia.elsevier.es/PublicationsMultimediaV1/item/epub/S0016716915000197?idApp=UINPBA00004N"
]

Article information

ISSN: 00167169
Original language: English

DOI:10.1016/j.gi.2015.04.016

The statistics are updated each day

Year/Month	Html	Pdf	Total

2024 November	5	0	5
2024 October	34	2	36
2024 September	37	6	43
2024 August	40	4	44
2024 July	34	11	45
2024 June	34	4	38
2024 May	22	7	29
2024 April	59	6	65
2024 March	60	10	70
2024 February	42	14	56
2024 January	77	23	100
2023 December	47	7	54
2023 November	55	11	66
2023 October	104	12	116
2023 September	78	5	83
2023 August	70	9	79
2023 July	71	6	77
2023 June	65	21	86
2023 May	84	12	96
2023 April	75	7	82
2023 March	69	11	80
2023 February	64	8	72
2023 January	59	21	80
2022 December	62	9	71
2022 November	69	20	89
2022 October	82	19	101
2022 September	58	26	84
2022 August	63	15	78
2022 July	64	19	83
2022 June	41	9	50
2022 May	48	17	65
2022 April	31	13	44
2022 March	53	19	72
2022 February	59	8	67
2022 January	111	8	119
2021 December	58	14	72
2021 November	52	26	78
2021 October	55	20	75
2021 September	57	27	84
2021 August	84	21	105
2021 July	51	14	65
2021 June	41	13	54
2021 May	67	13	80
2021 April	135	44	179
2021 March	83	30	113
2021 February	58	9	67
2021 January	44	16	60
2020 December	81	25	106
2020 November	55	17	72
2020 October	31	6	37
2020 September	34	17	51
2020 August	74	18	92
2020 July	54	15	69
2020 June	51	19	70
2020 May	48	18	66
2020 April	52	19	71
2020 March	55	25	80
2020 February	45	14	59
2020 January	56	7	63
2019 December	94	19	113
2019 November	62	21	83
2019 October	46	9	55
2019 September	43	16	59
2019 August	22	6	28
2019 July	31	8	39
2019 June	38	22	60
2019 May	106	58	164
2019 April	56	13	69
2019 March	17	3	20
2019 February	28	13	41
2019 January	18	6	24
2018 December	15	10	25
2018 November	30	17	47
2018 October	60	10	70
2018 September	58	21	79
2018 August	38	4	42
2018 July	28	8	36
2018 June	30	6	36
2018 May	43	18	61
2018 April	43	5	48
2018 March	25	3	28
2018 February	44	7	51
2018 January	39	8	47
2017 December	23	5	28
2017 November	21	10	31
2017 October	33	5	38
2017 September	33	9	42
2017 August	19	13	32
2017 July	21	6	27
2017 June	34	36	70
2017 May	32	14	46
2017 April	21	42	63
2017 March	21	42	63
2017 February	24	2	26
2017 January	13	2	15
2016 December	31	10	41
2016 November	34	10	44
2016 October	31	6	37
2016 September	25	12	37
2016 August	22	2	24
2016 July	22	4	26
2016 June	38	22	60
2016 May	15	20	35
2016 April	19	15	34
2016 March	31	21	52
2016 February	27	21	48
2016 January	20	23	43
2015 December	16	15	31
2015 November	14	15	29
2015 October	23	19	42
2015 September	10	8	18

Show all

Indexed in:

Follow us:

Statistics

Indexed in:

Follow us:

Statistics

Subscribe to our newsletter