Propuesta metodológica para la identificación de medidas de adaptación al cambio climático en sistemas de recursos hídricos

Esta tesis realiza una contribución metodológica en el estudio de medidas de adaptación potencialmente adecuadas a largo plazo, donde los sistemas de recursos hídricos experimentan fuertes presiones debido a los efectos del cambio climático. Esta metodología integra el análisis físico del sistema, basándose en el uso de indicadores que valoran el comportamiento de éste, y el análisis económico mediante el uso del valor del agua. El procedimiento metodológico inicia con la construcción de un conjunto de escenarios futuros, que capturan por un lado las características de variabilidad de las aportaciones de diversos modelos climáticos y, por otro, las características hidrológicas de la zona de estudio. Las zonas de estudio seleccionadas fueron las cuencas del Guadalquivir, Duero y Ebro y se utilizaron como datos observados las series de escorrentía en régimen natural estimadas por el modelo SIMPA que está calibrado en la totalidad del territorio español. Estas series observadas corresponden al periodo 1961-1990. Los escenarios futuros construidos representan el periodo 2071-2100. La identificación de medidas de adaptación se apoyó en el uso de indicadores que sean capaces de caracterizar el comportamiento de un sistema de recursos hídricos frente a los efectos del cambio climático. Para ello se seleccionaron los indicadores de calidad de servicio (I1) y de confiabilidad de la demanda (I2) propuestos por Martin-Carrasco et al. (2012). Estos indicadores valoran el comportamiento de un sistema mediante la identificación de los problemas de escasez de agua que presente, y requieren para su cuantificación el uso de un modelo de optimización. Para este estudio se ha trabajado con el modelo de optimización OPTIGES. La determinación de estos indicadores fue realizada para análisis a corto plazo donde los efectos del cambio climático no son de relevancia, por lo que fue necesario analizar su capacidad para ser usados en sistemas afectados por dichos efectos. Para este análisis se seleccionaron tres cuencas españolas: Guadalquivir, Duero y Ebro, determinándose que I2 no es adecuado para este tipo de escenarios. Por ello se propuso un nuevo indicador “Indicador de calidad de servicio bajo cambio climático” (I2p) que mantiene los mismos criterios de valoración que I2 pero que responde mejor bajo fuertes reducciones de aportaciones producto del cambio climático. La metodología propuesta para la identificación de medidas de adaptación se basa en un proceso iterativo en el cual se van afectando diversos elementos que conforman el esquema del sistema bajo acciones de gestión previamente identificadas, hasta llegar a un comportamiento óptimo dado por el gestor. Las mejoras de estas afectaciones son cuantificadas mediante los indicadores I1 e I2p, y de este conjunto de valores se selecciona la que se acerca más al comportamiento óptimo. Debido a la extensa cantidad de información manejada en este análisis, se desarrolló una herramienta de cálculo automatizada en Matlab. El proceso seguido por esta herramienta es: (i) Ejecución del modelo OPTIGES para las diferentes modificaciones por acciones de gestión; (ii) Cálculo de los valores de I1 e I2p para cada una de estas afectaciones; y (iii) Selección de la mejor opción. Este proceso se repite hasta llegar al comportamiento óptimo buscado, permitiendo la identificación de las medidas de adaptación mas adecuadas. La aplicación de la metodología para la identificación de medidas de adaptación se realizó en la cuenca del Guadalquivir, por ser de las tres cuencas analizadas bajo los indicadores I1 e I2p la que presenta los problemas más serios de escasez de agua. Para la identificación de medidas de adaptación se analizaron dos acciones de gestión: 1) incremento de los volúmenes de regulación y 2) reducción de las demandas de riego, primero bajo la valoración del comportamiento físico del sistema (análisis de sensibilidad) permitiendo identificar que la primera acción de gestión no genera cambios importantes en el comportamiento del sistema, que si se presentan bajo la segunda acción. Posteriormente, con la acción que genera cambios importantes en el comportamiento del sistema (segunda acción) se identificaron las medidas de adaptación más adecuadas, mediante el análisis físico y económico del sistema. Se concluyó que en la cuenca del Guadalquivir, la acción de reducción de las demandas de riego permite minimizar e incluso eliminar los problemas de escasez de agua que se presentarían a futuro bajo diferentes proyecciones hidrológicas, aunque estas mejoras implicarían fuertes reducciones en dichas demandas. Siendo las demandas más afectadas aquellas ubicadas en cabecera de cuenca. Los criterios para la reducción de las demandas se encuentran en función de las productividades y garantías con las que son atendidas dichas demandas. This thesis makes a methodological contribution to the study of potentially suitable adaptation measures in the long term, where water resource systems undergo strong pressure due to the effects of climate change. This methodology integrates the physical analysis of the system, by the use of indicators which assess its behavior, and the economic analysis by the use of the value of water. The methodological procedure begins with the building of a set of future scenarios that capture, by one hand, the characteristics and variability of the streamflow of various climate models and, on the other hand, the hydrological characteristics of the study area. The study areas chosen were the Guadalquivir, Ebro and Duero basins, and as observed data where used runoff series in natural regimen estimated by the SIMPA model, which is calibrated in the whole Spanish territory. The observed series are for the 1961-1990 period. The future scenarios built represent the 2071-2100 periods. The identification of adaptation measures relied on the use of indicators that were able of characterize the behavior of one water resource system facing the effects of climate change. Because of that, the Demand Satisfaction Index (I1) and the Demand Reliability Index (I2) proposed by Martin-Carrasco et al. (2012) were selected. These indicators assess the behavior of a system by identifying the water scarcity problems that it presents, and require in order to be quantified the use of one optimization model. For this study the OPTIGES optimization model has been used. The determination of the indicators was made for the short-term analysis where the climates change effect are not relevant, so it was necessary to analyze their capability to be used in systems affected by those these. For this analysis three Spanish basins were selected: Guadalquivir, Duero and Ebro. It was determined that the indicator I2 is not suitable for this type of scenario. It was proposed a new indicator called “Demand Reliability Index under climate change” (I2p), which keeps the same assessment criteria than I2, but responsive under heavy reductions of streamflow due to climate change. The proposed methodology for identifying adaptation measures is based on an iterative process, in which the different elements of the system´s schema are affected by previously defined management actions, until reach an optimal behavior given by the manager. The improvements of affectations are measured by indicators I1 e I2p, and from this set of values it is selected the affectation that is closer to the optimal behavior. Due to the large amount of information managed in this analysis, it was developed an automatic calculation tool in Matlab. The process followed by this tool is: Firstly, it executes the OPTIGES model for the different modifications by management actions; secondly, it calculates the values of I1 e I2p for each of these affectations; and finally it chooses the best option. This process is performed for the different iterations that are required until reach the optimal behavior, allowing to identify the most appropriate adaptation measured. The application of the methodology for the identification of adaptation measures was conducted in the Guadalquivir basin, due to this was from the three basins analyzed under the indicators I1 e I2p, which presents the most serious problems of water scarcity. For the identification of adaptation measures there were analyzed two management actions: 1) To increase the regulation volumes, and 2) to reduce the irrigation demands, first under the assessment of the physical behavior of the system (sensibility analysis), allowing to identify that the first management action does not generate significant changes in the system´s behavior, which there are present under the second management action. Afterwards, with the management action that generates significant changes in the system´s behavior (second management action), there were identified the most adequate adaptation measures, through the physical and economic analysis of the system. It was concluded that in the Guadalquivir basin, the action of reduction of irrigation demands allows to minimize or even eliminate the water scarcity problems that could exist in the future under different hydrologic projections, although this improvements should involve strong reductions of the irrigation demands. Being the most affected demands those located in basins head. The criteria for reducing the demands are based on the productivities and reliabilities with which such demands are meet.

