Modelización de servicios ecosistémicos en una cuenca hidrológica del área mediterránea = Modeling ecosystem services in a mediterranean area river basin

La modelización es un proceso por el que se obtienen modelos de los procesos del ´mundo real´ mediante la utilización de simplificaciones. Sin embargo, las estimaciones obtenidas con el modelo llevan implícitas incertidumbre que se debe evaluar. Mediante un análisis de sensibilidad se puede mejorar la confianza en los resultados, sin embargo, este paso a veces no se realiza debido básicamente al trabajo que lleva consigo este tipo de análisis. Además, al crear un modelo, hay que mantener un equilibrio entre la obtención de resultados lo más exactos posible mediante un modelo lo más sencillo posible. Por ello, una vez creado un modelo, es imprescindible comprobar si es necesario o no incluir más procesos que en un principio no se habían incluido. Los servicios ecosistémicos son los procesos mediante los cuales los ecosistemas mantienen y satisfacen el bienestar humano. La importancia que los servicios ecosistémicos y sus beneficios asociados tienen, junto con la necesidad de realizar una buena gestión de los mismos, han estimulado la aparición de modelos y herramientas para cuantificarlos. InVEST (Integrated Valuation of Ecosystem Services and Tradoffs) es una de estas herramientas específicas para calcular servicios eco-sistémicos, desarrollada por Natural Capital Project (Universidad de Stanford, EEUU). Como resultado del creciente interés en calcular los servicios eco-sistémicos, se prevé un incremento en la aplicación del InVEST. La investigación desarrollada en esta Tesis pretende ayudar en esas otras importantes fases necesarias después de la creación de un modelo, abarcando los dos siguientes trabajos. El primero es la aplicación de un análisis de sensibilidad al modelo en una cuenca concreta mediante la metodología más adecuada. El segundo es relativo a los procesos dentro de la corriente fluvial que actualmente no se incluyen en el modelo mediante la creación y aplicación de una metodología que estudiara el papel que juegan estos procesos en el modelo InVEST de retención de nutrientes en el área de estudio. Los resultados de esta Tesis contribuirán a comprender la incertidumbre involucrada en el proceso de modelado. También pondrá de manifiesto la necesidad de comprobar el comportamiento de un modelo antes de utilizarlo y en el momento de interpretar los resultados obtenidos. El trabajo en esta Tesis contribuirá a mejorar la plataforma InVEST, que es una herramienta importante en el ámbito de los servicios de los ecosistemas. Dicho trabajo beneficiará a los futuros usuarios de la herramienta, ya sean investigadores (en investigaciones futuras), o técnicos (en futuros trabajos de toma de decisiones o gestión ecosistemas). ABSTRACT Modeling is the process to idealize real-world situations through simplifications in order to obtain a model. However, model estimations lead to uncertainties that have to be evaluated formally. The role of the sensitivity analysis (SA) is to assign model output uncertainty based on the inputs and can increase confidence in model, however, it is often omitted in modelling, usually as a result of the growing effort it involves. In addition, the balance between accuracy and simplicity is not easy to assess. For this reason, when a model is developed, it is necessary to test it in order to understand its behavior and to include, if necessary, more complexity to get a better response. Ecosystem services are the conditions and processes through which natural ecosystems, and their constituent species, sustain and fulfill human life. The relevance of ecosystem services and the need to better manage them and their associated benefits have stimulated the emergence of models and tools to measure them. InVEST, Integrated Valuation of Ecosystem Services and Tradoffs, is one of these ecosystem services-specific tools developed by the Natural Capital Project (Stanford University, USA). As a result of the growing interest in measuring ecosystem services, the use of InVEST is anticipated to grow exponentially in the coming years. However, apart from model development, making a model involves other crucial stages such as its evaluation and application in order to validate estimations. The work developed in this thesis tries to help in this relevant and imperative phase of the modeling process, and does so in two different ways. The first one is to conduct a sensitivity analysis of the model, which consists in choosing and applying a methodology in an area and analyzing the results obtained. The second is related to the in-stream processes that are not modeled in the current model, and consists in creating and applying a methodology for testing the streams role in the InVEST nutrient retention model in a case study, analyzing the results obtained. The results of this Thesis will contribute to the understanding of the uncertainties involved in the modeling process. It will also illustrate the need to check the behavior of every model developed before putting them in production and illustrate the importance of understanding their behavior in terms of correctly interpreting the results obtained in light of uncertainty. The work in this thesis will contribute to improve the InVEST platform, which is an important tool in the field of ecosystem services. Such work will benefit future users, whether they are researchers (in their future research), or technicians (in their future work in ecosystem conservation or management decisions).

Efecto de la fertilización orgánica y mineral en las propiedades del suelo, la emisión de los principales gases de efecto invernadero y en las diferentes fases fenológicas del cultivo de café (Coffea arabica L.)

El cultivo de café es de gran importancia a nivel mundial (ICO, 2011), y en el Ecuador ha sido uno de los cultivos más importantes en la generación de divisas (COFENAC, 2011). Sin embargo en los sistemas productivos de este país se puede apreciar el uso inapropiado de fertilizantes, lo que conlleva a una pérdida de nutrientes, por lo que es importante estudiar las dosis adecuadas para la fertilización tanto mineral como orgánica. El objetivo del trabajo fue evaluar el efecto de la fertilización mineral y orgánica en diferentes dosis en un monocultivo de café en la provincia de Loja, sobre las propiedades del suelo, la emisión de los principales gases que provocan el efecto invernadero y la fenología y productividad del cultivo. En la provincia de Loja (Ecuador) se seleccionó un área de 2.520 m2 en la que se establecieron 21 parcelas de café arábigo (Coffea arabica L.) var. caturra y se aplicó tres tratamientos con tres repeticiones de fertilización mineral y tres orgánicos con dosis: bajas minerales (MIN 1= 157 Kg NPK ha-1 año-1 para el primer año y 425 Kg NPK ha-1 año-1 para el segundo año), medias minerales (MIN 2= 325 Kg NPK ha-1 año-1 para el primer año y 650 Kg NPK ha-1 año-1 en el segundo año) y altas minerales (MIN 3= 487 y 875 Kg NPK ha-1 año-1 para el primer y segundo año respectivamente), bajas orgánicas (ORG 1= 147 Kg NPK ha-1 año-1 en el primer año y 388 Kg NPK ha-1 año-1 en el año dos), medias orgánicas (ORG 2= 265 Kg NPK ha-1 año-1 para el primer año y 541 Kg NPK ha-1 año-1 en el segundo año), altas orgánicas (ORG 3= 368 Kg NPK ha-1 año-1 para el primer año y 727 Kg NPK ha-1 año-1 en el segundo año) y fertilización cero (TES = sin fertilización). Se usó urea, roca fosfórica y muriato de potasio en la fertilización mineral y humus (Bioabor) en la orgánica, más un tratamiento testigo, cada tratamiento tuvo tres repeticiones. El tiempo de evaluación de los fertilizantes aplicados fue de dos años consecutivos, la fertilización se la realizó dos veces por año y en base a análisis del suelo y demandas nutricionales del cultivo. para determinar las características del suelo se realizó muestreos de suelos en cada parcela a una profundidad de 20 cm de estas muestras los parámetro iniciales determinados fueron: color (Munsell), textura (método del hidrómetro), pH (relación 1:2,5 suelo-agua), Materia orgánica (Walkey y Black), Nitrógeno (Micro Kjendahl), Fósforo (Bray y Kurtz), Potasio (Olsen), estos procesos se repitieron cada seis meses para poder evaluar los cambios de que se producen debido a la fertilización mineral y orgánica en el cultivo. Las emisiones de gases efecto invernadero desde el suelo al ambiente se determinaron por el método de cámara cerrada (Rondón, 2000) y la concentración por cromatografía de gases. Las mediciones fisiológicas (altura de planta, ancho de copa, grosor de tallo y producción) se las evaluó cada dos meses, a excepción de la producción que fue anual al término de cada cosecha. Además se realizó el análisis económico de la productividad del cultivo. El análisis estadístico de datos se lo realizó con el programa SPSS v. 17.0. Las medias fueron comprobadas mediante ANOVAS de un factor con test de Tukey (P < 0,05). El beneficio económico se estimó en términos de ingresos y gastos totales que se presentaron en el ensayo. Los resultados obtenidos al término del ensayo indican que los tratamientos MIN 2 y MIN 3 produjeron cambios más significativos en comparación con los otros tratamientos establecidos en la mejora de fertilidad del suelo, el pH ha sido menos afectado en la acidificación en comparación con los tratamientos orgánicos que se han acidificado mayormente; la materia orgánica (MO) tuvo incrementos considerablemente bueno en estos dos tratamientos, sin embargo fueron superados por los tratamientos de fertilización orgánica; el nitrógeno total (Nt )y el potasio (K) también presentaron mejores valores al termino del ensayo y el fósforo (P) mostro incrementos buenos aunque un poco menores que los de los tratamientos ORG 2 y ORG 3. En lo que respecta a las emisiones de gases efecto invernadero, los flujos acumulados de óxido nitroso (N2O) en los dos años han aumentado en todos los tratamientos en comparación con el tratamiento Testigo, pero de manera considerable y con mayores flujos en el tratamiento MIN 3 y MIN 2 que se podrían considerarse los de mayor contaminación por N2O al ambiente lo que se le atribuye a las dosis de fertilización mineral aplicadas en el periodo de investigación, los tratamiento MIN 1 y todos los tratamientos orgánicos muestran menores emisiones al ambiente. Las emisiones de metano (CH4) no muestran mayores diferencias de emisiones entre tratamientos, siendo los mayores emisores los tratamientos ORG 3 y ORG 2 posiblemente debido al abono orgánico y añadido al suelo; para las emisiones de dióxido de carbono (CO2) de manera similar al CH4 el tratamiento ORG 3 fue el que presento mayores emisiones, los flujos de CO2 al ambiente de los otros tratamientos fueron menores y no presentaron diferencias significativas entre ellos. La variables fisiológicas en todos los casos apoyaron al desarrollo de las plantas de café, esto al ser comparadas con el tratamiento Testigo, sin embargo las que alcanzaron las mayores altitudes, anchos de copas y diámetro de tallo fueron las plantas del tratamiento MIN 3, seguido del MIN 3, no mostrando significancia entre ellos, y para los tratamientos orgánicos el que presento muy buenos resultados en estas variables ha sido el ORG 3, el cual no presento diferencias significativas con el MIN 2, lo cual comprueba que la fertilización mineral es más efectiva en este caso frente a la orgánica. Para el primer año de producción el tratamiento mineral con fertilización MIN 3 es el que obtuvo mayor producción no presentando diferencia estadística con el tratamiento con el MIN 2, no obstante fueron significativamente mayores que los otros tratamientos. Vale indicar que también el tratamiento MIN 1 y el tratamiento ORG 3 han presentado una producción considerable de café no mostrando diferencias estadísticas entre ellos. Para el segundo año la producción el cultivo mostró mayores rendimientos que el primer año de evaluación en todos los tratamientos, esto debido a la fisiología propia del cultivo y por otra parte se atribuye a la adición de fertilizantes que se ha realizado durante todo el ensayo; de manera similar al anterior los tratamientos MIN 3 y MIN 2 obtuvieron mejores rendimientos, no enseñando diferencias estadísticas significativas entre ellos, no obstante el tratamiento mineral dosis MEDIA no presentó significancia estadística con el ORG 3. El benéfico económico ha resultado mayor en el tratamiento MIN 3 y MIN 2, aunque el tratamiento MIN 2, es el que obtiene la mejor relación costo-beneficio; los tratamientos ORG 2 y ORG 3 y Testigo has producido beneficios negativos para el productor. En cuanto a la parte ambiental se considera que los mejores tratamientos en cuanto ha cuidado ambiental serían los tratamientos MIN 1 y ORG 1, sin embargo a nivel de producción y rentabilidad para el productor baja. ABSTRACT Coffee growing has great importance worldwide (ICO, 2011), and in Ecuador, it has been one of the most important crops to generate income (COFENAC, 2011). However, in the productive systems of this country, the inappropriate use of fertilizers has been observed which produces loss of nutrients, thus it is important to study suitable doses for mineral and organic fertilizing. The purpose of the study was to evaluate the effect of mineral and organic fertilizing at different doses in a coffee monoculture in the province of Loja on soil characteristics, emission of the main gasses that produce the greenhouse effect and the phenology and productivity of crops. In the province of Loja (Ecuador) an area of 2.520 m2 was chosen, where 21 plots of Arabica coffee (Coffea arabica L.), the caturra variety were cultivated and three treatments with three repetitions each one for mineral and organic fertilization were used with doses that ranged from: mineral low (MIN 1= 157 Kg NPK ha-1 año-1 for the first year y 425 Kg NPK ha-1 año-1 for the second year), mineral medium (MIN 2= 325 Kg NPK ha-1 año-1 for the first year y 650 Kg NPK ha-1 año-1 I the second year) y mineral high (MIN 3= 487 y 875 Kg NPK ha-1 año-1 for the first and second year respectively), organic low (ORG 1= 147 Kg NPK ha-1 año-1 in the first year y 388 Kg NPK ha-1 año-1 in the second year), organics medium (ORG 2= 265 Kg NPK ha-1 año-1 for the first year y 541 Kg NPK ha-1 año-1 in the second year), organics high (ORG 3= 368 Kg NPK ha-1 año-1 for the first year and 727 Kg NPK ha-1 año-1 in the second year) y fertilization zero (TES = no fertilization).; urea, phosphoric rock and muriate of potash were used in the mineral fertilization and humus (Bioabor) in the organic, plus a blank treatment. Time to evaluate the applied fertilizers was for two consecutive years, fertilization was done twice per year based on soil analysis and nutritional requirements of the crops. In order to determine the characteristics of the soil, samples of soil in each plot with a depth of 20 cm were done; from these samples, the determined initial parameters were: color (Munsell), texture (hydrometer method), pH (soil-water 1:2,5 relation), organic matter (Walkey y Black), nitrogen (Micro Kjendahl), phosphorus (Bray y Kurtz), potassium (Olsen); these processes were repeated each six months in order to evaluate the changes that are produced due to mineral and organic fertilization in the crops. The emissions of greenhouse gasses from the soil to the atmosphere were determined by using enclosure method (Rondón, 2000) and the concentration, by using gas chromatography during the whole testing. The physiological measures (plant height, width of the top of the tree, thickness of the stem and production) were evaluated each two months, except for production which was annual at the end of each harvest. Moreover, the economic analysis of the productivity of the crops was done. The statistical analysis of the data was done using SPSS v. 17.0. The means were proved by ANOVAS with a factor of a Tukey test (P < 0,05). The economic benefit was estimated in terms of incomes and total expenses which were presented in the essay. The results obtained at the end of the essay show that the MIN 2 and MIN 3 treatments produced more meaningful changes in comparison with the other treatments used to improve soil fertility; pH was less affected in the acidification compared with the organic treatments which were greatly acidified; organic matter (MO) had increased considerably in these two treatments; however, they were surpassed by the organic treatments of fertilization; total nitrogen (Nt) and potassium (K) also presented better results at the end of the essay and phosphorus (P) showed good increasing figures although a little lower compared with ORG 2 and ORG 3 treatments. Regarding the emission of the greenhouse gasses, the fluxes accumulated from nitrous oxide (N2O) in two years increased in all the treatments in comparison with the blank treatment, but in a greater form and with higher fluxes in the MIN 3 and MIN 2 treatments which can be considered as the ones with greater contamination of N2O in the atmosphere, this can be due to the applied mineral doses to fertilize during the process; MIN 1 treatments and all the organic ones showed lower emission to the atmosphere. Methane emissions (CH4) did not show major differences in emissions in the treatments, being the greater emissions the ORG 3 and ORG 2 treatments; this is possibly due to the organic compost added to the soil; regarding carbon dioxide (CO2) emissions, in a similar way to CH4, the ORG 3 treatment was the one that presented greater emissions, the CO2 emissions to the atmosphere in the other treatments were lower and did not present meaningful differences among them. The physiological variables in all the cases helped coffee crops grow, this was observed when compared with the blank treatment; however, plants that reached the greatest height, width of top and diameter of stem were the plants of the MIN 3 treatment, followed by MIN 3, which did not show much significance among them, and for the organic treatments, the one that presented great results in these variables was ORG 3, which did not show meaningful differences compared with MIN 2, which proves that mineral fertilization is more effective in this case compared with the organic. In the first year of production, the mineral treatment with MIN 3 fertilization obtained greater production and thus did not show statistical difference with MIN 2 treatment, although the other treatments were greater. It is worth mentioning that MIN 1 treatment and ORG 3 treatment presented a meaningful production of coffee, not showing statistical differences among them. For the second year, the production of the crops showed greater profits than in the first year of evaluation in all the treatments, this was due to the physiological properties of the crops and on the other hand, it might be due to the addition of fertilizers during the whole essay; in a similar way, MIN 3 and MIN 2 performed better, not showing greater statistical differences among them, although the mineral treatment MEDIUM doses did not show statistical difference compared with ORG 3. The economic benefit was greater in the MIN 3 and MIN 2 treatments, although MIN 2 treatment is the one that shows the best cost-benefit ratios; ORG 2 and ORG 3 treatments and the blank produced negative benefits for the producer. Regarding the environment, the best treatments to care for the atmosphere are considered to be MIN 1 and ORG 1 treatments; however, regarding production volume and profitability they were low for the producer.

