Mapa del stock de carbono orgánico del suelo en una cuenca de la Patagonia Andina, Neuquén, Argentina

FRUGONI, M. CRISTINA; GONZÁLEZ MUSSO, ROMINA; RABINO, ALBERTO

Puerto Madryn

Congreso; XXI Congreso Geológico Argentino; 2022

Asociación Geológica Argentina

El presente trabajo ha sido desarrollado en la cuenca del lago Lácar, próxima a San Martín de los Andes, en laprovincia del Neuquén, Argentina (71° 40´ 15.78"; 40° 6´ 23.49" 71° 20´ 16.15"; 40° 12´ 13.33"). El paisaje es demontañas modeladas por la acción glacial, presentando en su mayor parte pendientes de moderadamenteescarpadas a fuertemente escarpadas, a excepción de ciertas bases de pendiente y terrazas estructurales queposeen suaves inclinaciones. La geología de base corresponde a granitos de la Formación Huechulafquen, nosiendo éste el material de origen de los suelos, ya que toda el área ha sido recubierta por cenizas volcánicasholocénicas, aportadas por los volcanes presentes en la Cordillera de los Andes. Los suelos pertenecen al OrdenAndisoles. Poseen una secuencia de horizontes O-A-B w -C o multisecuencias de éstos debido a las sucesivastefras que dan lugar a suelos enterrados. Son profundos a muy profundos, a excepción de áreas próximas a cimaserosionadas o a arroyos de acción torrencial, en donde el manto de cenizas es más delgado. Son seudolimo-arenosos en superficie, tornándose seudolimosos en profundidad. Las estructuras son fluffy y granular ensuperficie, presentando bloques irregulares en profundidad. Son de color pardo muy oscuro a negro ensuperficie, y pardo-amarillentos en el subsuelo, a excepción de los perfiles que poseen horizontes enterrados quevuelven a ser oscuros. La consistencia es muy friable, no adhesivos, no plásticos, débil a fuertemente untuosos.Son alofánicos, con baja densidad aparente (<0,9 tn.m -3 ), bien drenados, con texturas dominadas por vidrio, porlo que satisfacen los criterios para ser considerados dentro del grupo de los andisoles menos alterados (SoilSurvey Staff, 2014). El clima es templado húmedo, con un importante gradiente de precipitaciones, siendo de3000mm de media anual en el extremo oeste de la cuenca y de 1300mm en el extremo este. La temperaturamedia anual oscila entre los 8° y 10°C. Bajo estas condiciones climáticas, el edafoclima es údico y mésico. Elsubgrupo más ampliamente distribuido en el área es el de Udivitrandes húmicos. También se identificaronUdivitrandes tápticos y líticos. La vegetación dominante es la del Bosque mixto de Nothofagus, especies de hojaancha tanto decidua como siempreverde. La cobertura del suelo es casi completa y en el sotobosque domina laChusquea culeou (caña colihue). Se relevaron 35 puntos de muestreo en donde se abrió una calicata y se realizóla descripción morfológica y toma de muestras del suelo por horizonte. Se determinó el contenido de carbonoorgánico (COS) por el método analítico de combustión húmeda (Walkley-Black) y la densidad aparente por elmétodo de muestra no disturbada con cilindro de volumen conocido. Se calculó el COS almacenado de 0 a 30cmde profundidad aplicando la ecuación que relaciona el porcentaje de carbono con la profundidad, la densidadaparente y el contenido de grava y piedras (FAO 2017, Gardi et al, 2014). Al ser suelos libres de gravas, esteúltimo componente no fue incluido. El cálculo se realizó considerando la información de concentración decarbono y densidad aparente de cada horizonte descripto y muestreado. Para caracterizar la concentración deCOS en porcentaje de 0 a 30cm se calculó la media ponderada de las profundidades parciales hasta alcanzar esaprofundidad. El mismo procedimiento se realizó para la densidad aparente (tn.m -3 ). Para el cálculo del COSalmacenado (tn.ha -1 ), si el horizonte superficial poseía una profundidad inferior a los 30cm, se lo calculaba hastaesa profundidad y se agregaban el o los horizontes subyacentes con sus respectivos COS y densidades aparentes(dado que esta última propiedad influye mucho en la reserva), sumando finalmente cada estrato para obtener elvalor total. Existen varios métodos de interpolación para el mapeo digital de las características del suelo, donde apartir de la ubicación espacial de los sitios de muestreo se predicen los valores de las regiones sin datos. Losmétodos geoestadísticos resultan muy eficientes para el mapeo de la distribución espacial del carbono orgánico.Particularmente el Kriging Ordinario (OK, por sus siglas en inglés) ha sido probado como uno de los métodosque arroja mejores resultados (Pham et al 2019, Liu et al 2006). OK incorpora las propiedades estadísticas de losdatos de provenientes del muestreo (autocorrelación espacial) y utiliza un semivariograma para expresar lacontinuidad espacial de los datos y la correlación espacial en función a la distancia entre sitios (Bhunia et al2016). El mapeo digital por el método de interpolación OK se llevó a cabo en R (R Team 2020). Para ello,primero se calculó un variograma experimental y posteriormente se ajustó a un variograma teórico que expliquela distribución de los datos. Una vez ajustado, se realizó la interpolación de datos para predecir los sitios nomuestreados sobre una grilla de 20 metros de resolución espacial. Finalmente, los resultados fueron exportadosen formato ráster e incorporados en un Sistema de Información Geográfica (SIG) usando el sofware QGIS(QGIS 2020) para su visualización, análisis y generación de cartografía. La densidad aparente espredominantemente de 0,8 tn.m -3 , aunque se han medido densidades de hasta 0,5 tn.m -3 , lo que indica una altaporosidad del suelo (70-80%) y por lo tanto una baja cantidad de fracción sólida. La concentración de carbonopresenta una importante variabilidad, característica muy frecuente en esta propiedad del suelo, aunque los valores obtenidos son similares a otros andisoles forestales (Takahashi et al, 2010) y a estudios regionalesrealizados por otros autores (Ferrer et al, 1990) (Tabla 1). La representación en una malla continua del contenidode carbono generada a partir de la interpolación espacial permitió calcular superficies, identificar sitios de mayorreserva y estimar la reserva total de carbono en el área de estudio. Es así entonces que el área cubierta por elmapa alcanza las 28.600 ha. La reserva total de carbono calculada en esa superficie es de 39,23 Mt. Además, elmapa grafica claramente la variabilidad en las reservas de COS presentes en la cuenca (Figura 1). Se observa unasubcuenca en la margen norte del lago Lácar que posee una reserva de 120 tn.ha -1 ; mientras que en otrossectores, como la cabecera este del lago acumula entre 50 y 70 tn.ha -1 . Cabe destacar que esta cuenca integra unárea de reserva del Parque Nacional Lanín, lo que ha permitido una importante intervención maderera desdemediados del Siglo XX y hasta la década del ?80. Por otra parte, existen áreas históricamente ocupadas porpobladores, quienes poseen ganado caprino, equino y ovino. Es uno de los principales atractivos turísticos queposee la ciudad, por lo que ha sido siempre importante la presión antrópica en toda la cuenca. Estos usos de latierra podrían explicar en alguna medida las importantes variaciones en las reservas del COS. Los andisoles, queposeen una capacidad de conservar importantes cantidades de COS (se han medido valores de 6%), al poseer tanbajas densidades aparentes y por lo tanto una baja proporción de fracción sólida, dan valores no tan elevados deCOS almacenado. Sería importante entonces, cuando no se cuenta con el dato, ajustar el valor de la densidadaparente con modelos confiables -particularmente en estos suelos con propiedades únicas- a fin de minimizar loserrores que se puedan cometer en el cálculo final de esta propiedad.