Encauzamiento de la cabecera del Arroyo del Chorreadero en el municipio de Villargordo

El objeto del presente proyecto es definir las obras para la implantación de colectores de pluviales en la población de Villargordo con el fin de evitar que los caudales, por falta de capacidad de la red de saneamiento existente, fluyan por las calles y ocasiones daños en las viviendas y otros edificios. Para la recogida de las aguas vertidas por las cuencas se proyectarán unos depósitos(o arquetas) de pluviales con rejas de dimensiones suficientes en la zona del Llano de las Viñas, en el lugar donde se encuentran en la actualidad. La inversión en las obras de construcción del encauzamiento serán amortizadas en un breve periodo de tiempo debido a que, de manera casi anual, los destrozos ocasionados por las riadas y por los grandes caudales que circulan por las calles de la población son considerablemente altos. Con esta actuación se pretende eliminar estos desperfectos y, sobre todo, las grandes afecciones que supone a una población, sobre todo económicas. También hay que evitar que la red de saneamiento existente entre en carga por la presencia de dichos caudales, lo que conlleva la inundación de las partes bajas de los edificios. Por ello esta conducción servirá únicamente para aliviar las puntas de daudal que se producen para el máximo caudal calculado. Esta medida se tomara dimensionando de manera que las conducciones proyectadas no entren en carga para el caudal de proyecto que corresponde al periodo de retorno tomado en este tipo de obras, que suele ser de 100 años.

Estimate of the total ice volume of Svalbard glaciers and their potential contribution to sea-level rise = Estimación del volumen de hielo de Svalbard y su contribución potencial a la subida del nivel del mar