Modelización de servicios ecosistémicos en una cuenca hidrológica del área mediterránea = Modeling ecosystem services in a mediterranean area river basin

La modelización es un proceso por el que se obtienen modelos de los procesos del ´mundo real´ mediante la utilización de simplificaciones. Sin embargo, las estimaciones obtenidas con el modelo llevan implícitas incertidumbre que se debe evaluar. Mediante un análisis de sensibilidad se puede mejorar la confianza en los resultados, sin embargo, este paso a veces no se realiza debido básicamente al trabajo que lleva consigo este tipo de análisis. Además, al crear un modelo, hay que mantener un equilibrio entre la obtención de resultados lo más exactos posible mediante un modelo lo más sencillo posible. Por ello, una vez creado un modelo, es imprescindible comprobar si es necesario o no incluir más procesos que en un principio no se habían incluido. Los servicios ecosistémicos son los procesos mediante los cuales los ecosistemas mantienen y satisfacen el bienestar humano. La importancia que los servicios ecosistémicos y sus beneficios asociados tienen, junto con la necesidad de realizar una buena gestión de los mismos, han estimulado la aparición de modelos y herramientas para cuantificarlos. InVEST (Integrated Valuation of Ecosystem Services and Tradoffs) es una de estas herramientas específicas para calcular servicios eco-sistémicos, desarrollada por Natural Capital Project (Universidad de Stanford, EEUU). Como resultado del creciente interés en calcular los servicios eco-sistémicos, se prevé un incremento en la aplicación del InVEST. La investigación desarrollada en esta Tesis pretende ayudar en esas otras importantes fases necesarias después de la creación de un modelo, abarcando los dos siguientes trabajos. El primero es la aplicación de un análisis de sensibilidad al modelo en una cuenca concreta mediante la metodología más adecuada. El segundo es relativo a los procesos dentro de la corriente fluvial que actualmente no se incluyen en el modelo mediante la creación y aplicación de una metodología que estudiara el papel que juegan estos procesos en el modelo InVEST de retención de nutrientes en el área de estudio. Los resultados de esta Tesis contribuirán a comprender la incertidumbre involucrada en el proceso de modelado. También pondrá de manifiesto la necesidad de comprobar el comportamiento de un modelo antes de utilizarlo y en el momento de interpretar los resultados obtenidos. El trabajo en esta Tesis contribuirá a mejorar la plataforma InVEST, que es una herramienta importante en el ámbito de los servicios de los ecosistemas. Dicho trabajo beneficiará a los futuros usuarios de la herramienta, ya sean investigadores (en investigaciones futuras), o técnicos (en futuros trabajos de toma de decisiones o gestión ecosistemas). ABSTRACT Modeling is the process to idealize real-world situations through simplifications in order to obtain a model. However, model estimations lead to uncertainties that have to be evaluated formally. The role of the sensitivity analysis (SA) is to assign model output uncertainty based on the inputs and can increase confidence in model, however, it is often omitted in modelling, usually as a result of the growing effort it involves. In addition, the balance between accuracy and simplicity is not easy to assess. For this reason, when a model is developed, it is necessary to test it in order to understand its behavior and to include, if necessary, more complexity to get a better response. Ecosystem services are the conditions and processes through which natural ecosystems, and their constituent species, sustain and fulfill human life. The relevance of ecosystem services and the need to better manage them and their associated benefits have stimulated the emergence of models and tools to measure them. InVEST, Integrated Valuation of Ecosystem Services and Tradoffs, is one of these ecosystem services-specific tools developed by the Natural Capital Project (Stanford University, USA). As a result of the growing interest in measuring ecosystem services, the use of InVEST is anticipated to grow exponentially in the coming years. However, apart from model development, making a model involves other crucial stages such as its evaluation and application in order to validate estimations. The work developed in this thesis tries to help in this relevant and imperative phase of the modeling process, and does so in two different ways. The first one is to conduct a sensitivity analysis of the model, which consists in choosing and applying a methodology in an area and analyzing the results obtained. The second is related to the in-stream processes that are not modeled in the current model, and consists in creating and applying a methodology for testing the streams role in the InVEST nutrient retention model in a case study, analyzing the results obtained. The results of this Thesis will contribute to the understanding of the uncertainties involved in the modeling process. It will also illustrate the need to check the behavior of every model developed before putting them in production and illustrate the importance of understanding their behavior in terms of correctly interpreting the results obtained in light of uncertainty. The work in this thesis will contribute to improve the InVEST platform, which is an important tool in the field of ecosystem services. Such work will benefit future users, whether they are researchers (in their future research), or technicians (in their future work in ecosystem conservation or management decisions).

Perfil de playas sustentadas : estudio hidrodinámico y experimental : directrices de diseño

Las playas sustentadas por medio de un pie sumergido son una atractiva alternativa de diseño de regeneración de playas especialmente cuando las condiciones física del emplazamiento o las características de la arena nativa y de préstamo producen perfiles de alimentación que no se intersectan. La observación y propuesta de este tipo de solución data de los años 1960’s, así como la experiencia internacional en la construcción de este tipo de playas. Sin embargo, a pesar de su utilización y los estudios en campo y laboratorio, no se dispone de criterios ingenieriles que apoyen el diseño de las mismas. Esta tesis consiste en un análisis experimental del perfil de playas sustentadas en un pie sumergido (o colgadas) que se concreta en una propuesta de directrices de diseño general que permiten estimar la ubicación y características geométricas del pie sumergido frente a un oleaje y material que constituye la playa determinados. En la tesis se describe el experimento bidimensional realizado en el modelo físico de fondo móvil, donde se combinan cinco tipos de oleaje con tres configuraciones del pie sumergido (“Sin estructura”, configuración baja o “Estructura 1” y configuración alta o “Estructura 2”), se presentan los resultados obtenidos y se realiza una discusión detallada de las implicaciones de los resultados desde el punto de vista hidrodinámico utilizando monomios adimensionales. Se ha realizado un análisis detallado del estado del arte sobre playas colgadas, presentando el concepto y las experiencias de realizaciones en distintos países. Se ha realizado una cuidadosa revisión de la literatura publicada sobre estudios experimentales de playas colgadas, modelos teóricos y otras cuestiones auxiliares, necesarias para la formulación de la metodología de la tesis. El estudio realizado se ha estructurado en dos fases. En la primera fase se ha realizado una experimentación en un modelo físico de fondo móvil construido en las instalaciones del Centro de Estudios de Puertos y Costas (CEPYC) del Centro de Estudios y Experimentación de Obras Públicas (CEDEX), consistente en un canal de 36 m de longitud, 3 m de anchura y 1.5 m de altura, provisto de un generador de oleaje de tipo pistón. Se ha diseñado una campaña de 15 ensayos, que se obtienen sometiendo a cinco tipos de oleaje tres configuraciones distintas de playa colgada. En los ensayos se ha medido el perfil de playa en distintos instantes hasta llegar al equilibrio, determinando a partir de esos datos el retroceso de la línea de costa y el volumen de sedimentos perdido. El tiempo total efectivo de ensayo asciende a casi 650 horas, y el número de perfiles de evolución de playa obtenidos totaliza 229. En la segunda fase se ha abordado el análisis de los resultados obtenidos con la finalidad de comprender el fenómeno, identificar las variables de las que depende y proponer unas directrices ingenieriles de diseño. Se ha estudiado el efecto de la altura de ola, del periodo del oleaje, del francobordo adimensional y del parámetro de Dean, constatándose la dificultad de comprensión del funcionamiento de estas obras ya que pueden ser beneficiosas, perjudiciales o inocuas según los casos. También se ha estudiado la respuesta del perfil de playa en función de otros monomios adimensionales, tales como el número de Reynolds o el de Froude. En el análisis se ha elegido el monomio “plunger” como el más significativo, encontrando relaciones de éste con el peralte de oleaje, la anchura de coronación adimensional, la altura del pie de playa adimensional y el parámetro de Dean. Finalmente, se propone un método de diseño de cuatro pasos que permite realizar un primer encaje del diseño funcional de la playa sustentada frente a un oleaje de características determinadas. Las contribuciones más significativas desde el punto de vista científico son: - La obtención del juego de resultados experimentales. - La caracterización del comportamiento de las playas sustentadas. - Las relaciones propuestas entre el monomio plunger y las distintas variables explicativas seleccionadas, que permiten predecir el comportamiento de la obra. - El método de diseño propuesto, en cuatro pasos, para este tipo de esquemas de defensa de costas. Perched beaches are an attractive beach nourishment design alternative especially when either the site conditions or the characteristics of both the native and the borrow sand lead to a non-intersecting profile The observation and suggestion of the use of this type of coastal defence scheme dates back to the 1960’s, as well as the international experience in the construction of this type of beaches. However, in spite of its use and the field and laboratory studies performed to-date, no design engineering guidance is available to support its design. This dissertation is based on the experimental work performed on a movable bed physical model and the use of dimensionless parameters in analyzing the results to provide general functional design guidance that allow the designer, at a particular stretch of coast - to estimate the location and geometric characteristics of the submerged sill as well as to estimate the suitable sand size to be used in the nourishment. This dissertation consists on an experimental analysis of perched beaches by means of a submerged sill, leading to the proposal of general design guidance that allows to estimate the location and geometric characteristics of the submerged sill when facing a wave condition and for a given beach material. The experimental work performed on a bi-dimensional movable bed physical model, where five types of wave conditions are combined with three configurations of the submerged sill (“No structure”, low structure or “Structure 1”, and high structure or “Structure 2”) is described, results are presented, and a detailed discussion of the results - from the hydrodynamic point of view – of the implications of the results by using dimensionless parameters is carried out. A detailed state of the art analysis about perched beaches has been performed, presenting the “perched beach concept” and the case studies of different countries. Besides, a careful revision of the literature about experimental studies on perched beaches, theoretical models, and other topics deemed necessary to formulate the methodology of this work has been completed. The study has been divided into two phases. Within the first phase, experiments on a movable-bed physical model have been developed. The physical model has been built in the Centro de Estudios de Puertos y Costas (CEPYC) facilities, Centro de Estudios y Experimentación de Obras Públicas (CEDEX). The wave flume is 36 m long, 3 m wide and 1.5 m high, and has a piston-type regular wave generator available. The test plan consisted of 15 tests resulting from five wave conditions attacking three different configurations of the perched beach. During the development of the tests, the beach profile has been surveyed at different intervals until equilibrium has been reached according to these measurements. Retreat of the shoreline and relative loss of sediment in volume have been obtained from the measurements. The total effective test time reaches nearly 650 hours, whereas the total number of beach evolution profiles measured amounts to 229. On the second phase, attention is focused on the analysis of results with the aim of understanding the phenomenon, identifying the governing variables and proposing engineering design guidelines. The effect of the wave height, the wave period, the dimensionless freeboard and of the Dean parameter have been analyzed. It has been pointed out the difficulty in understanding the way perched beaches work since they turned out to be beneficial, neutral or harmful according to wave conditions and structure configuration. Besides, the beach profile response as a function of other dimensionless parameters, such as Reynolds number or Froude number, has been studied. In this analysis, the “plunger” parameter has been selected as the most representative, and relationships between the plunger parameter and the wave steepness, the dimensionless crest width, the dimensionless crest height, and the Dean parameter have been identified. Finally, an engineering 4-step design method has been proposed, that allows for the preliminary functional design of the perched beach for a given wave condition. The most relevant contributions from the scientific point of view have been: - The acquisition of a consistent set of experimental results. - The characterization of the behavior of perched beaches. - The proposed relationships between the plunger parameter and the different explanatory variables selected, that allow for the prediction of the beach behavior. - The proposed design method, four-step method, for this type of coastal defense schemes.

Optimización de operaciones unitarias utilizando líquidos iónicos como agentes materiales de separación

En este trabajo se estudia la modelización y optimización de procesos industriales de separación mediante el empleo de mezclas de líquidos iónicos como disolventes. Los disolventes habitualmente empleados en procesos de absorción o extracción suelen ser componentes orgánicos muy volátiles y dañinos para la salud humana. Las innovadoras propiedades que presentan los líquidos iónicos, los convierten en alternativas adecuadas para solucionar estos problemas. La presión de vapor de estos compuestos es muy baja y apenas varía con la temperatura. Por tanto, estos compuestos apenas se evaporan incluso a temperaturas altas. Esto supone una gran ventaja en cuanto al empleo de estos compuestos como disolventes industriales ya que permite el reciclaje continuo del disolvente al final del proceso sin necesidad de introducir disolvente fresco debido a la evaporación del mismo. Además, al no evaporarse, estos compuestos no suponen un peligro para la salud humana por inhalación; al contrario que otros disolventes como el benceno. El único peligro para la salud que tienen estos compuestos es por tanto el de contacto directo o ingesta, aunque de hecho muchos Líquidos Iónicos son inocuos con lo cual no existe peligro para la salud ni siquiera a través de estas vías. Los procesos de separación estudiados en este trabajo, se rigen por la termodinámica de fases, concretamente el equilibrio líquido-vapor. Para la predicción de los equilibrios se ha optado por el empleo de modelos COSMO (COnductor-like Screening MOdel). Estos modelos tienen su origen en el empleo de la termodinámica de solvatación y en la mecánica cuántica. En el desarrollo de procesos y productos, químicos e ingenieros frecuentemente precisan de la realización de cálculos de predicción de equilibrios de fase. Previamente al desarrollo de los modelos COSMO, se usaban métodos de contribución de grupos como UNIFAC o modelos de coeficientes de actividad como NRTL.La desventaja de estos métodos, es que requieren parámetros de interacción binaria que únicamente pueden obtenerse mediante ajustes por regresión a partir de resultados experimentales. Debido a esto, estos métodos apenas tienen aplicabilidad para compuestos con grupos funcionales novedosos debido a que no se dispone de datos experimentales para llevar a cabo los ajustes por regresión correspondientes. Una alternativa a estos métodos, es el empleo de modelos de solvatación basados en la química cuántica para caracterizar las interacciones moleculares y tener en cuenta la no idealidad de la fase líquida. Los modelos COSMO, permiten la predicción de equilibrios sin la necesidad de ajustes por regresión a partir de resultados experimentales. Debido a la falta de resultados experimentales de equilibrios líquido-vapor de mezclas en las que se ven involucrados los líquidos iónicos, el empleo de modelos COSMO es una buena alternativa para la predicción de equilibrios de mezclas con este tipo de materiales. Los modelos COSMO emplean las distribuciones superficiales de carga polarizada (sigma profiles) de los compuestos involucrados en la mezcla estudiada para la predicción de los coeficientes de actividad de la misma, definiéndose el sigma profile de una molécula como la distribución de probabilidad de densidad de carga superficial de dicha molécula. Dos de estos modelos son COSMO-RS (Realistic Solvation) y COSMO-SAC (Segment Activity Coefficient). El modelo COSMO-RS fue la primera extensión de los modelos de solvatación basados en continuos dieléctricos a la termodinámica de fases líquidas mientras que el modelo COSMO-SAC es una variación de este modelo, tal y como se explicará posteriormente. Concretamente en este trabajo se ha empleado el modelo COSMO-SAC para el cálculo de los coeficientes de actividad de las mezclas estudiadas. Los sigma profiles de los líquidos iónicos se han obtenido mediante el empleo del software de química computacional Turbomole y el paquete químico-cuántico COSMOtherm. El software Turbomole permite optimizar la geometría de la molécula para hallar la configuración más estable mientras que el paquete COSMOtherm permite la obtención del perfil sigma del compuesto mediante el empleo de los datos proporcionados por Turbomole. Por otra parte, los sigma profiles del resto de componentes se han obtenido de la base de datos Virginia Tech-2005 Sigma Profile Database. Para la predicción del equilibrio a partir de los coeficientes de actividad se ha empleado la Ley de Raoult modificada. Se ha supuesto por tanto que la fracción de cada componente en el vapor es proporcional a la fracción del mismo componente en el líquido, dónde la constante de proporcionalidad es el coeficiente de actividad del componente en la mezcla multiplicado por la presión de vapor del componente y dividido por la presión del sistema. Las presiones de vapor de los componentes se han obtenido aplicando la Ley de Antoine. Esta ecuación describe la relación entre la temperatura y la presión de vapor y se deduce a partir de la ecuación de Clausius-Clapeyron. Todos estos datos se han empleado para la modelización de una separación flash usando el algoritmo de Rachford-Rice. El valor de este modelo reside en la deducción de una función que relaciona las constantes de equilibrio, composición total y fracción de vapor. Para llevar a cabo la implementación del modelado matemático descrito, se ha programado un código empleando el software MATLAB de análisis numérico. Para comprobar la fiabilidad del código programado, se compararon los resultados obtenidos en la predicción de equilibrios de mezclas mediante el código con los resultados obtenidos mediante el simulador ASPEN PLUS de procesos químicos. Debido a la falta de datos relativos a líquidos iónicos en la base de datos de ASPEN PLUS, se han introducido estos componentes como pseudocomponentes, de manera que se han introducido únicamente los datos necesarios de estos componentes para realizar las simulaciones. El modelo COSMO-SAC se encuentra implementado en ASPEN PLUS, de manera que introduciendo los sigma profiles, los volúmenes de la cavidad y las presiones de vapor de los líquidos iónicos, es posible predecir equilibrios líquido-vapor en los que se ven implicados este tipo de materiales. De esta manera pueden compararse los resultados obtenidos con ASPEN PLUS y como el código programado en MATLAB y comprobar la fiabilidad del mismo. El objetivo principal del presente Trabajo Fin de Máster es la optimización de mezclas multicomponente de líquidos iónicos para maximizar la eficiencia de procesos de separación y minimizar los costes de los mismos. La estructura de este problema es la de un problema de optimización no lineal con variables discretas y continuas, es decir, un problema de optimización MINLP (Mixed Integer Non-Linear Programming). Tal y como se verá posteriormente, el modelo matemático de este problema es no lineal. Por otra parte, las variables del mismo son tanto continuas como binarias. Las variables continuas se corresponden con las fracciones molares de los líquidos iónicos presentes en las mezclas y con el caudal de la mezcla de líquidos iónicos. Por otra parte, también se ha introducido un número de variables binarias igual al número de líquidos iónicos presentes en la mezcla. Cada una de estas variables multiplican a las fracciones molares de sus correspondientes líquidos iónicos, de manera que cuando dicha variable es igual a 1, el líquido se encuentra en la mezcla mientras que cuando dicha variable es igual a 0, el líquido iónico no se encuentra presente en dicha mezcla. El empleo de este tipo de variables obliga por tanto a emplear algoritmos para la resolución de problemas de optimización MINLP ya que si todas las variables fueran continuas, bastaría con el empleo de algoritmos para la resolución de problemas de optimización NLP (Non-Linear Programming). Se han probado por tanto diversos algoritmos presentes en el paquete OPTI Toolbox de MATLAB para comprobar cuál es el más adecuado para abordar este problema. Finalmente, una vez validado el código programado, se han optimizado diversas mezclas de líquidos iónicos para lograr la máxima recuperación de compuestos aromáticos en un proceso de absorción de mezclas orgánicas. También se ha usado este código para la minimización del coste correspondiente a la compra de los líquidos iónicos de la mezcla de disolventes empleada en la operación de absorción. En este caso ha sido necesaria la introducción de restricciones relativas a la recuperación de aromáticos en la fase líquida o a la pureza de la mezcla obtenida una vez separada la mezcla de líquidos iónicos. Se han modelizado los dos problemas descritos previamente (maximización de la recuperación de Benceno y minimización del coste de operación) empleando tanto únicamente variables continuas (correspondientes a las fracciones o cantidades molares de los líquidos iónicos) como variables continuas y binarias (correspondientes a cada uno de los líquidos iónicos implicados en las mezclas).

