MEMORIA DESCRIPTIVA DEL PROCEDIMIENTO … Deslinde Arrubal.pdf · memoria descriptiva del...

MEMORIA DESCRIPTIVA DEL PROCEDIMIENTO ADMINISTRATIVO DE DELIMITACIÓN DEL DOMINIO PÚBLICO HIDRÁULICO DE LAS PARCELAS CATASTRALES Nº 995, 996, 1018 Y 1019 DEL POLÍGONO 1 DEL TÉRMINO MUNICIPAL DE ARRÚBAL (FINCA REGISTRAL Nº 2)

MEMORIA DESCRIPTIVA DEL PROCEDIMIENTO

ADMINISTRATIVO DE DELIMITACIÓN DEL DOMINIO

PÚBLICO HIDRÁULICO DE LAS PARCELAS

CATASTRALES Nº 995, 996, 1018 Y 1019 DEL

POLÍGONO 1 DEL TÉRMINO MUNICIPAL DE ARRÚBAL

(FINCA REGISTRAL Nº 2)

CSV: MA0031481164583709E10336281604591907Verificación en https://sede.miteco.gob.es

DOCUMENTO FIRMADO ELECTRÓNICAMENTELa Jefa del Área de Hidrología y Cauces - Moreno Santaengracia Maria Luisa. Sello de tiempo: 06/11/2020 7:55:02

https://sede.miteco.gob.es


ÍNDICE

MEMORIA

1. INTRODUCCIÓN

2. OBJETO DEL DESLINDE

3. CARACTERISTICAS DEL TRAMO OBJETO DE DESLINDE

4. CARACTERISTICAS DE LA PROPIEDAD DE LOS TERRENOS

5. TRABAJOS REALIZADOS

6. FASES DEL PROCEDIMIENTO ADMINISTRATIVO EFECTUADAS HASTA LA FECHA

7. PROPUESTA DE DESLINDE

ANEXOS

ANEXO 1: RELACION DE PARCELAS AFECTADAS Y TITULARES

ANEXO 2: ESTUDIO HIDROLÓGICO

ANEXO 3: ESTUDIO HIDRÁULICO

ANEXO 4: ESTUDIO HISTORICO, GEOMORFOLÓGICO Y ECOLÓGICO

ANEXO 5: PLANOS PROPUESTA DE LA LINEA QUE DEFINE EL D.P.H.

ANEXO 6: REPORTAJE FOTOGRÁFICO





MEMORIA




MEMORIA


INDICE

1. INTRODUCCIÓN ............................................................................. 1

2. OBJETO DEL DESLINDE ................................................................... 1

3. CARACTERÍSTICAS DEL TRAMO OBJETO DE DESLINDE................. 2

4. CARACTERÍSTICAS DE LA PROPIEDAD DE LOS TERRENOS ........... 3

5. TRABAJOS REALIZADOS ................................................................. 3

6. FASES DEL PROCEDIMIENTO ADMINISTRATIVO EFECTUADAS HASTA LA FECHA .................................................................................... 5

7. PROPUESTA DE DESLINDE ............................................................... 7




MEMORIA


1

1. INTRODUCCIÓN

Se redacta la presente memoria descriptiva a partir de la información disponible en

este Organismo de Cuenca, y de la aportada hasta el momento en el procedimiento

administrativo de delimitación del dominio público hidráulico de las parcelas catastrales nº

995, 996, 1018 y 1019 del polígono 1 del término municipal de Arrúbal (finca registral nº 2),

conforme al Art. 242.3 del Real Decreto 606/2003, de 23 de mayo, por el que se modifica el

Real Decreto 849/1986, de 11 de abril, por el que se aprueba el Reglamento del Dominio

Público Hidráulico.

2. OBJETO DEL DESLINDE

Según el artículo 2. b) del Texto Refundido de la Ley de Aguas, constituyen el dominio

público hidráulico del Estado, entre otros, los cauces de corrientes naturales, continuas o

discontinuas, definiendo en su artículo 4 el álveo o cauce natural de una corriente continua

o discontinua como el terreno cubierto por las aguas en las máximas crecidas ordinarias.

El artículo 4.1 del Reglamento del Dominio Público Hidráulico, aprobado por Real Decreto 849/1986, de 11 de abril (en lo sucesivo RDPH), tras reproducir esa misma definición legal, prevé que la determinación de ese terreno se realizará atendiendo a sus características geomorfológicas, ecológicas y teniendo en cuenta las informaciones hidrológicas, hidráulicas, fotográficas y cartográficas que existan, así como las referencias históricas disponibles. El apartado 2 del artículo 4 del RDPH, en su primer párrafo, define el caudal de la máxima crecida ordinaria para los tramos de cauce donde exista información hidrológica suficiente, como la media de los máximos caudales instantáneos anuales en su régimen natural, calculada a partir de las series de datos existentes y seleccionando un período que incluirá el máximo número de años posible y será superior a diez años consecutivos. Dicho periodo será representativo del comportamiento hidráulico de la corriente y en su definición se tendrá en cuenta las características geomorfológicas, ecológicas y referencias históricas disponibles. Para los tramos de cauce en los que no exista información hidrológica suficiente, el párrafo segundo del apartado 2 del artículo 4 antes señalado, indica que el caudal de la máxima crecida ordinaria se establecerá a partir de métodos hidrológicos e hidráulicos alternativos, y, en especial, a partir de la simulación hidrológica e hidráulica de la determinación del álveo o cauce natural y teniendo en cuenta el comportamiento hidráulico de la corriente, las características geomorfológicas, ecológicas y referencias históricas disponibles.

La función de policía que la Ley atribuye a la Administración hidráulica sobre los cauces del dominio público hidráulico y sus zonas adyacentes, servidumbre y policía, se ejerce a través de diversas potestades entre las que se encuentra la de deslinde. Esta función protectora de los cauces contribuye igualmente a evitar y disminuir los riesgos asociados a las avenidas.




MEMORIA


2

Con fecha 9 de septiembre de 2018, la Letrada Dª Isabel Martínez García, en su

condición de miembro del consejo de administración judicial, nombrado en el

procedimiento de Ejecución de Títulos Judiciales ETJ 1247/16 del Juzgado de Primera

instancia nº 3 de Logroño, solicitó a este Organismo el comienzo a las operaciones precisas

para el deslinde del dominio público hidráulico en relación a las parcelas catastrales 995,

996, 1018 y 1019 del polígono nº 1, (finca registral nº 2), todas ellas en el Término Municipal de

Arrúbal (La Rioja), a fin de delimitar los terrenos que son propiedad del Grupo Sindical de

Colonización nº 24 de Arrúbal y los que pertenecen al Dominio Público Hidráulico, todo ello

de conformidad con el Auto judicial de 18 de mayo de 2017 del Juzgado antes referido, que

expresamente establecía en su fundamento de derecho IV lo siguiente:

“Es obligación del órgano de administración llevar a cabo todos los actos tendentes a

proceder al deslinde de las parcelas catastrales indicadas respecto a la parte que

pertenece al Dominio Público Hidráulico, debiendo para ello realizar las gestiones que sean

oportunas ante la Confederación Hidrográfica del Ebro, organismo administrativo

correspondiente”.

3. CARACTERÍSTICAS DEL TRAMO OBJETO DE DESLINDE

Las parcelas catastrales cuyo deslinde se solicita son las nº 995, 996, 1018 y 1019 del

polígono 1 del término municipal de Arrúbal, (Finca registral nº 2):

POL. PARCELA CATASTRAL

SUPERFICIE TOTAL CATASTRAL (M²)

1 995 86.534

1 996 36.319

1 1018 5.066

1 1019 1.184




MEMORIA


3

El tramo del río Ebro en el que se ubican las fincas a deslindar viene definido por las siguientes coordenadas ETRS89_UTM_H30:

SECCIÓN MARGEN COORDENADAS ETRS89_UTM_H30 LONGITUD (M)

INICIO Derecha X 562 524,8

1.485,83 Y 4 699 069

FINAL Derecha X 526 875,5

Y 4 697 705

4. CARACTERÍSTICAS DE LA PROPIEDAD DE LOS TERRENOS

En el Anexo I se adjunta la relación de los titulares de los terrenos afectados, elaborada a

partir de la cartografía digitalizada existente en este Organismo y la documentación

remitida al inicio del expediente de deslinde por la Gerencia Regional del Catastro de La

Rioja. La citada relación ha sido completada con la información aportada por el

Ayuntamiento de Arrúbal, en la fase procedimental correspondiente.

5. TRABAJOS REALIZADOS

Los estudios y trabajos principales, que han servido de base para la propuesta de deslinde del Dominio Público Hidráulico que se presenta en esta memoria, son los siguientes:

0. -Recopilación de la información existente

1.- Levantamiento topográfico de la zona obtenida mediante tecnología LIDAR.

2.- Estudio de la hidrología del tramo donde se ubican las fincas a deslindar, con base

en la información pluviométrica y foronómica disponible, mediante modelos

matemáticos que permiten deducir el caudal teórico de la máxima crecida ordinaria.

3.- Estudio hidráulico que permite fijar el nivel alcanzado por las aguas en la máxima

crecida ordinaria.

4.- Estudio de la geomorfología del terreno.




MEMORIA


4

5.- Estudio de los condicionantes ecológicos del tramo de río donde se encuentran

ubicadas las fincas objeto de deslinde y su entorno.

6.- Antecedentes administrativos existentes en los archivos de la Comisaría de aguas.

También se han considerado como elementos coadyuvantes para establecer la propuesta de deslinde del DPH otros criterios como son la observación del terreno y de las condiciones topográficas de la zona.

En el Anexo I se presenta la relación de parcelas afectadas y titulares.

En el Anexo II, que se adjunta a la presente memoria, se incluye un resumen de los estudios hidrológicos existentes y los realizados, aplicados al tramo del río donde se ubican las fincas a deslindar objeto de este expediente.

En el Anexos III se incluye un resumen del estudio hidráulico realizado para el tramo de río donde se ubican las fincas a deslindar objeto de este expediente.

En el Anexo IV se presenta el estudio histórico, geomorfológico y ecológico del tramo donde se ubican las fincas a deslindar objeto de este expediente.

En el Anexo V se adjuntan planos con la línea de propuesta del DPH a escala 1/1000.

En el Anexo VI se incluye un reportaje fotográfico de la visita de campo realizada.




MEMORIA


5

6. FASES DEL PROCEDIMIENTO ADMINISTRATIVO EFECTUADAS HASTA LA FECHA

6.1 TRABAJOS PREVIOS AL ACUERDO DE INCOACIÓN, RELATIVOS A LA INVESTIGACIÓN DE LA TITULARIDAD DE LOS TERRENOS AFECTADOS POR EL DESLINDE

Con fecha 9 de septiembre de 2018, la Letrada Dª Isabel Martínez García, en su condición

de miembro del consejo de administración judicial, nombrado en el procedimiento de

Ejecución de Títulos Judiciales ETJ 1247/16 del Juzgado de Primera instancia nº 3 de Logroño,

solicitó a este Organismo el comienzo de las operaciones precisas para el deslinde del

dominio público hidráulico, en relación a las parcelas catastrales 995, 996, 1018 y 1019 del

polígono nº 1, (finca registral nº 2), todas ellas en el Término Municipal de Arrúbal (La Rioja), a

fin de delimitar los terrenos que son propiedad del Grupo Sindical de Colonización nº 24 de

Arrúbal y los que pertenecen al Dominio Público Hidráulico, todo ello de conformidad con el

Auto judicial de 18 de mayo de 2017 del Juzgado antes referido, que expresamente

establecía en su fundamento de derecho IV lo siguiente:

“Es obligación del órgano de administración llevar a cabo todos los actos tendentes a

proceder al deslinde de las parcelas catastrales indicadas respecto a la parte que

pertenece al Dominio Público Hidráulico, debiendo para ello realizar las gestiones que

sean oportunas ante la Confederación Hidrográfica del Ebro, organismo administrativo

correspondiente”.

Con el fin de estudiar la titularidad de las fincas que, previsiblemente se encuentran

afectadas por el procedimiento, con fecha 3 de julio de 2020 se solicitó a la Gerencia

Regional de Catastro de La Rioja que, facilitase los planos, datos de identidad y domicilio de

los titulares catastrales de las fincas cuyo deslinde es objeto del presente procedimiento.

Con fecha de registro de entrada en este Organismo de 24 de julio de 2020, la Gerencia

Regional de Catastro de La Rioja remitió certificaciones catastrales, con la identificación y el

domicilio fiscal de los titulares catastrales; así como, un fichero en formato gml con la

representación gráfica georreferenciada de las parcelas.

De la misma manera, se remitió al Ayuntamiento de Arrúbal con fecha 3 de julio de 2020

información sobre las fincas a deslindar, con el objeto de que manifestase su conformidad o

adjuntara la información y observaciones que considerara oportunas. Con fecha 7 de

agosto de 2020, el Ayuntamiento de Arrúbal remitió certificaciones catastrales descriptivas y

gráficas de las fincas solicitadas, incluyéndose: titularidad, cultivo y relación de parcelas

colindantes, con sus respectivos titulares.




MEMORIA


6

6.2 ACUERDO DE INCOACIÓN

Para dar trámite a la solicitud del Juzgado de Primera instancia nº 3 de Logroño, con fecha 7

de julio de 2020 se acordó por parte de la Confederación Hidrográfica del Ebro, la incoación

del procedimiento administrativo de delimitación de dominio público hidráulico de las fincas

catastrales nº 995, 996, 1018 y 1019 del polígono 1 del término municipal de Arrúbal (finca

registral nº 2) ubicadas en un tramo del río Ebro de longitud 1.485,83 metros.

El 8 de julio de 2020 se comunica a las Áreas de Calidad de Aguas y de Gestión del Dominio

Público Hidráulico, la suspensión cautelar del otorgamiento de concesiones y autorizaciones

que pudiesen afectar al dominio público hidráulico o dificultasen los trabajos que deben

realizarse para su delimitación.

El Acuerdo de incoación de 7 de julio de 2020, se publicó en el Boletín Oficial de La Rioja con

fecha 5 de agosto de 2020 (BOR Núm. 97), y en el Diario de La Rioja (periódico de máxima

difusión en la zona) el 15 de julio de 2020. A su vez, fue comunicado a la Gerencia Territorial

de Catastro de La Rioja y al Ayuntamiento de Arrúbal que, con fecha de registro de entrada

en este Organismo de 11 de agosto de 2020, remitió la correspondiente diligencia de

exposición del acuerdo de incoación en el tablón de edictos de esa Corporación.

El Acuerdo de Incoación fue notificado también, de forma individualizada, a los titulares de

los terrenos y a los colindantes, previsiblemente afectados, mediante carta certificada con

acuse de recibo. Devuelta la notificación correspondiente al Grupo Sindical de Colonización

Nº 24 de Arrúbal, el interesado comunica la nueva dirección a efecto de notificaciones, así

como su dirección de correo electrónico, que han sido incorporadas al expediente.

Con fecha 23 de septiembre de 2020, de conformidad con el art. 242.3.b, se remitió la

relación de titulares previsiblemente afectados por el procedimiento de deslinde, al Registro

de la Propiedad nº 2 de Logroño, con el objeto de que manifestase su conformidad con la

misma o formulase las observaciones que estime pertinentes. Con fecha de registro de

entrada en este Organismo de 28 de octubre de 2020, el Registro de la Propiedad remitió

certificados señalando que las fincas sobre las que se solicitaba la información no aparecen

inscritas en dicho registro.

6.3 DOCUMENTACIÓN SOLICITADA Y APORTADA DURANTE LAS PRIMERAS FASES DEL PROCEMIENTO

Una vez notificado de forma individualizada el Acuerdo de Incoación a todos los

propietarios de las fincas afectadas por el procedimiento de deslinde, no se ha recibido

documentación hasta la fecha en este Organismo, a excepción de las certificaciones

catastrales remitidas por el Ayuntamiento de Arrúbal con fecha 7 de agosto de 2020, de la

que se ha hecho referencia en el apartado 6.1, y la comunicación de nueva dirección del




MEMORIA


7

Grupo Sindical de Colonización Nº 24 de Arrúbal, mencionado en el apartado 6.2 de la

presente memoria.

7. PROPUESTA DE DESLINDE

A partir de los resultados obtenidos en el estudio hidrológico, estudio hidráulico para la

Máxima Crecida Ordinaria, el estudio histórico, geomorfológico y ecológico, se ha definido

la propuesta de la línea que define el Dominio Público Hidráulico. También ha sido

considerada la información proporcionada por la topografía del terreno y la visita de

campo realizada, de conformidad con el artículo 4.1 del Reglamento del Dominio Público

Hidráulico, aprobado por Real Decreto 849/1986, de 11 de abril.

Se incluyen en el Anexo V los planos con la línea del Dominio Público Hidráulico definida por

los siguientes puntos referidos a coordenadas ETRS89/H30:

PARCELA DENOMINACIÓN ETRS89/H30

UTM_X UTM_Y 1018 H-1018\01 562 516,60 4 698 832,90 1018 H-1018\02 562 512,30 4 698 801,50 1018 H-1018\03 562 506,60 4 698 761,80 1018 H-1018\04 562 503,10 4 698 746,80 1018 H-1018\05 562 488,20 4 698 705,10 1018 H-1018\06 562 470,80 4 698 661,90 1019 -- -- -- 0995 H-0995\01 562 491,70 4 698 075,00 0995 H-0995\02 562 487,40 4 698 027,70 0995 H-0995\03 562 482,60 4 697 970,60 0995 H-0995\04 562 483,10 4 697 947,80 0995 H-0995\05 562 489,00 4 697 929,30 0995 H-0995\06 562 496,40 4 697 912,90 0995 H-0995\07 562 518,60 4 697 886,40 0995 H-0995\08 562 535,50 4 697 851,00 0995 H-0995\09 562 544,50 4 697 817,10 0995 H-0995\10 562 541,90 4 697 794,90 0995 H-0995\11 562 532,30 4 697 777,40 0995 H-0995\12 562 529,20 4 697 756,30 0995 H-0995\13 562 533,70 4 697 732,50 0996 H-0996\01 562 663,20 4 697 590,30 0996 H-0996\02 562 689,10 4 697 564,30




MEMORIA


8

De los estudios efectuados se desprende que, las fincas sobre las que se ha solicitado el

deslinde, se encuentran afectadas por el Dominio Público Hidráulico potencial de la

siguiente manera:

POL. PARCELA DPH / COLINDANCIA AFECCIÓN

1 995 DENTRO DEL DPH POTENCIAL PARCIAL



1 1019 FUERA SERVIDUMBRE

Resultado de la propuesta de deslinde, las siguientes fincas son colindantes con el Dominio

Público Hidráulico potencial, encontrándose fuera del mismo y sujetas a la zona de

servidumbre, de conformidad con el Art. 6 del RDPH:

POL. PARCELA DPH / COLINDANCIA AFECCIÓN



Zaragoza, a 5 de noviembre de 2020

LA JEFA DEL ÁREA DE HIDROLOGÍA Y CAUCES

D. María Luisa Moreno Santaengracia




ANEXO Nº 1. RELACIÓN DE PARCELAS AFECTADAS Y TITULARES


ANEXO Nº1.

RELACIÓN DE PARCELAS AFECTADAS Y TITULARES




ANEXO Nº 1. RELACIÓN DE PARCELAS AFECTADAS Y TITULARES


1. LISTADO DE PARCELAS

DPH / COLINDANCIA AFECCIÓN MARGEN NATURALEZA T.M. POLÍGONO PARCELA

CATASTRAL PARAJE NOMBRE

DENTRO DEL DPH POTENCIAL

PARCIAL DERECHA RÚSTICA ARRÚBAL 1 995 EL SOTO AYUNTAMIENTO DE ARRUBAL

DENTRO DEL DPH POTENCIAL PARCIAL DERECHA RÚSTICA ARRÚBAL 1 996 CABEZUELOS AYUNTAMIENTO DE ARRUBAL

DENTRO DEL DPH POTENCIAL PARCIAL DERECHA RÚSTICA ARRÚBAL 1 1018 EL SOTO AYUNTAMIENTO DE ARRUBAL

COLINDANTE CON EL DPH POTENCIAL SERVIDUMBRE DERECHA RÚSTICA ARRÚBAL 1 1019 EL SOTO AYUNTAMIENTO DE ARRUBAL

COLINDANTE CON EL DPH POTENCIAL SERVIDUMBRE DERECHA RÚSTICA ARRÚBAL 1 30371 PONTARRO CONSORCIO DE AGUAS Y RESIDUOS DE LA RIOJA

COLINDANTE CON EL DPH POTENCIAL SERVIDUMBRE DERECHA RÚSTICA ARRÚBAL 1 10371 PONTARRO GRUPO SINDICAL NUMERO 24




ANEXO Nº 2. ESTUDIO HIDROLÓGICO


ANEXO Nº2.

ESTUDIO HIDROLÓGICO






ÍNDICE

1. ESTUDIO HIDROLÓGICO. ........................................................................... 1

1.1. INTRODUCCIÓN. .................................................................................................. 1

1.2. DETERMINACIÓN DE LA MÁXIMA CRECIDA ORDINARIA ................................ 1

1.2.1. ESTUDIO DE LOS DATOS FORONÓMICOS. ................................................................. 1 1.2.1.1. Metodología. .......................................................................................................... 1 1.2.1.2. Selección de estación de aforos. ...................................................................... 2 1.2.1.3. Aplicación de la metodología a la estación seleccionada. ....................... 5

1.2.2. ESTUDIO MEDIANTE CAUMAX (CEDEX) ...................................................................... 8 1.2.2.1. Introducción. .......................................................................................................... 8 1.2.2.2. Metodología. .......................................................................................................... 8 1.2.2.3. Estimación de la Máxima Crecida Ordinaria con CAUMAX. ..................... 24 1.2.2.4. Resultados de CAUMAX. .................................................................................... 25

1.2.3. CAUDAL ADOPTADO PARA LA MCO. CONCLUSIONES. ....................................... 25






1

1. ESTUDIO HIDROLÓGICO.

1.1. INTRODUCCIÓN.

Estos estudios tienen por objeto la estimación, valoración o cálculo de la máxima crecida ordinaria en cada tramo fluvial, que es el volumen de agua que al circular ocupa el terreno que el texto refundido de la ley de aguas define como cauce de río y por tanto dominio público hidráulico.

En consecuencia, el objeto de este estudio hidrológico es, en base al análisis de la información existente sobre el comportamiento del río Ebro en el tramo que comprenden la zona objeto de estudio y su entorno, conocer la respuesta hidrológica de la cuenca, y en concreto el caudal asociado a su máxima crecida ordinaria.

Se han planteado dos metodologías para su estimación:

• Estudios mediante datos foronómicos.

• Estudios mediante el modelo CAUMAX.

1.2. DETERMINACIÓN DE LA MÁXIMA CRECIDA ORDINARIA

1.2.1. ESTUDIO DE LOS DATOS FORONÓMICOS.

1.2.1.1. Metodología.

Este estudio se ha basado en el análisis de los datos foronómicos que dispone la Comisaría de Aguas de la Confederación Hidrográfica del Ebro para las distintas estaciones de aforo emplazadas dentro del área de estudio, con el fin de establecer el caudal de la máxima crecida ordinaria de acuerdo con su descripción.

Estos datos foronómicos consisten en valores de niveles de agua en dichos puntos que pueden ser transformados en caudales a partir de la curva de gastos obtenida a partir de aforos directos en campo y de formulaciones teóricas.

El proceso realizado para determinar los caudales de la máxima crecida ordinaria ha sido el siguiente:

1. Recopilación de las series históricas correspondientes a las estaciones de aforo incluidas dentro del estudio.

2. Revisión de las series de niveles (limnogramas) y las correspondientes curvas de gasto. Revisión de los hidrogramas.

3. Estudio e identificación de las obras hidráulicas (presas, azudes, canales de derivación, etc.) que puedan influir en los registros de las estaciones de aforo. Estas infraestructuras modifican el régimen natural de flujo.






2

4. Identificación de las fechas en las cuales puede haberse modificado el sistema de gestión hidráulica del río (por obras o cambios en los criterios de gestión). Esta circunstancia puede hacer que las series no sean homogéneas.

5. Homogeneizar las series. Desde un punto de vista estadístico, deben pertenecer a la misma población; en caso contrario no sería correcto el análisis estadístico.

6. Obtención de la serie de caudales máximos anuales registrados en las estaciones de aforo con los criterios de gestión de las infraestructuras hidráulicas actuales.

7. Cálculo de los caudales, reales y restituidos a régimen natural, de máxima crecida ordinaria y de avenida.

Se ha obtenido el caudal de la máxima crecida ordinaria a partir de su descripción legal. En el apartado 1 del artículo 4 del Reglamento del Dominio Público Hidráulico se define el cauce natural de una corriente continua o discontinua como “el terreno cubierto por las aguas en las máximas crecidas ordinarias”. En el apartado 2 de ese mismo artículo explicita que en los cauces, o tramos de cauce, con información hidrológica suficiente “se considerará caudal de la máxima crecida ordinaria la media de los máximos caudales instantáneos anuales en su régimen natural, calculada a partir de las series de datos existentes y seleccionando un período que incluirá el máximo número de años posible y será superior a diez años consecutivos” y además que “dicho periodo será representativo del comportamiento hidráulico de la corriente y en su definición se tendrá en cuenta las características geomorfológicas, ecológicas y referencias históricas disponibles”. Igualmente, también refiere que, en los cauces, o tramos de cauce, en los que no se dispone de dicha información, “el caudal de la máxima crecida ordinaria se establecerá a partir de métodos hidrológicos e hidráulicos alternativos, y, en especial, a partir de la simulación hidrológica e hidráulica de la determinación del álveo o cauce natural”.

1.2.1.2. Selección de estación de aforos.

Las estaciones de aforos más cercanas son la 9280, “Río Ebro en Logroño” y la 9120 “Río Ebro en Mendavia”. En la figura 1 se presenta una ficha resumen con las características de la estación 9120 (A120, en SAIH Ebro), localizada en las inmediaciones de la zona de estudio, aguas abajo de la misma, y en la que se puede constatar que la serie hidrológica tiene calidad y extensión suficiente para cumplir con los requisitos exigidos por el reglamente para su empleo en el cálculo de la Máxima Crecida Ordinaria.

En la figura 2 se reproduce la ficha de la estación de aforos 9280 (A-280, en SAIH Ebro), situada aguas arriba de la zona de estudio. Se puede observar que, aunque la longitud de la serie, de once años, tiene la duración mínima requerida por la normativa vigente para su empleo en el cálculo del caudal de la máxima crecida ordinaria, la calidad de sus datos, con cinco años incompletos, no permite su empleo con las garantías mínimas por falta de información.






3

FIGURA 1: Ficha característica Estación de Aforos 9120 Ebro en Mendavia

FIGURA 2: Ficha característica Estación de Aforos 9280 Ebro en Logroño.

La siguiente estación de aforos en servicio localizada aguas arriba de la de Logroño es la denominada 9001, “Río Ebro en Miranda de Ebro”. Esta estación presenta una serie de larga duración con calidad sobrada para su empleo en el cálculo de la máxima crecida ordinaria (Figura 3). Dada su lejanía a la zona de estudio, existen importantes cursos que desembocan en el Ebro en el interfluvio desarrollado entre ella y la zona de estudio (Figura 4).






4

FIGURA 3: Ficha característica Estación de Aforos 9001 Ebro en Miranda de Ebro

FIGURA 4: Localización de las estaciones de aforo sobre el Ebro en servicio en relación con la zona de

estudio

Dada la mayor proximidad de la estación 9120 a la zona de estudio y para evitar la minoración de caudales que supondría utilizar la estación 9001 (que exluiría los caudales aportados por la intercuenca y los ríos afluentes Bayas y Zadorra por margen izquierda y Tirón, Najerilla, Iregua y Leza por margen derecha) se ha seleccionado como estación de aforos de referencia la estación 9120 “Río Ebro en Mendavia”.






5

1.2.1.3. Aplicación de la metodología a la estación seleccionada.

La serie de datos indicados es una serie en régimen regulado, estando principalmente afectados sus valores por los efectos de la gestión del embalse del Ebro y, en menor medida, la de los embalses de Ullibarri-Gamboa y Urrunaga.

El efecto de estos embalses sobre los caudales registrados implica una pequeña laminación de las avenidas y por lo tanto disminución de los caudales circulantes por la estación de aforos.

Esta disminución de los valores máximos registrados es inviable obtenerla de modo preciso, ya que la gestión de los embalses se ha ido optimizando con el paso de los años, principalmente debido a la información hidrológica disponible en tiempo real (proporcionada por el SAIH) y a que cada vez se ha dispuesto de un Sistema de Ayuda a la Decisión (SAD) más eficiente y completo compuesto por modelos de previsión meteorológica, hidrológica e hidráulica.

En la siguiente tabla quedan reflejados los principales datos, correspondiente a los registros de la estación 9120, una vez han sido analizados, validados y calculados los caudales (aplicando en cada caso la curva de gasto correspondiente a su periodo temporal). Es fundamental señalar que estos datos se han obtenidos del Anuario de Aforos de la Red Oficial de Estaciones de Aforo, publicados por el CEDEX, y por lo tanto son los datos oficiales, siendo los últimos datos disponibles los correspondientes al año hidrológico 2016/17.

Año hidrológico

Qci Caudal máximo instantáneo (m3/s)

Media móvil para 11 años consecutivos (m3/s)

Media para la serie más larga de años consecutivos con datos completos (m3/s)(*)

1948/49

1949/50

1950/51 703

912,7

1951/52 855 1952/53 775 1953/54 1282 1954/55 820 1955/56 755 1956/57 467 1957/58 775 1958/59 792 1959/60 1476 1960/61 1307 909,7 1961/62 1270 961,3 1962/63 1183 991,1 1963/64 672 981,7 1964/65 539 914,2 1965/66 994 930,0 1966/67 970,8 949,6 1967/68 1132 1010,1






6

Año hidrológico




1968/69 573 991,7 1969/70

1970/71

1971/72 989,8 1972/73 970,8 1973/74 986 1974/75 1092,4 1975/76 784,2 1976/77 1476 1977/78 1252 1978/79 1247,2 1979/80

1980/81 1903,2 1981/82 478,55 1982/83 1406,4 1983/84 753,05 1984/85 623,2 1985/86 614,4 1986/87 509 1987/88 1239,2 1988/89 298,8 1989/90

1990/91 898,4 1991/92 887,25 1992/93 825,4 1993/94 1659,5 1994/95 1159,2 1995/96 494,2 1996/97 808,6 1997/98 744,2 1998/99 449,1 1999/00 631,08 2000/01 1213 888,2 2001/02 469,5 849,2 2002/03

2003/04

2004/05 423,932 2005/06 671,502 2006/07 944,981 2007/08 1212,702 2008/09 939,22 2009/10 968,4






7

Año hidrológico




2010/11

2011/12 485,8 2012/13 963 2013/14 727 2014/15 1283,4 2015/16 1125,9 2016/17 861,3

(*) La serie debe tener una longitud temporal no inferior a 11 años

Debido a la importante incertidumbre en la restitución a régimen natural de los datos (exigida por la definición de la MCO), se han aplicado los análisis estadísticos directamente a los datos registrados en la estación de aforos, viendo como han ido evolucionando temporalmente los valores, siendo conscientes de que se va a obtener como resultado un umbral inferior para el caudal correspondiente a la MCO.

Hay que tener en cuenta que el embalse del Ebro, se pone en explotación en 1946, por lo que, a pesar de que su efecto es muy reducido, toda la serie se ve afectada por su laminación.

Como se ha indicado anteriormente, el apartado 2 del artículo 4 del Reglamento del Dominio Público Hidráulico indica que “se considerará caudal de la máxima crecida ordinaria la media de los máximos caudales instantáneos anuales en su régimen natural, calculada a partir de las series de datos existentes y seleccionando un período que incluirá el máximo número de años posible y será superior a diez años consecutivos”

En base a esta disposición y toda la información disponible, se toma como periodo representativo el comprendido entre 1950/51 y 1968/69, al que le corresponde un caudal de 912,7 m3/s.

Año hidrológico


Media para la serie más larga de años consecutivos con datos

completos (m3/s) 1950/51 703

912,7

1951/52 855 1952/53 775 1953/54 1282 1954/55 820 1955/56 755 1956/57 467 1957/58 775 1958/59 792 1959/60 1476 1960/61 1307 1961/62 1270






8

Año hidrológico


Media para la serie más larga de años consecutivos con datos

completos (m3/s) 1962/63 1183 1963/64 672 1964/65 539 1965/66 994 1966/67 970,8 1967/68 1132 1968/69 573

Estos 19 años hidrológicos consecutivos, en régimen natural, son representativos del comportamiento hidráulico de la corriente y, además, al estar lo más alejado temporalmente de la actualidad, previsiblemente el efecto laminador del embalse fuera todavía menor, así como el nivel de antropización del cauce.

Ahora bien, teniendo en cuenta la longitud de la serie y que ésta presenta un total de 60 años hidrológicos con información completa, parece lógico que la media de los caudales máximos instantáneos registrados sea tan representativa o más que la serie de 19 años reseñada para el cálculo de la máxima crecida ordinaria. Considerando todos los años con dato disponible, el caudal máximo instantáneo medio resultante asciende a 913 m3/s, valor idéntico en la práctica al obtenido para el periodo representativo seleccionado.

1.2.2. ESTUDIO MEDIANTE CAUMAX (CEDEX)

1.2.2.1. Introducción.

Dentro del ámbito del Sistema Nacional de Cartografía de Zonas Inundables (SNCZI), el entonces Ministerio de Agricultura, Alimentación y Medio Ambiente, encargó al CEDEX mediante un Convenio de colaboración la elaboración de un mapa de los caudales máximos asociados a distintas probabilidades de recurrencia en la red de ríos que gestiona a través de las distintas Confederaciones Hidrográficas.

Este trabajo se presenta en la aplicación CAUMAX, integrada en un sistema de información geográfica, en la que es posible consultar los caudales máximos instantáneos en régimen natural asociados a distintos periodos de retorno para los cauces con una cuenca superior a 50 km2 y calcular estos caudales mediante el método racional modificado para cauces con cuencas inferiores a 50 km2.

Esta herramienta presenta un valor especial en cuencas de gran tamaño, afectadas por condicionantes regionales, ya que tiene en cuenta la sincronización de las avenidas en las zonas de confluencia como la que nos ocupa. No obstante, los resultados de este trabajo son una referencia inicial de los caudales de avenidas existentes, que en cada caso podrán ser completados a partir de estudios específicos y más detallados.

1.2.2.2. Metodología.

La metodología aplicada en el estudio ha partido de la selección de estaciones de aforos y embalses con suficiente información sobre caudales máximos como para



http://hercules.cedex.es/caumax/caumax_v21.rar




9

realizar una estimación suficientemente fiable de los valores correspondientes a distintas probabilidades o periodos de retorno (cuantiles). A las series de datos seleccionadas se les han aplicado diversas técnicas estadísticas para verificar la calidad de la información, así como la coherencia en la información de las distintas estaciones entre sí y se les ha ajustado una función de distribución con la finalidad de determinar sus cuantiles.

En el análisis de las leyes de frecuencia se han incorporado las siguientes técnicas estadísticas con objeto de mejorar la precisión de la estimación de los cuantiles, fundamentalmente de los correspondientes a los más altos periodos de retorno (100 y 500 años) para cuya estimación es necesario extrapolar la ley de frecuencia fuera del rango de probabilidad cubierto por los datos presentes en la serie:

• Estudio de la combinación función de distribución–procedimiento de ajuste con un mejor comportamiento para el tipo de datos de la cuenca.

• Empleo de información regional en la determinación del valor del coeficiente de sesgo. Este estadístico tiene una gran influencia en la extrapolación de la función de distribución a altos periodos de retorno y, sin embargo, su estimación a partir de muestras de pequeña longitud presenta una gran imprecisión.

• Incorporación de información histórica al análisis de las leyes de frecuencia en aquellos lugares donde ese tipo de información está disponible. La información histórica correspondiente a las mayores avenidas ocurridas en la zona, es decir, a eventos de pequeña probabilidad de ocurrencia, proporciona información directa sobre el valor de los caudales en el rango de probabilidad no cubierto por los datos de aforo sistemático, mejorando de forma importante la estimación de los cuantiles de alto periodo de retorno.