El objetivo final de las investigaciones recogidas en esta tesis doctoral es la estimación del volumen de hielo total de los ms de 1600 glaciares de Svalbard, en el Ártico, y, con ello, su contribución potencial a la subida del nivel medio del mar en un escenario de calentamiento global. Los cálculos más exactos del volumen de un glaciar se efectúan a partir de medidas del espesor de hielo obtenidas con georradar. Sin embargo, estas medidas no son viables para conjuntos grandes de glaciares, debido al coste, dificultades logísticas y tiempo requerido por ellas, especialmente en las regiones polares o de montaña. Frente a ello, la determinación de áreas de glaciares a partir de imágenes de satélite sí es viable a escalas global y regional, por lo que las relaciones de escala volumen-área constituyen el mecanismo más adecuado para las estimaciones de volúmenes globales y regionales, como las realizadas para Svalbard en esta tesis. Como parte del trabajo de tesis, hemos elaborado un inventario de los glaciares de Svalbard en los que se han efectuado radioecosondeos, y hemos realizado los cálculos del volumen de hielo de más de 80 cuencas glaciares de Svalbard a partir de datos de georradar. Estos volúmenes han sido utilizados para calibrar las relaciones volumen-área desarrolladas en la tesis. Los datos de georradar han sido obtenidos en diversas campañas llevadas a cabo por grupos de investigación internacionales, gran parte de ellas lideradas por el Grupo de Simulación Numérica en Ciencias e Ingeniería de la Universidad Politécnica de Madrid, del que forman parte la doctoranda y los directores de tesis. Además, se ha desarrollado una metodología para la estimación del error en el cálculo de volumen, que aporta una novedosa técnica de cálculo del error de interpolación para conjuntos de datos del tipo de los obtenidos con perfiles de georradar, que presentan distribuciones espaciales con unos patrones muy característicos pero con una densidad de datos muy irregular. Hemos obtenido en este trabajo de tesis relaciones de escala específicas para los glaciares de Svalbard, explorando la sensibilidad de los parámetros a diferentes morfologías glaciares, e incorporando nuevas variables. En particular, hemos efectuado experimentos orientados a verificar si las relaciones de escala obtenidas caracterizando los glaciares individuales por su tamaño, pendiente o forma implican diferencias significativas en el volumen total estimado para los glaciares de Svalbard, y si esta partición implica algún patrón significativo en los parámetros de las relaciones de escala. Nuestros resultados indican que, para un valor constante del factor multiplicativo de la relacin de escala, el exponente que afecta al área en la relación volumen-área decrece según aumentan la pendiente y el factor de forma, mientras que las clasificaciones basadas en tamaño no muestran un patrón significativo. Esto significa que los glaciares con mayores pendientes y de tipo circo son menos sensibles a los cambios de área. Además, los volúmenes de la población total de los glaciares de Svalbard calculados con fraccionamiento en grupos por tamaño y pendiente son un 1-4% menores que los obtenidas usando la totalidad de glaciares sin fraccionamiento en grupos, mientras que los volúmenes calculados fraccionando por forma son un 3-5% mayores. También realizamos experimentos multivariable para obtener estimaciones óptimas del volumen total mediante una combinación de distintos predictores. Nuestros resultados muestran que un modelo potencial simple volumen-área explica el 98.6% de la varianza. Sólo el predictor longitud del glaciar proporciona significación estadística cuando se usa además del área del glaciar, aunque el coeficiente de determinación disminuye en comparación con el modelo más simple V-A. El predictor intervalo de altitud no proporciona información adicional cuando se usa además del área del glaciar. Nuestras estimaciones del volumen de la totalidad de glaciares de Svalbard usando las diferentes relaciones de escala obtenidas en esta tesis oscilan entre 6890 y 8106 km3, con errores relativos del orden de 6.6-8.1%. El valor medio de nuestras estimaciones, que puede ser considerado como nuestra mejor estimación del volumen, es de 7.504 km3. En términos de equivalente en nivel del mar (SLE), nuestras estimaciones corresponden a una subida potencial del nivel del mar de 17-20 mm SLE, promediando 19_2 mm SLE, donde el error corresponde al error en volumen antes indicado. En comparación, las estimaciones usando las relaciones V-A de otros autores son de 13-26 mm SLE, promediando 20 _ 2 mm SLE, donde el error representa la desviación estándar de las distintas estimaciones. ABSTRACT The final aim of the research involved in this doctoral thesis is the estimation of the total ice volume of the more than 1600 glaciers of Svalbard, in the Arctic region, and thus their potential contribution to sea-level rise under a global warming scenario. The most accurate calculations of glacier volumes are those based on ice-thicknesses measured by groundpenetrating radar (GPR). However, such measurements are not viable for very large sets of glaciers, due to their cost, logistic difficulties and time requirements, especially in polar or mountain regions. On the contrary, the calculation of glacier areas from satellite images is perfectly viable at global and regional scales, so the volume-area scaling relationships are the most useful tool to determine glacier volumes at global and regional scales, as done for Svalbard in this PhD thesis. As part of the PhD work, we have compiled an inventory of the radio-echo sounded glaciers in Svalbard, and we have performed the volume calculations for more than 80 glacier basins in Svalbard from GPR data. These volumes have been used to calibrate the volume-area relationships derived in this dissertation. Such GPR data have been obtained during fieldwork campaigns carried out by international teams, often lead by the Group of Numerical Simulation in Science and Engineering of the Technical University of Madrid, to which the PhD candidate and her supervisors belong. Furthermore, we have developed a methodology to estimate the error in the volume calculation, which includes a novel technique to calculate the interpolation error for data sets of the type produced by GPR profiling, which show very characteristic data distribution patterns but with very irregular data density. We have derived in this dissertation scaling relationships specific for Svalbard glaciers, exploring the sensitivity of the scaling parameters to different glacier morphologies and adding new variables. In particular, we did experiments aimed to verify whether scaling relationships obtained through characterization of individual glacier shape, slope and size imply significant differences in the estimated volume of the total population of Svalbard glaciers, and whether this partitioning implies any noticeable pattern in the scaling relationship parameters. Our results indicate that, for a fixed value of the factor in the scaling relationship, the exponent of the area in the volume-area relationship decreases as slope and shape increase, whereas size-based classifications do not reveal any clear trend. This means that steep slopes and cirque-type glaciers are less sensitive to changes in glacier area. Moreover, the volumes of the total population of Svalbard glaciers calculated according to partitioning in subgroups by size and slope are smaller (by 1-4%) than that obtained considering all glaciers without partitioning into subgroups, whereas the volumes calculated according to partitioning in subgroups by shape are 3-5% larger. We also did multivariate experiments attempting to optimally predict the volume of Svalbard glaciers from a combination of different predictors. Our results show that a simple power-type V-A model explains 98.6% of the variance. Only the predictor glacier length provides statistical significance when used in addition to the predictor glacier area, though the coefficient of determination decreases as compared with the simpler V-A model. The predictor elevation range did not provide any additional information when used in addition to glacier area. Our estimates of the volume of the entire population of Svalbard glaciers using the different scaling relationships that we have derived along this thesis range within 6890-8106 km3, with estimated relative errors in total volume of the order of 6.6-8.1% The average value of all of our estimates, which could be used as a best estimate for the volume, is 7,504 km3. In terms of sea-level equivalent (SLE), our volume estimates correspond to a potential contribution to sea-level rise within 17-20 mm SLE, averaging 19 _ 2 mm SLE, where the quoted error corresponds to our estimated relative error in volume. For comparison, the estimates using the V-A scaling relations found in the literature range within 13-26 mm SLE, averaging 20 _ 2 mm SLE, where the quoted error represents the standard deviation of the different estimates.

