Trabajos de restauración de Zostera noltii en la costa vasca: actividades de 2010 Informe para: Agencia Vasca del Agua – Uraren Euskal Agentzia Albert Einstein E7 46-3º (Parque Tecnológico de Álava), 01510 Miñano (Álava) Pasaia, 21 de diciembre de 2010
106
Embed
Trabajos de restauración de Zostera noltii en la costa vasca ......fanerógamas marinas en la costa vasca, concretamente de Zostera noltii. En el primer año de estudio, llevado a
This document is posted to help you gain knowledge. Please leave a comment to let me know what you think about it! Share it to your friends and learn new things together.
Transcript
Trabajos de restauración de Zostera noltii en la costa vasca:
actividades de 2010
Informe para:
Agencia Vasca del Agua – Uraren Euskal Agentzia
Albert Einstein E7 46-3º (Parque Tecnológico de Álava), 01510 Miñano (Álava)
Tipo documentoTipo documentoTipo documentoTipo documento Informe de actividades de 2010
Titulo documentoTitulo documentoTitulo documentoTitulo documento Trabajos de restauración de Zostera noltii en la costa vasca: actividades de 2010
FechaFechaFechaFecha 21/12/2010
ProyectoProyectoProyectoProyecto Red de Calidad de Aguas
CódigoCódigoCódigoCódigo ATM2010RedCalidad
ClienteClienteClienteCliente Agencia Vasca del Agua
Equipo de proyecto:Equipo de proyecto:Equipo de proyecto:Equipo de proyecto: Dr. Joxe Mikel Garmendia
Mireia Valle
Dr. Guillem Chust
Dr. Javier Franco
Dr. Ángel Borja
Responsable Responsable Responsable Responsable proyectoproyectoproyectoproyecto Dr. Joxe Mikel Garmendia
Revisado porRevisado porRevisado porRevisado por Dr. Javier Franco
FechaFechaFechaFecha 21-12-2010
Si procede, este documento deberá ser citado del siguiente modo:
Garmendia, J.M., M. Valle, G. Chust, J. Franco y Á. Borja, 2010. Trabajos de restauración de Zostera noltii en la costa vasca: actividades de 2010. Informe elaborado por AZTI-Tecnalia para la Agencia Vasca del Agua – Uraren Euskal Agentzia, Departamento de Medio Ambiente, Planificación Territorial, Agricultura y Pesca del Gobierno Vasco, 106 pp.
4.2.1.4.2.1.4.2.1.4.2.1. Seguimiento del trasplante en sedimento arenosoSeguimiento del trasplante en sedimento arenosoSeguimiento del trasplante en sedimento arenosoSeguimiento del trasplante en sedimento arenoso............................................... 15
4.2.2.4.2.2.4.2.2.4.2.2. Seguimiento del trasplante en sedimento fangosoSeguimiento del trasplante en sedimento fangosoSeguimiento del trasplante en sedimento fangosoSeguimiento del trasplante en sedimento fangoso............................................... 24
4.2.3.4.2.3.4.2.3.4.2.3. Comparación de trasplantes en arenas y fangosComparación de trasplantes en arenas y fangosComparación de trasplantes en arenas y fangosComparación de trasplantes en arenas y fangos.................................................. 31
4.3.2.4.3.2.4.3.2.4.3.2. Material y métodosMaterial y métodosMaterial y métodosMaterial y métodos................................................................................................. 40
4.4.1.4.4.1.4.4.1.4.4.1. Labores de trasplanteLabores de trasplanteLabores de trasplanteLabores de trasplante............................................................................................. 55
4.4.2.4.4.2.4.4.2.4.4.2. Caracterización de las zonas donante y receptoraCaracterización de las zonas donante y receptoraCaracterización de las zonas donante y receptoraCaracterización de las zonas donante y receptora ............................................... 62
4.4.3.4.4.3.4.4.3.4.4.3. Seguimiento de lSeguimiento de lSeguimiento de lSeguimiento de los trasplantesos trasplantesos trasplantesos trasplantes.............................................................................. 64
4.5 RECUPERACIÓN DE LA ZONA DONANTE ................................................72
4.6 TALLER DE TRABAJO SOBRE PRADERAS MARINAS..............................77
4.6.1.4.6.1.4.6.1.4.6.1. I European Seagrass Restoration WorksI European Seagrass Restoration WorksI European Seagrass Restoration WorksI European Seagrass Restoration Workshophophophop........................................................ 77
4.6.2.4.6.2.4.6.2.4.6.2. Conclusiones del tallerConclusiones del tallerConclusiones del tallerConclusiones del taller ........................................................................................... 78
Tal y como se comentó en el informe del año 2009, a la hora de la planificación de
actividades de restauración, resulta interesante y necesario disponer del conocimiento de los
niveles de variabilidad genética, diferenciación y flujo genético de las poblaciones de Zostera
noltii existentes en la actualidad en el País Vasco.
El análisis realizado el pasado año fue efectuado a partir de 7, 54 y 70 muestras de los
estuarios del Lea, Bidasoa y Oka respectivamente. Esta diferente contribución de muestras
por parte de cada estuario pudo influir en el resultado, por lo que se decidió aumentar el
número de muestras procedentes del Lea y realizar de nuevo el análisis genético. Por ello,
con la intención de equilibrar el número de muestras procedentes de cada estuario, se
añadieron otras 40 muestras a las del estuario del Lea (recogidas el 18 de octubre de 2009) y
se volvió a proceder de la misma manera.
De esta manera, se recogieron 171 muestras de los tres estuarios del País Vasco
(Figura 1) que cuentan con poblaciones de Zostera noltii: Lea, Bidasoa y Oka (Garmendia et
al., 2010). En cada punto seleccionado (47 en Lea, 54 en Bidasoa y 70 en Oka) se tomó un haz
o pie (conjunto de 2-4 hojas) de la planta y, bajo las necesarias condiciones de conservación,
todas las muestras fueron enviadas a la Universidad de Algarve (Portugal) para su análisis
genético.
SAN SEBASTI
BILBAO
MAR CANT ÁBRICO
0 10 20
Escacale(km)
N
BARBADUN
NERBIOI
BUTROE
OKA LEA
ARTIBAI
DEBAUROLA ORIA
URUMEAOIARTZUN
BIDASOA
2º W3º W
43º 20’N
SAN SEBASTIÁN
BILBAO
MAR CANTÁBRICO
0 10 200 10 20 km
N
BARBADUN
NERBIOI
BUTROE
OKA LEA
ARTIBAI
DEBAUROLA ORIA
URUMEAOIARTZUN
BIDASOA
2º W3º W
43º 20’N
�
SAN SEBASTI
BILBAO
MAR CANT ÁBRICO
0 10 20
Escacale(km)
0 10 20
Escacale(km)
NNN
BARBADUN
NERBIOI
BUTROE
OKA LEA
ARTIBAI
DEBAUROLA ORIA
URUMEAOIARTZUN
BIDASOA
2º W3º W
43º 20’N
SAN SEBASTIÁN
BILBAO
MAR CANTÁBRICO
0 10 200 10 20 km
N
BARBADUN
NERBIOI
BUTROE
OKA LEA
ARTIBAI
DEBAUROLA ORIA
URUMEAOIARTZUN
BIDASOA
2º W3º W
43º 20’N
�
Figura Figura Figura Figura 1111.... Estuarios de la CAPV. Se resaltan con un cuadro rojo aquellos que presentan poblaciones de Zostera noltii en la actualidad (no se consideran la poblaciones de Butroe producto de los trabajos de trasplante).
Las hojas fueron recogidas al azar con una separación mínima de 1 m entre cada
punto. Las muestras fueron caracterizadas mediante el contenido genético (genotipo)
utilizando 7 loci microsatélites (Coyer et al., 2004). Tras la eliminación de genotipos
multilocus duplicados1 (ejemplares clónicos o ramet pertenecientes al mismo individuo
genético o genet), se calculó la diferenciación genética (Fst); la heterocigosidad (HE y HO); y el
coeficiente de endogamia (FIS). Todo ello se realizó mediante el uso del software Genetix.
También se calculó la variación genética como la riqueza alélica (Â) estandarizada por el
tamaño de muestra. Dicha estandarización se llevó a cabo de acuerdo con el número de
genets y también por el número de ramets (Tabla 1). Por otro lado, también se recurrió al
Análisis Factorial de Correspondencias (AFC), disponible en la versión 4.04 de Genetix
(Belkhir et al., 2004) con el fin de visualizar individuos en un espacio multidimensional
basado en datos alélicos. Este procedimiento no agrupa a priori individuos en grupos, y por
tanto proporciona una ilustración más objetiva de la estructura de la población.
Se llevó a cabo un “test de asignación” utilizando Genclass 2.0 (Piry et al., 2004)
siguiendo el método bayesiano y el cálculo de probabilidad (Rannala y Mountain, 1997). Este
test nos informa sobre el origen de la población.
Tabla Tabla Tabla Tabla 1111.... Número de individuos o muestras recolectadas (N), Número de genets (G), Diversidad genotípica (R) y Diversidad genética (Â), Coeficiente de endogamia (FIS), Heterocigosidad esperada (HE) y observada (HO) para las tres poblaciones de Zostera noltii de la CAPV. Nc: No calculado. (*) Desviación significativa (p<0,05).
PoblaciónPoblaciónPoblaciónPoblación NNNN GGGG RRRR Â Â Â Â
La diversidad genotípica varía entre 0,23 en Lea y 0,73 en Oka (Tabla 1). La variación
en Bidasoa es intermedia (0,55). La heterocigosidad observada fue de 0,31 en Bidasoa y 0,59
y 0,49 en Lea y Oka respectivamente. En el caso de Lea, se observó una desviación
significativa frente al equilibrio de Hardy Weinberg como resultado de la excesiva
heterocigosidad. La riqueza alélica fue de 2,78 en Bidasoa, siendo mucho más elevada en Oka
1 El sistema de crecimiento clónico de las fanerógamas marinas hace difícil la caracterización individual de las plantas. Por ello, los métodos genéticos son necesarios para diferenciar los genets (individuos genéticos) de los ramets (individuos morfológicos generados a partir de la reproducción vegetativa). Una vez realizado esto, se puede llevar a cabo un adecuado análisis de la diversidad genotípica y flujo de genes.
y Lea tras la estandarización por el tamaño de muestra de Lea. En la comparación entre Oka
y Bidasoa (estandarizado por el tamaño de muestra del Bidasoa), la población de Oka
presenta una riqueza alélica dos veces superior a la del Bidasoa.
La representación del análisis de correspondencias factorial muestra que dos ejes
pueden explicar toda la variación de los datos (Eje 1 = 72% y Eje 2 = 28%). Además, la figura
muestra que las poblaciones están claramente diferenciadas, formando distintos grupos. Este
hecho también se aprecia en los valores de Fst o diferenciación de la población (Figura 2).
Altos valores de Fst (dentro del rango 0-1) se obtienen entre las poblaciones más
distanciadas. Oka y Lea se encuentran geográficamente más próximas y presentan un valor
de Fst menor (0,14) que los de Oka y Lea con Bidasoa (0,32 y 0,40 respectivamente).
Figura Figura Figura Figura 2222.... Representación de los dos primeros ejes (factores) del análisis de correspondencias factorial (AFC) basado en la variación alélica en 7 loci de microsatélite para 94 genets de Zostera noltii de los estuarios del Lea (en blanco), Oka (en azul) y Bidasoa (en amarillo). Fst=diferenciación de la población.
Al analizar cada población a la escala de clones, se observa que los clones son de gran
tamaño. En el caso de Bidasoa, algunos individuos recolectados en diferentes parches o zonas
muy distanciadas (A-G) pertenecen al mismo genotipo multilocus. También se ha podido
CRECIMIENTO DE LOS PARCHES DE LA ESTACICRECIMIENTO DE LOS PARCHES DE LA ESTACICRECIMIENTO DE LOS PARCHES DE LA ESTACICRECIMIENTO DE LOS PARCHES DE LA ESTACIÓÓÓÓN ARENOSA B1N ARENOSA B1N ARENOSA B1N ARENOSA B1
CRECIMIENTO DE LOS PARCHES DE LA ESTACICRECIMIENTO DE LOS PARCHES DE LA ESTACICRECIMIENTO DE LOS PARCHES DE LA ESTACICRECIMIENTO DE LOS PARCHES DE LA ESTACIÓÓÓÓN ARENOSA B1N ARENOSA B1N ARENOSA B1N ARENOSA B1
Figura Figura Figura Figura 3333.... Desarrollo del parche de Zostera noltii B1-1 perteneciente a la zona arenosa. Los números del interior de los recuadros corresponden al factor de crecimiento (veces que ha crecido el parche respecto a su superficie inicial). A partir de diciembre se observan repetidos fenómenos de sedimentación (reflejado en un enterramiento parcial de los parches) y en junio se produce una importante riada que provoca la erosión del sedimento y la pérdida de gran cantidad de arena de esta zona.
CRECIMIENTO DE LOS PARCHES DE LA ESTACICRECIMIENTO DE LOS PARCHES DE LA ESTACICRECIMIENTO DE LOS PARCHES DE LA ESTACICRECIMIENTO DE LOS PARCHES DE LA ESTACIÓÓÓÓN ARENOSA B1N ARENOSA B1N ARENOSA B1N ARENOSA B1
CRECIMIENTO DE LOS PARCHES DE LA ESTACICRECIMIENTO DE LOS PARCHES DE LA ESTACICRECIMIENTO DE LOS PARCHES DE LA ESTACICRECIMIENTO DE LOS PARCHES DE LA ESTACIÓÓÓÓN ARENOSA B1N ARENOSA B1N ARENOSA B1N ARENOSA B1
Figura Figura Figura Figura 4444.... Desarrollo del parche de Zostera noltii B1-2 perteneciente a la zona arenosa. Los números del interior de los recuadros corresponden al factor de crecimiento (veces que ha crecido el parche respecto a su superficie inicial). A partir de diciembre se observan repetidos fenómenos de sedimentación (reflejado en un enterramiento parcial de los parches) y en junio se produce una importante riada que provoca la erosión del sedimento y la pérdida de gran cantidad de arena de esta zona.
CRECIMIENTO DE LOS PARCHES DE LA ESTACICRECIMIENTO DE LOS PARCHES DE LA ESTACICRECIMIENTO DE LOS PARCHES DE LA ESTACICRECIMIENTO DE LOS PARCHES DE LA ESTACIÓÓÓÓN ARENOSA B1N ARENOSA B1N ARENOSA B1N ARENOSA B1
CRECIMIENTO DE LOS PARCHES DE LA ESTACICRECIMIENTO DE LOS PARCHES DE LA ESTACICRECIMIENTO DE LOS PARCHES DE LA ESTACICRECIMIENTO DE LOS PARCHES DE LA ESTACIÓÓÓÓN ARENOSA B1N ARENOSA B1N ARENOSA B1N ARENOSA B1
Figura Figura Figura Figura 5555.... Desarrollo del parche de Zostera noltii B1-3 perteneciente a la zona arenosa. Los números del interior de los recuadros corresponden al factor de crecimiento (veces que ha crecido el parche respecto a su superficie inicial). A partir de diciembre se observan repetidos fenómenos de sedimentación (reflejado en un enterramiento parcial de los parches) y en junio se produce una importante riada que provoca la erosión del sedimento y la pérdida de gran cantidad de arena de esta zona.
CRECIMIENTO DE LOS PARCHES DE LA ESTACICRECIMIENTO DE LOS PARCHES DE LA ESTACICRECIMIENTO DE LOS PARCHES DE LA ESTACICRECIMIENTO DE LOS PARCHES DE LA ESTACIÓÓÓÓN ARENOSA B1N ARENOSA B1N ARENOSA B1N ARENOSA B1
CRECIMIENTO DE LOS PARCHES DE LA ESTACICRECIMIENTO DE LOS PARCHES DE LA ESTACICRECIMIENTO DE LOS PARCHES DE LA ESTACICRECIMIENTO DE LOS PARCHES DE LA ESTACIÓÓÓÓN ARENOSA B1N ARENOSA B1N ARENOSA B1N ARENOSA B1
FigFigFigFigura ura ura ura 6666.... Desarrollo del parche de Zostera noltii B1-4 perteneciente a la zona arenosa. Los números del interior de los recuadros corresponden al factor de crecimiento (veces que ha crecido el parche respecto a su superficie inicial). A partir de diciembre se observan repetidos fenómenos de sedimentación (reflejado en un enterramiento parcial de los parches) y en junio se produce una importante riada que provoca la erosión del sedimento y la pérdida de gran cantidad de arena de esta zona.
Figura 7.Figura 7.Figura 7.Figura 7. Aspecto general de la zona de trasplante antes (mayo) y después (junio) de la riada de mediados de junio 2010 que alteró la zona de los sedimentos arenosos (señalada con una flecha).
Figura 8.Figura 8.Figura 8.Figura 8. Aspecto general de la zona de trasplante arenosa (B1) antes (febrero, abril y mayo) de la riada ocurrida en junio de 2010. Puede observarse que la zona es prácticamente llana.
Figura 9.Figura 9.Figura 9.Figura 9. Aspecto general de la zona de trasplante arenosa (B1) después (finales de junio) de la riada ocurrida a mediados de junio de 2010. Puede observarse que la zona ahora presenta gran cantidad de montículos y pozas, lo cual es consecuencia del movimiento de grandes cantidades de sedimento, que ha supuesto un gran cambio en la morfología de la zona.
CRECIMIENTO DE LOS PARCHES DE LA ESTACICRECIMIENTO DE LOS PARCHES DE LA ESTACICRECIMIENTO DE LOS PARCHES DE LA ESTACICRECIMIENTO DE LOS PARCHES DE LA ESTACIÓÓÓÓN FANGOSA B2N FANGOSA B2N FANGOSA B2N FANGOSA B2
CRECIMIENTO DE LOS PARCHES DE LA ESTACICRECIMIENTO DE LOS PARCHES DE LA ESTACICRECIMIENTO DE LOS PARCHES DE LA ESTACICRECIMIENTO DE LOS PARCHES DE LA ESTACIÓÓÓÓN FANGOSA B2N FANGOSA B2N FANGOSA B2N FANGOSA B2
Figura 10.Figura 10.Figura 10.Figura 10. Desarrollo del parche de Zostera noltii B2-1 perteneciente a la zona fangosa. Los números del interior de los recuadros corresponden al factor de crecimiento (veces que ha crecido el parche respecto a su superficie inicial).
CRECIMIENTO DE LOS PARCHES DE LA ESTACICRECIMIENTO DE LOS PARCHES DE LA ESTACICRECIMIENTO DE LOS PARCHES DE LA ESTACICRECIMIENTO DE LOS PARCHES DE LA ESTACIÓÓÓÓN FANGOSA B2N FANGOSA B2N FANGOSA B2N FANGOSA B2
CRECIMIENTO DE LOS PARCHES DE LA ESTACICRECIMIENTO DE LOS PARCHES DE LA ESTACICRECIMIENTO DE LOS PARCHES DE LA ESTACICRECIMIENTO DE LOS PARCHES DE LA ESTACIÓÓÓÓN FANGOSA B2N FANGOSA B2N FANGOSA B2N FANGOSA B2
Figura 11.Figura 11.Figura 11.Figura 11. Desarrollo del parche de Zostera noltii B2-2 perteneciente a la zona fangosa. Los números del interior de los recuadros corresponden al factor de crecimiento (veces que ha crecido el parche respecto a su superficie inicial).
CRECIMIENTO DE LOS PARCHES DE LA ESTACICRECIMIENTO DE LOS PARCHES DE LA ESTACICRECIMIENTO DE LOS PARCHES DE LA ESTACICRECIMIENTO DE LOS PARCHES DE LA ESTACIÓÓÓÓN FANGOSA B2N FANGOSA B2N FANGOSA B2N FANGOSA B2
CRECIMIENTO DE LOS PARCHES DE LA ESTACICRECIMIENTO DE LOS PARCHES DE LA ESTACICRECIMIENTO DE LOS PARCHES DE LA ESTACICRECIMIENTO DE LOS PARCHES DE LA ESTACIÓÓÓÓN FANGOSA B2N FANGOSA B2N FANGOSA B2N FANGOSA B2
Figura 12.Figura 12.Figura 12.Figura 12. Desarrollo del parche de Zostera noltii B2-3 perteneciente a la zona fangosa. Los números del interior de los recuadros corresponden al factor de crecimiento (veces que ha crecido el parche respecto a su superficie inicial).
CRECIMIENTO DE LOS PARCHES DE LA ESTACICRECIMIENTO DE LOS PARCHES DE LA ESTACICRECIMIENTO DE LOS PARCHES DE LA ESTACICRECIMIENTO DE LOS PARCHES DE LA ESTACIÓÓÓÓN FANGOSA B2N FANGOSA B2N FANGOSA B2N FANGOSA B2
CRECIMIENTO DE LOS PARCHES DE LA ESTACICRECIMIENTO DE LOS PARCHES DE LA ESTACICRECIMIENTO DE LOS PARCHES DE LA ESTACICRECIMIENTO DE LOS PARCHES DE LA ESTACIÓÓÓÓN FANGOSA B2N FANGOSA B2N FANGOSA B2N FANGOSA B2
Figura 13.Figura 13.Figura 13.Figura 13. Desarrollo del parche de Zostera noltii B2-4 perteneciente a la zona fangosa. Los números del interior de los recuadros corresponden al factor de crecimiento (veces que ha crecido el parche respecto a su superficie inicial).
Figura Figura Figura Figura 14141414.... Ulva sp. en la zona central intermareal cercana a los puntos de trasplante de los parches de Zostera noltii en sedimento fangoso del estuario del Butroe. Abril de 2010.
Figura Figura Figura Figura 15151515.... El alga Ulva sp. cubriendo casi completamente el parche trasplantado de Zostera noltii B2-3 (correspondiente al sedimento fangoso) en abril (a y b) y mayo (c y d) de 2010.
4.2.3. Comparación de trasplantes en arenas y fango s
A continuación, para facilitar la comparación del desarrollo de los parches en ambos
tipos de sedimento, se presentan algunas figuras representando conjuntamente los datos
para ambas zonas de trasplante (arenosa y fangosa).
Crecimiento de los parches de Crecimiento de los parches de Crecimiento de los parches de Crecimiento de los parches de Zostera noltiiZostera noltiiZostera noltiiZostera noltii
0,0
1,0
2,0
3,0
4,0
5,0
6,0
7,0
Mr
09
Ab My Jn Jl Ag Sp Oc No Di En
10
Fe Mr Ab My Jn Ag Oc
mesmesmesmes
Fa
ctor
de
crec
imie
nto
Fa
ctor
de
crec
imie
nto
Fa
ctor
de
crec
imie
nto
Fa
ctor
de
crec
imie
nto
B1-1
B1-2
B1-3
B1-4
B2-1
B2-2
B2-3
B2-4
sedimentación
riada / erosión
Figura Figura Figura Figura 16161616.... Incremento (factor de crecimiento) de la superficie de cada uno de los parches trasplantados de Zostera noltii. B1 en sedimento arenoso y B2 en sedimento fangoso. A partir de diciembre se observan repetidos fenómenos de sedimentación (reflejado en un enterramiento parcial de los parches) especialmente en la zona arenosa, y en junio se produce una importante riada que provoca la erosión del sedimento y la pérdida de gran cantidad de arena de la zona B1.
En la figura 16 y tabla 2 se observa claramente el mejor crecimiento inicial
experimentado por los parches de la zona arenosa, y su posterior disminución hasta su
desaparición. Mientras tanto, los parches fangosos, aunque no ofrezcan un crecimiento
rápido, sobreviven durante un tiempo más prolongado, aunque al final parece intuirse cierta
Tabla Tabla Tabla Tabla 2222.... Factor de crecimiento (media y desviación típica, d.t.) de la superficie de los parches trasplantados de Zostera noltii. B1 en zona arenosa y B2 en zona fangosa.
Figura Figura Figura Figura 17171717.... Evolución de la densidad promedio (barras: desviación típica) a lo largo del seguimiento post-trasplante en las zonas arenosa (B1) y fangosa (B2) del estuario del Butroe. A partir de diciembre se observan repetidos fenómenos de sedimentación (reflejado en un enterramiento parcial de los parches) especialmente en la zona arenosa, y en junio se produce una importante riada que provoca la erosión del sedimento y la pérdida de gran cantidad de arena de la zona B1.
4.2.4. Discusión
Teniendo en cuenta los resultados observados a lo largo de los primeros meses de
seguimiento, el trasplante se presumía exitoso (Garmendia et al., 2010). Sin embargo, en la
segunda parte de este seguimiento (la correspondiente al año 2010) se han podido identificar
dos presiones que han resultado cruciales para la disminución y desaparición de alguno de
los parches trasplantados. En primer lugar, a partir de diciembre de 2009 se observaron unos
eventos de sedimentación reflejados en un mayor o menor grado de enterramiento de los
parches. La aparición y desaparición (o mejor dicho, desenterramiento y enterramiento) de
las estacas de referencia fijadas en los bordes de los parches en la zona arenosa reflejaron el
elevado dinamismo que sufre el sedimento en esta zona (hay que recordar que se encuentra
muy cerca del canal principal del estuario). Como consecuencia de este enterramiento, la
visualización y medida de los parches visitados a partir de este invierno resultó muy difícil y
laboriosa. Por tanto, aunque no se descarta el efecto negativo directo producido por el
enterramiento en la supervivencia de las plantas, también se considera que las superficies
calculadas subestiman la extensión real que ocupa cada parche debido a esta dificultad a la
podrán tenerse en cuenta los efectos causados por este proceso y adecuar la selección de las
zonas a recuperar a este nuevo escenario ambiental.
El contenido del presente capítulo ha sido publicado en la Revista de Investigación
Marina (Valle et al., 2010).
4.3.2. Material y métodos
Área de estudioÁrea de estudioÁrea de estudioÁrea de estudio
El área de estudio correspondiente al presente capítulo se centra en el estuario del
Oka, ubicado en la zona centro-occidental de la costa del País Vasco. Sus características
hidromorfológicas más relevantes pueden verse en Borja et al. (2006). Este estuario es uno de
los biológicamente más diversos y mejor conservados del País Vasco, junto con el del Bidasoa
(Borja et al., 2004). Debido a ello, se encuentra amparado bajo diferentes figuras de
protección y conservación de ámbito autonómico e internacional (Castro et al., 2004). En sus
marismas se desarrolla un número considerable de especies incluidas en el Catálogo Vasco de
especies amenazadas (Otxoa et al., 2007), entre ellas Z. noltii, especie objeto del presente
estudio.
Recopilación de datosRecopilación de datosRecopilación de datosRecopilación de datos
En este apartado se describen los diferentes tipos de datos recopilados para la
elaboración del trabajo y sus fuentes.
a) Localización geográfica de las poblaciones de a) Localización geográfica de las poblaciones de a) Localización geográfica de las poblaciones de a) Localización geográfica de las poblaciones de Z. noltiiZ. noltiiZ. noltiiZ. noltii
La distribución geográfica de la especie se obtuvo a partir del Informe ‘Trabajos de
restauración de Zostera noltii en la Costa Vasca’ (Garmendia et al., 2008) y partir del trabajo
de cartografiado a escala 1:5.000 de las poblaciones vizcaínas de plantas incluidas en el
Catálogo Vasco de Especies Amenazadas (Otxoa et al., 2007).
En junio de 2009 se realizó un muestreo complementario para comprobar el correcto
posicionamiento de las manchas y georreferenciar manchas nuevas. Durante este muestreo
se cartografió un nuevo polígono de Z. noltii. Las coordenadas X e Y (en Proyección UTM y
Datum WGS84) del contorno de los polígonos se recogieron mediante un GPS modelo
Integración de datosIntegración de datosIntegración de datosIntegración de datos
Para el correcto funcionamiento del ENFA es preciso enmascarar cada uno de los
mapas ecogeográficos obtenidos a partir de cada una de las variables, ya que el área cubierta
por cada mapa debe ser exactamente la misma para poder realizar el análisis. Por lo tanto, se
creó una máscara delimitando la zona intermareal del estuario, zona potencial de
distribución de la especie (Garmendia et al., 2008). El proceso de enmascarado se realizó
mediante la opción ‘Raster calculator’ de la extensión ‘Spatial Analyst’, implementada en el
programa ArcGis 9.2.
Análisis Factorial de Nicho Ecológico y Mapa de Hábitat IdóneoAnálisis Factorial de Nicho Ecológico y Mapa de Hábitat IdóneoAnálisis Factorial de Nicho Ecológico y Mapa de Hábitat IdóneoAnálisis Factorial de Nicho Ecológico y Mapa de Hábitat Idóneo
Para elaborar el modelo de HI y obtener el MHI, se utilizó el software BioMapper
versión 4 (http://www2.unil.ch/biomapper/), dentro del cual está implementada la
herramienta ENFA, desarrollada por Hirzel et al. (2002). El ENFA calcula funciones de
idoneidad de hábitat comparando la distribución de la especie en el área cubierta por las
variables ecogeográficas. Este método selecciona, por un procedimiento de análisis factorial,
la información relevante del hábitat mediante la obtención de los factores de Marginalidad
(M) y de Especialización (E). La M representa la distancia ecológica entre el óptimo de la
especie y el promedio del hábitat en el área de referencia. Su valor global se sitúa
generalmente entre 0 y 1, aunque puede alcanzar valores mayores (Hirzel et al., 2002). Un
valor de M alto indica que la especie tiene requerimientos de hábitat que difieren de las
condiciones medias disponibles. Por otro lado, el factor de E, que varía entre 0 e infinito,
mide cuán estrecho es el rango de la variable para la especie. Su inverso, la Tolerancia, con
valores entre 0 y 1, permite interpretar la amplitud del nicho de la especie. Una especie que
presenta una Tolerancia cercana a 1 muestra un nicho más amplio, unos requerimientos
ecológicos menos estrictos que una especie con una Tolerancia cercana a 0. Aparte de los
valores globales de cada factor, el ENFA produce una tabla de varianza donde otorga un
coeficiente a cada variable ecogeográfica. Los coeficientes del factor M expresan la
marginalidad de la especie en cada variable ecogeográfica. La variable que presenta el valor
absoluto más alto será aquella que se aleja más de las condiciones medias disponibles en el
área de estudio. Valores negativos indican que la especie prefiere valores más bajos que la
media, valores positivos indican que la especie aparece en aquellas celdas que presentan
valores más altos que la media en el área de estudio. Respecto a los coeficientes del resto de
factores, la variable con el valor absoluto más alto será aquella que es más restrictiva para la
aparición de la especie. Con los factores extraídos, aquellos que expliquen la mayor parte de
la varianza, el programa construye el Mapa predictivo de HI.
En un principio se añadieron todos los mapas que representaban las características
ambientales. Mediante las diferentes opciones del software BioMapper, estos mapas fueron
verificados y se construyeron las matrices de covarianza y correlación entre todas las
variables ecogeográficas, junto con su correspondiente árbol de correlación. Para la
construcción del modelo se excluyeron aquellas variables altamente correlacionadas debido a
que resultan redundantes, no aportan información adicional al modelo y pueden provocar
problemas de multicolinearidad. Una vez seleccionadas las variables, se procedió a ejecutar
el ENFA. A partir del modelo obtenido se elaboró el MHI para Z. noltii. Este mapa se generó
en base al algoritmo de medianas (Hirzel et al., 2002; Hirzel y Arlettaz, 2003) incluyendo los
factores que explicaban la mayor parte de la varianza. Una vez realizado el mapa, se evaluó
su eficacia predictiva por medio del procedimiento de la validación cruzada aplicable en
BioMapper, siguiendo el método descrito por Boyce et al. (2002). Este tipo de validación
genera un intervalo de confianza (generalmente entre 0 y 1) acerca de la eficacia predictiva
del modelo (Skov et al., 2008). El índice de Boyce produce una curva de proporción entre lo
predicho y lo esperado, esta curva ofrece más información sobre la calidad del modelo:
robustez, resolución HI y desviación de la aleatoriedad. Dicha información permite la
reclasificación del mapa predicho en clases o intervalos significativos de HI.
Extrapolación del modelo al escenario de cambio climáticoExtrapolación del modelo al escenario de cambio climáticoExtrapolación del modelo al escenario de cambio climáticoExtrapolación del modelo al escenario de cambio climático
Se generó el modelo de extrapolación a partir del modelo de HI ya obtenido y por otro
lado se transformó el mapa de la variable altura topográfica con el objeto de simular el
escenario de cambio climático. Esta transformación consistió en disminuir las cotas del mapa
49 cm, ascenso del nivel medio del mar previsto para el Golfo de Vizcaya a lo largo del siglo
XXI (Caballero et al., 2009). Una vez transformado este mapa, se procedió a la extrapolación
del modelo en la misma área de estudio con la variable de altura topográfica modificada.
4.3.3. Resultados
El área total caracterizada dentro del estuario fue de 4,16 km², con 1.040.644 celdas.
De ella, el área ocupada por los polígonos de Z. noltii fue de 0,21 km², con 53.503 celdas
(Figura 18). La capa de puntos aleatorios creados sobre el área ocupada por los polígonos
presentó una superficie de 0,014 km², compuesta por 5.081 celdas.
Figura Figura Figura Figura 18181818. (a) Distribución de los polígonos de presencia de Zostera noltii en el estuario del Oka. Detalle de la zona de Arketas (b) y Murueta (c).
Análisis Factorial de Nicho EcológicoAnálisis Factorial de Nicho EcológicoAnálisis Factorial de Nicho EcológicoAnálisis Factorial de Nicho Ecológico
Para aplicar el ENFA de los 15 mapas ecogeográficos (generados a partir de cada una
de las variables antes mencionadas) se excluyeron aquellos que procedían de variables con
una correlación mayor al 0,724. Así, para la construcción del modelo se utilizaron finalmente
Tabla Tabla Tabla Tabla 3333. Listado de las variables ecogeográficas. Abreviaturas: IPT_fino, Índice de Posición Topográfico a escala fina; IPT_grueso, Índice de Posición Topográfico a escala gruesa.
Tipo de variableTipo de variableTipo de variableTipo de variable Variables ecogeográficas Variables ecogeográficas Variables ecogeográficas Variables ecogeográficas ComentariosComentariosComentariosComentarios
Características del sustrato Media granulométrica (phi) Utilizada en el análisis
Selección sedimentaria (phi) Utilizada en el análisis
Contenido en gravas (%) Utilizada en el análisis
Contenido en arenas (%) Eliminada (r= 0,973 con media)
Contenido en pelitas (%) Eliminada (r= 0,977 con media)
Contenido en materia orgánica (%) Utilizada en el análisis
Potencial redox (mV) Utilizada en el análisis
Características del terreno Altura (m) Utilizada en el análisis
Orientación Utilizada en el análisis
IPT_fino Utilizada en el análisis
IPT_grueso Eliminada (r= 1 con IPT_fino)
Pendiente (º) Utilizada en el análisis
Rugosidad Utilizada en el análisis
Condiciones hidrográficas Velocidad corriente de marea en llenante (m.s-1) Eliminada (r= 0,804 con vaciante)
Velocidad corriente de marea en vaciante (m.s-1) Utilizada en el análisis
El análisis ENFA realizado para Z. noltii en el estuario del Oka dio como resultado un
valor de M global de 1,008, lo que indica que el hábitat de la especie difiere
considerablemente de la media de las condiciones ambientales en el área de estudio; un valor
de E global de 3,455, valor que se traduce en un nicho ecológico estrecho; y por último, un
valor de Tolerancia global de 0,289, valor cercano a cero, lo cual confirma la especialización
de la especie, siendo ésta restrictiva en el rango de condiciones en las que vive.
A partir de los resultados obtenidos con el ENFA, se retuvieron cuatro factores para la
producción del MHI, los cuales explican el 90% de la especialización. El primer eje
seleccionado o de Marginalidad (que maximiza la diferencia absoluta entre la media global
ambiental y la media de la especie) explica el 19% de la variabilidad. Los otros tres factores
seleccionados (factores de Especialización) explican, respectivamente, el 36%, el 27% y el 8%
(Tabla 4).
Las variables con coeficientes de M mayores (que determinan la presencia de Z. noltii)
fueron: media granulométrica (0,71), potencial redox (-0,55), altura topográfica (0,26),
selección sedimentaria (0,24), pendiente (-0,15), contenido en gravas (-0,13) e IPT a escala
fina (0,11). Las variables ecogeográficas materia orgánica, orientación, rugosidad del terreno
y velocidad de la corriente vaciante presentan coeficientes de M inferiores, por lo que su
importancia es menor a la hora de determinar el nicho de la especie (Tabla 4).
Tabla Tabla Tabla Tabla 4444. . . . Varianza explicada por los cuatro primeros factores ecológicos y, entre paréntesis, valores de los coeficientes de las diferentes variables. Abreviaturas: Pot. redox, Potencial redox; IPT_fino, Índice de Posición Topográfica a escala fina; Mat.org., porcentaje de contenido en materia orgánica.
Con el objetivo de interpretar el valor de los coeficientes de las diferentes variables que
explican el factor de M, se comparó la distribución de los valores de las variables en aquellas
zonas donde la especie está presente y en el área de estudio (Tabla 5).
Tabla Tabla Tabla Tabla 5555. . . . Distribución de los valores de cada una de las variables ecogeográficas utilizadas para la elaboración del modelo de distribución potencial de Zostera noltii en el estuario del Oka. Para cada variable se calculó, mínimo, máximo, media y desviación estándar para las áreas de presencia de la especie y para el área de estudio. Abreviaturas: D.E., desviación estándar; IPT_fino: Índice de Posición Topográfico a escala fina.
Áreas de presencia Áreas de presencia Áreas de presencia Áreas de presencia Área de estudioÁrea de estudioÁrea de estudioÁrea de estudio
Mínimo Máximo Media Media Media Media D.E. Mínimo Máximo Media Media Media Media D.E.
Figura Figura Figura Figura 19191919. (a) Mapa de Hábitat Idóneo actual para Zostera noltii en el estuario del Oka. Detalle de la zona de Arketas (b) y Murueta (c).
Extrapolación del modelo al escenario de ascenso del nivel medio del Extrapolación del modelo al escenario de ascenso del nivel medio del Extrapolación del modelo al escenario de ascenso del nivel medio del Extrapolación del modelo al escenario de ascenso del nivel medio del marmarmarmar
Como resultado de la extrapolación del modelo se obtuvo un nuevo MHI, reclasificado
también en tres intervalos de idoneidad (Figura 20): de 0 a 33, (80 % del área total; 808.477
celdas ocupadas); de 33 a 67 (13,8% del área total; 141.651 celdas ocupadas); y de 67 a 100
(7,2% del área total; 73.564 celdas ocupadas). Comparando los nuevos porcentajes de área
ocupada por cada intervalo con los obtenidos en el MHI actual, se predijo la modificación que
sufrirá el HI actual bajo el escenario de ascenso del nivel medio del mar. Así, para las zonas
que actualmente presentan un grado de idoneidad de hábitat bajo (intervalo de 0 a 33), las
nuevas condiciones implicarán un aumento leve de su superficie (de 77,5% a 80%); para las
zonas clasificadas con un grado de idoneidad medio el aumento será mucho mayor (de 8,0% a
13,8%); mientras que las zonas de grado de idoneidad alta se reducirán drásticamente
pasando de ocupar un 14,5% del área total a un 7,2% (Figura 20).
Figura Figura Figura Figura 20202020. (a) Mapa de Hábitat Idóneo para Zostera noltii bajo condiciones de ascenso del nivel medio del mar de 49 cm en el estuario del Oka. Detalle de la zona de Arketas (b) y Murueta (c).
Tal y como se observa en la Figura 20 (b), en la zona de Arketas desaparecen las áreas
de alto grado de idoneidad y la superficie ocupada por las áreas de grado de idoneidad medio
se reduce. En la Figura 20 (c) se aprecia el aumento de las áreas de grado de idoneidad medio
a ambos lados del canal. En la margen derecha del canal la disminución de las áreas de alta
idoneidad es casi completa, aunque cabe destacar que en la zona baja aparece una mancha de
este grado que en el MHI actual se clasificaba como de grado de idoneidad bajo. Por último,
en la margen izquierda del estuario se aprecia un desplazamiento del área de idoneidad alta
hacia el interior, lo que provoca una disminución de la superficie ocupada por este tipo de
hábitat, y permite un aumento de la superficie del hábitat de grado de idoneidad media.
Tras extraerse los 12 parches se seleccionaron tres puntos a una distancia aproximada
de 200 m desde la zona donante, que fueron identificados como hábitat idóneo para Zostera
noltii. En cada uno de ellos se trasplantaron 4 parches en una disposición "en cruz", que
fueron ubicados a una distancia aproximada de 2 parches con respecto al punto central de la
zona (Figura 22). El traslado de los parches (dentro de cajas de madera) desde el punto
donante hasta los puntos receptores se hizo a pie usando unas planchas ligeras a modo de
trineo (Figuras 23, 24 y 25). Los puntos receptores fueron georreferenciados y señalizados
mediante estacas de madera para facilitar su localización en las posteriores visitas de
seguimiento, cuyo fin es evaluar el desarrollo de los ejemplares trasplantados y el nivel de
éxito de la actuación ejecutada.
� N
4 2
3
1
Canal Canal Canal Canal principal principal principal principal del rdel rdel rdel rííííoooo
� N
� N
4 2
3
1
Canal Canal Canal Canal principal principal principal principal del rdel rdel rdel rííííoooo
Figura 22.Figura 22.Figura 22.Figura 22. Disposición y distancia relativa entre parches (en verde) con las estacas de señalización de madera (en naranja) utilizadas en cada uno de los puntos. Cada uno de los 4 parches fue numerado (1, 2, 3 y 4) teniendo en cuenta su orientación, siguiendo la misma pauta en cada uno de los tres puntos (OK1, OK2 y OK3).
Figura 2Figura 2Figura 2Figura 23333.... Diversas fases de la extracción de los parches: puede verse la mezcla de Zostera noltii con el alga Gracilaria sp., parches extraídos y preparados para ser transportados, y transporte de los parches mediante arrastre.
Figura Figura Figura Figura 22224444.... Estado de la zona donante tras la extracción de los 12 parches. Colocando las estacas de referencia de madera para posteriores visitas.
Como consecuencia de esta actuación no se generaron residuos y lo único que se
introdujo en el medio fueron unas estacas de madera de unos 20 cm de longitud (2 por
parche) y unas tablillas de madera fácilmente biodegradables que quedaron enterradas en el
fondo de cada parche.
Esta actuación requerirá visitas posteriores, de periodicidad mensual-bimensual, que
consistirán en una inspección de la zona, con toma de fotografías de los parches
trasplantados y de la zona donante, con el fin de poder estimar el desarrollo de los ejemplares
trasplantados y la evolución de la recuperación de la mancha donante de Zostera.
Figura 2Figura 2Figura 2Figura 26666.... Zostera noltii en la zona donante, una muestra de sedimento, medida del potencial redox en el sedimento y tamizado de la muestra para la extracción de la macrofauna.
4.4.2. Caracterización de las zonas donante y recep tora
A continuación se exponen los resultados obtenidos el 30 de marzo de 2010 en los
puntos donante (se tomaron dos réplicas al inicio y final de la zanja, OK-D1 y OK-D2) y
receptores (OK1, OK2 y OK3) para facilitar la comparación y visualizar mejor sus similitudes
y diferencias.
Tanto los datos de la Tabla 6 como los de la Figura 27 confirman la gran similitud
existente entre las zonas en lo referente a la composición sedimentológica. La estación
receptora OK3 es la más distinta al resto: muestra una composición granulométrica
claramente más arenosa que el resto, lo cual no se considera un impedimento para que
Zostera noltii pueda sobrevivir, ya que esta planta presenta poblaciones tanto en sedimentos
más arenosos como más fangosos que los de OK-D (Garmendia et al., 2008). Aunque se ha
visto en Butroe que las plantas en zona arenosa están sujetas a la erosión por riadas, en esta
zona, en principio, no parece que esto pueda ser un problema, ya que el Oka tiene un caudal
mucho menor que el Butroe y, a diferencia de las zonas de trasplante del Butroe, las del Oka
están ubicadas lejos del canal principal del río. La estación receptora OK3 también se
distingue de las demás en su relación C/N, contenido en Cu, Cr y Mn. En general, las
características sedimentológicas de las zonas donante y las zonas receptoras son muy
similares, siendo quizás lo más destacable el mayor contenido en Pb de la zona donante.
Tabla Tabla Tabla Tabla 6666.... Parámetros sedimentológicos correspondientes a las zonas donante (OK-D1 y OK-D2) y receptoras (OK1, OK2 y OK3) del trasplante de Zostera noltii: Gravas, Arenas y Finos; contenido en materia orgánica (m.o.); temperatura y potencial redox (Eh); Compuestos orgánicos (COP: carbono orgánico particulado; NOP: nitrógeno orgánico particulado; C/N: relación carbono/nitrógeno) y metales en la fracción inferior a 0,063 mm del sedimento.
Figura 2Figura 2Figura 2Figura 27777.... Composición granulométrica de las zonas donante (OK-D1 y OK-D2) y receptoras (OK1, OK2 y OK3) del trasplante: porcentajes de cada fracción granulométrica y porcentajes acumulados.
Tabla Tabla Tabla Tabla 7777.... Parámetros estructurales de la fauna en las zonas donantes (OK-D1 y OK-D2) y receptoras (OK1, OK2 y OK3) del trasplante: Riqueza Específica (RE) en número de especies, Abundancia (Abund) en número de individuos en la muestra, Densidad, Diversidad de Shannon (H) y Equitatividad de Pielou (J).
Tabla Tabla Tabla Tabla 8888.... Listado de la fauna encontrada en diferentes muestras correspondientes a las zonas donantes (OK-D1 y OK-D2) y receptoras (OK1, OK2 y O3) del trasplante de Zostera noltii. En naranja se resaltan las especies coincidentes.
Figura 29.Figura 29.Figura 29.Figura 29. Valores promedio de la temperatura (T) registrada dentro (D) y fuera (F) de los parches ubicados en las estaciones de trasplante OK1 y OK2 en el estuario del OKA.
Figura 29 (Continuación).Figura 29 (Continuación).Figura 29 (Continuación).Figura 29 (Continuación). Valores promedio de la temperatura (T) registrada dentro (D) y fuera (F) de los parches de la estación de trasplante OK3 y de la zona donante (OKD) en el estuario del OKA.
Figura Figura Figura Figura 30303030.... Valores promedio del potencial de óxido-reducción (Eh) registrado dentro (D) y fuera (F) de los parches de las estaciones de trasplante OK1 y OK2 en el estuario del OKA.
FigurFigurFigurFigura 30 (Continuación).a 30 (Continuación).a 30 (Continuación).a 30 (Continuación). Valores promedio del potencial de óxido-reducción (Eh) registrado dentro (D) y fuera (F) de los parches de la estación de trasplante OK3 y de la zona donante (OKD) en el estuario del OKA.
En este apartado se aborda el trabajo correspondiente a la evaluación del impacto
generado en el medio físico de la zona donante de los esquejes trasplantados y del tiempo de
recuperación necesario tras este tipo de actuaciones.
4.5.1. Metodología
La extracción de los 12 parches fue realizada de manera ininterrumpida y siguiendo
una misma dirección para alterar la menor superficie posible. Como consecuencia de ello, se
provocó una zanja de unos 15 cm de profundidad y 2,5 m de largo por 1 m de ancho
(superficie total de parche útil extraído de 1,23 m²) y la consecuente alteración del sedimento
a causa del pisoteo. Al finalizar la extracción, se introdujeron unas tablillas de madera para
señalizar los bordes de la zanja y facilitar su visualización en las posteriores visitas de
seguimiento. Así mismo, se introdujeron 4 varas verticales (alineadas por pares) para
realizar las medidas de la progresión de la sedimentación o relleno natural de la zanja
originada (Figura 32). Estas varas sobresalen 15 cm de la superficie del sedimento y
permiten el acondicionamiento de una regla para poder dibujar el perfil de la zona a lo largo
de dos transectos.
Figura 32.Figura 32.Figura 32.Figura 32. Esquema de la extracción de los esquejes en la zona donante: cada rectángulo azul corresponde a un parche, los puntos naranjas son las varillas verticales fijas que servirán de apoyo a las reglas de medición de la profundidad de la zanja (que se llevarán en cada visita): perfiles OK-D1 y OK-D2.
En la figura 33 puede verse el estado en el que quedó la zona donante tras la
extracción de los esquejes para el trasplante.
Figura 3Figura 3Figura 3Figura 33333.... Zona donante tras la extracción de los esquejes, con las tablillas de señalización marcando el perímetro de la zanja. La vara blanca de la derecha mide 1 m de longitud.
Con la finalidad de evaluar el tiempo necesario para la recuperación de esta
perturbación, se planificó realizar visitas bimensuales para inspeccionar visualmente el
estado de la zona alterada, obtener fotografías, tomar medidas de la densidad de la población
natural de Zostera noltii y de la profundidad de la zanja.
Figura 34.Figura 34.Figura 34.Figura 34. Regla apoyada sobre las varillas preparada para la medición de la profundidad de la zanja.
4.5.2. Resultados
Las visitas realizadas hasta el momento han dado como resultado las fotografías de la
figura 35. En ellas se observa la rápida ocupación de la zanja por el alga Ulva sp.
Figura 3Figura 3Figura 3Figura 35555.... Estado de la zona donante durante las primeras 4 visitas del seguimiento: abril, junio, agosto y octubre de 2010.
Los perfiles OK-D1 y OK-D2 (Figura 36) muestran variaciones como consecuencia del
relleno de la zanja que se está produciendo de manera natural. A su vez, los bordes de la
zanja han sufrido cierto allanamiento por efecto de la erosión y derrumbe de las paredes de la
zanja.
OK-D1
0
5
10
15
20
25
30
0 10 20 30 40 50 60 70 80 90 100 110 120 130 140
cm
Pro
f. (c
m)
27/04/2010 28/06/2010 26/08/2010 26/10/2010
OK-D2
0
5
10
15
20
25
30
0 10 20 30 40 50 60 70 80 90 100 110 120 130 140
cm
Pro
f. (c
m)
27/04/2010 28/06/2010 26/08/2010 26/10/2010
Figura Figura Figura Figura 36363636.... Evolución de cada uno de los perfiles o profundidades correspondientes a cada sección de la zanja provocada (OK-D1 y OK-D2) en la zona donante.
(http://www.ccmar.ualg.pt/biomares/index_en.html), y como colofón a los cuatro años de su
duración, a finales de septiembre de 2010 se organizó un taller de trabajo en Portugal con el
objetivo de reunir a los distintos investigadores, principalmente europeos, que trabajan con
praderas marinas. El contenido del taller giró fundamentalmente en torno a la grave
situación actual de las praderas marinas europeas, a la preocupación por la tendencia
negativa que éstas están mostrando en cuanto a su superficie de ocupación en los fondos
marinos, y los grandes esfuerzos que, durante los últimos años, se vienen realizando en
distintos países para frenar su declive y contribuir a su recuperación. Este taller se celebró
durante los días 30 de septiembre y 1 de octubre de 2010 en Portinho da Arrábida (Portugal)
bajo el título “I European Seagrass ReI European Seagrass ReI European Seagrass ReI European Seagrass Restoration Workshopstoration Workshopstoration Workshopstoration Workshop”.
Dado el interés que presentaba esta reunión, y la oportunidad que ofrecía para
establecer contactos con otros expertos en la materia, representantes de AZTI-Tecnalia
acudieron a dicho taller donde expusieron sus experiencias e informaron sobre la situación de
las praderas marinas en el País Vasco mediante la presentación de la comunicación “Zostera Zostera Zostera Zostera
noltii noltii noltii noltii restoration essays in the estuaries of the Basque Countryrestoration essays in the estuaries of the Basque Countryrestoration essays in the estuaries of the Basque Countryrestoration essays in the estuaries of the Basque Country”.
En el mencionado taller de intercambio de experiencias participaron unos 30
investigadores y expertos procedentes de varios países y pertenecientes a diferentes
organismos o instituciones: Portugal (Instituto da Conservação da Natureza e da
Biodiversidade-ICNB/PNA, Centro de Ciências do Mar do Algarve-CCMAR, Universidade de
Algarve, Universidade de Évora, IMAR), Estados Unidos (Cornell University Cooperative
Extension), España (Instituto Mediterráneo de Estudios Avanzados-IMEDEA, Consejo
Superior de Investigaciones Científicas-CSIC, Instituto de Física Interdisciplinar y Sistemas
Complejos-IFISC (Universidad de las Islas Baleares-UIB/CSIC), AZTI-Tecnalia, Instituto
Canario de Ciencias Marinas-ICCM, ECOs-Environmental Consulting St., Universidad de
Las Palmas, OCEANA), Francia (Groupement d'Intérêt Public pour la Réhabilitation de
l'étang de Berre-GIPREB), y Holanda (Radboud University Nijmegen).
La primera jornada se dedicó a la exposición de las distintas ponencias por parte de
cada uno de los grupos participantes, junto con sus respectivas rondas de preguntas y
discusión. La segunda jornada consistió en una sesión de debate y análisis de la situación
actual de las praderas marinas en Europa, del avance y estado de las investigaciones
relativas a la restauración y recuperación de dichas praderas, y del grado de éxito obtenido
con las diferentes alternativas metodológicas empleadas. Esta jornada concluyó con el
consenso y redacción de una serie de conclusiones y con la idea generalizada de que deberían
repetirse este tipo de reuniones con el objeto de aunar conocimientos y mejorar metodologías
en favor de la recuperación de las praderas marinas europeas.
4.6.2. Conclusiones del taller
En este taller se alcanzaron algunas conclusiones cuya divulgación, tanto en el entorno
de las distintas autoridades gestoras del dominio acuático como del público en general, fue
considerada muy interesante e incluso necesaria dada la delicada situación por la que
atraviesan las praderas marinas, especialmente en el ámbito europeo.
A continuación se presentan las conclusiones mencionadas que han servido de base
para realizar una comunicación en el 9th International Seagrass Biology Workshop (ISBW),
Restoration of coastal ecosystems in Southeastern Asian regions: “Seagrass decline;
requisites and successful restoration” celebrado en Trang (Tailandia) en noviembre de 2010
(Conclusions and recommendations from the I European Seagrass Restoration WorkshopConclusions and recommendations from the I European Seagrass Restoration WorkshopConclusions and recommendations from the I European Seagrass Restoration WorkshopConclusions and recommendations from the I European Seagrass Restoration Workshop. . . .
Cunha, A.H., Marbá, N., van Katwijk, M., Pickerell, C., Henriques, M., Bernard, G., Ferreira, Cunha, A.H., Marbá, N., van Katwijk, M., Pickerell, C., Henriques, M., Bernard, G., Ferreira, Cunha, A.H., Marbá, N., van Katwijk, M., Pickerell, C., Henriques, M., Bernard, G., Ferreira, Cunha, A.H., Marbá, N., van Katwijk, M., Pickerell, C., Henriques, M., Bernard, G., Ferreira,
En base a los resultados presentados en este taller se recomiendan las siguientes
acciones antes de comenzar con la ejecución de un proyecto de restauración:
1. Establecer unas metas y objetivos claros antes de iniciar la restauración.
2. Definir la metodología de seguimiento o supervisión y fijar los criterios de éxito
antes de comenzar.
3. Llevar a cabo convenios para asegurar una supervisión a largo plazo si ésta no forma
parte del proyecto inicial.
4. Realizar los esfuerzos necesarios para asegurar que las amenazas locales (por
ejemplo, la bioturbación, el herbivorismo, la hidrología, los movimientos del sedimento, los
impactos humanos, etc.), que puede que influyan en la restauración de las praderas marinas,
están identificadas y son bien conocidas antes de comenzar con el proyecto.
5. Considerar todo posible impacto basado en una revisión de la literatura y la medida
de las condiciones físicas de la zona, todo ello combinado con el conocimiento local.
6. Los proyectos de restauración solo deben llevarse a cabo cuando las amenazas que
provocan la regresión de las praderas se hayan eliminado.
7. Antes de implicarse en un proyecto de restauración a gran escala deben realizarse
ensayos o pruebas piloto a pequeña escala, aunque en algunos casos, la realización directa a
gran escala podría ser necesaria (van Katwijk et al., 2009).
2 COUNCIL REGULATION (EC) No 1967/2006 of 21 December 2006 concerning management measures for the sustainable exploitation of fishery resources in the Mediterranean Sea, amending Regulation (EEC) No 2847/93 and repealing Regulation (EC) No 1626/94.
14. Se deben publicar los resultados y compartir las experiencias, particularmente en
los comienzos de los nuevos proyectos. Sería bueno también entrevistar a los científicos y
expertos en restauración que ya realizaron pruebas de trasplante en los años setenta, como
Meinesz (Francia), Thorhaug (Florida, Caribe)… y/o estudiar su trabajo.
15. Casi todos los participantes muestran frustraciones ya sea en praderas naturales
que han sido alteradas o en recuperaciones naturales que han sido impedidas (pesca de
arrastre, marisqueo, recolección de cebo, actividades turísticas...). En parte debido a la
ausencia de aplicación de la ley, en parte debido al limitado estado de regulación o
protección, o incluso a la modificación del estado de protección cuando prevalece el interés
económico. Antes de empezar cualquier esfuerzo de restauración debe asegurarse la
identificación de todas estas trabas y presiones, su magnitud y frecuencia.
Importantes recomendaciones para la conservación de las praderas marinasImportantes recomendaciones para la conservación de las praderas marinasImportantes recomendaciones para la conservación de las praderas marinasImportantes recomendaciones para la conservación de las praderas marinas
1. Los esfuerzos para la conservación de las praderas existentes actualmente deben ser
una prioridad europea.
2. En ningún caso la restauración debe ser considerada como la primera alternativa
para la mitigación de los efectos en la planificación de los proyectos de desarrollo costeros. Al
trabajar en áreas donde existen praderas marinas, los gestores deben seguir una adecuada
secuencia de opciones de mitigación. Se debe evitar el daño a las praderas existentes a través
de la adecuada selección del lugar (es decir, trasladar el proyecto a un área que no dañe a la
pradera). Cuando sea imposible trasladarlo a otra área, se debe reducir la huella del impacto
para minimizar la perturbación sobre la pradera. Finalmente, cuando las dos primeras
alternativas no sean posibles, se podrían considerar acciones de restauración como
mitigación del impacto. En tal caso, e insistiendo en que debe ser el último recurso, la
restauración debería ser concebida en el siguiente orden: mismo lugar y mismo tipo, distinto
lugar y mismo tipo, distinto lugar y distinto tipo.
3. El éxito de la restauración es generalmente demasiado bajo como para justificarla
como una adecuada medida de compensación por actividades económicas. En caso de
intereses primordiales (como la seguridad) podría emplearse la restauración, pero con el
objetivo adicional de desarrollarla en base a una investigación científica y aprovechar los
resultados para ampliar el conocimiento sobre el funcionamiento de las praderas marinas.
Cada medio es único en sí y la conservación de todas las praderas marinas existentes es
crucial para mantener y restaurar la degradación de los mares.
4. Se debe dar prioridad a la potencial restauración natural.
5. Se debe reconocer que la dinámica de cada zona depende de fenómenos naturales y
las actuaciones deben enfocarse desde una perspectiva de paisaje y conectividad entre
poblaciones.
6. Se debe presionar a los gobiernos para catalogar las especies de las praderas
marinas como especies prioritarias o protegidas.
7. Se debe presionar a los gobiernos para reconsiderar las áreas de arrastre con el fin
de evitar el daño a las zonas de pradera y/o praderas históricas para permitir la
recuperación, así como todas las actividades humanas que causan la regresión de las
praderas europeas. La Directiva de la Estrategia Marina Europea3 de 2008 plantea que los
estados miembros tienen la obligación de lograr un “buen estado ambiental” para el año
2020, por tanto, todas las actividades humanas que tienen un impacto sobre los ecosistemas
marinos deben ser abordadas con el objetivo de preservar el medio marino.
8. Se deben contratar especialistas en relaciones públicas para desarrollar las acciones
de mercado sobre la conservación de las praderas marinas necesarias con el objeto de
alcanzar un adecuado intercambio de información y comunicación entre las partes
implicadas.
9. Se debe considerar e implicar a las ONGs en los proyectos.
10. Se debe evaluar el interés de desarrollar un proyecto INTERREG (estudios pilotos
de restauración) o LIFE (evaluación de la distribución histórica basada en el conocimiento
local, etc.) sobre praderas marinas europeas.
11. Se deben buscar fuentes alternativas de financiación para proyectos de
investigación, conservación y restauración, a saber, fondos privados de compañías o
instituciones interesadas en las áreas donde esos proyectos tendrían lugar (es decir, la
implicación de los agentes sociales en la financiación).
3 EC (2008). Directive 2008/56/EC of the European Parliament and of The Council of 17 June 2008 establishing a framework for community action in the field of marine environmental policy (Marine Strategy Framework Directive). Official Journal of the European Union. 25th of June of 2008. 164/19-40.
En el transcurso de los distintos trabajos que se han llevado a cabo a lo largo de este
año 2010 se han ido obteniendo una serie de resultados y conclusiones fundamentales para la
posterior toma de decisiones necesarias para ir avanzando en las tareas a ejecutar. Estas
conclusiones, tanto intermedias como finales, son las siguientes:
1.1.1.1. En base a un nuevo análisis genético con mayor cantidad de muestras, se confirma la
idoneidad del estuario destuario destuario destuario del Okael Okael Okael Oka como el más apropiado para actuar como donantedonantedonantedonante de
ejemplares de Zostera noltii para el estuario del Butroe. Asimismo, se confirma que las
poblaciones de esta especie en cada uno de los estuarios estudiados presenta una identidad
genética propia y diferente de los otros, con una mayor diferencia entre los estuarios más
alejados.
2.2.2.2. Los parchesparchesparchesparches trasplantados en los sedimentos arenosossedimentos arenosossedimentos arenosossedimentos arenosos del Butroe han desaparecido como
consecuencia del elevado hidrodinamismo provocado por una riadariadariadariada (junio 2010) que ha
erosionado el sedimento que servía de sustrato a los ejemplares de Zostera noltii y ha
cambiando sustancialmente la morfología del lugar. A la hora de analizar la idoneidad de
una zona receptora e identificar las posibles presiones existentes, se debe recabar
información histórica o plurianual y tener en cuenta también aquellos agentes
meteorológicos extremos, de baja frecuencia relativa y de carácter imprevisible.
3.3.3.3. Los parchesparchesparchesparches trasplantados en los sedimentos sedimentos sedimentos sedimentos fangosofangosofangosofangosossss del Butroe siguensiguensiguensiguen desarrollándose
lentamente tras 18 meses18 meses18 meses18 meses. Dada su ubicación más alejada del canal principal del río, sus
sedimentos no han sufrido el efecto devastador de la riada de junio. No obstante, no se
descarta el posible impacto de futuros episodios de elevado crecimiento de macroalgas que se
ha repetido en primavera de los últimos años y que refleja la existencia de eutrofización de
las aguas.
4.4.4.4. La interpretación definitiva del resultado de un trasplante debe realizarse tras un
prolongadoprolongadoprolongadoprolongado periodo de seguimientoseguimientoseguimientoseguimiento, que debe incluir al menos dos épocas invernales. El buen
crecimiento y desarrollo inicial del trasplante no garantiza su éxito. A pesar de un
espectacular desarrollo de los trasplantes durante los primeros meses o época estival,
cualquier evento puntual puede afectar la salud y supervivencia de los parches e invertir las
tendencias observadas hasta entonces.
5.5.5.5. El Análisis Factorial de Nicho Ecológico (ENFA), que produce modelos de hábitat idóneo a
partir de la caracterización del nicho ecológico de la especie, ha identificado como principales principales principales principales
variables ambientalesvariables ambientalesvariables ambientalesvariables ambientales que determinan la distribución de Z. noltii en el estuario del Oka:
media granulométrica, potencial redox, altura topográfica, selección del sedimento,
pendiente, porcentaje de contenido en gravas del sedimento e Índice Topográfico de Posición
a escala fina.
6.6.6.6. Considerando los escenarios futuros en cuanto a ascenso del nivel del mar en la costa
vasca previsto para finales de este siglo XXI (49 cm), existe una elevada vulnerabilidadelevada vulnerabilidadelevada vulnerabilidadelevada vulnerabilidad de
estas praderas, ya que se estima que la superficie actual del hábitat más idóneo de las
praderas de Z. noltii en el estuario del Oka se reducirá hasta un 50%se reducirá hasta un 50%se reducirá hasta un 50%se reducirá hasta un 50%; se debe tener presente
en todo momento que este valor está sobreestimado, ya que no contempla la probable
creación de nuevos hábitats adecuados para la especie durante el periodo que transcurra
hasta que se produzca dicho ascenso de nivel del mar. Esto hace que sea importante
emprender acciones de conservación de las praderas existentes y restauración de los
estuarios donde haya desaparecido en el pasado.
7777.... Se observa un desarrollo diferencialdesarrollo diferencialdesarrollo diferencialdesarrollo diferencial de los trasplantestrasplantestrasplantestrasplantes realizados en el año 2010201020102010 en zonas
identificadas por el modelo como idóneas dentro del estuario del Oka. Uno de los trasplantes
parece que va a desaparecer, otro se mantiene por ahora y el tercero parece que se
desarrollará correctamente. No obstante, tal y como se comenta en una conclusión anterior,
se deben esperar algunos meses más para obtener unas conclusiones más definitivas y
ajustadas a la realidad.
8888.... A falta de una confirmación estadística, la presencia de Zostera noltii influyeinfluyeinfluyeinfluye en las
característicascaracterísticascaracterísticascaracterísticas físico-químicas del sedimentosedimentosedimentosedimento: en la zona ocupada por esta planta la
temperatura es menor y el potencial de óxido-reducción registra valores más negativos que
en la zona desprovista de ella.
9.9.9.9. La recuperación natural de la alteración provocada por la actuación de trasplante
requiere más de 7 meses, ya que, tras ese periodo, aún permanecen señales de la zanja
originada tras la extracción de los parches. Debe continuarse con el seguimiento de dicha
zona para determinar el momento de su plena recuperación.
10101010.... Se deben continuar con las reuniones con otros gruposreuniones con otros gruposreuniones con otros gruposreuniones con otros grupos de investigación para
intercambiar experiencias, adquirir nuevos conocimientos en beneficio de un mayor
aprovechamiento del esfuerzo realizado y generar un mayor respaldo científico en la defensa
Auby, I. y P.J. Labourg, 1996. Seasonal dynamics of Zostera noltii Hornem. in the Bay of Arcachon (France). Journal of Sea Research, 35(4): 269-277.
Baeta, A., I. Valiela, F. Rossi, R. Pinto, P. Richard, N. Niquil y J.C. Marques, 2009. Eutrophication and trophic structure in response to the presence of the eelgrass Zostera noltii. Marine Biology, 156: 2107-2120.
Balestri, E., L. Piazzi, y F. Cinelli, 1998. Survival and growth of transplanted and natural seedlings of Posidonia oceanica (L.) Delile in a damaged coastal area. Journal of Experimental Marine Biology and Ecology, 228: 209-225.
Belkhir, K., Borsa, P., Chikhi, L., Raufaste, N. y P. Bonhomme, 2004. 1996-2004 GENETIX 4.05, logiciel sous Windows TM pour la génétique des populations. Laboratoire Génome, Populations, Interactions, CNRS UMR 5171, Université de Montpellier II, Montpellier (France). http://www.genetix.univ-montp2.fr/genetix/genetix.htm
Boese, B.L. y B.D. Robbins, 2008. Effects of erosion and macroalgae on intertidal eelgrass (Zostera marina) in a northeastern Pacific estuary (USA). Botanica Marina, 51: 247-257.
Boese, B.L., J.E. Kaldy, P.J. Clinton, P.M. Eldridge y C.L. Folger, 2009. Recolonization of intertial Zostera marina L. (eelgrass) following experimental shoot removal. Journal of Experimental Marine Biology and Ecology, 374: 69-77.
Borja, A., F. Aguirrezabalaga, J. Martínez, J.C. Sola, L. García-Arberas y J.M. Gorostiaga, 2004. Benthic communities, biogeography and resources management. In: Borja, A., Collins, M. (Eds.), Oceanography and Marine Environment of the Basque Country. Elsevier Oceanography Series, 70: 455-492.
Borja, A., I. Galparsoro, O. Solaun, I. Muxika, E. M. Tello, A. Uriarte y V. Valencia, 2006. The European Water Framework Directive and the DPSIR, a methodological approach to assess the risk of failing to achieve good ecological status. Estuarine, Coastal and Shelf Science, 66: 84-96.
Borja, A., J. Bald, M.J. Belzunce, J. Franco, J.M. Garmendia, J. Larreta, I. Muxika, M. Revilla, G. Rodríguez, O. Solaun, A. Uriarte, V. Valencia, I. Adarraga, F. Aguirrezabalaga, I. Cruz, A. Laza, M.A. Marquiegui, J. Martínez, E. Orive, J.Mª Ruiz, S. Seoane, J.C. Sola, J.Mª Trigueros y A. Manzanos, 2008. Red de seguimiento del estado ecológico de las aguas de transición y costeras de la Comunidad Autónoma del País Vasco. Informe de AZTI-Tecnalia para Agencia Vasca del Agua, Gobierno Vasco. 20 Tomos, 624 pp.
Borum, J., C.M. Duarte, D. Krause-Jensen y T.M. Greve (Eds.), 2004. European seagrasses: an introduction to monitoring and management. Monitoring and Managing of European seagrasses (M&MS) EU project. 88 pp. http://www.seagrasses.org
Bos, A.R. y M.M. van Katwijk, 2007. Planting density, hydrodynamics exposure and mussel beds affect survival of transplanted intertidal eelgrass. Marine Ecology Progress Series, 336: 121-129.
Bos, A.R., N. Dankers, A.H. Groeneweg, D.C.R. Hermus, Z. Jager, D.J. de Jong, T. Smit, J. de Vlas, M. van Wieringen y M.M.M. van Katwijk, 2005. Eelgrass (Zostera marina L.) in
the western Wadden Sea: monitoring, habitat suitability model, transplantations and communication. En: Herrier, J.L., J. Mees, A. Salman, J. Syes, H. van Nieuwenhuyse y I. Dobbelaere (Eds.). Proceedings “Dunes and Estuaries 2005” - International Conference on Nature Restoration Practices in European Coastal Habitats, VLIZ Special Publication, 19: 95:109.
Boyce, M., P.R. Vernier, S.E. Nielsen y F.K. Schmiegelow, 2002. Evaluating resource selection functions, Ecological Modelling, 157: 281–300.
Braunisch, V., K. Bollmann, R.F. Graf y A.H. Hirzel, 2008. Living on the edge- Modelling habitat suitability for species at the edge of their fundamental niche. Ecological Modelling, 214 (2-4): 153-167.
Brun, F.G., J.L. Pérez-Llorens, I. Hernández y J.J. Vergara, 2003a. Patch distribution and within-patch dynamics of the seagrass Zostera noltii Hornem. in Los Toruños salt-marcs, Cádiz Bay, Natural Park, Spain. Botanica Marina, 46: 513-524.
Brun, F.G., J.J. Vergara, G. Navarro, I. Hernández y J.L. Pérez-Llorens, 2003b. Effect of shading by Ulva rigida canopies on growth and carbon balance of the seagrass Zostera noltii. Marine Ecology Progress Series, 265: 85-96.
Burkholder, J.A.M., D.A. Tomasko y B.W. Touchette, 2007. Seagrasses and eutrophication. Journal of Experimental Marine Biology and Ecology, 350: 46-72.
Cabaço, S. y R. Santos, 2007. Effects of burial and erosiono on the seagrass Zostera noltii. Journal of Experimental Marine Biology and Ecology, 340: 204-212.
Cabaço, S., R. Machás y R. Santos, 2009. Individual and population plasticity of the seagrass Zostera noltii along a vertical intertidal gradient. Estuarine, Coastal and Shelf Science, 82: 301-308.
Cabaço, S., R. Machás, V. Vieira y R. Santos, 2008a. Impacts of urban wastewater discharge on seagrass meadows (Zostera noltii). Estuarine, Coastal and Shelf Science, 78: 1-13.
Cabaço, S., R. Santos y C.M. Duarte, 2008b. The impact of sediment burial and erosion on seagrasses: a review. Estuarine, Coastal and Shelf Science, 79: 354-366.
Caballero, A., M. Marcos y G. Chust, 2009. Regional sea level change: projections and impacts in the Basque coast. 33rd International Symposium on Remote Sensing of Environment (ISRSE), 4-8 May 2009, Stresa (Italy).
Cardoso, P.G., D. Raffaelli y M.A. Pardal, 2008. The impact of extreme weather events on the seagrass Zostera noltii and related Hydrobia ulvae population. Marine Pollution Bulletin, 56(3): 483-492.
Cardoso, P.G., S. Leston, T.F. Grilo, M.D. Bordalo, D. Crespo, D. Raffaelli y M.A. Pardal, 2010. Implications of nutrient decline in the seagrass ecosystem success. Marine Pollution Bulletin, 60: 601-608.
Castro, R., A. Uriarte, A. Martínez de Murguía, A. Borja, 2004. Biodiversity and conservation of wildlife and natural habitats. Borja, A. and Collins, M. (Eds.). Oceanography and Marine Environment of the Basque Country, Elsevier Oceanography Series, 70: 531-547.
Chust, G., A. Borja, P. Liria, I. Galparsoro, M. Marcos, A. Caballero y R. Castro, 2009a. Human impacts overwhelm the effects of sea-level rise on Basque coastal habitats (N Spain) between 1954 and 2004. Estuarine, Coastal and Shelf Science, 84(4): 453-462.
Chust, G., A. Borja, J. Franco, P. Liria, I. Galparsoro y A. Caballero, 2009b. Aplicación de sensores de alta resolución para la evaluación de cambios recientes y escenarios futuros del cambio climático en los hábitats de la Reserva de la Biosfera de Urdaibai (Vizcaya). Informe relativo a las actividades realizadas durante 2008. Informe inédito
elaborado por AZTI-Tecnalia para la Secretaría de Estado de Medio Rural y Agua (Ministerio de Medio Ambiente y Medio Rural y Marino), 53 pp.
Coyer, J.A., T.B.H. Reusch, W.T. Stam, E.S Serrao, G. Pearson, G. Procaccini y J.L. Olsen, 2004. Characterization of microsatellite loci in the dwarf eelgrass Zostera noltii (Zosteraceae) and cross–reactivity with Z. japonica. Moolecular Ecology Notes, 4: 497-499.
Cunha, A.H., R.P. Santos, A.P. Gaspar y M.F. Barrios, 2005. Seagrass landscape-scale changes in response to disturbance created by the dynamics of barrier-islands: a case study from Ria Formosa (Southern Portugal). Estuarine, Coastal and Shelf Science, 64: 636-644.
Curiel, D., F. Scarton, A. Rismondo y M. Marzocchi, 2003. Trasplanting seagrasses in the Lagoon of Venice: results and perspectives. En: E. Özhan (Ed.), Proceedings of the Sixth International Conference on the Mediterranean Coastal Environment, MEDCOAST 03, Ravena, Italy, Vol. 2: 853-864.
de Jonge, V.N., D.J. de Jong y M.M. van Katwijk, 2000. Policy plans and management measures to restore eelgrass (Zostera marina L.) in the Dutch Wadden Sea. Helgoland Marine Research, 54: 151-158.
Dolbeth, M., P.G. Cardoso, S.M. Ferreira, T. Verdelhos, D. Raffaelli y M.A. Pardal, 2007. Anthropogenic and natural distubance effects on a macrobenthic estuarine community over a 10-year period. Marine Pollution Bulletin, 54(5): 576-585.
Duarte, C.M., 2002. The future of seagrass meadows. Environmental Conservation, 29(2): 192–206.
Duarte, C.M., J. Terrados, N.S.W. Agawin, M.D. Fortes, S. Bach y W.J. Kenworthy, 1997. Response of a mixed Philippine seagrass meadow to experimental burial. Marine Ecology Progress Series, 147: 285-294.
Erftemeijer, P.L.A. y R.R.R. Lewis, 2006. Environmental impacts of dredging on seagrasses: a review. Marine Pollution Bulletin, 52: 1553-1572.
Estrada-Peña, A. y J.M. Venzal, 2007. Climate niches of tick species in the Mediterranean region: modelling of occurrence data, distributional constrains, and impact of climate change. Journal of Medical Entomology. 44, 1130-1138.
Foden, J. y D.P. Brazier, 2007. Angiosperms (seagrass) within the EU water framework directive: A UK perspective. Marine Pollution Bulletin, 55: 181-195.
Foden, J. y D.J. de Jong, 2007. Assessment metrics for littoral seagrass under the European Water Framework Directive; outcomes of UK intercalibration with the Netherlands. Hydrobiologia, 579: 187-197.
Fonseca, M.S. y W.J. Kenworthy, 1987. Effects of current on photosynthesis and distribution of seagrasses. Aquatic Botany, 27: 59-78.
Fonseca, M.S., W.J. Kenworthy, B.E. Julius, S. Shutler y S. Fluke, 2002. Seagrasses. En: Perrow, M.R. y A.J. Davy (Eds.), Handbook of Ecological Restoration. University Press, Cambridge, pp: 149-170.
Galparsoro, I., A. Borja, J. Bald, P. Liria y G. Chust, 2008. Predicting habitat suitability for European Lobster (Homarus gammarus), on the Basque continental shelf (Bay of Biscay), using Ecological-Niche Factor Analysis. Ecological modelling, 220: 556-567.
García, P., E. Zapico y A. Colubi, 2009. An angiosperm quality index (AQI) for Cantabrian estuaries. Ecological Indicators, 9: 856-865.
Garmendia, J.M., A. Borja y J. Franco, 2008. Trabajos de restauración de Zostera noltii en la costa vasca. Informe elaborado por AZTI-Tecnalia para la Agencia Vasca del Agua. 94 pp.
Garmendia, J.M., Á. Borja y J. Franco, 2010. Trabajos de restauración de Zostera noltii en la costa vasca: actividades de 2009. Informe elaborado por AZTI-Tecnalia para la Agencia Vasca del Agua – Uraren Euskal Agentzia, Departamento de Medio Ambiente, Planificación Territorial, Agricultura y Pesca del Gobierno Vasco, 95 pp.
Golden, R.R., K.E. Busch, L.P. Karrh, T.A. Parham, M.J. Lewandowski y M.D. Naylor, 2010. Large-scale Zostera marina (eelgrass) restoration in Chesapeake Bay, Maryland, USA. Part II: a comparison of restoration methods in the Patuxet and Potomac Rivers. Restoration Ecology, 18(4): 501-513.
González, M., 1994. Un modelo numérico en Elementos Finitos para la corriente inducida por la marea. Aplicaciones al Estrecho de Gibraltar. Tesina de especialidad. Escuela Técnica de Ingenieros de Caminos, Canales y Puertos de Barcelona. UPC.
González, M., M.A. García, M. Espino y A.S. Arcilla, 1995. Un modelo numérico en Elementos Finitos para la corriente inducida por la marea. Aplicaciones al Estrecho de Gibraltar. Revista Internacional de Métodos Numéricos para Cálculo y Diseño en Ingeniería, 11(3): 383-400.
Harwell M.C. y R.J. Orth, 2002. Long-distance dispersal potential in a marine macrophyte. Ecology, 83(12): 3319-3330.
Hauxwell, J., J. Cebrian, C. Furlong y I. Valiela, 2001. Macroalgal canopies contribute to eelgrass (Zostera marina) decline in temperate estuarine ecosystems. Ecology, 82(4): 1007-1022.
Hawkins, S. J., J. R. Allen, S. Bray, 1999. Restoration of temperate marine and coastal ecosystems: nudging nature. Aquatic Conservation-Marine and Freshwater Ecosystems, 9: 23-46.
Hirzel, A.H. y R. Arlettaz, 2003. Modelling habitat suitability for complex species distributions by the environmental-distance geometric mean. Environmental Management, 32: 614-623.
Hirzel, A.H., J. Hausser, D. Chessel y N. Perrin, 2002. Ecological Niche Factor Analysis: How to compute habitat suitability maps without absence data? Ecology, 83(7): 2027-2036.
Jenness, J., 2006. Topographic Position Index (tpi_jen.avx) extension for ArcView 3.x 1.3a. Jenness Enterprises. www.jennessent.com/arcview/tpi.htm
Kairis, P.A. y J.M. Rybczyk, 2009. Sea level rise and eelgrass (Zostera marina) production: A spatially explicit relative elevation model for Padilla Bay, WA. Ecological Modelling. DOI:10.1016/j.ecolmodel.2009.01.025.
Källström, B., A. Nyqvist, P. Aberg, M. Bodin y C. André, 2008. Seed rafting as dispersal strategy for eelgrass (Zostera marina). Aquatic Botany, 88: 148-153.
Krause Jensen, D., T.M. Greve, K. Nielsen, 2005. Eelgrass as a bioindicator under the European Water Framework Directive. Water Resources Management, 19: 63-75.
Lee, K.S., Park, J.I. Kim, Y.K. Park, S.R. y J.H. Kim, 2007. Recolonization of Zostera marina following destruction caused by a red tide algal bloom: the role of new shoot recruitment from seed banks. Marine Ecology Progress Series, 342: 105-115.
Leschen, A.S., K.H. Ford y N.T. Evans, 2010. Successful eelgrass (Zostera marina) restoration in a formerly eutrophic estuary (Boston Harbor) supports the use of a multifaceted watershed approach to mitigating eelgrass loss. Estuaries and Coasts, 33: 1340-1354.
Martins, I., J.M. Neto, M.G. Fontes, J.C. Marques y M.A. Pardal, 2005. Seasonal variation in short-term survival of Zostera noltii transplants in a declining meadow in Portugal. Aquatic Botany, 82: 132-142.
Meinesz, A., G. Caye, F. Locques y S. Macaux, 1990. Analyse bibliographique sur la culture des phanerogames marines. Posidonia Newsletter, 3: 1-67.
Molenaar, H. y A. Meinesz, 1995. Vegetative reproduction in Posidonia oceanica: survival and development of transplanted cuttings according to different spacings, arrangements and substrates. Botanica Marina, 38: 313-322.
Moore, K.A. y F.T. Short, 2006. Zostera: biology, ecology and management. Cap. 16, pp.361-386. In. Larkum, A.W.D., Orth, R.J. & Duarte, C.M. (Eds), Seagrasses: biology, ecology and conservation. Springer, Dordrecht, The Netherlands, 648 pp.
Oetjen, K., S. Ferber, I. Dankert y T. Reusch, 2010. New evidence for habitat-specifi selection in Wadden Sea Zostera marina populations revealed by genome scanning using SNP and microsatellite markers. Marine Biology, 157: 81-89.
Orth, R.J., S.R. Marion, K.A. Moore y D.J. Wilcox, 2010. Eelgrass (Zostera marina L.) in the Chesapeae Bay region of Mid-Atlantic coast of the USA: challenges in conservation and restoration. Estuaries and Coasts, 33: 139-150.
Otxoa, E., T. Oyanarte, S. Patino, J. Valencia, E. Miguel, A. Prieto y J. Elorza, 2007. Cartografi ado a escala 1:5000 de las poblaciones vizcaínas de plantas incluidas en el catálogo vasco de especies amenazadas. Departamento de Botánica, Sociedad de Ciencias Naturales de Sestao, 71 pp.
Oviedo, L., 2007. Dolphin sympatric ecology in tropical ford: habitat bathymetry and topography as a strategy to coexist. Journal of the Marine Biological Association of the UK., 87: 1327-1335.
Park, J.I. y K.S. Lee, 2007. Site-specific success of three transplanting methods and the effect of planting time on the establishment of Zostera marina transplants. Marine Pollution Bulletin, 54: 1238-1248.
Peralta, G., F.G. Brun, J.L. Pérez-Lloréns y T.J. Bouma, 2006. Direct effects of current velocity on the growth, morphometry and architecture of seagrasses: a case study on Zostera noltii. Marine Ecology Progress Series, 327: 135-142.
Philippart, C.J.M., 1995. Seasonal variation in growth and biomass of an intertidal Zostera noltii stand in the Dutch Wadden Sea. Netherlands Journal of Sea Research, 28: 287-299.
Piry, S., A. Alapetit, J.M. Cornuet, D. Paetkau, L. Baudouin y A. Estoup, 2004. GeneClass 2: a software for genetic assignment and first generation migrants detection. Journal of Heredity, 95: 536-539.
Praca, E. y A. Gannier, 2008. Ecological niches of three teuthophageous odontocetes in the northwestern Mediterranean Sea. Ocean Science, 4(1): 49-59.
Rannala, B. y J.L. Mountain, 1997. Detecting immigration by using multilocus genotypes. Proceedings of the National Academy of Sciences USA 94: 9197-9201.
Rodríguez-Salinas, P., R. Riosmena-Rodríguez, G. Hinojosa-Arango y R. Muñiz-Salazar, 2010. Restoration experiment of Zostera marina L. in a subtropical lagoon. Ecological Engineering, 36: 12-18.
Schanz, A. y H. Almus, 2003. Impact of hydrodynamics on development and morphology of intertidal seagrasses in the Wadden Sea. Marine Ecology Progress Series, 261: 123-134.
Selig, U., A. Eggert, D. Schories, M. Schubert, C. Blumel y H. Schubert, 2007. Ecological classification of macroalgae and angiosperm communities of inner coastal waters in the southern Baltic Sea. Ecological Indicators, 7: 665-678.
Sheppard, C. y R. Rioja-Nieto, 2005. Sea surface temperature 1871-2099 in 38 cells in the Caribbean region. Marine Environmental Research, 60: 389-396.
Short, F.T. y H.A. Neckles, 1998. The effects of global climate change on seagrasses. Aquatic Botany, 63(3-4): 169-196.
Short, F.T., E.W. Koch, J.C. Creed, K.M. Magalhaes, E. Fernández y J.L. Gaeckle, 2006. SeagrassNet monitoring across the Americas: case studies of seagrass decline. Marine Ecology, 27: 277-289.
Short, F.T., R.C. Davis, B.S. Koop, C.A. Short y D.M. Burdick, 2002. Site selection model for optimal transplantation of eelgrass Zostera marina in the northeastern US. Marine Ecology Progress Series, 227: 253-267.
Silván, F. y J.A. Campos, 2002. Estudio de la fl ora vascular amenazada de los estuarios de la Comunidad Autónoma del País Vasco. Dpto. de Medio Ambiente y Ordenación del Territorio, Gobierno Vasco, 100 pp.
Skov, H., E. Humphreys, S. Garthe, K. Geitner, D. Gremillet, K.C. Hamer, J. Hennicke, H. Parner y S. Wanless, 2008. Application of habitat suitability modelling to tracking data of marine animals as a means of analyzing their feeding. habitats. Ecological Modelling, 212(3-4): 504-512.
Solaun, O., M.J. Belzunce, J. Franco, V. Valencia y A. Borja, 2009a. Estudio de la contaminación en los sedimentos de los estuarios del País Vasco (1998-2001). Revista de Investigación Marina. 10: 47 pp.
Solaun, O., A. Borja, A. del Campo, J. Franco, J. Mader y M. González, 2009b. Estudio de impacto ambiental del dragado de la zona contigua a Astilleros de Murueta. Informe inédito elaborado por AZTI-Tecnalia para Astilleros de Murueta. 100 pp.
Tanner, C., S. Hunter, J. Reel, T. Parham, M. Naylor, L. Karrh, K. Busch, R.R. Golden, M. Lewandowski, N. Rybicki y E. Schenk, 2010. Evaluating a large-scale eelgrass restoration project in the Chesapeake Bay. Restoration Ecology, 18(4): 538-548.
Thom, R.M, 1990. A review of eelgrass (Zostera marina L.) transplanting projects in the Pacific Northwest. The Northwest Environmental Journal, 6: 121-137.
Tullrot, 2009. Background document for Zostera beds, seagrass beds. Biodiversity Series. Publication Number 426/2009. OSPAR Commission, London, United Kingdom, 37 pp.
van der Heide, T., E.T.H.M. Peeters, D.C.R. Hermus, M.M. van Katwijk, J.G.M. Roelofs, A.J.P. Smolders, 2009. Predicting habitat suitability in temperate seagrass ecosystems. Limnology and Oceanography, 54: 2018-2024.
van Katwijk, M.M. y D.C.R. Hermus, 2000. Effects of water dynamics on Zostera marina: transplantation experiments in the intertidal Dutch Wadden Sea. Marine Ecology Progress Series, 208: 107-118.
van Katwijk, M.M., D.C.R. Hermus, D.J. de Jong, R.M. Asmus y V.N. de Jonge, 2000. Habitat suitability of the Wadden Sea for restoration of Zostera marina beds. Helgoland Marine Research, 54: 117-128.
van Katwijk, M.M., L.J.M. Wijgergangs, 2004. Effects of locally varying exposure, sediment type and low-tide water cover on Zostera marina recruitment from seed. Aquatic Botany, 80: 1-12.
van Katwijk, M.M., A.R. Bos, V.N. de Jonge, L.S.A.M. Hanssen, D.C.R. Hermus, D.J. de Jong, 2009. Guidelines for seagrass restoration: importance of habitat selection and donor population, spreading of risks, and ecosystem engineering effects. Marine Pollution Bulletin, 58: 179-188.
Vermaat, J.E. y F.C.A. Verhagen, 1996. Seasonal variation in the intertidal Zostera noltii Hornem.: coupling demographic and physiological patterns. Aquatic Botany, 52: 259-281.
Vermaat, J.E., M.J.M. Hootsmans y P.H. Nienhuis, 1987. Seasonal dynamics and leaf growth of Zostera noltii Hornem., a perennial intertidal seagrass. Aquatic Botany, 28: 287-299.
Vermaat, J.E., N.S.R. Agawin, M.D. Fortes y J.S. Uri, 1997. The capacity of seagrasses to survive increased turbidity and siltation: the significance of growth form an light use. Ambio, 25(2): 499-504.
Viña, A., S. Bearer, H.M. Zhang, Z.Y. Ouyang y J.G. Liu, 2008. Evaluating MODIS data for mapping wildlife habitat distribution. Remote Sensing of Environment, 112(5): 2160-2169.
Waycott, M., B.J. Longstaff y J. Mellors, 2005. Seagrass population dynamics and water quality in the Great Barrier Reef region: a review and future research directions. Marine Ecology Progress Series, 51: 343-350.
West, R.J., N.E. Jacobs y D.E. Roberts, 1990. Experimental transplanting of seagrasses in Botany Bay, Australia. Marine Pollution Bulletin, 21(4): 197-203.
Wright, D.J., E.R. Lundblad, R.W. Rinehart, E.M. Larkin, J. Murphy y L. Cary-Kothera, 2005. Benthic Terrain Modeller Toolbar. Oregon State University Davey Jones Locker Seafl oor Mapping/Marine GIS Lab. http://dusk.geo.orst.edu/esri04/p1433_ron.html
Zipperle, A., J.A. Coyer, K. Reise, E. Gitz, W.T. Stam y J.N. Olsen, 2009. Clonal architecture in an intertidal bed of the dwarf eelgrass Zostera noltii in the Northern Wadden Sea: persistence through extreme physical perturbation and the importance of a seed bank. Marine Biology, 156: 2139-2148.
Como consecuencia de los trabajos realizados en relación a Zostera noltii en el País
Vasco se han elaborado y efectuado las siguientes comunicaciones y documentos:
�COMUNICACIONES A CONGRESOS COMUNICACIONES A CONGRESOS COMUNICACIONES A CONGRESOS COMUNICACIONES A CONGRESOS O REUNIONES O REUNIONES O REUNIONES O REUNIONES ((((5555))))
A.A.A.A.----PANELES o POSTERS (2):
1111....----"Selección de un estuario adecuado para recibir trasplantes de Zostera
noltii: ejemplo del País Vasco". Garmendia, J.M., Borja, Á., Rodríguez, J.G., Franco,
J. "XVI Simposio Ibérico de Estudio de Biología Marina""XVI Simposio Ibérico de Estudio de Biología Marina""XVI Simposio Ibérico de Estudio de Biología Marina""XVI Simposio Ibérico de Estudio de Biología Marina" (SIEBM), Alicante, 6-10
septiembre 2010.
2222....----"Restauración de praderas marinas: experimentos de trasplante de Zostera
noltii en estuarios del País Vasco". Garmendia, J.M., Borja, Á., Valle, M., Menchaca,
I., Franco, J. "XVI Simposio Ibérico de Estudio de Biología Marina""XVI Simposio Ibérico de Estudio de Biología Marina""XVI Simposio Ibérico de Estudio de Biología Marina""XVI Simposio Ibérico de Estudio de Biología Marina" (SIEBM),
Alicante, 6-10 septiembre 2010.
B.B.B.B.----ORALES (3):
1.1.1.1.---- "Sea level rise in the Basque coast from 1950 to 2010: observed trends,
future projections and impacts". Chust, G., Caballero, A., Marcos, M., Liria, P.,
Borja, Á., Galparsoro, I., Valle, M., Garmendia, J.M., Franco, J. "12 International "12 International "12 International "12 International
Symposium on Oceanography of the Bay of Biscay (ISOBAY)"Symposium on Oceanography of the Bay of Biscay (ISOBAY)"Symposium on Oceanography of the Bay of Biscay (ISOBAY)"Symposium on Oceanography of the Bay of Biscay (ISOBAY)", Brest (Francia), 3-6
mayo 2010.
2.2.2.2.---- "Predicting suitable habitat for Zostera noltii in the Oka estuary (Basque
Country) and its modification under mean sea level rise scenario". Valle, M., Borja,
Á., Galparsoro, I., Garmendia, J.M., Chust, G. "12 Internati"12 Internati"12 Internati"12 International Symposium on onal Symposium on onal Symposium on onal Symposium on
Oceanography of the Bay of Biscay (ISOBAY)"Oceanography of the Bay of Biscay (ISOBAY)"Oceanography of the Bay of Biscay (ISOBAY)"Oceanography of the Bay of Biscay (ISOBAY)", Brest (Francia), 3-6 mayo 2010.
3.3.3.3.---- "Zostera noltii restoration essays in the estuaries of the Basque Country".
Garmendia, J.M., Valle, M., Borja, Á. "I European Seagrass Restoration Workshop""I European Seagrass Restoration Workshop""I European Seagrass Restoration Workshop""I European Seagrass Restoration Workshop",
Portinho da Arrábida (Portugal), 30 septiembre-1 octubre 2010.
10. ANEXO I.I.I.I.---- Seguimiento mensual de la temperatura en el sedimento (Temp.) y potencial redox (Eh) Seguimiento mensual de la temperatura en el sedimento (Temp.) y potencial redox (Eh) Seguimiento mensual de la temperatura en el sedimento (Temp.) y potencial redox (Eh) Seguimiento mensual de la temperatura en el sedimento (Temp.) y potencial redox (Eh) medidos a unos 4 cm de profundidad y los parámetros de crecimiento de los parches de medidos a unos 4 cm de profundidad y los parámetros de crecimiento de los parches de medidos a unos 4 cm de profundidad y los parámetros de crecimiento de los parches de medidos a unos 4 cm de profundidad y los parámetros de crecimiento de los parches de Zostera noltiiZostera noltiiZostera noltiiZostera noltii trasplantados a sedimentos arenosos (ZT trasplantados a sedimentos arenosos (ZT trasplantados a sedimentos arenosos (ZT trasplantados a sedimentos arenosos (ZT----B1) B1) B1) B1) y fangosos (ZTy fangosos (ZTy fangosos (ZTy fangosos (ZT----B2) del estuario B2) del estuario B2) del estuario B2) del estuario del Butroe en marzo del 2009. d.t.=desviación típica; n.m.= no medido.del Butroe en marzo del 2009. d.t.=desviación típica; n.m.= no medido.del Butroe en marzo del 2009. d.t.=desviación típica; n.m.= no medido.del Butroe en marzo del 2009. d.t.=desviación típica; n.m.= no medido.