Determinación y caracterización de los indicadores que intervienen en un sistema integral de iluminación natural y protección solar : aplicado a la incorporación de bandejas en los huecos de fachadas existentes

En esta tesis se estudia cómo lograr el equilibrio entre dos fenómenos naturales que afectan a los huecos de fachadas: la iluminación natural y la ganancia solar. Es decir, el cómo, conseguir la optimización de la iluminación natural que se introduce a través de las ventanas existentes, sin realizar una laboriosa intervención de sustitución de las mismas y el cómo, conseguir la protección de la radiación solar directa de la zona acristalada, para evitar tanto las ganancias térmicas innecesarias como el deslumbramiento que afecta a la calidad lumínica de los recintos. Para el desarrollo esta investigación se ha propuesto una metodología de estudio dividida en dos fases: La primera, de Análisis y Diagnóstico, en la que se han de definir los estándares a cumplir y las variables con las que se evaluarán las bandejas. Y, la segunda, de Evaluación y Comprobación en la que se han de establecer los criterios de valoración y ponderación de cada variable. En la primera fase, se definirán las variables físico-ambientales, para lo que se seleccionarán algunas ciudades a estudiar, a las cuales se les estudiará las necesidades de confort térmico, se determinarán las dimensiones que deben tener las protecciones solares en esas ciudades, además se determinarán las actuaciones a realizar en los huecos de fachada según sea su posición en el plano vertical de la misma. Así mismo, se plantea hacer un análisis de casos reales, para lo que se caracterizarán desde el punto de vista lumínico algunas viviendas, realizando medidas “in situ” y comparando resultados con los que se obtienen de los programas de simulación, para seleccionar las viviendas en las que la realidad y la simulación se aproximen más, una de estas viviendas servirá de modelo en las simulaciones que se realizará en la segunda fase. También, en esta primera fase, mediante un modelo neutro, se estudiará el comportamiento térmico y lumínico del tamaño del hueco en el que se insertará la bandeja posteriormente, para luego estudiar la posición de la bandeja en el plano vertical de la ventana, desde el punto de vista ergonómico. Y finalmente se estudiará, el comportamiento térmico y lumínico del modelo con la bandeja ubicada a 40, 50 y 60cm del techo. En la segunda fase, se establecerá la valoración y ponderación de las variables con las que seleccionar la bandeja que mejor equilibre los aspectos térmicos y lumínicos, teniendo en cuenta estrategias pasivas de acondicionamiento ambiental, como favorecer las ganancias solares en invierno en horas diurnas y evitar las pérdidas de calor en horas nocturnas; y en verano implementar sistemas de sombreamiento en la zona acristalada para evitar las ganancias de calor; y, tanto en verano como en invierno, aprovechar la iluminación natural, para favorecer la iluminancia útil y evitar el deslumbramiento. Una vez definidos los criterios de valoración y ponderación se aplicará a la evaluación térmica y lumínica del modelo neutro con la bandeja, consiguiendo seleccionar la bandeja con mejor comportamiento. Posteriormente se comprobará la metodología de estudio desarrollada en el modelo seleccionado, se evaluará el comportamiento térmico y lumínico, con la incorporación de algunas alternativas de bandeja. Con esta investigación se quiere demostrar que mediante la aplicación de esta metodología de estudio, es posible evaluar y seleccionar bandejas que respondan a las necesidades requeridas en distintos casos de estudio, por lo que se considera que, la bandeja puede ser un elemento arquitectónico aplicable tanto en rehabilitación como en nueva construcción, de espacios en los que sea necesario mejorar sus condiciones lumínicas y térmicas simultáneamente. ABSTRACT This thesis studies how to balance two natural events that affect the window opening of facades: daylighting and solar gain. That is to say, how to achieve optimization of natural light that gets in through the existing windows, without making a laborious intervention of replacing them and how to get protection from direct solar radiation from the glass area, to avoid unnecessary heat gain and glare affecting the light quality of the enclosures. To develop this research, it has been proposed a methodology of study divided into two phases: First phase, Analysis and Diagnostics, in which the variables with which the light shelf are evaluated will be defined along with the standards the light shelves will meet. The second phase, Assessment and Verification, in which the assessment criteria and weighting of each variable will be established. In the first phase, the physical and environmental variables shall be defined, various cities will be selected to be studied, and in each the needs of thermal comfort will be determined along with the dimensions of shading devices in the cities. In addition the actions to be taken in the window opening of the façade will be determined, depending on their position in the vertical plane. An analysis of real cases will be undertaken, which will be characterized from the luminous point of view, performing "in situ" measurements and comparing results with those obtained from simulation programs, to select places/dwellings where reality and simulation are closer, one of these places/dwellings will be a model, in the simulations to perform at the second phase. Also, in this first phase, by a neutral model, the thermal and light behavior of the size of the window opening will be studied, in which the light shelf is inserted later, the position of the light shelf in the vertical plane of the window is studied, from an ergonomic point of view. And finally to study the thermal and light behavior of the model with the light shelf located at 40, 50 and 60cm from the ceiling. In the second phase, the evaluation and weighting of the variables will be established selecting the light shelf that best balances the thermal and daylighting aspects, taking into account passive environmental conditioning strategies; such as getting solar gains in winter during daylight hours, and preventing heat loss during the night hours; and in summer implementing shading systems in the glazing area to avoid heat gains. And in both summer and winter, taking advantage of natural lighting, to improve useful illuminance and avoid glare. Once defined, the evaluation criteria and weighting will be applied to thermal and daylighting evaluation to the neutral model with the light shelf, the best performing light shelf will be selected. The study methodology developed in the selected model will be verified the thermal and daylighting performance with the addition of some light shelf alternative will also be studied. With this research, we want to show that by applying this study methodology it is possible to evaluate and select the light shelf that meets the needs required in different case studies, so it is considered that the light shelf may be an applicable architectural element in both refurbishment and new construction of spaces where necessary to improve their daylighting and thermal conditions simultaneously.