E.D.A.R. Renedo de Piélagos (Cantabria)

Las obras definidas en este proyecto están integradas dentro del “Plan Director de Saneamiento y Depuración de Aguas Residuales de Cantabria”, que contempla la actuación en materia de saneamiento por determinadas zonas o cuencas hidrográficas, de acuerdo con la Directiva 91/271/CEE. La Estación Depuradora de Aguas Residuales (E.D.A.R) de Renedo de Piélagos, perteneciente al “Saneamiento de la Cuenca Media del Pas-Pisueña” tiene por objeto el tratamiento de los vertidos de aguas residuales de los municipios de Santa María de Cayón, Castañeda, Puente Viesgo y Piélagos de modo que se garantice una adecuada calidad de las aguas de los ríos Pisueña y Pas. El objeto del presente proyecto es la completa definición del conjunto de obras necesarias para el correcto funcionamiento del tratamiento de aguas residuales procedentes de los municipios de la cuenca media del Pas-Pisueña. De acuerdo con los objetivos generales perseguidos por las actuaciones a realizar, se pretende dotar a la zona de una Estación Depuradora de Aguas Residuales en Renedo de Piélagos necesaria para asegurar el correcto funcionamiento de la red. De esta forma se intenta evitar los vertidos directos a los ríos y arroyos de la zona, y de este modo, garantizar la calidad ambiental del medio receptor, que una vez finalizadas las actuaciones será el río Pas. El ámbito de actuación de este proyecto abarca los términos municipales de Puente Viesgo, Castañeda y Piélagos, aunque el colector interceptor principal, que dará servicio a la E.D.A.R. de Renedo de Piélagos, se diseña con capacidad suficiente para recoger una futura incorporación de los caudales de aguas residuales del municipio de Sta. María de Cayón.

La habitación subterránea en la Península Ibérica : investigación geológica

