2.1 Muestras de crudo sometidas a meteorización en suelo

El crudo fue recolectado en el yacimiento del Mosquito (S 51°51'51'' y O 68°57'7''), perteneciente a la cuenca Austral (formación Springhill), provincia de Santa Cruz, Argentina. El mismo se caracterizó por una gravedad API de 32° a 15 °C (crudo liviano) y una salinidad de 52 g.m-3. Se transportó al laboratorio en una botella limpia de vidrio (un litro de capacidad, color ámbar) para evitar la fotooxidación. En el llenado de la botella se evitó la presencia de una cámara de aire, para minimizar el impacto del oxígeno en la estabilidad del crudo y, además, prevenir su biodegradación. Luego fue almacenada en un lugar oscuro y seco a temperatura ambiente hasta la realización de la experiencia dentro de las siguientes 72 h ^{(Tomas y Acuña, 2022}).

2.2 Ensayo de meteorización en suelo

Se tomó una muestra de suelo patagónico de aproximadamente 5 kg en las inmediaciones de la ciudad de Río Gallegos, Santa Cruz, Argentina, a una profundidad comprendida entre los 0 y 15 cm (horizonte A) respecto de la superficie, tras eliminar la hojarasca presente. Dicha muestra se tamizó con una malla de 2 mm de poro y se almacenó en refrigeración hasta el momento de su estudio. El suelo se caracterizó por tener una textura franco-arenosa, y las siguientes propiedades físico-químicas: pH 7.6; cloruro 37.4 ppm; sulfato 48.7 ppm; bicarbonato 45.7 ppm; carbonato < 1 ppm; calcio 40.8 ppm; magnesio 12.2 ppm; nitrito 0.52 ppm; nitrato 18.1 ppm; amonio 0.05 ppm; fosfato < 1 ppm; conductividad 84.5 μS.cm-1; humedad 1.07 %; materia orgánica 0.98 %; materia inorgánica 99.02 %; densidad aparente 1,31 g.cm-3; densidad real 2,46 g.cm-3; porosidad 46 %; hidrocarburos totales < 1 mg.kg-1; bacterias aerobias totales (BAT) 5.30 x 105 UFC.g-1; bacterias degradadoras de hidrocarburo (BDH) 8.90 x 103 UFC.g-1 (^{Acuña y Pucci, 2022}).

El estudio de meteorización se llevó a cabo en 30 reactores a escala de laboratorio (frascos de vidrio transparentes de 300 ml de capacidad) monitoreados por el período de un año y repartidos equitativamente a lo largo de diez etapas. En cada uno de ellos se añadieron 10 g de suelo y de forma progresiva 1 ml de crudo; con la ayuda de una espátula se distribuyó de forma homogénea el mismo, con el fin de garantizar el mayor esparcimiento de los hidrocarburos a lo largo de las partículas del suelo. Los reactores se dejaron abiertos, es decir, sin tapa a la temperatura ambiente promedio del laboratorio de 20 °C, para que fenómenos tales como la fotooxidación, evaporación y biodegradación modificaran la composición química del hidrocarburo. El tiempo cero (T0) se consideró para las muestras de crudo agregadas a los tres reactores iniciales y recuperadas inmediatamente para su análisis, es decir, a los cero días de comenzada la experiencia. El resto de las etapas integradas también por tres reactores se detalla a continuación: 7, 14, 21, 28, 60, 120, 180, 270 y 365 días constituyen los tiempos de análisis T1, T2, T3, T4, T5, T6, T7, T8 y T9, respectivamente.

2.3 Obtención del crudo meteorizado

Cumplidos los intervalos de tiempo de cada etapa, se procedió a extraer el crudo agregado a los reactores mediante el empleo de la extracción acelerada con disolventes (ASE) con el equipo Dionex ASE 150 de Thermo Scientific, con n-pentano como solvente de extracción. Las condiciones de extracción fueron: temperatura de celda de 175 ºC con una presión de 1500 psi y un tiempo de extracción de 5 min en un único ciclo. Los extractos obtenidos se concentraron bajo corriente de nitrógeno hasta un volumen de 10 ml; de cada una de ellas se recogió una alícuota de un 1 ml en viales independientes y se refrigeraron (extractos ASE). Los volúmenes restantes siguieron evaporándose hasta alcanzar los 0.5 ml, luego se sometieron a una separación por cromatografía de adsorción sólido-líquido en columna de vidrio (20 cm x 1.2 cm) empacada con 3 g de gel de sílice (activado a 150 °C durante un período de 24 h en la estufa), 50 mg de sulfato de sodio activado y 50 mg de alúmina activada. Para esta separación se sembraron los 0.5 ml de los extractos del suelo en la columna y posteriormente se eluyeron con 10 ml de n-pentano y 10 ml de diclorometano, para obtener las fracciones de hidrocarburos alifáticos y aromáticos, respectivamente. Estos extractos, a su vez, se concentraron hasta 0.5 ml bajo corriente de nitrógeno, se transfirieron a los viales correspondientes, y se almacenaron a -15 °C hasta su análisis (Tomas y ^{Acuña, 2022}).

2.4 Identificación de n-alcanos, isoprenoides y HAP por GC/MS

El análisis por GC/MS se llevó a cabo teniendo en cuenta las publicaciones de Tomas y ^{Acuña (2022}). Por cada muestra se inyectó en un cromatógrafo de gases, un μl del extracto ASE en modo split y un μl de los extractos alifáticos y aromáticos en modo splitless. La separación se realizó en un equipo Agilent modelo 7890A, con un detector de espectrometría de masas de la misma marca (modelo 5975C). Se utilizó una columna HP5ms de 30 m de largo, con un diámetro interno de 0.32 mm y un espesor de película de 0.25 μm. La temperatura del inyector se ajustó a 290 °C y se utilizó helio como gas acarreador con una velocidad de flujo de 1.2 ml.min-1. El programa de temperatura utilizado fue el siguiente: temperatura inicial de 55 °C durante 2 min, seguido de una rampa de 6 °C.min-1 hasta alcanzar los 270 °C, pasando directamente a otra rampa de 3 °C.min-1 hasta alcanzar 300 °C, temperatura que se mantuvo durante 17 min. El tiempo total de corrida fue de 65 min. El detector de masas se utilizó con una temperatura de la fuente de iones y de la línea de transferencia de 230 °C y 180 °C, respectivamente, y una energía de impacto de 70 eV.

El escaneo de masas entre 30 y 400 uma se realizó para obtener el cromatograma de iones totales (TIC) con el fin de determinar la distribución de n-parafinas y los isoprenoides acíclicos pristano y fitano. Por su parte, el monitoreo de iones seleccionados (modo SIM/Scan) se utilizó para analizar los iones m/z = 128 (naftaleno), 142 (metilnaftalenos), 156 (dimetilnaftalenos), 170 (trimetilnaftalenos), 178 (fenantreno), 184 (tetrametilnaftalenos y dibenzotiofeno), 192 (metilfenantrenos), 198 (metildibenzotiofenos) en la fracción aromática, y los iones m/z = 191 (terpanos) y 217 (esteranos), en el extracto alifático (Cuadro1). De esta manera se obtuvieron los fragmentogramas de interés correspondientes a cada muestra, haciendo uso del software MSD ChemStation Data Analysis. A partir de los tiempos de retención, bibliografía de referencia e integración manual de los picos se identificaron los compuestos antes mencionados.

2.5. Análisis de los datos

