Repositorio Institucional de Documentos

APLICACIÓN DEL MÉTODO LANDSCAPE FUNCTION ANALYSIS (LFA) PARA EVALUAR EL ESTADO DE LAS RESTAURACIONES MINERAS. EL CASO DE LA MINERÍA DE CARBÓN A CIELO ABIERTO EN UTRILLAS (TERUEL) Realizado por: Juan Alejandro Oyarzábal Mateo Dirigido por: José Manuel Nicolau Ibarra Ciencias Ambientales Octubre 2013 Escuela Politécnica Superior de Huesca [“APLICACIÓN DEL MÉTODO LANDSCAPE FUNCTION ANALYSIS (LFA) PARA EVALUAR EL ESTADO DE LAS RESTAURACIONES MINERAS. EL CASO DE LA MINERÍA DE CARBÓN A CIELO ABIERTO EN UTRILLAS (TERUEL)”.] 2013 2 [“APLICACIÓN DEL MÉTODO LANDSCAPE FUNCTION ANALYSIS (LFA) PARA EVALUAR EL ESTADO DE LAS RESTAURACIONES MINERAS. EL CASO DE LA MINERÍA DE CARBÓN A CIELO ABIERTO EN UTRILLAS (TERUEL)”.] 2013 3 APLICACIÓN DEL MÉTODO LANDSCAPE FUNCTION ANALYSIS (LFA) PARA EVALUAR EL ESTADO DE LAS RESTAURACIONES MINERAS. EL CASO DE LA MINERÍA DE CARBÓN A CIELO ABIERTO EN UTRILLAS (TERUEL) Realizado por: Juan Alejandro Oyarzábal Mateo Dirigido por: José Manuel Nicolau Ibarra Ciencias Ambientales Octubre 2013 Escuela Politécnica Superior de Huesca [“APLICACIÓN DEL MÉTODO LANDSCAPE FUNCTION ANALYSIS (LFA) PARA EVALUAR EL ESTADO DE LAS RESTAURACIONES MINERAS. EL CASO DE LA MINERÍA DE CARBÓN A CIELO ABIERTO EN UTRILLAS (TERUEL)”.] 2013 4 Resumen La minería es una actividad necesaria para la sociedad que produce un intenso impacto ambiental y cierto rechazo social. Para poder compatibilizar minería con medio ambiente la restauración es una herramienta necesaria aunque no siempre suficiente. Las restauraciones mineras se encuentran reguladas por la ley desde 1982 en España. Sin embargo, las normativas no contemplan ningún método empírico para la evaluación del resultado de las labores de restauración. Este es un aspecto importante desde el punto de vista ambiental, pero también administrativo, ya que para la devolución del aval económico depositado por las empresas mineras se requiere un informe técnico sobre el éxito de la restauración. Para cubrir la falta de una metodología de evaluación de restauraciones, en este Trabajo de Fin de Grado se evalúa la validez del método LFA (Landscape Function Analysis), desarrollado en Australia, como herramienta de aplicación en la cuenca minera de Utrillas (Teruel). LFA es un método de evaluación del estado de degradación de ecosistemas en ambientes secos, que se ha aplicado con éxito en zonas mineras restauradas. LFA se fundamenta en tres tipos de índices: Índice de Estabilidad, Índice de Infiltración e Índice de Reciclado de Nutrientes. El método se basa en los procesos hidrológicos de superficie, teniendo que realizar una serie transectos en el terreno que se quiere evaluar. Sobre estos últimos, hay que medir: a) La organización espacial de los parches de vegetación (midiendo longitud, anchura y distancia entre ellos) y b) Evaluación de las características de la superficie del suelo mediante 11 indicadores que se agrupan en los 3 índices enunciados anteriormente. El método se ha aplicado en un conjunto de 6 taludes restaurados que representan un gradiente de degradación-agradación, es decir un gradiente de restauraciones fallidas a exitosas. Sobre estas laderas existe información publicada relativa a su funcionamiento ecohidrológico, lo que ha permitido comparar las estimaciones del LFA con las mediciones empíricas publicadas sobre las mismas. Adicionalmente, se ha aplicado también sobre una ladera natural, un encinar tomado como ecosistema de referencia. Los resultados han mostrado las diferencias –estadísticamente significativasentre los taludes según su estado de conservación y, también, que tales diferencias se corresponden con las encontradas mediante la aproximación empírica de tipo ecohidrológico. Así se ha diferenciado: taludes con regueros y escasa cobertura vegetal con estado de degradación; taludes que se encuentran en una situación de incertidumbre en su evolución, con presencia de regueros discontinuos y cobertura vegetal más estable. Y finalmente un talud que se encuentra en un estado de agradación, con ausencia de regueros y vegetación continua (dominada por Genista scorpius) mostrando un estado funcional mejor que el propio ecosistema de referencia. Por último se han identificado 3 tipos de comportamientos ecohidrológicos de los parches de vegetación: a) Generadores de escorrentía (exportadores), b) Sumideros (importadores) y c) Según las condiciones y circunstancias pueden actuar como sumideros o exportadores. [“APLICACIÓN DEL MÉTODO LANDSCAPE FUNCTION ANALYSIS (LFA) PARA EVALUAR EL ESTADO DE LAS RESTAURACIONES MINERAS. EL CASO DE LA MINERÍA DE CARBÓN A CIELO ABIERTO EN UTRILLAS (TERUEL)”.] 2013 5 [“APLICACIÓN DEL MÉTODO LANDSCAPE FUNCTION ANALYSIS (LFA) PARA EVALUAR EL ESTADO DE LAS RESTAURACIONES MINERAS. EL CASO DE LA MINERÍA DE CARBÓN A CIELO ABIERTO EN UTRILLAS (TERUEL)”.] 2013 6 Agradecimientos Me gustaría dar las gracias a José Manuel Nicolau, por la oportunidad que me ha dado de hacer este proyecto, por todo lo que me ha enseñado, por solucionar todas mis dudas y ayudarme en todo lo posible durante el muestreo en Utrillas con su enorme amabilidad. Gracias a mi familia y amig@s, por su ánimo y apoyo. Gracias a al grupo de Ecología de Alcalá de Henares (Mariano y Tíscar) por ayudarme en las dudas surgidas durante el muestreo y por el buen rato pasado. [“APLICACIÓN DEL MÉTODO LANDSCAPE FUNCTION ANALYSIS (LFA) PARA EVALUAR EL ESTADO DE LAS RESTAURACIONES MINERAS. EL CASO DE LA MINERÍA DE CARBÓN A CIELO ABIERTO EN UTRILLAS (TERUEL)”.] 2013 7 Índice del Trabajo fin de Grado 1. Introducción 13 1.1 La necesidad de las restauraciones en la minería 13 1.2 La evaluación del éxito de las restauraciones mineras 14 1.3 Objetivos del TFG (Trabajo Fin de Grado) 16 2. Método LFA 17 2.1 Marco conceptual 17 2.1.1 Ecohidrología en ambientes Mediterráneos secos 17 2.1.2 La determinación de la heterogeneidad en los ambientes semiáridos 18 2.1.3 Modelo TTRP: Marco conceptual del método LandScape Function Analysis 19 2.2 Descripción del método 20 3. Material y métodos 29 3.1 Área de estudio 29 3.1.1 Descripción de las restauraciones mineras en el área de estudio 33 3.1.2 Comportamiento ecohidrológico de las laderas restauradas 36 3.2 Protocolo de muestreo 38 3.3 Análisis de datos 44 4. Resultados 46 4.1 Resultados a nivel de ladera 46 4.1.1 Caracterización de las laderas 46 4.1.2 Caracterización del paisaje 50 4.1.3 Validación del método LFA para la cuenca de Utrillas 53 4.2 Resultados a nivel de parche 56 4.2.1 Contribución de cada tipo de parche/interparche a los índices 57 4.2.2 Caracterización funcional de cada tipo de parche/interparche 58 4.2.3 Comparación cualitativa con la caracterización de parches empírica 62 5. Discusión 63 5.1 Escala de ladera 63 5.2 Escala de parche 65 [“APLICACIÓN DEL MÉTODO LANDSCAPE FUNCTION ANALYSIS (LFA) PARA EVALUAR EL ESTADO DE LAS RESTAURACIONES MINERAS. EL CASO DE LA MINERÍA DE CARBÓN A CIELO ABIERTO EN UTRILLAS (TERUEL)”.] 2013 8 6. Conclusiones 68 7. Bibliografía 69 8. Anexos 74 8.1 Anexo 1 74 8.1.1 La minería del carbón y su deterioro ambiental 74 8.2.1 Evolución histórica de la minería del carbón en la provincia de Teruel 75 8.2 Anexo 2 79 8.2.1 Marco legal y normativa aplicable a la restauración 79 8.2.2 Disposiciones específicas medioambientales y particulares de aplicación a la Minería. 80 8.2.3 Normativa específica para la restauración en explotaciones de carbón a cielo abierto. 82 8.3 Anexo 3 85 8.3.1Metodología LFA 85 8.3.1.1 Ejemplo y ficha de trabajo para la caracterización de la organización del paisaje. 85 8.3.1.2 Ficha de trabajo para la evaluación de los indicadores del suelo 86 8.3.1.3 Indicadores de las características del suelo 86 8.3.2 Fotografías de cada tipo de parche/interparche 89 8.3.3 Tablas de referencia de los indicadores del suelo 94 8.4.4 Coordenadas de los transectos 97 8.4 Anexo 4 98 8.4.1 Datos de los resultados del Índice de Infiltración, Estabilidad y Reciclaje de Nutrientes para cada transecto y ladera. 98 8.4.2 Datos de los resultados de la organización del paisaje para cada transecto y ladera. 100 8.4.3 Datos y resultados del comportamiento individual de cada tipo de parche /interparche respecto al IE, II Y IRN. 103 8.4.4 Contribución (%) de cada parche/interparche individual al IE, II, IRN. 106 8.4.5 Datos de los resultados de la caracterización de cada tipo de parche/ interparche basados en los índices funcionales del paisaje 111 8.4.6 Datos y resultados del análisis Kruskal-Wallis. 112 [“APLICACIÓN DEL MÉTODO LANDSCAPE FUNCTION ANALYSIS (LFA) PARA EVALUAR EL ESTADO DE LAS RESTAURACIONES MINERAS. EL CASO DE LA MINERÍA DE CARBÓN A CIELO ABIERTO EN UTRILLAS (TERUEL)”.] 2013 9 8.4.7 Resultados del análisis estadístico de correlación de Spearman para la validación del método LFA. 122 [“APLICACIÓN DEL MÉTODO LANDSCAPE FUNCTION ANALYSIS (LFA) PARA EVALUAR EL ESTADO DE LAS RESTAURACIONES MINERAS. EL CASO DE LA MINERÍA DE CARBÓN A CIELO ABIERTO EN UTRILLAS (TERUEL)”.] 2013 16 diversidad similar y estructura de la comunidad en comparación con sitios de referencia, (2) la presencia de especies autóctonas, (3) la presencia de grupos funcionales necesarios para la estabilidad a largo plazo, (4) capacidad del medio físico para mantener poblaciones reproductoras, (5) el funcionamiento normal, (6) la integración con el paisaje, (7) la eliminación de las amenazas potenciales; (8) capacidad de resistencia a las perturbaciones naturales, y (9) autosostenibilidad. Aunque las mediciones de estos atributos proporcionan una excelente evaluación de éxito de la restauración, algunos estudios no cuentan con los recursos financieros para supervisar todos estos atributos. Por otra parte, las estimaciones de muchos atributos a menudo requieren estudios detallados a largo plazo, pero la fase de seguimiento de la mayoría de los proyectos de restauración rara vez dura más de 5 años. Para resumir este apartado, se ha observado que en España existe una carencia en cuanto a los métodos de evaluación de restauraciones mineras, las cuales son de obligatoriedad por parte de la legislación nacional después de haber realizado una actividad extractiva, siendo la Generalitat de Catalunya la única que posee un método adecuado, claro y conciso, que se ha desarrollado en base a otros métodos del ámbito internacional como los que han sido enunciados anteriormente que permiten evaluar adecuadamente un restauración minera. 1.3 Objetivos del TFG El objetivo general de este TFG es comprobar la validez del método LFA en la evaluación de las restauraciones de la minería del carbón a cielo abierto en Utrillas (Teruel). Para ello se han perseguido dos objetivos específicos: 1. Comparar el estado de las restauraciones obtenido a través del método LFA en un conjunto de siete taludes que representan un gradiente del estado ecológico de las restauraciones (fracaso-éxito), con los resultados obtenidos en ellos por otros autores mediante el monitoreo de diversos procesos ecohidrológicos. 2. Comparar el estado de las restauraciones obtenido a través del método LFA con el de un encinar próximo, considerado como ecosistema de referencia, aplicando el mismo método LFA. [“APLICACIÓN DEL MÉTODO LANDSCAPE FUNCTION ANALYSIS (LFA) PARA EVALUAR EL ESTADO DE LAS RESTAURACIONES MINERAS. EL CASO DE LA MINERÍA DE CARBÓN A CIELO ABIERTO EN UTRILLAS (TERUEL)”.] 2013 17 2.Método LFA 2.1 Marco conceptual 2.1.1 Ecohidrología en ambientes Mediterráneos secos Hay que destacar que en los ambientes mediterráneos secos, a los cuales pertenece nuestra área de estudio, se originan precipitaciones inferiores a los 500 mm, lo que origina que los procesos de drenaje sean muy limitados, por lo que la mayoría del agua de las precipitaciones se pierde a través de procesos de evapotranspiración. Por debajo de esos niveles de disponibilidad hídrica, la presencia de una cubierta vegetal estable y continua es casi imposible (Cortina et al 2004). Si pasamos al ámbito del clima mediterráneo-continental, presenta una serie de condiciones que dificultan de forma considerable el proceso de desarrollo vegetal, que estará limitado a dos periodos vegetativos correspondientes, la primavera y el otoño, los cuales son cortos y están separados temporalmente. Por una parte tenemos la escasez hídrica estival y los fríos y duros inviernos que originan que las plantas estén en parada biológica durante un gran parte del año. Si a esto, le sumamos la discontinuidad de las precipitaciones y en ocasiones de carácter torrencial de estas, se producen procesos erosivos muy intensos que seguidos de grandes periodos sin lluvia, originan que el establecimiento de la vegetación sea dificilísimo. Hay que tener en cuenta que la erosión además limita la disponibilidad hídrica ya que modifica las características superficiales del suelo (profundidad del suelo funcional para el desarrollo vegetal, disminución del materia orgánica, encostramiento superficial…) lo que produce un empobrecimiento de la comunidad vegetal, la cual disminuye la capacidad de esta para protegerse de la erosión hídrica (Nicolau 2003, Moreno de las Heras, 2009). Debemos tener en cuenta que la erosión alcanza el valor máximo en zonas con 300 mm de precipitación anual (García-Fayos, 2004). El factor de la vegetación es fundamental sobre el control de la intensidad de los efectos de la erosión hídrica, ya que incrementa la propia coherencia del suelo y la estabilidad de los agregados del mismo, protección frente al daño producido por el impacto de las gotas de lluvia, además de mejorar la capacidad de infiltración y actuando como obstáculo para la escorrentía. Si la vegetación no es capaz de controlar estos factores, la erosión va disminuyendo la cantidad efectiva de recursos disponibles para su correcto desarrollo, que puede deberse a el arrastre y el transporte a otros lugares o por un deterioro de las características físico-químicas del suelo que permiten captar y almacenar los propios recursos, lo que dificulta considerablemente el tan importante proceso de sucesión ecológica (Nicolau, 2003). La erosión hídrica afecta a la capa superficial del suelo, donde se produce la acumulación de restos vegetales, nutrientes, agua y propágulos. En este horizonte superficial se produce la germinación de las semillas, por lo que la pérdida por erosión del mismo supone una importante perturbación y el consecuente reforzamiento del estrés hídrico de la cubierta vegetal. Por eso la medición de las características superficiales del suelo es un importante indicador del estado funcional del ecosistema. [“APLICACIÓN DEL MÉTODO LANDSCAPE FUNCTION ANALYSIS (LFA) PARA EVALUAR EL ESTADO DE LAS RESTAURACIONES MINERAS. EL CASO DE LA MINERÍA DE CARBÓN A CIELO ABIERTO EN UTRILLAS (TERUEL)”.] 2013 18 Esta escorrentía superficial se comporta como un factor clave en los ambientes restaurados y su control es vital para el devenir de la restauración. Ya que puede actuar como redistribuidor de un recurso escaso, el agua; o como un importante agente perturbador del ecosistema. El dominio y control experto de estos procesos de escorrentía dentro de sistemas que presentan estados degradados, resulta una herramienta prioritaria para emprender procesos de restauración de estos ecosistemas (Nicolau, 2003). Existe una relación determinante la cual afecta a la dinámica erosiva del sistema entre las formas de superficie, el suelo y la vegetación, teniendo en cuenta el régimen climático. Cada suelo presenta un umbral crítico entre la erosión que puede soportar y la vegetación que puede sustentar. Si este umbral es sobrepasado, la disminución de la cobertura vegetal e incremento de la erosión debido a la falta de vegetación, puede llegar a alcanzar un estado de desertificación irreversible (García-Fayos, P. 2004). La composición de la vegetación y el patrón espacial de ésta depende y es influida por la erosión hídrica (Espigares et al, 2009), pero también el patrón espacial vegetal actúa modificando la dinámica hídrica del sistema de forma reversible. Los procesos erosivos son claves tanto en la estructura de la vegetación (cobertura o biomasa) como en la riqueza de especies, ya que las especies vegetales que están presentes en estos ambientes tan duros presentan determinadas características funcionales que les debe permitir sobrevivir (García-Fayos, P. 2004). Para resumir este apartado, podemos establecer que hay dos características de los ambientes mediterráneos secos:  Presencia de una cubierta vegetal discontinua debido a la baja disponibilidad de recursos y heterogeneidad ambiental a pequeña escala. Cabe destacar que la cobertura vegetal en estos ambientes, va desde un 100% a menos de 20% según los condiciones que se den (Cortina et al, 2004)  Cobertura escasa, generando “islas de fertilidad” rodeadas de huecos desnudos. De estas dos características, concluimos que la cubierta vegetal está limitada por el recurso hídrico, lo que impide una cubierta frondosa y continua, que sumado a procesos erosivos originados por precipitaciones intensas a los cuales no pueden hacer frente o mitigar su intensidad, y unas características del suelo que no permiten una rápida y eficiente infiltración del exceso de agua, el problema se agrava. 2.1.2 La determinación de la heterogeneidad ambiental en los ambientes semiáridos A pequeña escala, si se presenta una heterogeneidad espacial generada, se produce una variación importante en el comportamiento de la superficie del suelo (Tongway et al, 2005) como sucede en la hidrogeomorfología, que varía notablemente debido a la disponibilidad de sedimentos, infiltración o rugosidad. Estos cambios generan un diferenciación clave en los tipos de superficies o cubiertas del terreno, en las que habrá elementos que actúen como fuente o como sumidero según la escorrentía superficial, sedimentos y propágulos. [“APLICACIÓN DEL MÉTODO LANDSCAPE FUNCTION ANALYSIS (LFA) PARA EVALUAR EL ESTADO DE LAS RESTAURACIONES MINERAS. EL CASO DE LA MINERÍA DE CARBÓN A CIELO ABIERTO EN UTRILLAS (TERUEL)”.] 2013 19 Según Tongway et al (2004a) la producción de biomasa es mayor en los medios semiáridos si la lluvia se concentra en algunas zonas, en vez de que está se distribuya uniformemente por la superficie del suelo. Esta afirmación es clave en la ecología de paisaje en las zonas áridas (Maestre, 2003). Es decir, si se producen precipitaciones capaces de generar escorrentía, y estas, se concentran en determinadas zonas, poseerán una mayor disponibilidad de condiciones hídricas generándose las llamadas “islas de fertilidad”. Estas zonas presentan mayores concentraciones de otros recursos (nutrientes, materia orgánica) o propágulos (semillas, esporas de hongos, etc.), ya que la escorrentía superficial transporta los nutrientes disueltos, restos de materia orgánica y sedimentos que acaban siendo retenidos en estos puntos de máxima fertilidad local. Por eso en estos ambientes áridos o semiáridos, estas manchas o parches de vegetación perennes son tan importantes, actuando como colectores de recursos que son escasos (Martín et al, 1999). Otra interacción muy importante es la de tipo “facilitación” de planta a planta, que es vital también para formar estas “islas de fertilidad” o micrositios con mejores condiciones para la dispersión, germinación, crecimiento y supervivencia de las plántulas y semillas. Un ejemplo de ello es la facilitación entre plántulas de especies arbóreas y arbustivas en ambientes mediterráneos, siendo un fenómeno generalizado (Zamora, 2001; Maestre, 2003a; Maestre et al, 2003d; Marañón et al, 2004). Podemos concluir, que la vegetación mejora las condiciones edáficas y microclimáticas en sus inmediaciones en comparación a los espacios de suelo desnudo contiguos, formándose lo que hemos enunciado antes, las “islas de fertilidad”, donde la relación facilitación-competencia entre plantas presenta un balance positivo a favor de la facilitación (Maestre, 2003a). El modelo de funcionamiento de este tipo de ecosistemas semiáridos se caracteriza en estar formados por dos unidades o zonas, las zonas fuente (suelo desnudo) ya que permiten el transporte a través de ella de esos recursos tan importantes que hemos enunciado antes y las zonas sumidero (plantas vasculares) capaces de retener esos recursos procedentes de los sumideros (Cortina et al, 2004). La función de la planta en estas zonas llamadas sumidero, la cual mejora las condiciones microclimáticas, es aprovechar el recurso hídrico de forma más eficiente, de forma que las tasas de evapotranspiración disminuyan, mejorando las propiedades físicas y químicas del suelo, además de la habilidad para la retención de los propágulos debido a un transporte anemocoro. En suma, estas manchas pueden ser empleadas por las especies vegetales para avanzar etapas de la sucesión al emplearlas como puntos de refugio en áreas degradadas (Barbera et al, 2006). Este efecto nodriza o niñera ha sido ampliamente comprobado en el ámbito semiárido mediterráneo (Barbera et al. 2006; Gómez-Aparicio et al. 2004, Maestre, 2003a). 2.3.1 Modelo TTRP: El marco conceptual del método LFA (Landscape Function Analysis) Ludwig y Tongway (1995) proponen un modelo conceptual que se ajusta muy bien al concepto de heterogeneidad ambiental y procesos ecohidrológicos en ambientes áridos y semiáridos. La heterogeneidad espacial está formada por la alternancia de suelo desnudo con ausencia de vegetación (claros) y por manchas de vegetación con variaciones en su tamaño (matas) como consecuencia de un flujo de recursos que recorre el sistema. En este tipo de ecosistemas [“APLICACIÓN DEL MÉTODO LANDSCAPE FUNCTION ANALYSIS (LFA) PARA EVALUAR EL ESTADO DE LAS RESTAURACIONES MINERAS. EL CASO DE LA MINERÍA DE CARBÓN A CIELO ABIERTO EN UTRILLAS (TERUEL)”.] 2013 20 intervienen procesos como el transporte, la utilización, el ciclo de los escasos recursos y limitaciones como el agua, el suelo, la materia orgánica y propágulos en un contexto temporal y espacial (Ludwig et al. 1995). Este enfoque explica el funcionamiento del ecosistema y las diferencias que existen entre la composición biológica y la estructura de este. El marco conceptual TTRP (Desencadenante-Transferencia-Reserva-Pulso) se basa en lo enunciado anteriormente, es decir, en el análisis de los procesos involucrados en el correcto funcionamiento del paisaje que son capaces de controlar la forma de utilizar eficientemente los recursos limitantes. Si consideramos, según este TFG, una ladera como elemento de estudio, entendiendo que actúa como un sistema biofísico, estaría formado por “parches” (islas vegetales) y zonas intermedias de suelo desnudo (interparches o claros), donde las zonas de suelo desnudo ejercen una función de fuente de recursos o inputs (agua, sedimentos y propágulos) en dirección a hacia los parches, matas o islas, las cuales generan una serie de procesos de retroalimentación dentro del sistema y a escala de ladera. A continuación se presenta un esquema del modelo TTRP (Figura 1): Figura 1. Esquema del Modelo conceptual Desencadenante-Transferencia-Reserva-Pulso (Trigger– transfer–reserve–pulse, TTRP). Tomado de Tongway et al., 1995 Este modelo explica, que las precipitaciones actúan como “Desencadenante”, ya que con una determinada intensidad genera una escorrentía superficial que va redistribuir (“Transferencia”) los recursos a través del sistema. Los recursos que se encuentran en las zonas intermedias o claros son transportados por la escorrentía hasta ser retenidos y acumulados en otros puntos del sistema, que son los parches, matas o islas. El resultado de este proceso genera una distribución desigual de los recursos y por tanto una diferenciación en las características como la infiltración y disponibilidad hídrica, materia orgánica, propágulos…). Si se dan las condiciones se genera un “Pulso” en la productividad vegetal y la mineralización de nutrientes en las zonas receptoras de estos aportes. Parte de esa productividad vegetal se puede perder o verse disminuida debido a procesos perturbadores como el fuego, la herbivoría o puede quedarse acumulada en forma de “Reserva”. Este “Pulso” puede general procesos de retroalimentación y modificar el proceso de “Transferencia” al modificarse las [“APLICACIÓN DEL MÉTODO LANDSCAPE FUNCTION ANALYSIS (LFA) PARA EVALUAR EL ESTADO DE LAS RESTAURACIONES MINERAS. EL CASO DE LA MINERÍA DE CARBÓN A CIELO ABIERTO EN UTRILLAS (TERUEL)”.] 2013 21 condiciones físicas del suelo. Es este proceso de retroalimentación o autorregulación es la clave de la dinámica sucesional del ecosistema, y por lo tanto un factor determinante que caracteriza el estado en que se encuentra el sistema. Por ejemplo, si el sistema posee más pérdidas de los recursos que reservas, la dinámica del sistema va a ser negativa, mientras que si es al revés, es decir, hay un mayor control de los recursos y disminuyen las fugas, el sistema está en proceso de una rehabilitación exitosa. 2.2 Descripción del método A continuación se va explicar en qué consiste el método que se ha utilizado en este trabajo, Landscape Function Analysis, desarrollado por D J Tongway and N L Hindley (2004). Se trata de un método mediante el cual podemos evaluar el estado funcional de un sistema biofísico, basándonos principalmente en los procesos implicados en la hidrología de superficie: precipitación, infiltración, escorrentía, erosión, crecimiento de las plantas y ciclos de nutrientes. La metodología va consistir en:  Caracterización de organización del paisaje (Distribución espacial de los parches fértiles e interparches): El objetivo de este paso es caracterizar y cartografiar el sitio monitoreado en términos pérdida o acumulación de recursos en un patrón espacial. Es decir, analizar cómo es la distribución de los diferentes tipos de parches fértiles e interparches, en este caso a lo largo de las laderas de estudio. Para poder realizar esta caracterización se realizan una serie de transectos con una cinta métrica (Figura 2). Figura 2. Foto durante la realización de un transecto con la cinta métrica [“APLICACIÓN DEL MÉTODO LANDSCAPE FUNCTION ANALYSIS (LFA) PARA EVALUAR EL ESTADO DE LAS RESTAURACIONES MINERAS. EL CASO DE LA MINERÍA DE CARBÓN A CIELO ABIERTO EN UTRILLAS (TERUEL)”.] 2013 22 Los transectos se deben iniciar en el borde superior de la pendiente de la cuenca hidrográfica local, es decir, en la parte superior de la ladera (Figura 3). Estos transectos se harán a favor de la pendiente máxima (dirección preferente del flujo de recursos y aguan en el ecosistema) y en dirección de los vientos más fuertes. Figura 3. Transecto realizado en una ladera En cada transecto habrá que realizar una serie de mediciones, para ello deberemos definir claramente que es un parche e interparche. Un parche es elemento físico perdurable en el tiempo (perenne) que tenga capacidad para la acumulación de recursos de agua y sedimentos arrastrados por la erosión superficial, actuando como un área o elemento de captación y retención de recursos dentro de la ladera (Figura 4). Figura 4. Ejemplo gráfico de un parche [“APLICACIÓN DEL MÉTODO LANDSCAPE FUNCTION ANALYSIS (LFA) PARA EVALUAR EL ESTADO DE LAS RESTAURACIONES MINERAS. EL CASO DE LA MINERÍA DE CARBÓN A CIELO ABIERTO EN UTRILLAS (TERUEL)”.] 2013 23 Una vez explicado el parche nos disponemos a explicar un interparche, que es el espacio comprendido entre dos parches, el cual se caracteriza por ser una zona donde los recursos como el agua, los materiales del suelo y la materia orgánica son libremente transportados pendiente abajo cuando el transporte hídrico y eólico están activos. En la figura 5 se presenta la sección de un transecto donde se muestran los límites de parches e interparches. Figura 5. Sección de un transecto donde se muestran los límites de los parches e interparches Las mediciones que hay que realizar en cada transecto son las siguientes: 1. Nº y tipología de las obstrucciones (parches) del flujo superficial por unidad de longitud del transecto. Nota: Se considera como parche plantas superiores a 1cm y rocas de más de 10 cm. 2. Anchura de las obstrucciones por unidad de longitud del transecto 3. Distancia media, y relación, entre obstrucciones (longitud interparche) por unidad de longitud del transecto. En las figuras 6 y 7 se presentan una serie de imágenes de lo explicado anteriormente: [“APLICACIÓN DEL MÉTODO LANDSCAPE FUNCTION ANALYSIS (LFA) PARA EVALUAR EL ESTADO DE LAS RESTAURACIONES MINERAS. EL CASO DE LA MINERÍA DE CARBÓN A CIELO ABIERTO EN UTRILLAS (TERUEL)”.] 2013 24 Figura 6. Ilustración de las mediciones que hay que realizar en los parches/interparches a lo largo de transecto Figura 7. Ilustración esquemática que muestra la manera de medir la y la anchura de los parches. Los procesos identificados a esta escala son muy informativos sobre la función y el estado de salud de la vegetación, que en estos ambientes realizan una regulación de los escasos recursos. En particular, el papel de la disposición espacial de las plantas perennes, es fundamental a la hora de detener el flujo de agua de escorrentía, filtrar los sedimentos y materia orgánica, procesos de vital importancia. Mediante este método es posible medir rápidamente la disposición espacial de la vegetación, reflejando el control que tiene ésta sobre la escorrentía y erosión.  Evaluación de la superficie del suelo de cada tipo de parche/interparche En los mismos transectos medidos anteriormente en la caracterización de organización del paisaje, se realizarán una toma de muestras con la finalidad de medir unos parámetros que caractericen el comportamiento de cada tipo de parche e interparche presentes en el transecto. Para ello se tomarán 5 réplicas de cada tipo de [“APLICACIÓN DEL MÉTODO LANDSCAPE FUNCTION ANALYSIS (LFA) PARA EVALUAR EL ESTADO DE LAS RESTAURACIONES MINERAS. EL CASO DE LA MINERÍA DE CARBÓN A CIELO ABIERTO EN UTRILLAS (TERUEL)”.] 2013 25 parche/interparche con el objetivo de conseguir una mayor formalidad estadística de los resultados obtenidos. En este muestreo a través del transecto se evitará una selección de los sitios más interesantes, distribuyéndose a lo largo de todo el transecto. Los datos que se toman, son de una serie de indicadores de las características superficiales del suelo (Tabla 1). Indicadores/variable tomada Significado Índice en que se emplea 1. Protección frente al golpeo de gotas de lluvia Estima la vulnerabilidad a la erosión por golpeo de las gotas de lluvia Estabilidad 2. Cobertura vegetal perenne Evalúa la contribución de la biomasa de raíces a los procesos de reciclaje de nutrientes Infiltración Reciclado de nutrientes 3. Restos vegetales depositados Indica la disponibilidad de materia orgánica superficial para la descomposición y el reciclado de nutrientes Infiltración Reciclado de nutrientes 4. Cobertura de criptógamas Indicador de la estabilidad de superficie del suelo, de su resistencia a la erosión y disponibilidad de nutrientes Reciclaje de nutrientes 5. Grado de fragmentación de la costra física Mide la cantidad de costra superficial disponible para la erosión hídrica o eólica Estabilidad 6. Tipo y grado de la erosión del suelo Estima la naturaleza y severidad de los procesos erosivos actuales Estabilidad 7. Deposición de materiales (sedimentos) Evalúa la cantidad de depósitos aluviales Estabilidad 8. Rugosidad de la superficie del suelo Indicador de la rugosidad de la superficie del suelo en base a su capacidad de retener agua, sedimentos y semillas Infiltración Reciclado de nutrientes 9. Naturaleza del suelo (resistencia a la perturbación) Estima la probabilidad de perder suelo por una perturbación mecánica Estabilidad 10. Test de humectación Evalúa la estabilidad/dispersión de los agregados del suelo cuando está húmedo Estabilidad Infiltración 11. Textura Indicador de la capacidad de infiltración y almacenamiento de agua Infiltración Tabla 1. Identificación de los indicadores de la superficie del suelo y relación con los índices de estabilidad [“APLICACIÓN DEL MÉTODO LANDSCAPE FUNCTION ANALYSIS (LFA) PARA EVALUAR EL ESTADO DE LAS RESTAURACIONES MINERAS. EL CASO DE LA MINERÍA DE CARBÓN A CIELO ABIERTO EN UTRILLAS (TERUEL)”.] 2013 32 carácter excepcional. (Merino et al., 2009). A continuación se presenta un diagrama ombrotérmico de la zona (Figura 12): Figura 12. Diagrama ombrotérmico (Climodiagrama de Walter-Gausen) para los datos de estación de Montalbán (León, 1991). En cuanto al paisaje vegetal de la zona se caracteriza por un mosaico de cultivos de cereal en las zonas planas, terrazas de cultivo abandonadas y vegetación discontinua en laderas, dominada sobre todo por matorrales de Genista scorpius, Rosmarinus officinalis y Thymus vulgaris entre otros, así como pequeños espacios reforestados con Pinus halepensis. La vegetación potencial varía desde la dominancia de Quercus ilex hacia la de Quercus faginea (Monserrat, 1990), Geológicamente hablando, la zona se encuentra en el interior de la cuenca sedimentaria de Teruel. Los materiales carboníferos se encuentran compuestos por hulla subituminosa de edad Albiense. También decir que las escombreras mineras están formadas por estériles, los cuales se componen de arcillas, arenas, arcillas arenosas, alternando con areniscas pardoamarillentas, argilitas y limonitas. Los suelos de la región varían entre Typic y Lithic Xerorthent hacia Calcic Xerochrept (Arranz, 2004), con un pH neutro o moderadamente básico (7.3 a 8.6), generalmente pobres en materia orgánica (0.5-2.1%), de textura franco arcillosa y pedregosidad elevada en algunas ocasiones (20-47% en peso). Una vez explicado las condiciones climáticas, geológicas, edáficas etc. Vamos a profundizar en una serie de características ambientales que dificultan bastante el establecimiento y desarrollo de la cubierta vegetal y por lo tanto del éxito de la restauración de los ecosistemas [“APLICACIÓN DEL MÉTODO LANDSCAPE FUNCTION ANALYSIS (LFA) PARA EVALUAR EL ESTADO DE LAS RESTAURACIONES MINERAS. EL CASO DE LA MINERÍA DE CARBÓN A CIELO ABIERTO EN UTRILLAS (TERUEL)”.] 2013 33 restaurados ejecutados como consecuencia de la minería de carbón a cielo abierto en Utrillas (Teruel):  Los sustratos están formados por estériles procedentes de la actividad minera que poseen unas características físico-químicas adversas para el establecimiento de la vegetación.  Deficiencia en la llegada de propágulos debido a la lejanía de las fuentes productoras.  Diseño topográfico poco adecuado que en ocasiones aumenta la intensidad de la escorrentía notablemente y por lo tanto la erosión.  Sustratos impermeables que dificultan enormemente la infiltración y aumentan la escorrentía.  Pastoreo intensivo y no ordenado Además de estas características, hay que tener en cuenta las condiciones climáticas de la zona enunciadas anteriormente como la escasez de precipitaciones con episodios lluviosos intensos añadido a grandes valores de evapotranspiración. Cabe destacar que la erosividad climática se ve incrementada como hemos dicho antes por los sustratos poco permeables con poca materia orgánica, sumado a un diseño topográfico que favorece la entrada de escorrentía por la cabecera de la ladera ( Moreno de las Heras et al., 2008). 3.1.1 Descripción de las restauraciones mineras realizadas en el área de estudio Mina de “El Murciélago” La restauración que se realizó en esta mina (Figura 13) , pertenece al grupo de las de “Primera Generación”, basadas en el modelo “platarforma-talud-cuenca”, donde se dedicaban al cultivo de secano en las plataformas, aunque en este caso no ha sido así. Sobre los taludes se extendió tierra vegetal (20-50 cm) realizándose labores de revegetación con mezclas comerciales de semillas herbáceas (entre estas cabe destacar el uso de las siguientes especies: Lolium perenne, Medicago sativa, Poa pratensis, F. arundinacea…). El relieve abrupto de los taludes, pistas y cunetas limita notablemente el desarrollo de la vegetación y constituye una fuente de sedimentos que contaminan los cauces naturales (cubierta vegetal 15-25 %, riqueza de especies 10-15 y tasa de erosión 10-20 t/ha año; Moreno de las Heras et al., 2008). Este modelo no es sostenible a largo plazo. [“APLICACIÓN DEL MÉTODO LANDSCAPE FUNCTION ANALYSIS (LFA) PARA EVALUAR EL ESTADO DE LAS RESTAURACIONES MINERAS. EL CASO DE LA MINERÍA DE CARBÓN A CIELO ABIERTO EN UTRILLAS (TERUEL)”.] 2013 34 Figura 13. Fotografía de la antigua mina de “El Murciélago” “Antigua mina de “El Moral” La restauración de “El Moral” (Figura 16) representa una transición entre las de “Primera Generación” y las de “Segunda Generación” en cuencas, desarrolladas por la empresa MFUSA. Se caracteriza por un sistema de taludes de pendientes más suaves con un mejor desarrollo de la vegetación (Figura 15) con bermas y plataformas de cultivo intercaladas (Figura 14). Figura 14. Imagen de una plataforma aprovechada por cultivos en “El Moral” [“APLICACIÓN DEL MÉTODO LANDSCAPE FUNCTION ANALYSIS (LFA) PARA EVALUAR EL ESTADO DE LAS RESTAURACIONES MINERAS. EL CASO DE LA MINERÍA DE CARBÓN A CIELO ABIERTO EN UTRILLAS (TERUEL)”.] 2013 35 Figura 15. Imagen de una ladera de “El Moral” suavizada con buen desarrollo de vegetación La construcción de las laderas tenía las siguientes características:  Pendiente entre el 18-21 %  Aislamiento de la ladera de las escorrentías procedentes de las plataformas, pistas o bermas superiores.  Extendido de una capa de tierra vegetal de unos 50cm de espesor.  Labrado transversal a la pendiente. Labor profunda  Aplicación de purines .Labrado para evitar malos olores.  Siembra de mezcla de herbáceas a finales de febrero tras el frío invernal  Labor superficial para el enterramiento de las semillas  Dos años después se realiza la plantación de especies leñosas Figura 16. Foto de la antigua mina del “El Moral” [“APLICACIÓN DEL MÉTODO LANDSCAPE FUNCTION ANALYSIS (LFA) PARA EVALUAR EL ESTADO DE LAS RESTAURACIONES MINERAS. EL CASO DE LA MINERÍA DE CARBÓN A CIELO ABIERTO EN UTRILLAS (TERUEL)”.] 2013 36 3.1.2 Comportamiento ecohidrológico de las laderas restauradas En un estudio a escala de la provincia de Teruel se han identificado los principales factores que controlan la colonización vegetal de las escombreras mineras (Moreno de las Heras et al 2008):  En condiciones de sustratos ácidos la colonización vegetal se ve muy limitada, dado que la flora del entorno es basófila y no se puede desarrollar adecuadamente en este tipo de condiciones tan adversas.  Si nos centramos en la zona oriental de la provincia, donde está más influenciada por el clima mediterráneo, y donde los taludes fueron rellanados de tierra vegetal aunque no hubo revegetación, se ha producido una colonización vegetal a través de propágulos procedentes de vegetación natural.  En la zona occidental, con carácter más continental, el conjunto de herbáceas sembradas producen un impedimento para el establecimiento de las plantas del entorno. Aunque este se puede ver modificado si hay alguna perturbación tipo sobrepastoreo o fungicidas que generan huecos y una posible colonización. Como se ha comentado antes en el marco conceptual, unos de los principales factores limitantes es la erosión superficial provocada por la escorrentía. Este efecto de la erosión se ve incrementado por la impermeabilidad del sustrato y/o por las morfologías convexas o abruptas en la parte superior de los taludes, o por conexiones con flujos hídricos procedentes de pistas o bermas externas a las laderas. Estos flujos hídricos son muy negativos y perjudiciales si se concentran en regueros. Estos regueros (Figura 17) realizan una evacuación de la escorrentía que produce una disminución de la humedad edáfica en las zonas interreguero (que se encuentran entre dos regueros) (Figura 18) siendo un efecto restrictivo producido por la erosión. Según varias investigaciones realizadas sobre el efecto de la erosión en regueros sobre el desarrollo del suelo y la vegetación, se ha determinado que la escorrentía exógena procedente de la parte superior de los taludes produce una fuerte inestabilidad que se debe evitar. (Moreno de las Heras et al., 2008; Espigares et al., 2009). Figura 17. Foto de un reguero localizado en una ladera de estudio en el “Cerro del Moral” [“APLICACIÓN DEL MÉTODO LANDSCAPE FUNCTION ANALYSIS (LFA) PARA EVALUAR EL ESTADO DE LAS RESTAURACIONES MINERAS. EL CASO DE LA MINERÍA DE CARBÓN A CIELO ABIERTO EN UTRILLAS (TERUEL)”.] 2013 37 La intensificación del déficit hídrico propiciada por la erosión en regueros condiciona notablemente la colonización vegetal de las laderas. En concreto se ha comprobado que limita y dificulta la germinación, incorporación de la materia orgánica, la supervivencia, el desarrollo de estabilidad de agregados, mineralización y reciclado de nutrientes (N y P) así como el tamaño y la actividad de las poblaciones microbianas y la producción de semillas por parte de las plantas (Espigares et al 2009). Figura 18. Foto de una zona interreguero en una ladera de la zona de estudio de “El Moral” Según Moreno de las Heras (2009) se han identificado tres posibles direcciones de evolución de los sistemas de ladera restaurados: 1) Hacia la degradación, donde hay presencia y formación de regueros y una bajísima productividad vegetal con presencia de alfalfa. Este sería el caso de las laderas de ”Primera generación”. 2) Hacia la gradación, donde el sistema está formado por comunidades mixtas de matorral y gramíneas (dominadas por Genista scorpius). Este sería el caso de laderas de “Segunda generación”. 3) Situación intermedia donde hay incertidumbre en la evolución que va a seguir el sistema. Este tipo de dirección de evolución también correspondería a laderas de “Segunda generación”. [“APLICACIÓN DEL MÉTODO LANDSCAPE FUNCTION ANALYSIS (LFA) PARA EVALUAR EL ESTADO DE LAS RESTAURACIONES MINERAS. EL CASO DE LA MINERÍA DE CARBÓN A CIELO ABIERTO EN UTRILLAS (TERUEL)”.] 2013 38 A continuación se presenta la tabla 4 con datos indicativos de la trayectoria evolutiva del sistema de laderas restaurado: % de cobertura vegetal Tasa de erosión de regueros (t ha-1 año -1) Trayectoria de evolución del sistema <30 >20 Situación de degradación >50 <5 Situación de desarrollo vegetal 30-50 5-20 Cubierta vegetal discontinua en un patrón de matas y claros Tabla 4. Datos referentes a la trayectoria evolutiva del sistema de laderas restaurado en la zona de estudio. Espigares et al (2012) confirmaron que existe un flujo de escorrentía desde los claros hacia las matas, favoreciendo a las manchas de vegetación que disponen de más agua desarrollando una mayor biomasa y sufriendo menos estrés. Si las condiciones que se dan en el sistema son favorables, debido sobre todo a un menos volumen de escorrentía, se generan islas de fertilidad que controlan el balance hídrico, con unas ventajas muy importantes: Se establece en los claros y contribuye al cerramiento del dosel arbustivo-herbáceo, además de optimizar su control del uso de agua en la ladera. La escorrentía no es siempre la principal fuerza directriz de la sucesión en los ambientes continentales de Teruel. En laderas restauradas donde los recursos edáficos están estabilizados frente a la erosión, la distancia de las fuentes de propágulos o el efecto inhibidor de las herbáceas sembradas controlan la colonización vegetal. (Moreno de las Heras et al., 2008). 2.2 Protocolo de muestreo Para abordar los objetivos del TFG, aplicando el método LFA, se recogió información de seis laderas restauradas y de una ladera natural, tomada como referencia. La selección de las seis laderas restauradas responde a los siguientes criterios: a) Recoger un gradiente de degradación-agradación de los ecosistemas restaurados en los taludes, que represente el éxito-fracaso de las actuaciones. b) Aprovechar las laderas de “El Moral” de las que existe información ecohidrológica obtenida mediante una aproximación experimental (Merino et al. 2009). Para la ladera de referencia se seleccionó un encinar ubicado en la misma litología que los taludes (Formación Escucha), en la misma orientación (Norte) y similar pendiente. La identificación de cada ladera y la mina a la que pertenece se indican en la tabla 5. Igualmente se presentan las fotografías de las laderas (Figuras 19, 20, 21, 22, 23, 24, 25). [“APLICACIÓN DEL MÉTODO LANDSCAPE FUNCTION ANALYSIS (LFA) PARA EVALUAR EL ESTADO DE LAS RESTAURACIONES MINERAS. EL CASO DE LA MINERÍA DE CARBÓN A CIELO ABIERTO EN UTRILLAS (TERUEL)”.] 2013 39 Identificación Lugar L1 “El Moral” L2 “El Moral” L3 “El Moral” L4 “El Moral” L5 “El Murciélago” L6 “El Murciélago” L7(Ecosistema de Referencia) Carretera las “Parras de Martín) Tabla 5: Identificación de las laderas estudiadas La ladera 1 (L1) se caracteriza por presentar una densa red de regueros y una escasa cobertura vegetal. La ladera 2 (L2) presenta una red de regueros discontinua y un patrón de vegetación discontinuo con manchas de gramíneas remanentes de la siembra inicial en zonas de importación de agua y claros exportadores, colonizados por Santolina chamaecyparissus. La ladera 3 (L3) no presenta red de regueros, con un patrón de vegetación continuo, donde domina Genista Scorpius. La ladera 4 (L4) presenta una red de regueros discontinúa, pero con menor grado que L2, con un patrón de vegetación discontinuo con manchas de gramíneas remanentes de la siembra inicial en zonas de importación de agua y claros exportadores, colonizados por Brachypodium ramosum. La ladera 5 (L5) y ladera 6 (L6) presenta algún reguero discontinuo de forma ocasional, con mínimo establecimiento de la cubierta vegetal, predominado por suelo desnudo, donde las matas suelen ser de Santolina chamaecyparissus y Medicago Sativa. La ladera 7 (L7) ha sido explicada en el área de estudio. Figura 19. Fotos de la ladera 1 [“APLICACIÓN DEL MÉTODO LANDSCAPE FUNCTION ANALYSIS (LFA) PARA EVALUAR EL ESTADO DE LAS RESTAURACIONES MINERAS. EL CASO DE LA MINERÍA DE CARBÓN A CIELO ABIERTO EN UTRILLAS (TERUEL)”.] 2013 40 Figura 20. Fotos de la ladera 2 Figura 21. Fotos de la ladera 3 [“APLICACIÓN DEL MÉTODO LANDSCAPE FUNCTION ANALYSIS (LFA) PARA EVALUAR EL ESTADO DE LAS RESTAURACIONES MINERAS. EL CASO DE LA MINERÍA DE CARBÓN A CIELO ABIERTO EN UTRILLAS (TERUEL)”.] 2013 41 Figura 22. Fotos de la ladera 4 Figura 23. Fotos de la ladera 5 Figura 24. Fotos de ladera 6 [“APLICACIÓN DEL MÉTODO LANDSCAPE FUNCTION ANALYSIS (LFA) PARA EVALUAR EL ESTADO DE LAS RESTAURACIONES MINERAS. EL CASO DE LA MINERÍA DE CARBÓN A CIELO ABIERTO EN UTRILLAS (TERUEL)”.] 2013 48 (B) (C) Figura 28. Test Kruskal-Wallis. IE (A), II (B) e IR (C) [“APLICACIÓN DEL MÉTODO LANDSCAPE FUNCTION ANALYSIS (LFA) PARA EVALUAR EL ESTADO DE LAS RESTAURACIONES MINERAS. EL CASO DE LA MINERÍA DE CARBÓN A CIELO ABIERTO EN UTRILLAS (TERUEL)”.] 2013 49 Según los resultados de la figura 27 y 28, se deduce en cuanto al Índice de Estabilidad (%), la existencia de 3 grupos bien diferenciados que serían los siguientes:  Un primer grupo formado por la ladera L1, L2, L4 y L6 que no presentan diferencias significativas entre ellas. Aunque cabe destacar que L2 y L4 poseen valores muy similares y siendo superiores a L1 y L6.  Un segundo grupo formado por L3 y L7. Poseen los mejores valores de este índice, siendo las más estables y sin diferencias significativas entre ellas, siendo L3 un poco más estable que L7 (ecosistema de referencia).  Un tercer grupo formado únicamente por L5, que se diferencia significativamente de los grupos enunciados anteriormente, presentando el valor más bajo del Índice de Estabilidad. En cuanto al Índice de Infiltración (%), al contrario que el Índice de Estabilidad no hay presencia de tres grandes grupos diferenciados. Por una parte, L3 se diferencia de todas significativamente, siendo esta la que mayor valores de este índice presenta. Seguida de L7, que con la única ladera que no tiene diferencias significativas es con L4. Además L4 no tiene diferencias significativas con L2. Finalmente L5 se diferencia de todas significativamente presentado de los valores más bajos, mientras que L1 y L6 presentan el mismo comportamiento. Podríamos establecer que L3 posee la mayor capacidad de infiltración, seguido por L7, L2 y L4, que presentan un comportamiento frente al control de la escorrentía parecido. Le seguiría L5 y finalmente el último grupo con el peor comportamiento en cuanto este índice que son L1 y L6. Con el Índice de Reciclado de Nutrientes (%) se vuelven a formar 3 grupos bien diferenciados:  L3 y L7 formarían un primer grupo caracterizado por tener los valores más altos, y por lo tanto mejores.  Un segundo grupo formado por L2 y L4, con valores intermedios entre el último grupo que es el más deficiente.  El tercer grupo está compuesto por L1, L5 y L6, que presentan los valores más bajos. Haciendo un balance global de los 3 índices, L3 y L7 son las laderas cuyo estado funcional de su ecosistema se encuentra mejor, siendo ambas muy parecidas, incluso teniendo L3, un mejor estado que el propio ecosistema de referencia (L7) lo cual es muy significativo. Ya que L3 es una ladera que alcanzado el grado de comportamiento de un ecosistema de referencia local, y como hemos dicho incluso está en mejor estado. [“APLICACIÓN DEL MÉTODO LANDSCAPE FUNCTION ANALYSIS (LFA) PARA EVALUAR EL ESTADO DE LAS RESTAURACIONES MINERAS. EL CASO DE LA MINERÍA DE CARBÓN A CIELO ABIERTO EN UTRILLAS (TERUEL)”.] 2013 50 Los peores estados funcionales son los de las laderas L1, L5 y L6, que aunque L1 pertenezca al grupo de la restauraciones de segunda generación se comportan de forma muy similar a las de primera (L5 y L6) por lo que el resultado de la restauración de esta ladera L1 se aleja notablemente de las otras laderas de 2º generación. Y finalmente L2 y L4, están en esa situación intermedia, en que su estado funcional no es muy deficiente pero no alcanza los niveles de L3 y L7. 4.1.2 Caracterización del paisaje Una vez presentados los resultados del los índices a escala de ladera, se presentan a continuación los resultados procedentes de la caracterización de la organización del paisaje (Nº de parche/10 m; Área total de parches por transecto; Índice de organización del paisaje*; y Distancia media entre parches) con el fin determinar la función y estado de la salud de la vegetación, ya que la disposición espacial de esta, refleja el control que tiene sobre la escorrentía y los recursos. Además de ayudar a identificar en qué estado sucesional se encuentran las respectivas laderas. Para determinar las diferencias significativas estadísticas entre las laderas se ha realizado un test estadístico Kruskal-Wallis y un “Post hoc” de Bonferroni que aparece en la figura 29. (A) [“APLICACIÓN DEL MÉTODO LANDSCAPE FUNCTION ANALYSIS (LFA) PARA EVALUAR EL ESTADO DE LAS RESTAURACIONES MINERAS. EL CASO DE LA MINERÍA DE CARBÓN A CIELO ABIERTO EN UTRILLAS (TERUEL)”.] 2013 51 (B) (C) [“APLICACIÓN DEL MÉTODO LANDSCAPE FUNCTION ANALYSIS (LFA) PARA EVALUAR EL ESTADO DE LAS RESTAURACIONES MINERAS. EL CASO DE LA MINERÍA DE CARBÓN A CIELO ABIERTO EN UTRILLAS (TERUEL)”.] 2013 52 (D) Figura 29. Análisis de Kruskal-Wallis para las distintas laderas en base a los datos de organización del paisaje. Nº total de parches/10 m (A); Área total de parches (B); Distancia media entre parches (C); Índice de organización del paisaje (D) En cuanto a los resultados expuestos de las variables de la figura 29:  “Nº de parches /10m de transecto”: L1 se diferencia significativamente de todas la laderas menos L5 y L6, mientras que L2 (con el valor más alto) únicamente no tiene diferencias significativas con L3. Por otra parte L3 y L4 presentan diferencias significativas solo con L6. Cuanto mayor número de parches, mayor control de la escorrentía y erosión. Lo que no implica que los índices que muestran el estado funcional de las laderas sean mayores cuanto mayores sean el nº de parches. Hay que tener en cuenta el área de los parches, como ocurre en L3 y L7 (presentan las mayores áreas) que poseen los mayores valores de IE, II y IRN aunque en las variables de nº de parches/10m de transecto y distancia media entre parches no tengan lo valores más óptimos.  “Área total de parches”: La ladera 3 posee el valor más alto y se diferencia significativamente de todas menos de la ladera L2 y L7. La L2 y L7 respectivamente no se diferencian significativamente de las demás. Cuanta más área de parches haya mayor control habrá sobre la escorrentía y la erosión.  “Distancia media entre parches”: La ladera 1 se diferencia significativamente de todas las demás. Mientras hay otro grupo formado por la ladera L2, L3, L4, L5 y L7 que no [“APLICACIÓN DEL MÉTODO LANDSCAPE FUNCTION ANALYSIS (LFA) PARA EVALUAR EL ESTADO DE LAS RESTAURACIONES MINERAS. EL CASO DE LA MINERÍA DE CARBÓN A CIELO ABIERTO EN UTRILLAS (TERUEL)”.] 2013 53 poseen diferencias significativas, aunque L2, L3 y L7 poseen valores muy similares e inferiores a L5. Y finalmente L6 no tiene diferencias con L4, L5 y L7. Cuanta más distancia entre los parches menos control sobre la escorrentía y la erosión hay. Siendo L1, L5 y L6 las que menos distancias poseen, lo que implica bajos resultados en los índices (IE, II, IRN) como hemos observado antes. Mientras que el resto de laderas ocurre lo contrario.  “Índice de organización del paisaje”: Lo que nos indica esta variable es la cantidad de espacio que ocupa los parches respecto a la distancia del transecto. Es decir, una relación entre la longitud total de todos los parches respecto a longitud del transecto. L2, L3 y L7 no tienen diferencias significativas entre ellas, lo que indica un comportamiento parecido. Las laderas 5 y 6, se diferencias de todas las demás, siendo estas homogéneas. L2 se diferencia aparte de L5 y L6, de L1. Y L2, L4 y L7 no tienen diferencias significativas entre ellas. Más distancia que ocupan los parches mejor para el estado funcional de la ladera, ya que retendrá mayor recurso hídrico y se frenara la erosión. Comentados los resultados de las variables de la caracterización de la organización del paisaje, podemos establecer que L3 y L7, y seguida de L2 en algunas ocasiones, poseen los valores más óptimos, que nos indican que controlan mejor la escorrentía y erosión. L1, L5 y L6 poseen los peores, lo que se traduce en peores estados funcionales que se observan claramente en los resultados de los índices de Estabilidad, Infiltración y Reciclado de Nutrientes. 4.1.3 Validación del método LFA para la cuenca de Utrillas En cuanto la validación del LFA para las condiciones locales de la zona de estudio se ha realizado mediante un análisis de correlación de Spearman con parámetros medidos empíricamente por el grupo de investigación de los departamentos de ecología de la Universidad de Alcalá de Henares y Zaragoza. Para la correlación, de tanto el Índice de Infiltración como el de Estabilidad para las laderas L5, L6 y L7 (Ecosistema de referencia), no hay datos empíricos disponibles ya que no se han realizado aún la medida de las siguientes variables: Coeficiente Escorrentía (%), Tasa de erosión regueros (Tn/ha/año), Densidad de regueros (m/m2) y Tasa de erosión (Tn/ha/año). Tampoco hay datos de nitrógeno total y materia orgánica de L7 (Ecosistema de referencia). Una vez dicho esto, cabe resaltar que no supone ningún inconveniente a la hora de realizar la correlación, ya que L3 es muy similar a L7 y L5 y L6 muy parecidas a la L1, habiendo dos valores extremos: uno máximo correspondiente a L7 y otro mínimo correspondiente a L1, de tal forma, que la correlación no queda desvirtuada. A continuación en la figura 30 se presentan los resultados de la correlación de Spearman: [“APLICACIÓN DEL MÉTODO LANDSCAPE FUNCTION ANALYSIS (LFA) PARA EVALUAR EL ESTADO DE LAS RESTAURACIONES MINERAS. EL CASO DE LA MINERÍA DE CARBÓN A CIELO ABIERTO EN UTRILLAS (TERUEL)”.] 2013 54 (A) [“APLICACIÓN DEL MÉTODO LANDSCAPE FUNCTION ANALYSIS (LFA) PARA EVALUAR EL ESTADO DE LAS RESTAURACIONES MINERAS. EL CASO DE LA MINERÍA DE CARBÓN A CIELO ABIERTO EN UTRILLAS (TERUEL)”.] 2013 55 (B) (C) Figura 30. Relación entre el Índice de Estabilidad y la tasa de erosión, la densidad de regueros y la tasa de erosión en regueros (A). Relación entre el Índice de Infiltración y el coeficiente de escorrentía (B). Relación entre el Índice de Reciclado de Nutrientes y al contenido en materia orgánica y el contenido en Nitrógeno (C). [“APLICACIÓN DEL MÉTODO LANDSCAPE FUNCTION ANALYSIS (LFA) PARA EVALUAR EL ESTADO DE LAS RESTAURACIONES MINERAS. EL CASO DE LA MINERÍA DE CARBÓN A CIELO ABIERTO EN UTRILLAS (TERUEL)”.] 2013 56 Para la relación del Índice de estabilidad con la tasa de erosión (R2=0.75, p= 0.01) y tasa de erosión de regueros (R2=0.84, p= 0.01) la correlación es positiva alta, mientras que para la relación del Índice de Estabilidad con la densidad de regueros (R2=0.898, p=0.01) es muy alta. Por lo que podemos establecer que el método se ajusta a las condiciones locales del área de estudio debido a la alta relación observada entre el Índice de Estabilidad y los datos de validación con los que ha sido correlacionado. Para la relación entre el Índice de Infiltración y el Coeficiente de escorrentía (R2=0.856, p=0.05) la correlación es positiva alta, ajustándose a los datos de validación. Finalmente la validación del Índice de Reciclado Nutrientes correlacionándolo con el contenido de materia orgánica (R2= 0.877, p=0.05) y contenido de nitrógeno total (R2=0.921, p=0.05) se observa una correlación positiva alta y muy alta respectivamente. 4.2 Resultados a nivel de parche Lo primero se presenta el resultado de la identificación de la tipología de cada parche/interparche realizado en el muestro (Tabla 7) Presencia en ladera Tipo de cubierta de suelo Clave L1, L2, L3, L4, L5,L6 y LE Calvero o Suelo desnudo BARE L1, L2, L3, L4, L5,L6 y LE Roca ROC L1, L2, L3, L4 y LE Tapiz Herbáceo TAPH L2, L4 y LE Lolium perenne LOL L2, L4,L6 y LE Dactylis glomerata DAC L4 y LE Brachypodium ramosum BRACH L3, L4,L5 y LE Genista scorpius GEN L1,L2,L3,L4,L5,L6 y LE Santonina chamaecyparissus SAN L1,L3,L4 y LE Thymus vulgaris THY L5 y L6 Medicago sativa MED LE Quercus ilex QUER LE Poa bulbosa POA LE Rosal sp ROS L1 Pinus halepensis PIN LE Asociación Genista scorpius y Quercus ilex GEN-QUER Tabla 7. Identificación de los tipos de parches/interparches Cabe destacar que hay varios tipos de parches que no han sido utilizados en el análisis estadístico debido a su poca presencia (n<3) como son: Spartium junceum, Juniperus communis, Erinacea anthyllis, Carduus acanthoides y una asociación producida entre Genista scorpius y Juniperus communis. [“APLICACIÓN DEL MÉTODO LANDSCAPE FUNCTION ANALYSIS (LFA) PARA EVALUAR EL ESTADO DE LAS RESTAURACIONES MINERAS. EL CASO DE LA MINERÍA DE CARBÓN A CIELO ABIERTO EN UTRILLAS (TERUEL)”.] 2013 57 4.2.1 Contribución de cada tipo parche/interparche a los índices Para observar el comportamiento individual que tiene cada tipo parche/interparche, se presentan, por una parte el resultado de cada uno de los índices (IE, II e IRN) que proporciona cada tipo de parche/interparche y la contribución en (%) de cada tipo de parche, es decir, la cantidad con la que contribuye cada parche al resultado de los índices, que dependerá de la presencia o la distancia y el espacio que ocupan cada parche con respecto a la distancia total del transecto. En la figura 31, se presentan los resultados de la contribución (%) de cada tipo de parche/interparche a los índices. (A) (B) [“APLICACIÓN DEL MÉTODO LANDSCAPE FUNCTION ANALYSIS (LFA) PARA EVALUAR EL ESTADO DE LAS RESTAURACIONES MINERAS. EL CASO DE LA MINERÍA DE CARBÓN A CIELO ABIERTO EN UTRILLAS (TERUEL)”.] 2013 64 Según este modelo hay 3 posibles direcciones de evolución de los sistemas de ladera restaurados, como se ha identificado igualmente por el método LFA y concuerdan con los resultados del propio modelo que fue desarrollado en el área de estudio.  Las laderas L1, L5 y L6 pertenecen al tipo 1 según el modelo, es decir, su dirección de evolución es hacia la degradación, ya que hay presencia de regueros y pobres establecimientos de parches de vegetación por lo que hay un escaso control de la escorrentía y grandes tasas de erosión que impiden el establecimiento adecuado de la vegetación, como queda reflejado en estas laderas.  Las laderas L3 y L7, pertenecen al tipo 3, siguen un camino hacia la agradación, formado por comunidades mixtas de matorral y gramíneas, donde domina Genista scorpius. Estas laderas se caracterizan por ausencia de regueros, con presencia de importantes parches que controlan la escorrentía de forma adecuada permitiendo esta evolución positiva del sistema. En cuanto a L3 cabe destacar que posee valores incluso un poco mejores que el propio ecosistema de referencia, lo cual es muy significativo, ya que una restauración minera ya ha alcanzado un estado funcional adecuado para la zona de estudio y se puede tomar de guía o ejemplo para futuras restauraciones.  Y finalmente L2 y L4, son del tipo 2, donde hay incertidumbre sobre la evolución que va a seguir el sistema, ya que se trata de una situación intermedia entre el tipo 1 y 3. Hay presencia de parches bien establecidos con presencia de regueros. La dirección que sigan estos sistemas dependerá de el control que tengan los parches sobre la cantidad de escorrentía, y por tanto que haya un control del sistema sobre la erosión y no al revés, ya que esto último produciría una disminución de la cobertura vegetal que en un futuro se vería incapaz de controlar la erosión y alcanzaría un estado hacia la degradación. Cabe destacar que estos umbrales de desarrollo de la restauración a escala de ladera, además de basarse en la tipología de laderas de Moreno de las Heras, funcionan según el modelo enunciado anteriormente en la introducción, el TTRP (figura 1), base conceptual del método LFA, ya que el desencadenante o trigger (lluvia), genera una escorrentía superficial que realiza una transferencia de los recursos a través del sistema, los cuales son retenidos por los parches, pudiéndose generar un pulso en la productividad vegetal. Es decir, al haber una presencia de parches estables que controlan esta escorrentía que posee los recursos de forma adecuada, como sucede en las laderas de tipo 3, permite que haya buenas productividades vegetales controlando la erosión y evitando las pérdidas de estos recursos y siguiendo el sistema una dirección de evolución hacia la gradación. En cambio, en las laderas de tipo 1, no se controlan esos recursos, que fluyen a través del sistema y se pierden debido la mínima presencia de parches, formándose regueros. Ello limita notablemente que algún parche por pequeño que sea pueda almacenar este flujo hídrico. Además este efecto es imposible de mitigar si no hay una presencia de unos parches adecuados, por lo que estos sistemas llevan caminos hacia la degradación. [“APLICACIÓN DEL MÉTODO LANDSCAPE FUNCTION ANALYSIS (LFA) PARA EVALUAR EL ESTADO DE LAS RESTAURACIONES MINERAS. EL CASO DE LA MINERÍA DE CARBÓN A CIELO ABIERTO EN UTRILLAS (TERUEL)”.] 2013 65 Esto se ve reflejado también en los resultados obtenidos en la figura 29, relacionados con las variables de organización del paisaje (Nº de parche/10 m; Área total de parches por transecto; Índice de organización del paisaje*; y Distancia media entre parches). Las laderas que poseen el mejor estado funcional L3 y L7, seguidas de L2 y L4, poseen los mejores valores en el computo global de todas estas variables: nº de parches por transecto, área de parches, longitud de parches/longitud de transecto, lo que facilita funcionar según el modelo TTRP (figura 1). En cambio las laderas L1, L5 y L6, con los peores valores de los índices, controlan peor la escorrentía, de manera que no operan los procesos TTRP. 5.2 Escala de parche En cuanto a los resultados obtenidos a nivel de parche, en la figura 31, se han presentado los datos de contribución de cada tipo de parche a los índices del LFA que caracterizan el estado funcional del sistema. En las laderas 1, 5 y 6 (las que se encuentran en peor estado), el interparche BARE (suelo desnudo), de carácter exportador de agua y nutrientes, es el más abundante y hegemónico. En estas laderas la escorrentía generada aquí no es interceptada por parches importadores como en las laderas que funcionan según el TTRP. Las laderas 2 y 4 de carácter intermedio; presentan altos valores del parche BARE (Suelo desnudo) y de SAN (Santonina chamaecyparissus), ambos de carácter exportador. Pero también de parches interceptores, aunque aparecen también otros parches que poseen niveles significativos de contribución como TAPH (Tapiz herbáceo) y LOL (Lolium perenne); este último remanente de las siembras de vegetación que realizó MFU. Se ha demostrado que estos dos últimos parches herbáceos constituyen islas de fertilidad hidrológicamente activadas (Espigares et al., 2012). La evolución de estas laderas se considera incierta. Probablemente dependerá de la magnitud de las perturbaciones que reciba en los próximos años: tormentas erosivas intensas, sequías acusadas o sobre-pastoreo podrían llevar a la degradación. La entrada de propágulos y/o introducción de plantones podría llevarlo a la agradación. Las laderas con un mejor estado funcional (3 y 7) presentan una distribución más homogénea de los parches. BARE (Suelo desnudo) y TAPH (Tapiz herbáceo) son los más abundantes, aunque en el caso de la ladera 3, el parche GEN (Genista scorpius), alcanza el mayor valor. Ello hace que sean las laderas con mayor control biológico de los recursos hídricos. Además en la ladera 7, encinar de referencia, el parche GEN-QUER (Genista scorpius – Quercus ilex) adquiere bastante importancia, lo que de la una resilencia frente a las perturbaciones (fuego por ejemplo). Del resultado obtenido para el parche GEN (Genista scorpius), podemos deducir, que esta especie es interesante en las restauraciones mineras. Ya que subsana la dificultad para la colonización y establecimiento de especies arbóreas y arbustivas, además de realizar funciones de especie nodriza, al proteger a otras especies como se ve reflejado en parches como GEN- [“APLICACIÓN DEL MÉTODO LANDSCAPE FUNCTION ANALYSIS (LFA) PARA EVALUAR EL ESTADO DE LAS RESTAURACIONES MINERAS. EL CASO DE LA MINERÍA DE CARBÓN A CIELO ABIERTO EN UTRILLAS (TERUEL)”.] 2013 66 QUER (Genista scorpius - Quercus ilex) o GEN-JUN (Genista Scorpius-Junniperus communis), observado durante la identificación de la tipología de parches durante el muestreo. Otros resultados obtenidos a escala de parche son los de la figura 33, relacionados con la caracterización de cada tipo de parche en función de los índices que caracterizan el estado funcional del sistema. Se observa que se cumple un patrón para los 3 índices, en que hay 3 grupos que se comportan de forma diferente:  Un primer grupo formado por el parche MED (Medicago Sativa), ROC (Roca), BARE (Suelo desnudo) y DAC (Dactylis glomerata). Este grupo tiene los peores valores en los 3 índices, ya que tienen un mal comportamiento en cuanto a la capacidad de retener el agua de escorrentía, que se traduce en un menor crecimiento vegetal lo que conlleva a una menor contribución a la hora del reciclaje de nutrientes que finalmente disminuye la estabilidad que confiere el parche al sistema.  Un segundo grupo formado por GEN (Genista scorpius) y GEN-QUER (Genista scorpiusQuercus ilex), que poseen los mejores comportamientos en los índices, como se ve reflejado en los resultados. Al contrario que el primer grupo poseen una alta capacidad para la retención de agua que mejora los procesos que se producen según el modelo TTRP.  Un último tercer grupo formado por BRA (Brachypodium ramosum), LOL (Lolium perenne), SAN (Santolina chamaecyparissus), TAPH (Tapiz herbáceo) y THY (Thymus vulgaris), que poseen un comportamiento intermedio entre los dos grupos anteriores. La comparación cualitativa entre la caracterización funcional de los parches en base a los índices del LFA con la caracterización de unidades eco-hidrológicas identificadas por Merino (2008) muestra una elevada similitud entre ambas. Podemos decir, que los distintos tipos de parche identificados muestran que se produce un ordenamiento en función de los valores de los índices de Estabilidad, Infiltración y Reciclado de Nutrientes que se corresponde, en gran medida, con el comportamiento hidrológico descubierto por el grupo de investigación en Ecohidrología de áreas restauradas y Restauración de minería a cielo abierto de la universidad de Alcalá para las unidades ecohidrológicas (figura 35) (Merino, 2008) que se presenta a continuación: [“APLICACIÓN DEL MÉTODO LANDSCAPE FUNCTION ANALYSIS (LFA) PARA EVALUAR EL ESTADO DE LAS RESTAURACIONES MINERAS. EL CASO DE LA MINERÍA DE CARBÓN A CIELO ABIERTO EN UTRILLAS (TERUEL)”.] 2013 67 Figura 35. Resultado del test de Kruskal-Wallis para determinar el comportamiento hidrológico de las comunidades vegetales. (Tomado de Merino, 2008) En esta gráfica de Merino (2008) identifica 7 unidades ecohidrológicas, que aunque no es el mismo número realizado en análisis Kruskal-Wallis que se ha realizado en este trabajo (11), nos da una gran idea de que el resultado a escala de parche obtenido por LFA, es válido al igual que a escala de ladera, ya que se ha comparado con datos empíricos de investigaciones. Este comportamiento ecohidrológico de cada tipo de parche establecido por Merino, y por tanto los parches identificados en este trabajo funcionan a través del modelo TTRP. Ya que hay parches como MED (Medicago Sativa), DAC (Dactylis glomerata), ROC (Roca) y un interparche BARE (Suelo desnudo) que son generadores de escorrentía cuando la lluvia actúa de desencadenante, mientras que GEN (Genista scorpius) y GEN-QUER (Genista scorpius-Quercus ilex), actúan como sumideros, reteniendo esta escorrentía, lo que les convierte en parches fértiles, produciéndose una distribución desigual entre los generadores y sumideros que produce una diferenciación en las características (infiltración y disponibilidad de agua, de materia orgánica, de propágulos, etc.). El resto de parches como SAN(Santolina chamaecyparissus), THY (Thymus vulgaris), LOL (Lolium perenne), BRA (Brachypodium ramosum) y el interparche TAPH se pueden comportar según las circunstancias y condiciones como sumideros o generadores. [“APLICACIÓN DEL MÉTODO LANDSCAPE FUNCTION ANALYSIS (LFA) PARA EVALUAR EL ESTADO DE LAS RESTAURACIONES MINERAS. EL CASO DE LA MINERÍA DE CARBÓN A CIELO ABIERTO EN UTRILLAS (TERUEL)”.] 2013 68 6.Conclusiones Los resultados del trabajo han permitido alcanzar las siguientes conclusiones:  El método LFA permite realizar la caracterización del estado funcional en la que se encuentran las restauraciones mineras de carbón a cielo abierto para las condiciones locales de la zona de estudio en Utrillas (Teruel). Por lo que se recomienda utilizarla como herramienta para la evaluación de los resultados obtenidos de las restauraciones.  Todos los índices de la caracterización del estado funcional de las laderas presentan una correlación entre alta y muy alta con las variables cuantitativas utilizadas para la verificación del método para las condiciones locales.  La caracterización de los tipos de parches realizada mediante el LFA ha permitido identificar diferentes comportamientos ecohidrológicos que se corresponden bien con los establecidos por Merino (2008).  Se refuerza la aplicabilidad y efectividad del método después de los resultados obtenidos en este TFG, así como la confirmación de que se trata de un método rápido y práctico, que supone un ahorro en esfuerzo y tiempo.  La estructura y dinámica de la vegetación en claros y matas es clave en el comportamiento de las laderas mineras estudiadas. Cuanto mayor sea el número de parches y la superficie que ocupan se dispondrá de un mayor control de la escorrentía y de los efectos erosivos que son la principal casusa del descenso de la calidad del estado funcional.  El método LFA ha identificado las laderas con los estados sucesionales propuestos por Moreno de las Heras (2010): Las laderas bien colonizadas por aliaga (Genista scorpius) se asemejan funcionalmente al ecosistema de referencia (encinar intervenido); en situación intermedia con incertidumbre sobre su futura evolución se encuentran laderas con regueros discontinuos y parches exportadores e importadores equilibrados; en estado de degradación se han identificado laderas reguerizadas con predominio de parches exportadores.  “Genista scorpius” se muestra como una especie interesante a la hora de utilizarla en las restauraciones mineras, ya que puede ser gran parte de la solución al proceso sucesional de este tipo de ecosistemas, minimizando la dificultad para la colonización y ayudando al establecimiento de especies arbóreas y arbustivas como “Quercus ilex” y “Juniperus communis” [“APLICACIÓN DEL MÉTODO LANDSCAPE FUNCTION ANALYSIS (LFA) PARA EVALUAR EL ESTADO DE LAS RESTAURACIONES MINERAS. EL CASO DE LA MINERÍA DE CARBÓN A CIELO ABIERTO EN UTRILLAS (TERUEL)”.] 2013 69 7.Bibliografía David Tongway. (2005). Landscape function analysis: A systems approach to assessing rangeland condition. David Tongway. (2011). Using landscape function principles to design and monitor active restoration. Moreno de las Heras, M. 2008. Efectos ecológicos de la erosión en laderas derivadas de la minería del carbón a cielo abierto. Tesis Doctoral. Universidad de Alcalá. Departamento de Ecología. Martín, R. Aguilar and Osvaldo E. Sala. 1999. Patch strucure, dynamics and implications for the functioning of arid ecosistems. Tree vol 14, pp 273-277 Rivas-Martínez S. (1987). Memoria del mapa de series de vegetación de España. Ministerio de Agricultura, Pesca y Alimentación. Madrid Espigares, T., Moreno-de las Heras, M., Nicolau, J.M., 2009. Performance of vegetation in reclaimed slopes affected by soil erosion. Restoration Ecology. José Manuel Nicolau. (2010). In Avelino García Álvarez (Ed.), Plan director para la restauración edafopaisajística y recuperación ambiental de los espacios degradados por la minería de carbón en la comarca de bierzo. José Manuel Nicolau Ibarra, José Francisco Martín Duque, Silvia Pérez Domingo, Francisco Comín Sebastián, Cristina Martínez Moreno, Ana Lucía Vela. La reconstrucción del relieve en la restauración de la minería a cielo abierto (2º Curso Internacional de Restauración Ecológica de Áreas Degradadas por la Minería a Cielo Abierto). José Manuel Nicolau, Tíscar Espigares, Mariano Moreno de las Heras y Luis Merino-Martín. (2011). Ecohidrología: Erosión hídrica y dinámica de la vegetación en laderas artificiales. Restauración ecológica de áreas afectadas por infraestructuras de transporte. (2011). Juan Pedro Rodríguez-López, Nieves Meléndez, Ana Rosa Soria, Poppe L. de Boer. (2009). Reinterpretación estratigráfica y sedimentológica de las formaciones Escucha y Utrillas de la cordillera ibérica. 22(Revista de la Sociedad Geológica) L. Merino-Martín, M. Moreno-de las Heras, S. Pérez-Domingo1, T. Espigares, and J. M. Nicolau. (2012). Hydrological heterogeneity in mediterranean reclaimed slopes: Runoff and sediment yield at the patch and slope scales along a gradient of overland flow. Hydrol. Earth Syst. Sci., 16, 1305–1320, 2012. Merino, L., Nicolau, J.M, Espigares, M.T., y De Las Heras, M., 2009. Dinámica de la vegetación y de los procesos hidrológicos en laderas restauradas de la minería del carbón a cielo abierto en Teruel. (Informe de proyecto para IET no publicado) Luis Merino-Martín, David D. Breshears, Mariano Moreno-de las Heras, Juan Camilo Villegas, Silvia Pérez-Domingo, Tíscar Espigares, and José M. Nicolau. (2012). Ecohydrological source- [“APLICACIÓN DEL MÉTODO LANDSCAPE FUNCTION ANALYSIS (LFA) PARA EVALUAR EL ESTADO DE LAS RESTAURACIONES MINERAS. EL CASO DE LA MINERÍA DE CARBÓN A CIELO ABIERTO EN UTRILLAS (TERUEL)”.] 2013 70 sink interrelationships between vegetation patches and soil hydrological properties along a disturbance gradient reveal a restoration threshold. 20 (Restoration Ecology). M. Moreno de las Heras, & J. M. Nicolau y T. Espigares. Interacción entre la erosión en regueros, contenido de humedad edáfica y colonización vegetal en laderas restauradas de la minería a cielo abierto del carbón en ambiente mediterráneo-continental (Teruel). Estudios De La Zona no Saturada Del Suelo Vol. VII. F.J. Samper Calvete y A. Paz González, 2005. M. Moreno-de las Heras, T. Espigares, L. Merino-Martín, J.M. Nicolau. (2011). Water-related ecological impacts of rill erosion processes in mediterranean-dry reclaimed slopes.84 (114124)(Catena). Mª Remedios Zamora Roselló. (2012). La transformación del territorio a través de la restauración de los espacios afectados por actividades extractivas. (Congreso Nacional del Medio Ambiente) Marcos Mejuto Mendieta. (2011). Afectación de la minería del carbón en las propiedades físicas y químicas de los suelos de la cuenca hidrográfica del río Rodrigatos (el bierzo, león) Maria C. Ruiz-Jaen and T. Mitchell Aide. (2005). Restoration success: How is it being measured? Restoration Ecology Vol. 13, no. 3, Pp. 569–577. Mariano Moreno‐de las Heras, José M. Nicolau, Luis Merino‐Martín, and Bradford P. Wilcox. (2010). Plot‐scale effects on runoff and erosion along a slope degradation gradient.46 (WATER RESOURCES RESEARCH) Real Decreto 975/2009, De 12 De Junio, Sobre Gestión De Los Residuos De Las Industrias Extractivas y De Protección y Rehabilitación Del Espacio Afectado Por Actividades Mineras. (2009). Tíscar Espigares, Mariano Moreno-de las Heras, and José M. Nicolau. (2011). Performance of vegetation in reclaimed slopes affected by soil erosion. Restoration Ecology. Tíscar Espigares, Luis Merino-Martín, Mariano Moreno-de las Heras and José-Manuel Nicolau. (2012). Intensity of ecohydrological interactions in reclaimed mediterranean slopes: Effects of run-off redistribution on plant performance. Published Online in Wiley Online Library (Wileyonlinelibrary.Com) DOI: 10.1002/eco.1307, ECOHYDROLOGY . VICENÇ CARABASSA – MARC VIZCANO– ESTEVE SERRA – ORIOL ORTIZ –JOSEP MARIA ALCAÑIZ. (2010). Getting real: A procedure for self-evaluation of quarry restorations. Avignon, France, 23-27/08/2010, (Proceedings 7th European Conference on Ecological Restoration). Vicente Manzano. (2006). Análisis de datos en psicología I. Guía práctica para el trabajo con SPSS. Montserrat, P. 1990. Vegetación. En: Mapa forestal del España: Daroca (ed. de la Torre R). Ministerio de Agricultura, Pesca y Alimentación, ICONA Madrid, pp 47-68. Tongway, David and Hindley, Norman. (1995) Assessment of Soil Condition of Tropical Grasslands. CSIRO Division of Wildlife and Ecology, Canberra [“APLICACIÓN DEL MÉTODO LANDSCAPE FUNCTION ANALYSIS (LFA) PARA EVALUAR EL ESTADO DE LAS RESTAURACIONES MINERAS. EL CASO DE LA MINERÍA DE CARBÓN A CIELO ABIERTO EN UTRILLAS (TERUEL)”.] 2013 71 Barberá, G. C., Navarro-Cano, J.A. y Castillo, V.M. 2006. Seedling recruitment in a semiarid steppe: The role of microsite and post-dispersal seed predation. Journal of Arid Environments 67 (2006) 701-714. Beeby, A. 1993. Applying Ecology. Chapman & Hall, London, 441 pp. Contreras, S. y Solé-Benet, A. 2003. Hidrofobia en suelos mediterráneos semiáridos: Implicaciones hidrológicas para una pequeña cuenca experimental en el este ibérico. Revista C&G. 17 (1-2), pp. 19-45. García-Fayos, P. 2004. Interacciones entre la vegetación y la erosión hídrica. En: Valladares, F (Ed.), Ecología del Bosque Mediterráneo en un Mundo Cambiante. Ministerio de Medio Ambiente, EGRAF, S.A., Madrid, pp. 309-334. Gómez-Aparicio, L. et al. 2004. Appliying plant facillitation to forest restoration: A meta analysis of the use of shrubs as nurse plants. Ecological Applications, 14 (4), 2004, pp. 11281138. Gómez Sal, A. y Nicolau, J.M. 1999. Hacia una caracterización ecológica de los paisajes degradados. Bases para su restauración. En: González Alonso, S. (ed) Libro homenaje a D. Ángel Ramos: 629-639. E.T.S.I. Montes. Madrid. Gobierno de Aragón, Manual de restauración de explotaciones mineras a cielo abierto de Aragón. Gutiérrez, J.R. y Squeo, F.A. 2004. Importancia de los arbustos en los ecosistemas áridos de Chile. Ecosistemas 13 (1): 36-45. Kapolka, N.M., Dollholpf, D.J. 2001. Effect of slope gradient and plant growth on soil loss on reconstructed steep slopes. International Journal of Surface Mining, eclamation and Environment, 15: 86-99. Ludwig, J.A., Wilcox, B.P., Breshears, D.D., Tongway D.J. and Imerson A.C. 2005. Vegetation patches and runoff-erosion as interacting ecohydrological processes in semiarid landscapes. Ecology, 86. Pp 288-297. Marañón, T., Camarero, J.J., Castro, J., Díaz, M., Espelta, J.M., Hampe, A., Jordano, P,Valladares, F., Verdú, M., Zamora, R., 2004. Heterogeneidad ambiental y nicho de regeneración. En: Valladares, F. 2004. Ecología del bosque mediterráneo en un mundo cambiante. Páginas 69-99. Ministerio de Medio Ambiente, EGRAF, S. A., Madrid. Martín, R. Aguilar and Osvaldo E. Sala. 1999. Patch structure, dynamics and implications for the functioning of arid ecosistems. Tree vol 14, pp 273-277. Palmer, J. 1992. Nutrient Cycling: Key to Reclamation Success. In: Chambers, J. & Wade, J. (eds): Evaluating Reclamation Success: The Ecological Consideration: 27-36. USDA, General Technical Report EN-164. [“APLICACIÓN DEL MÉTODO LANDSCAPE FUNCTION ANALYSIS (LFA) PARA EVALUAR EL ESTADO DE LAS RESTAURACIONES MINERAS. EL CASO DE LA MINERÍA DE CARBÓN A CIELO ABIERTO EN UTRILLAS (TERUEL)”.] 2013 72 Maestre, F.T. 2003 (b). Variaciones en el patrón especial a pequeña escala de los componentes de la costra biológica en un ecosistema mediterráneo semiárido. Revista Chilena de Historia Natural. No. 76, pp. 35-46 Tongway, D.J., Cortina, J. y Maestre F.T., 2004 (b). Heterogeneidad espacial y gestión de medios semiáridos. En Revista Ecosistemas 13 (1). Enero 2004. Pp 2-7. Zamora, R., Castro, J., Gómez, J. M., García, D., Hódar, J. A., Gómez, L. y Baraza, E. 2001. El papel de los matorrales en la regeneración forestal. Quercus 187: 40-47. Noy-Meir I. (1973) Desert Ecosystems: Environment and Producers. Annual Review of Ecology and Systematics, 4, 25-51 Haigh, M.J. 2000. The aims of land reclamation. In Haigh, M.J. (ed) Reclaimed land. Erosion control, soils and ecology. Pp: 1-20 [“APLICACIÓN DEL MÉTODO LANDSCAPE FUNCTION ANALYSIS (LFA) PARA EVALUAR EL ESTADO DE LAS RESTAURACIONES MINERAS. EL CASO DE LA MINERÍA DE CARBÓN A CIELO ABIERTO EN UTRILLAS (TERUEL)”.] 2013 73 Anexos [“APLICACIÓN DEL MÉTODO LANDSCAPE FUNCTION ANALYSIS (LFA) PARA EVALUAR EL ESTADO DE LAS RESTAURACIONES MINERAS. EL CASO DE LA MINERÍA DE CARBÓN A CIELO ABIERTO EN UTRILLAS (TERUEL)”.] 2013 80 8.2.2 Disposiciones específicas medioambientales y particulares de aplicación a la Minería Una vez normalizada la situación política del país, en la década de los 80, y con el horizonte de nuestra integración en la Comunidad Europea, se toma conciencia de la importante necesidad de paliar los impactos ambientales de todo tipo y en especial aquellos que son susceptibles de influir en la competencia y el desarrollo económico, como es en el caso de la minería. Debido a esta situación, se empiezan a emitir diversos decretos y normativas que van mejorando la regulación de los aspectos procedimentales para la conseguir los objetivos adecuados en restauración ambiental. En principio son normativas emanadas de la Administración Central, pero siempre con mención y reserva de su aplicación a las posibles competencias que puedan ser asumidas por las Comunidades Autónomas en proceso de creación. A pesar de esta dirección en líneas positivas que toma la normativa, no queda claro el posible carácter retroactivo de estas normas y su aplicación en actividades anteriores a la entrada en vigor de esta legislación, por lo que hay centenares de explotaciones, la mayoría inactivas en la actualidad, que se encuentran en un vació legal en cuanto a la obligación de efectuar la restauración e integración ambiental. En todo caso se deberán tener en cuenta estos procedimientos para su aplicación, si se considera procedente, antes de la ejecución de la oportuna restauración ambiental en las escombreras abandonadas. En otros aspectos, como es el caso de la responsabilidad de los titulares de la concesión, resulta complicado, cuando no inviable, la aplicación de modo estricto de la legalidad vigente en estas materias. A continuación se hace mención expresa a los Reales Decretos en el ámbito de la restauración para ser tenidos en cuenta ante una eventual actuación en los espacios degradados: - Real Decreto 2994/1982, de 15 de octubre, sobre Restauración del Espacio Natural Afectado por Actividades Mineras. (BOE nº 274 de 15 de noviembre de 1982, Mº de Industria y Energía) Uno de los problemas ambientales causados en ocasiones por la minería es el del deterioro de los terrenos circundantes a la zona de actividad, circunstancia que se manifiesta de modo especial en las explotaciones a cielo abierto, pudiéndose provocar perjuicios, no sólo de orden estético, sino también geomorfológico, como la erosión. La Ley de Minas de mil novecientos setenta y tres, se halla inmersa en una filosofía conservadora del medio ambiente, en la idea de que “la obtención de un recurso natural, como es el producto minero, sólo debe comprometer en la menor cuantía posible la utilización y conservación de otros bienes, como el espacio en el que se sitúan las explotaciones, procurando, al mismo tiempo, que las legítimas medidas de protección de dichos bienes han de evitar ser excesivamente maximalistas, de tal forma que no hagan económica o técnicamente inviable el desarrollo de las actividades extractivas”, ya que ello supondría un importante [“APLICACIÓN DEL MÉTODO LANDSCAPE FUNCTION ANALYSIS (LFA) PARA EVALUAR EL ESTADO DE LAS RESTAURACIONES MINERAS. EL CASO DE LA MINERÍA DE CARBÓN A CIELO ABIERTO EN UTRILLAS (TERUEL)”.] 2013 81 deterioro social y económico que el país no puede permitirse. La necesidad de guardar el preciso equilibrio entre los dos fines indicados, obliga al estudio particular de cada uno, con objeto de ponderar las numerosas y muy diversas circunstancias -algunas de ellas difícilmente posibles de preverque concurren en cada explotación y los requerimientos que se derivan de las características de su entorno natural, que presentan grandes diferencias de un lugar a otro. Por ello, el presente Real Decreto configura un sistema mediante el cual, en primer lugar, el titular de una solicitud de las previstas en la Ley de Minas, debe presentar un Plan de Restauración del Espacio Natural, afectado por las labores. El Plan tiene dos partes: - La primera dedicada a suministrar información sobre la descripción del lugar previsto para las labores mineras y su entorno, con información acerca del medio socioeconómico, ya que todo ello es necesario para ponderar la mayor o menor intensidad del Plan. - La segunda parte, contiene el proyecto de restauración propiamente dicho, incluyendo las medidas previstas para la protección del paisaje, acondicionamiento de la superficie del terreno, prevención de la erosión y otros. El Plan, una vez aprobado por la Administración, se convierte en obligatorio para el titular del derecho minero, quien puede ejecutarlo por sí o confiar la realización a la Administración, mediante la entrega de una cantidad periódica, con la cual aquélla dota un fondo destinado al efecto. Con ello se otorga flexibilidad al sistema, ya que en muchas ocasiones el titular del aprovechamiento carece de las posibilidades técnicas para acometer con garantía la realización del Plan. Otras veces la restauración sólo es posible una vez finalizada la explotación, por lo que sería muy difícil conseguir que su titular emprenda aquélla, debiendo ser pues la Administración la responsable de la ejecución del Plan con las cantidades periódicamente obtenidas. Por lo que se refiere a las explotaciones en marcha, dispone el Real Decreto que sus titulares presenten, cuando sean requeridos para ello, un proyecto de restauración. Entre los aspectos más destacables en este apartado cabría señalar: - Como se puede constatar, este Real Decreto, a pesar de tener casi tres décadas, mantiene su plena vigencia, no solo jurídicamente, sino también en sus aspectos técnicos y aplicados, siendo un claro referente para acometer las tareas de restauración ambiental. En todo caso es ya más cuestionable si se aplica con eficacia o no se obliga con la debida autoridad a los concesionarios y si no se hace con el debido rigor o, por otro lado, si las técnicas de revegetación, las especies y labores de cultivo que se usan para acometer estas restauraciones son las más adecuadas o se debe investigar y mejorar todavía mucho más. [“APLICACIÓN DEL MÉTODO LANDSCAPE FUNCTION ANALYSIS (LFA) PARA EVALUAR EL ESTADO DE LAS RESTAURACIONES MINERAS. EL CASO DE LA MINERÍA DE CARBÓN A CIELO ABIERTO EN UTRILLAS (TERUEL)”.] 2013 82 8.2.3 Normativa específica para la restauración en explotaciones de carbón a cielo abierto En 1984, con la inminente creación de las Comunidades Autónomas y el horizonte de la plena incorporación al concierto de las Naciones Europeas, se hace necesaria una re-adaptación de las normativas, dándole mayor concreción y regulando de modo especifico las restauraciones ambientales de la minería del carbón a cielo abierto. Así aparece el Real Decreto 1116/1984, de 9 de mayo, sobre Restauración del Espacio Natural Afectado por las Explotaciones de Carbón a cielo abierto y el Aprovechamiento Racional de estos Recursos Energéticos. (BOE nº 141 de 13 de junio de 1984, Mº de Industria y Energía). Los aspectos más destacables de este decreto, en lo que contempla las tareas de restauración, es la obligatoriedad de presentar un plan de restauración, donde hay una identificación del área de explotación y su entorno, descripción del medio físico y socioeconómico, además de un programa de restauración (acciones preventivas y correctivas). El Real Decreto 1116/1984 que hemos comentado, dispone que por el Ministro de Industria y Energía se dictarán las normas a las que deben ajustarse la elaboración y aprobación de los planes de explotación y restauración que deberán ser presentados por las explotaciones, dada la necesidad de adaptar los planes actuales, en su caso, y ajustar los nuevos a estos condicionantes. Dichas normas se recogen en la Orden de 13 de junio de 1984 sobre normas de elaboración de los planes de explotación y restauración del espacio natural afectado por las explotaciones de carbón a cielo abierto. (BOE nº 143 de 15 de junio de 1984, Mº de Industria y Energía). En la Orden se establecen las garantías o fianzas exigibles para el cumplimiento del plan de restauración que, pudiendo ser variables en el transcurso de la explotación, tengan en cuenta los daños pendientes de corrección en cada momento. El Real Decreto 2994/1982, de 15 de octubre, que estableció la obligación, a cargo de quienes realicen el aprovechamiento de recursos regulado por la Ley de Minas, de realizar trabajos de restauración del espacio natural afectado por las labores mineras, en los términos señalados en la misma disposición, permite optar porque los trabajos de restauración previstos en el proyecto o Plan de explotación sean realizados por el propio explotador, o que sea la Administración la encargada de ejecutar el Plan de restauración, obligándose al titular de la explotación minera a entregar a la Administración una cantidad periódica suficiente para cubrir el coste de ejecución del mismo. Por ello teniendo en cuenta que la disposición final del Real Decreto citado autoriza al Ministerio de Industria y Energía a adoptar las disposiciones necesarias para el desarrollo del mismo, en cuanto al sistema para llevar a cabo los ingresos y su posterior aplicación a financiar los Planes de restauración, así como en lo relativo a las garantías que pueden exigirse pará asegurar la ejecución de las labores de restauración cuando sea el titular de aprovechamiento o explotación quien deba llevar a efecto el Plan, aspecto que se desarrolla en la: [“APLICACIÓN DEL MÉTODO LANDSCAPE FUNCTION ANALYSIS (LFA) PARA EVALUAR EL ESTADO DE LAS RESTAURACIONES MINERAS. EL CASO DE LA MINERÍA DE CARBÓN A CIELO ABIERTO EN UTRILLAS (TERUEL)”.] 2013 83 - Orden de 20 de noviembre de 1984 desarrollo del RD 2994/1982 de 15 de octubre, sobre restauración del espacio natural afectado por actividades mineras. (BOE nº 285 de 28 de noviembre de 1984, Mº de Industria y Energía) Con posterioridad se ha legislado en relación con la seguridad que deben ofrecer las estructuras con el transcurso de la explotación minera. En concreto, con la promulgación del decreto que se menciona a continuación, se abordan aspectos relacionados con el tratamiento de las escombreras. Real Decreto 863/1985, de 2 de abril. Reglamento General de Normas Básicas sobre Seguridad Minera. Existen otras disposiciones en materia de Evaluación de Impacto Ambiental y protección del Dominio Público Hidráulico, que en el ámbito de las actuaciones en materia de rehabilitación de escombreras, es conveniente tener en cuenta: - Real Decreto Legislativo 1302/1986, de 28 de junio, de Evaluación del Impacto Ambiental. BOE nº 155, de 30 de junio de 1986. - Real Decreto 1131/1988, de 30 de septiembre, por el que se aprueba el Reglamento para la ejecución del Real Decreto Legislativo 1302/1986, de 28 de junio, de evaluación de impacto ambiental. - Ley 6/2001, de 8 de mayo, de modificación del Real Decreto Legislativo 1302/1986, de 28 de junio, de evaluación de impacto ambiental. BOE nº 111 de 9 de mayo de 2001. - Real Decreto 1/2008, de 11 de enero, por el que se aprueba el texto refundido de la Ley de Evaluación de Impacto Ambiental de proyectos. BOE nº 23 de 26/1/2008. - Real Decreto 9/2008, de 11 de enero, por el que se modifica el Reglamento del Dominio Público Hidráulico, aprobado por el Real Decreto 849/1986, de 11 de abril. BOE nº 14 miércoles 16 enero 2008 Y en cuanto a la disposición en materia de avales para la restauración minera: - Orden 5282/2002 de 25 de julio. Avales para responder de la restauración del espacio natural de explotaciones mineras y depósitos de lodo. Recientemente se ha publicado en el BOE el Real Decreto 975/2009 de 12 de junio sobre gestión de los residuos de las industrias extractivas y de protección y rehabilitación del espacio afectado por actividades mineras, que transpone la directiva 2006/21/ce del Parlamento Europeo y del Consejo de la UE de 15 de marzo de 2006 sobre la gestión de los residuos de industrias extractivas y por la que se modifica la Directiva 2004/35/CE que obliga a los Estados miembros a tomar todas las medidas necesarias para garantizar que los residuos de extracción se gestionen de un modo que no suponga peligro para la salud de las personas y sin utilizar procesos o métodos que puedan dañar el medio ambiente y, en particular, suponer riesgos para el agua, el aire, el suelo, la fauna o la flora, sin causar molestias debidas al ruido o los malos olores y sin afectar negativamente al paisaje ni a lugares que presenten un interés especial. [“APLICACIÓN DEL MÉTODO LANDSCAPE FUNCTION ANALYSIS (LFA) PARA EVALUAR EL ESTADO DE LAS RESTAURACIONES MINERAS. EL CASO DE LA MINERÍA DE CARBÓN A CIELO ABIERTO EN UTRILLAS (TERUEL)”.] 2013 84 Como resumen podemos establecer de forma general que existe una legislación abundante y minuciosa en el ámbito de la rehabilitación de espacios degradados por las actividades mineras y en particular por la minería del carbón. Sin embargo, existen numerosas zonas afectadas que se encuentran en una situación de vacío legal y en la actualidad nadie se considera responsable o competente para intervenir en su rehabilitación. [“APLICACIÓN DEL MÉTODO LANDSCAPE FUNCTION ANALYSIS (LFA) PARA EVALUAR EL ESTADO DE LAS RESTAURACIONES MINERAS. EL CASO DE LA MINERÍA DE CARBÓN A CIELO ABIERTO EN UTRILLAS (TERUEL)”.] 2013 85 8.3 Anexo 3 8.3.1 Metodología LFA 8.3.1.1 Ejemplo y ficha de trabajo para la caracterización de la organización del paisaje [“APLICACIÓN DEL MÉTODO LANDSCAPE FUNCTION ANALYSIS (LFA) PARA EVALUAR EL ESTADO DE LAS RESTAURACIONES MINERAS. EL CASO DE LA MINERÍA DE CARBÓN A CIELO ABIERTO EN UTRILLAS (TERUEL)”.] 2013 86 8.3.1.2 Ficha de trabajo para la evaluación de los indicadores del suelo: 8.3.1.3 Indicadores de las características del suelo  Protección frente al golpeo de gotas de lluvia: Evaluar el porcentaje de la cobertura de vegetación perenne proyectada a una altura de 0,5 m. Además de rocas> 2 cm y material leñoso> 1 cm de diámetro o más, así como objetos inamovibles de larga vida. Estos objetos interceptan y rompen las gotas de lluvia, por lo que se produce menos erosión y costra física del suelo.  Cobertura vegetal perenne: El objetivo es estimar la "cobertura basal" de la hierba perenne y / o la densidad del dosel de los árboles y arbustos.  Restos vegetales depositados: Se refiere a los pastos anuales y herbáceas efímeras (tanto de pie como individual), así como hojas desprendidas, tallos, ramas, frutas, estiércol etc. La posición de la basura en el paisaje global también ayuda en la definición de parches fértiles. La acumulación de hojarasca está fuertemente relacionada con los elementos carbono, nitrógeno y otros almacenados en la capas superficiales del suelo y adquiridos por procesos de descomposición. [“APLICACIÓN DEL MÉTODO LANDSCAPE FUNCTION ANALYSIS (LFA) PARA EVALUAR EL ESTADO DE LAS RESTAURACIONES MINERAS. EL CASO DE LA MINERÍA DE CARBÓN A CIELO ABIERTO EN UTRILLAS (TERUEL)”.] 2013 87 Hay tres propiedades de basura que deben ser evaluadas en el orden siguiente:  La cubierta (en 10 clases) de acuerdo con la tabla.  El origen: Local de (l) (ejemplo) = deriva de las plantas que crecen en una proximidad muy cercana a la zona de consulta y no muestra señales de transporte / depósito por el viento o el agua que fluye y transportado (t) (ejemplo) = tiene claros signos de haber sido desplazada.  El grado de descomposición / incorporación en 4 clases: 1. Descomposición Ninguna (n) (ejemplo): La basura se esparce libremente en la superficie con pocos signos de descomposición e incorporación. 2. Descomposición leve (s) (ejemplo): Basura se descompone en pequeños fragmentos y están íntimamente en contacto con el suelo; algunos fragmentos pueden estar parcialmente enterrados. 3. Descomposición moderada (m) (ejemplo): Aparición de varias capas distintas; algún ataque de los hongos es visible; la capa próxima al suelo está algo humificada; algunos oscurecimientos de la tierra a una profundidad de menos de 10 mm. 4. Descomposición extensa (e) (ejemplo): Aparecen al menos 3 capas o etapas descomposición que van desde material fresco en la parte superior hasta las capas más humificadas (muy oscuro, sin fragmentos identificables).  Cobertura de criptógamas: "Cryptogama" es un término genérico que incluye algas, hongos, líquenes, musgos y hepáticas. La fructificación de cuerpos de micorrizas se incluirían. Cuando éstos están presentes, indican que hay estabilidad en la superficie del suelo y niveles elevados de nutrientes disponibles en las capas superficiales del suelo.  Grado de fragmentación de la costra física: Los suelos con costras físicas en buen estado son costras que son suaves y se ajusta a las suaves ondulaciones en la superficie del suelo. Este tipo de cortezas pierden poco material del suelo en un evento de escorrentía. Sin embargo, las cortezas se vuelven inestables, frágiles y fácilmente perturbadas por los animales de pastoreo y los materiales que estén disponibles para la erosión hídrica y eólica.  Tipo y grado de erosión de suelo: Erosión en este contexto se refiere a la erosión acelerada causada por la interacción de la gestión y eventos climáticos. [“APLICACIÓN DEL MÉTODO LANDSCAPE FUNCTION ANALYSIS (LFA) PARA EVALUAR EL ESTADO DE LAS RESTAURACIONES MINERAS. EL CASO DE LA MINERÍA DE CARBÓN A CIELO ABIERTO EN UTRILLAS (TERUEL)”.] 2013 88 Cinco formas principales se describen a continuación:  Surcos y cárcavas: Son canales cortados por el agua que fluye. Arroyuelos con menos de 300 mm de profundidad mientras que los barrancos son mayores de 300 mm de profundidad (McDonald et al). Pueden ser iniciadas por el agua que fluye y pastoreo. Su presencia es un signo seguro de que el agua fluye rápidamente fuera del paisaje. A menudo lleva a la vez hojarasca y suelo con él. Ellos están alineados aproximadamente con la pendiente máxima.  Terrazas: Son abruptas paredes de 1 a 10cm de altura, alineados con el contorno local . El material erosionado se deposita en un abanico aluvial.  Erosión laminar: Es la eliminación progresiva de capas muy finas de suelo a través de extensas áreas. A veces se confunde con el “scalding”, pero típicamente se asocia con suelos de texturas gradacionales o uniforme. Muchas superficies laminadas están cubiertas por capas de grava o piedra que quedan después de la erosión de material más fino, cuando en una etapa avanzada.  Scalding: Es el resultado de una pérdida masiva del material del horizonte A u horizonte B que son típicamente muy duros, cuando se seca y tienen muy bajas tasas de infiltración. Los “scalding” tienen un potencial productivo cero, y capacidad para incrementar la escorrentía. A menudo se encuentran en terrenos planos, aunque no exclusivamente.  Pedestalling: Es el resultado de la eliminación de suelo por la erosión de un área a una profundidad de por lo menos varios cm, que salen de las culatas de las plantas supervivientes sobre una columna de suelo por encima del nuevo nivel general de la paisaje. Raíces expuestas son un sello distintivo de esta forma la erosión. Esta es una señal de que el propio tipo de suelo es muy erosionable y que la pérdida de la vegetación en el paisaje fue precedida por la erosión, y al revés. A menudo se asocia con ambientes mineros.  Deposición de materiales (sedimentos): El objetivo es evaluar la naturaleza y la cantidad de limo transportado y depositado. La presencia de suciedad y materiales en la zona de consulta indica la disponibilidad para el transporte de recursos procedentes de pendientes arriba en el paisaje y supone una cierta inestabilidad. Limos, arenas y gravas por lo general comprenden el aluvión. Los abanicos aluviales pueden llegar a ser bastante estables y productivos, en función de la tensión y la perturbación que se produce en la superficie. Un abanico aluvial puede convertirse en un parche productivo en poco tiempo si se producen las condiciones adecuadas de temporada. [“APLICACIÓN DEL MÉTODO LANDSCAPE FUNCTION ANALYSIS (LFA) PARA EVALUAR EL ESTADO DE LAS RESTAURACIONES MINERAS. EL CASO DE LA MINERÍA DE CARBÓN A CIELO ABIERTO EN UTRILLAS (TERUEL)”.] 2013 89  Rugosidad de la superficie: El objetivo es evaluar la rugosidad de la superficie por su capacidad para capturar y retener móvil recursos como el agua, propágulos, materia vegetal y orgánica.  Naturaleza del suelo (resistencia a la perturbación): El objetivo es evaluar la facilidad con la cual el suelo puede ser perturbado mecánicamente para producir material disponible para la erosión por el viento o el agua.  Test de humectación: La prueba se lleva a cabo mediante la inmersión en agua de fragmentos de suelo seco al aire de aproximadamente 1 cm y se observa la respuesta en un período de un minuto más o menos. Fragmentos del suelo que mantienen estables su cohesión cuando están mojados, implica un bajo potencial de erosión hídrica. - La calidad del agua es importante. - El agua salina es inadecuada. - La corteza del suelo debe permanecer más alta después de la inmersión  Textura 8.3.2 Fotografías de cada tipo de parche/interparche Calvero o suelo desnudo (BARE) [“APLICACIÓN DEL MÉTODO LANDSCAPE FUNCTION ANALYSIS (LFA) PARA EVALUAR EL ESTADO DE LAS RESTAURACIONES MINERAS. EL CASO DE LA MINERÍA DE CARBÓN A CIELO ABIERTO EN UTRILLAS (TERUEL)”.] 2013 96 Deposición de materiales (sedimentos) Rugosidad de la superficie del suelo Naturaleza del suelo (Resistencia a la perturbación) Test de humectación [“APLICACIÓN DEL MÉTODO LANDSCAPE FUNCTION ANALYSIS (LFA) PARA EVALUAR EL ESTADO DE LAS RESTAURACIONES MINERAS. EL CASO DE LA MINERÍA DE CARBÓN A CIELO ABIERTO EN UTRILLAS (TERUEL)”.] 2013 97 8.4.4 Coordenadas de los transectos Las coordenadas de los transectos están descritas en el datum ETRS89 en la zona 30 Ladera Transecto Longitud Latitud X Y L1 1 0o50´23.12´´ W 40o47´44.26`` N 682.246 4.518.315 2 0o50´22.83´´ W 40o47´44.14`` N 682.254 4.518.312 3 0o50´22.45´´ W 40o47´44.01`` N 682.262 4.518.308 4 0o50´22.11´´ W 40o47´44.84`` N 682.270 4.518.303 5 0o50´21.77´´ W 40o47´43.70`` N 682.278 4.518.299 6 0o50´21.47´´ W 40o47´43.53`` N 682.286 4.518.294 L2 1 0o50´22.81´´ W 40o47´46.41`` N 682.252 4.518.382 2 0o50´22.48´´ W 40o47´46.22`` N 682.260 4.518.376 3 0o50´22.12´´ W 40o47´46.01`` N 682.268 4.518.370 4 0o50´21.82´´ W 40o47´45.84`` N 682.276 4.518.365 5 0o50´21.46´´ W 40o47´45.65`` N 682.284 4.518.359 6 0o50´20.69´´ W 40o47´45.29`` N 682.302 4.518.349 L3 1 0o50´21.24´´ W 40o47´53.06`` N 682.284 4.518.588 2 0o50´20.92´´ W 40o47´53.06`` N 682.291 4.518.588 3 0o50´20.56´´ W 40o47´53.09`` N 682.300 4.518.589 4 0o50´20.21´´ W 40o47´53.15`` N 682.308 4.518.591 5 0o50´19.85´´ W 40o47´53.23`` N 682.316 4.518.594 6 0o50´19.51´´ W 40o47´53.32`` N 682.324 4.518.597 L4 1 0o50´25.66´´ W 40o47´47.18`` N 682.185 4.518.404 2 0o50´25.29´´ W 40o47´47.17`` N 682.193 4.518.404 3 0o50´24.97´´ W 40o47´47.16`` N 682.201 4.518.404 4 0o50´24.62´´ W 40o47´47.12`` N 682.209 4.518.403 5 0o50´24.27´´ W 40o47´47.05`` N 682.217 4.518.401 6 0o50´23.91´´ W 40o47´46.96`` N 682.226 4.518.398 L5 1 0o51´40.86´´ W 40o48´58.74`` N 680.368 4.520.568 2 0o51´40.56´´ W 40o48´59.09`` N 680.376 4.520.578 3 0o51´40.24´´ W 40o48´59.63`` N 680.382 4.520.595 4 0o51´39.9´´ W 40o48´59.93`` N 680.390 4.520.605 5 0o51´39.56´´ W 40o49´0.53`` N 680.398 4.520.624 6 0o51´39.2´´ W 40o40´0.97`` N 680.406 4.520.637 L6 1 0o51´48.48´´ W 40o48´57.23`` N 680.192 4.520.517 2 0o51´48.75´´ W 40o48´57.75`` N 680.184 4.520.533 3 0o51´49.01´´ W 40o48´58.12`` N 680.178 4.520.544 4 0o51´49.35´´ W 40o48´58.43`` N 680.170 4.520.553 5 0o51´49.67´´ W 40o48´58.82`` N 680.162 4.520.565 6 0o51´49.98´´ W 40o48´59.22`` N 680.154 4.520.577 L7 1 0o53´21.73´´ W 40o48´29.78`` N 678.027 4.519.617 2 0o53´22.06´´ W 40o48´29.79`` N 678.019 4.519.617 3 0o53´22.40´´ W 40o48´29.81`` N 678.011 4.519.618 4 0o53´22.75´´ W 40o48´29.83`` N 678.003 4.519.618 5 0o53´23.1´´ W 40o48´29.82`` N 677.995 4.519.618 6 0o53´23.43´´ W 40o48´29.81`` N 677.987 4.519.617 [“APLICACIÓN DEL MÉTODO LANDSCAPE FUNCTION ANALYSIS (LFA) PARA EVALUAR EL ESTADO DE LAS RESTAURACIONES MINERAS. EL CASO DE LA MINERÍA DE CARBÓN A CIELO ABIERTO EN UTRILLAS (TERUEL)”.] 2013 98 8.4 Anexo 4 8.4.1 Datos de los resultados del Índice de Infiltración, Estabilidad y Reciclaje de nutrientes para cada transecto y ladera:  Ladera 1 Nº de transecto II (%) IE (%) IRN (%) 1 18,4 52,1 16,5 2 19,8 45,2 13,9 3 18,5 38,9 11,8 4 17,5 50,2 15,1 5 19,6 49,1 14 6 18 45,5 14,6 Media 18,633 46,833 14,316  Ladera 2 Nº de transecto II (%) IE (%) IRN (%) 1 27,4 58,4 19,5 2 26,8 57,4 18,2 3 28,3 56,2 20 4 27,2 50,7 20,4 5 29,2 52,8 22,3 6 26,1 40,2 17,7 Media 27,5 52,616 19,683  Ladera 3 Nº de transecto II (%) IE (%) IRN (%) 1 40,3 61,1 29,7 2 37,6 59,3 25,1 3 39 60,9 27,3 4 40,7 62 29,7 5 38,3 58,8 26,4 6 33,6 56,9 21,7 Media 38,25 59,833 26,65 [“APLICACIÓN DEL MÉTODO LANDSCAPE FUNCTION ANALYSIS (LFA) PARA EVALUAR EL ESTADO DE LAS RESTAURACIONES MINERAS. EL CASO DE LA MINERÍA DE CARBÓN A CIELO ABIERTO EN UTRILLAS (TERUEL)”.] 2013 99  Ladera 4 Nº de transecto II (%) IE (%) IRN (%) 1 27,6 51,2 23,2 2 29,4 53,7 21,9 3 29,8 50,4 18,1 4 28,6 54,7 22,1 5 26,9 52,7 19,5 6 26,5 54,5 21,8 Media 28,133 52,866 21,1  Ladera 5 Nº de transecto II (%) IE (%) IRN (%) 1 25,6 39 12,9 2 24,6 37,4 11,3 3 23,2 36,5 11,1 4 21,9 40,7 12,5 5 21,6 40,5 11,9 6 16,5 42,5 11,4 Media 22,233 39,433 11,85  Ladera 6 Nº de transecto II (%) IE (%) IRN (%) 1 14,8 46,4 11,3 2 16,5 47,9 11,7 3 15,3 48,3 11 4 16,4 45,1 11,5 5 17,6 48,4 15,4 6 16,7 46,7 11,9 Media 16,216 47,133 12,133 [“APLICACIÓN DEL MÉTODO LANDSCAPE FUNCTION ANALYSIS (LFA) PARA EVALUAR EL ESTADO DE LAS RESTAURACIONES MINERAS. EL CASO DE LA MINERÍA DE CARBÓN A CIELO ABIERTO EN UTRILLAS (TERUEL)”.] 2013 100  Ladera 7 Nº de transecto II (%) IE (%) IRN (%) 1 33,9 59,2 28,7 2 29,2 58,9 25,4 3 30,2 61,6 28,3 4 32,5 61,5 28,3 5 31,8 58,2 26,1 6 30,7 58,6 26,7 Media 31,38333333 59,66666667 27,25  Laderas Ladera II (%) IE (%) IRN (%) L1 18,633 46,833 14,316 L2 27,5 52,616 19,683 L3 38,25 59,833 26,65 L4 28,133 52,866 21,1 L5 22,233 39,433 11,85 L6 16,216 47,133 12,133 L7 31,383 59,666 27,25 8.4.2 Datos de los resultados de la organización del paisaje para cada transecto y ladera:  Ladera 1: Transecto 1 2 3 4 5 6 Media Nº de parches / 10m transecto 3 1,3 2 1 2 0,7 1,666 Distancia media de parches (m) 2,89 5,9 3,66 5,52 3,27 7,3 4,756 Índice Org Paisaje 0,04 0,02 0,27 0,45 0,02 0,03 0,138 Área total parches (s.q m) 0,1 0,1 0,2 3,7 0,1 0,5 0,783 [“APLICACIÓN DEL MÉTODO LANDSCAPE FUNCTION ANALYSIS (LFA) PARA EVALUAR EL ESTADO DE LAS RESTAURACIONES MINERAS. EL CASO DE LA MINERÍA DE CARBÓN A CIELO ABIERTO EN UTRILLAS (TERUEL)”.] 2013 101  Ladera 2 Transecto 1 2 3 4 5 6 Media Nº de parches / 10m transecto 10,7 8 9 10,7 8,3 3,7 8,4 Distancia media de parches (m) 0,37 0,81 0,65 0,6 0,86 2,22 0,918 Índice Org Paisaje 0,59 0,41 0,42 0,32 0,32 0,11 0,361 Área total parches (s.q m) 25,9 27 22,8 6,9 7 0,9 15,083  Ladera 3 Transecto 1 2 3 4 5 6 Media Nº de parches / 10m transecto 6,7 6,1 6,7 8 4,7 7,7 6,65 Distancia media de parches (m) 0,62 1,07 0,92 0,8 1,2 1,07 0,946 Índice Org Paisaje 0,58 0,39 0,51 0,47 0,48 0,32 0,458 Área total parches (s.q m) 20,9 37 47,3 12,7 25,4 6 24,883  Ladera 4 Transecto 1 2 3 4 5 6 Media Nº de parches / 10m transecto 6,7 4 4,4 4,5 5,6 7,3 5,416 Distancia media de parches (m) 1,18 1,68 2,25 1,5 1,19 0,91 1,451 Índice Org Paisaje 0,18 0,33 0,09 0,16 0,25 0,31 0,22 Área total parches (s.q m) 3,1 6,5 1,3 4,2 9 16,2 6,716 [“APLICACIÓN DEL MÉTODO LANDSCAPE FUNCTION ANALYSIS (LFA) PARA EVALUAR EL ESTADO DE LAS RESTAURACIONES MINERAS. EL CASO DE LA MINERÍA DE CARBÓN A CIELO ABIERTO EN UTRILLAS (TERUEL)”.] 2013 102  Ladera 5 Transecto 1 2 3 4 5 6 Media Nº de parches / 10m transecto 3,7 4,3 3,7 3,7 4,9 4 4,05 Distancia media de parches (m) 2,3 2,22 2,36 2,34 2,68 2,2 2,35 Índice Org Paisaje 0,08 0,04 0,05 0,06 0,09 0,05 0,06 Área total parches (s.q m) 0,4 0,2 0,5 0,6 0,7 0,5 0,483  Ladera 6 Transecto 1 2 3 4 5 6 Media Nº de parches / 10m transecto 2,7 4,7 3,3 4 2,7 2,7 3,35 Distancia media de parches (m) 3,24 1,85 2,64 2,21 3,66 2,91 2,751 Índice Org Paisaje 0,02 0,07 0,03 0,04 0,03 0,02 0,035 Área total parches (s.q m) 0,1 0,4 0,1 0,2 0,1 0,2 0,1833  Ladera 7 Transecto 1 2 3 4 5 6 Media Nº de parches / 10m transecto 5,7 6 8,3 4 6,1 3,7 5,633 Distancia media de parches (m) 1,06 1,1 1,06 1,05 0,99 1,26 1,086 Índice Org Paisaje 0,29 0,19 0,33 0,48 0,37 0,37 0,338 Área total parches (s.q m) 4 2 3,1 18,6 15,2 30 12,15  Laderas Ladera 1 2 3 4 5 6 7 Nº de parches / 10m transecto 1,666 8,4 6,65 5,416 4,05 3,35 5,633 Distancia media de parches (m) 4,756 0,9183 0,946 1,451 2,35 2,751 1,086 Índice Org Paisaje 0,138 0,362 0,458 0,22 0,062 0,035 0,338 Área total parches (s.q m) 0,783 15,08 24,88 6,717 0,483 0,183 12,15 [“APLICACIÓN DEL MÉTODO LANDSCAPE FUNCTION ANALYSIS (LFA) PARA EVALUAR EL ESTADO DE LAS RESTAURACIONES MINERAS. EL CASO DE LA MINERÍA DE CARBÓN A CIELO ABIERTO EN UTRILLAS (TERUEL)”.] 2013 103 8.4.3 Datos y resultados del comportamiento individual de cada tipo de parche/interparche respecto al IE, II Y IRN: [“APLICACIÓN DEL MÉTODO LANDSCAPE FUNCTION ANALYSIS (LFA) PARA EVALUAR EL ESTADO DE LAS RESTAURACIONES MINERAS. EL CASO DE LA MINERÍA DE CARBÓN A CIELO ABIERTO EN UTRILLAS (TERUEL)”.] 2013 104 [“APLICACIÓN DEL MÉTODO LANDSCAPE FUNCTION ANALYSIS (LFA) PARA EVALUAR EL ESTADO DE LAS RESTAURACIONES MINERAS. EL CASO DE LA MINERÍA DE CARBÓN A CIELO ABIERTO EN UTRILLAS (TERUEL)”.] 2013 105 [“APLICACIÓN DEL MÉTODO LANDSCAPE FUNCTION ANALYSIS (LFA) PARA EVALUAR EL ESTADO DE LAS RESTAURACIONES MINERAS. EL CASO DE LA MINERÍA DE CARBÓN A CIELO ABIERTO EN UTRILLAS (TERUEL)”.] 2013 112 8.4.6 Datos y resultados del análisis Kruskal-Wallis Índices - Laderas  Índice Estabilidad-Laderas [“APLICACIÓN DEL MÉTODO LANDSCAPE FUNCTION ANALYSIS (LFA) PARA EVALUAR EL ESTADO DE LAS RESTAURACIONES MINERAS. EL CASO DE LA MINERÍA DE CARBÓN A CIELO ABIERTO EN UTRILLAS (TERUEL)”.] 2013 113  Índice Infiltración-Laderas Estadísticos de contrastea,b Indice estabilidad Chi-cuadrado 34,239 gl 6 Sig. asintót. ,000 a. Prueba de Kruskal-Wallis b. Variable de agrupación: Ladera Estadísticos de contrastea,b Indice estabilidad Chi-cuadrado 34,239 gl 6 Sig. asintót. ,000 a. Prueba de Kruskal-Wallis b. Variable de agrupación: Ladera [“APLICACIÓN DEL MÉTODO LANDSCAPE FUNCTION ANALYSIS (LFA) PARA EVALUAR EL ESTADO DE LAS RESTAURACIONES MINERAS. EL CASO DE LA MINERÍA DE CARBÓN A CIELO ABIERTO EN UTRILLAS (TERUEL)”.] 2013 114  Índice de Reciclaje de Nutrientes-Laderas [“APLICACIÓN DEL MÉTODO LANDSCAPE FUNCTION ANALYSIS (LFA) PARA EVALUAR EL ESTADO DE LAS RESTAURACIONES MINERAS. EL CASO DE LA MINERÍA DE CARBÓN A CIELO ABIERTO EN UTRILLAS (TERUEL)”.] 2013 115 Organización del Paisaje-Laderas  Área total de parches-Laderas [“APLICACIÓN DEL MÉTODO LANDSCAPE FUNCTION ANALYSIS (LFA) PARA EVALUAR EL ESTADO DE LAS RESTAURACIONES MINERAS. EL CASO DE LA MINERÍA DE CARBÓN A CIELO ABIERTO EN UTRILLAS (TERUEL)”.] 2013 116  Distancia media entre parches-Laderas [“APLICACIÓN DEL MÉTODO LANDSCAPE FUNCTION ANALYSIS (LFA) PARA EVALUAR EL ESTADO DE LAS RESTAURACIONES MINERAS. EL CASO DE LA MINERÍA DE CARBÓN A CIELO ABIERTO EN UTRILLAS (TERUEL)”.] 2013 117  Índice Organización del Paisaje-Laderas [“APLICACIÓN DEL MÉTODO LANDSCAPE FUNCTION ANALYSIS (LFA) PARA EVALUAR EL ESTADO DE LAS RESTAURACIONES MINERAS. EL CASO DE LA MINERÍA DE CARBÓN A CIELO ABIERTO EN UTRILLAS (TERUEL)”.] 2013 118  Nº de parches/10m de transecto [“APLICACIÓN DEL MÉTODO LANDSCAPE FUNCTION ANALYSIS (LFA) PARA EVALUAR EL ESTADO DE LAS RESTAURACIONES MINERAS. EL CASO DE LA MINERÍA DE CARBÓN A CIELO ABIERTO EN UTRILLAS (TERUEL)”.] 2013 119 Tipo de parche /interparche – Índices  Índice de EstabilidadTipo de parche [“APLICACIÓN DEL MÉTODO LANDSCAPE FUNCTION ANALYSIS (LFA) PARA EVALUAR EL ESTADO DE LAS RESTAURACIONES MINERAS. EL CASO DE LA MINERÍA DE CARBÓN A CIELO ABIERTO EN UTRILLAS (TERUEL)”.] 2013 120  Índice de infiltración-Tipo de parche [“APLICACIÓN DEL MÉTODO LANDSCAPE FUNCTION ANALYSIS (LFA) PARA EVALUAR EL ESTADO DE LAS RESTAURACIONES MINERAS. EL CASO DE LA MINERÍA DE CARBÓN A CIELO ABIERTO EN UTRILLAS (TERUEL)”.] 2013 121  Índice de Reciclaje de NutrientesTipo de parche [“APLICACIÓN DEL MÉTODO LANDSCAPE FUNCTION ANALYSIS (LFA) PARA EVALUAR EL ESTADO DE LAS RESTAURACIONES MINERAS. EL CASO DE LA MINERÍA DE CARBÓN A CIELO ABIERTO EN UTRILLAS (TERUEL)”.] 2013 128