La habitación rupestre en la Península Ibérica conforma un amplio conjunto de manifestaciones cuyos orígenes son difíciles de dilucidar. Existen conjuntos primitivos, posiblemente horadados durante la Antigüedad Clásica y Tardía, y otros que probablemente se originaron y excavaron en la Alta Edad Media, en los que se mezclan grupos de origen religioso y otros de probable uso defensivo. El conjunto peninsular es sin duda el más variado de Europa Occidental, pues recibió influencias árabes y usos de la cristiandad antigua, importados de Oriente Próximo, que por tanto relacionan estas manifestaciones con otras del arco mediterráneo. La supervivencia del uso de las cuevas a lo largo de la Baja Edad Media -una vez afianzada la Reconquista- es una incógnita, pero los usos rupestres volvieron a generalizarse en ciertos sectores de la Península durante la Edad Moderna, toda vez que el hábitat cuevero manifestó una eclosión relacionada con la peripecia de los moriscos, primero desterrados de sus habitaciones en el Sur y Este peninsular y luego expulsados en los albores del siglo XVII. Los que quedaron, nominalmente conversos, debieron habitar de nuevo cuevas en sus lugares de origen. Esos núcleos son los más abundantes, particularmente en la provincia de Granada y aledañas. Este substrato fue seguramente el punto de apoyo para la proliferación de las cuevas de habitación a partir del siglo XIX, en el cual un creciente proletariado agrícola y urbano necesitó de alojamientos baratos y no hizo sino imitar usos preexistentes, que se pueden rastrear en muchos de los núcleos rupestres que sobrevivieron mayoritariamente hasta bien entrado el siglo XX, y que se abandonaron gradualmente a partir de los años 60 de esa centuria. Para entonces, existían barrios de cuevas extensos en muchas provincias, destacando, aparte de las andaluzas, ciertas zonas de la Cuenca del Ebro (aragonesa, navarra y riojana), del arco periurbano de Valencia, del Sur de Madrid, de la Mancha toledana, o de las provincias de Albacete, Guadalajara, Murcia e incluso de Palencia. Los núcleos antiguos se excavaron -en razón de su origen dedicado a defensa y refugio- en lugares poco accesibles, que mayoritariamente se dan en relieves anfractuosos, en las orlas marginales detríticas y carbonatadas de las cuencas terciarias y en terrenos más antiguos de la geológicamente denominada Cuenca Vasco Cantábrica, en las cuales florecieron centros de eremitismo del primitivo condado de Castilla. También son lugares inaccesibles los riscos asomados a cantiles fluviales, cuya regularización morfológica natural ha sido causa de la ruina de múltiples hipogeos que se labraron con las mismas intenciones de refugio, defensa o retiro espiritual. Los núcleos modernos se han excavado ya en terrenos más propicios (los que componen las cuencas terciarias o "España arcillosa", mayoritariamente), y por ello observamos cómo abundan las litologías sedimentarias, que son aquellas en las que se horadaron casi todas las "colmenas" de habitación moderna en la Península. En unos casos y otros, existen rasgos comunes en lo relativo a la litología y comportamiento de los materiales excavados, y también en lo relativo a su evolución, meteorización y conservación. Se han estudiado por ello estas pautas comunes -como un posible avance para el establecimiento futuro de estudios de geoconservación del patrimonio rupestre habitado-, que se traducen en la determinación de los procesos de meteorización más característicos en los antros de la Península Ibérica -sean modernos o antiguos- y en la determinación de las relaciones más habituales entre geomorfología y tipología de las cuevas de habitación. También se exponen algunas conclusiones relativas a la resistencia de los tipos pétreos en relación a la antigüedad de los emplazamientos. Esta relación se explica bien si tenemos en cuenta que las cuevas antiguas trataban de ser lugares apartados en los que dominan rocas más resistentes -como se ha explicado- y que en las modernas se ha buscado la habitación permanente, en los materiales más blandos ocupan las depresiones terciarias del Centro, Este y Sur de la Península, colonizados no ya como refugios sino de modo seguro, y ya en arrabales "extramuros" de las ciudades. Geomorfológicamente, esta razón histórica tiene consecuencias sobre la posición de las cuevas, sobre su organización, y sobre su conservación. La extensión del dominio estudiado obliga prácticamente a bosquejar algunas de estas conclusiones geológicas, pero permite a su vez proporcionar una visión global acerca del patrimonio troglodítico desde una perspectiva geológica, y en ello radica la principal novedad de la investigación. Cave dwelling in the Iberian Peninsula comprises a great deal of examples whose origins are sometimes difficult to elucidate. There are primitive groups of caves, probably belonging to Classical and Late Antiquity, and other settlements that appear to have been created and excavated in the Early Middle Ages. Some of them are due to religious reasons and some others may probably have served for defensive uses. The Peninsular group is very likely the most diverse in Western Europe, for it was not only influenced by Arabs, but it also assimilated uses from the Antique Christendom, imported from the Middle East. In this sense, Iberian cave dwellings connect with those of the Mediterranean area. There is not total certainty about the survival of caves serving for dwelling throughout the Late Middle Ages, once the Reconquista was a fact. However, underground excavations for human habitation were once again dispread in certain zones of the Peninsula during the Early Modern period. This growth of underground habitats appears in connexion with Morisco’s vicissitudes; first of all, their removal from their settlements in South and Eastern Peninsula, and finally their expulsion from Spanish territory at the very beginning of XVII th century. Those of them who rested in Spain –and that were nominally “converts”- seem to have returned to cave dwelling in their places of origin, particularly in the province of Granada and its neighbouring zones. This substrate may have been the toehold for a new spreading of cave dwellings since XIX th century, when the increasing rural and urban proletariat returned to the pre-existent uses of caves in order to solve the necessity of affordable housing. This fact can still be detected in many of the rock settlements that have survived during a great part of the XX th century and which were gradually abandoned from the 60´s onwards. There were important cave dwelling districts in many Spanish territories by that time, and not only in Andalusia. We also find them in certain areas of the Ebro basin (those of Aragon, Navarra and La Rioja), in the peri-urban arc of Valencia, in the South of Madrid province and also in the provinces of Toledo (the so called “Mancha toledana”), Albacete, Guadalajara, Murcia or even Palencia as well. Due to their defensive and refuge uses, primitive underground habitats were dug in hardly accessible places. The majority of them are located in mountainous and rough areas, when not in the marginal borders of Tertiary basins, where coarse detritic and carbonate formations outcrop. Cave dwellings can also be found in more ancient rock masses, such as those of the Basque Cantabrian Mesozoic Basin, which is the area where hermit centres of the primitive County of Castile first flourished. Cliffs surrounding fluvial valleys are as well inaccessible places, but here we find that geological evolution has caused the destruction of many rock sanctuaries and cliff dwellings that were originally dug with the same purposes of defence, refuge and spiritual retreat. Later modern cave settlements were dug in quite more favourable terrains, mainly in the soils that compose the Tertiary basins, generally known as “España arcillosa” or “Clayey Spain”. Therefore, we find abundant sedimentary fine and medium grained lithologies, which are the ones that have hosted the majority of Modern Era warren cave dwellings in the Iberian Peninsula. Actually, both types of cave dwelling share some standards regarding the lithology of the excavated materials, and they share as well certain patterns that affect to their evolution, weathering and preservation. These common patterns have been studied here in order to determine the most characteristic weathering processes that affect the majority of the Iberian caves, both Antique and Modern. And also with a view to establish the most habitual relationships between geomorphology and typology of cave dwellings. The study may as well provide a first basis for future studies on geo-preservation of cave dwellings heritage. We also reach some conclusions about the strength of different rocks concerning the antiquity of the sites. As we have already pointed out, this relation comes from the fact that ancient caves were placed in remote or isolated locations, where harder rocks outcrop, while more modern ones result from people´s search of permanent dwelling. In this sense, the softer rocks of the Tertiary Basins of Middle, East and South Peninsula provided a secure colonization to this second and modern group. And moreover, considering geomorphological features, this historical reason has had an effect not only on the position and location of the caves, but also on their organisation/structure/distribution and preservation. The huge extension of our domain of interest almost forces to sketch out some of these geological conclusions. But at the same time it gives a global panorama of Spanish troglodyte heritage, seen from a geological perspective. And here is the main novelty of this research.

MSB: Una aplicacion WWW para generación de resúmenes de comportamiento

La disponibilidad en Internet de grandes volúmenes de información en los que se basan decisiones relacionadas con el funcionamiento de sistemas dinámicos hace que cada vez sean más necesarios métodos efectivos de generación automáti-ca de resumen de información de comportamiento. En este artículo se describe la aplicación MSB (Multimedia Summari-zer of Behavior) como una herramienta software que genera presentaciones multimedia de resumen de comportamiento de sistemas dinámicos para ser presentadas en la Web. MSB usa como entrada datos de comportamiento de fuentes remo-tas (medidos por sensores) que pueden estar disponibles en diferentes servidores Web. Estas fuentes de información se interpretan y resumen usando métodos del campo de la Inteligencia Artificial simulando formas de razonamiento que realizan las personas al resumir información de interés. En el artículo se describe la arquitectura basada en el conocimien-to del sistema MSB y la aplicación desarrollada en el dominio de hidrología sobre comportamiento de cuencas hidrográ-ficas.

Espacios del agua en el territorio urbanizado. Los Angeles, California