El análisis estadístico de los datos foronómicos permite conocer los cuantiles en aquellos puntos en los que existen medidas de caudal, pero no puede aplicarse en aquellos puntos donde no están disponibles dichas medidas. Puesto que el objetivo del trabajo es elaborar un mapa que dé información sobre los caudales máximos de forma casi continua a lo largo de la red fluvial, es necesario estimar también los cuantiles en aquellos puntos no aforados. Dicha estimación debe realizarse forzosamente mediante extrapolación de los resultados obtenidos previamente para los puntos aforados.

La extrapolación de resultados puede realizarse básicamente mediante dos procedimientos:

• Modelos hidrometeorológicos calibrados de tal forma que reproduzcan adecuadamente los resultados, desde un punto de vista estadístico, en los puntos aforados, es decir, que para las hipótesis de cálculo adoptadas den como resultado los cuantiles previamente calculados mediante el análisis estadístico, asumiendo que el modelo así calibrado da resultados correctos en los puntos no aforados. La obtención de los cuantiles de caudal es fundamentalmente un problema estadístico y no un problema de modelación hidrometeorológica. La interpretación estadística de los resultados del modelo sólo podría llevarse a cabo si se conociese la probabilidad conjunta de las distintas variables que intervienen en la determinación del caudal (cantidad total de precipitación






10

durante la tormenta, distribución espacial y temporal de la precipitación, humedad antecedente del suelo en la cuenca, etc.). En la práctica dichas probabilidades no se conocen, por lo que hay que recurrir a adoptar determinadas hipótesis de cálculo (precipitación uniforme sobre la cuenca, distribución temporal según hietogramas sintéticos, etc.) y a calibrar determinados parámetros de la modelación (habitualmente los parámetros del modelo de infiltración) para forzar el ajuste de los resultados del modelo hidrometeorológico a los del análisis estadístico en las estaciones de aforo. Esto conlleva que estos modelos sólo puedan aplicarse a cuencas de tamaño no muy grande en las que puedan asumirse las hipótesis realizadas.

• Otra posibilidad es el empleo de modelos estadísticos basados en relacionar mediante ecuaciones de regresión múltiple, bien los propios cuantiles calculados a partir de las series temporales de las estaciones de aforo, bien determinados estadísticos de dichas series temporales, con determinadas características fisiográficas y climáticas de las cuencas, lo que permite realizar estimaciones de los cuantiles en aquellos puntos no aforados simplemente conociendo el valor en sus cuencas de las características fisiográficas y climáticas empleadas en la regresión.

El esquema de cálculo a adoptar es una combinación de ambos procedimientos. Modelos hidrometeorológicos en las cuencas de menor tamaño (hasta unos 500 km2 de cuenca) y modelos estadísticos en aquellos puntos con mayor cuenca vertiente.

FIGURA 5: Imagen de la aplicación CAUMAX en el eje del río Ebro en las inmediaciones de la zona de

estudio.

Identificación de regiones homogéneas.

Previamente al análisis, se ha realizado un estudio para la identificación de regiones estadísticamente homogéneas. Las regiones homogéneas se han identificado siguiendo






11

criterios geográficos, a partir de la orografía de la cuenca, la altitud, el régimen de lluvias y el valor de los cuantiles de precipitación más altos. Este método tiene la ventaja de facilitar la asignación de una región a un punto determinado, ya que las regiones están delimitadas perfectamente por los accidentes geográficos de las cuencas.

La identificación de una región como homogénea debe ser comprobada mediante medidas que analizan la homogeneidad estadística de las muestras. Se han utilizado dos tests para cuantificar el grado de homogeneidad de las regiones, el test de Wiltshire, basado en los momentos ordinarios, y el test de Hosking y Wallis, basado en los L-momentos. En función de los resultados obtenidos, se ha decidido adoptar un valor regional únicamente para el coeficiente de sesgo, determinando el coeficiente de variación a partir de la información local.

FIGURA 6: Regiones identificadas por CAUMAX en la cuenca del Ebro.

La zona objeto de estudio, correspondiente al eje del Ebro, queda incluida en la región 96.

Selección y análisis de las series temporales de caudales máximos anuales.

Criterios de selección de los puntos aforados.

En la selección de las estaciones de aforo a emplear en el estudio estadístico es necesario tener en cuenta dos aspectos contrapuestos. Por una parte, es importante disponer para el análisis de series temporales con el mayor número de datos posible puesto que, como se indicó anteriormente, la precisión en la estimación de los cuantiles, en mayor medida en los de mayor periodo de retorno, depende de forma crítica del número de datos disponibles. Esto nos llevaría a seleccionar aquellas estaciones de medida con una mayor longitud de registro o, al menos, a exigir a las series temporales un número de datos mínimo.






12

Por otra parte, teniendo en cuenta que los puntos de aforo seleccionados se emplearán para extrapolar resultados a los puntos no aforados, permitiendo así estimar los cuantiles en puntos sin información, es conveniente que cubran de forma más o menos homogénea todo el territorio de la zona de estudio, evitando dejar zonas sin información donde la extrapolación de resultados a los puntos de la red fluvial no aforados sería menos precisa. Esta circunstancia nos llevaría a seleccionar el mayor número de estaciones de medida posible.

Obviamente, la selección final se realizará como un equilibrio entre los dos condicionantes anteriores, imponiendo a las estaciones seleccionadas una longitud mínima de registro que permita realizar estimaciones de los cuantiles suficientemente precisas, y que permita seleccionar el mayor número de estaciones posible de tal manera que los cuantiles obtenidos en el análisis estadístico tengan suficiente representatividad espacial.

También con el objetivo de aumentar el número de estaciones disponible, y aunque el dato realmente necesario para el análisis es el caudal instantáneo, es decir, el máximo caudal alcanzado durante las crecidas, se ha decidido seleccionar también estaciones en las que en todo o en parte de su registro solo se disponía del caudal medio diario. En estos casos se han transformado los caudales medios diarios en sus correspondientes caudales instantáneos mediante el procedimiento que se expondrá más adelante. Esta decisión ha llevado a considerar también las series temporales medidas en los embalses (en los que únicamente se mide el caudal medio diario) y no sólo las registradas en las estaciones de aforo.

Otro aspecto a considerar es la posible alteración del funcionamiento natural de la cuenca por causas antrópicas, lo que provoca que el caudal medido sea el resultado de dichas alteraciones y que, por tanto, difiera del que habría circulado por el río de no existir dichas causas. Dado que el objeto del análisis estadístico son las series de caudales máximos anuales, es decir, el máximo caudal que ha circulado por el río cada año, los caudales de mayor interés para el estudio son caudales importantes cuya alteración significativa no es, en general, fácil a pesar de la importante infraestructura hidráulica que pudiera existir en las cuencas. En general, la única infraestructura hidráulica capaz de alterar de forma significativa los caudales de avenida son las presas de cierta importancia, por lo que es ese tipo de infraestructura la que hay que tener en cuenta para analizar la posible alteración del régimen natural de la cuenca.

Finalmente, el criterio empleado para seleccionar las estaciones de aforo y los embalses a considerar en el análisis estadístico ha sido el siguiente:

• Estaciones de aforo o embalses cuyo registro no afectado por la presencia de embalses aguas arriba tenga una longitud mínima de 20 datos anuales de caudal máximo medio diario o de caudal instantáneo.

• En zonas con poca densidad espacial de estaciones se ha rebajado la longitud mínima anterior, seleccionando también estaciones con una longitud mínima no alterada de entre 15 y 20 datos.






13

En aquellas estaciones con presencia de embalses aguas arriba se ha considerado que la alteración del registro de caudales máximos anuales no es significativa si se cumplía al menos una de las dos condiciones siguientes:

• El porcentaje de cuenca controlada por los embalses es menor del 10 %.

• El volumen total de embalse es menor del 10 % del volumen medio de las avenidas.

FIGURA 7: Estaciones y embalses seleccionados en una primera fase en la cuenca del Ebro.

Revisión de la calidad y representatividad de los datos de las series temporales.

Un aspecto que condiciona de forma crítica, junto con la longitud de la serie temporal, la precisión de los cuantiles obtenidos mediante el análisis estadístico es la calidad de los datos.

Por este motivo, se han llevado a cabo diversos procesos de revisión y contraste de la información para detectar aquellas estaciones con mayores problemas (en la totalidad o en parte de su registro) y proceder a su corrección o, de no ser posible, a eliminar la estación del análisis.

A continuación se detallan los procesos de revisión llevados a cabo:

• Identificación de tendencias temporales en las series

o Una de las hipótesis básicas de la que parten los modelos estadísticos de máximos es suponer que cada uno de los datos de la serie corresponden a resultados independientes del fenómeno aleatorio objeto de estudio. De esta forma, los datos de las series temporales de máximos deberían ser independientes entre sí y en consecuencia, no deberían presentar ningún tipo de tendencia temporal. De hecho, la construcción de la serie






14

temporal seleccionando únicamente el mayor caudal de cada año tiene como uno de sus objetivos garantizar (o casi garantizar) la independencia de los datos, puesto que los máximos caudales ocurridos en años hidrológicos diferentes, separados generalmente por intervalos temporales importantes, muy difícilmente serán dependientes entre sí.

o Para identificar estas posibles tendencias temporales se ha aplicado el test estadístico de Mann-Kendall. Esto ha permitido detectar algunas estaciones de aforo con problemas en sus datos, habitualmente motivados por la existencia dentro de la serie temporal de dos o más periodos de cierta extensión con caudales de magnitud muy diferente entre sí. Este proceso ha llevado a corregir parte del registro temporal de algunas estaciones o a descartarlo en aquellos casos que esa corrección no ha sido posible.

o La existencia de periodos temporales marcadamente diferentes dentro de las series tiene generalmente su origen en cambios en la curva de gasto a lo largo del tiempo, realizados con intención de ajustarla a nueva información sobre aforos directos. Estas revisiones son imprescindibles para adaptar las curvas a los pequeños cambios que se producen en los cauces como consecuencia de los procesos de erosión y sedimentación, pero en algunos casos esta mejora del ajuste en la parte baja de la curva provoca, al extrapolar, cambios importantes en la parte alta que en general no están justificados.

• Identificación de inconsistencias entre los datos de distintas estaciones

o Este proceso ha consistido en contrastar los datos de caudal máximo anual (tanto medio diario como instantáneo) de cada estación de aforos con los de aquellas otras estaciones situadas aguas arriba de la misma para los años de registro común. Aquellos casos en los que se han detectado medidas de caudal aguas arriba superiores a las registradas aguas abajo han sido objeto de un estudio detallado posterior para determinar las posibles causas de las inconsistencias.

• Identificación de datos anómalamente altos (outliers)

o Este proceso no pretende detectar deficiencias en la calidad de la información sino localizar aquellos datos singulares dentro de las series temporales poco representativos desde un punto de vista estadístico a la hora de realizar el ajuste de la ley de frecuencia.

o De forma general, la frecuencia con la que se mide en una estación de aforos un determinado valor de caudal está relacionada con su probabilidad de ocurrencia, no siendo probable, por tanto, registrar dentro de un periodo temporal más o menos corto (lo habitual en la mayoría de las estaciones de aforo) caudales de elevada magnitud que sólo ocurren durante avenidas de alto periodo de retorno. Sin embargo, es posible que en algún caso, por azar, se registren caudales extraordinariamente altos cuya probabilidad de ocurrencia dentro del periodo temporal de registro






15

sea muy pequeña. Lógicamente, la inclusión de estos datos en el análisis no representa adecuadamente las características estadísticas del proceso aleatorio, y supone alterar significativamente la magnitud de los estadísticos de la serie.

o Para detectar los datos anómalamente altos se ha empleado el procedimiento recomendado por el Water Resources Council (WRC) de EE.UU. (USWRC, 1981). Los datos identificados como anómalos mediante la aplicación del test anterior han sido eliminados de la serie sistemática.

• Identificación de estaciones discordantes

o Tras la identificación de tendencias temporales y outliers en las series sistemáticas de datos observados, se identifican aquellas muestras que tienen unas características estadísticas significativamente diferentes al resto de las estaciones. La identificación se realiza mediante la medida de discordancia (D), basada en el análisis de los valores de los L-momentos (Hosking y Wallis, 1997). El análisis se debe realizar sobre un grupo considerable de estaciones, por ejemplo, una cuenca entera.

o La medida D se fundamenta en la representación de los valores de los coeficientes de LMomentos (t2, t3 y t4) de cada estación en un gráfico en tres dimensiones, obteniendo una nube de puntos. Se consideran como discordantes, aquellos puntos que se alejen de la nube o que tengan una distancia excesiva respecto de los valores medios de la misma.

Selección de puntos aforados

Como consecuencia del proceso anteriormente descrito, de depuración y filtrado de posibles errores y datos anómalos en las estaciones que cumplían las condiciones de partida, se obtiene un conjunto de puntos de control formado por estaciones de aforo y embalses, cuyos registros constituyen la información seleccionada para el presente trabajo. En la siguiente figura aparecen representados los puntos de control seleccionados localizados en el territorio.






16

FIGURA 8: Puntos de aforo seleccionados en la cuenca del Ebro.

Función de distribución y procedimiento de ajuste seleccionado

Como resultado de los estudios y análisis anteriores, se ha decidido realizar los ajustes de las leyes de frecuencia empleado la función de valores extremos generalizada (GEV) ajustada mediante el método de los L-momentos (L-Mom), asumiendo un valor regional para el Lcoeficiente de sesgo en cada una de las regiones estadísticamente homogéneas, y tomando el valor local para el L-coeficiente de variación.

Como caso especial, se ha decidido utilizar una función Gumbel (caso particular de la GEV) ajustada por el procedimiento de los L-momentos en el eje principal del Ebro hasta la confluencia con el Segre.

Estimación de los cuantiles en los puntos de la red fluvial no aforados.

La estimación de los cuantiles en los puntos de la red fluvial no aforados, es decir, en ausencia de medidas locales, puede ser abordada mediante la aplicación de los modelos estadísticos regionales, y puede ser la única opción viable en aquellos puntos con una importante área vertiente donde el uso de los modelos hidrometeorológicos puede no ser adecuado. Estos métodos suelen constar de tres fases principales en su aplicación:

• Identificación de la región con comportamiento estadístico homogéneo donde debe incluirse el punto de estudio.

• Establecimiento de relaciones entre distintas características (fisiográficas, climáticas, edafológicas, etc.) de la cuenca vertiente al punto (variables explicativas) y los estadísticos o cuantiles a estimar, siendo frecuente la aplicación de ecuaciones de regresión logarítmicas.

• Aplicación de las relaciones anteriores en los puntos sin registro.






17

Los distintos métodos propuestos pueden clasificarse en dos grupos según el tipo de estadístico que se seleccione para realizar la extrapolación de información a los puntos no aforados:

• Métodos que relacionan directamente las estimaciones locales de los cuantiles realizadas en los puntos aforados con las variables características de las cuencas vertientes. Este modelo ha sido ampliamente utilizado en EE.UU. por el U.S. Geological Survey (USGS).

• Modelos que asumen una determinada función de distribución regional que puede ser particularizada al punto de estudio a partir de determinados estadísticos (como la media o el coeficiente de variación) que son los que se relacionan con las variables características de las cuencas. Como caso particular se encuentran los modelos de “índice de avenida”, que toman como único estadístico la media de la serie. Para la elaboración del mapa de caudales máximos se ha empleado el primer tipo de métodos por presentar las siguientes ventajas frente a los modelos que asumen una función de distribución regional:

o Capacidad para considerar en el análisis la información sobre avenidas históricas.

o Capacidad para considerar la correlación entre estaciones.

o Amplia experiencia en su aplicación a escala nacional. Como se ha dicho anteriormente, es el método empleado por el USGS para el cálculo de caudales máximos en EE.UU.

Bases teóricas de los modelos regionales basados en ecuaciones de regresión múltiple

Los modelos regionales basados en el análisis de regresiones lineales múltiples tratan de determinar qué variables hidrológicas características de la cuenca vertiente explican la variabilidad del estadístico que se quiere regionalizar, en este caso los cuantiles de caudal máximo. Para ello desarrollan ecuaciones de regresión lineal múltiple entre los cuantiles de un determinado periodo de retorno y las variables seleccionadas. Dichas ecuaciones pueden ser aplicadas a puntos no aforados, de tal forma que permiten estimar los cuantiles a partir de las variables características de la cuenca, variables que son conocidas en cualquier punto.

Selección y ajuste del modelo de regresión múltiple

El ajuste de las ecuaciones de regresión múltiple se ha realizado para cada una de las zonas con similar comportamiento estadístico, las mismas que se obtuvieron anteriormente como resultado de la aplicación de los test de homogeneidad.

Estas regiones han demostrado tener, en general, un buen comportamiento a la hora de ser empleadas para determinar los modelos regionales de regresión, permitiendo la definición de modelos con errores moderados y elevados coeficientes de determinación. De forma complementaria, se ha realizado un análisis de la distribución espacial de los residuos de la estimación, es decir, de las diferencias entre los cuantiles estimados por el modelo y los reales, con la finalidad de detectar posibles sesgos que indicasen la






18

necesidad de subdividir alguna de las regiones consideradas. De esta forma, se ha decidido utilizar la subdivisión de la región 94 para obtener un mejor ajuste de las ecuaciones de regresión.

Antes de abordar el análisis de los modelos regionales de regresión lineal múltiple que nos permitan estimar los cuantiles de caudal máximo anual en los puntos no aforados, es necesario seleccionar y calcular aquellas variables descriptoras de las características de las cuencas que puedan tener mayor influencia en el valor de los cuantiles y, por tanto, que sean más adecuadas para explicar la variabilidad de los mismos entre los distintos puntos de medida seleccionados para el análisis.

Las variables seleccionadas en una primera fase fueron las siguientes:

• Características físiográficas de la cuenca: Área de la cuenca, área drenante con altitud mayor de 1500 m, perímetro de la cuenca, longitud del río principal, pendiente del río principal, altitud media de la cuenca, pendiente media de la cuenca.

