Ver el código de «Crítica al embalse del Ciguiñuela»

{{ TrabajoSIG | Crítica al embalse del Ciguiñuela | Álvaro Cáceres Alonso <br /> Javier Díez Ávila <br /> Gonzalo Marinas Sanz | [[:Categoría:SIGAIC_14/15|Curso 14/15]] }}

Tras la propuesta del presidente de la Confederación Hidrográfica del Duero
de la construcción del embalse en el río Ciguiñuela (afluente del Eresma)
para solucionar los problemas que existen de abastecimiento en la provincia
de Segovia, se quiere analizar la viabilidad del proyecto y su impacto en el
territorio.

Para ello, se procederá a generar el embalse conocido el volumen propuesto
y la cerrada. Una vez hecho esto, se detallarán los terrenos, municipios y construcciones
afectadas, así como posibles zonas de protección medioambiental
dañadas. Además, se calculará el coste de las expropiaciones y se comentará
la litología encima de la cual se va a embalsar. Finalmente, se hará un somero
estudio hidrológico como complemento a su estudio.

Durante el análisis, se usará cartografía del IGN (MTN25, MDT05, BTN25...)
así como capas del Catastro y mapas del IGME. También se utilizarán archivos
del Geoportal de la IDE de Castilla y León, de la Diputación Provoncial
de Segovia, y datos del INE, de AEMET y de la CHD.

Se esperan obtener mapas con el embalse y las superficies expropiadas,
infraestructuras y zonas protegidas afectadas. . . con los cuales sea posible
llevar a cabo una crítica a grandes rasgos de la necesidad o no de la obra.

== Introducción ==
[[Archivo:SIG Marinas 01.png|500x353px|marco|derecha|Situación geográfica]]
La ciudad de Segovia, a día de hoy, se abastece gracias a dos embalses situados en sus alrededores:
* Pontón Alto (en el río Eresma y en los municipios de La Granja de San Ildefonso y Palazuelos de Eresma): 7,41hm3 [1]
*  Puente Alta (en el río de la Acebeda en el municipio de Revenga): 2,50hm3 [2]
En el proyecto de “Regulación de caudales para el abastecimiento de Segovia y otras poblaciones de la cuenca del Eresma”, planteado por el Ministerio de Agricultura, Alimentación y Medio Ambiente, se propone la construcción de un embalse de 29hm3 en
el río Ciguiñuela (tributario del Eresma, tal y como se ve en la figura 1), que cuadruplicará
la cantidad de agua disponible para abastecimiento (de los 9,91hm3 actuales a
los 38,91hm3 futuros). Evidentemente, tal incremento para una ciudad y sus municipios
circundantes que han tenido un crecimiento poblacional muy discreto en los últimos
años, está injustificado, por lo que buena parte del agua del nuevo embalse se destinará
a reducir la presión sobre algunos acuíferos de Castilla y León, muy castigados hoy día
por las dotaciones requeridas para regadío). [3]

La crítica de esta infraestructura hidráulica se articulará en diferentes apartados, que
considerarán su emplazamiento geográfico y las posibles repercusiones a su entorno, así
como la viabilidad económica y técnica del proyecto. Para ello, se empleará principalmente
el programa de código libre QuantumGIS.
Por último, es necesario mencionar los datos iniciales disponibles a partir de los cuales
se realizarán todas las operaciones:

1. Posición aproximada de la presa [4]

2. Volumen del embalse: 29 hm3[3, 5]

3. Superficie ocupada por la lámina libre: 200 ha [3, 5]

== Metodología ==

=== Digitalización de la presa y tanteo de inundaciones ===

La primera operación a realizar es digitalizar la presa mediante una línea. Con ese fin,
a partir del MTN25 hoja 0483-1 y el MDT05-LIDAR hoja 0483 [6], se dibuja de forma
grosera el eje de la presa con el primer mapa, y se realiza un ajuste más fino una vez
extraídas las curvas de nivel metro a metro del segundo.

Para limitar las superficies de inundación a la presa y la cuenca del Ciguiñuela aguas
arriba de ella (y no obtener las inundaciones en todo la extensión de la malla LIDAR), se
genera un polígono coincidente exactamente con la cerrada y las crestas que delimitan las
hondonadas del río y sus afluentes más próximos. Después, se cargan en GRASS las capas
adecuadas y se rasteriza el polígono para poder operar con él. Mediante la herramienta
r.mapcalc se opera para introducir en el polígono ráster las elevaciones asociadas a los
píxeles equivalentes del MDT05.

Por último, se tantean distintas inundaciones con r.lake.seed para obtener el embalse
dato y se vectoriza y exporta el resultado para operar con él más cómodamente a
posteriori.

=== Zonificación del espacio fluvial ===

Para calcularlas, basta con realizar sendos buffers de cinco

=== Municipios afectados ===