Cada una de las diez etapas del ensayo estuvo integrada por tres reactores que se analizaron independientemente y cuyos valores se promediaron para eliminar así el efecto de la variabilidad analítica. Se calcularon 12 RD consideradas a partir de la literatura publicada por ^{Fernández-Varela et al. (2010}), derivadas de isoprenoides cíclicos y comúnmente utilizadas en la identificación de fuentes de contaminación por hidrocarburos para evaluar su recalcitrancia frente a la meteorización. Además, en base al trabajo de ^{Wang y Fingas (2003}) se consideraron otros índices de diagnóstico generados a partir de n-alcanos, pristano, fitano y HAP. Estos parámetros se presentan en el apartado de resultados, más precisamente en el inciso 3.3.

La desviación estándar relativa (DER) se empleó como indicador para evaluar la estabilidad de las moléculas frente a las condiciones de meteorización. Para esto se analizaron las RD en función del tiempo de estudio a través de sus DER obtenidas. Para que una RD se considere apta para su uso a lo largo del tiempo definido su DER debe ser menor al 5 % (^{Zhang et al., 2015}). El cálculo de este parámetro estadístico se llevó a cabo mediante el programa Excel Microsoft Office utilizando la Ecuación 1:

DER = (DE(T0:T9)/P(T0:T9))100 (E.1)

Donde:

DE = Desviación estándar de las diez etapas comprendidas en el ensayo.

P = Promedio de las diez etapas comprendidas en el ensayo.

3.1. Cromatogramas de iones totales (TIC)

El primer paso para determinar si se produjeron alteraciones en las muestras de hidrocarburos sometidas a las condiciones de laboratorio por el período de un año en reactores con suelo es fundamental observar los TIC correspondientes a cada una de las etapas del análisis (Figura1 ). En el tiempo cero (T0) se observó una distribución bimodal de los n-alcanos, desde el undecano (n-C11) hasta el triacontano (n-C30; Figura1A). Durante el primer mes (T1-T4) no se visualizó un cambio apreciable en la composición de los n-alcanos (Figuras1B-E); por otro lado, se evidenció un ligero levantamiento de la línea de base, lo cual es indicativo de un proceso de biodegradación incipiente sobre las muestras de crudo (^{Sutton et al., 2005}).

Luego de dos meses (T5) se consideró a la biodegradación y a la evaporación (Figura1F) como las responsables de la disminución total de los n-alcanos livianos undecano (n-C₁₁) y dodecano (n-C12) y, a los cuatro meses (T6), del tridecano (n-C13). Además, se produjo la formación de una mezcla compleja no resuelta (MCNR) producto del aumento paulatino y progresivo de la línea de base. La presencia de MCNR o ''joroba'' en el cromatograma es una típica evidencia de biodegradación (Figura1G), y consiste en compuestos orgánicos que no se pueden identificar en las condiciones de análisis utilizadas para la separación de los n-alcanos a través de los TIC. En otras palabras, consta de compuestos biorresistentes como saturados cíclicos, aromáticos, naftenoaromáticos, y compuestos polares que no pueden separarse mediante la técnica de cromatografía (^{Sutton et al., 2005}).

Entre los seis y nueve meses del ensayo (T7-T8) tuvo lugar la disminución del tetradecano (n-C₁₄) y el incremento de la MCNR (Figuras1H-I). Es importante considerar que la fracción más liviana de n-alcanos generalmente es la primera en degradarse, por ser la de menor complejidad y de rápida metabolización por los microorganismos (^{Peters et al., 2005}). Finalmente, a los 12 meses de la experiencia (T9) se observó un decrecimiento en las abundancias del heptadecano (n-C₁₇) y del octadecano (n-C₁₈) con respecto a los isoprenoides acíclicos pristano (P) y fitano (F), respectivamente (Figura1J). Este comportamiento es característico de un proceso mediado por microorganismos que produce la degradación de los alcanos lineales, más fáciles de metabolizar que el P y el F, pues estos presentan una estructura ramificada ^{(Peters et al., 2005}).