• Características climáticas: Cuantiles de precipitación máxima diaria, precipitación media anual.

• Características del suelo y de los usos del suelo: Parámetro P0, tasa mínima de infiltración.

De entre las variables inicialmente consideradas, es necesario seleccionar aquéllas que finalmente se van a incluir en el análisis de los modelos de regresión como variables explicativas, y deben ser aquellas que presenten una mayor correlación con la variable dependiente y, por tanto, permitan explicar mejor su variabilidad, y que al mismo tiempo no tengan una correlación significativa entre sí. Dichas correlaciones se establecerán, lógicamente, entre los logaritmos de las variables, al ser logarítmicos los modelos de regresión múltiple considerados.

Las variables explicativas finalmente seleccionadas teniendo en cuenta las consideraciones anteriores han sido las siguientes:

• Área drenante: A

• Área drenante con altitud mayor de 1500 m: A1500

• Cuantiles de precipitación máxima diaria: Px (x= periodo de retorno del cuantil de precipitación).

• Altitud media de la cuenca: H

• Pendiente media de la cuenca: Pcuenca

• Umbral de escorrentía: P0

• Tasa mínima de infiltración: Inf






19

El cálculo de las variables se ha realizado a partir de la información contenida en las siguientes coberturas:

• Modelo Digital del Terreno del Servicio Geográfico del Ejército y mapa de direcciones de drenaje elaborado por el Centro de Estudios Hidrográficos del CEDEX a partir del mismo con resolución de 500mx500m (figura 9).

FIGURA 9: Modelo digital del terreno

• Mapas de precipitaciones elaborados a partir del trabajo de Máximas Lluvias Diarias en la España Peninsular con resolución de 1000mx1000 m (figura 10).

FIGURA 10: Mapa de precipitaciones máximas diarias

• Mapa de umbral de escorrentía (P0), elaborado siguiendo la metodología utilizada en la Tesis doctoral “Análisis de nuevas fuentes de datos para la estimación del parámetro número de curva del modelo hidrológico del SCS: Datos de perfiles de suelos y teledetección” (Ferrer, 2003), a partir del mapa de






20

usos de suelo Corine Land Cover 2000 y datos de infiltración actualizados, con resolución de 500mx500 m (figura 11).

FIGURA 11: Mapa de umbral de escorrentía en condiciones medias de humedad.

De entre todas las variables características de las cuencas que presentan una correlación significativa con la variable dependiente y una correlación pequeña entre sí, es necesario seleccionar aquél conjunto de variables que mejor explican la variabilidad de los cuantiles. La decisión sobre cuál de entre todos los modelos posibles presenta un mejor comportamiento se toma comparando el valor que adoptan distintos estadísticos en cada uno de los modelos. Los estadísticos utilizados para realizar ese análisis son habitualmente los siguientes:

• Coeficiente de determinación (R2): este estadístico mide la proporción de variación de la variable dependiente que es explicada por el modelo.

• Coeficiente de determinación ajustado (R2 corr): representa una corrección del estadístico anterior teniendo en cuenta el número de estaciones incluidas en el análisis y el número de variables explicativas adoptadas por el modelo.

• Error estándar de regresión (e): Es una medida de la diferencia entre los valores reales de los cuantiles en cada estación y los valores estimados por los modelos.

• Suma de los errores cuadráticos de predicción (PRESS): Este estadístico se calcula realizando el ajuste de la ecuación de regresión con todas las estaciones menos una, aplicando la ecuación deducida a la estación no incluida en el análisis y calculando el error cuadrático respecto al valor real del cuantil en dicha estación. Este proceso se repite para cada una de las estaciones, calculando finalmente la suma de los errores.

• Coeficiente de Mallow (Cp): El coeficiente de Mallow representa un compromiso entre maximizar la proporción de varianza de la variable dependiente que es explicada por el modelo incluyendo todas las variables relevantes, y minimizar el error de la estimación reduciendo el número de variables empleadas.






21

El modelo finalmente seleccionado debe ser aquél con un menor valor del error cuadrático medio, del coeficiente de Mallow y del error cuadrático de predicción, y un mayor valor del coeficiente de determinación ajustado. En cualquier caso, el modelo seleccionado no sólo debe tener un buen comportamiento desde un punto de vista matemático, sino que, además, la ecuación debe tener sentido físico, reflejando relaciones funcionales entre los cuantiles y las variables que representen el comportamiento real de la cuenca.

Los modelos finalmente seleccionados para cada una de las regiones y cada uno de los periodos de retorno se muestran en la tabla adjunta.






22

Modelo de regresión para los ejes principales

Para obtener los caudales en los puntos de las corrientes principales correspondientes a grandes cuencas vertientes no se han empleado las ecuaciones de regresión múltiple descritas en el apartado anterior, ya que las relaciones entre los caudales y las características medias de las cuencas podrían diferir de las correspondientes a cuencas de menor tamaño e, incluso, el valor medio en la cuenca de alguna de las variables podría no ser significativo a la hora de determinar el caudal cuando se trata de cuencas






23

de gran tamaño. De esta manera, y puesto que se dispone de varias estaciones de aforos con series suficientemente largas y fiables a lo largo de los ejes principales, algunas de ellas con información histórica, se ha decidido establecer simplemente unas relaciones sencillas en función del área vertiente al punto que permitan interpolar y dar continuidad a los cuantiles calculados en las distintas estaciones de aforo presentes a lo largo del río.

En la figura 12 se muestran las relaciones entre los logaritmos de los cuantiles de caudales y el logaritmo del área de la cuenca. Se observa cómo los ajustes son muy buenos, con coeficientes de correlación próximos a uno, por lo que se ha considerado que estas relaciones sencillas en función del área son suficientes para dar continuidad a los valores de los cuantiles de caudal máximo a lo largo de la corriente principal del río.

FIGURA 12: Relaciones entre los cuantiles de caudal máximo y el área vertiente a lo largo de la corriente

principal de los grandes ríos para T= 100 años






24

FIGURA 13: Modelos seleccionados para la región 96 (eje del Ebro) y factor de ajuste obtenido

1.2.2.3. Estimación de la Máxima Crecida Ordinaria con CAUMAX.

Se puede realizar una estimación de la máxima crecida ordinaria, si se conoce la ley de frecuencia de caudales máximos en situación natural, aplicando la metodología expuesta en el informe “Aspectos prácticos de la definición de la máxima crecida ordinaria” (CEDEX, 1994).

En el citado informe, se realizó un estudio experimental en 20 puntos distribuidos a lo largo de la geografía nacional, en los que se disponía de la serie histórica de máximos caudales anuales, así como de información sobre el caudal de desbordamiento del cauce. El caudal de desbordamiento se considera dentro de la hidráulica fluvial como representativo del comportamiento de la corriente, y su línea de agua marca los límites de lo que comúnmente se entiende por cauce. Además, se pudo comprobar en distintos puntos de la geografía española que dicho caudal está dentro del rango de los valores compatibles con los términos en que el Reglamento del Dominio Público Hidráulico define como máxima crecida ordinaria, y además es indicativo del posicionamiento central de los mismos, motivos por los cuales se propuso identificar la máxima crecida ordinaria con el caudal de desbordamiento.

Una expresión, para la que es necesario conocer la ley de frecuencia de caudales máximos, es la siguiente en la que se relaciona el periodo de retorno correspondiente al caudal de la máxima crecida ordinaria (TMCO) con el coeficiente de variación de la serie temporal:

Se ha realizado una regionalización estadística de la cuenca de estudio en la que se conocen los valores regionales de los principales estadísticos, entre ellos el coeficiente de variación para estimar el periodo de retorno y con él, el caudal en los distintos puntos de la red fluvial.

En la figura 14 se recogen los valores regionales del coeficiente de variación en las regiones estadísticas adoptadas en el estudio, así como los valores correspondientes para el periodo de retorno aproximado de la máxima crecida ordinaria.






25

Región CV T

91 0,47 2,5

92 0,70 3,5

93 1,36 7,0

94 1,04 5,0

95 0,69 3,5

96 0,50 2,5

FIGURA 14: Coeficientes de variación regionales y periodos de retorno para la máxima crecida ordinaria en las regiones estadísticas. En rojo, región estadística de la zona de estudio.

1.2.2.4. Resultados de CAUMAX.

A partir de la aplicación de la metodología expuesta se obtienen los siguientes resultados en régimen natural para el tramo de estudio:

PERIODO DE RETORNO (Años) 2 2.5 (M.C.O.) 10 50 100 500

Q (m3/s) 911 1.206 1.748 2.308 2.553 3.105

Como se observa en la tabla el caudal correspondiente a la máxima crecida ordinaria corresponde a un período de retorno de 2.5 años y un caudal de 1.206 m3/s.

Es importante señalar que el período de retorno asociado a la MCO, en este caso 2,5 años, es el resultado del exhaustivo análisis descrito anteriormente y conceptualmente significa que como media se produce cada 2,5 años un caudal de magnitud similar (o superior) a 1.206 m3/s.

Este caudal, por otra parte, es superior al que se obtiene con la definición legal de MCO, media de 19 años consecutivos de los caudales máximos instantáneos anuales, que en nuestro para la estación de aforos de Mendavia se corresponde con 912,7 m3/s, lo cual queda perfectamente justificado por el efecto laminador que introducen los embalses en el caudal circulante, ya que como se ha indicado, los caudales registrados en la estación de aforos están en régimen alterado (no natural).

1.2.3. CAUDAL ADOPTADO PARA LA MCO. CONCLUSIONES.

De todo cuanto antecede se puede concluir lo siguiente:

1. Para el cálculo del caudal correspondiente a la MCO se ha aplicado rigurosamente la definición incluida en el reglamento y se ha contrastado con el CAUMAX, metodología desarrollada por el CEDEX (Centro de Estudios y Experimentación de Obras Públicas) a nivel nacional, que está validada y contrastada y es una referencia en esta materia.






26

2. Los caudales obtenidos por ambas metodologías (estudios foronómicos y CAUMAX) son coherentes si se tiene en cuenta la restitución a régimen natural de CAUMAX, ya que con régimen alterado el valor es algo menor debido al efecto laminador de las presas situadas aguas arriba, efecto más significativo cuanto de menor periodo de retorno sea la avenida.

3. Se ha adoptado el caudal de 913 m3/s para la máxima crecida ordinaria, inferior al proporcionado por el CAUMAX en régimen natural y equivalente tanto a la media de los caudales máximos instantáneos anuales de la serie representativa de diecinueve años adoptada para la estación de aforos del Ebro en Mendavia (1950/51 a 1968/69) como a la media de los caudales máximos instantáneos anuales registrados en dicha estación. Este valor en cualquier caso es un umbral inferior para la Máxima Crecida Ordinaria; es decir, las zonas inundables obtenidas en el modelo hidráulico serán en todo caso iguales o inferiores a las probables.




ANEXO Nº 3. ESTUDIO HIDRÁULICO


ANEXO Nº3.

ESTUDIO HIDRÁULICO






ÍNDICE

1. INTRODUCCIÓN. ........................................................................................ 1

2. ESTUDIOS CARTOGRÁFICOS. ..................................................................... 1

2.1. INTRODUCCIÓN. .................................................................................................. 1

2.2. METODOLOGÍA. ................................................................................................... 2

2.2.1. ADQUISICIÓN DE LA INFORMACIÓN LiDAR ............................................................... 2 2.2.1.1. Introducción a la tecnología LiDAR. .................................................................. 3 2.2.1.2. Resolución y precisión del LiDAR. ....................................................................... 5 2.2.1.3. Clasificación, edición, corrección y tratamiento de datos LiDAR .............. 6 2.2.1.4. Control de calidad de la información LiDAR suministrada ........................... 6 2.2.1.5. Clasificación automática .................................................................................... 7 2.2.1.6. Clasificación manual ............................................................................................ 7 2.2.1.7. Elaboración de productos finales ...................................................................... 8 2.2.1.8. Control de calidad de los productos .............................................................. 11 2.2.1.9. Elaboración de la cartografía vectorial ......................................................... 11

2.2.2. CROQUIS DE OBRAS DE FÁBRICA .............................................................................. 11 2.2.3. PERFILES BATIMÉTRICOS ............................................................................................... 12

3. ESTUDIOS HIDRÁULICOS. .......................................................................... 14

3.1. INTRODUCCIÓN ................................................................................................. 14

3.2. METODOLOGÍA .................................................................................................. 14

3.2.1. SELECCIÓN DEL MODELO. ......................................................................................... 14 3.2.2. INFORMACIÓN DE PARTIDA ....................................................................................... 14 3.2.3. MODELOS HIDRÁULICOS BIDIMENSIONALES ........................................................... 15 3.2.4. ESTIMACIÓN INICIAL DE LAS PÉRDIDAS DE CARGA ............................................... 20 3.2.5. CALIBRACIÓN Y VALIDACIÓN DEL MODELO. ......................................................... 23 3.2.6. RESULTADOS OBTENIDOS. ........................................................................................... 24

4. LAMINAS ................................................................................................... 26






1

1. INTRODUCCIÓN.

Este anexo se ha dividido en dos apartados. En primer lugar, se describen los trabajos cartográficos, ya que es fundamental disponer de una buena cartografía del terreno, que junto con los caudales son los datos a introducir en el modelo hidráulico. A continuación, se realiza una descripción de la metodología utilizada para realizar el estudio hidráulico.

2. ESTUDIOS CARTOGRÁFICOS.

2.1. INTRODUCCIÓN.

Consiste en obtener la topografía del terreno del tramo fluvial a modelizar en una anchura suficiente para representar toda la zona inundable. Es necesario disponer de un modelo del terreno que incluya no sólo la zona del cauce (dominio público hidráulico) por donde circulan las crecidas ordinarias, sino también hay que tener bien representada topográficamente la llanura de inundación.

Para la realización de estos trabajos, el Ministerio de Agricultura, Alimentación y Medio Ambiente suscribió convenio con el Instituto Geográfico Nacional. El Instituto Geográfico (dependiente del Ministerio de Fomento) realizó un vuelo con tecnología LiDAR (altimetría láser, datos: x y z), obteniendo coordenadas cartográficas de, al menos, 0.5 puntos por cada metro cuadrado que ha permitido elaborar un modelo digital del terreno de detalle. La obtención de cada punto de información (posición de coordenadas x, y, además de su altitud, z) se realiza mediante una tecnología láser que permite desde un avión, cuya posición está en todo momento georreferenciada, emitir rayos láser cuyo eco es (tras rebotar con la superficie terrestre) capturado de nuevo. De este modo se obtienen datos del terreno, vegetación, edificios, etc., presentes en la superficie terrestre. Esta cartografía digital ha servido de base posteriormente para representar el terreno (cauce de los ríos y sus márgenes) y fue completada con visitas de reconocimiento de campo para su adecuada calibración y representación de los detalles.

En la zona de estudio, entre los años 2016 y 2018 se ha llevado a cabo un segundo vuelo LiDAR que ha permitido elaborar un modelo digital de terreno de mayor resolución (segunda cobertura). Ello ha sido consecuencia de que, en esta segunda cobertura, fueron 2 puntos por cada metro cuadrado los obtenidos en territorio riojano (margen derecha del Ebro) y 14 en el territorio foral de Navarra (margen izquierda del Ebro).

La precisión de la cartografía obtenida tras la realización de cualquiera de estos trabajos supera a la topografía clásica de escala 1:1.000 exigida en la tramitación de un procedimiento de deslinde según el artículo 242.3 del Reglamento de Dominio Público Hidráulico, por lo que es válida para los trabajos incluidos en este estudio.






2

FIGURA 1: Localización del área estudiada

2.2. METODOLOGÍA.

Una de las informaciones base para la delimitación del Dominio Público Hidráulico es la cartografía del tramo de estudio, es decir, el Modelo Digital del Terreno (MDT), que represente fielmente la realidad del terreno en el tramo de estudio.

Para ello, ha sido utilizado un modelo digital del terreno generado mediante la tecnología LiDAR (Laser Imaging Detection and Ranging), el cual debe ser tratado posteriormente en gabinete para eliminar los valores correspondientes a elementos diferenciados del terreno: vegetación, edificaciones, puentes, etc.

En los siguientes apartados se realiza una breve descripción de los procedimientos seguidos para la generación de la cartografía.

2.2.1. ADQUISICIÓN DE LA INFORMACIÓN LiDAR

La información LIDAR empleada para la elaboración de la cartografía necesaria para la ejecución de los trabajos ha procedido de los vuelos realizados por el Plan Nacional de Ortofotografía Aérea (PNOA) entre los años 2016-2018.

El PNOA es un proyecto cofinanciado y cooperativo entre la Administración General del Estado (AGE) y las comunidades autónomas que se enmarca dentro del Plan Nacional de Observación del Territorio (PNOT), siendo coordinado por el Instituto Geográfico Nacional (IGN) y el Centro Nacional de Información Geográfica (CNIG). Tiene como objetivo la obtención de productos fotogramétricos comunes para todo el territorio nacional, con especificaciones técnicas consensuadas entre todas las administraciones implicadas y cuyo resultado es el proceso y obtención de ortofotografías digitales, incluyendo: el vuelo fotogramétrico, apoyo de campo, aerotriangulación y el modelo digital de elevaciones.






3

Desde el año 2004, el Plan Nacional de Ortofotografía Aérea proporciona imágenes aéreas, ortofotos y modelos digitales de elevaciones (MDE) de todo el territorio, con una periodicidad de actualización de dos años. El proyecto se encuentra en continua evolución, adaptándose a las necesidades de los usuarios y al desarrollo de nuevas tecnologías. En 2009, se planteó la necesidad, por parte de una serie de usuarios, de la obtención de Modelos Digitales de alta precisión, obtenidos por tecnología LiDAR, para la realización de cartografía de áreas de inundación, proyectos de carreteras, inventarios forestales, etc.

Tal y como se observa en la siguiente ilustración, la disponibilidad de datos LIDAR dentro de la demarcación hidrográfica del Ebro no sólo es completa, sino que en buena parte de ella, incluida la zona estudiada, ya se dispone de la segunda cobertura LiDAR:

FIGURA 2: Localización del área estudiada El vuelo y la obtención de cartografía para la zona de estudio se llevó a cabo entre los meses de agosto y septiembre de 2016 para el área perteneciente a la Comunidad Autónoma de La Rioja (margen derecha del Ebro) y de septiembre a noviembre de 2017 en los terrenos adscritos a la Comunidad Foral de Navarra (margen izquierda del Ebro).

2.2.1.1. In troducc ión a la tecnología LiDAR.

El LiDAR (Light Detection And Ranging) es un sistema que contiene un emisor de impulsos y un receptor de detección remota basado en un sensor láser. El instrumental se instala en un avión desde el cual se procede a la emisión de un pulso láser y a la medida del tiempo que tarda dicho pulso en llegar a un objeto, en el cual rebota, y volver al avión. De este modo, a cada tiempo de respuesta se le asigna una cota. El sistema instalado en el avión no trabaja solamente en la vertical del mismo. También se realizan medidas transversalmente a la dirección de avance, con un cierto ángulo respecto a la vertical.

Según las características del objeto al que llega el pulso láser se pueden dar diferentes respuestas:

• En una superficie sólida (edificios, suelo, etc.) el rayo se refleja y vuelve al avión.

• En el agua, el rayo láser es absorbido y no vuelve al avión, por lo que no se obtiene información.






4

• En vegetación, al chocar con la copa del árbol parte del rayo se refleja y vuelve al avión, pero parte del rayo la atraviesa hasta llegar al suelo. En este caso, el sistema guarda el primer y último pulso.

FIGURA 3: Toma de datos de LiDAR

El registro del primer pulso permite describir la superficie (parte más alta) de los objetos, mientras que el registro del último eco se usa para conocer la superficie del suelo. La capacidad multirretorno se utiliza para detector vegetación y su altura.

Las principales características de los datos LiDAR facilitados por el PNOA son los siguientes:

• Los datos se facilitan en ficheros digitales (formato LAS) con información altimétrica de la nube de puntos LiDAR, distribuidos en ficheros de 2x2 km de extensión.






5

• Las nubes de puntos han sido capturadas mediante vuelos con sensor LiDAR con una densidad de 2 (La Rioja) o 14 (Navarra) puntos/m2, y posteriormente clasificadas de manera automática y coloreadas mediante RGB obtenido a partir de ortofotos del Plan Nacional de Ortofotografía Aérea (PNOA), con tamaño de pixel de 25 o 50 cm.

• Distancia nominal entre puntos: 1,4 m.

• Alturas ortométricas con precisión altimétrica mejor de 20 cm RMSE Z.

• Sistema geodésico de referencia ETRS89 y proyección UTM en el huso correspondiente a cada fichero.

Una vez recibida la información LiDAR correspondiente a los vuelos del PNOA, se procede a la recopilación, organización y clasificación de la documentación según las áreas de trabajo a las que pertenezcan, para facilitar el acceso a la misma en las siguientes fases de procesado de la información.

2.2.1.2. Res oluc ión y prec is ión de l LiDAR.

La resolución del vuelo LIDAR se caracteriza por el número de puntos medio por metro cuadrado del que se dispone de información.

Para el SNCZI se recomiendo emplear densidades medias, y, en especial, la densidad que aporta el LIDAR PNOA del Instituto Geográfico Nacional.

En este sentido, el Instituto Geográfico Nacional comenzó en el año 2009 los trabajos necesarios para obtener el vuelo LIDAR de gran parte del territorio nacional, de densidad promedio 0.5 puntos/m2.

En el año 2010, dentro del marco del proyecto XERRAC, se realizó una comparación de los resultados entre el LiDAR PNOA y un LIDAR convencional con vuelo a 4 puntos/m2. Los resultados, como se puede ver en la figura 4, muestran que, pese a la menor precisión, la definición de elementos es suficiente para caracterizar el terreno. Se puede comprobar que un vuelo a 0,5 puntos/m2 ya sería suficientemente representativo de la realidad.

Para poder referenciar correctamente el punto que se ha medido en el terreno, se utiliza la combinación de dos técnicas diferentes. De esta forma se obtienen las alturas del terreno con un error de +/- 15 cm en altura.






6

FIGURA 4: Comparativa entre LiDAR con vuelo convencional y LiDAR con vuelo PNOA

2.2.1.3. Clas ificac ión , ed ic ión , correcc ión y tra tamiento de da tos LiDAR

El procesado de la información LIDAR del PNOA requiere una serie de tareas hasta obtener los productos finales:

• Control de calidad de la información LIDAR suministrada.

• Clasificación automática.

• Clasificación manual.

• Elaboración de productos finales

• Control de Calidad de los productos

2.2.1.4. Contro l de ca lidad de la in formac ión LiDAR s uminis trada

En algunos casos, la información LIDAR PNOA presenta ciertas anomalías que han tenido que ser subsanadas previamente al procesado de la información.

Los principales controles elaborados en esta fase inicial han sido:

• Control del Sistema de Referencia, verificando que toda la información está en el sistema de referencia designado para el proyecto (ETRS89 UTM 30N).

• Corrección de las líneas de vuelo, ya que en ocasiones aparecen anomalías por una mala unión entre las diferentes líneas de vuelo que dan como resultado saltos de entre 40 cm y 1 m en el terreno.

• Corrección de dichas anomalías previa a la clasificación de la información LAS.






7

• Revisión de los valores de elevación y, en caso de ser necesario, realización de las correcciones oportunas de la misma empleando para ello el geoide.

2.2.1.5. Clas ificac ión automática

En primera instancia, se ha realizado un procesado automático de la nube de puntos que permite obtener su clasificación con un alto porcentaje de efectividad.

Este nivel de detalle en la clasificación puede no ser suficiente para obtener Modelos Digitales del Terreno con la calidad mínima requerida para su uso en simulación hidráulica, razón por la que ha sido necesario aplicar métodos adicionales. A modo de ejemplo, en la figura 5 se muestra un bloque de edificios que no ha sido correctamente clasificado por el proceso de clasificación automática.

FIGURA 5: De izquierda a derecha: imagen de la nube de puntos en 2D, comparativa 3D. Ambas permiten

determinar que se trata de una edificación y no de vegetación (clase asignada por el proceso de clasificación automática)

2.2.1.6. Clas ificac ión manual

La clasificación manual de datos LiDAR se ha realizado tras revisar el proceso e clasificación automático, mediante la edición de las diferentes clases de interés obtenidas en dicho proceso (suelo, edificaciones, vegetación media, alta, baja, etc.).

FIGURA 6: Error de clasificación en edificación subsanado tras aplicar el proceso de clasificación manual






8

2.2.1.7. Elaborac ión de produc tos fina les

Durante esta fase se han utilizado herramientas y soluciones específicas para elaborar los MDT de suelo desnudo y los MDS con todos los elementos clasificados.

Todos los productos raster generados tienen una resolución de 1x1 m.

1. MODELOS DIGITALES DEL TERRENO (MDT): Los MDT básicos contienen únicamente la información del suelo desnudo. Dentro de los productos considerados MDT han sido elaborados tres que contienen diferente información con respecto al terreno.

a) MDT terreno: se corresponde con un MDT que contiene toda la información clasificada como suelo desnudo. Este modelo no conlleva ninguna manipulación especial.

b) MDT cauce natural: consiste en un modelo representativo de la superficie desnuda, pero eliminando del modelo todas las intervenciones antrópicas, como por ejemplo motas, puentes, azudes, presas, canales, etc. Este modelo digital tiene como objetivo poder realizar la simulación hidráulica en régimen natural. La figura 7 muestra los cambios aportados por este modelo en comparación con el MDT de suelo desnudo.

FIGURA 7: Comparativa entre MDT terreno (izquierda), con una protección o mota en una de las márgenes,

y MDT cauce natural (derecha), una vez eliminada la protección gracias al proceso de clasificación empleado

c) MDT terreno 2: se corresponde al MDT de suelo desnudo modificado e incluye también las infraestructuras situadas sobre el cauce, pero eliminando los azudes y presas localizadas dentro del cauce procesado. La figura 8 muestra cómo se ha eliminado un alto de agua del perfil longitudinal del río, mediante procedimientos de mejora del cauce, pero conservando el resto de rasgos existentes.






9

FIGURA 8: Ortoimagen de cauce con azud (arriba). MDT de suelo desnudo con perfil longitudinal del

cauce reseñado (abajo, izquierda). MDT 2 con el salto de agua corregido (abajo, derecha).

2. MODELOS DIGITALES DE SUPERFICIES (MDS): Los MDS recogen diferentes elementos presentes en el terreno y que han sido excluidos de los MDT. Han sido elaborados diferentes tipos de MDS gracias al exhaustivo proceso de clasificación al que fueron sometidos los datos LAS, diferenciándose dos tipos de vegetación (arbustiva y arbórea), edificaciones, puentes, infraestructuras localizadas sobre los cauces (puentes, presas…), motas, etc.

a) MDS arboles: Para un correcto funcionamiento de estos modelos, la vegetación se mantiene en todo el territorio a excepción del cauce, para facilitar la simulación hidráulica. Uno de los aspectos tenidos en cuenta durante la fase de control de calidad de los modelos ha sido la comprobación de la no existencia de vegetación sobre el cauce.

b) MDS arbustos: Al igual que los modelos con árboles, los modelos con arbustos tienen las mismas características y funcionalidades que el MDS árboles. La diferencia existente entre los dos modelos es que éste (MDS arbustos) recoge la vegetación de media y baja altura.

c) MDS edificios: Este modelo de superficies es el resultado de la combinación del MDT de suelo desnudo con los elementos clasificados como edificios y presentan dos tipos de usos y ventajas diferentes:

• Permite simular el efecto de la presencia de edificios, especialmente en zonas muy urbanizadas como núcleos urbanos o zonas industriales.






10

• Son de gran utilidad para calcular la vulnerabilidad de determinadas zonas con alta densidad de edificios, pudiendo valorar el grado en que la inundación le afecta: altura de la lámina de inundación, determinación de zonas con velocidades de flujo inaceptables para el uso residencial, etc.

FIGURA 9: Ejemplo de MDS edificios

d) MDS edificios + infraestructuras: Se trata de un modelo derivado del anterior, que combina el MDT terreno 2 con los elementos clasificados como edificios.

FIGURA 10: Ejemplo de MDS edificios + infraestructuras

3. OTROS PRODUCTOS:

a) Imagen de intensidades: a partir de la amplitud de la señal recepcionada desde el avión después de rebotar en la superficie terrestre, se obtiene una imagen de intensidades que permite realizar distinciones entre superficies, identificando carreteras, agrupaciones industriales, etc. La geodatabase de productos finales incluye la imagen de intensidades del primer pulso.






11

b) Capa de densidad de puntos: ha sido obtenida una capa que permite conocer la densidad de puntos en cada píxel de la zona de estudio. La capa estará en formato raster.

c) Líneas de rotura: se trata de la capa vectorial que recoge la información básica de las infraestructuras existentes en el terreno. Permite generar los modelos generación de los modelos finales con o sin puentes y obtener datos imprescindibles a la hora de acometer la modelización hidráulica.

2.2.1.8. Contro l de ca lidad de los produc tos

Una vez elaborados los productos entregables, se lleva a cabo un exhaustivo control de calidad de los productos obtenidos mediante análisis espaciales, comprobación de perfiles longitudinales y transversales del cauce y verificación del proceso de clasificación.

2.2.1.9. Elaborac ión de la ca rtogra fía vec toria l

Usando como información base los productos LIDAR y sus subproductos derivados (modelos digitales de elevaciones, imagen de intensidades, etc) ha sido elaborada una cartografía vectorial que recoge aquellos elementos vectorizables que son de interés en la generación de modelos de simulación y que son necesarios para su correcto funcionamiento en la mayoría de los programas de simulación hidráulica existentes actualmente en el mercado.

Las capas vectoriales elaboradas han sido:

• Los límites de los cauces naturales presentes en la zona de estudio, representados por dos líneas que recorren ambas márgenes del río.

• Edificaciones presentes en la zona de estudio. Siendo de gran interés aquellas que se encuentran dentro de las potenciales láminas de inundación.

• Los elementos significativos que puedan ser vectorizados y que han sido necesarios para la generación de la simulación hidráulica (barreras naturales o artificiales que actúan como líneas de rotura): viaductos, puentes, pasarelas, terraplenes de carreteras, líneas de coronación y pies de base de las motas del cauce…

2.2.2. CROQUIS DE OBRAS DE FÁBRICA

La caracterización de obstáculos en cauce tiene por objeto obtener los datos necesarios de puentes y azudes con el fin de poder introducirlos en los modelos hidráulicos de cálculo.

Esta caracterización no precisa un detalle a nivel de perfil taquimétrico de la obra de paso, pero sí debe inluir los datos más relevantes del obstáculo que los modelos hidráulicos demandan para su caracterización: anchos libres y anchos “ocupados” del cauce (pilas, estribos), alturas libres para el paso del agua, forma de los vanos, etc.






12

En el caso concreto de este estudio no existe ninguna obra de fábrica en el río que requiera ser incluida para la modelización hidráulica.

2.2.3. PERFILES BATIMÉTRICOS

El objetivo de realizar batimetrías es conseguir información adicional con el fin de completar la información del LiDAR donde esta tecnología no informa adecuadamente.

Las batimetrías proporcionarán datos de perfiles transversales del cauce en forma de puntos con coordenadas (x, y, z). La ubicación de estos perfiles transversales se decide de forma que su elección proporcione una visión lo más cercana a la realidad de la topografía del río. Para ello, la ubicación de los perfiles batimétricos contempla los siguientes criterios generales:

• Caracterización del cauce: los perfiles batimétricos se ubican a una distancia adecuada para proporcionar la topografía más representativa del cauce. La distancia variará entre 100-500 m y los perfiles se situarán perpendiculares al cauce en la medida de lo posible.

• Caracterización de elementos en el cauce: se dispondrán perfiles batimétricos aguas arriba y aguas abajo de azudes y de islas.

• Caracterización de canales secundarios: los perfiles batimétricos se realizarán en el cauce principal y en derivaciones de importancia del mismo que puedan suponer un trasvase de volumen de agua no despreciable.

• Caracterización de otras zonas singulares: se dispondrán perfiles batimétricos en confluencias de importancia u otras particularidades que se puedan dar.

Los datos batimétricos contienen la información de puntos del cauce constreñidos en los perfiles antes señalados y también de puntos de la línea de agua. Es decir, cada perfil dispone una coordenada en cada margen del punto donde se produce la transición agua-tierra. Estos puntos agua/tierra han de obtenerse siempre que sea posible puesto que, en ocasiones, su caracterización será dificultosa por coincidir con zonas de alta vegetación, rocas, etc. Estos puntos sirven para calibrar los datos batimétricos, puesto que se podrán comparar, en la mayoría de los casos, con los puntos obtenidos de LiDAR directamente.

Las técnicas para obtener los datos batimétricos varían en función de las características del cauce, pudiendo tratarse de perfiles obtenidos directamente con jalón cuando puede vadearse el tramo o con ecosonda y embarcación si las características de la corriente lo exigen.

Una vez obtenidos los datos batimétricos, se interpola una superficie batimétrica que se incorpora directamente a los modelos digitales del terreno obtenidos mediante LiDAR. Para ello, en primer lugar, se realiza una la calibración de los datos (con los puntos tierra-agua ya mencionados) con el fin de detectar posibles desfases sistemáticos en la coordenada z de los puntos.

En segundo lugar, se trata la topografía del cauce de forma aislada al resto de puntos del terreno. De esta forma, se trabaja sólo en las zonas contenidas dentro de la zona






13

cauce y se interpolan los datos obtenidos de los perfiles transversales sin tener en cuenta islas o azudes. Con este procedimiento, se obtiene un modelo digital del terreno “limpio” del cauce.

Posteriormente se incorporan elementos elevados dentro del cauce, como pueden ser azudes o islas. Por último, este MDT del cauce se integra en el MDT general de la zona a estudiar, siendo el MDT final homogéneo y adecuado para tratar en los modelos hidráulicos.

Para el caso que nos ocupa el MDT final se ha obtenido utilizando los vuelos LiDAR de segunda cobertura llevados a cabo entre los meses de agosto y septiembre de 2016 para el área perteneciente a la Comunidad Autónoma de La Rioja (margen derecha del Ebro) y de septiembre a noviembre de 2017 en los terrenos adscritos a la Comunidad Foral de Navarra (margen izquierda del Ebro). Este producto se ha corregido en la zona del cauce a partir del MDT proporcionado por la Confederación Hidrográfica del Ebro, el cual no incorpora información batimétrica.

FIGURA 11: Modelo digital del terreno sin batimetría, sin edificios, sin puentes y sin vegetación de la zona

de estudio.






14

3. ESTUDIOS HIDRÁULICOS.

3.1. INTRODUCCIÓN

Una vez caracterizada la geometría de los cauces y de las llanuras de inundación a través del MDT generado a partir de los datos LiDAR, apoyados por el estudio de fotografías aéreas y la comprobación en campo, y conocida la hidrología del tramo de cauce, el estudio hidráulico persigue la obtención de la extensión de la zona inundable correspondiente a la Máxima Crecida Ordinaria.

3.2. METODOLOGÍA

3.2.1. SELECCIÓN DEL MODELO.

La selección del modelo hidráulico a emplear en cada tramo de estudio depende fundamentalmente de las particularidades de la zona.

El uso de modelos bidimensionales se hace indispensable allí donde las dos componentes de la velocidad, en el sentido del flujo y en el sentido transversal, tienen un peso importante. Estas circunstancias suelen darse en grandes llanuras aluviales o zonas en las que se producen desbordamientos laterales de importancia.

Los modelos unidimensionales se aplican en las situaciones en las que la componente transversal de la velocidad sea nula o despreciable. Estos modelos presentan la ventaja con respecto a los bidimensionales de una mayor rapidez de computación, lo que facilita enormemente el procedimiento de calibración de los parámetros del modelo.

En este caso se ha utilizado un modelo hidráulico bidimensional, ya que es más adecuado para la simulación del flujo transversal en un tramo de río meandriforme y, por lo tanto, proporciona unos resultados más precisos.

3.2.2. INFORMACIÓN DE PARTIDA

La información de partida de cualquier modelo hidráulico debe incluir al menos los siguientes elementos:

• Datos de caudales punta o hidrogramas que resulten del estudio hidrológico.

• Información básica de caracterización física del cauce.

• Información sobre elementos localizados aguas abajo de la zona de estudio que ayude a definir las condiciones de contorno en el modelo hidráulico, como nivel, azudes, zonas en las que se produzca el calado crítico, etc.

En este trabajo, se ha contado con la siguiente información de partida para la confección del modelo hidráulico:

• El modelo digital del terreno (MDT) del tramo fluvial a estudiar, procedente del tratamiento de los datos LiDAR del PNOA y la cartografía disponible.






15

• Ortofotos actuales y fotografías aéreas históricas de la zona de estudio a la mejor resolución posible.

• Información sobre elementos localizados aguas arriba y abajo del tramo de estudio que ayude a definir las condiciones de contorno o de borde de la simulación, como por ejemplo niveles de embalses, niveles de marea, azudes, puentes, etc.

• Información sobre usos del suelo para determinar las pérdidas de energía del agua.

3.2.3. MODELOS HIDRÁULICOS BIDIMENSIONALES

Como principal herramienta para la modelización hidráulica bidimensional se ha utilizado InfoWorks ICM. El modelo InfoWorks ICM ha sido desarrollado por HR Wallingford (Reino Unido) y supone de una potente herramienta para el análisis de riesgos de inundabilidad.

A diferencia de otros modelos bidimensionales, InfoWorks ICM ofrece la posibilidad de incorporar elementos unidimensionales dentro del cálculo y permite la simulación de coberturas, incluso a presión, en combinación con el flujo en superficie; es decir, permite, de manera rápida y sencilla, simular avenidas en zonas urbanas altamente antropizadas, así como valorar la aplicación de medidas estructurales que reduzcan el riesgo de inundación y mejoren la capacidad hidráulica de desagüe del río de una forma eficiente.

FIGURA 12: Vista 3D del mallado con InfoWorks ICM de una zona a estudiar

Una de las principales ventajas de InfoWorks ICM respecto de sus competidores es el algoritmo de mallado propio que incorpora el programa. Una vez definido por zonas el tamaño mínimo y máximo de los elementos triangulares que definirán el terreno, incorporadas las líneas de rotura que identifican, por ejemplo, los taludes del cauce o los terraplenes de infraestructuras, y establecidas las áreas exentas de cálculo, como las edificaciones, el programa construye de manera automática la malla que mejor






16

compatibiliza el coste computacional con la precisión en el cálculo y la representación de la geometría del terreno.

Otra de las grandes ventajas que presenta InfoWorks ICM es que, una vez ejecutada la simulación, el acceso a los resultados hidrodinámicos es sencillo y su interfaz permite obtener:

• Planta de calados y velocidades vectoriales.

• Gráficas de evolución puntual de variables hidráulicas a lo largo del periodo de simulación.

• Integración de resultados sobre una línea definida por el usuario; por ejemplo, flujo a través de una sección.

FIGURA 13: Ejemplos de resultados gráficos de una simulación con InfoWorks ICM

Las características principales de este modelo y la forma en que resuelve las ecuaciones son las que se describen a continuación:

• Discretización espacial: Malla no estructurada triangular.

• Integración de las ecuaciones de Saint-Venant mediante volúmenes finitos asumiendo:

o Flujo predominantemente horizontal.

o Variación vertical de la velocidad despreciable.

El proceso de integración es un proceso iterativo que parte de una situación inicial (n) de partida.

• Se integra para la situación n+1, a lo largo de toda la geometría.

• Se comprueba que la precisión alcanzada es suficiente.

• En caso contrario, se itera incluso variando el .

• Si la precisión es suficiente, se inicia un nuevo proceso iterativo para el siguiente momento temporal.






17

FIGURA 14: Esquema de ecuaciones 2D resueltas por InfoWorks ICM

La estabilidad numérica del modelo depende entre otros de los siguientes factores:

• Situación inicial.

• .

• Gradientes de velocidad, calado y caudal.

InfoWorks almacena todos los datos hidráulicos del modelo en una base de datos maestra. Como datos de entrada, esta base requiere al menos de los siguientes:

• Modelo digital del terreno (MDT) en alguno de los siguientes formatos:

o TIN (Triangulated Irregular Network o Red Irregular de Triángulos).

o GRID (Red regular tipo rejilla).

• Definición de la zona de mallado sobre el MDT anterior (contorno 2D).

• Condiciones de contorno a la entrada y a la salida.

• Archivo de caudales a simular.

• Archivo de simulación con los parámetros temporales.

El listado anterior constituye una condición de mínimos que se puede ampliar opcionalmente con alguno de los siguientes elementos, si el modelo lo precisa:

• Condición inicial del contorno de simulación.

• Puentes (se pueden modelizar como elementos unidimensionales conectados a la malla 2D o como elementos 2D puros).

• Infraestructuras (azudes, compuertas, etc.).






18

• Zonas de mallado específicas que precisen un mayor o menor detalle.

• Zonas de rugosidad.

La metodología general que se sigue a la hora de crear un modelo bidimensional en InfoWorks ICM se puede sintetizar con el siguiente esquema de trabajo:

• Crear las carpetas de trabajo y carpetas de resultados.

• Generar la base de datos maestra que contendrá la información de trabajo y los resultados. Este fichero puede ser intercambiado con otros usuarios o archivado y contiene todos los datos asociados a la simulación.

• El MDT se importa a dicha base de datos en formato ASCII.

FIGURA 15: Mallado sobre ortofoto (arriba) y sobre MDT (abajo) con InfoWorks ICM

• Definición de la zona 2D donde se realizarán los cálculos hidrodinámicos bidimensionales. Esta zona se crea directamente en la interfaz del programa mediante la definición manual de un polígono, o bien se importa desde un CAD o shape entre otros formatos disponibles. Los parámetros de discretización de la






19

malla se definen en esta capa, donde se aplica para todos los modelos la opción de mallado sensible al terreno1.

• La distribución de rugosidades generada con el procedimiento descrito anteriormente a partir de los datos del SIOSE se importa en formato shape. Este archivo debe coincidir exactamente (es recomendable la aplicación de reglas topológicas) con la extensión de la zona 2D anterior, o bien cubrir un área mayor que dicho polígono.

• Se incorporan a la red del modelo las líneas de rotura (elemento tipo “línea general”) teniendo en cuenta la morfología del terreno y prestando especial atención a:

o El pie de talud y la coronación en el cauce principal.

o Los escarpes en las llanuras de inundación.

o Los terraplenes de infraestructuras lineales.

o Muros de encauzamiento.

• Antes de ejecutar el proceso de mallado, se incorporan las estructuras. Para crear una estructura tipo puente dentro del mallado 2D se deben definir dos elementos:

o Tablero del puente: (Estructura Base Lineal 2D - Base Linear Structure 2D).

o Vanos del puente: (Estructura Lineal de Puente 2D - Bridge Linear Structure 2D), en el que la forma del vano se puede definir con tantos puntos como se desee.

• Se definen las líneas (o nodos en su caso) sobre las que se introducen las condiciones de entrada y salida del tramo analizado. La condición de entrada aplicada en todos los casos es una línea donde se introducen los datos de caudal correspondientes, mientras que la de salida es siempre la opción de Normal Depth o condición de calado normal (este calado se calcula considerando que la pendiente final del tramo se prolonga indefinidamente).

1 El algoritmo de mallado de InfoWorks ICM permite definir un error cordal y un ángulo mínimo objetivo. Igualmente, el algoritmo de triangulación se apoya en las líneas de rotura. Se entiende por error cordal la máxima altura admisible en los vértices de un elemento.






20

FIGURA 16: Ejemplo de introducción de puentes en InfoWorks ICM

• En los tramos que confluyen en el mar, se impone además una condición de marea y oleaje en la modelización. Estos datos proceden del estudio realizado por el CEDEX en el marco de los trabajos de desarrollo del Sistema Nacional de Cartografía de Zonas Inundables.

• Finalmente se ejecuta el mallado del modelo con la opción de mallado sensible y los parámetros de discretización mencionados anteriormente.

• Con la geometría completamente definida, ya se pueden introducir los datos de caudal en las condiciones de entrada al modelo, definir los parámetros temporales de la simulación (paso temporal, tiempo de simulación, etc.) y los archivos que la conforman (red del modelo, archivos de caudal o de condiciones iniciales, etc.) y ejecutar el modelo para proceder al análisis de sensibilidad y calibración previos a la validación del mismo.

3.2.4. ESTIMACIÓN INICIAL DE LAS PÉRDIDAS DE CARGA

Siguiendo las recomendaciones de la Guía Metodológica para el desarrollo del Sistema Nacional de Zonas Inundables (en adelante Guía Metodológica), editada en el año 2011 por el Ministerio de Agricultura, Alimentación y Medio Ambiente, para la estimación de las pérdidas de carga por fricción se siguen dos procedimientos diferenciados con aplicación en dos zonas también diferenciadas del medio fluvial: el cauce y la llanura de inundación.

La rugosidad en la llanura de inundación depende fundamentalmente de los usos del suelo en la llanura de inundación. Existen numerosos estudios que estiman el coeficiente de rugosidad en función de la actividad desarrollada sobre la zona inundable y estas estimaciones están referidas normalmente al número de Manning y, por lo general, toman valores que pueden variar dentro de unos rangos.

Como información de partida sobre los usos del suelo, se ha empleado el SIOSE aplicando la siguiente equivalencia a las coberturas simples que permite ponderar en cada polígono el resultado final en caso de ser éste compuesto:






21

Valores de Manning asignados a las coberturas del SIOSE

Cobertura simple Etiqueta Manning adoptado Arroz CHA 0,040

Cultivos herbáceos distintos de arroz CHL 0,040 Cítricos LFC 0,060

Frutales no cítricos LFN 0,060 Viñedo LVI 0,050 Olivar LOL 0,060

Otros cultivos leñosos LOC 0,040 Prados PRD 0,045

Pastizales PST 0,045 Perennifolias FDP 0,070 Caducifolias FDC 0,070 Coníferas CNF 0,070 Matorral MTR 0,060

Playas, dunas y arenales PDA 0,025 Ramblas RMB 0,035 Roquedo RQD No existe

Acantilados marinos ACM 0,025 Afloramientos rocosos y roquedos ARR 0,035

Coladas lávicas cuaternarias CLC 0,030 Canchales CCH 0,035

Suelo desnudo SDN 0,025 Zonas quemadas ZQM 0,040

Glaciares y nieves permanentes GNP 0,025

Edificación EDF 0,04 para modelización con

edificios y 0,12 para modelización sin edificios

Zona verde artificial y arbolado urbano ZAU 0.090

Lámina de agua artificial LAA 0,015 Vial, aparcamiento, zona peatonal sin

vegetación VAP 0,040

Suelo no edificado SNE 0,040 Otras construcciones OCT 0,040

Zonas de extracción o vertido ZEV 0,040 Zonas pantanosas HPA 0,040

Turberas HTU 0,040 Salinas continentales HSA 0,040

Marismas HMA 0,040 Salinas marinas HSM 0,040 Cursos de agua ACU 0,040 Lagos y lagunas ALG 0,025

Embalses AEM 0,025 Lagunas costeras ALC 0,030

Estuarios AES 0,030 Mares y océanos AMO 0,025

En la zona del cauce, para la obtención de unos valores lo más aproximados a la realidad, a partir de los datos recopilados en las visitas de campo y de información






22

gráfica disponible, se caracterizarán los distintos tramos a partir de su naturaleza y propiedades morfológicas. En concreto, se ha utilizado la clasificación de tramos fluviales propuesta por Cowan, que permite, mediante una sencilla formulación, obtener el número de rugosidad de Manning. De esta forma:

n = (n0 + n1 + n2 + n3 + n4).m5

donde:

• n0 depende del material que conforma el lecho.

• n1 depende del grado de irregularidad del lecho.

• n2 depende del tipo de variación de la sección transversal.

• n3 depende del efecto relativo de obstrucciones.

• n4 depende de tipo de vegetación existente.

• m5 es función de la cantidad de meandros.

Valores de los parámetros de la formulación de Cowan

Material n0 Hormigón/Asfalto 0,015

Tierra 0,020 Roca cortada 0,025

Grava fina 0,024 Grava gruesa 0,028

Grado de irregularidad n1 Bajo 0,000

Menor 0,005 Medio 0,007

Moderado 0,010 Alto 0,020

Variaciones de sección transversal n2 Gradual 0,000

Ocasionalmente alternante 0,005 Frecuentemente alternante 0,010

Muy frecuentemente alternante 0,015 Efecto de las obstrucciones n3

Pequeño <5% 0,000 Menor 5-10% 0,010

Menor 10-15% 0,015 Apreciable 15-25% 0,020 Apreciable 25-40% 0,025 Apreciable 40-50% 0,030

Alto 50-75% 0,050






23

Valores de los parámetros de la formulación de Cowan

Alto >75% 0,060 Vegetación n4

Nula 0,000 Baja 0,008

Media 0,017 Alta 0,034

Muy alta 0,075 Cantidad de meandros m5

Menor 2 1,000 Apreciable 1,150

Alta 2 1,300

En el caso de modelos bidimensionales, los valores de n2 y m5 adoptan siempre un valor de 0 y 1 respectivamente.

Una vez definidos los coeficientes de rugosidad del cauce y de las llanuras de inundación, se integran en una única capa shape que sirve directamente como dato de entrada para el modelo hidráulico.

FIGURA 17: Ejemplo de distribución espacial del coeficiente de rugosidad de Manning a partir de la

cartografía de usos del suelo del SIOSE y su revisión con la ortofoto PNOA de Máxima Actualidad

3.2.5. CALIBRACIÓN Y VALIDACIÓN DEL MODELO.

Una vez confeccionado el modelo hidráulico, resulta imprescindible realizar un proceso de calibración del mismo, consistente básicamente en ajustar la superficie de inundación que obtiene el modelo hidráulico con las producidas en la realidad. En el proceso de calibración se determinarán los valores asignados a los coeficientes de rugosidad, mejorando las estimaciones realizadas inicialmente.






24

En aquellos casos en los que se disponga de información de detalle de inundaciones históricas se llevará a cabo una calibración del modelo bastante precisa. En este sentido resulta de gran ayuda disponer de fotografías aéreas de la inundación en las que se pueda apreciar la extensión de la misma. Se trata de una información de amplia extensión, por lo que no se ve influida por efectos locales.

En su defecto, se debe recopilar información fotográfica de la extensión de la inundación que han alcanzado las avenidas recientes, que dadas las condiciones de contorno de la zona de estudio se pueden dar como válidas al no apreciarse por efectos locales, como obstrucciones del cauce, rotura de motas, etc.

FIGURA 18: Fotografía panorámica de la avenida de 2016 en un sector del Ebro localizado aguas arriba

de la zona de estudio. 28 de febrero de 2016. Hora. 10.30 a.m.

3.2.6. RESULTADOS OBTENIDOS.

Se ha realizado la simulación hidráulica del tramo de estudio, utilizando como caudal circulante 913 m3/s, caudal correspondiente al adoptado para la Máxima Crecida Ordinaria. El área inundada se recoge en la figura 19.

Adicionalmente se ha simulado un caudal de 1.200 m3/s, máximo caudal esperable correspondiente al obtenido del programa CAUMAX para T=2,5 años. Este caso se muestra en la figura 20, pudiéndose comprobar que, con dicho volumen de agua, el desbordamiento es generalizado, especialmente por la margen derecha.

En el apartado siguiente se adjuntan las láminas de detalle de la superficie afectada por la Máxima Crecida Ordinaria en el tramo de estudio.






25

FIGURA 19: Resultados de la inundación obtenida para un caudal de 913 m3/s, correspondiente al caudal

adoptado para la MCO.

FIGURA 20: Resultados de la inundación obtenida para un caudal de 1200 m3/s, correspondiente al

caudal obtenido por el programa CAUMAX para T=2,5 años






26

4. LAMINAS

Localización

Delimitación del D.P.H.. Criterios Hidráulicos




ANEXO Nº 4. ESTUDIO HISTÓRICO, GEOMORFOLÓGICO Y ECOLÓGICO


ANEXO Nº4.

ESTUDIO HISTÓRICO, GEOMORFOLÓGICO Y ECOLÓGICO






ÍNDICE

1. ESTUDIOS HISTÓRICO, GEOMORFOLÓGICO Y MEDIOAMBIENTAL. ........ 1

2. ANÁLISIS HISTÓRICO. ................................................................................. 1

3. ANÁLISIS GEOMORFOLÓGICO ................................................................. 4

4. ANÁLISIS AMBIENTAL ................................................................................. 6

5. LAMINAS ................................................................................................... 11






1

1. ESTUDIOS HISTÓRICO, GEOMORFOLÓGICO Y ECOLÓGICO.

Dentro de estos estudios, un hito fundamental para el registro histórico de la evolución de los cauces fluviales lo constituyen la realización de los denominados vuelos americanos A (1945/46) y B (1956/57), en especial el segundo por la mayor calidad de los fotogramas aéreos y su disponibilidad. Sin desdeñar el vuelo americano A, de georreferenciación menos precisa y no siempre disponible, los fotogramas del vuelo de 1956/57 permiten delimitar un cauce histórico suficientemente antiguo con precisión y definir su evolución hasta su estado actual mediante comparación con fotografías aéreas de fechas intermedias, identificar los cambios habidos y localizar las zonas más activas e inundables del medio fluvial.

El trabajo de recopilación de esta información gráfica se encuentra disponible en los visores y servicios WMS de esta Confederación (SITEbro - Sistema de Información Territorial de la cuenca del Ebro), del Ministerio de Transición Ecológica y el Reto Demográfico, del IGN y de las Comunidades Autónomas.

Debe indicarse que, conforme al artículo 4 del Reglamento de Dominio Público Hidráulico, las referencias históricas son uno de los medios de determinación del terreno cubierto por las aguas en las máximas crecidas ordinarias (cauce de río).

El análisis geomorfológico, necesario para delimitar cada uno de los elementos a definir cartográficamente se ha de basar en trabajos de campo, análisis sobre los resultados del vuelo LiDAR o en su defecto, fotografías aéreas con visión estereoscópica, ortofotos y las topografías más detalladas disponibles. Los objetivos de este análisis son:

a) Identificación de zonas inundables en llanuras aluviales y en zonas susceptibles de sufrir avenidas.

b) Delimitación de la máxima crecida ordinaria. Se analiza el trazado del cauce y su evolución mediante fotografías históricas (al menos, desde el año 1957 o, si no estuviera disponible, la más próxima a esta fecha). El análisis combinado de esta información mostrará la delimitación del Dominio Público Hidráulico según criterios geomorfológicos.

c) Zonificación del área inundable.

El análisis ambiental o ecológico se basa en trabajos de campo y en ortofotos históricas para determinar las zonas de vegetación de ribera y sotos y caracterizar, si los hubiere, los sectores afectados de las zonas incluidas dentro de la red NATURA 2000, que es una red que incluye las áreas de conservación de la biodiversidad en la Unión Europea. Consta de Zonas Especiales de Conservación (ZEC), designadas de acuerdo con la Directiva Hábitat, así como de Zonas de Especial Protección para las Aves (ZEPA) establecidas en virtud de la Directiva Aves y Lugares de Interés Comunitario (LIC). Su finalidad es asegurar la supervivencia a largo plazo de las especies y los hábitats más amenazados de Europa, contribuyendo a detener la pérdida de biodiversidad ocasionada por el impacto adverso de las actividades humanas. Es el principal instrumento para la conservación de la naturaleza en la Unión Europea.

2. ANÁLISIS HISTÓRICO.

La base de datos disponible para evaluar los episodios de inundación en la cuenca del Ebro contiene más de 1330 eventos de este tipo. Desgraciadamente, aunque extensa, dicha relación no puede considerarse exhaustiva, pues muchos episodios, especialmente los antiguos, no han






2

dejado referencias explícitas en archivos o registros por haber afectado a áreas poco pobladas, acontecer en tiempos de los que apenas se conservan reseñas escritas o por no haber causados daños significativos.

De los episodios registrados, muchos de ellos, en especial los anteriores al siglo XX, hacen referencia general al río Ebro como causante de la inundación, pero sin definir áreas concretas o citando sólo los principales núcleos urbanos afectados, tales como Zaragoza, Logroño, Tudela, etc. En muchos de estos casos, es difícil imaginar una afección simultánea a las localidades de, por ejemplo, Logroño y Zaragoza, cuya ubicación les hace menos propensa a sufrir la acción de las aguas que otras poblaciones cercanas, y que no haya constancia de daños en el tramo intermedio. En estas situaciones, lo lógico es pensar que gran parte de la ribera del río entre ambas localidades se ha visto afectada por el episodio de inundación, pero los daños causados en esas zonas no tienen reflejo en las fuentes bibliográficas al resultar menos llamativos que los de la gran población: para magnificar el evento ante la autoridad, tiene mayor impacto la reseña de daños o la pérdida de algún edificio o infraestructura en Logroño que la de una parte de la cosecha en Arrúbal.

Con esta premisa, se ha de colegir que la ribera del Ebro en Arrúbal ha debido sufrir un número considerable de eventos de inundación, pero de los que no se tiene constancia por falta de registros específicos. En tiempos modernos sigue pasando algo parecido, ya que sólo hay una referencia concreta de inundación en Arrúbal, que data del año 2003 y se trata de una avenida de muy bajo periodo de retorno, inferior a 5 años: los caudales registrados en la cercana estación de aforos de Mendavia fueron ligeramente superiores a los de la máxima crecida ordinaria que sirve para definir el DPH hidráulico. En este caso, además se dispone de una estimación del área inundada en el paraje de El Soto (Figura 1), sensiblemente concordante con el área objeto de deslinde.

Figura 1. Estimación a partir de fotointerpretación del área inundada en el paraje de El Soto, en Arrúbal,

por la avenida del Ebro de febrero de 2003. FUENTE: SITEbro

Además de esta referencia, el inventario antes citado recoge otros dos eventos del Ebro acontecidos en diciembre de 1993 que afectaron al vecino término municipal de Mendavia (Navarra) y que, lógicamente, también debieron repercutir en las riberas del río en Arrúbal pero de los que no hay reseña explícita, probablemente por falta de daños significativos.






3

Por otro lado, en las imágenes aéreas disponibles correspondientes al denominado Vuelo Americano A, de 1945/46, se pueden observar daños en campos de cultivo sólo achacables a una avenida acontecida poco antes (meses o algún año) de su toma y de la que no hay referencia.

FIGURA 2: Evolución del río Ebro en el paraje de El Soto en Arrúbal desde 1945 a la actualidad

En cualquier caso, el seguimiento de las imágenes aéreas disponibles pone de manifiesto la evolución del Ebro, consistente en la consolidación del cauce de aguas bajas histórico como cauce ordinario actual y la transformación de sus brazos secundarios activos en parte del cauce de aguas altas del río; es decir, se ha simplificado su trazado reduciéndose el número de brazo activos, lo que supone caudales circulantes menores y, por tanto, una reducción aparente del área afectada por la que debería ser la máxima crecida ordinaria.






4

3. ANÁLISIS GEOMORFOLÓGICO

Las modificaciones geomorfológicas ligadas a la actividad fluvial tienden a concentrarse en las inmediaciones del cauce, en las zonas afectadas por las crecidas ordinarias y por avenidas de bajo periodo de retorno. Sin embargo, en zonas más alejadas del cauce, en especial las situadas a cotas sensiblemente superiores al lecho, su repercusión es mucho menor, incluso despreciable. Es decir, en avenidas y crecidas, las afecciones por el agua incluyen la modificación de los rasgos morfológicos preexistentes en el entorno del cauce en cuestión, mientras que en zonas alejadas la lámina de agua tiende a adaptarse a los rasgos existentes sin modificarlos en grado significativo.

Estas modificaciones tienen lugar mientras los ríos mantengan un cierto grado de naturalidad, ya que son entidades dinámicas que evolucionan en el tiempo modificando sus trazados y principales rasgos morfológicos definitorios. Por este motivo y bajo dichas circunstancias, tanto el río como sus riberas son cambiantes, resultando imposible definir, a partir de una imagen tomada en un momento dado, límites perpetuos para los rasgos morfoestructurales que delimitan los distintos espacios fluviales.

Como consecuencia de ello, el análisis de los rasgos geomorfológicos del entorno fluvial permite definir con precisión adecuada las zonas de inundabilidad elevada (las que se inundan frecuentemente por el desbordamiento de los ríos) disminuyendo la fiabilidad del criterio geomorfológico a medida que su grado de inundabilidad se reduce. Por este motivo, el análisis geomorfológico se ha centrado fundamentalmente en definir el dominio correspondiente al cauce, entendiendo como tal lo recogido en el artículo 4 del Texto Refundido de la Ley de Aguas: “el terreno cubierto por las aguas en las máximas crecidas ordinarias”. Además, el Texto Refundido de la Ley de Aguas, en su artículo 6, diferencia dentro del cauce la zona de ribera, definiendo las riberas como “las fajas laterales de los cauces situadas por encima del nivel de aguas bajas”.

Por todo ello se han efectuado dos series de cartografías geomorfológicas, cuyas características básicas son las siguientes:

• Cartografía geomorfológica del cauce en su situación actual, que refleja el estado y rasgos del cauce bajo condiciones alteradas, derivadas de las obras de regulación y derivación de caudales y de los tratamientos de márgenes y encauzamientos efectuados en las últimas décadas. El análisis se ha realizado fundamentalmente a partir de la interpretación de las ortoimágenes del PNOA en combinación con los Modelos Digitales del Terreno (derivados de los registros LIDAR tratados que han sido utilizados en la modelación hidráulica) y el apoyo de imágenes estereoscópicas. Se ha digitalizado esta cartografía directamente sobre las ortoimágenes del PNOA.

• Cartografía geomorfológica de los cauces a mediados del siglo XX, correspondiente a una situación seminatural o natural. Se ha efectuado a partir del análisis de los fotogramas aéreos de la época georreferenciados a tal efecto. Como consecuencia de las deformaciones resultantes a pesar de la ortorrectificación de los fotogramas, se han tenido que realizar pequeños ajustes finales manuales.

Para la definición cartográfica del cauce en el sentido recogido por el Texto Refundido de la Ley de Aguas se estudia la situación actual alterada y la histórica a partir de imágenes aéreas. En este caso concreto, las imágenes empleadas en estas cartografías han sido, para la situación actual, la ortoimagen del PNOA de máxima actualidad y, para la situación histórica, los fotogramas aéreos del denominado Vuelo Americano A (1945-46), el más antiguo existente con calidad suficiente para efectuar con detalle este tipo de trabajo. La comparación de los resultados obtenidos en ambos casos permite comprobar el dinamismo del medio fluvial y, con ello, la superficie cambiante identificable como DPH en el paraje de El Soto en Arrúbal.






5

La evolución geomorfológica de este paraje se recoge en las figuras 3 y 4, y en ellas se constata, como ya se ha manifestado en el apartado precedente, una simplificación del trazado del río y una reducción aparente de la superficie afectada por la máxima crecida ordinaria, crecida que marca el límite teórico del DPH.

FIGURA 3: Análisis geomorfológico del paraje de El Soto en Arrúbal a partir de las imágenes aéreas de

Vuelo Americano A, de 1945-46. Leyenda.1: Terrenos no inundables. 2: Terrenos inundables por avenidas extraordinarias asociadas a

diferentes periodos de retorno. 3. Terrenos inundables afectados por una avenida reciente de bajo periodo de retorno. 4: Terrenos habitualmente inundables adscribibles a la máxima crecida ordinaria. 5: Canales secundarios funcionales en situación ordinaria. 6: Canal de aguas bajas (estiaje) siempre inundado. Líneas negras: Canales secundarios funcionales sólo en situación de crecida.

Esta reducción aparente sugiere una disminución en los caudales circulantes, lo que constituye una circunstancia obvia desde el momento en que se construyeron los grandes embalses regulatorios localizados en la cuenca del Ebro aguas arriba del paraje objeto de la solicitud de deslinde (embalses del Ebro, Ullibarri-Gamboa, Urrunaga, etc.). Estos embalses contribuyen a la regulación del caudal del río y permiten una cierta laminación de sus avenidas. Es decir, con su construcción, el régimen del Ebro ha pasado de natural, cuya máxima crecida ordinaria marcaría el límite del DPH hidráulico, a alterado, situación en la que ese caudal circulando correspondería a una situación de avenida para el régimen natural.

Todo ello ha facilitado, además de la simplificación ya mencionada del cauce del Ebro, la progresión de las áreas cultivadas en detrimento de los terrenos inundables asociados a la máxima crecida ordinaria.






6

FIGURA 4: Análisis geomorfológico del paraje de El Soto en Arrúbal a partir de la ortoimagen del PNOA de

máxima actualidad. Leyenda.1: Terrenos no inundables. 2: Terrenos inundables por avenidas extraordinarias asociadas a

diferentes periodos de retorno. 3. Terrenos inundables afectados por avenidas de muy bajo periodo de retorno. 4: Terrenos habitualmente inundables adscribibles a la máxima crecida ordinaria, incluyendo cauces abandonados habitualmente no funcionales o galachos (G). 5: Cauce ordinario, incluyendo canales secundarios funcionales en situación de crecida. 6: Canal de aguas bajas (estiaje) siempre inundado.

4. ANÁLISIS ECOLÓGICO

En esta zona, el río Ebro conforma el límite norte de La Rioja con la Comunidad Foral de Navarra y por su longitud y caudal, así como por la abundancia y singularidad de las especies y hábitats que alberga, constituye el curso fluvial más importante tanto de La Rioja como, en este caso junto al sistema Aragón-Arga, de Navarra. Esta importancia queda de manifiesto por la existencia en sus márgenes o proximidades de importantes áreas sujetas a protección ambiental: LIC’s ES2200040 “Río Ebro” y ES2200031 “Yesos de la Ribera Estellesa”, en Navarra, y LIC ES2300008 "Sotos y Riberas del Ebro”, en La Rioja.

No obstante, el paraje de El Soto, donde se ubican las parcelas objeto de solicitud de deslinde, no está incluido en ninguna de estas áreas protegidas (FIGURA 5). A pesar de ello, la potencialidad ambiental del paraje a largo plazo puede ser equivalente a la de dichas zonas protegidas, pudiendo constituirse, en caso de recuperación de su población, en una potencial área de expansión del visón europeo (Mustela lutreola), especie declarada “En Peligro Crítico” y que es el mamífero en mayor peligro de extinción de Europa por la competencia que le supone el invasor y más agresivo visón americano (Neovison vison).






7

FIGURA 5: Relación espacial entre el área de estudio y los espacios protegidos más cercanos

El paraje de El Soto de Arrúbal queda incluido, de acuerdo con el Atlas Faunístico de España, en la cuadrícula de 10x10 km denominada 30TWM69. De acuerdo con dicho Atlas, en la cuadrícula mencionada se han inventariado un total de 170 especies de vertebrados, desglosados en 24 especies de mamíferos, 119 de aves, 10 de reptiles, 6 de anfibios y 11 de peces. La relación de estas especies, la mayoría de las cuales son susceptibles de localizarse en un momento u otro en la zona de estudio, se recoge en el cuadro siguiente:

ESPECIE ORDEN CLASE Rattus norvegicus Rata parda o común Roedores Mamíferos

Vulpes vulpes Zorro Carnívoros Mamíferos Crocidura russula Musaraña gris Soricomorfos Mamíferos

Mus musculus Ratón doméstico Roedores Mamíferos Lutra lutra Nutria Carnívoros Mamíferos

Arvicola sapidus Rata de agua Roedores Mamíferos Pipistrellus kuhlii Murciélago de bosque claro Quirópteros Mamíferos

Mustela nivalis Comadreja Carnívoros Mamíferos Pipistrellus nathusi Murciélago enano de bosque Quirópteros Mamíferos

Mustela putorius Turón Carnívoros Mamíferos Erinaceus europaeus Erizo europeo Erinaceomorfos Mamíferos Pipistrellus pygmaeus Murciélago de Cabrera Quirópteros Mamíferos

Lepus granatensis Liebre ibérica Lagomorfos Mamíferos Genetta genetta Gineta Carnívoros Mamíferos

Martes foina Garduña Carnívoros Mamíferos Miniopterus schreibersii Murciélago de cueva Quirópteros Mamíferos

Oryctolagus cuniculus Conejo Lagomorfos Mamíferos Rhinolophus euryale Murciélago de herradura Quirópteros Mamíferos

Talpa europaea Topo europeo Erinaceomorfos Mamíferos Pipistrellus pipistrellus Murciélago común o enano Quirópteros Mamíferos

Mustela lutreola Visón europeo Carnívoros Mamíferos






8

ESPECIE ORDEN CLASE Sus scrofa Jabalí Artiodáctilos Mamíferos

Felis silvestris Gato montés Carnívoros Mamíferos Meles meles Tejón Carnívoros Mamíferos

Remiz pendulinus Pájaro moscón Paseriformes Aves Porzana pusilla Polluela chica Galliformes Aves

Apus melba Vencejo real Apodiformes Aves Hippolais pallida Zarcero pálido Paseriformes Aves

Ardea cinerea Garza real Ciconiformes Aves Accipiter nisus Gavilán Falconiformes Aves

Falco tinnunculus Cernícalo vulgar Falconiformes Aves Pterocles orientalis Ortega Columbiformes Aves

Sylvia hortensis Curruca mirlona Paseriformes Aves Ptyonoprogne rupestris Avión roquero Paseriformes Aves

Alectoris rufa Perdiz roja Galliformes Aves Columba livia/domestica Paloma bravía Columbiformes Aves

Phoenicurus ochruros Colirrojo tizón Paseriformes Aves Emberiza calandra Triguero Paseriformes Aves

Motacilla cinerea Lavandera cascadeña Paseriformes Aves Buteo buteo Ratonero común Falconiformes Aves

Columba oenas Paloma zurita Columbiformes Aves Gyps fulvus Buitre leonado Falconiformes Aves

Cisticola juncidis Buitrón Paseriformes Aves Burhinus oedicnemus Alcaraván Charadriformes Aves

Passer montanus Gorrión molinero Paseriformes Aves Otus scops Autillo Strigiformes Aves

Galerida theklae Cogujada montesina Paseriformes Aves Falco subbuteo Alcotán Falconiformes Aves Ardea purpurea Garza imperial Ciconiformes Aves

Certhia brachydactyla Agateador común Paseriformes Aves Nycticorax nycticorax Martinete Ciconiformes Aves

Cuculus canorus Cuco Cuculiformes Aves Turdus philomelos Zorzal común Paseriformes Aves

Fringilla coelebs Pinzón vulgar Paseriformes Aves Sturnus unicolor Estornino negro Paseriformes Aves

Carduelis cannabina Pardillo común Paseriformes Aves Hippolais polyglotta Zarcero común Paseriformes Aves

Phylloscopus bonelli Mosquitero papialbo Paseriformes Aves Lanius excubitor Alcaudón real norteño Paseriformes Aves Falco peregrinus Halcón peregrino Falconiformes Aves

Motacilla alba Lavandera blanca Paseriformes Aves Pica pica Urraca Paseriformes Aves

Anthus campestris Bisbita campestre Paseriformes Aves Columba palumbus Paloma torcaz Columbiformes Aves

Sylvia undata Curruca rabilarga Paseriformes Aves Sylvia cantillans Curruca carrasqueña Paseriformes Aves

Cettia cetti Ruiseñor bastardo Paseriformes Aves Emberiza cirlus Escribano soteño Paseriformes Aves

Streptopelia decaocto Tórtola turca Columbiformes Aves Parus caeruleus Herrerillo común Paseriformes Aves

Parus major Carbonero común Paseriformes Aves Delichon urbica Avión común Paseriformes Aves

Streptopelia turtur Tórtola común Columbiformes Aves Tachybaptus ruficollis Zampullín chico Gaviiformes Aves

Acrocephalus scirpaceus Carricero común Paseriformes Aves Passer domesticus Gorrión común Paseriformes Aves Charadrius dubius Chorlitejo chico Charadriformes Aves

Acrocephalus arundinaceus Carricero tordal Paseriformes Aves






9

ESPECIE ORDEN CLASE Sylvia atricapilla Curruca capirotada Paseriformes Aves

Circaetus gallicus Águila culebrera Falconiformes Aves Tyto alba Lechuza común Strigiformes Aves