Lo único que hay que hacer en este punto es cargar la capa con los municipios de
Segovia del CartoCiudad [6] e intersecarla dos veces:

1. Previamente con un rectángulo que delimita la zona de estudio a la extensión de
la hoja 0483-1 del MTN25

2. Posteriormente con la capa vectorial del embalse

El área afectada se obtiene añadiendo la columna apropiada en la tabla de atributos de
la última capa hallada.

=== Catastro ===

A la hora de analizar las parcelas a expropiar para establecer el embalse, se recurre a los
archivos .shp que proporciona el Catastro con las parcelas de cada municipio afectado [7].
Una vez cargadas, se añade una columna con el nombre del pueblo al que pertenecen, se
intersecan con la lámina del embalse y se exporta la tabla de atributos a Excel. Después,
estableciendo la condición de que aquellas parcelas ocupadas en una superficie mayor al
30 % serán expropiadas completamente, se calcula la superficie expropiada y su coste
total.

=== Geología ===

Lo primero que se hace es descargar la hoja 0483 georreferenciada del MAGNA 50 [14],
para después superponer encima el embalse obtenido en el apartado 2.1. Finalmente se
añaden unas anotaciones de texto para suplir la ausencia de la leyenda original.

=== Superficie vegetal ===

Para analizar la cubierta terrestre solo es preciso cargar las dos capas vectoriales
(puntual y lineal) con esa información [8].

=== Medio ambiente ===

Las afecciones medioambientales se han estudiado desde varios puntos de vista que se
detallan a continuación.

==== LIC y ZEPA ====

Para comprobar si la infraestructura afecta o no a estos lugares, se abren las capas que
los incluyen y se eliminan todos los que no sean cercanos a la zona de estudio. Existen
varias fuentes de las cuales se pueden obtener, pero aquí se ha optado por el Ministerio
de Agricultura, Alimentación y Medio Ambiente (MAGRAMA) por tener la información
más actualizada [9].

==== Zona forestal de valor ====

La Diputación de Segovia define en su cartografía provincial (hoja 45) zonas forestales
de valor [10] fácilmente digitalizables a partir de la hoja 0483 del PNOA [6]. Estas zonas
se observan bien en la fotografía georreferenciada y cargada en QGIS con el complemento
Photo2Shape.

==== Vías pecuarias ====

Cerca de la ubicación del embalse se encuentra la Cañada Real Soriana Occidental,
por lo que se considera preciso comprobar la interrupción o no de las sendas para el
ganado. De la IDE de Castilla y León se toman los caminos por la zona del embalse
[12], y se consideran aquellos pertenecientes a la Red Nacional de Vías Pecuarias según
el MAGRAMA [11]. Por otra parte, de la BTN25 [6] se abre el .shp perteneciente a este
ámbito.

==== Aguas subterráneas ====

Se hace uso de dos capas, obtenidas del MAGRAMA [13] superpuestas al mapa de
España [6]:
* Aguas subterráneas
* Acuíferos sobreexplotados

=== Estudio hidrológico ===

Calculada la extensión del embalse, lo primero que se necesita es su cuenca para
hacer el estudio pluviométrico. Para ello, partiendo de la malla LIDAR y con GRASS, se
realizan las cuencas con el comando r.watershed, que se convierten a capa vectorial y
se exportan a QGIS. A partir de esta capa, se seleccionan las subcuencas que drenan en el
embalse a partir de la capa de ríos (obtenida de la BTN25) y del embalse. Seleccionando
como objetos espaciales las cuencas necesarias, se realiza la diferencia entre la capa de
subcuencas y estos objetos, obteniendo la cuenca del embalse.

Creada la cuenca, se calcula la precipitación media anual a partir de los polígonos
de Thiessen o de Voronoi. Los pluviómetros de la provincia de Segovia (aquellos que
dan datos mensuales) [15] se descargan en una tabla Excel, que se guarda con formato
.txt para, acto seguido, crear una capa desde un fichero de texto delimitado con sus
coordenadas. Con el centroide (obtenido con las herramientas de geometría de objetos
vectoriales) del embalse se realiza un buffer de 15 km, que interseca a la de pluviómetros.
Con esta nueva capa, se calculan los polígonos de Vornoi, y, realizando la intersección
con la cuenca del embalse, se tienen los datos necesarios para el cálculo de la lluvia de
proyecto y el volumen de lluvia anual, conociendo la superficie de las subcuencas y la
precipitación media de cada zona (a partir de los datos pluviométricos obtenidos de la
AEMET).

Para hallar los perfiles de los ríos, se hace uso de la herramienta Profile Tool seleccionando
como polilínea la correspondiente a los ríos que atraviesan la subcuenca. Para
el umbral de escorrentía, se carga la capa del MAGRAMA [16].

=== Eutrofización ===