Figura 1 Cromatogramas de iones totales para los crudos extraídos de suelo a los tiempos T0 (A), T1 (B), T2 (C), T3 (D), T4 (E), T5 (F), T6 (G), T7 (H), T8 (I), T9 (J). 1: tridecano (n-C₁₃), 2: heptadecano (n-C₁₇), 3: octadecano (n-C18), 4: docosano (n-C22), 5: pentacosano (n-C25), 6: octacosano (n-C28), P: pristano, F: fitano. 

Por otro lado, es importante añadir que por encima de la MCNR se observaron señales de alta intensidad para los n-alcanos comprendidos entre el pentadecano (n-C₁₅) y el tridecano (n-C₃₀), lo que sugiere en primera instancia que estas moléculas no se han biodegradado. Estos resultados se pueden comparar con los obtenidos en otros estudios; por ejemplo, en uno de ellos se evaluó la presencia y distribución de n-alcanos y biomarcadores alifáticos en muestras de suelo recogidas en Owaza, en el delta del Níger (Nigeria) durante el lapso de un año; en sus resultados expuso que la fracción alifática estaba dominada por los n-alcanos de n-C₁₉ a n-C₃₁ (^{Faboya et al., 2016}). Otra investigación llevada a cabo en el suelo del Amazonas peruano, vinculada al vertido de hidrocarburos provenientes de instalaciones en esa región, evidenció un estado muy erosionado de los crudos luego de varios años de producido el siniestro, en el cual la distribución de n-alcanos se caracterizó por la ausencia de sus componentes ligeros y por el predominio de los compuestos de alto peso molecular (^{Rosell-Melé et al., 2018}).

3.2 Fragmentogramas de isoprenoides cíclicos

Los fragmentogramas correspondientes a los isoprenoides conocidos como terpanos (ion m/z = 191) se caracterizaron por no presentar cambios evidentes en cada una de las etapas comprendidas dentro de este análisis, con la salvedad de un leve incremento de la línea de base que no modificó la composición ni la distribución de estos biomarcadores (Figuras2A-J). En detalle se detectaron las señales de los terpanos tricíclicos (C23 y C26), pentacíclicos (hopanos C29 y C30) y los hopanos extendidos (homohopanos C31, C32 y C33). Cabe destacar que no se visualizó una disminución de los picos vinculados a los homohopanos (C₃₁ a C₃₃), compuestos susceptibles a la biodegradación; y tampoco del resto de isoprenoides cíclicos pertenecientes al ion m/z = 191. En los yacimientos petrolíferos, la degradación de estos biomarcadores mediada por microorganismos se desencadena por la pérdida de un grupo metilo, que da lugar a la formación de los 25-norhopanos u hopanos desmetilados (^{Peters et al., 2005}), los cuales pueden identificarse a través del fragmentograma de masas asociado al ion m/z = 177. Estas moléculas no se observaron durante este ensayo, lo que refuerza los resultados obtenidos para los homohopanos (datos no presentados).

Por otro lado, en los fragmentogramas de masas para el ion m/z = 217 se observó un descenso de las abundancias relativas de los esteranos (C27, C28 y C29) y de los diasteranos C27 en función del pregnano C21 y homopregnano C22. Sin embargo, no se modificó la intensidad relativa del patrón de distribución de los esteranos regulares (C27, C28, C29) en el transcurso de la experiencia (Figuras3A-J). Estos biomarcadores no presentaron un orden de alteración específico bajo las condiciones de ensayo al cabo de 12 meses; sin embargo, la intensidad de las señales entre C₂₁ y C₂₂ respecto a C27, C28 y C29, sugiere una posible alteración por biodegradación (^{López e Infante, 2021}).

Es importante añadir que experiencias llevadas a cabo en el laboratorio sobre crudos, estudiando el efecto de la meteorización, mostraron que ni los terpanos ni los esteranos fueron afectados por el fenómeno de evaporación, ya que estos biomarcadores se concentraron proporcionalmente en relación con el aumento de los porcentajes de degradación de los hidrocarburos (^{Wang y Fingas, 2003}). Los resultados obtenidos muestran, en mayor o menor medida, la estabilidad de estos isoprenoides cíclicos bajo las condiciones de laboratorio en las cuales transcurrió este ensayo.

Figura 2 Fragmentogramas de la relación masa/carga (m/z) correspondiente al ion 191 para los crudos extraídos de suelo a los tiempos T0 (A), T1 (B), T2 (C), T3 (D), T4 (E), T5 (F), T6 (G), T7 (H), T8 (I), T9 (J). 1: terpano tricíclico C23, 2: terpano tricíclico C26, 3: trisnorneohopano (Ts), 4: trisnorhopano (Tm), 5: norhopano C29, 6: hopano C30, 7-8: homohopanos C31 (R y S), 9-10: homohopanos C32 (R y S), 11-12: homohopanos C33 (R y S). 

La capacidad de los terpanos y de los esteranos de permanecer inalterados cuando los hidrocarburos son sometidos a procesos de meteorización se ha documentado en numerosas investigaciones. Por ejemplo, un estudio llevado a cabo por ^{Rosell-Melé et al. (2018}) en el norte del Amazonas peruano permitió determinar los responsables de la contaminación petrolífera debido a que las muestras de suelo analizadas mantenían los perfiles de terpanos y esteranos constantes luego de varios años de acontecido el vertido de hidrocarburos a partir de la ruptura de un oleoducto. Estos isoprenoides son muy fiables debido a su especificidad, naturaleza hidrofóbica y largo tiempo de residencia en el medio ambiente (^{Seguel et al., 2001}).

Otra contingencia de esta naturaleza tuvo lugar en sistemas litosféricos de Nigeria y fue estudiada a partir de los 12 meses de producida. En este caso se evidenció que el crudo responsable de la contaminación provenía de un yacimiento localizado en el delta del Níger debido a la integridad que mantuvieron las distribuciones de estos biomarcadores (^{Faboya et al., 2016}). La capacidad de los terpanos y esteranos para indicar la fuente responsable del petróleo derramado en la intemperie se debe a su naturaleza refractaria y su alta resistencia a la biodegradación (^{Peters et al., 2005}). Además, el estudio de un derrame en el Golfo de México luego de un año y medio mostró un enriquecimiento relativo de estas moléculas respecto a los alcanos (^{Aeppli et al., 2012}).

Por otro lado, la degradación biológica de estos biomarcadores se ha logrado en condiciones agresivas de laboratorio utilizando cultivos de enriquecimiento aeróbico (^{Douglas et al., 2012}). Además, se ha observado la alteración de estos compuestos en varios estudios de campo que involucraban períodos de tiempo superiores a los cinco años (^{Peters et al., 2005}). Por ejemplo, ocho años después de un experimento en el cual se derramaron hidrocarburos en un manglar de Guadalupe (Antillas Menores), los esteranos y hopanos totales se redujeron en más de un 25 % (^{Munoz et al., 1997}). Se hicieron observaciones similares 20 años después a partir de un vertido de crudo en el Ártico (^{Prince et al., 2003}).

Figura 3 Fragmentogramas de la relación masa/carga (m/z) correspondiente al ion 217 para los crudos extraídos de suelo a los tiempos T0 (A), T1 (B), T2 (C), T3 (D), T4 (E), T5 (F), T6 (G), T7 (H), T8 (I), T9 (J). 1: esterano C20, 2: pregnano C21, 3: homopregnano C22, 4: diasteranos C27, 5: colestanos C27, 6: ergostanos C28, 7: estigmastanos C29. 

3.3 RD de n-alcanos e isoprenoides

La mayoría de las RD presentadas en el Cuadro2 mostraron valores de DER por debajo del 5 % en la mayoría de los casos; por consiguiente, se considera que estas relaciones no han sido afectadas por los factores fisicoquímicos y biológicos inherentes al ambiente en el que se desarrolló el ensayo (^{Zhang et al., 2015}). A partir de las áreas se calcularon los cocientes P/n-C₁₇ y F/n-C₁₈ para ser empleados como índices de biodegradación, ya que los mismos consideran moléculas con grados de resistencia biológica diferente. En este sentido, el pristano y el fitano son más resistentes a la biodegradación que sus pares n-C17 y n-C18 (^{Peters et al., 2005}). Como ambas relaciones sobrepasaron el 5 % de su DER, se puede decir que ha tenido lugar el agotamiento preferencial de los hidrocarburos fácilmente degradables (n-C₁₇ y n-C₁₈) frente a los isoprenoides acíclicos P y F, indicio de que las muestras de crudo en suelo se biodegradaron.

Este comportamiento también se ha observado en numerosos estudios, como los presentados por ^{Aeppli et al. (2012}), ^{Faboya et al. (2016}) y ^{Rosell-Melé et al. (2018}). Por otro lado, la relación (n-C13 + n-C14)/(n-C25 + n-C26) asociada al proceso de evaporación exhibió un valor de DER próximo al 60 % (Cuadro2). En detalle se observó una disminución progresiva en el valor de la RD a medida que el tiempo del ensayo transcurría, lo cual puede explicarse por una mayor velocidad de evaporación de los alcanos n-C13 y n-C14 respecto a sus pares n-C25 y n-C26.

Las RD determinadas para los terpanos y esteranos (Cuadro2) presentaron DER inferiores al 5 % en todos los casos; por lo tanto, se sugiere que estos indicadores no han sido afectados por la evaporación y la biodegradación que tuvo lugar durante la experiencia. Respecto a la primera, cuando se verifica la evaporación de un hidrocarburo, compuestos de elevada masa molecular como estos biomarcadores tienden a aumentar su concentración en el crudo residual sin modificar sus RD (^{Wang et al., 2006}). Por otro lado, la biodegradación no ha tenido el tiempo suficiente bajo las condiciones de ensayo para afectar pronunciadamente las abundancias relativas de estos isoprenoides cíclicos. La recalcitrancia, tanto de terpanos como de esteranos, fue documentada en investigaciones ambientales asociadas a vertidos de hidrocarburos in situ o en simulaciones llevadas a cabo en el laboratorio. Por ejemplo, un estudio publicado por ^{Zhang et al. (2015}), que evaluó las 12 RD de terpanos y esteranos presentadas en el Cuadro2, mostró que los cocientes no superaron el 5 % de la DER luego de tres meses de monitoreo. Por otra parte, la contaminación de suelos con hidrocarburos en el delta del Níger pudo comprobarse a través de la determinación de parámetros tales como Ts/Tm que no se modificaron luego de un año de acontecidos los derrames (^{Faboya et al., 2016}). Además, no se observó ningún cambio en la relación Ts/Tm luego de 24 años de producido un derrame en el estrecho de Magallanes (Wang y Fingas, 2003), y 25 años después del accidente ocurrido en Nipisi, los cocientes H₂₉/H₃₀ y Ts/ Tm todavía permanecían estables.

Cuadro 2 Desviaciones estándar relativas (DER) de las relaciones de diagnóstico (RD) para las muestras de crudo sometidas a meteorización en suelo considerando los diez tiempos definidos para el ensayo. 

3.4 RD e histogramas de HAP

Los HAP se han empleado para conocer el origen y destino ambiental del crudo y sus derivados, ya que son moléculas relativamente estables. Por ejemplo, en numerosas investigaciones vinculadas al derrame de hidrocarburos del Exxon Valdez se ha observado que algunas proporciones de los HAP permanecen constantes en estadios tempranos e intermedios de meteorización (^{Douglas et al., 1996}). Respecto a las mismas, podemos decir que el cociente entre el naftaleno y sus derivados metilados (N₀ + N₁)/N₂ en este trabajo sugiere que un proceso de evaporación ha afectado las muestras de crudo sometidas a meteorización. En el transcurso del año se observó una disminución ligera y constante de esta relación que derivó en una DER cercana al 40 % (Cuadro2); por lo tanto, el naftaleno y los metilnaftalenos se han evaporado a una mayor velocidad que los dimetilnaftalenos, asociándose este comportamiento con el menor peso molecular de los mismos (^{Zhang et al., 2015}).