Apus apus Vencejo común Apodiformes Aves Corvus corax Cuervo Paseriformes Aves

Corvus monedula Grajilla Paseriformes Aves Jynx torquilla Torcecuello Piciformes Aves

Clamator glandarius Críalo Cuculiformes Aves Monticola solitarius Roquero solitario Paseriformes Aves Dendrocopos minor Pico menor Piciformes Aves Carduelis carduelis Jilguero Paseriformes Aves

Carduelis chloris Verderón común Paseriformes Aves Oriolus oriolus Oropéndola Paseriformes Aves

Bubo bubo Búho real Strigiformes Aves Sylvia conspicillata Curruca tomillera Paseriformes Aves

Ciconia ciconia Cigüeña común Ciconiformes Aves Pyrrhocorax pyrrhocorax Chova piquirroja Paseriformes Aves

Corvus corone Corneja Paseriformes Aves Sylvia borin Curruca mosquitera Paseriformes Aves

Phylloscopus collybita/ibericus Mosquitero común Paseriformes Aves Hirundo rustica Golondrina común Paseriformes Aves Turdus merula Mirlo común Paseriformes Aves

Sturnus vulgaris Estornino pinto Paseriformes Aves Vanellus vanellus Avefría Charadriformes Aves

Asio otus Búho chico Strigiformes Aves Upupa epops Abubilla Coraciformes Aves

Gallinula chloropus Polla de agua Galliformes Aves Saxicola torquata Tarabilla común Paseriformes Aves

Picus viridis Pito real Piciformes Aves Luscinia megarhynchos Ruiseñor común Paseriformes Aves

Hieraaetus pennatus Águila calzada Falconiformes Aves Alauda arvensis Alondra común Paseriformes Aves Rallus aquaticus Rascón Galliformes Aves

Fulica atra Focha común Galliformes Aves Motacilla flava Lavandera boyera Paseriformes Aves

Aegithalos caudatus Mito Paseriformes Aves Coturnix coturnix Codorniz Galliformes Aves

Riparia riparia Avión zapador Paseriformes Aves Falco naumanni Cernícalo primilla Falconiformes Aves

Actitis hypoleucos Andarríos chico Charadriformes Aves Phylloscopus brehmii (ibericus) Mosquitero iberico Paseriformes Aves

Anas platyrhynchos Ánade real Anseriformes Aves Athene noctua Mochuelo Strigiformes Aves

Alcedo atthis Martín pescador Coraciformes Aves Caprimulgus europaeus Chotacabras gris Piciformes Aves

Petronia petronia Gorrión chillón Paseriformes Aves Emberiza hortulana Escribano hortelano Paseriformes Aves

Troglodytes troglodytes Chochín Paseriformes Aves Lanius collurio Alcaudón dorsirrojo Paseriformes Aves

Muscicapa striata Papamoscas gris Paseriformes Aves Calandrella brachydactyla Terrera común Paseriformes Aves

Serinus serinus Verdecillo Paseriformes Aves Merops apiaster Abejaruco común Coraciformes Aves

Oenanthe oenanthe Collalba gris Paseriformes Aves Lanius senator Alcaudón común Paseriformes Aves

Galerida cristata Cogujada común Paseriformes Aves Melanocorypha calandra Calandria Paseriformes Aves






10

ESPECIE ORDEN CLASE Sylvia melanocephala Curruca cabecinegra Paseriformes Aves

Oenanthe hispanica Collalba rubia Paseriformes Aves Circus aeruginosus Aguilucho lagunero Falconiformes Aves

Milvus migrans Milano negro Falconiformes Aves Circus cyaneus Aguilucho pálido Falconiformes Aves

Circus pygargus Aguilucho cenizo Falconiformes Aves Hieraaetus fasciatus Águila perdicera Falconiformes Aves

Emberiza schoeniclus Escribano palustre Paseriformes Aves Neophron percnopterus Alimoche Falconiformes Aves

Malpolon monspessulanus Culebra bastarda Squamata Reptiles Natrix maura Culebra viperina Squamata Reptiles

Anguis fragilis Lución Squamata Reptiles Chalcides striatus Eslizón tridáctilo ibérico Squamata Reptiles

Podarcis hispanica Lagartija ibérica Squamata Reptiles Lacerta lepida (Timon lepidus) Lagarto ocelado Squamata Reptiles

Rhinechis (Elaphe) scalaris Culebra de escalera Squamata Reptiles Natrix natrix Culebra de collar Squamata Reptiles

Psammodromus hispanicus Lagartija cenicienta Squamata Reptiles Coronella girondica Culebra lisa meridional Squamata Reptiles

Rana perezi Rana común Anura Anfibios Pelodytes punctatus Sapillo moteado Anura Anfibios

Pelobates cultripes Sapo de espuelas Anura Anfibios Triturus marmoratus Tritón jaspeado Caudata Anfibios

Hyla arborea Ranita de San Antón Anura Anfibios Bufo calamita Sapo corredor Anura Anfibios Gobio lozanoi Gobio ibérico Ciprínidos Peces

Barbus graellsii Barbo de Graells Ciprínidos Peces Chondrostoma arcasii Bermejuela Ciprínidos Peces

Barbatula barbatula Lobo de río Ciprínidos Peces Chondrostoma miegii Madrilla Ciprínidos Peces

Cyprinus carpio Carpa común Ciprínidos Peces Barbus haasi Barbo colirrojo Ciprínidos Peces

Carassius auratus Carpa dorada Ciprínidos Peces Phoxinus phoxinus Piscardo Ciprínidos Peces

Cobitis calderoni Lamprehuela Ciprínidos Peces Salmo trutta Trucha común Salmónidos Peces

Los elementos faunísticos no constituyen factores claramente determinantes para la definición del límite del DPH, ya que, salvo la ictiofauna, asociada necesariamente a masas de agua, la presencia de ninguno de ellos es exclusiva de dicho ámbito: como consecuencia de su movilidad, ni siquiera la identificación continuada de individuos clasificables como propios de ecosistemas fluviales garantiza que ese contacto tenga lugar en dicho ámbito territorial.

Por el contrario, ciertas especies florísticas, las riparias, son características de las riberas fluviales, entendiendo éstas en el sentido que recoge el artículo 6 del Texto Refundido de la Ley de Aguas. Por dicho motivo, la presencia natural de estas especies riparias en el ámbito fluvial se considera suficiente para incluir las superficies ocupadas por ellas en el DPH (inclusión por motivo ecológico).

En el paraje de El Soto, en Arrúbal, se identifican manchas de vegetación riparia con diferentes estados de conservación, con sectores de alta calidad ambiental coexistiendo con otros peor conservados o, incluso, con vegetación riparia inexistente. Esto no es óbice para poder asignar al DPH las zonas ocupadas por ellas.

En general, las áreas ocupadas por este tipo de vegetación son compatibles con los límites definidos con criterios geomorfológicos, desbordando los mismos a lo largo del pie del escarpe de terraza que actúa de límite entre los terrenos inundables y los no inundables, donde progresan en






11

dirección al núcleo urbano de Arrúbal (FIGURA 6). Esta zona, esta vegetación no se puede ligar al DPH asociado al río Ebro, sino con un pequeño cauce secundario que desemboca en el Ebro en el límite sur de la zona de estudio.

FIGURA 6: Cartografía de la vegetación riparia asociable a DPH en el paraje de El Soto, en Arrúbal.

5. LAMINAS

Localización

Delimitación del dominio público hidráulico. Criterios geomorfológicos. Orografía

Delimitación del dominio público hidráulico. Criterios geomorfológicos. Vuelo Americano A. 1945-46.

Delimitación del dominio público hidráulico. Criterios geomorfológicos. Vuelo Americano B. 1956-57

Delimitación del dominio público hidráulico. Criterios geomorfológicos. Vuelo SIGPAC. 1997-2003

Delimitación del dominio público hidráulico. Criterios geomorfológicos. PNOA 2012

Delimitación del dominio público hidráulico. Criterios geomorfológicos. PNOA Máxima Actualidad. 2017

Delimitación del dominio público hidráulico. Criterios geomorfológicos. Cartografía sobre Vuelo Americano A. 1945-46.

Delimitación del dominio público hidráulico. Criterios geomorfológicos. Cartografía sobre PNOA Máxima Actualidad.






12

Delimitación del dominio público hidráulico. Criterios ecológicos. Vegetación riparia.




ANEXO Nº5. PLANOS DE LA PROPUESTA DE DESLINDE


ANEXO Nº5.

PLANOS DE LA PROPUESTA DE DESLINDE




ANEXO Nº5. PLANOS DE LA PROPUESTA DE DESLINDE


1. LAMINAS

Propuesta de Deslinde. Plano Guía Propuesta de deslinde. Escala 1:1000 Parcelas afectadas




DDDD

DD

DDDD

DDDD

DDDD

DDDD

DDDD

DDDD

DD

DDDD

DDDD

DD

DD

DD

DD

DD

DD

DD

DD

DD

DD

DDDD

DDDD

DDDD

DDDD

DDDD

DDDD

DDDD

DD

DD

DD

DDDD

DD

DD

DD

DD

DD

DD

DD

DD

DD

DD

DD

DD

DD

DD

DD

DD

DD

DD

DD

DD

DD

DD

DD

DD

DD

DD

DD

DDDD

DDDD

DDDD

DD

DD

DD

DD

DD

DD

DD

DDDD

DDDD

DDDD

DDDD

DD

DD

DDDD

DDDD

DDDD

DDDD

DDDD

DDDD

DDDD

DDDD

DDDD

DDDD

DDDD

DDDD

DDDD

DD

DD

DD

DD

DD

DD

DD

DD

DD

DDDD

DDDD

DDDD

DDDD

DDDD

DDDD

DDDD

DDDD

DDDD

DDDD

DDDD

DDDD

DDDD

DD

DD

DD

DD

DD

DD

DDDD

DDDD

DDDD

DDDD

DDDD

DDDD

DDDD

DD

DD

DDDD

DD

DD

DD

DD

DD

DD

DDDD

DDDD

DDDD

DDDD

DDDD

DDDDDD

DD

DD

DD

DD

DD

DD

DDDD

DD DD

DDDD

DDDD

DD

DD

DD

DD

DD

995996

1018

1019

1

5

4

3

2

Arrubal

ESCALAS: FECHA:

SEPTIEMBRE 20201:6.500 DESIGNACIÓN:EMPRESA CONSULTORA:

HOJA:

ORIGINAL A-3

Nº DE PLANO:

1.10 60 120Meters

1 de 1

PROPUESTA DE DESLINDE. ESCALA 1:1000.PLANO GUÍA.

¯

LEYENDA:DD DD Límites Municipales Cuadricula E 1:1000 Límite D.P.H Área Parcelas Afectadas Parcelario

TITULO DEL PROYECTO:

MEMORIA DESCRIPTIVA DEL PROCEDIMIENTO ADMINISTRATIVO DE DELIMITACIÓN DEL DOMINIO PÚBLICO HIDRÁULICO DE LAS

PARCELAS CATASTRALES Nº 995, 996, 1018 Y 1019 DEL POLÍGONO 1 DEL TÉRMINO MUNICIPAL DE ARRÚBAL (FINCA REGISTRAL Nº 2)




A

A

A

A

A

A

A

A

A

A

A

A

A

A

A

A

A

A

A

A

A

06

05

04

03

02

01

429

1016

594

1018

608

607

600

592

604

589

590

606593

603

588

602

601

591

611

605

609

612

610

575

599

576

587

20579

573

577

574

578

581

10579

582

583

584

10021002

596

ESCALAS: FECHA:


HOJA:

ORIGINAL A-3

Nº DE PLANO:

1.20 10 20Meters PROPUESTA DE DESLINDE. ESCALA 1:1000.

¯

LEYENDA:A Hitos Deslinde Límite D.P.H Área de Parcela Afectada Parcelario

1 de 5







A

A

A

A

A

A

A

A

A

A

A

A

A

A

A

A

A

A

A

A

A

06

429

618

612

584594

614

596

77

1018

615

619

78

611

1019

79

712

610

609

621 620

622

608

999

783

976

1001

1000

623

583

645

ESCALAS: FECHA:


HOJA:

ORIGINAL A-3

Nº DE PLANO:


¯


2 de 5







A

A

A

A

A

A

A

A

A

A

A

A

A

A

A

A

A

A

A

A

A

06

05

04

03

02

01

428

995

709

1017

710

708707

700

699704698

996

1027

703

706

705

696697

977

ESCALAS: FECHA:


HOJA:

ORIGINAL A-3

Nº DE PLANO:


¯


3 de 5







A

A

A

A

A

A

A

A

A

A

A

A

A

A

A

A

A

A

A

A

A

13

12

11

10

09

08

07

06

995

428

10428

373

996

375

374

376

372

10371

ESCALAS: FECHA:


HOJA:

ORIGINAL A-3

Nº DE PLANO:


¯


4 de 5







A

A

A

A

A

A

A

A

A

A

A

A

A

A

A

A

A

A

A

A

A02

01

20371

428

10371

1015

10428

996

373

997

463

30371

372

997

370

374

740

ESCALAS: FECHA:


HOJA:

ORIGINAL A-3

Nº DE PLANO:


¯


5 de 5







DDDD

DDDD

DDDD

DDDD

DDDD

DDDD

DDDD

DDDD

DDDD

DDDD

DDDD

DDDD

DD

DD

DD

DD

DD

DD

DD

DD

DD

DD

DD

DD

DDDD

DDDD

DDDD

DDDD

DDDD

DDDD

DDDD

DD

DD

DD

DD

DD

DD

DD

DD

DD

DD

DD

DD

DD

DD

DD

DD

DD

DD

DD

DD

DD

DD

DD

DD

DD

DD

DD

DD

DD

DD

DD

DD

DD

DD

DD

DD

DD

DD

DD

DD

DDDD

DDDD

DDDD

DD

DD

DDDD

DD

DD

DD

DD

DD

DD

DDDD

DDDD

DDDD

DDDD

DD

DD

DD

DDDD

DDDD

DDDD

DDDD

DDDD

DDDD

DDDD

DDDD

DDDD

DDDD

DDDD

DDDD

DDDD

DDDD

DDDD

DD

DD

DD

DD

DD

DD

DD

DD

DD

DD

DD

DDDD

DDDD

DDDD

DDDD

DDDD

DDDD

DDDD

DDDD

DDDD

DDDD

DDDD

DDDD

DDDD

DDDD

DDDD

DDDD

DD

DD

DD

DD

DD

DD

DDDD

DDDD

DDDD

DDDD

DDDD

DDDD

DDDD

DDDD

DD

DD

DD

DDDD

DD

DD

DD

DD

DD

DD

DD

DD

DDDD

DDDD

DDDD

DDDD

DDDD

DDDD

DDDD

DD

DD

DD

DD

DD

DD

DD

DD

DDDD

DD

DD

DDDD

DDDD

DD

DD

DD

DD

DD

DD

DD

DD

995996

1018

1019

Arrubal

ESCALAS: FECHA:


HOJA:

ORIGINAL A-3

Nº DE PLANO:

20 50 100Meters

1 de 1PARCELAS AFECTADAS.

¯

LEYENDA:DD DD Límites Municipales Área Parcelas Afectadas Límite D.P.H Parcelas_estudio Parcelario

POLIGONO N_PARCELA SUP. AFECTADA (ha) AFECCIÓN (%) AFECCIÓN1 1018 0,204 40,27% Parcial en D.P.H.Técnico1 995 1,994 23,04% Parcial en D.P.H.Técnico1 996 0,03 0,83% Parcial en D.P.H.Técnico




* La parcela 1019 se encuentra Fuera del D.P.H. Técnico.




ANEXO Nº 6. REPORTAJE FOTOGRÁFICO


ANEXO Nº6.

REPORTAJE FOTOGRÁFICO






1

1. LOCALIZACIÓN GENERAL DE LAS IMÁGENES






2

2. REPORTAJE FOTOGRÁFICO

FIGURA 03: Parcela 1018. Panorámica hacia aguas arriba desde el punto 01. El camino correspondería

aproximadamente a la zona de servidumbre del DPH, el cual ocuparía todo el espacio desarrollado entre el camino y el río Ebro, localizado tras los árboles visibles a la derecha de la imagen

FIGURA 04: Parcela 1018. Panorámica hacia el río Ebro desde el punto 01. El registro visible, asociado a la

conducción del bombeo de la central eléctrica de ciclo combinado de Arrúbal, se localiza en DPH y dentro de la parcela. El límite entre el DPH y su zona de servidumbre coincide con la cuneta del camino.






3

FIGURA 05: Parcela 1018. Panorámica hacia aguas abajo desde el punto 02. La cuneta de la derecha del

camino coincide aproximadamente con el límite interior de la parcela, la cual se extiende más allá del camino ocupando parte del DPH del Ebro.

FIGURA 06: Parcelas 1018 y 1019. Límite entre las parcelas 1018, izquierda, y 1019, derecha, en el punto

03. El camino desde el que se toma la imagen separa ambas parcelas.






4

FIGURA 07: Parcela 995. Panorámica hacia el Ebro desde el punto 04. El límite del DPH se localizaría

aproximadamente a la altura de la primera línea de vegetación arbustiva; al fondo el soto asociado al galacho que se conserva en la zona, ya en DPH.

FIGURA 08: Parcela 995. Panorámica hacia la parcela 996 desde el punto 04. Al fondo, soto desarrollado

entre ambas parcelas, fuera del DPH asociado al río Ebro.






5

FIGURA 09: Parcela 996. Vista de la parcela 996 desde el punto 05. El límite del DPH se situaría

aproximadamente sobre el estribo del paso fotografiado.

FIGURA 10: Acequia, con dominio público asociado, localizada entre las parcelas 995 y 996. Vista hacia

aguas abajo desde el punto 06.






6

FIGURA 11: Parcela 996. Escarpe que define la linde sur de la parcela visto desde el punto 05.




MEMORIA DESCRIPTIVA DEL PROCEDIMIENTO … Deslinde Arrubal.pdf · memoria descriptiva del...

Documents

Transcript of MEMORIA DESCRIPTIVA DEL PROCEDIMIENTO … Deslinde Arrubal.pdf · memoria descriptiva del...