Para este apartado solo se calcula la cota media del embalse a partir de la herramienta
ráster Estadísticas de zona.

== Resultados ==

=== Digitalización de la presa y tanteo de inundaciones ===
[[Archivo:SIG Marinas 02.png|300x212px|marco|derecha|Figura 2. Presa]]
En la figura 2 se ve la posición final de la presa una vez encajada adecuadamente. Es
necesario mencionar que, para la cota 1093, se obtiene una longitud de cerrada de casi
1500m.

Para generar las inundaciones, se emplea el polígono ráster con las elevaciones que se
muestra en la figura 3, tal y como se indica en el apartado 2.1.

Después, se realizan varios tanteos con diferentes cotas y los resultados expuestos por
r.lake.seed se comparan con los datos iniciales:

* Inundación a la cota 1085:

<small> 
::Profundidad del lago de 0.00 a 43.8149 metros
::Superficie del lago: 1429761.82 metros cuadrados
::Volumen del lago: 19646003 metros cubicos</small>

* Inundación a la cota 1090:

<small> 
::Profundidad del lago de 0.00 a 48.8149 metros
::Superficie del lago: 1867055.83 metros cuadrados
::Volumen del lago: 27842255 metros cubicos</small>

[[Categoría:SIGAIC_14/15]]

* Inundación a la cota 1091:

<small> 
::Profundidad del lago de 0.00 a 49.8149 metros
::Superficie del lago: 1974197.51 metros cuadrados
::Volumen del lago: 29761121 metros cubicos</small>

Aunque, llegados a este punto, está bastante claro por dónde debe andar el Nivel Máximo
Normal (NMN), es necesario observar (en la figura ) cómo se produce el desbordamiento
del agua por la margen izquierda del río al aumentar la cota de la superficie libre a 1092.
[[Archivo:SIG Marinas 03.png|500x353px|marco|centro|Figura 4. Inundación de la presa]]
Para ver este fenómeno, es necesario considerar otro polígono de referencia con las
elevaciones del terreno.

Como se señala en el apartado 4, se toma como Nivel de Avenida de Proyecto la cota
1090 (lo cual será usado para los siguientes cálculos), aunque eso suponga no alcanzar
las condiciones de contorno establecidas.
[[Archivo:SIG Marinas 04.png|500x353px|marco|centro|Figura 5. Embalse a la cota 1090]]

=== Zonificación del espacio fluvial ===

Antes de nada, es necesario recordar los conceptos asociados a este punto.

:'''Dominio público''' hidráulico es el conjunto de las aguas continentales superficiales y subterráneas renovables, los cauces de corrientes naturales continuas y discontinuas, los lechos de lagos y lagunas o embalses superficiales en el cauce público y los acuíferos subterráneos, delimitados a la máxima crecida ordinaria (en este caso, 1090).

:'''Zona de servidumbre''' es la franja de protección de cinco metros de anchura al dos márgenes del cauce, que son de uso público y que tiene la finalidad de proteger el ecosistema fluvial y el dominio público hidráulico, de paso público peatonal y para el desarrollo de los servicios de vigilancia, conservación y salvamento, y para el amarre de embarcaciones de manera ocasional o en caso de necesidad.

:'''Zona de policía''' es la constituida por una franja lateral de cien metros de anchura a cada lado, contados a partir de la línea que delimita el cauce, en las que se condiciona el uso del suelo y las actividades que en él se desarrollen.

Todos estas superficies se delimitan en la figura 6.

[[Archivo:SIG Marinas 06.png|500x352px|marco|centro|Figura 6: Zonificacióon del espacio fluvial]]

=== Municipios afectados ===
[[Archivo:SIG Marinas 07.png|200x141px|marco|izquierda|Figura 7: Discrepancias entre límites municiapales]]
A la hora de analizar este punto, lo primero que llama la atención es la discrepancia
existente entre las diversas entidades oficiales a la hora de delimitar las poblaciones,
llamando la atención el solape de superficies según el Catastro. En la figura 7 se detallan
los lindes según el IGN (Cartociudad y MTN25) y el Catastro, y se observa que se
producen desfases de hasta 50m.

Entrando ahora en materia, en la figura se ve cómo la lámina de agua a la cota 1090,
solapándola a la capa del Cartociudad, inunda tres términos municipales :

:1. Espirdo: 86,02 ha

:2. Trescasas: 74,95 ha

:3. Torrecaballeros: 25,73 ha

[[Archivo:SIG Marinas 08.png|500x353px|marco|centro|Figura 8: Municipios afectados y superficies [ha]]]

=== Catastro ===

Para la construcción del embalse habrá que expropiar una serie de parcelas, que se
muestran en la figura 9. Aquellas fincas no completamente cubiertas por el embalse,
se expropian en toda su extensión si están ocupadas en más del 30 %, obteniendo los
siguientes resultados: