Proyecto de corrección y repoblación forestal de la Cuenca de la Rambla de la Hoya de los Álamos (provincia de Almería)

La Sección objeto de este Proyecto, comprendiendo también la llamada "Rambla del Gallo”, que de cuenca muy reducida, corre paralela a la de la Hoya de los Álamos, se encuentra toda ella en los términos municipales de Fondón y Fuente- Victoria, correspondiendo a éste último toda la cuenca de la Rambla de la Hoya de los Álamos y una buena parte de la del Gallo, y comprendiendose en el término de Fondón silo la parte restante de esta. En resuman, de las 1637,50 has, que totalizan ambas cuencas,1.380 has. son de Fuente Victoria y el resto de Fondón, Términos ambos del partido judicial de Canjáyar. Esta situación administrativa puede considerarse como estable, por estar dichos términos desandados por el Instituto Geográfico

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.

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.

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.

Metodología para la gestión integral de las presas de laminación de avenidas y su aplicación a la cuenca hidrográfica del Segura

En el levante español se producen periódicamente intensos aguaceros (gota fría), que generan grandes riadas, desbordándose los cauces y asolando los campos y las ciudades, ocasionando graves daños en la economía local y con frecuencia, por desgracia, víctimas personales. La Confederación Hidrográfica del Segura ha puesto en marcha un vasto sistema de protección, conocido como Plan de Defensa contra Avenidas, del que ya se encuentran en explotación un conjunto de presas cuyo objetivo prioritario es la laminación de las avenidas. Las presas de laminación ubicadas en esta área mediterránea son estructuras muy singulares, al estar ubicadas en cuencas extremadamente secas, con alta pendiente, rápida respuesta, picos altos de caudal, y gran poder erosivo del flujo. Las aguas arrastran en suspensión un alto porcentaje de aporte sólido y flotantes, que quedan retenidos en el embalse, aterrándolo e inutilizando los desagües. En esta tesis doctoral se analizan los criterios de actuación para la gestión de estas presas de laminación, con objeto de que se puedan introducir en las Normas de Explotación innovaciones destinadas a mejorar la eficiencia del conjunto, de manera que se consiga la máxima protección posible y al mismo tiempo se puedan aprovechar las aguas, sin menoscabo de su objetivo prioritario, para otros fines complementarios. Como consecuencia del trabajo realizado se indica la conveniencia de establecer un conjunto de normas de actuación comunes a todas las presas del sistema, como es la estrategia de recibir siempre a los hidrogramas de la riada con las compuertas de los desagües de fondo cerrados, o la conveniencia de que en todas las presas se dispongan compuertas en estos órganos. También se dictan recomendaciones en la gestión del vaciado, para conseguir el máximo aprovechamiento del agua, introduciendo el concepto innovador de los trasvases virtuales. Asimismo, se analiza la conveniencia de fijar normas diferenciales específicas atendiendo a las peculiaridades de cada presa, en cuanto se refiere al nivel máximo de explotación, al cumplimiento de los condicionantes ambientales, al mantenimiento de la calidad del agua, y a otros aspectos singulares que caracterizan a cada uno de estos embalses. The Spanish Mediterranean region hydrology is characterised by low rainfall, but concentrated in extreme convective events. Those heavy rains and storms lead to large, flash-floods which cause devastation of fields, villages and, unfortunately, loss of lives. Segura River Basin Authority has developed an ambitious plan for flood protection and control that is constituted by a system of flood protection dams, river channelling, river diversions diversions and levees. Flood protection dams placed on this area, are distinct structures, as they are located on extremely dry basins, with high slopes and limited vegetal covering. These facts, along with the very high rainfall intensity, are determinants of a rapid response and a very erosive flow. The flow transports a large rate of solids in suspension, and also debris, which many times block and clog the outlets intakes, disabling them. In such a case, the dam would be inoperative, in an emergency state and remedial measures should be taken. This PhD thesis analyses the flood protection dams operation criteria, for providing managers a methodology to analyse each individual case and also the flood protection system as a whole. Additionally practical rules and recommendations, to cope with the objectives of dam safety, flood protection and water resources use, are given. The rules consider the specific differences and risks that may be faced by the dams. To that end, they were grouped, depending on their characteristics: allowance for permanent water storage, environmental restrictions, water quality, geotechnical problems, type of structure, etc. As a result of this research, it is recommended to design these dams with valves for regulating the bottom outlets. The availability of valves is linked with the routing strategy of receiving flood with the outlets closed. This has been shown as the best measure to prevent clogging and opens the possibility of storing temporarily the water, making then a scheduled emptying of the reservoir, which could be computed as a virtual transfer to the regulating dams. Emptying should be done in accordance with the situation of the other dams of the system and of the rivers downstream. This detention time would make possible also the beginning of the sedimentation process which is necessary for using the water with the modern irrigation systems that are in use in this area.