Mecánica de las bóvedas de fábrica: el enfoque del equilibrio

El padre Vicente tosca, arquitecto, filósofo, matemático, y astrónomo, Y comienza su Tratado de Arquitectura (que forma parte de su Compendio matemático, 9 vols. 1701-1715) abordando el tema de la traza o proyecto de bóvedas: Lo más sutil y primoroso de la Arquitectura... es la construcción de todo género de arcos y bóvedas, cortando sus piedras, y ajustándolas con tal artificio, que la misma gravedad y peso que las había de precipitar hácia la tierra, las mantenga constantes en el ayre, sustentándose las unas á las otras en virtud de la mutua complicación que las enlaza, con lo que cierran por arriba las fábricas con toda seguridad y firmeza. El equilibrio se alcanza a través de la geometría y, de este modo es posible la construcción de edificios de fábrica seguros. Las antiguas reglas tradicionales para el proyecto de bóvedas y estribos de fábrica tienen un carácter geométrico, ya que establecen ciertas relaciones entre las dimensiones de los elementos estructurales. Regulan, por ejemplo, que el grosor de un estribo debe ser una cierta fracción entera de la luz de las bóvedas. De hecho en la afirmación del padre Tosca encontramos la esencia misma del proyecto de estructuras de fábrica (Huerta 2004). Sin embargo para nosotros, arquitectos e ingenieros del siglo XXI, todo esto nos parece demasiado ingenuo, la demostración de la ignorancia de los antiguos maestros; de hecho, hasta el siglo XVIII no se desarrolló una teoría científica de las estructuras, basada en la Resistencia de Materiales y las leyes de la Mecánica (Huerta 1996). De cualquier forma estos *ignorantes+ maestros constructores levantaron el Pantheon de Roma, Santa Sofía y las catedrales góticas. Por tanto, puede que, después de todo, el enfoque geométrico tradicional no esté demasiado equivocado. Es posible que los antiguos maestros tuvieran una teoría, distinta de nuestra teoría científica, pero basada en un profundo conocimiento de la naturaleza y comportamiento de las estructuras de fábrica. Si fuera así, sería interesante conocer algo sobre esta teoría, la cual, si juzgamos por los resultados, era extraordinaria. Sin embargo, aunque podemos ampliar nuestros conocimientos, no podemos sustraernos a los que ya tenemos. Estamos forzados a abordar el tema del análisis de arcos y bóvedas desde un enfoque científico, dentro del marco de la teoría de estructuras actual. Antes de proseguir se deben hacer dos observaciones. La primera se refiere al objetivo de la Teoría de Estructuras. El objetivo de la teoría estructural es el proyecto de construcciones seguras o el análisis de la seguridad de las que ya existen. Es una ciencia aplicada no una ciencia pura. Como Rankine (1858) señaló, si la pregunta que se hace el científico es *qué quiero saber+, la del arquitecto o ingeniero es *qué quiero hacer+. Todas las consideraciones teóricas están condicionadas por la necesidad de una respuesta a un problema sin demora. La segunda observación hace referencia a nuestra propia ignorancia respecto al tema. Ya no se construyen bóvedas y toda la tradición de la construcción en fábrica se ha perdido en el mundo occidental. Esta parte de la arquitectura, que fue considerada *lo más sutil y primoroso+, nos es ajena. La mayor parte de los arquitectos e ingenieros nunca han visto construir una sencilla bóveda. Carecemos del oficio y experiencia del antiguo constructor, que seleccionaba la piedra, dibujaba las plantillas para cortarla, trazaba la cimbra, dirigía el proceso de construcción y, finalmente, supervisaba el descimbrado. Rodrigo Gil de Hontañón (Tratado de arquitectura, ca. 1540, copiado en García, 1681) después de describir la construcción de una bóveda gótica de crucería advierte que: . . . estas cosas, podran ser difiçiles de comprehender faltando en quien las procura la experiencia, la practi¬ca, la profesion de la canteria, y la execuçion, o el aberse allado presente a algunos çierres de cruçeria, para haçerse capaz en el asiento de ella. Esta es precisamente la situación de cualquier arquitecto o ingeniero de la actualidad. A partir de aquí trataremos de mostrar que los enfoques tradicional y moderno del cálculo de estructuras de fábrica llevan a la misma conclusión fundamental: la absoluta importancia de la geometría. Se demostrará que la moderna teoría del Análisis Límite de las estructuras de Fábrica, que ha sido desarrollada por el profesor Heyman, es la herramienta más apropiada para comprender y analizar las construcciones de fábrica. Esta teoría conduce al *enfoque del equilibrio+; el calculista necesita sólo estudiar posible estados de equilibrio con la fábrica trabajando a compresión. La existencia de estos posibles estados de equilibrio depende de la geometría. Una construcción segura es una construcción en equilibrio: la teoría moderna conduce a las mismas disposiciones geométricas que la tradicional. No podía ser de otro modo, los espectaculares logros de los antiguos arquitectos no pudieron ser fruto de la casualidad. Se harán numerosas referencias históricas con la intención de dejar claro que existe una antigua tradición en el cálculo científico de las estructuras de fábrica según el enfoque del equilibrio. Tenemos mucho que aprender de los arquitectos e ingenieros de la antigüedad. Pudieron no tener una profunda comprensión de la teoría, pero poseían los conocimientos esenciales, ésos que proporciona la reflexión sobre la práctica. Ars sine scientia nihil est, la práctica no es nada sin la teoría, pero la teoría sin la práctica es simplemente peligrosa. La experiencia, en nuestro caso, debemos buscarla en los edificios construidos y en lo que podemos extraer de la lectura crítica de los antiguos tratados de arquitectura e ingeniería.