Análisis estructural de modelos hidrológicos y de sistemas de recursos hídricos en zonas semiáridas

espanolAnalisis estructural de modelos hidrologicos y de sistemas de recursos hidricos en zonas semiaridas RESUMEN El desarrollo de la modelacion hidrologica ha experimentado un impulso extraordinario en las ultimas decadas. Desde los primeros modelos conceptuales de la decada de los 60 del pasado siglo hasta hoy no solo el numero de modelos desarrollados es incontable, sino que sus tipologias se han multiplicado, sus aplicaciones practicas se han extendido, y la propia teoria subyacente a la modelacion ha sido objeto de avances sustanciales. En esta tesis nos centraremos, en primer lugar, en un tipo especifico de esos modelos, los denominados modelos conceptuales de balance hidrico o de lluvia-escorrentia, que simulan de forma agregada la respuesta hidrologica de una cuenca a partir de las precipitaciones producidas, con el objetivo de realizar una aplicacion masiva a numerosas cuencas espanolas y estudiar su comportamiento estructural. Para ello, y tras una presentacion general del ciclo hidrologico y sus procesos subyacentes desde la perspectiva de estos modelos, se describen algunos comunmente utilizados mostrando en detalle sus principios basicos, ecuaciones descriptivas, parametrizacion empleada y condiciones de aplicacion. Seguidamente se aborda el problema del analisis estructural, entendiendo por tal los aspectos relativos a la sensibilidad de los modelos, analisis de datos y perturbaciones, incertidumbres, problemas de identificabilidad, y validacion de su funcionamiento. La descripcion es autocontenida e incluye los desarrollos matematicos necesarios para la comprension de los conceptos y tecnicas introducidos. Para poder aplicar estas tecnicas se ha desarrollado una aplicacion informatica, que incorpora muchos de los modelos y tecnicas descritas y permite aplicarlas a gran escala. Brevemente se describe esta aplicacion, sus interfaces y los resultados proporcionados. Con este codigo se ha procedido a la aplicacion practica de la teoria a un extenso conjunto de cuencas espanolas peninsulares, centrando la perspectiva en las singularidades y diferencias de los modelos, desde el punto de vista estructural, segun el grado de aridez de las cuencas a las que se aplican. Se introducen los conceptos de regionalizacion y estabilidad estructurales dependientes de indicadores, adoptando como tal indicador la aridez de las cuencas. Las aplicaciones practicas se completan con el estudio de perturbaciones aplicado a la evaluacion de efectos del cambio climatico mediante modelos conceptuales agregados. Tras estos analisis, puramente hidrologicos, la investigacion se completa con una consideracion teorica del posible efecto de la aridez y torrencialidad de las cuencas sobre el funcionamiento de los modelos de optimizacion y simulacion aplicados a los sistemas de recursos hidricos, aplicandolo a un caso real concreto. Por ultimo, se concluye con algunas consideraciones y posibles extensiones de las tecnicas mostradas, orientadas a la modelacion semidistribuida. El analisis ha mostrado la viabilidad de la aplicacion practica de estas metodologias, la posibilidad de mejorar la parametrizacion de los modelos, y que muchos resultados estructurales relativos a modelos de amplio uso son, en efecto, dependientes de la aridez de las cuencas a las que se aplican. La aridez no solo condiciona el valor esperable y rango admisible de sus parametros, informacion de utilidad para acotar informacion requerida a priori, sino que tambien influye en sus niveles de correlacion y, por tanto, su identificabilidad estructural. Asimismo se ha verificado que la aridez o irregularidad de flujos un territorio puede afectar a los resultados de los modelos de optimizacion y simulacion de sistemas de recursos hidricos, al invalidar el supuesto de paso y captacion mensual completa asumido por estos modelos. La importancia de esta afeccion depende de la configuracion de estos sistemas y las caracteristicas de sus caudales a escalas inframensuales como la diaria. Con igual objetivo de valorar efectos de la aridez, se ha analizado el efecto de la concentracion subdiaria de las intensidades de lluvia sobre la infiltracion a escala diaria, observando que el efecto puede ser significativo dependiendo de la textura del suelo y su humedad inicial. Finalmente se han apuntado algunos criterios o bases para la extension de las tecnicas mostradas a otras tipologias de modelos como los diarios o mensuales semidistribuidos, linea de gran interes practico al permitir realizar evaluaciones de recursos a gran escala, en extensos territorios con numerosas estaciones, y contrastadas mediante calibracion con datos observados en subcuencas aforadas, cuya extrapolacion puede llevarse a cabo con las tecnicas descritas. EnglishABSTRACT The development of hydrological modeling has experienced an extraordinary boost in recent decades. From the first conceptual models of the 60s of the last century until today, not only the number of developed models is countless, but their typologies have multiplied, their practical applications have been extended, and underlying theory to modeling has been the subject of substantial progress. This thesis will focus, firstly, in a specific type of those models, the so-called rainfall-runoff conceptual models of continuous water balance, simulating the hydrologic response of a basin, in order to carry out a massive application of several models to numerous Spanish basins, and study their structural behaviour. To do this, and after a general presentation of the hydrological cycle and its underlying processes from the perspective of these models, some of them commonly used are described in detail showing its basic principles, describing equations, parameterization used and application conditions. Then, the problem of structural analysis is discussed, including sensitivity of the models, data analysis and perturbations, uncertainties, identifiability problems, and validation of its operation. The description is self-contained and includes the mathematical developments needed for an understanding of the introduced concepts and techniques. To apply these techniques we developed a software that incorporates many of the models and techniques described, and allows a large-scale execution. The software interfaces and the results provided are briefly described. Using this software we proceeded to the practical application of the theory to an extensive set of Spanish basins, focusing perspective in singularities and differences of models, from a structural point of view, according to the degree of aridity of the considered basins. The concepts of regionalization and structural stability dependence on indicators are introduced, adopting as indicator the basin aridity. The practical applications are completed with the study through conceptual models of perturbations applied to the evaluation of effects of climate change on water resources. Following these purely hydrological issues, research is completed with a theoretical consideration of the possible effect of highly irregulars daily flow on the results of optimization and simulation models for water resources systems, with a monthly time step, and applying it to real case. Finally, we conclude with some considerations and possible extensions of the exposed techniques to semidistributed and daily new model under development. The analysis has shown the feasibility of the practical application of these methodologies, the possibility of improving the parametrization of the models, and many structural results for widely used models (structural identification matrix, covariances of parameter, eigenvalues, conditioned sensitivity, sensitivity ratio, etc.) that are in fact dependent on the aridity of the basins. Aridity not only determines the expected value and permissible range of parameters, useful information to narrow the a priori information required, but also it affects their levels of correlation and therefore structural identifiability. It has also been verified that the irregularity of daily flows, typical in semiarid areas, can affect the results of optimization and simulation models of water resources systems, due to invalidate the assumption of full monthly step diversion taken by these models. The importance of this effect depends on the configuration of the system and the flow characteristics at submonthly scales as daily. With the same objective of evaluating the effects of aridity, we analyzed the effect of intraday concentration of rainfall intensities on a daily scale infiltration, noting that the effect can be significant depending on the soil texture and initial water content. Finally we include some considerations for extend the techniques to other types of models such as daily or monthly semi-distributed using tessellation, line of great practical interest to enable water resource assessments on a large scale, in extensive territories with numerous gage stations. Model results can be contrasted by calibration with gauged watersheds, and parameter extrapolation can be carried out with the described techniques.