Cabe resaltar que el resto de las relaciones vinculadas con los HAP no experimentó cambios sustanciales, es decir, se mantuvieron por debajo del 5 % de su DER. Por ejemplo, el cociente 2-MP/1-MP (metilfenantrenos) tiene la capacidad de indicar si los procesos fotooxidativos han afectado a los HAP (^{Olson et al., 2017}), pero como la DER de esta relación no superó el 5 % -y considerando que los reactores del ensayo se colocaron durante todo el experimento fuera de la incidencia directa de luz-, se sugiere que la luz no ha promovido cambios en la composición de las muestras de crudo.

En esta misma línea se han utilizado otras RD que provienen de los HAP, tanto el 4MeDBT como el 2 + 3-MeDBT (metildibenzotiofenos) son menos resistentes a la degradación biológica que el 1-MeDBT y se utilizan para indicar los efectos de la biodegradación (^{Wang y Fingas, 2003}). Las DER del 4-MeDBT/1-MeDBT y del 2 + 3-MeDBT/1-MeDBT fueron inferiores al 3%, lo cual sugiere que estas relaciones aún no han sido afectadas por un proceso de biodegradación considerando el tiempo y condiciones del ensayo. Además, es importante mencionar que no se estimuló la comunidad microbiana del suelo con el uso de nutrientes. Una investigación que empleó microorganismos heterótrofos en un microambiente controlado de laboratorio para estudiar la biorremediación en suelos contaminados con crudo mostró el mismo número de bacterias 55 días después de iniciado el ensayo (^{Mariano et al., 2007}). Esto puede asociarse a una ralentización de la biodegradación vinculada a un cambio en la función microbiana, más que a una reducción de la población bacteriana (^{Greenwood et al., 2008}). Sin embargo, existen antecedentes de estudios que presentaron cambios importantes en las RD de los HAP. Por ejemplo, muestras tomadas en el curso de ocho meses luego del derrame del Hebei Spirit en Corea del Sur evidenciaron una progresiva disminución tanto de 4-MeDBT/1-MeDBT como de 2 + 3-MeDBT/1-MeDBT debido procesos degradativos mediados por microorganismos (^{Yim et al., 2011}). En este sentido, en una experiencia llevada a cabo en una playa de China en la cual se vertieron 50 cm3 de crudo, 95 días después se observó que los indicadores 4-MeDBT/1-MeDBT y 2-MP/1-MP presentaban DER por encima del 5 %: un indicio de biodegradación y fotólisis sobre los hidrocarburos derramados (^{Zhang et al., 2015}).

Los porcentajes relativos de los HAP variaron en cada una de las etapas de este estudio (Figura4). Inicialmente (Figura4A) se observó la presencia de naftaleno (N₀) y sus metilderivados (N1, N2, N3, N4), fenantreno (P0) y metilfenantrenos (P1), dibenzotiofeno (D0) y metildibenzotiofenos (D₁), siendo los trimetilnaftalenos (N₃) -con una abundancia relativa inferior al 30 %- los compuestos mayoritarios. Entre T1 y T4 (Figuras4B-E) se produjo una sostenida disminución del naftaleno y de sus metilderivados que pueden evaporarse de la superficie del suelo u oxidarse por la acción de microorganismos (^{An et al., 2005}). Este agotamiento decreciente con la sustitución de grupos alquilo se asoció a un proceso de evaporación que se ha observado en otros estudios llevados a cabo en el laboratorio (^{Brakstad et al., 2014}; ^{Olson et al., 2017}), o de accidentes como el caso de la plataforma petrolífera Deepwater Horizon, en las que la evaporación de alcanos ligeros y HAP de bajo peso molecular se concluyó como el principal proceso de agotamiento del crudo superficial (^{Brown et al., 2011}).

