nueva página del texto (beta)
Español (pdf)
Artículo en XML
Referencias del artículo
Traducción automática
Enviar artículo por email
Citado por SciELO 
Similares en
SciELO 
Permalink
1. Introducción
El agua meteórica reside en la superficie continental en forma de glaciares, mantos acuíferos, ríos y lagos. Se origina por la descarga de nubes en forma de lluvia, granizo o nieve al interior de los continentes (Mook, 2002; Sharp, 2007). El ~90 % de la precipitación pluvial en la Tierra se produce sobre los océanos, y tan solo un ~10 % se moviliza en masas de aire al interior de los continentes.
Esta fuente de agua genérica no posee una composición isotópica fija, pues se modifica en su distribución a escala global mediante procesos de fraccionamiento isotópico (Craig, 1961) que controlan la composición isotópica en la precipitación de una masa dada, principalmente en la fase de vapor remanente que queda en dicha masa de aire (Poage y Chamberlain, 2001; Mook, 2002; Sharp, 2007). El fraccionamiento isotópico se inicia con la evaporación de agua oceánica (la principal fuente de agua meteórica), de forma que las moléculas de agua evaporada incorporan preferentemente los isótopos más ligeros de O y H (1H y 16O). Por el contrario, las moléculas de agua constituidas por átomos más pesados de O y H (2H y 18O) se fraccionan principalmente hacia la fase líquida (Dansgaard, 1964). Los principales factores de fraccionamiento isotópico que controlan las variaciones en la concentración de isótopos estables de H y O en el agua meteórica son:
La temperatura: cuando la radiación solar incide sobre la superficie oceánica genera evaporación en la superficie de la misma, generando un cambio de fase que propicia fraccionamiento isotópico en las moléculas de agua. El fraccionamiento isotópico producido por las variaciones en la temperatura obedece a que la energía de vibración de las moléculas de agua es mayor a medida que incrementa la temperatura, propiciando la incorporación de las moléculas más pesadas hacia la fase vapor (Mook, 2002).
Latitud: está relacionada con el grado de radiación solar sobre la superficie oceánica. La composición isotópica disminuye con el incremento en la latitud, de forma que a mayores latitudes el agua en la fase de vapor será isotópicamente más ligera que en zonas ecuatoriales. Tal fenómeno es debido esencialmente a la disminución de la temperatura con la latitud, de forma que en zonas ecuatoriales tiende a atenuarse el fraccionamiento entre el isótopo ligero y el pesado de un elemento, pudiéndose incorporar una mayor cantidad relativa del isótopo pesado a la fase vapor (Craig, 1961; Poage y Chamberlain, 2001; Mook, 2002, Sharp, 2007).
Distancia que recorren las nubes al internarse al continente: asociado dinámicamente a los cambios en la temperatura atmosférica y a las barreras topográficas con las que se encuentre configurado el relieve. Las variaciones en la temperatura atmosférica generan el punto de rocío y, consigo, las precipitaciones que dan lugar a cuerpos de agua meteórica en la superficie. Por su parte, las barreras topográficas impiden la descarga homogénea de agua a lo largo de las cordilleras (Craig, 1961; Poage y Chamberlain, 2001; Mook, 2002, Sharp, 2007).
Elevación topográfica: se registran composiciones isotópicas menores a medida que se incrementa la elevación. De acuerdo con Mook (2002), este efecto se relaciona con la temperatura, debido a que la condensación se produce a partir de la disminución de la misma con el aumento de la altitud. Al disminuir la presión con el incremento en la altitud, el sistema requiere un descenso de la temperatura mayor que en el caso de la condensación por presión isobárica para poder alcanzar la presión de vapor de agua saturada (Dansgaard, 1964; Molnar, 2010). El fraccionamiento isotópico debido a variaciones topográficas tiene lugar en las cordilleras montañosas cuando las nubes viajan al interior de los continentes a lo largo de una pendiente topográfica en incremento: una parte de la fase vapor se condensa y precipita a determinada altura al llegar al punto de rocío. Esto continúa produciéndose a medida que la elevación topográfica se incrementa en su recorrido por la cordillera (Ambach et al., 1968; Rowley y Grazione, 2007; Sharp, 2007; Campani et al., 2012; Poulsen y Jeffery, 2012), en el mismo proceso de fraccionamiento isotópico. De tal forma, las moléculas de agua constituidas por los isótopos ligeros de O y H permanecerán preferentemente en la fase vapor respecto a la fase líquida que precipita a partir de ella. A partir de precipitaciones sucesivas, el agua de lluvia devendrá isotópicamente más ligera con cada precipitación; dicho proceso se conoce como “destilación de tipo Rayleigh”.
Efecto de la orografía: parte de la problemática que diversos autores han encontrado al realizar estudios sobre las variaciones isotópicas en agua meteórica radica en la configuración de barreras en el relieve que obstruyen el libre paso del vapor atmosférico en las cordilleras. Éstas provocan un fraccionamiento isotópico dinámico al canalizar y concentrar el vapor atmosférico, y forzando con ello la descarga de nubes hacia los valles conformados de manera local.
Efecto de pseudoaltitud: consiste en un enriquecimiento en 18O y 2H durante la evaporación de las gotas de lluvia por debajo de la base de las nubes, principalmente en las zonas de poca altitud, tales como valles profundos en zonas montañosas (Moser y Stichler, 1974; Schemmel et al., 2013).
El estudio de la composición isotópica de agua meteórica se ha aplicado en los últimos años con especial énfasis para determinar la paleoaltimetría de algunas cordilleras orogénicas (Chamberlain et al., 1999; Garzione et al., 2000; Rowley 2007; Rowley y Grazione, 2007; Campani et al., 2012). Por ello, la necesidad de entender con mayor detalle los factores que modifican la composición isotópica del agua meteórica que descarga en los continentes ha adquirido una importancia especial. De ahí parte la necesidad de estudios relacionados con la elevación topográfica y la distancia que recorren las nubes a lo largo de la pendiente topográfica hasta encontrar barreras topográficas que puedan promover la generación de precipitaciones. Hasta el momento, sólo existe un estudio sistemático en México enfocado en determinar la línea de agua meteórica local en asociación con un orógeno (Pérez Quezadas et al., 2015). Los estudios existentes sobre la composición isotópica del agua meteórica son de escaso alcance, y están mayoritariamente relacionados con estudios de origen del agua en campos geotérmicos (e.g., Truesdell et al., 1978; Tello, 1982; González-Partida et al., 2005; Birkle et al., 2016).
En el presente estudio se analizan las variaciones en la composición isotópica de agua meteórica en una sección del Cinturón de Pliegues y Cabalgaduras Mexicano (CPCM), a fin de obtener sólo las variaciones debidas al incremento en la elevación topográfica. Se trata de establecer la línea del agua meteórica local (LAML) a partir del análisis de las concentraciones isotópicas de agua de lluvia, con la finalidad de fijarla como grupo de valores de referencia que permita discriminar entre fuentes de agua distintas en cualquier estudio isotópico de O y H. La LAML permitirá identificar la presencia de agua meteórica en estudios sobre geoquímica de fluidos en yacimientos minerales, geología estructural, hidrogeología, etc. La sección del CPCM estudiada se ubica geográficamente desde el centro del país hacia el noreste, hasta las inmediaciones de la Planicie Costera del Golfo de México (Figura 1). Por sus características morfológicas y ubicación geográfica, el CPCM en su parte central constituye un laboratorio natural para este tipo de estudio.
2. Características de la región estudiada
2.1. Geología
El Cinturón de Pliegues y Cabalgaduras Mexicano (CPCM) es un rasgo orogénico en forma de cuña que se emplazó por deformación de piel delgada thin-skinned, producto de la exhumación y acortamiento de cuatro elementos paleogeográficos (Suter, 1980, 1984, 1987; Eguiluz de Antuñano et al., 2000; Fitz-Díaz, 2010; Fitz-Díaz et al., 2011). Tales elementos son dos plataformas marinas (El Doctor y Valles-San Luis Potosí) y dos cuencas sedimentarias (Zimapán y Tampico-Misantla; Figuras 2 y 3). Éstos corresponden a la configuración predeposicional del basamento en pilares y fosas, respectivamente, generadas durante el proceso de apertura del Golfo de México. En las rocas de la Plataforma El Doctor se desarrollaron principalmente cabalgaduras. La parte oriental de dicha plataforma cabalga sobre rocas de la Cuenca de Zimapán, donde la deformación propició el emplazamiento de anticlinales y sinclinales asimétricos con sus flancos frontales invertidos, inclinados hacia el suroeste (Carrillo, 1971; Carrillo-Martínez, 1990). El borde oriental de la Cuenca de Zimapán cabalga, a su vez, sobre la Plataforma Valles-San Luis Potosí, la cual muestra deformación similar a la de la Plataforma El Doctor. Asimismo, la Plataforma Valles-San Luis Potosí cabalga en su borde oriental sobre la Cuenca Tampico-Misantla, que presenta desarrollo de pliegues transportados en dirección ENE. Es decir, las cuencas y plataformas se invirtieron, pasando a ser elementos topográficos positivos, en forma de escamas tectónicas imbricadas hacia el oriente, y registran una propagación en la intensidad de la deformación con esfuerzos orientados de SW a NE (Campa-Uranga, 1983; Davis et al., 1983; Eguiluz de Antuñano et al., 2000; Fitz-Díaz et al., 2010, 2011, 2014; Ortega-Flores, 2011; ver Figura 3).
Figura 2 Mapa geológico de la región estudiada. La sección estudiada va de A hacia A’, en círculos rojos se indican las localidades de muestreo. Modificado de Fitz-Díaz et al. (2011).
Figura 3 Perfil geológico de la sección estudiada, que muestra la disposición espacial de los elementos paleogeográficos que la componen, la geometría de la cuña orogénica en función de los mismos y una representación del relieve con las elevaciones topográficas máximas y mínimas. Modificado de Fitz-Díaz et al. (2011).
2.2. Geomorfología y clima
La sección de muestreo se ubica en el centro del país y se extiende hacia el noreste del país atravesando los estados de Querétaro, Hidalgo y San Luis Potosí. El conjunto de características morfoclimáticas que la caracterizan son:
Se ubica en una latitud prácticamente fija entre 21º y 22ºN, por lo que no se esperan variaciones en la composición isotópicas de agua meteórica debidas al efecto latitudinal.
Presenta un relieve en forma de cuña orogénica con una pendiente continua entre 1º y 1.5º (Fitz-Díaz et al., 2011), con elevaciones máximas de aproximadamente 3100 msnm, descendiendo hasta la Planicie Costera del Golfo de México (Figura 4), y sólo existen algunos valles y depresiones del orden de pocas decenas de metros. Tales características permiten obviar el efecto de pseudoaltitud.
No existen fuentes relevantes adicionales de agua meteórica que puedan modificar significativamente la composición isotópica de agua precipitada por el acarreo de nubes desde el Golfo de México.
No existen barreras topográficas que impidan la libre circulación de los vientos alisios, lo cuales son responsables del transporte de nubes que descargan al interior del continente. Ello permitió la captación del agua meteórica derivada de cada destilación de tipo Rayleigh producida a consecuencia de la continentalización de las masas de vapor de agua.
Los vientos que transportan el vapor atmosférico producto de la evaporación de agua oceánica proveniente del Golfo de México son los alisios, los cuales son regulares durante todo el año.
Figura 4 Mapa digital de elevaciones de la región estudiada: los tonos cálidos corresponden a las elevaciones mayores y los fríos a las de menor cota topográfica. La localización de los colectores de agua de lluvia instalados se indica con triángulos azules, los de escorrentía se indican en rectángulos ocres.
Dichas características deben de permitir establecer una LAML) perfectamente representativa de su composición isotópica, en ausencia de perturbaciones debidas a (a) factores de fraccionamiento distintos al que es debido solamente a la altitud, ni (b) fuentes de agua distintas a las masas de vapor derivadas de evaporación marina en el Golfo de México.
A medida que la elevación topográfica se incrementa, la temperatura de las masas de aire reportadas y las precipitaciones disminuyen. De acuerdo con el Servicio Meteorológico Nacional (2012), existen variaciones climáticas contrastantes relacionadas con la elevación topográfica, teniendo un clima cálido sub-húmedo en la Planicie Costera del Golfo con temperaturas promedio anuales entre 22° y 26 °C y una precipitación promedio anual de 974.4 mm. En el estado de Hidalgo, al centro de la sección, el clima es templado húmedo, con temperaturas de 18° a 22 °C y una precipitación promedio anual de 652.4 mm. En el centro del país, en las inmediaciones de los estados de Querétaro e Hidalgo, el clima es templado sub-húmedo y presenta temperaturas de 18° a 22 °C y una precipitación promedio anual de 589.9 mm.
3. Metodología
La metodología llevada a cabo para cumplir los propósitos del trabajo se dividió en 3 etapas: análisis de relieve, muestreo de campo y análisis de laboratorio, que se detallan a continuación.
3.1. Análisis de relieve
Se desarrolló un análisis morfométrico del relieve con el software ArcGis versión 9.2, utilizando la topografía de 8 cartas a escala 1:50000 de INEGI en formato shapefile: F14C48, F14C49, F14C58, F14C59, F14D31, F14D32, F14D41 y F14D22. Se obtuvo, en primer lugar, un MDE (Modelo
Digital de Elevación) que se reclasificó con la finalidad de obtener un mapa de elevaciones a la misma escala (Figura 4). El mapa de elevaciones obtenido a partir de un MDE muestra en colores cálidos las elevaciones topográficas más elevadas de la región, cambiando a tonos verdosos conforme disminuye la cota. Posteriormente, se determinó la sección de estudio en función del continuo decremento en la cota topográfica que se obtuvo en el mapa, eligiendo localidades para el muestreo en las zonas topográficamente más altas que representaran la pendiente de la cuña orogénica. Se aprecia que las localidades donde se colocaron gradualmente disminuyen en elevación, siguiendo un relieve que decrece hacia el noreste (Figura 4). Con ello se evitó instalar colectores de agua pluvial en depresiones topográficas mostradas en el análisis de elevaciones; dichas depresiones representan valles o cuerpos de agua importantes, como lagos o lagunas. Igualmente, se pudieron evitar zonas con microclimas o evapotranspiración de zonas boscosas, favoreciendo con ello la captación específicamente de agua de lluvia no afectada por procesos secundarios de evaporación. Finalmente, se eligieron 12 localidades de muestreo, con elevación máxima de 2740 msnm y mínima de 56 msnm (ver Tabla 1).
Tabla 1 Relación de puntos de muestreo de agua, su localización geográfica de los mismos, y composición isotópica en H y O del agua respecto a VSMOW.
| Muestra | Localidad | Altura (msnm) | X (m) | Y (m) | δ2H(‰) | δ18O(‰) |
|---|---|---|---|---|---|---|
| Agua de lluvia: | ||||||
| TRA-12 | El Doctor | 2743 | 521439 | 2351174 | -80.34 | -11.75 |
| TRA-11 | San Joaquín | 2430 | 511266 | 2346815 | -68.87 | -10.28 |
| TRA-10 | Apartadero | 1640 | 510007 | 2339356 | -60.06 | -9.05 |
| TRA-9 | Jiliapan | 1525 | 501441 | 2335912 | -52.42 | -8.33 |
| TRA-8 | Pacula | 1313 | 492144 | 2331613 | -57.18 | -8.84 |
| TRA-7 | Cerro Prieto | 1589 | 469229 | 2327523 | -48.15 | -8.01 |
| TRA-6 | Santa Ana de Allende | 1387 | 464454 | 2318887 | -42.09 | -7.03 |
| TRA-5 | Chapulhuacán | 961 | 449449 | 2312756 | -42.83 | -7.15 |
| TRA-4 | Taman (Los Amigos) | 390 | 441149 | 2312756 | -32.55 | -5.91 |
| TRA-3 | Tamazunchale | 153 | 438715 | 2305551 | -36.72 | -6.38 |
| TRA-2 | San Martín Chalchicuautla | 201 | 425475 | 2294723 | -30.76 | -5.63 |
| TRA-1 | Tanquián | 56 | 401437 | 2292969 | -25.35 | -4.89 |
| Escorrentía: | ||||||
| Agua de río | ||||||
| E-1 | San Martín Chalchicuautla | 282 | 506912 | 2370178 | -31 | -5.71 |
| E-2 | Xilitla | 694 | 500914 | 2364554 | -35.8 | -6.28 |
| E-3 | Xilitla | 696 | 500914 | 2364554 | -35.6 | -6.25 |
| E-4 | Pacula | 1818 | 466176 | 2321347 | -53.3 | -8.04 |
| E-5 | Apartadero | 1865 | 446306 | 2314238 | -64.5 | -9.79 |
| Agua de manantiales y cuevas | ||||||
| E-6 | Apartadero | 1885 | 446306 | 2314238 | -63.3 | -9.66 |
| E-7 | Xilitla | 529 | 501012 | 2365438 | -34.5 | -6.19 |
Nota: la altura y coordenadas utilizadas fueron las registradas por medio de GPS, con datum WGS84, zona 14.
3.2. Muestreo de campo
Para realizar el muestreo de agua meteórica se construyeron 12 colectores de agua pluvial de tal forma que garantizaran el adecuado almacenaje del agua captada, teniendo particular cuidado en evitar la evaporación dentro del colector y el escape de vapor a la atmósfera. Para ello se utilizaron bidones de 20 litros de capacidad, forrados con materiales reflectantes (hacia el exterior) y aislantes térmicos (hacia el interior), como hule bajo-alfombra. Se colocó una cubierta reflectante de aluminio como recubrimiento final, con el propósito de evitar la radiación solar incidente. En la parte superior del bidón se colocó un embudo de 15 centímetros de diámetro para captar el agua pluvial, cubierto por una malla de alambre que sirvió como filtro de partículas con diámetros mayores a 2 mm. Por último se agregó una capa de aproximadamente 3 mm de aceite mineral en el interior del bidón, a modo de tapón para evitar el escape de vapor que pudiera ser generado dentro del bidón.
Los colectores se instalaron en las localidades seleccionadas, en azoteas con espacios libres de cualquier escurrimiento y de acceso restringido (preferencialmente en edificios como escuelas o presidencias municipales) y sobre bases pintadas de blanco para evitar al máximo efectos indeseables debidos al calentamiento del piso por radiación solar. La época seleccionada para colectar el agua pluvial fue la temporada de lluvias del año 2012, entre los meses de julio y octubre, En el mismo periodo se tomaron muestras de agua de manantiales y ríos de flujo constante en 7 localidades ubicadas en las inmediaciones de los puntos de muestreo donde se instalaron los colectores, con la finalidad de obtener composiciones isotópicas de control en el sistema isotópico de la región (Figura 4). Durante la temporada de lluvia los colectores fueron constantemente monitoreados para evitar que se derramaran o contaminaran.
3.3. Análisis isotópicos
Una vez concluida la temporada de lluvia se recuperaron intactos los colectores de los sitios donde se instalaron, pudiendo recuperar el 100 % de muestra de agua en ellos, en condiciones óptimas para su análisis en laboratorio. El análisis de la composición isotópica de H y O se efectuó en el Laboratorio de Isótopos Estables del Instituto de Geología de la UNAM. En primer lugar, se separaron las fases de aceite y sólidas de las muestras de agua mediante embudos de separación y papel filtro. Las muestras que contenían materia orgánica se agitaron con carbón activado para adsorberla y eliminarla de la fase líquida. Una vez obtenidos aproximadamente 100 mL de muestra de agua de cada bidón, se filtraron para eliminar cualquier residuo sólido. Para las muestras tomadas de escorrentía se siguió el mismo procedimiento.
Los análisis de δ2H se efectuaron mediante un espectrómetro de absorción láser Los Gatos Research modelo 908-0008-3001, y para los análisis de δ18O se utilizó un espectrómetro de masas Thermo MAT 253 con una interfaz acoplada Thermo Finnigan de GASBENCH II. La calidad de los análisis de cada muestra fue controlada mediante el uso de los estándares de trabajo VSMOW, GISP, SLAP (ver Coplen, 1988) y agua tipo I FLEX 2010 grado HPLC. Los datos se reportaron mediante la notación convencional delta como desviaciones por mil respecto al estándar del promedio del agua oceánica (VSMOW, Vienna Standard Mean Oceanic Water).
4. Resultados
Los resultados de δ18O y δ2H obtenidos de agua pluvial y escorrentía se encuentran contenidos en la Tabla 1 y en las Figuras 5 y 6. La relación entre las composiciones isotópicas de las muestras de agua de lluvia y la elevación de las localidades estudiadas se muestra en las Figuras 7 y 8 (para δ2H y δ18O, respectivamente). Dicha relación es lineal e inversa respecto al incremento en la cota topográfica. El valor de R2 obtenido de δ2H y δ18O es de 0.90 y 0.91, respectivamente.
Figura 5 Línea de Agua Meteórica Local obtenida por el análisis de 12 colectores de agua pluvial dentro de la zona estudiada.
Figura 6 Comparativa de la Línea de Agua Meteórica Local en negro y Línea de Agua Meteórica de escorrentía en rojo, esta muestra la diferencia entre las ecuaciones de la recta que definen a ambas.
Figura 7 Variaciones en los valores de la composición isotópica de H en el perfil geológico y topográfico de la sección estudiada. Clave: PD = Plataforma El Doctor; CZ = Cuenca de Zimapán; PV-SLP = Plataforma Valles-San Luis Potosí; CTM = Cuenca Tampico-Misantla.
Figura 8 Variaciones en los valores de la composición isotópica de O en el perfil geológico y topográfico de la sección estudiada. Clave: PD = Plataforma El Doctor; CZ = Cuenca de Zimapán; PV-SLP = Plataforma Valles-San Luis Potosí; CTM = Cuenca Tampico-Misantla.
Los valores de δ2H obtenidos para muestras de agua de lluvia fluctúan entre -80 y -25 ‰ y los de δ18O varían entre -12 y -4 ‰. Los valores de δ2H para las muestras de agua de escorrentía varían entre -60 y -36 ‰ y los de δ18O varían entre -10 y -6 ‰. Estos valores, tanto de agua pluvial como de agua de escorrentía en la sección estudiada son progresivamente menores a medida que se incrementa la elevación topográfica (Tabla 1, Figuras 9 y 10). Cabe hacer notar que el agua de escorrentía presenta composiciones isotópicas similares a la del agua de lluvia en muestras obtenidas a elevaciones topográficas similares (Figura 6).
Figura 9 Correlación entre la composición isotópica de H obtenida en el análisis del agua de lluvia de la sección estudiada respecto a las variaciones de altitud.
Figura 10 Correlación entre la composición isotópica de H obtenida en el análisis del agua de lluvia de la sección estudiada respecto a las variaciones de altitud.
Con los resultados obtenidos de las concentraciones isotópicas en agua de lluvia se definió la ecuación de regresión lineal que representa la LAML, que define la variación de δ2H y δ18O (Figura 5) del agua meteórica a lo largo de la sección estudiada en el Cinturón de Pliegues y Cabalgaduras Mexicano. La ecuación de la recta que la define es:
Ésta es similar a la Línea de Agua Meteórica Mundial (Craig, 1961):
La diferencia de la ordenada al origen que exhibe la LAML respecto a la línea de agua meteórica correspondiente al promedio mundial definido por Craig (1961) es de +5; el valor de R2 obtenido en la LAML es de 0.99, mientras que la variación en la pendiente es de 0.15.
En comparación, la Línea de Agua Meteórica que se obtiene del análisis de agua de escorrentía se define por la ecuación:
La diferencia de la ordenada al origen que exhibe dicha correlación lineal respecto a la LAML definida en el presente estudio es de +1; el valor de R2 es de 0.99, mientras que la variación en la pendiente es de 0.17 respecto a la misma (Figura 6).
La tasa de fraccionamiento isotópico calculada por mínimos cuadrados por cada kilómetro de elevación obtenido a partir de los valores de muestras de agua de lluvia es de -2.19 ‰ km-1 para δ18O y de -17.75 ‰ km-1 para δ2H (Figuras 9 y 10).
5. Discusión
La sección elegida para el análisis de la variación en la composición isotópica de agua meteórica en este trabajo forma parte de una porción central del Cinturón de Pliegues y Cabalgaduras Mexicano, el cual tiene un rasgo morfológico de cuña orogénica que decrece en elevación topográfica de manera regular al noreste del país (Figura 3). La zona de estudio presenta una pendiente prácticamente constante, con diferencias de latitud muy escasas y en ausencia de barreras orográficas relevantes, a diferencia de otros sitios en el mundo donde se han realizado estudios con propósitos similares, tales como la Sierra Nevada en la Alta California (Ingraham y Taylor, 1991), los Alpes (Campani et al., 2012), los Táurides en Turquía (Schemmel et al., 2013), o una sección que parte del Citlaltépetl hasta la planicie del Golfo de México (Pérez Quezadas et al., 2015), entre otros. Las ventajas de contar con una sección de estudio como la utilizada en este trabajo son múltiples, pues ello permite (1) calcular la tasa de fraccionamiento en función de la altura con mucha mayor precisión que en los casos antes mencionados, y (2) eliminar los efectos de los diversos factores de fraccionamiento adicionales al que es exclusivamente debido a la elevación topográfica.
Los valores que se obtuvieron durante el análisis de la composición isotópica de agua meteórica de la temporada de lluvias de 2012 muestran valores de composiciones isotópicas progresivamente menores conforme se incrementa la elevación topográfica (Figuras 7 y 8). Dicha tendencia ha sido discutida en diversos estudios, destacando el de Poage y Chamberlain (2001), quienes compilaron datos de 68 regiones del mundo ubicadas en Norteamérica, Centroamérica, Sudamérica, Europa, el Himalaya (hasta los 5000 msnm) y localidades en extremos latitudinales terrestres norte y sur. Poage y Chamberlain (2001) reportan la misma tendencia en el empobrecimiento de isótopos pesados de agua meteórica en todos los casos estudiados, debido a que el vapor atmosférico que se interna al continente es afectado por la destilación tipo Rayleigh y la obstrucción de barreras topográficas.
La relación entre δ2H y δ18O para el agua meteórica es proporcional y fue definida por Craig (1961) como la Línea de Agua Meteórica Mundial (LAMM), en un estudio que integra datos de diversas partes del mundo. En éste se considera la posición latitudinal de algunas regiones como principal factor de empobrecimiento o enriquecimiento de isótopos pesados. Con base en dicha observación y a lo discutido en trabajos posteriores (v.g., Friedman et al., 1964; Mook, 2002; Sharp, 2007) se puede afirmar que, en la sección objeto del presente estudio, la posición latitudinal controla la temperatura que incide sobre la superficie oceánica, siendo éste el primer factor de fraccionamiento isotópico hacia la fase de vapor. Lo anterior se encuentra representado por el valor de la ordenada al origen de la ecuación de la LAML que es +15 en las muestras colectadas con agua de lluvia, respecto al de la LAMM que es +10, así como la ligera desviación en la pendiente, de 8.15, respecto al valor del promedio global, que es de 8. Ello refleja condiciones de mayor evaporación al promedio mundial en la latitud comprendida entre los 21º y 22º norte, representada por la sección orogénica estudiada. El valor de R2 de la LAML es de 0.99, lo que indica un buen ajuste hacia la recta, ocasionado por la adecuada captura in situ, la ausencia de fuentes de vapor atmosférico residual en el muestreo y el correcto proceso de análisis de agua de lluvia. Por otra parte, las muestras analizadas de agua de escorrentía forman parte del mismo sistema isotópico y no parecen tener una fuente secundaria que aporte agua a ríos y manantiales, reflejado en la ecuación lineal que la define (Ecuación 3) la cual muestra una ligera desviación en la pendiente (de 8.32), así como en la ordenada al origen (+16) respecto a la LAML (Figura 6). Esto es un comportamiento esperable originado por el transporte que ha tenido sobre la superficie terrestre y la evaporación que se origina como consecuencia de ello.
En comparación, en un estudio similar en la Sierra Nevada de la Alta California, Ingraham y Taylor (1991) reportaron tasas de fraccionamiento en función del incremento en la cota topográfica para isótopos de O entre -3.1 y -2.1 ‰ km-1, correspondientes respectivamente a latitudes entre 37º y 39º N. En cambio, a latitudes entre 33º y 37º S, Rozanski y Araguás (1995) obtuvieron tasas de fraccionamiento de isótopos de O entre -2.1 y -2.2‰ km-1. Las diferencias en el valor de las tasas de fraccionamiento isotópico reportados en dichos trabajos no son ocasionadas sustancialmente por la posición latitudinal, de lo cual es posible deducir que el empobrecimiento en isótopos pesados en las secciones anteriormente mencionadas es más sensible a cambios en el incremento de la topografía y a la obstrucción propiciada por las barreras orogénicas que a otros factores de fraccionamiento. La pendiente regular del área de estudio propicia la entrada homogénea de vapor atmosférico derivado de la evaporación de agua de mar, por medio de los vientos alisios que entran por el Golfo de México. Como consecuencia de ello, el valor de R2 (de 0.99) de la LAML ajusta a la recta, lo que indica ausencia de fuentes secundarias de vapor que modificarían la composición isotópica de la fuente primordial. Poage y Chamberlain (2001) y Schemmel et al. (2013) discuten la participación de barreras topográficas y diversas masas de aire entrando en una cordillera orogénica como factores principales que modifican la composición isotópica. Schemmel et al. (2013) determinaron tres Líneas de Agua Meteórica Local en distintas secciones estudiadas en las montañas Taurus y Pontiac en Turquía, con diferencias considerables respecto a la LAMM obtenida por Craig (1961):
Las Ecuaciones 4 a 6 (Schemmel et al., 2013) muestran diferencias sustanciales entre sí, tanto en la pendiente como en la ordenada al origen. En las Ecuaciones 4 y 5 no existe una fuente secundaria de vapor que modifique la concentración de isótopos de agua meteórica. Sin embargo, en la Ecuación 6, determinada en el valle formado por ambas cordilleras, existe una notoria aportación de vapor originado por procesos secundarios, lo cual puede identificarse tanto en la pendiente como en la ordenada al origen. Así, en comparación, la similitud entre la LAML obtenida para este trabajo (Ecuación 1) respecto a las ecuaciones 2 (LAMM de Craig, 1961), 4 y 5 (Schemmel et al., 2013) sugiere la ausencia de fuentes secundarias de vapor, lo cual indica que el vapor del que derivaron las precipitaciones analizadas en este estudio es debido exclusivamente a la evaporación de agua marina en el Golfo de México.
Las variaciones en los valores de δ2H y δ18O establecen, en ambos casos, correlaciones lineales con la elevación topográfica, de modo que ésta constituye el factor de fraccionamiento esencial para la sección estudiada. Este resultado era el esperado por el modo en que fue seleccionada la zona de estudio y en que fue diseñado el muestreo. Las tasas de fraccionamiento obtenidas en el presente estudio son debidas a la pendiente continua que permite llevar a cabo la destilación Rayleigh de manera ideal y sin barreras que afecten dicho proceso. La buena correlación ente valores isotópicos y altitud obtenida indica que no existen algunas otras fuentes de vapor residual (e.g., concentrado por algún alto topográfico) las cuales pudieran modificar el valor de la composición isotópica en la región. Es decir, el valor de las tasas de fraccionamiento isotópico para H y O analizadas en este trabajo sólo es debido al incremento de la elevación topográfica.
6. Conclusiones
En este trabajo se determinaron las variaciones en δ2HVSMOW y δ18O VSMOW, y el fraccionamiento isotópico en el agua meteórica a lo largo de una sección del Cinturón de Pliegues y Cabalgaduras Mexicano (CPCM) con características óptimas para ello. Los resultados obtenidos reflejan un descenso en δ2H y δ18O a medida que la elevación topográfica aumenta, desde δ2H = -25.35 ‰ y δ18O = -4.89 ‰ en la cota topográfica más baja (56 msnm), hasta δ2H = -80.34 ‰ y δ18O = -11.75 ‰ en la cota topográfica más alta (2740 msnm).
La ecuación de la recta que define la Línea del Agua Meteórica Local (LAML; δ2H = δ18O * 8.15 + 15) tiene una ordenada al origen de 15, a diferencia del promedio mundial definido en la Línea del Agua Meteórica Mundial, con un origen de 10. Ello es debido al efecto de una mayor evaporación a los 21º N en el Golfo de México que en el promedio global.
Las tasas de fraccionamiento en función de la elevación topográfica que se obtuvieron en la sección estudiada fueron de -2.19 ‰ km-1 para δ18O, y de -17.75 ‰ km-1 para δ2H.
La ecuación de la recta que define la Línea de Agua Meteórica en agua de escorrentía muestra una ligera variación respecto a la LAML, tanto en la pendiente, como en la ordenada al origen, lo cual es consecuencia de la evaporación que ha sufrido durante su recorrido desde la zona de recarga, hacia el punto donde se colectó.
