<?xml version="1.0" encoding="utf-8"?>
<!DOCTYPE article
  PUBLIC "-//NLM//DTD JATS (Z39.96) Journal Publishing DTD v1.1 20151215//EN" "https://jats.nlm.nih.gov/publishing/1.1/JATS-journalpublishing1.dtd">
<article article-type="research-article" dtd-version="1.1" specific-use="sps-1.8" xml:lang="es" xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink">
	<front>
		<journal-meta>
			<journal-id journal-id-type="publisher-id">rge</journal-id>
			<journal-title-group>
				<journal-title>Revista Geográfica de América Central</journal-title>
				<abbrev-journal-title abbrev-type="publisher">Rev. Geog. Amér. Central</abbrev-journal-title>
			</journal-title-group>
			<issn pub-type="ppub">1011-484X</issn>
			<issn pub-type="epub">2215-2563</issn>
			<publisher>
				<publisher-name>Universidad Nacional, Costa Rica</publisher-name>
			</publisher>
		</journal-meta>
		<article-meta>
			<article-id pub-id-type="doi">10.15359/rgac.63-2.5</article-id>
			<article-categories>
				<subj-group subj-group-type="heading">
					<subject>Artículo</subject>
				</subj-group>
			</article-categories>
			<title-group>
				<article-title>Vulnerabilidad hidrogeológica en la zona costera de Jacó, Pacífico Central, Costa Rica</article-title>
			</title-group>
			<contrib-group>
				<contrib contrib-type="author">
					<name>
						<surname>Gómez-Cruz</surname>
						<given-names>Alicia</given-names>
					</name>
					<xref ref-type="aff" rid="aff1"><sup>1</sup></xref>
				</contrib>
				<contrib contrib-type="author">
					<name>
						<surname>Núñez-Solís</surname>
						<given-names>Christian</given-names>
					</name>
					<xref ref-type="aff" rid="aff2"><sup>2</sup></xref>
				</contrib>
				<contrib contrib-type="author">
					<name>
						<surname>Madrigal-Solís</surname>
						<given-names>Helga</given-names>
					</name>
					<xref ref-type="aff" rid="aff3"><sup>3</sup></xref>
				</contrib>
				<contrib contrib-type="author">
					<name>
						<surname>Calderón-Sánchez</surname>
						<given-names>Hazel</given-names>
					</name>
					<xref ref-type="aff" rid="aff4"><sup>4</sup></xref>
				</contrib>
				<contrib contrib-type="author">
					<name>
						<surname>Fonseca-Sánchez</surname>
						<given-names>Alicia</given-names>
					</name>
					<xref ref-type="aff" rid="aff5"><sup>5</sup></xref>
				</contrib>
				<contrib contrib-type="author">
					<name>
						<surname>Jiménez-Gavilán</surname>
						<given-names>Pablo</given-names>
					</name>
					<xref ref-type="aff" rid="aff6"><sup>6</sup></xref>
				</contrib>
			</contrib-group>
			<aff id="aff1">
				<label>1</label>
				<institution content-type="original"> Laboratorio de Hidrología Ambiental. Escuela de Ciencias Biológicas. Universidad Nacional de Costa Rica. agomezcruz@gmail.com.</institution>
				<institution content-type="normalized">Universidad Nacional</institution>
				<institution content-type="orgdiv2">Laboratorio de Hidrología Ambiental</institution>
				<institution content-type="orgdiv1">Escuela de Ciencias Biológicas</institution>
				<institution content-type="orgname">Universidad Nacional de Costa Rica</institution>
				<country country="CR">Costa Rica</country>
				<email>agomezcruz@gmail.com</email>
			</aff>
			<aff id="aff2">
				<label>2</label>
				<institution content-type="original"> Laboratorio de Hidrología Ambiental. Escuela de Ciencias Biológicas. Universidad Nacional de Costa Rica. christian.nunez.solis@una.cr.</institution>
				<institution content-type="normalized">Universidad Nacional</institution>
				<institution content-type="orgdiv2">Laboratorio de Hidrología Ambiental</institution>
				<institution content-type="orgdiv1">Escuela de Ciencias Biológicas</institution>
				<institution content-type="orgname">Universidad Nacional de Costa Rica</institution>
				<country country="CR">Costa Rica</country>
				<email>christian.nunez.solis@una.cr</email>
			</aff>
			<aff id="aff3">
				<label>3</label>
				<institution content-type="original"> Laboratorio de Hidrología Ambiental. Escuela de Ciencias Biológicas. Universidad Nacional de Costa Rica. helga.madrigal.solis@una.cr.</institution>
				<institution content-type="normalized">Universidad Nacional</institution>
				<institution content-type="orgdiv2">Laboratorio de Hidrología Ambiental</institution>
				<institution content-type="orgdiv2">Escuela de Ciencias Biológicas</institution>
				<institution content-type="orgname">Universidad Nacional de Costa Rica</institution>
				<country country="CR">Costa Rica</country>
				<email>helga.madrigal.solis@una.cr</email>
			</aff>
			<aff id="aff4">
				<label>4</label>
				<institution content-type="original"> Laboratorio de Hidrología Ambiental. Escuela de Ciencias Biológicas. Universidad Nacional de Costa Rica. hazelcalderon@gmail.com.</institution>
				<institution content-type="normalized">Universidad Nacional</institution>
				<institution content-type="orgdiv2">Laboratorio de Hidrología Ambiental</institution>
				<institution content-type="orgdiv1">Escuela de Ciencias Biológicas</institution>
				<institution content-type="orgname">Universidad Nacional de Costa Rica</institution>
				<country country="CR">Costa Rica</country>
				<email>hazelcalderon@gmail.com</email>
			</aff>
			<aff id="aff5">
				<label>5</label>
				<institution content-type="original"> Laboratorio de Hidrología Ambiental. Escuela de Ciencias Biológicas. Universidad Nacional de Costa Rica. alicia.fonseca.sanchez@una.cr</institution>
				<institution content-type="normalized">Universidad Nacional</institution>
				<institution content-type="orgdiv2">Laboratorio de Hidrología Ambiental</institution>
				<institution content-type="orgdiv1">Escuela de Ciencias Biológicas</institution>
				<institution content-type="orgname">Universidad Nacional de Costa Rica</institution>
				<country country="CR">Costa Rica</country>
				<email>alicia.fonseca.sanchez@una.cr</email>
			</aff>
			<aff id="aff6">
				<label>6</label>
				<institution content-type="original"> Departamento de Ecología y Geología. Facultad de Ciencias. Universidad de Málaga. pgavilan@uma.es.</institution>
				<institution content-type="normalized">Universidad de Málaga</institution>
				<country country="ES">Spain</country>
			</aff>
			<pub-date pub-type="collection">
				<season>Jul-Dec</season>
				<year>2019</year>
			</pub-date>
			<issue>63</issue>
			<fpage>123</fpage>
			<lpage>145</lpage>
			<history>
				<date date-type="received">
					<day>10</day>
					<month>09</month>
					<year>2017</year>
				</date>
				<date date-type="accepted">
					<day>08</day>
					<month>09</month>
					<year>2018</year>
				</date>
			</history>
			<permissions>
				<license license-type="open-access" xlink:href="https://creativecommons.org/licenses/by-nc-sa/4.0/" xml:lang="es">
					<license-p>Este es un artículo publicado en acceso abierto bajo una licencia Creative Commons</license-p>
				</license>
			</permissions>
			<abstract>
				<title>Resumen</title>
				<p>En las últimas dos décadas, el Pacífico Central de Costa Rica experimentó un fuerte aumento en el desarrollo urbanístico y turístico, lo que representa una mayor amenaza a la calidad del agua subterránea. El objetivo de este estudio fue el generar una propuesta del Mapa de Vulnerabilidad Hidrogeológica de los acuíferos de Jacó, utilizando la metodología GOD. Para esto, se elaboraron perfiles geológicos y se realizó comprobación de campo. Se evidencia que el área de estudio está conformada por a) un acuífero costero detrítico, poroso y libre, albergado en la planicie aluvial de Jacó y, b) un acuífero fracturado en las rocas basálticas del Complejo de Nicoya, libre en la zona montañosa y subyaciendo la formación detrítica en la zona central, donde el agua subterránea se integra a la del acuífero detrítico, conformándose en un mismo acuífero. La zona del acuífero aluvial obtuvo una vulnerabilidad moderada, en tanto que el fracturado una vulnerabilidad baja. Se propusieron recomendaciones de manejo.</p>
			</abstract>
			<trans-abstract xml:lang="en">
				<title>Abstract</title>
				<p>In the last two decades, the Central Pacific coast of Costa Rica experienced a strong growth in urban and touristic development. If not regulated, this development can result in a pressure to groundwater resources and, therefore, a threat in terms of quality and quantity. In this study, a hydrogeological vulnerability map of Jacó aquifers was performed, by the GOD method. Geological profiles were elaborated, and field verification was performed. It is confirmed that, in the study area, two aquifers are found: a) a detritic, porous, unconfined coastal aquifer, settled in the Jacó alluvial plain and, b) an aquifer in the basaltic fractured rocks of the Nicoya Complex, unconfined in the mountainous zone and underlying the detritic formation in the central area, where the groundwater is integrated to that of the detrital aquifer. The zone of the alluvial aquifer resulted in a medium vulnerability index, whereas the mountainous areas in a low vulnerability. Management recommendations were proposed.</p>
			</trans-abstract>
			<kwd-group xml:lang="es">
				<title>Palabras clave:</title>
				<kwd>acuífero costero</kwd>
				<kwd>vulnerabilidad hidrogeológica</kwd>
				<kwd>GOD</kwd>
				<kwd>acuífero fracturado</kwd>
				<kwd>protección del agua subterránea.</kwd>
			</kwd-group>
			<kwd-group xml:lang="en">
				<title>Keywords:</title>
				<kwd>coastal aquifer</kwd>
				<kwd>hydrogeological vulnerability</kwd>
				<kwd>GOD index</kwd>
				<kwd>fissured aquifer</kwd>
				<kwd>groundwater protection.</kwd>
			</kwd-group>
			<counts>
				<fig-count count="16"/>
				<table-count count="0"/>
				<equation-count count="0"/>
				<ref-count count="23"/>
				<page-count count="23"/>
			</counts>
		</article-meta>
	</front>
	<body>
		<sec sec-type="intro">
			<title>Introducción</title>
			<p>A nivel mundial, el acceso al agua potable se ha visto amenazado debido a un aumento en la demanda para satisfacer una mayor cantidad y variedad de actividades antropogénicas, una carga contaminante mayor y un uso no racional del recurso (<xref ref-type="bibr" rid="B20">Llamas y Custodio, 2003</xref>). Con base en los modelos de desarrollo (<xref ref-type="bibr" rid="B13">Foster, 1992</xref>), el Banco Mundial (BM) ha recomendado la elaboración de las cartografías de amenaza, vulnerabilidad y riesgo de los acuíferos a la contaminación como una medida para la protección de la calidad de las aguas subterráneas (<xref ref-type="bibr" rid="B17">Foster <italic>et al</italic>., 2007</xref>). El concepto de vulnerabilidad se remonta desde la década de los 70 (<xref ref-type="bibr" rid="B1">Albinet y Margat, 1970</xref>), de forma intuitiva; posteriormente, su uso fue ampliamente utilizado en los años 80 (<xref ref-type="bibr" rid="B15">Foster e Hirata 1988</xref>). De acuerdo con <xref ref-type="bibr" rid="B17">Foster <italic>et al</italic>. (2007)</xref>, el concepto de vulnerabilidad del agua subterránea está en función de las características propias de los estratos que separan la zona no saturada de la saturada. Estas características determinan la tendencia o probabilidad de que algún contaminante alcance una posición especifica en el sistema acuífero, después de su introducción en algún punto sobre el terreno.</p>
			<p>En Costa Rica, durante las últimas décadas, han surgido y se han consolidado políticas ambientales orientadas a reforzar la protección y uso racional de los recursos naturales. No obstante, la gestión del recurso hídrico enfrenta grandes desafíos, los cuales deben superarse con acciones eficientes para garantizar la calidad ambiental y la sostenibilidad del recurso hídrico subterráneo (<xref ref-type="bibr" rid="B23">Reynolds y Fraile, 2002</xref>). Por tal motivo, a partir del 2008, el Estado dictaminó la elaboración de estudios de vulnerabilidad hidrogeológica para ser incluidos en la normativa urbanística de cada cantón del país, según el método GOD (Groundwater Hidraulic Confinement, Overlaying Strata and Deph to Groundwater Table). El método GOD considera dos factores básicos: el grado de inaccesibilidad hidráulica de la zona saturada y la capacidad de atenuación de los estratos suprayacentes a la zona saturada (<xref ref-type="bibr" rid="B14">Foster e Hirata, 1991</xref>). Esta metodología utiliza los siguientes parámetros (<xref ref-type="bibr" rid="B17">Foster <italic>et al</italic>. 2007</xref>): grado de confinamiento hidráulico (G), características de los sustratos de la zona no saturada o capas confinantes, tomando en cuenta aspectos litológicos y el grado de consolidación (O) y; distancia al agua, considerada como la profundidad al nivel freático o al techo del acuífero confinado (D). Se entiende por grado de confinamiento hidráulico si el acuífero es libre, confinado o semiconfinado (<xref ref-type="bibr" rid="B17">Foster <italic>et al.</italic> 2007</xref>).</p>
			<p>El abastecimiento de agua potable ha sido un factor clave para el desarrollo de las comunidades costeras. Sin embargo, el crecimiento desordenado en la zona marítimo-terrestre en el litoral pacífico de Costa Rica pone en evidencia el seguimiento limitado a las regulaciones existentes (<xref ref-type="bibr" rid="B12">Estado de la Nación, 2007</xref>). Además, en los últimos años, se ha generado un rápido crecimiento urbano y turístico en esta zona (<xref ref-type="bibr" rid="B8">Barrantes, 2004</xref>). Según el Instituto Nacional de Estadísticas y Censos (<xref ref-type="bibr" rid="B18">INEC, 2012</xref>), el cantón de Garabito, al cual pertenece la comunidad de Jacó, ha presentado la mayor tasa de crecimiento poblacional del país durante los últimos 30 años, lo que podría conllevar a un aumento en la carga de contaminantes que se liberan al ambiente (<xref ref-type="bibr" rid="B21">Orozco, 2015</xref>). </p>
			<p>Por su parte, lo somero del nivel freático, en algunos sectores (<xref ref-type="bibr" rid="B6">Arias y Morera, 2002</xref>), implica que las actividades antrópicas representen una amenaza a la calidad del agua subterránea. Debido al rápido crecimiento urbano y turístico de la ciudad de Jacó, es necesario contar con una herramienta que contribuya con el ordenamiento territorial para la protección de los recursos hídricos subterráneos, los cuales son ampliamente utilizados para consumo humano y otras actividades productivas. Por esta razón, el objetivo del presente estudio es elaborar una propuesta de Mapa de Vulnerabilidad Hidrogeológica (MVH) en este sector.</p>
		</sec>
		<sec>
			<title>Área de estudio</title>
			<p>El sitio de estudio (<xref ref-type="fig" rid="f1">Figura 1</xref>) está localizado en el distrito de Jacó, cantón de Garabito de Puntarenas, Costa Rica, entre las coordenadas geográficas 9°35'40,66&quot; - 9°40'34,647&quot; latitud norte y 84°38'48,466&quot; - 84°34'27,373&quot; longitud oeste, y entre las coordenadas 1 061 000 m - 1 070 000 m norte y 429 000 m - 437 000 m este, del sistema de Proyección Transversal de Mercator para Costa Rica (CRTM). La zona de estudio está delimitada por una cuenca hidrográfica de nombre Jacó, nombre que ha sido asignado en este estudio para el área en análisis, con una extensión de 47 km<sup>2</sup>; compuesta por tres microcuencas que pertenecen a los ríos Pueblo Nuevo (Mona), Copey y Naranjal. En la parte baja y media de estas microcuencas se localiza un acuífero albergado en depósitos no consolidados, con un área aproximada de 15 km<sup>2</sup>. El área restante está compuesta, principalmente, por basaltos y rocas sedimentarias, en donde se encuentra un acuífero fracturado. Estas rocas se ubican en los cerros (<xref ref-type="fig" rid="f2">Figura 2</xref>).</p>
			<p>
				<fig id="f1">
					<label>Figura 1</label>
					<caption>
						<title>Mapa de ubicación de la zona de estudio, Jacó, Puntarenas, Costa Rica</title>
					</caption>
					<graphic xlink:href="2215-2563-rge-63-123-gf1.jpg"/>
				</fig>
			</p>
			<p>
				<fig id="f2">
					<label>Figura 2</label>
					<caption>
						<title>Mapa geológico de la zona de estudio, Jacó, Puntarenas, Costa Rica. Modificado de <xref ref-type="bibr" rid="B5">Arias (2003</xref>). Redefinición de la Formación Tulín (Maastrichtiano-Eoceno Inferior) del Pacífico Central de Costa Rica</title>
					</caption>
					<graphic xlink:href="2215-2563-rge-63-123-gf2.jpg"/>
				</fig>
			</p>
		</sec>
		<sec>
			<title>Características generales</title>
			<p>La época seca se extiende por cinco meses, comprendidos de diciembre a abril, siendo enero el más seco, en tanto que la época lluviosa comprende los meses de mayo a noviembre y el mes más lluvioso corresponde con octubre. Según la base de datos del Instituto Meteorológico Nacional de Costa Rica del 2014, en la Estación Jacó, la precipitación máxima promedio, entre 1997 y 2013, fue de 616,6 mm en los meses de octubre; el mínimo de precipitación se obtuvo en los meses de enero, con un promedio de 25 mm. </p>
		</sec>
		<sec>
			<title>Aspectos geológicos</title>
			<p>El área de estudio está conformada por tres tipos de unidades litoestatigráficas claramente diferenciables entre sí: a) basaltos del Complejo de Nicoya, dentro del cual se han incluido rocas tipo radiolaritas, calizas, areniscas y lutitas, b) basaltos asociados a la Formación Tulín y, c) depósitos no consolidados (<xref ref-type="fig" rid="f2">Figura 2</xref>). </p>
			<p><italic>Complejo de Nicoya.</italic> Dengo (1962) define el Complejo de Nicoya como una serie de rocas intensamente plegadas, consistentes en basaltos, unidades sedimentarias e intrusivas. <xref ref-type="bibr" rid="B19">Kuypers (1980, citado en Sprechmann et al., (1984)</xref> indica que existe ambigüedad debido a que las series de rocas suprayacentes se encuentran, localmente, muy plegadas y a la presencia de extrusiones basálticas. Por lo tanto, se incluyen diferentes litologías dentro de la unidad y, en consecuencia, varía la edad que se le atribuye. Debido a lo anterior, <xref ref-type="bibr" rid="B19">Kuypers (1980, citado en Sprechmann <italic>et al</italic>., 1984)</xref> y <xref ref-type="bibr" rid="B9">Baumgartner (1984</xref>) proponen una subdivisión estratigráfica para el Complejo de Nicoya, la que consiste en dos unidades, Matapalo y Esperanza. <xref ref-type="bibr" rid="B5">Arias (2003</xref>) denomina a los basaltos que han sido asociados al Complejo de Nicoya como basaltos del <italic>Plateau</italic> Oceánico Caribe (POC); estos basaltos son las rocas más antiguas. </p>
			<p>La Unidad Matapalo posee 150 m de espeso; se caracteriza por estar conformada por basaltos tholeiíticos oceánicos cubiertos por radiolaritas, con espesor mínimo de 40 m. La edad asignada es del Calloviense hasta el Aptiense (<xref ref-type="bibr" rid="B19">Kuypers, 1980, citado en Sprechmann <italic>et al</italic>., 1984)</xref>. El contacto con la suprayacente Esperanza no es claramente visible (<xref ref-type="bibr" rid="B19">Kuypers, 1980, citado en Sprechmann <italic>et al</italic>., 1984</xref>). En la Unidad Esperanza, conformada por basaltos tholeiíticos y diabasas ofiolíticas, son comunes las rocas gabroicas, stocks plagiograníticos y algunas intercalaciones de radiolaritas de hasta 3 m en la parte superior de la unidad. La edad asignada es del Santoniense (<xref ref-type="bibr" rid="B19">Kuypers, 1980, citado en Sprechmann <italic>et al</italic>., 1984</xref>). </p>
			<p>En la zona de estudio, la geología en los cerros está compuesta de basaltos y rocas sedimentarias tipo areniscas, lutitas, radiolaritas y calizas que se asocian a las rocas del Complejo de Nicoya (<xref ref-type="bibr" rid="B6">Arias y Morera, 2002</xref>). Estas lavas, de composición basáltica, son de una coloración negra y, por su origen, tienen una textura en almohadilla, lo que las caracteriza. En superficie, están meteorizadas y son fácilmente confundibles con arenas gruesas. </p>
			<p>Las radiolaritas del Complejo de Nicoya forman parte de la Unidad Matapalo, del Jurásico Superior - Cretácico Inferior. Son rocas de origen pelágico, de grano fino a muy fino, estratificadas, bastante duras y versicolor (<xref ref-type="bibr" rid="B2">Agudelo, 2005</xref>). Están relacionadas a pedernales y jaspes y se encuentran muy a menudo asociadas a mineralizaciones de hierro-manganeso, resultado de la alteración hidrotermal en fondos oceánicos a partir de vulcanismo exhalativo, sea en forma diseminada o en nódulos. Sobreyace disconformemente la Unidad de basaltos y es atravesada por numerosos cuerpos intrusivos gabroides a diabásicos. Presenta una fuerte deformación tectónica (pliegues isoclinales-chevron) y un alto grado de fallamiento <xref ref-type="bibr" rid="B2">(Agudelo, 2005</xref>). </p>
			<p><italic>Formación Tulín.</italic><xref ref-type="bibr" rid="B5">Arias (2003</xref>) sugiere que la Formación Tulín, del (Maastrichtiano-Eoceno Inferior, está conformada en un 95 % por basaltos, los cuales se presentan como coladas en almohadillas; estos basaltos presentan intercalaciones de sedimentos epiclásicos, piroclasto, gabros y cumulitas de olivino. Es una formación más joven que el POC, cuyas rocas constituyen el basamento de la Formación Tulín. </p>
			<p><italic>Depósitos no consolidados.</italic> Pertenecen a esta unidad todos los depósitos recientes asociados a los depósitos por gravedad al pie de monte o coluvios y a aquellos relacionados con procesos aluviales. También aquellos depósitos de materiales arenosos con influencia marina. Los depósitos asociados al arrastre por acción fluvial o por gravedad, por lo general, son fragmentos de rocas de diversos tamaños y composición inmersos en una matriz limo-arcillo arenosa. Hacia el pie de los cerros, los bloques de fragmentos son de mayor tamaño y sub-redondeados a redondeados, en tanto que hacia la parte baja de la cuenca los fragmentos son más pequeños, ya que han sido arrastrados a mayor distancia, lo que se evidencia por su menor grado de angulosidad (<xref ref-type="bibr" rid="B6">Arias y Morera, 2002</xref>). </p>
		</sec>
		<sec>
			<title>Marco teórico-conceptual</title>
			<p>La estimación del índice de vulnerabilidad GOD (<xref ref-type="bibr" rid="B15">Foster e Hirata, 1988</xref>) involucra una serie de etapas concretas. Primera, identificar el grado de confinamiento hidráulico del acuífero y asignarle un índice a este parámetro en una escala de 0.0 a 1.0. Segunda, especificar las características del sustrato suprayacente a la zona saturada del acuífero en términos de: a) grado de consolidación (teniendo en cuenta la probable presencia o ausencia de permeabilidad por fisuras) y, b) tipo de litología (considerando indirectamente porosidad efectiva, permeabilidad de la matriz y contenido de humedad en la zona no saturada o retención específica) y, asignar un índice a este parámetro en una escala de 0.4 a 1.0. Tercera, estimar la distancia o profundidad al nivel del agua (en acuíferos no confinados) o profundidad al techo del primer acuífero confinado, con la consiguiente asignación de un índice en una escala de 0.6 a 1.0. En este sentido, el peligro de contaminación del agua subterránea se definiría luego como la probabilidad de que el agua subterránea, en la parte superior de un acuífero, sea contaminada en un nivel inaceptable por las actividades que se desarrollan en la superficie del terreno suprayacente (<xref ref-type="bibr" rid="B15">Foster e Hirata, 1988</xref>; <xref ref-type="bibr" rid="B17">Foster <italic>et al.</italic> 2007</xref>). </p>
			<p>Así, en lugar de aplicar controles universales sobre los usos del territorio potencialmente contaminantes y la descarga de efluentes, es más efectivo (y menos perjudicial para el desarrollo económico) modificar el tipo y nivel de control de acuerdo con esta capacidad de atenuación. Esta es la premisa básica detrás del concepto de vulnerabilidad a la contaminación de acuíferos y lo que conduce a la necesidad de contar con el mapeo de la misma (<xref ref-type="bibr" rid="B17">Foster <italic>et al</italic>., 2007</xref>).</p>
		</sec>
		<sec sec-type="methods">
			<title>Marco metodológico</title>
			<p>Perfiles geológicos. En la zona de estudio existen registrados 57 pozos privados y cuatro pozos de abastecimiento público pertenecientes al AyA. De estos 57, solamente 21 cuentan con la descripción litológica en sus informes de perforación. Estos informes se recuperaron en octubre del 2016 del archivo institucional de la Dirección de Investigación y Gestión Hídrica del Servicio Nacional de Aguas Subterráneas, Riego y Avenamiento (SENARA) y del Archivo de Concesiones de la Dirección de Aguas del Ministerio de Ambiente y Energía (MINAE). Con base en esta información, se trazaron al menos 11 perfiles geológicos, de los que se escogieron los cuatro perfiles más representativos para determinar la relación entre las unidades geológicas. Un perfil se trazó con una dirección de suroeste a noreste, el segundo perfil de noroeste- sureste, atravesando el sector de la cuenca baja y el tercero de noroeste-sureste, atravesando el sector de la cuenca alta.</p>
			<p>Aforos diferenciales. Para definir la relación hidráulica entre las aguas superficiales y las subterráneas se planificaron dos campañas de aforos diferenciales, una en época seca y otra en la época lluviosa, en los principales ríos y quebradas de la cuenca. Los aforos diferenciales han sido ampliamente utilizados en diversas situaciones climáticas para realizar estimaciones de caudal, por ejemplo, en <xref ref-type="bibr" rid="B4">Andreo <italic>et al</italic>. (2004</xref>) y <xref ref-type="bibr" rid="B7">Barberá y Andreo (2015</xref>). Para la selección de los sitios se consideró que éstos no tuvieran influencia por entradas o salidas de aguas provenientes de otros cursos de agua, y si las había, se aforaban estos ingresos o egresos de agua para valorar la relación hidráulica entre el río y el acuífero. En noviembre del 2013, durante la época lluviosa, se realizaron aforos en 18 sitios; en marzo del 2014, época seca, en 12 sitios. La diferencia entre ambas campañas se debe a que seis de los puntos aforados en época lluviosa se encontraban secos para la época de estiaje. Una vez calculado el caudal en cada punto aforado, se definieron los tramos donde el río es influente y los tramos en donde el río es efluente.</p>
			<p>Mapa de vulnerabilidad. Se utilizó el método GOD (Foster <italic>et al</italic>., 2002), para caracterizar la vulnerabilidad de los acuíferos en la zona de estudio. Para la obtención de los valores de vulnerabilidad en cada una de las capas de las que se compone el método GOD, se hizo una revisión del Mapa Geológico de Costa Rica (<xref ref-type="bibr" rid="B10">Denyer y Alvarado, 2007</xref>), del Atlas tectónico de Costa Rica (<xref ref-type="bibr" rid="B11">Denyer <italic>et al</italic>., 2003</xref>), del estudio de <xref ref-type="bibr" rid="B5">Arias (2003</xref>) y de la Cartografía digital y ortofoto (<xref ref-type="bibr" rid="B22">PRCR, 2007</xref>). Además, se hizo corroboración de campo de la geología de la zona, correlacionando los afloramientos observados con las descripciones realizadas por los autores mencionados. También, se recopiló la información litológica de los pozos en la zona a partir de los informes de perforación archivados en el SENARA (Servicio Nacional de Aguas Subterráneas Riego y Avenamiento) y en la Dirección de Aguas del Ministerio de Ambiente y Energía, y de los datos de perforación de los pozos del AyA, para estimar el espesor de la zona no saturada y la profundidad del techo de los acuíferos. Finalmente, se utilizaron los resultados de los perfiles geológicos elaborados durante esta investigación para mejorar el entendimiento del modelo hidrogeológico conceptual. Se utilizó un Sistema de Información Geográfica (SIG) para integrar la información requerida para la estimación de los valores en cada capa del método GOD. Una vez definidas las vulnerabilidades en el área de estudio, se elaboró una matriz con las actividades sugeridas en cada zona, según su vulnerabilidad, considerando las actividades desarrolladas en la región. </p>
		</sec>
		<sec sec-type="results">
			<title>Resultados</title>
			<sec>
				<title>-Hidrometría</title>
				<p>Durante las campañas de medición de caudal, se corroboró en el campo la ausencia de tomas de agua para riego, embalses o cualquier otro aprovechamiento a lo largo del río o quebrada. De esta forma, se confirmó el carácter influente-efluente de los ríos estudiados. La quebrada Mona, en época lluviosa, tuvo un tramo en el sector central de su recorrido en donde el caudal pasó de 79,9 a 14,7 l/s. Esta reducción en el caudal indica que en este sector la quebrada se comportó como influente y el tramo que sigue hasta su desembocadura como efluente (<xref ref-type="fig" rid="f3">Figura 3</xref>). El río Copey, en época lluviosa, se comportó como efluente en la mayor parte del trayecto, excepto el sector medio, en donde el caudal pasó de 133,2 a 112,9 l/s. La quebrada Seca, también en época lluviosa, tuvo un comportamiento efluente, excepto en su desembocadura, donde su caudal cambió de 99,6 a 52,4 l/s, convirtiéndose en una quebrada influente en ese sitio (<xref ref-type="fig" rid="f3">Figura 3</xref>). En época de lluvia, el comportamiento influente-efluente de los ríos a lo largo de su cauce evidenció que el acuífero tiene una profundidad al nivel de agua subterránea por encima del lecho del río, de tal manera que el acuífero puede alimentar los ríos y quebradas, mientras que en otros sectores la posición del nivel del agua está por debajo del lecho del río, invirtiendo el comportamiento. En época seca, el comportamiento de los ríos se mantuvo igual, y lo único que varió fue el caudal. Así mismo, los tramos en los que el caudal disminuyó drásticamente en época lluviosa son los mismos que se encontraron completamente secos en época seca debido al descenso del nivel del agua en el acuífero (<xref ref-type="fig" rid="f3">Figura 3</xref>), lo que corrobora el carácter influente de los ríos y quebradas en estos sectores. </p>
				<p>
					<fig id="f3">
						<label>Figura 3</label>
						<caption>
							<title>Mapa de aforos diferenciales en la zona de estudio, Jacó, Puntarenas, Costa Rica</title>
						</caption>
						<graphic xlink:href="2215-2563-rge-63-123-gf3.jpg"/>
					</fig>
				</p>
			</sec>
			<sec>
				<title>-Modelo hidrogeológico</title>
				<p>En la zona de Jacó se han desarrollado dos acuíferos costeros; uno superficial, detrítico, libre, albergado en el depósito no consolidado y, otro fracturado y albergado en las rocas del Complejo de Nicoya (<xref ref-type="fig" rid="f4">Figura 4</xref>). Este acuífero fracturado es libre y profundo en los cerros y, cuando subyace al acuífero detrítico en la llanura aluvial se puede mostrar confinado en al menos un sector de la zona central de la llanura, donde la carga de presión alcanza hasta 1 m por encima del nivel del terreno (registros de perforación HE-134 y Pozo 5 del AyA) (perfil 4, <xref ref-type="fig" rid="f6">Figura 6</xref>). El agua subterránea en las zonas montañosas viaja a través de los basaltos y rocas sedimentarias; luego, probablemente, una parte se integra al agua subterránea del depósito no consolidado y conforma un único acuífero (perfiles 1 y 2, <xref ref-type="fig" rid="f5">Figuras 5</xref> y 6). Aun cuando este acuífero estuviera conformado por dos formaciones geológicas, materiales no consolidados y basaltos por debajo de ellos, con distintos valores de conductividad hidráulica y permeabilidad, es probable que el agua subterránea mantenga un mismo nivel por la ausencia de una capa impermeable que separe ambas formaciones, lo que se deduce a partir de los informes de perforación obtenidos del SENARA. En el informe del pozo HE-134 se indicó, que cuando el pozo alcanzó los basaltos fracturados, a los 22 m de profundidad, se produjo una surgencia (<xref ref-type="fig" rid="f6">Figura 6</xref>, perfil 4). La surgencia reportada en ese informe y en el Pozo 5 del AyA podría deberse a la diferencia en permeabilidad de la formación fracturada con respecto a la detrítica, a una mayor presión del agua en dicho material fracturado debido al efecto de algunos lentes de arcilla cerca de la superficie (<xref ref-type="fig" rid="f6">Figura 6</xref>, perfil 4), o a que los pozos que alcanzaron los basaltos están expuestos a un mayor potencial hidráulico, debido a que alcanzaron profundidades con distintas líneas equipotenciales a las de los pozos que extraen agua de la formación detrítica, sin que estas tres razones sean mutuamente excluyentes. Además, en los perfiles 1 y 2 (<xref ref-type="fig" rid="f5">Figura 5</xref>) se evidencia que la topografía de los basaltos es muy irregular y el nivel del agua subterránea en estos basaltos es afectado por las fuertes pendientes.</p>
				<p>
					<fig id="f4">
						<label>Figura 4</label>
						<caption>
							<title>Mapa hidrogeológico de la zona de estudio, Jacó, Puntarenas, Costa Rica</title>
						</caption>
						<graphic xlink:href="2215-2563-rge-63-123-gf4.jpg"/>
					</fig>
				</p>
				<p>
					<fig id="f5">
						<label>Figura 5</label>
						<caption>
							<title>Perfiles hidrogeológicos 1 y 3 en la zona de estudio, Jacó, Puntarenas, Costa Rica</title>
						</caption>
						<graphic xlink:href="2215-2563-rge-63-123-gf5.jpg"> </graphic>
					</fig>
				</p>
				<p>
					<fig id="f6">
						<label>Figura 6</label>
						<caption>
							<title>Perfiles hidrogeológicos 2 y 4 en la zona de estudio, Jacó, Puntarenas, Costa Rica</title>
						</caption>
						<graphic xlink:href="2215-2563-rge-63-123-gf6.jpg"/>
					</fig>
				</p>
			</sec>
			<sec>
				<title>-Acuífero detrítico Jacó</title>
				<p>Se trata de un acuífero, libre, detrítico, albergado en materiales no consolidados tipo arenas, limos y gravas con una capa de suelo de poco espesor. En los sectores con información de nivel de agua, este se encuentra entre 0,5 y 5 m de profundidad, con equipotenciales entre los 0,5 y 15 m.s.n.m. (<xref ref-type="fig" rid="f4">Figura 4</xref>). La dirección hidráulica de este acuífero es hacia el suroeste (SW) y el gradiente es de 0,01. Otro mecanismo de recarga se da por infiltración desde el lecho de los ríos, como se ha demostrado con los aforos diferenciales realizados (<xref ref-type="fig" rid="f3">Figura 3</xref>). La descarga se produce de manera natural, como descarga a los ríos, en algunos sectores (<xref ref-type="fig" rid="f3">Figura 3</xref>), y hacia al mar (<xref ref-type="fig" rid="f4">Figura 4</xref>), y producto de la extracción de los pozos. En la planicie aluvial no se encontraron manantiales durante las visitas al campo ni tampoco existen registros. Los caudales de explotación indican una productividad variable desde 1 hasta 50 l/s. Los pozos que reportan un caudal de extracción mayor son los pozos que alcanzaron los basaltos fracturados.</p>
			</sec>
			<sec>
				<title>-Acuífero fracturado</title>
				<p>Este acuífero fracturado, con una permeabilidad secundaria originada por las fisuras, se ha desarrollado en lavas y rocas sedimentarias mencionadas anteriormente, las que conforman los cerros que rodean la cuenca. En este acuífero fracturado, existían 19 pozos registrados hasta la fecha; únicamente 15 con información litológica. La profundidad del nivel del agua subterránea es variable debido a que los pozos tienen profundidades variables y el nivel del agua resultante corresponde con un nivel de equilibrio entre las distintas capas productoras. Por ejemplo, el pozo HE-109 posee 105 m de profundidad y un nivel estático reportado a 70 m de profundidad, en tanto que el pozo HE-58 tiene 180 m de profundidad y un nivel estático a 100 m de profundidad. En los informes de los pozos HE-18, HE-19, HE-47, HE-68 y HE-134, en la zona central, se reportó surgencia entre los 22 y los 38 m, en donde se ubican los basaltos fracturados (perfil 4, <xref ref-type="fig" rid="f6">Figura 6</xref>).</p>
				<p>En cuanto a la productividad del acuífero, los informes de perforación obtenidos de los archivos del SENARA indican caudales muy bajos, entre 0,5 a 7 l/s, con descensos muy fuertes. A pesar de no contarse con pruebas de bombeo, esta condición, posiblemente, refleja un acuífero de baja productividad debido a una baja permeabilidad secundaria. Aun cuando no se pudo trazar equipotenciales para este acuífero por falta de información, es probable que el flujo se dirija hacia las zonas aluviales, ya que las rocas que contienen este acuífero se encuentran en los cerros bordeando la llanura. Posteriormente, en el sector subyacente al acuífero detrítico, el agua se integra con el agua subterránea en el acuífero detrítico, conformándose un único nivel y un único acuífero. Luego, el flujo se dirige hacia el SW, es decir, hacia el mar. Sin embargo, esta es una suposición basada en las características topográficas y geomorfológicas de la zona (<xref ref-type="fig" rid="f4">Figura 4</xref>). </p>
			</sec>
			<sec>
				<title>-Vulnerabilidad hidrogeológica</title>
				<p>Una vez definidas las condiciones hidrogeológicas, se estimó la vulnerabilidad del acuífero poroso y del acuífero fracturado en la zona de los cerros. </p>
				<p><italic>Grado de confinamiento hidráulico (Valor G).</italic> Al acuífero detrítico se le asignó un valor de 0,9, por ser un acuífero con poca cobertura, es decir, con una zona vadosa de poco espesor, según lo observado mediante los perfiles realizados. Para el acuífero fracturado se estimó un valor de 0,7, correspondiente a un acuífero libre con cobertura de acuerdo a los perfiles 1, 2 y 3. (<xref ref-type="fig" rid="f7">Figura 7</xref>). <italic>Ocurrencia del sustrato suprayacente (Valor O).</italic> Con base en la descripción litológica en los reportes de perforación de los pozos, se identificó una capa de suelo arcillo-limoso en la zona del acuífero aluvial, de espesor variable entre los 0 y 14 m, por lo que se estimó un valor de 0,5. En los cerros, el valor es de 0,6, corresponde a rocas basálticas fracturadas (<xref ref-type="fig" rid="f7">Figura 7</xref>). <italic>Distancia al nivel del agua subterránea (Valor D)</italic>. En el acuífero detrítico, la profundidad del nivel de agua con respecto al nivel de suelo se encuentra entre 1 y 5 m, por lo que se le asignó un valor de 0.9, mientras que en el acuífero libre fracturado, en la zona de los cerros, el agua se encontró entre 20 y 50 m de profundidad, correspondiéndole un valor de 0,7, aunque hay niveles en algunos pozos a mayor profundidad (<xref ref-type="fig" rid="f7">Figura 7</xref>). </p>
				<p>
					<fig id="f7">
						<label>Figura 7</label>
						<caption>
							<title>Grado de confinamiento (G), Ocurrencia del sustrato (O) y Distancia del agua (D) de los acuíferos en la zona de estudio, Jacó, Puntarenas, Costa Rica</title>
						</caption>
						<graphic xlink:href="2215-2563-rge-63-123-gf7.jpg"/>
					</fig>
				</p>
				<p>Como resultado, se obtuvo un nivel medio de vulnerabilidad hidrogeológica en la zona en donde se ubica el acuífero detrítico, con un valor de 0,41 (<xref ref-type="fig" rid="f8">Figura 8</xref>). Este sector representó un 30,5 % del total del área de estudio, mientras que el área correspondiente del acuífero fracturado se clasificó con un nivel bajo de vulnerabilidad hidrogeológica, con un valor de 0,29.</p>
				<p>
					<fig id="f8">
						<label>Figura 8</label>
						<caption>
							<title>Propuesta de Vulnerabilidad Hidrogeológica de los acuíferos en la zona de estudio, Jacó, Puntarenas, Costa Rica</title>
						</caption>
						<graphic xlink:href="2215-2563-rge-63-123-gf8.jpg"/>
					</fig>
				</p>
			</sec>
		</sec>
		<sec sec-type="discussion">
			<title>Discusión</title>
			<p>Si bien la vulnerabilidad a la contaminación depende de la capacidad de atenuación natural del perfil de suelo a los contaminantes, esta varía ampliamente con las condiciones geológicas cercanas a la superficie del terreno. Entre estas condiciones se encuentra la profundidad a la que está el agua subterránea. Esto tuvo implicaciones importantes al designarle una vulnerabilidad baja al acuífero fracturado, ubicado en los cerros de la zona de Jacó, debido a que el nivel del agua se encontró a profundidades superiores a los 50 m. Los acuíferos con vulnerabilidad baja son vulnerables a contaminantes conservativos cuando son descargados o lixiviados a largo plazo, con regularidad (<xref ref-type="bibr" rid="B17">Foster <italic>et al</italic>., 2007</xref>). Según Foster <italic>et al.</italic> (2013), un acuífero con vulnerabilidad baja puede ser vulnerable a la lixiviación de nitratos, líquidos altamente salinos y compuestos organoclorados en fase densa no acuosa.</p>
			<p>Por su parte, aun cuando el nivel del agua subterránea es somero, el acuífero sobre la llanura aluvial obtuvo una clasificación de vulnerabilidad a la contaminación moderada. Esto debido a que los informes de los pozos en esta zona indican una capa de arcillas de varios metros de espesor en superficie, e inclusive indican intercalaciones de arcillas con las gravas dentro del paquete aluvial, lo que le ha dado a este acuífero cierto grado de protección ante ciertas sustancias contaminantes. Según <xref ref-type="bibr" rid="B17">Foster <italic>et. al.</italic> (2007</xref>), los acuíferos con vulnerabilidad moderada son aquellos que son vulnerables a algunos contaminantes cuando son descargados o lixiviados con regularidad. Aun cuando hay un pequeño sector dentro del acuífero en donde se reporta surgencia, los informes de los pozos no señalan la profundidad a la que ocurrió esta condición y los niveles de agua subterránea indicados en dichos informes corresponden con niveles de equilibrio entre el sector confinado y el acuífero detrítico más superficial. Por tal motivo, este sector se clasificó bajo la misma clase de vulnerabilidad que el resto del acuífero detrítico. </p>
			<p>La inclusión de la vulnerabilidad hidrogeológica en el Plan Regulador de los gobiernos locales en Costa Rica resultaría de gran provecho para definir las actividades a desarrollar dentro del territorio municipal, lo que permitiría implementar procesos de gestión integrada del recurso hídrico. Para esto, la implementación de la propuesta del mapa de vulnerabilidad debe ir asociada a la aplicación de una matriz de vulnerabilidad. En el caso de Jacó, la mayor parte de los proyectos urbanísticos se ubican en el sector de vulnerabilidad moderada, así como las principales actividades productivas y socioeconómicas del área de estudio. Sin embargo, el plan regulador para la zona costera de Jacó, el cual no ha sido aprobado a la fecha, no contempla aún la vulnerabilidad del acuífero y permite un área de construcción por predio del 70 %, no siendo necesariamente la medida adecuada para un sector con vulnerabilidad moderada, tanto por razones de recarga como de carga hidráulica potencialmente contaminante. Tampoco existe una normativa que regule el tamaño de los predios a futuro en zonas no urbanizadas dentro del distrito. Además, se debe evaluar qué cantidad de predios actuales cumple la norma y cuáles serán las medidas de control a futuro. </p>
			<p>La elaboración de los mapas de vulnerabilidad en la zona evidenció una necesidad de mayor inversión por parte de las instituciones gubernamentales para generar datos hidrogeológicos, principalmente a través de la instalación de más piezómetros en las áreas del Complejo de Nicoya. Esta indicación es necesaria para cualquier método de evaluación de la vulnerabilidad hidrogeológica (<xref ref-type="bibr" rid="B16">Foster <italic>et al</italic>., 2013</xref>). Debido a la escasez de información hidrogeológica en la zona de estudio, se debió extrapolar parte de la información a otros sectores. Esto ocurrió en sectores, que geológicamente se asocian a las rocas del Complejo de Nicoya, con pocos datos sobre la profundidad del nivel del agua y litología. Sin embargo, afortunadamente, la mayor proporción de estos terrenos con poca información se ubican en zonas con pendientes mayores al 40 %, recibiendo una protección ya definida en la Ley Forestal y, por lo tanto, no son susceptibles a ser desarrolladas. Por su parte, debido a la escaza información disponible en la zona, la metodología que se ajustó mejor a este estudio fue la del método GOD. De igual forma, se recomienda ejecutar, en la medida en que se genere más información, otros métodos de mayor precisión, como DRASTIC, propuesto por <xref ref-type="bibr" rid="B3">Aller <italic>et al.</italic> (1987</xref>), entre otros. Sin embargo, este método tiende a subestimar la vulnerabilidad de acuíferos fracturados (Foster <italic>et al.</italic> 2007), y otros métodos más complejos requieren de una mayor inversión en cuanto a la generación de insumos, como la tasa de percolación o de recarga acuífera. Además, según Foster <italic>et al.</italic> (2013), no necesariamente los métodos más complejos son los que generan mapas de vulnerabilidad más cercanos a la realidad. </p>
			<p>En Costa Rica, los planes reguladores de ordenamiento territorial no toman en cuenta la amenaza y el riesgo a la contaminación y solamente incorporan la cartografía de vulnerabilidad hidrogeológica como una herramienta de gestión del recurso hídrico subterráneo. El mapa de riesgo de contaminación, al integrar la amenaza de contaminación y la vulnerabilidad de los acuíferos, permitiría regular el territorio, tomando en cuenta las amenazas actuales y, contribuiría con la identificación de zonas de conflicto en donde existe alguna vulnerabilidad, con el propósito de diseñar o corregir medidas regulatorias para la protección del agua subterránea.</p>
		</sec>
		<sec>
			<title>Recomendaciones</title>
			<p>Tomando en cuenta el nivel de vulnerabilidad estimada, se generó una serie de recomendaciones de manejo y protección del suelo en la zona de estudio. Aun cuando se cuenta con información litológica limitada y se carece de pruebas de infiltración en el acuífero fracturado, se recomienda tomar varias medidas de protección de manera precautoria hasta que se conozca mejor la zona no saturada sobre el acuífero fracturado. </p>
			<p>Es así como, en la llanura aluvial, la cual se clasificó como zona de Vulnerabilidad Hidrogeológica Moderada, se justifica la necesidad de la instalación de un sistema de alcantarillado sanitario, el cual cuente con un sistema de tratamiento y un emisario submarino, debido a la lixiviación potencial de sustancias contaminantes desde los tanques sépticos y letrinas. Se recomienda la introducción de estas medidas de manera paulatina con apoyo del gobierno local, del AyA y del Consejo Cantonal de Coordinación Interinstitucional (CCCI) de Jacó. Por su parte, las plantas de tratamiento de aguas residuales mejoradas son obligatorias para las actividades hoteleras, turísticas e industriales de alto impacto, en concordancia con el concepto de armonía con el entorno, dirigido al ecoturismo y la sostenibilidad ambiental.</p>
			<p>Por su parte, en el sector de vulnerabilidad moderada, es decir, sobre el acuífero detrítico, se recomienda que el cambio de uso de suelo de pasto y cultivos a urbano sea permitido, mientras se cumpla con las indicaciones sobre la disposición de aguas servidas mencionadas anteriormente. Además, se recomienda que los predios tengan un área mínima de 300 m<sup>2</sup>, con un área de construcción máxima del 50 %, para evitar la impermeabilización de los terrenos y favorecer la recarga del acuífero. De conservarse las zonas de cultivo en la zona del acuífero detrítico, se recomienda incentivar prácticas de manejo agrícola mejoradas, como el uso racional de agroquímicos, la agricultura orgánica, cosechas controladas para el descanso de los suelos y proyectos de reforestación. En esta zona, no se recomiendan los rellenos sanitarios, canteras y/o tajos. </p>
			<p>Así mismo, en el sector de vulnerabilidad baja, sobreyaciendo el acuífero fracturado, las regulaciones pueden ser más flexibles y permisivas. En cuanto al desarrollo urbanístico, se recomiendan predios de un mínimo de 250 m<sup>2</sup> y con un área de construcción máxima del 50 %, siempre y cuando se encuentren en terrenos que cuenten con una pendiente menor al 40 %. Al igual que en las zonas de Vulnerabilidad Moderada, la agricultura y el pastoreo deben de seguir prácticas de manejo de suelos adecuadas y en lo posible incentivar a la reforestación y a la agricultura orgánica. Las actividades turísticas y hoteleras deberían encaminarse al ecoturismo. Por último, se recomienda que sean permitidos: los rellenos sanitarios (con un diseño adecuado a las condiciones), los tajos, canteras y la industria con plantas de tratamientos mejoradas. Para los terrenos con una pendiente mayor al 40 %, se recomienda un uso de suelo forestal, el cual se ha designado como tal en la Ley Forestal de Costa Rica para, entre otras cosas, minimizar la erosión.</p>
			<p>Finalmente, se recomienda implementar o seguir ejecutando otras herramientas de gestión integral del recurso hídrico: incentivos económicos para promover cambios hacia actividades con menor potencial contaminante en zonas con vulnerabilidad moderada y elevada; redes de monitoreo para la calidad del agua subterránea; e implementación de programas de concientización social derivados de un manejo de cuencas conjunto entre cantones y localidades, apoyados por universidades, las municipalidades, el gobierno y la empresa privada.</p>
		</sec>
	</body>
	<back>
		<ref-list>
			<title>Referencias</title>
			<ref id="B1">
				<mixed-citation>Albinet, M. y Margat, J. (1970) Cartographie de la vulnerabilite a la pollution des nappes d’eau souterraine (Contamination vulnerability mapping of groundwater). <italic>Bulletin de la Bureau de Recherches Géologiques et Minieres</italic>. Segunda serie 3(4):13- 22.</mixed-citation>
				<element-citation publication-type="journal">
					<person-group person-group-type="author">
						<name>
							<surname>Albinet</surname>
							<given-names>M.</given-names>
						</name>
						<name>
							<surname>Margat</surname>
							<given-names>J.</given-names>
						</name>
					</person-group>
					<year>1970</year>
					<article-title>Cartographie de la vulnerabilite a la pollution des nappes d’eau souterraine (Contamination vulnerability mapping of groundwater)</article-title>
					<source>Bulletin de la Bureau de Recherches Géologiques et Minieres</source>
					<comment>Segunda serie</comment>
					<volume>3</volume>
					<issue>4</issue>
					<fpage>13</fpage>
					<lpage> 22</lpage>
				</element-citation>
			</ref>
			<ref id="B2">
				<mixed-citation>Agudelo, C. (2005). <italic>Diagnóstico de las aguas subterráneas en las hojas Tárcoles, Herradura y Candelaria</italic>. Servicio Nacional de Aguas Subterráneas, Riego y Avenamiento (SENARA), San José. pp. 35 (Informe interno).</mixed-citation>
				<element-citation publication-type="book">
					<person-group person-group-type="author">
						<name>
							<surname>Agudelo</surname>
							<given-names>C.</given-names>
						</name>
					</person-group>
					<year>2005</year>
					<source>Diagnóstico de las aguas subterráneas en las hojas Tárcoles, Herradura y Candelaria</source>
					<publisher-name>Servicio Nacional de Aguas Subterráneas, Riego y Avenamiento (SENARA)</publisher-name>
					<publisher-loc>San José</publisher-loc>
					<fpage>35</fpage>
					<lpage>35</lpage>
					<comment>Informe interno</comment>
				</element-citation>
			</ref>
			<ref id="B3">
				<mixed-citation>Aller, L.; Bennett, T.; Lehr, J.; Petty, R. y Hackett, G. (1987). DRASTIC: <italic>A standardized system for evaluating geroundwater pollution potential using hydrogeologic settings</italic>. Environmental Protection Agency Report 600/2-87-035. Washington, D.C.</mixed-citation>
				<element-citation publication-type="book">
					<person-group person-group-type="author">
						<name>
							<surname>Aller</surname>
							<given-names>L.</given-names>
						</name>
						<name>
							<surname>Bennett</surname>
							<given-names>T.</given-names>
						</name>
						<name>
							<surname>Lehr</surname>
							<given-names>J.</given-names>
						</name>
						<name>
							<surname>Petty</surname>
							<given-names>R.</given-names>
						</name>
						<name>
							<surname>Hackett</surname>
							<given-names>G.</given-names>
						</name>
					</person-group>
					<year>1987</year>
					<source>DRASTIC: <italic>A standardized system for evaluating geroundwater pollution potential using hydrogeologic settings</italic></source>
					<comment>Environmental Protection Agency Report 600/2-87-035</comment>
					<publisher-loc>Washington, D.C.</publisher-loc>
				</element-citation>
			</ref>
			<ref id="B4">
				<mixed-citation>Andreo, B.; Vadillo, I.; Carrasco, F.; Neukum, C.; Jiménez, P.; Goldscheider, N.; Hötzl, H.; Vías, J.M.; Pérez, I. y Göppert, N. (2004). Precisiones sobre el funcionamiento hidrodinámico y la vulnerabilidad a la contaminación del acuífero kárstico de la Sierra de Líbar (provincias de Málaga y Cádiz, Sur de España) a partir de un ensayo de trazadores. <italic>Revista de la Sociedad Geológica de España</italic>, 17(3-4): 187-197.</mixed-citation>
				<element-citation publication-type="journal">
					<person-group person-group-type="author">
						<name>
							<surname>Andreo</surname>
							<given-names>B.</given-names>
						</name>
						<name>
							<surname>Vadillo</surname>
							<given-names>I.</given-names>
						</name>
						<name>
							<surname>Carrasco</surname>
							<given-names>F.</given-names>
						</name>
						<name>
							<surname>Neukum</surname>
							<given-names>C.</given-names>
						</name>
						<name>
							<surname>Jiménez</surname>
							<given-names>P.</given-names>
						</name>
						<name>
							<surname>Goldscheider</surname>
							<given-names>N.</given-names>
						</name>
						<name>
							<surname>Hötzl</surname>
							<given-names>H.</given-names>
						</name>
						<name>
							<surname>Vías</surname>
							<given-names>J.M.</given-names>
						</name>
						<name>
							<surname>Pérez</surname>
							<given-names>I.</given-names>
						</name>
						<name>
							<surname>Göppert</surname>
							<given-names>N.</given-names>
						</name>
					</person-group>
					<year>2004</year>
					<article-title>Precisiones sobre el funcionamiento hidrodinámico y la vulnerabilidad a la contaminación del acuífero kárstico de la Sierra de Líbar (provincias de Málaga y Cádiz, Sur de España) a partir de un ensayo de trazadores</article-title>
					<source>Revista de la Sociedad Geológica de España</source>
					<volume>17</volume>
					<issue>3-4</issue>
					<fpage>187</fpage>
					<lpage>197</lpage>
				</element-citation>
			</ref>
			<ref id="B5">
				<mixed-citation>Arias, O. (2003). Redefinición de la Formación Tulín (Maastrichtiano-Eoceno Inferior) del Pacífico Central de Costa Rica. <italic>Revista Geológica de Centroamérica</italic>. 28. 47-68.</mixed-citation>
				<element-citation publication-type="journal">
					<person-group person-group-type="author">
						<name>
							<surname>Arias</surname>
							<given-names>O.</given-names>
						</name>
					</person-group>
					<year>2003</year>
					<article-title>Redefinición de la Formación Tulín (Maastrichtiano-Eoceno Inferior) del Pacífico Central de Costa Rica</article-title>
					<source>Revista Geológica de Centroamérica</source>
					<volume>28</volume>
					<fpage>47</fpage>
					<lpage>68</lpage>
				</element-citation>
			</ref>
			<ref id="B6">
				<mixed-citation>Arias, M. y Morera, S. (2002). <italic>Evaluación del acuífero de Jaco, cantón de Garabito, Provincia de Puntarenas</italic>. Servicio Nacional de Aguas Subterráneas, Riego y Avenamient (SENARA), San José (Informe interno). pp.15.</mixed-citation>
				<element-citation publication-type="book">
					<person-group person-group-type="author">
						<name>
							<surname>Arias</surname>
							<given-names>M.</given-names>
						</name>
						<name>
							<surname>Morera</surname>
							<given-names>S.</given-names>
						</name>
					</person-group>
					<year>2002</year>
					<source>Evaluación del acuífero de Jaco, cantón de Garabito, Provincia de Puntarenas</source>
					<publisher-name>Servicio Nacional de Aguas Subterráneas, Riego y Avenamient (SENARA)</publisher-name>
					<publisher-loc>San José</publisher-loc>
					<comment>(Informe interno)</comment>
					<fpage>15</fpage>
					<lpage>15</lpage>
				</element-citation>
			</ref>
			<ref id="B7">
				<mixed-citation>Barberá, J. A. y Andreo, B. (2015). Hydrogeological processes in a fluviokarstic area inferred from the analysis of natural hydrogeochemical tracers. The case study of eastern Serranía de Ronda (S Spain). <italic>Journal of Hydrology</italic>, 523, 500-514.</mixed-citation>
				<element-citation publication-type="journal">
					<person-group person-group-type="author">
						<name>
							<surname>Barberá</surname>
							<given-names>J. A.</given-names>
						</name>
						<name>
							<surname>Andreo</surname>
							<given-names>B.</given-names>
						</name>
					</person-group>
					<year>2015</year>
					<article-title>Hydrogeological processes in a fluviokarstic area inferred from the analysis of natural hydrogeochemical tracers. The case study of eastern Serranía de Ronda (S Spain)</article-title>
					<source>Journal of Hydrology</source>
					<volume>523</volume>
					<fpage>500</fpage>
					<lpage>514</lpage>
				</element-citation>
			</ref>
			<ref id="B8">
				<mixed-citation>Barrantes, G. (2004). Disponibilidad del recurso hídrico y sus implicaciones para el desarrollo en Costa Rica. <italic>Undécimo Informe Estado de la Nación en Desarrollo Humano Sostenible, Programa Estado de la Nación</italic>. San José, Costa Rica. (Informe interno).</mixed-citation>
				<element-citation publication-type="journal">
					<person-group person-group-type="author">
						<name>
							<surname>Barrantes</surname>
							<given-names>G.</given-names>
						</name>
					</person-group>
					<year>2004</year>
					<article-title>Disponibilidad del recurso hídrico y sus implicaciones para el desarrollo en Costa Rica</article-title>
					<source>Undécimo Informe Estado de la Nación en Desarrollo Humano Sostenible, Programa Estado de la Nación</source>
					<publisher-loc>San José, Costa Rica</publisher-loc>
					<comment>(Informe interno)</comment>
				</element-citation>
			</ref>
			<ref id="B9">
				<mixed-citation>Baumgartner, P. (1984). El Complejo Ofiolitico de Nicoya (Costa Rica): Modelos estructurales analizados en función de las edades de los radiolarios (Calloviense a Santoniense). Citado en: Sprechmann, P. (1984). Estratigrafía del Complejo de Nicoya. En: Sprechmann, P. (Ed.). <italic>Manual de geología de Costa Rica</italic>, (115-123). Volúmen1: Estratigrafía. Editorial de la Universidad de Costa Rica. San José, Costa Rica.</mixed-citation>
				<element-citation publication-type="book">
					<person-group person-group-type="author">
						<name>
							<surname>Baumgartner</surname>
							<given-names>P.</given-names>
						</name>
					</person-group>
					<year>1984</year>
					<chapter-title>El Complejo Ofiolitico de Nicoya (Costa Rica): Modelos estructurales analizados en función de las edades de los radiolarios (Calloviense a Santoniense). Citado en: Sprechmann, P. (1984). Estratigrafía del Complejo de Nicoya</chapter-title>
					<person-group person-group-type="editor">
						<name>
							<surname>Sprechmann</surname>
							<given-names>P</given-names>
						</name>
					</person-group>
					<source>Manual de geología de Costa Rica</source>
					<fpage>115</fpage>
					<lpage>123</lpage>
					<volume>Volúmen1</volume>
					<comment>Estratigrafía</comment>
					<publisher-name>Editorial de la Universidad de Costa Rica. San José, Costa Rica</publisher-name>
				</element-citation>
			</ref>
			<ref id="B10">
				<mixed-citation>Denyer, P. y Alvarado, G. (2007), <italic>Mapa Geológico de Costa Rica</italic>. - Escala 1: 400000, DGM-ECG-UCR, San José, Costa Rica.</mixed-citation>
				<element-citation publication-type="book">
					<person-group person-group-type="author">
						<name>
							<surname>Denyer</surname>
							<given-names>P.</given-names>
						</name>
						<name>
							<surname>Alvarado</surname>
							<given-names>G.</given-names>
						</name>
					</person-group>
					<year>2007</year>
					<source>Mapa Geológico de Costa Rica</source>
					<comment>- Escala 1: 400000</comment>
					<publisher-name>DGM-ECG-UCR</publisher-name>
					<publisher-loc>San José, Costa Rica</publisher-loc>
				</element-citation>
			</ref>
			<ref id="B11">
				<mixed-citation>Denyer, P.; Montero, W. y Alvarado, G. (2003). <italic>Atlas tectónico de Costa Rica</italic>. - Escala 1: 500000, EUCR, San José, Costa Rica.</mixed-citation>
				<element-citation publication-type="book">
					<person-group person-group-type="author">
						<name>
							<surname>Denyer</surname>
							<given-names>P.</given-names>
						</name>
						<name>
							<surname>Montero</surname>
							<given-names>W.</given-names>
						</name>
						<name>
							<surname>Alvarado</surname>
							<given-names>G.</given-names>
						</name>
					</person-group>
					<year>2003</year>
					<source>Atlas tectónico de Costa Rica</source>
					<comment>- Escala 1: 500000</comment>
					<publisher-name>EUCR</publisher-name>
					<publisher-loc>San José, Costa Rica</publisher-loc>
				</element-citation>
			</ref>
			<ref id="B12">
				<mixed-citation>Estado de la Nación (2007). <italic>Decimotercer Informe Estado de la Nación en Desarrollo Humano Sostenible. Programa Estado de la Nación</italic>. San José, Costa Rica (Inf. interno). pp. 437.</mixed-citation>
				<element-citation publication-type="book">
					<person-group person-group-type="author">
						<collab>Estado de la Nación</collab>
					</person-group>
					<year>2007</year>
					<source>Decimotercer Informe Estado de la Nación en Desarrollo Humano Sostenible. Programa Estado de la Nación</source>
					<publisher-loc>San José, Costa Rica</publisher-loc>
					<comment>(Inf. interno)</comment>
					<fpage>437</fpage>
					<lpage>437</lpage>
				</element-citation>
			</ref>
			<ref id="B13">
				<mixed-citation>Foster, S. (1992). Unsustainable development and irrational exploitation of groundwater resources in developing nations an overview. <italic>IAH Hydrogeology Selected Papers</italic>. Heise, Hannover. 3: 321-336. </mixed-citation>
				<element-citation publication-type="journal">
					<person-group person-group-type="author">
						<name>
							<surname>Foster</surname>
							<given-names>S.</given-names>
						</name>
					</person-group>
					<year>1992</year>
					<article-title>Unsustainable development and irrational exploitation of groundwater resources in developing nations an overview</article-title>
					<source>IAH Hydrogeology Selected Papers</source>
					<publisher-name>Heise</publisher-name>
					<publisher-loc>Hannover</publisher-loc>
					<volume>3</volume>
					<fpage>321</fpage>
					<lpage>336</lpage>
				</element-citation>
			</ref>
			<ref id="B14">
				<mixed-citation>Foster, S. e Hirata, R. (1991). <italic>Determinación del riesgo de contaminación de aguas subterráneas, una metodología basada en datos existentes</italic> (2a ed.). - 90 págs. CEPIS, Centro Panamericano de Ingeniería Sanitaria y Ciencias del Ambiente, Lima, Perú.</mixed-citation>
				<element-citation publication-type="book">
					<person-group person-group-type="author">
						<name>
							<surname>Foster</surname>
							<given-names>S.</given-names>
						</name>
						<name>
							<surname>Hirata</surname>
							<given-names>R.</given-names>
						</name>
					</person-group>
					<year>1991</year>
					<source>Determinación del riesgo de contaminación de aguas subterráneas, una metodología basada en datos existentes</source>
					<comment>(2a ed.).</comment>
					<fpage>90</fpage>
					<lpage>90</lpage>
					<publisher-name>CEPIS, Centro Panamericano de Ingeniería Sanitaria y Ciencias del Ambiente</publisher-name>
					<publisher-loc>Lima, Perú</publisher-loc>
				</element-citation>
			</ref>
			<ref id="B15">
				<mixed-citation>Foster, S. e Hirata, R. (1988). <italic>Groundwater pollution risk assessment: a methodology using available data</italic>. WHO-PAHO-CEPIS, Lima.</mixed-citation>
				<element-citation publication-type="book">
					<person-group person-group-type="author">
						<name>
							<surname>Foster</surname>
							<given-names>S.</given-names>
						</name>
						<name>
							<surname>Hirata</surname>
							<given-names>R.</given-names>
						</name>
					</person-group>
					<year>1988</year>
					<source>Groundwater pollution risk assessment: a methodology using available data</source>
					<publisher-name>WHO-PAHO-CEPIS</publisher-name>
					<publisher-loc>Lima</publisher-loc>
				</element-citation>
			</ref>
			<ref id="B16">
				<mixed-citation>Foster, S.; Hirata, R. y Andreo, B. (2013). The aquifer pollution vulnerability concept: aid or impediment in promoting groundwater protection? <italic>Hydrogeology Journal</italic>, 21, 1389-3922 p.</mixed-citation>
				<element-citation publication-type="journal">
					<person-group person-group-type="author">
						<name>
							<surname>Foster</surname>
							<given-names>S.</given-names>
						</name>
						<name>
							<surname>Hirata</surname>
							<given-names>R.</given-names>
						</name>
						<name>
							<surname>Andreo</surname>
							<given-names>B.</given-names>
						</name>
					</person-group>
					<year>2013</year>
					<article-title>The aquifer pollution vulnerability concept: aid or impediment in promoting groundwater protection?</article-title>
					<source>Hydrogeology Journal</source>
					<volume>21</volume>
					<fpage>1389</fpage>
					<lpage>3922</lpage>
				</element-citation>
			</ref>
			<ref id="B17">
				<mixed-citation>Foster, S.; Hirata, R.; Gomes, D.; D´Elia, M. y Paris, M. (2007). <italic>Protección de la Calidad del Agua Subterránea</italic>. Banco Mundial, Washington, D.C. (pp. 115). </mixed-citation>
				<element-citation publication-type="book">
					<person-group person-group-type="author">
						<name>
							<surname>Foster</surname>
							<given-names>S.</given-names>
						</name>
						<name>
							<surname>Hirata</surname>
							<given-names>R.</given-names>
						</name>
						<name>
							<surname>Gomes</surname>
							<given-names>D.</given-names>
						</name>
						<name>
							<surname>D´Elia</surname>
							<given-names>M.</given-names>
						</name>
						<name>
							<surname>Paris</surname>
							<given-names>M.</given-names>
						</name>
					</person-group>
					<year>2007</year>
					<source>Protección de la Calidad del Agua Subterránea</source>
					<publisher-name>Banco Mundial</publisher-name>
					<publisher-loc>Washington, D.C.</publisher-loc>
					<fpage>115</fpage>
					<lpage>115</lpage>
				</element-citation>
			</ref>
			<ref id="B18">
				<mixed-citation>INEC (2012). <italic>Resultados generales. Censo Nacional de Población y VI de Vivienda 2011</italic>. Instituto Nacional de Estadística y Censos. San José, Costa Rica. Recuperado: <ext-link ext-link-type="uri" xlink:href="www.inec.go.cr">www.inec.go.cr</ext-link> (Consulta 4/03/2015).</mixed-citation>
				<element-citation publication-type="book">
					<person-group person-group-type="author">
						<collab>INEC</collab>
					</person-group>
					<year>2012</year>
					<source>Resultados generales. Censo Nacional de Población y VI de Vivienda 2011</source>
					<publisher-name>Instituto Nacional de Estadística y Censos</publisher-name>
					<publisher-loc>San José, Costa Rica</publisher-loc>
					<ext-link ext-link-type="uri" xlink:href="www.inec.go.cr">www.inec.go.cr</ext-link>
					<comment>(Consulta 4/03/2015)</comment>
				</element-citation>
			</ref>
			<ref id="B19">
				<mixed-citation>Kuypers, E.P. (1980). The geologic history of the Nicoya Ophiolite Complex, Costa rica, and its geotectonic significance. Amsterdam, The Netherlands. Tectonophysics 68(3/4): 233-255. Citado en: Sprechmann, P. (1984). Estratigrafía del Complejo de Nicoya. En: Peter Sprechmann (Ed.). <italic>Manual de geología de Costa Rica</italic>. Volúmen1: Estratigrafía. Editorial de la Universidad de Costa Rica. San José, Costa Rica. pp.320.</mixed-citation>
				<element-citation publication-type="book">
					<person-group person-group-type="author">
						<name>
							<surname>Kuypers</surname>
							<given-names>E.P.</given-names>
						</name>
					</person-group>
					<year>1980</year>
					<chapter-title>The geologic history of the Nicoya Ophiolite Complex, Costa rica, and its geotectonic significance. Amsterdam, The Netherlands. Tectonophysics 68(3/4): 233-255. Citado en: Sprechmann, P. (1984). Estratigrafía del Complejo de Nicoya</chapter-title>
					<person-group person-group-type="editor">
						<name>
							<surname>Sprechmann</surname>
							<given-names>Peter</given-names>
						</name>
					</person-group>
					<source><italic>Manual de geología de Costa Rica</italic>.</source>
					<volume>Volúmen1</volume>
					<comment>Estratigrafía</comment>
					<publisher-name>Editorial de la Universidad de Costa Rica</publisher-name>
					<publisher-loc>San José, Costa Rica</publisher-loc>
					<fpage>320</fpage>
					<lpage>320</lpage>
				</element-citation>
			</ref>
			<ref id="B20">
				<mixed-citation>Llamas, M.R. y Custodio, E. (2003). Intensive use of groundwater: a new situation which demands proactive actions. En: Llamas, M.R.; Custodio, E. (Eds.), <italic>Intensive Use of Groundwater, Challenges and Opportunities</italic> (13-34). A.A. Balkema Publishers. Lisse, Países Bajos.</mixed-citation>
				<element-citation publication-type="book">
					<person-group person-group-type="author">
						<name>
							<surname>Llamas</surname>
							<given-names>M.R.</given-names>
						</name>
						<name>
							<surname>Custodio</surname>
							<given-names>E.</given-names>
						</name>
					</person-group>
					<year>2003</year>
					<chapter-title>Intensive use of groundwater: a new situation which demands proactive actions</chapter-title>
					<person-group person-group-type="editor">
						<name>
							<surname>Llamas</surname>
							<given-names>M.R.</given-names>
						</name>
					</person-group>
					<source>Intensive Use of Groundwater, Challenges and Opportunities</source>
					<fpage>13</fpage>
					<lpage>34</lpage>
					<publisher-name>A.A. Balkema Publishers</publisher-name>
					<publisher-loc>Lisse, Países Bajos</publisher-loc>
				</element-citation>
			</ref>
			<ref id="B21">
				<mixed-citation>Orozco, R. (2015). <italic>Propuesta de manejo de uso y cobertura de la tierra para la reducción del riesgo de contaminación del acuífero costero Jacó, Pacífico Central, Costa Rica</italic>. 56 p. (Tesis de Licenciatura Universidad Nacional), Heredia, Costa Rica. </mixed-citation>
				<element-citation publication-type="thesis">
					<person-group person-group-type="author">
						<name>
							<surname>Orozco</surname>
							<given-names>R.</given-names>
						</name>
					</person-group>
					<year>2015</year>
					<source>Propuesta de manejo de uso y cobertura de la tierra para la reducción del riesgo de contaminación del acuífero costero Jacó, Pacífico Central, Costa Rica</source>
					<fpage>56</fpage>
					<lpage>56</lpage>
					<comment>(Tesis de Licenciatura Universidad Nacional)</comment>
					<publisher-loc>Heredia, Costa Rica</publisher-loc>
				</element-citation>
			</ref>
			<ref id="B22">
				<mixed-citation>PRCR (2007). <italic>Proyecto de Regularización Catastral y Registral. Cartografía digital y ortofoto</italic>. Escala 1:50000. Instituto Geográfico Nacional. San José, Costa Rica.</mixed-citation>
				<element-citation publication-type="book">
					<person-group person-group-type="author">
						<collab>PRCR</collab>
					</person-group>
					<year>2007</year>
					<source>Proyecto de Regularización Catastral y Registral. Cartografía digital y ortofoto</source>
					<comment>Escala 1:50000</comment>
					<publisher-name>Instituto Geográfico Nacional</publisher-name>
					<publisher-loc>San José, Costa Rica</publisher-loc>
				</element-citation>
			</ref>
			<ref id="B23">
				<mixed-citation>Reynolds, J. y Fraile, J. (2002): Presente y futuro de las aguas subterráneas en el Valle Central. En: Reynolds, J. (ed.): <italic>Manejo Integrado de Aguas Subterráneas</italic>. EUNED. San José, Costa Rica. pp. 313.</mixed-citation>
				<element-citation publication-type="book">
					<person-group person-group-type="author">
						<name>
							<surname>Reynolds</surname>
							<given-names>J.</given-names>
						</name>
						<name>
							<surname>Fraile</surname>
							<given-names>J.</given-names>
						</name>
					</person-group>
					<year>2002</year>
					<chapter-title>Presente y futuro de las aguas subterráneas en el Valle Central</chapter-title>
					<person-group person-group-type="editor">
						<name>
							<surname>Reynolds</surname>
							<given-names>J</given-names>
						</name>
					</person-group>
					<source>Manejo Integrado de Aguas Subterráneas</source>
					<publisher-name>EUNED</publisher-name>
					<publisher-loc>San José, Costa Rica</publisher-loc>
					<fpage>313</fpage>
					<lpage>313</lpage>
				</element-citation>
			</ref>
		</ref-list>
	</back>
	<sub-article article-type="translation" id="s1" xml:lang="en">
		<front-stub>
			<article-categories>
				<subj-group subj-group-type="heading">
					<subject>Articles</subject>
				</subj-group>
			</article-categories>
			<title-group>
				<article-title>Hydrogeological vulnerability in Jacó coastal aquifers, Central Pacific, Costa Rica</article-title>
			</title-group>
			<abstract>
				<title>Abstract</title>
				<p>In the last two decades, the Central Pacific of Costa Rica experienced a strong growth in urban and tourism development. If not regulated, this development can result in a pressure to groundwater resources and, therefore, a threat in terms of quality and quantity. In this study, a hydrogeological vulnerability map of Jacó aquifers was performed, by the GOD method. Geological profiles were elaborated, and field verification was performed. It is confirmed that, in the study area, two aquifers are found: a) a detritic, porous, unconfined coastal aquifer, settled in the Jacó alluvial plain and, b) an aquifer in the basaltic fractured rocks of the Nicoya Complex, unconfined in the mountainous zone and underlying the detritic formation in the central area, where the groundwater is integrated to that of the detrital aquifer. The zone of the alluvial aquifer resulted in a medium vulnerability index, whereas the mountainous areas in a low vulnerability. Management recommendations were proposed.</p>
			</abstract>
			<kwd-group xml:lang="en">
				<title>Keywords:</title>
				<kwd>Coastal aquifer</kwd>
				<kwd>hydrogeological vulnerability</kwd>
				<kwd>GOD index</kwd>
				<kwd>fractured aquifer</kwd>
				<kwd>groundwater protection.</kwd>
			</kwd-group>
		</front-stub>
		<body>
			<sec sec-type="intro">
				<title>Introduction</title>
				<p>At the global level, access to drinking water has been threatened due to an increase in demand to satisfy a greater number and variety of anthropogenic activities, greater pollution load and an irrational use of the resource (<xref ref-type="bibr" rid="B20">Llamas and Custodio, 2003</xref>). Based on development models (<xref ref-type="bibr" rid="B13">Foster, 1992</xref>), the World Bank (WB) has recommended the mapping of aquifer threats, vulnerability and risk cartographies, as a measure for protecting groundwater quality (<xref ref-type="bibr" rid="B17">Foster et al., 2007</xref>). The intuitive concept of aquifer vulnerability dates back to the 1970s (<xref ref-type="bibr" rid="B1">Albinet and Margat, 1970</xref>); subsequently, it found widespread use in the 1980s (Foster and Hirata 1988). According to Foster et al. (2007), the concept of groundwater vulnerability depends on the characteristics of the strata that separates the unsaturated zone from the saturated zone. These characteristics determine the tendency or probability of a given pollutant reaching a specific position within the aquifer system, after its introduction on the ground.</p>
				<p>In Costa Rica, in the past last decades, environmental policies aimed at reinforcing the protection and rational use of natural resources have emerged and consolidated. Nevertheless, water resource management faces great challenges, which must be overcome with efficient actions to guarantee environmental quality and sustainability of the groundwater resource (<xref ref-type="bibr" rid="B23">Reynolds and Fraile, 2002</xref>). For this reason, as of 2008, the government of Costa Rica ruled for the development of hydrogeological vulnerability studies and their inclusion in urban development regulations of each canton of the country, according to the GOD method (Groundwater Hydraulic confinement - Overlaying strata - Depth to groundwater table). The GOD method considers two basic factors: the degree of hydraulic inaccessibility of the saturated zone and the attenuation capacity of the overlying strata to the saturated zone (<xref ref-type="bibr" rid="B14">Foster and Hirata, 1991</xref>). This methodology uses the following parameters (<xref ref-type="bibr" rid="B17">Foster et al., 2007</xref>): degree of hydraulic confinement (G), characteristics of substrates of the unsaturated zone or confining layers, taking into account lithological aspects and the degree of consolidation (O) and; distance to water, considered as the depth at the water table or the roof of the confined aquifer (D). The degree of hydraulic confinement refers to whether the aquifer is free, confined or semi-confined (<xref ref-type="bibr" rid="B17">Foster et al., 2007</xref>).</p>
				<p>The supply of drinking water has been a key factor for the development of coastal communities. However, the disorderly growth in the maritime-terrestrial zone of the Central Pacific region of Costa Rica reveals the limited follow-up to existing regulations (<xref ref-type="bibr" rid="B12">Estado de la Nación, 2007</xref>). In addition, in recent years, rapid urban and tourism development has spurred in this area (<xref ref-type="bibr" rid="B8">Barrantes, 2004</xref>). According to the National Institute of Statistics and Censuses (<xref ref-type="bibr" rid="B18">INEC, 2012</xref>), the canton of Garabito, to which the community of Jacó belongs, has presented the highest population growth rate in the country during the last 30 years, which could lead to an increase in the load of pollutants released into the environment (<xref ref-type="bibr" rid="B21">Orozco, 2015</xref>).</p>
				<p>On the other hand, the shallowness of the phreatic level, in some sectors (<xref ref-type="bibr" rid="B6">Arias and Morera, 2002</xref>), implies that anthropic activities represent a threat to the quality of groundwater. Due to the rapid urban and tourism growth in the city of Jacó, it is necessary to have a tool that contributes to territorial ordering for the protection of groundwater resources, which are widely used for human consumption and other productive activities. For this reason, the objective of the present study is to elaborate a Hydrogeological Vulnerability Map (HVM) proposal for this geographical region.</p>
			</sec>
			<sec>
				<title>Study area</title>
				<p>The study site (<xref ref-type="fig" rid="f9">Figure 9</xref>) is located in the district of Jacó, canton of Garabito, province of Puntarenas, Costa Rica, between the geographic coordinates 9°35'40.66&quot;-9°40'34,647&quot; north latitude and 84°38'48,466&quot;-84°34'27,373&quot; west longitude, and between coordinates 1 061 000 m - 1 070 000 m north and 429 000 m - 437 000 m east in the Mercator Transversal Projection System for Costa Rica (CRTM). The study area is delimited by a hydrographic basin labeled <italic>Jacó</italic>, a name given to the study area under analysis, with an extension of 47 km<sup>2</sup>; it consists of three micro-basins that belong to the Pueblo Nuevo (Mona), Copey and Naranjal rivers. An aquifer is located in the lower and middle part of these micro-basins, settled in unconsolidated deposits, with an approximate extension area of 15 km<sup>2</sup>. The remaining area is mainly composed of basalts and sedimentary rocks, where a fractured aquifer is found. These rocks are located in the mountains (<xref ref-type="fig" rid="f10">Figure 10</xref>).</p>
				<p>
					<fig id="f9">
						<label>Figure 9</label>
						<caption>
							<title>Study area location map, Jacó, Puntarenas, Costa Rica. </title>
						</caption>
						<graphic xlink:href="2215-2563-rge-63-123-gf9.png"/>
					</fig>
				</p>
				<p>
					<fig id="f10">
						<label>Figure 10</label>
						<caption>
							<title>Study area geological map, Jacó, Puntarenas, Costa Rica. Modified from <xref ref-type="bibr" rid="B5">Arias (2003</xref>). Redefinition of the Tulín Formation (Maastrichtian-Lower Eocene) of the Central Pacific of Costa Rica.</title>
						</caption>
						<graphic xlink:href="2215-2563-rge-63-123-gf10.png"/>
					</fig>
				</p>
			</sec>
			<sec>
				<title>General characteristics</title>
				<p>The dry season extends for five months, from December to April, January being the driest, while the rainy season includes the months of May to November and the rainiest month corresponds to October. According to the published 2014 database of the <italic>National Meteorological Institute of Costa Rica</italic>, between 1997 and 2013, the Jacó metereological station registered an average maximum rainfall of 616.6 mm in the months of October; the minimum precipitation occurred in the month of January, with an average of 25 mm.</p>
			</sec>
			<sec>
				<title>Geological aspects</title>
				<p>The study area is comprised of three types of lithostratigraphic units that are clearly distinguishable from each other: a) basalts from the Nicoya Complex, including radiolarite, limestone, sandstone and lulite rocks, b) basalts associated with the Tulín Formation, and c) unconsolidated deposits (<xref ref-type="fig" rid="f10">Figure 10</xref>).</p>
				<p><bold>
 <italic>Nicoya Complex</italic>.</bold> Dengo (1962) defines the Nicoya Complex as a series of intensely folded rocks, consisting of basalts, sedimentary and intrusive units. <xref ref-type="bibr" rid="B19">Kuypers (1980, cited in Sprechmann et al., 1984)</xref> indicates that ambiguity in its definition exists because the series of overlying rocks are locally folded and also due to the presence of basaltic extrusions. Therefore, different lithologies are included within the unit and, consequently, the attributed age varies. Based on the aforementioned arguments, <xref ref-type="bibr" rid="B19">Kuypers (1980, cited in Sprechmann et al., 1984)</xref> and <xref ref-type="bibr" rid="B9">Baumgartner (1984</xref>) propose a stratigraphic subdivision for the Nicoya Complex, which consists of two units, Matapalo and Esperanza. <xref ref-type="bibr" rid="B5">Arias (2003</xref>) refers to the basalts that have been associated with the Nicoya Complex as basalts of the Caribbean Oceanic Plateau (COP); these basalts are the oldest rocks.</p>
				<p>The Matapalo Unit is 150 m thick; it is characterized by being formed by oceanic tholeiitic basalts covered by radiolarites, with a minimum thickness of 40 m. The assigned age spans from Callovian to the Aptian (<xref ref-type="bibr" rid="B19">Kuypers, 1980, cited in Sprechmann et al., 1984</xref>). Its contact with the overlying Esperanza unit is not clearly visible (<xref ref-type="bibr" rid="B19">Kuypers, 1980, cited in Sprechmann et al., 1984</xref>). In the Esperanza unit, consisting of tholeiitic basalts and ophiolitic diabases, gabbro rocks, plagiogranite stocks and some radiolarite, intercalations up to 3 m in the upper part of the unit are common. The assigned age is Santoniense (<xref ref-type="bibr" rid="B19">Kuypers, 1980, cited in Sprechmann et al., 1984</xref>).</p>
				<p>In the study area, the geology in mountains is composed of basalts and sedimentary rocks of the types: sandstone, lutite/shale, radiolarites and limestone, that are associated with the rocks of the Nicoya Complex (<xref ref-type="bibr" rid="B6">Arias and Morera, 2002</xref>). These lavas, of basaltic composition, are of a black coloration and, by their origin, have a texture of pillowy folds, which characterizes them. On the surface, they are weathered and easily confused with coarse sands.</p>
				<p>The radiolarites of the Nicoya Complex are part of the Matapalo Unit, from the Upper Jurassic-Lower Cretaceous. They are rocks of pelagic origin, from fine to very fine grain, stratified, quite hard and versicolor (<xref ref-type="bibr" rid="B2">Agudelo, 2005</xref>). They are related to flints and jaspers and are very often associated with iron-manganese mineralization, the result of hydrothermal alteration in ocean beds from exhalative volcanism, either in disseminated form or in nodules. The radiolarites overlap disconformingly the basalts and is crossed by numerous intrusive gabbroid to diabasic bodies. It presents a strong tectonic deformation (isoclinal-chevron folds) and a high degree of faulting (<xref ref-type="bibr" rid="B2">Agudelo, 2005</xref>).</p>
				<p><bold>
 <italic>Tulín Formation</italic>.</bold><xref ref-type="bibr" rid="B5">Arias (2003</xref>) suggests that the Tulín Formation, (Maastrichtian-Lower Eocene), is composed of 95% basalts, which are presented as pillow castings; these basalts have interspersed epiclassic sediments, pyroclastic, gabbros, and olivine cumulites. It is a younger formation than COP, whose rocks constitute the basement of the Tulín Formation.</p>
				<p><bold>
 <italic>Unconsolidated deposits</italic>.</bold> All recent deposits associated with gravity deposits at the foot of the mountain or colluviums and those related to alluvial processes belong to this unit. Also, those deposits of sandy materials with marine influence. Deposits associated to materials dragged by fluvial action or gravity action, are generally fragments of rocks of diverse sizes and composition immersed in a silt-clay sandy matrix. Towards the foot of the mountains, the blocks of fragments are larger and sub-rounded to rounded, while toward the lower part of the basin the fragments are smaller, since they have been dragged a greater distance, as evidence by their lower degree of angularity (<xref ref-type="bibr" rid="B6">Arias and Morera, 2002</xref>).</p>
			</sec>
			<sec>
				<title>Theoretical-conceptual framework</title>
				<p>The estimation of the GOD vulnerability index (<xref ref-type="bibr" rid="B15">Foster and Hirata, 1988</xref>) involves a series of specific stages. First, identify the degree of hydraulic confinement of the aquifer and assign an index to this parameter on a scale from 0.0 to 1.0. Second, specify the characteristics of the substrate overlying the saturated zone of the aquifer in terms of: a) degree of consolidation (taking into account the probable presence or absence of permeability by fissures) and, b) type of lithology (indirectly considering effective porosity, matrix permeability and moisture content in the unsaturated zone or specific retention) and, assigning to this parameter an index on a scale from 0.4 to 1.0. Third, estimate the distance or depth to the water level (in unconfined aquifers) or depth to the ceiling of the first confined aquifer and subsequently assigning an index on a scale from 0.6 to 1.0. In this regard, the danger of groundwater contamination is later defined as the probability of groundwater, in the upper part of an aquifer, becoming contaminated to an unacceptable degree relative to water quality standards as imposed by the activities developed on the surface of the overlying terrain (<xref ref-type="bibr" rid="B15">Foster and Hirata, 1988</xref>, <xref ref-type="bibr" rid="B17">Foster et al., 2007</xref>).</p>
				<p>Thus, instead of applying universal controls over the potentially contaminating activities taking place on the territory and effluent discharges, it is more effective (and less harmful to economic development) to modify the type and level of control in accordance with the attenuation capacity of the aquifer. This is the basic premise behind the concept of contamination vulnerability in aquifers and what leads to the necessity of its mapping (<xref ref-type="bibr" rid="B17">Foster et al., 2007</xref>).</p>
			</sec>
			<sec sec-type="methods">
				<title>Methodological framework</title>
				<p>Geological profiles. In the study area there are 57 private wells registered and four public supply wells belonging to <italic>AyA</italic>. Of 57 private wells, only 21 include a lithological description in their drilling reports. These reports were retrieved in October 2016 from the institutional archive of the Directorate of Water Research and Management of the <italic>National Groundwater, Irrigation and Drainage Service</italic> (<italic>SENARA</italic>) and the Water Directorate of the <italic>Ministry of Environment and Energy</italic> (<italic>MINAE</italic>). Based on this information, at least 11 geological profiles were drawn, from which the four most representative profiles were chosen to determine the relationship between the geological units. A profile was drawn with direction from southwest to northeast, the second profile from northwest to southeast, crossing the lower basin sector, and the third from northwest-southeast, crossing the sector of the upper basin.</p>
				<p>Differential gauging. To define the hydraulic relationship between surface and groundwater, two differential gauging campaigns were planned, one during the dry season and the other in the rainy season, in the main rivers and streams of the basin. Differential gauging has been widely used in various climatic situations to perform flow estimates, for example, in <xref ref-type="bibr" rid="B4">Andreo et al. (2004</xref>) and <xref ref-type="bibr" rid="B7">Barberá and Andreo (2015</xref>). For the selection of gauging sites, consideration was taken not to include those influenced by inflows or outflows of water from other water bodies, and if so, such water flows were quantified to assess the hydraulic relationship between the water bodies (rivers and streams) and the aquifer. In November 2013, during the rainy season, gauging was carried out in 18 sites; in March 2014, dry season, in 12 sites. The difference between the two campaigns resides in that six of points gauged during the rainy season were dry at base flow. Once the flow at each gauging point was calculated, the influent and effluent sections of the river were defined.</p>
				<p>Vulnerability map. The GOD method (Foster et al., 2002) was used to characterize the vulnerability of the aquifers in the study area. To obtain the vulnerability values in each of the layers comprising the GOD method, a review was made of the Geological Map of Costa Rica (<xref ref-type="bibr" rid="B10">Denyer and Alvarado, 2007</xref>), Tectonic Atlas of Costa Rica (<xref ref-type="bibr" rid="B11">Denyer et al. al., 2003</xref>), and a study carried out by <xref ref-type="bibr" rid="B5">Arias (2003</xref>) as well as available digital cartography and orthophotography (<xref ref-type="bibr" rid="B22">PRCR, 2007</xref>) data. In addition, field corroboration of the area’s geology was made, correlating the observed outcrops with the descriptions made by the mentioned authors. Also, lithological information of the wells present in the area was compiled from the drilling reports filed in <italic>SENARA</italic> (<italic>National Service of Groundwater Irrigation and Drainage</italic>) and in the Water Directorate of the <italic>Ministry of Environment and Energy (MINAE)</italic>, as well as drilling/perforation data from <italic>AyA</italic> wells, to estimate the thickness of the unsaturated zone and the depth of the aquifer roofs. Finally, the resultant geological profiles elaborated in this investigation were used to improve the understanding of the conceptual hydrogeological model. A Geographic Information System (GIS) was used to integrate the information required to estimate the values in each layer of the GOD method. Once the vulnerabilities in the study area were defined, a matrix was elaborated with the suggested activities for each zone in accordance with its particular vulnerability and consideration of the activities developed in the region.</p>
			</sec>
			<sec sec-type="results">
				<title>Results</title>
				<sec>
					<title>-Hydrometry</title>
					<p>During flow measurement campaigns, the absence of water intakes for irrigation, reservoirs or any other water use along the the stream was corroborated in the field. In this way, the influent-effluent character of the studied streams was confirmed. In the rainy season, the Mona Stream had a section in the central sector of its route where the flow decreased from 79.9 to 14.7 l/s. This flow reduction indicates that in this sector the stream behaved as an influent and the section that continues to its outlet as an effluent (<xref ref-type="fig" rid="f11">Figure 11</xref>). The Copey River, during the rainy season, behaved as an effluent for most of its course, except for the middle sector, where the flow decreased from 133.2 to 112.9 l/s. The Seca Stream, also during the rainy season, had an effluent behavior, except in its mouth, where its flow changed from 99.6 to 52.4 l/s, becoming an influent stream in that site (<xref ref-type="fig" rid="f11">Figure 11</xref>). During the rainy season, the influent-effluent behavior of the rivers along its course showed that the aquifer has a depth at the same level as the groundwater above the river bed thus the aquifer is capable of feeding rivers and streams, while in other sectors the position of the aquifer water level is below the riverbed, reversing the behavior (i.e. rivers and streams feed aquifer). During the dry season, the behavior of the rivers remained the same, varying only in flow. Likewise, the sections in which the flow decreased drastically during the rainy season are the same ones that were completely dry in the dry season due to the decrease in the aquifer water level (<xref ref-type="fig" rid="f11">Figure 11</xref>), which corroborates the influent character of the rivers and streams in these sectors.</p>
					<p>
						<fig id="f11">
							<label>Figure 11</label>
							<caption>
								<title>Mapping of differential gauging sites in the study area, Jacó Puntarenas, Costa Rica.</title>
							</caption>
							<graphic xlink:href="2215-2563-rge-63-123-gf11.png"/>
						</fig>
					</p>
				</sec>
				<sec>
					<title>- Hydrogeological model</title>
					<p>In the Jacó area, two coastal aquifers have developed; one superficial, detrital, free, settled in the unconsolidated deposit, and a second fractured aquifer settled in the rocks of the Nicoya Complex (<xref ref-type="fig" rid="f12">Figure 12</xref>). This fractured aquifer is free and deep in the mountains and, when it underlies the detrital aquifer in the alluvial plain, it can be confined in at least one sector of the central area of the plain, where the pressure reaches up to 1 m above ground level (drilling record of HE-134 and well 5-AyA) (profile 4, <xref ref-type="fig" rid="f14">Figure 14</xref>). Groundwater in mountainous areas travels through basalts and sedimentary rocks; then, probably, a part is integrated into the groundwater of the unconsolidated deposit forming a single aquifer (profiles 1 and 2, <xref ref-type="fig" rid="f13">Figures 13</xref> and 14). Even if this aquifer was made up of two geological formations, unconsolidated materials and basalts below them, with different hydraulic conductivity and permeability values, it is likely for groundwater to maintain the same level due to the absence of an impermeable layer separating both formations, which is deduced from the drilling reports obtained from <italic>SENARA</italic>. The HE-134 well report indicated that when the well reached the fractured basalts, at 22 m depth, an upwelling occurred (profile 4, <xref ref-type="fig" rid="f14">Figure 14</xref>). The declared upwelling in this, and in the well 5-AyA report, could be due to: a) the difference in permeability between the fractured and detrital formations, b) higher water pressure within the fractured material due to the effect of some clay lenses near the surface (profile 4, <xref ref-type="fig" rid="f14">Figure 14</xref>) and/or, c) because the wells that reached the basalts are exposed to a greater hydraulic potential, for they reached depths with different equipotential lines to those wells that extract water from the detrital formation. In addition, profiles 1 and 2 (<xref ref-type="fig" rid="f13">Figure 13</xref>) show that the topography of the basalts is very irregular and the groundwater level in these basalts is affected by steep slopes.</p>
					<p>
						<fig id="f12">
							<label>Figure 12</label>
							<caption>
								<title>Study area hydrogeological map, Jacó, Puntarenas, Costa Rica.</title>
							</caption>
							<graphic xlink:href="2215-2563-rge-63-123-gf12.jpg"/>
						</fig>
					</p>
					<p>
						<fig id="f13">
							<label>Figure 13</label>
							<caption>
								<title>Hydrogeological profiles 1 and 3 located in the study area, Jacó, Puntarenas, Costa Rica.</title>
							</caption>
							<graphic xlink:href="2215-2563-rge-63-123-gf13.jpg"/>
						</fig>
					</p>
					<p>
						<fig id="f14">
							<label>Figure 14</label>
							<caption>
								<title>Hydrogeological profiles 2 and 4 in the study area, Jacó, Puntarenas, Costa Rica.</title>
							</caption>
							<graphic xlink:href="2215-2563-rge-63-123-gf14.jpg"/>
						</fig>
					</p>
				</sec>
				<sec>
					<title>- Jacó detrital aquifer</title>
					<p>It is a free, detrital aquifer settled in unconsolidated materials such as sands, silts and gravels with a thin layer of soil. For those sectors where water level information is available, it reaches from 0.5 to 5 m deep, with equipotentials between 0.5 and 15 m.a.s.l (<xref ref-type="fig" rid="f12">Figure 12</xref>). The hydraulic direction of this aquifer is to the southwest (SW) and the gradient is 0.01. Another mechanism of recharge is via infiltration from river beds, as demonstrated by differential gauging assessment (<xref ref-type="fig" rid="f11">Figure 11</xref>). </p>
					<p>The aquifer discharge occurs naturally into rivers, and in some sectors (<xref ref-type="fig" rid="f11">Figure 11</xref>), into the sea (<xref ref-type="fig" rid="f12">Figure 12</xref>), as well as artificially by the water extraction from wells. In the alluvial plain no springs were found during field visits, nor are there records of their existence. Well exploitation flows indicate a variable productivity from 1 to 50 l/s; those reporting the greatest extraction flows are wells that reached the fractured basalts.</p>
				</sec>
				<sec>
					<title>-Fractured aquifer</title>
					<p>This fractured aquifer, with a secondary permeability caused by fissures, has developed in lavas and the sedimentary rocks mentioned above, which compose the mountains surrounding the basin. To date, the fractured aquifer registered 19 wells, of which only 15 declared lithological data in their corresponding reports. The depth of groundwater level is variable because the wells themselves show variable depths and the resulting water level corresponds to an equilibrium level between the different producing layers. For example, well HE-109 has a depth of 105 m and a static level reported at 70 m depth, while well HE-58 has a depth of 180 m and a static level at 100 m depth. In the reports of wells HE-18, HE-19, HE-47, HE-68 and HE-134, all located in the central zone, upwelling was declared between 22 and 38 m, where the fractured basalts are located (profile 4, <xref ref-type="fig" rid="f14">Figure 14</xref>).</p>
					<p>As for the productivity of the aquifer, the drilling reports obtained from the <italic>SENARA</italic> archives indicate very low flows, between 0.5 to 7 l/s, with very strong declines. In spite of not having performed pumping tests, this condition possibly reflects a low productivity aquifer as a result of low secondary permeability. Even though it was not possible to map equipotentials for this aquifer due to lack of information, it is likely that the flow is directed towards the alluvial zones, since the rocks that contain this aquifer are found in the hills bordering the plain. Subsequently, in the sector underlying the detrital aquifer, the water is integrated with the groundwater in the detrital aquifer, forming a single level and a single aquifer. Then, the flow is directed southwest towards the sea. However, this is an assumption based on topographic and geomorphological characteristics of the area (<xref ref-type="fig" rid="f12">Figure 12</xref>).</p>
				</sec>
				<sec>
					<title>- Hydrogeological vulnerability</title>
					<p>Once the hydrogeological conditions were defined, the vulnerabilities of the porous and fractured aquifers in the hill areas was estimated.</p>
					<p>Degree of hydraulic confinement (G Value)</p>
					<p>The detrital aquifer was assigned a value of 0.9, as it is an aquifer with little coverage, that is, with a shallow vadose zone, as observed in the created profiles. For the fractured aquifer a value of 0.7 was estimated, corresponding to a free aquifer with coverage according to profiles 1, 2 and 3. (<xref ref-type="fig" rid="f15">Figure 15</xref>). </p>
					<p>Occurrence of the overlying substrate (Value O). </p>
					<p>Based on the lithological description in the drilling/perforation reports of the wells, a layer of clay-silt soil was identified in the alluvial aquifer zone, with a variable thickness between 0 and 14 m, thus a value of 0.5 was estimated. In the mountains, the value is 0.6 and corresponds to fractured basaltic rocks (<xref ref-type="fig" rid="f15">Figure 15</xref>).</p>
					<p>Distance to the groundwater level (D Value). </p>
					<p>In the detrital aquifer, the depth of the water level with respect to surface level is between 1 and 5 m, so a value of 0.9 was assigned, while in the free fractured aquifer, in the mountainous areas, water was found between 20 and 50 m depth, corresponding to a value of 0.7, although there are levels in some wells at greater depth (<xref ref-type="fig" rid="f15">Figure 15</xref>).</p>
					<p>
						<fig id="f15">
							<label>Figure 15</label>
							<caption>
								<title>Degree of hydraulic confinement (G), Occurrence of overlying substrate (O) and Distance to groundwater level (D) of aquifers in the study area, Jacó, Puntarenas, Costa Rica.</title>
							</caption>
							<graphic xlink:href="2215-2563-rge-63-123-gf15.png"/>
						</fig>
					</p>
					<p>As a result, a medium degree of hydrogeological vulnerability was obtained in the area where the detrital aquifer is located, with a value of 0.41 (<xref ref-type="fig" rid="f16">Figure 16</xref>). This sector represented 30.5% of the total study area, while the corresponding area of the fractured aquifer was classified with a low degree of hydrogeological vulnerability, with a value of 0.29.</p>
					<p>
						<fig id="f16">
							<label>Figure 16</label>
							<caption>
								<title>Proposed hydrogeological vulnerability of aquifers in the study area, Jacó, Puntarenas, Costa Rica.</title>
							</caption>
							<graphic xlink:href="2215-2563-rge-63-123-gf16.png"/>
						</fig>
					</p>
				</sec>
			</sec>
			<sec sec-type="discussion">
				<title>Discussion</title>
				<p>Although the vulnerability to contamination depends on the natural attenuation capacity of the soil᾿s profile to pollutants, it varies widely with the geological conditions that prevail near land surface. Among these conditions is the depth at which groundwater resides. This had important implications in designating a low vulnerability to the fractured aquifer, located in the mountains around the Jacó area, because the water level was found at depths greater than 50 m. Low vulnerability aquifers are vulnerable to conservative contaminants, when these are regularly discharged or leached for long periods of time (<xref ref-type="bibr" rid="B17">Foster et al., 2007</xref>). According to <xref ref-type="bibr" rid="B16">Foster et al. (2013)</xref>, an aquifer with low vulnerability can be vulnerable to the leaching of nitrates, highly saline liquids and organochlorine compou nds in dense, non-aqueous phase.</p>
				<p>With respect to the aquifer on the alluvial plain, even though groundwater level is shallow, it was classified as moderate towards vulnerability to contamination. This is because the reports of the wells in this area indicate a surface layer of clays several meters thick, and even reveal intercalations of clays with the gravels inside the alluvial pack, which has given this aquifer a certain degree of protection against certain polluting substances. According to <xref ref-type="bibr" rid="B17">Foster et. (2007</xref>), aquifers with moderate vulnerability are those that are vulnerable to some pollutants when they are discharged or leached regularly. Even though there is a small sector within the aquifer where upwelling is reported, well reports do not indicate the depth at which this condition occurred and the groundwater levels pointed out in these reports represent the equilibrium levels between the confined sector and the more superficial detrital aquifer. For this reason, this sector was classified under the same vulnerability class as the rest of the detrital aquifer.</p>
				<p>The inclusion of hydrogeological vulnerability in local planning programs (<italic>Plan Regulador,</italic> in Costa Rican legal jargon) by local governments in Costa Rica, would be a very useful tool for defining adequate activities to be developed within a municipal territory, and would allow the implementation of integrated water resource management processes. The implementation of the proposed vulnerability map must be accompanied with the application of a vulnerability matrix. In the case of Jacó, most of the urban developments, as well as the main productive and socioeconomic activities, are in the sector of moderate vulnerability of the study area. However, the local plan (<italic>Plan Regulador</italic>) for the development of the coastal zone of Jacó, which has not been approved to date, does not yet contemplate the vulnerability of the aquifer and allows a construction area per property equal to 70%, which is not necessarily the appropriate measure for a moderate vulnerability sector, both for reasons of recharge and potentially polluting hydraulic load. Nor is there regulation that governs the size of future properties in non-urbanized areas within the district. In addition, it must be evaluated how many of the current properties meet the current construction area regulation and the required future control measures.</p>
				<p>The elaboration of vulnerability maps in the area evidenced a need for greater government investment to generate hydrogeological data, mainly through the installation of more piezometers around the Nicoya Complex area. This indication is necessary for any method of assessing hydrogeological vulnerability (<xref ref-type="bibr" rid="B16">Foster et al., 2013</xref>). Due to the scarcity of hydrogeological data in the study area, part of the information must be extrapolated from other sectors. This occurred in sectors that are geologically associated with the rocks of the Nicoya Complex, with little data on the depth of the water level and lithology. Fortunately, the largest proportion of zones with little information are located in areas with slopes greater than 40% hence receiving protection already typified in the current Forestry Law and, therefore, are not susceptible to being developed. On the other hand, due to the scarcity of information available in the area, the methodology that is best suited for this study was that of the GOD method. Similarly, it is recommended to utilize, if more information is generated, other methods of greater accuracy, such as DRASTIC, proposed by <xref ref-type="bibr" rid="B3">Aller et al. (1987</xref>), among others. However, DRASTIC method tends to underestimate the vulnerability of fractured aquifers (Foster et al., 2007), and other more complex methods require greater investment in terms of data input generation, such as percolation rate or aquifer recharge rate. In addition, according to Foster et al. (2013), the most complex methods do not necessarily afford the most realistic vulnerability maps.</p>
				<p>In Costa Rica, land use plans for territorial ordering do not take into account the threat and risk of contamination but only incorporate hydrogeological vulnerability mapping as a tool for managing the groundwater resource. The pollution risk map, by integrating the threat of contamination and aquifer vulnerability, would allow the identification of conflict zones and the territory regulation, with the purpose of designing or taking regulatory corrective measures for the protection of groundwater.</p>
			</sec>
			<sec>
				<title>Recommendations</title>
				<p>Taking into account the estimated level of vulnerability, a series of soil management and protection recommendations were generated for the study area. Even with limited availability of lithological information and lack of infiltration testing for the fractured aquifer, it is recommended to pursue several protective actions as precautionary measures until a better understanding exist of the unsaturated zone over the fractured aquifer.</p>
				<p>This is how, for the alluvial plain, which was classified as a moderate hydrogeological vulnerability zone, the installation of a sanitary sewer system having a treatment system is justified due to potential leaching of polluting substances from septic tanks and latrines. It is recommended that this measure be introduced gradually with the support of the local government, <italic>AyA</italic> and the <italic>Cantonal Council for Interinstitutional Coordination</italic> (<italic>CCCI</italic>) of Jacó. Moreover, improved wastewater treatment plants are mandatory for high-impact hotel, tourism and industrial activities, in accordance with the concept of harmony with the environment, aimed at ecotourism and environmental sustainability.</p>
				<p>Furthermore, it is recommended that for the sector of moderate vulnerability (i.e. detrital aquifer), land use status be changed from pasture and cropland to urban if the previously mentioned sewage water disposal indications are met. In addition, it is recommended that the properties possess a minimum area of 300 m<sup>2</sup>, with a maximum construction area equivalent to 50% of the total area to avoid land surface impermeabilization and favor aquifer recharge. If the cultivation zones in the detrital aquifer area are conserved, it is recommended to encourage improved agricultural management practices, such as the rational use of agrochemicals, organic agriculture, crop rotation and controlled cropping systems to decrease soil utilization burden (soil relaxation) as well as reforestation projects. In this zone, landfills, quarries and/or pits are not recommended.</p>
				<p>Likewise, in the sector of low vulnerability, overlying the fractured aquifer, regulations can be more flexible and permissive. As for urban development, we recommend properties with a minimum 250 m<sup>2</sup> and with a maximum construction area of 50%, as long as the terrain has a slope of less than 40%. As in the areas of moderate vulnerability, agriculture and grazing, must follow appropriate soil management practices, and where possible, reforestation and organic agriculture should be encouraged. Tourism and hotel activities should be conceptually focused towards ecotourism. It is recommended that the following infrastructure and activities be allowed: landfills (of suitable design), pits, quarries and industry with improved treatment plants. Land use recommendation for terrains possessing a slope greater than 40% is as forest land as is currently mandated under Costa Rica Forestry Law to among other things this minimizes erosion.</p>
				<p>Finally, it is recommended to implement or continue executing other integral water resource management tools: economic incentives to promote a shift towards activities of less polluting potential in zones with moderate and high vulnerability; groundwater quality monitoring networks; and implementation of social awareness programs derived from joint watershed management agreements between cantons and localities, supported by universities, municipalities, government and private enterprise.</p>
			</sec>
		</body>
	</sub-article>
</article>