Retomando el presente estudio, entre los dos y cuatro meses (T5 y T6) se evidenció en los histogramas (Figuras4F-G) una reducción en los valores del naftaleno y sus derivados metilados y, además, un incremento en el fenantreno y en los metilfenantrenos. La adsorción de los HAP al suelo ralentiza su volatilización, efecto que se incrementa cuando hay abundancia de materia orgánica en la matriz del sedimento. El suelo de la Patagonia utilizado en esta experiencia se caracteriza por su contenido de materia orgánica reducido, que le confiere poca capacidad de sorción para los contaminantes orgánicos (^{Toledo et al., 2022}). Además, la hidrofobicidad del naftaleno, el fenantreno y sus metilcompuestos determina que su capacidad para adsorberse en las partículas del suelo sea lenta; por lo tanto, sus masas moleculares constituyen el principal aspecto que condiciona su comportamiento frente a la evaporación (^{An et al., 2005}). En base a lo expuesto, los resultados exhiben un comportamiento esperado, ya que los HAP más pesados permanecen en el suelo por más tiempo debido a una estructura molecular más estable e hidrófoba ^{(Abdulazeez y Fantke, 2017}).

Luego de seis y nueve meses (T7 y T8), los histogramas permitieron observar la continua disminución del naftaleno y sus metilmoléculas y un sostenido incremento del fenantreno y de los metilfenantrenos (P₁) con porcentajes cercanos al 20 y 50 %, respectivamente (Figuras4H-I). En muchos casos, las proporciones relativas de naftalenos en relación con otras series de HAP disminuyen de manera pronunciada por su marcada volatilidad (^{Kao et al., 2015}). Finalmente, luego de un año (T9) tanto el naftaleno como sus derivados metilados presentaron valores muy bajos en relación con las muestras de crudo iniciales (T0) y, por otro lado, el fenantreno y los metilfenantrenos siguieron aumentando sus abundancias relativas (Figura4J). Es interesante añadir que en la Figura4 se observó un patrón tipo escalera desde el naftaleno a los metilfenantrenos. Este comportamiento fue reportado en un estudio sobre crudos derramados en la Antártida, en el cual se atribuyó que la volatilización había sido el principal mecanismo de eliminación de los hidrocarburos en las zonas afectadas (^{Aislabie et al., 1999}). Con respecto a los HAP azufrados, se puede mencionar que el dibenzotiofeno (D₀) permaneció prácticamente inalterado a lo largo de los 12 meses con una baja abundancia relativa y que los metildibenzotiofenos (D₁) sufrieron un paulatino descenso en el mismo lapso del 3 al 2 %.

Figura 4 Histogramas de hidrocarburos aromáticos policíclicos (HAP) para los crudos extraídos del suelo a los tiempos T0 (A), T1 (B), T2 (C), T3 (D), T4 (E), T5 (F), T6 (G), T7 (H), T8 (I), T9 (J). NX: naftaleno y metilderivados, PX: fenantreno y metilderivados, DX: dibenzotiofeno y metilderivados. 

En resumen, conocer el perfil de biomarcadores de un crudo determinado y la estabilidad ambiental de esas moléculas bajo determinadas condiciones de meteorización contribuye de manera significativa al entendimiento, manejo y conservación de los ambientes afectados por este recurso no renovable (^{López e Infante, 2021}). Esto se debe a que los biomarcadores encontrados y su comportamiento durante la intemperización constituyen una clara evidencia de los fenómenos fisicoquímicos que han alterado los hidrocarburos vertidos al ambiente (^{Tomas y Acuña, 2023}). A partir de este análisis se pueden establecer las metodologías necesarias para la remediación de las zonas afectadas (^{Zhang et al., 2015}).