Los Ángeles, en California, muestra de forma excepcional la hidrología de las grandes metrópolis en zonas de paisaje árido, basadas en grandes obras de ingeniería civil que concentran en la ciudad el agua de ámbitos geográficos muy extensos mediante aportes elevados de energía, al tiempo que canalizan la evacuación como residuo de los recursos locales. La desertización creciente derivada de este modelo sugiere la necesidad de alternativas de baja energía, basadas en la operatividad de la forma y de los sistemas vivos, para devolver al espacio y al territorio competencias de depuración, almacenaje y control de inundaciones. La construcción del suelo como infraestructura extensiva de captación, filtrado y reutilización modifica el funcionamiento hidrológico de la ciudad, reduciendo deslizamiento y escorrentía para conservar el valioso patrimonio de lluvia local y convierte la fábrica urbana en un gran sistema polivalente de gestión hidráulica. La recuperación de las cuencas locales de menor escala devuelve al agua y sus espacios su papel principal como articuladores de la forma urbana, marcando la localización de áreas húmedas y verdes. La naturaleza urbana se aleja del referente decorativo de pradera inglesa a favor de una concepción más utilitarista, entendida como una máquina orgánica integrada que resuelve la depuración, contribuye a la construcción de suelo y a la climatización a nivel metropolitano.

Estudio de la alteración hidrológica actual en el río Eresma por el embalse del Pontón Alto y de la alteración futura por la construcción del nuevo embalse del Ciguinuela.

El trabajo que se presenta se encuadra en el área de los trabajos científico-técnicos relacionados con la evaluación de impactos. En concreto, se incluye en el campo de la valoración de la alteración hidrológica en ríos. La zona de estudio en la que se centra el trabajo se localiza en la provincia de Segovia (Castilla y León). Se estudia el tramo del río Eresma aguas abajo de la Presa del Pontón Alto hasta la confluencia con el río Ciguiñuela. En su recorrido el río deja atrás la Granja de San Ildefonso y pasa por los municipios de Palazuelos del Eresma y San Cristóbal de Segovia hasta llegar a la propia ciudad de Segovia. El motivo de la redacción de este estudio se basa en que La Sociedad Estatal de Aguas de las Cuencas de España, S.A. ha planteado la construcción de una nueva presa en el río Ciguiñuela, afluente del río Eresma, que sería llenada con aguas procedentes del río Eresma derivadas a través de una conducción que uniría el futuro embalse del Ciguiñuela con el actualmente existente embalse del Pontón Alto. Este trabajo pretende profundizar en los impactos hidrológicos que se producirían en el río Eresma tras la creación de esta nueva obra.

El programa ERHIN y la gestión de embalses en áreas con alta precipitación nival = The ERHIN programme and reservoir management in areas of high nival precipitation

El programa de Evaluación Recursos Hídricos procedentes de la Innivación (ERHIN) ha desarrollado un modelo hidrológico con subrutinas que contemplan el proceso de acumulación/fusión de nieve que permite el seguimiento del manto de nieve y su fusión , esencial para la previsión de avenidas en embalses cuyas cuencas vertientes recogen altos porcentajes de precipitación solida. En determinados sectores del territorio español (Pirineos, Cordilleras Cantábrica e Ibérica, Sistema Central y Sierra Nevada) la precipitación nival alcanza valores significativos durante el invierno, que conducen a espesores de nieve persistentes en las cuencas alimentadoras de numerosos embalses y a la fusión rápida de esta nieve en períodos cortos de tiempo, directamente relacionados con la elevación de la isoterma de cero grados en primavera. Estos procesos condicionan la gestión de los mencionados embalses.

Evaluación de la protección conseguida en la cuenca de un torrente de montaña ante los daños causados por eventos torrenciales mediante su restauración hidrológico-forestal. Cuenca del torrente de Arratiecho en el Pirineo Aragonés, España.

La demanda continua e insistente de innovación en el ámbito cotidiano, nos hace perder la perspectiva de análisis cuando los procesos requieren un periodo prolongado para resolverse. Esto ocurre con los bosques ubicados en cuencas vertientes a torrentes de montaña, que tienen su origen en la restauración hidrológico-forestal llevada a cabo en las mismas, como complemento de los trabajos de corrección hidráulica efectuados en los torrentes, para controlar los procesos de erosión, inundaciones y aterramientos, desencadenados en diferentes zonas de la cuenca, cuando se dinamizan en ella los ciclos del agua y de los sedimentos a causa de los aguaceros o de la repentina fusión del manto de nieve. La cuenca vertiente al torrente de Arratiecho, situada en el municipio de Biescas, Pirineo Aragonés, fue restaurada y el barranco corregido entre 1903-05; para evitar su desertificación ante la erosión generalizada que sufría la misma, frenar el deslizamiento de una ladera morrénica y reducir los efectos de las inundaciones y aterramientos en su cono de sedimentación, donde se ubican los terrenos más fértiles y le atraviesa la carretera nacional N-260. Se comentan los trabajos realizados y los resultados obtenidos un siglo después, atendiendo especialmente a la evolución de las reforestaciones, siguiendo los criterios utilizados en el análisis de los bosques protectores.

Evaluación de modelos multivariados para la simulación estocástica de lluvia a escala diaria