Desarrollo de metodologías para mejorar la estimación de los hidrogramas de diseño para el cálculo de los órganos de desagüe de las presas

En la presente Tesis se ha llevado a cabo el contraste y desarrollo de metodologías que permitan mejorar el cálculo de las avenidas de proyecto y extrema empleadas en el cálculo de la seguridad hidrológica de las presas. En primer lugar se ha abordado el tema del cálculo de las leyes de frecuencia de caudales máximos y su extrapolación a altos periodos de retorno. Esta cuestión es de gran relevancia, ya que la adopción de estándares de seguridad hidrológica para las presas cada vez más exigentes, implica la utilización de periodos de retorno de diseño muy elevados cuya estimación conlleva una gran incertidumbre. Es importante, en consecuencia incorporar al cálculo de los caudales de diseño todas la técnicas disponibles para reducir dicha incertidumbre. Asimismo, es importante hacer una buena selección del modelo estadístico (función de distribución y procedimiento de ajuste) de tal forma que se garantice tanto su capacidad para describir el comportamiento de la muestra, como para predecir de manera robusta los cuantiles de alto periodo de retorno. De esta forma, se han realizado estudios a escala nacional con el objetivo de determinar el esquema de regionalización que ofrece mejores resultados para las características hidrológicas de las cuencas españolas, respecto a los caudales máximos anuales, teniendo en cuenta el numero de datos disponibles. La metodología utilizada parte de la identificación de regiones homogéneas, cuyos límites se han determinado teniendo en cuenta las características fisiográficas y climáticas de las cuencas, y la variabilidad de sus estadísticos, comprobando posteriormente su homogeneidad. A continuación, se ha seleccionado el modelo estadístico de caudales máximos anuales con un mejor comportamiento en las distintas zonas de la España peninsular, tanto para describir los datos de la muestra como para extrapolar a los periodos de retorno más altos. El proceso de selección se ha basado, entre otras cosas, en la generación sintética de series de datos mediante simulaciones de Monte Carlo, y el análisis estadístico del conjunto de resultados obtenido a partir del ajuste de funciones de distribución a estas series bajo distintas hipótesis. Posteriormente, se ha abordado el tema de la relación caudal-volumen y la definición de los hidrogramas de diseño en base a la misma, cuestión que puede ser de gran importancia en el caso de presas con grandes volúmenes de embalse. Sin embargo, los procedimientos de cálculo hidrológico aplicados habitualmente no tienen en cuenta la dependencia estadística entre ambas variables. En esta Tesis se ha desarrollado un procedimiento para caracterizar dicha dependencia estadística de una manera sencilla y robusta, representando la función de distribución conjunta del caudal punta y el volumen en base a la función de distribución marginal del caudal punta y la función de distribución condicionada del volumen respecto al caudal. Esta última se determina mediante una función de distribución log-normal, aplicando un procedimiento de ajuste regional. Se propone su aplicación práctica a través de un procedimiento de cálculo probabilístico basado en la generación estocástica de un número elevado de hidrogramas. La aplicación a la seguridad hidrológica de las presas de este procedimiento requiere interpretar correctamente el concepto de periodo de retorno aplicado a variables hidrológicas bivariadas. Para ello, se realiza una propuesta de interpretación de dicho concepto. El periodo de retorno se entiende como el inverso de la probabilidad de superar un determinado nivel de embalse. Al relacionar este periodo de retorno con las variables hidrológicas, el hidrograma de diseño de la presa deja de ser un único hidrograma para convertirse en una familia de hidrogramas que generan un mismo nivel máximo en el embalse, representados mediante una curva en el plano caudal volumen. Esta familia de hidrogramas de diseño depende de la propia presa a diseñar, variando las curvas caudal-volumen en función, por ejemplo, del volumen de embalse o la longitud del aliviadero. El procedimiento propuesto se ilustra mediante su aplicación a dos casos de estudio. Finalmente, se ha abordado el tema del cálculo de las avenidas estacionales, cuestión fundamental a la hora de establecer la explotación de la presa, y que puede serlo también para estudiar la seguridad hidrológica de presas existentes. Sin embargo, el cálculo de estas avenidas es complejo y no está del todo claro hoy en día, y los procedimientos de cálculo habitualmente utilizados pueden presentar ciertos problemas. El cálculo en base al método estadístico de series parciales, o de máximos sobre un umbral, puede ser una alternativa válida que permite resolver esos problemas en aquellos casos en que la generación de las avenidas en las distintas estaciones se deba a un mismo tipo de evento. Se ha realizado un estudio con objeto de verificar si es adecuada en España la hipótesis de homogeneidad estadística de los datos de caudal de avenida correspondientes a distintas estaciones del año. Asimismo, se han analizado los periodos estacionales para los que es más apropiado realizar el estudio, cuestión de gran relevancia para garantizar que los resultados sean correctos, y se ha desarrollado un procedimiento sencillo para determinar el umbral de selección de los datos de tal manera que se garantice su independencia, una de las principales dificultades en la aplicación práctica de la técnica de las series parciales. Por otra parte, la aplicación practica de las leyes de frecuencia estacionales requiere interpretar correctamente el concepto de periodo de retorno para el caso estacional. Se propone un criterio para determinar los periodos de retorno estacionales de forma coherente con el periodo de retorno anual y con una distribución adecuada de la probabilidad entre las distintas estaciones. Por último, se expone un procedimiento para el cálculo de los caudales estacionales, ilustrándolo mediante su aplicación a un caso de estudio. The compare and develop of a methodology in order to improve the extreme flow estimation for dam hydrologic security has been developed. First, the work has been focused on the adjustment of maximum peak flows distribution functions from which to extrapolate values for high return periods. This has become a major issue as the adoption of stricter standards on dam hydrologic security involves estimation of high design return periods which entails great uncertainty. Accordingly, it is important to incorporate all available techniques for the estimation of design peak flows in order to reduce this uncertainty. Selection of the statistical model (distribution function and adjustment method) is also important since its ability to describe the sample and to make solid predictions for high return periods quantiles must be guaranteed. In order to provide practical application of previous methodologies, studies have been developed on a national scale with the aim of determining a regionalization scheme which features best results in terms of annual maximum peak flows for hydrologic characteristics of Spanish basins taking into account the length of available data. Applied methodology starts with the delimitation of regions taking into account basin’s physiographic and climatic characteristics and the variability of their statistical properties, and continues with their homogeneity testing. Then, a statistical model for maximum annual peak flows is selected with the best behaviour for the different regions in peninsular Spain in terms of describing sample data and making solid predictions for high return periods. This selection has been based, among others, on synthetic data series generation using Monte Carlo simulations and statistical analysis of results from distribution functions adjustment following different hypothesis. Secondly, the work has been focused on the analysis of the relationship between peak flow and volume and how to define design flood hydrographs based on this relationship which can be highly important for large volume reservoirs. However, commonly used hydrologic procedures do not take statistical dependence between these variables into account. A simple and sound method for statistical dependence characterization has been developed by the representation of a joint distribution function of maximum peak flow and volume which is based on marginal distribution function of peak flow and conditional distribution function of volume for a given peak flow. The last one is determined by a regional adjustment procedure of a log-normal distribution function. Practical application is proposed by a probabilistic estimation procedure based on stochastic generation of a large number of hydrographs. The use of this procedure for dam hydrologic security requires a proper interpretation of the return period concept applied to bivariate hydrologic data. A standard is proposed in which it is understood as the inverse of the probability of exceeding a determined reservoir level. When relating return period and hydrological variables the only design flood hydrograph changes into a family of hydrographs which generate the same maximum reservoir level and that are represented by a curve in the peak flow-volume two-dimensional space. This family of design flood hydrographs depends on the dam characteristics as for example reservoir volume or spillway length. Two study cases illustrate the application of the developed methodology. Finally, the work has been focused on the calculation of seasonal floods which are essential when determining the reservoir operation and which can be also fundamental in terms of analysing the hydrologic security of existing reservoirs. However, seasonal flood calculation is complex and nowadays it is not totally clear. Calculation procedures commonly used may present certain problems. Statistical partial duration series, or peaks over threshold method, can be an alternative approach for their calculation that allow to solve problems encountered when the same type of event is responsible of floods in different seasons. A study has been developed to verify the hypothesis of statistical homogeneity of peak flows for different seasons in Spain. Appropriate seasonal periods have been analyzed which is highly relevant to guarantee correct results. In addition, a simple procedure has been defined to determine data selection threshold on a way that ensures its independency which is one of the main difficulties in practical application of partial series. Moreover, practical application of seasonal frequency laws requires a correct interpretation of the concept of seasonal return period. A standard is proposed in order to determine seasonal return periods coherently with the annual return period and with an adequate seasonal probability distribution. Finally a methodology is proposed to calculate seasonal peak flows. A study case illustrates the application of the proposed methodology.