La obtención de leyes de frecuencia de caudales de avenida mediante simulación hidrometeorológica requiere de la extensión de las series observadas de lluvia mediante series sintéticas que conserven sus propiedades extremales. Para cuencas pequeñas, cuya lluvia media pueda asimilarse a la registrada en un único pluviómetro, existen varios modelos puntuales que pueden realizar eficazmente la tarea. En el caso de cuencas medianas y grandes, deben utilizarse modelos multivariados, cuya calibración es más difícil que en el caso univariado. El objetivo de este trabajo es la simulación estocástica de lluvias a escala a través del empleo de un modelo multivariado en el tiempo y en el espacio. Para ello se ha empleado el programa Rainsim, que se basa en el modelo espacio-temporal de Neyman-Scott de pulsos rectangulares (e.g. Coperwait et al., 2002; Fowler et al., 2005; Burton et al., 2008). A partir de series históricas de lluvias se obtienen sus estadísticos, con los que se ajustan los parámetros del modelo para poder generar largas series de precipitaciones. En el trabajo se ha empleado el modelo Rainsim en un caso univariado y otro multivariado con tres pluviómetros, y se han comparado los resultados. El caso univariado es el embalse de Alloz, situado en la cuenca del Ebro, que cuenta con un pluviómetro con una serie histórica de 43 años con paso diario. En el caso multivariado se ha trabajado con el embalse de Fuensanta, situado en la cuenca del Segura, donde se escogieron tres pluviómetros (Gontar, Arguellite y Siles), que cuentan con series de datos de precipitación diaria de 31, 42 y 22 años. En los trabajos realizados se han obtenido resultados muy satisfactorios para el modelo puntual de un pluviómetro estudiado en la cuenca de Alloz, el comparar el comportamiento extremal de los datos observados y los datos simulados. En cuanto al modelo multivariado, los datos simulados presentan un comportamiento no tan acertado como en el caso univariado.

Aplicación del análisis de evaluación de riesgo a una cuenca hidrográfica en Panamá

Fueron evaluados los problemas relativos al manejo de los recursos hídricos y la gestión del agua en Panamá. A partir del comportamiento del clima, de la evaluación físico química de las principales cuencas hidrográficas de Panamá y del uso del suelo y del agua se aplicó el análisis de evaluación de riesgo.

Los nuevos regímenes ecológicos de caudales como herramienta para la mejora ambiental de los ecosistemas acuáticos y riparios.

Dos estrategias fundamentales para una gestión correcta de los regímenes fluviales, en el ámbito de las exigencias impuestas por la Directiva Marco del Agua (2000/60/CE) y por la propia Directiva Hábitat (92/43/CE) son el análisis de la alteración hidrológica de las masas de agua, y el establecimiento de un régimen ambiental de caudales que permita alcanzar un estado ecológico adecuado a los objetivos medioambientales de cada una de esas masas. Por lo que respecta al primero de estos aspectos, el análisis de la alteración hidrológica de las masas superficiales, existen diferentes metodologías para su cálculo e interpretación, entre las que destaca el modelo español IAHRIS, basado en la utilización de Indicadores de Alteración Hidrológica. En cuanto al establecimiento de caudales ambientales, las últimas modificaciones de la Ley de Aguas española (RDL 1/2001) y del Reglamento de Planificación Hidrológica (R.D. 907/2007), que incluyen cambios sustanciales en los procedimientos de cálculo de los requerimientos hídricos mínimos de los ecosistemas fluviales, ha motivado un cambio de estrategia en cuanto a su cálculo, concertación y seguimiento adaptativo. El Ministerio de Medio Ambiente y Medio Rural y Marino ha desarrollado a lo largo de los últimos años, con el apoyo técnico del CEDEX en la coordinación y desarrollo de los trabajos, un amplio conjunto de estudios técnicos destinados a incorporar todos estos requerimientos normativos en los Planes de Gestión de Cuenca. Se han realizado estudios de detalle en más de 400 masas de agua, correspondientes a las cuencas hidrográficas intercomunitarias. Los resultados de los estudios permiten extraer importantes consecuencias científicas y de gestión sobre la compatibilización de la protección ambiental y el aprovechamiento de los recursos hídricos nacionales, cuyas implicaciones en la conservación y restauración de los ecosistemas acuáticos y riparios resulta del máximo interés para la gestión de las masas de agua españolas. Los resultados muestran la necesidad de adaptar las metodologías existentes en el ámbito internacional a las singularidades hidrológicas, geomorfológicas y ecológicas de los ríos españoles, con el fin de facilitar su integración con los objetivos medioambientales de la Directiva Marco del Agua, y de asegurar una protección efectiva de los ecosistemas fluviales. Es preciso, asimismo, de cara a los siguientes ciclos de planificación, incrementar el conocimiento sobre los requerimientos hídricos de las comunidades animales y vegetales, los efectos sobre ellas de diferentes escenarios de gestión, y los mecanismos de protección de los ecosistemas fluviales en condiciones de sequías prolongadas.

Modelos hidráulicos bidimensionales de lluvia directa aplicados al cálculo de eventos hidrológicos

Una de las principales preocupaciones a la hora de realizar cualquier proyecto de ingeniería es la adecuada determinación de los hidrogramas de avenida de los cursos de agua que lo afecten. Hoy en día es posible la simulación de modelos hidrológicos complejos de lluvia directa sobre malla 2D de elementos finitos, empleando los avances en cálculo vectorial que ofrece la tecnología CUDA (Arquitectura Unificada de Dispositivos de Cálculo), sin la necesidad de determinar subcuencas y tránsitos (como en HEC-HMS). Con los datos que ofrecen los Modelos Digitales del Terreno GRID del PNOA español (Plan Nacional de Ortofotografía Aérea) con precisión de 5 x 5 m y las Coberturas del Terreno SIG CORINE (Coordinación de Información del Medio Ambiente), que permite la evaluación de la rugosidad del suelo, se tiene la información necesaria para construir fácilmente este tipo de modelos. En este documento se presenta un procedimiento para delimitar fácilmente los cursos de agua principales y disponer en estas zonas un tamaño más pequeño de los elementos de la malla. El método propuesto permite una mejor definición de estas áreas después de unas pocas iteraciones, acelerando por tanto el proceso de construcción del modelo sin comprometer la calidad de los resultados. Para demostrar el método, se modelizaron dos cuencas completas mediante lluvia directa sobre una malla de elementos 2D. Se ensayaron diferentes niveles de precisión (variando el número de elementos) y se incluyeron dos presas, usando los datos del Modelo Digital del Terreno del PNOA. También se completó el cálculo hidrológico de ambas mediante HEC-HMS. En primer lugar se estudió la cuenca del río Zapardiel (cuenca del Duero, con superficie total de 1.450 km2, no presenta regulación artificial, pero sí algunas áreas de almacenamiento naturales). La segunda cuenca ensayada correspondió al río Zújar (cuenca del Guadiana, con superficie total de 8.500 km2, que incluye dos grandes presas: La Serena, con una capacidad de 3.200 hm3 y el Embalse del Zújar, situado aguas abajo de la primera, con un volumen útil de 300 hm3). Los resultados confirmaron que la forma de modelizar clásica con enfoque de subcuencas y tránsitos (tipo HEC-HMS) puede ser reemplazada por simulaciones de lluvia directa sobre malla 2D, por lo tanto reduciendo la incertidumbre inherente en la estimación de tiempos de concentración de subcuencas y parámetros de tránsito. Gracias a la disponibilidad de datos públicos de GRID PNOA y cobertura del terreno CORINE, actualmente es posible construir modelos de mallas de elementos 2D para simular lluvia directa sobre cuenca sin necesidad de costosos levantamientos topográficos. Con el hardware disponible hoy en día en ordenadores de uso general, se puede establecer un límite superior razonable para mallas 2D de 4 millones de elementos, lo que permite la correcta simulación de la lluvia directa en cuencas de hasta 10.000 km2. Para cuencas más grandes, este proceso se podría aplicar de forma repetida a varias zonas más pequeñas o en sucesivas secciones de la cuenca.

Utilización de modelos hidráulicos bidimensionales en la determinación del tiempo de concentración

El tiempo de concentración de una cuenca sigue siendo relativamente desconocido para los ingenieros. El procedimiento habitual en un estudio hidrológico es calcularlo según varias fórmulas escogidas entre las existentes para después emplear el valor medio obtenido. De esta media se derivan los demás resultados hidrológicos, resultados que influirán en el futuro dimensionamiento de las infraestructuras. Este trabajo de investigación comenzó con el deseo de conseguir un método más fiable y objetivo que permitiera obtener el tiempo de concentración. Dada la imposibilidad de poner en práctica ensayos hidrológicos en una cuenca física real, ya que no resulta viable monitorizar perfectamente la precipitación ni los caudales de salida, se planteó llevar a cabo los ensayos de forma simulada, con el empleo de modelos hidráulicos bidimensionales de lluvia directa sobre malla 2D de volúmenes finitos. De entre todos los disponibles, se escogió InfoWorks ICM, por su rapidez y facilidad de uso. En una primera fase se efectuó la validación del modelo hidráulico elegido, contrastando los resultados de varias simulaciones con la formulación analítica existente. Posteriormente, se comprobaron los valores de los tiempos de concentración obtenidos con las expresiones referenciadas en la bibliografía, consiguiéndose resultados muy satisfactorios. Una vez verificado, se ejecutaron 690 simulaciones de cuencas tanto naturales como sintéticas, incorporando variaciones de área, pendiente, rugosidad, intensidad y duración de las precipitaciones, a fin de obtener sus tiempos de concentración y retardo. Esta labor se realizó con ayuda de la aceleración del cálculo vectorial que ofrece la tecnología CUDA (Arquitectura Unificada de Dispositivos de Cálculo). Basándose en el análisis dimensional, se agruparon los resultados del tiempo de concentración en monomios adimensionales. Utilizando regresión lineal múltiple, se obtuvo una nueva formulación para el tiempo de concentración. La nueva expresión se contrastó con las formulaciones clásicas, habiéndose obtenido resultados equivalentes. Con la exposición de esta nueva metodología se pretende ayudar al ingeniero en la realización de estudios hidrológicos. Primero porque proporciona datos de manera sencilla y objetiva que se pueden emplear en modelos globales como HEC-HMS. Y segundo porque en sí misma se ha comprobado como una alternativa realmente válida a la metodología hidrológica habitual. Time of concentration remains still fairly imprecise to engineers. A normal hydrological study goes through several formulae, obtaining concentration time as the median value. Most of the remaining hydrologic results will be derived from this value. Those results will determine how future infrastructures will be designed. This research began with the aim to acquire a more reliable and objective method to estimate concentration times. Given the impossibility of carrying out hydrological tests in a real watershed, due to the difficulties related to accurate monitoring of rainfall and derived outflows, a model-based approach was proposed using bidimensional hydraulic simulations of direct rainfall over a 2D finite-volume mesh. Amongst all of the available software packages, InfoWorks ICM was chosen for its speed and ease of use. As a preliminary phase, the selected hydraulic model was validated, checking the outcomes of several simulations over existing analytical formulae. Next, concentration time values were compared to those resulting from expressions referenced in the technical literature. They proved highly satisfactory. Once the model was properly verified, 690 simulations of both natural and synthetic basins were performed, incorporating variations of area, slope, roughness, intensity and duration of rainfall, in order to obtain their concentration and lag times. This job was carried out in a reasonable time lapse with the aid of the parallel computing platform technology CUDA (Compute Unified Device Architecture). Performing dimensional analysis, concentration time results were isolated in dimensionless monomials. Afterwards, a new formulation for the time of concentration was obtained using multiple linear regression. This new expression was checked against classical formulations, obtaining equivalent results. The publication of this new methodology intends to further assist the engineer while carrying out hydrological studies. It is effective to provide global parameters that will feed global models as HEC-HMS on a simple and objective way. It has also been proven as a solid alternative to usual hydrology methodology.