La geoquímica de los sedimentos de llanuras de inundación como guía para la valoración ambiental en España

En la presente tesis se estudian los contenidos geoquímicos de los sedimentos de llanura de inundación en diversas cuencas fluviales seleccionadas con el objetivo de contribuir a un mejor conocimiento de sus condiciones medioambientales. En cada cuenca, se ha muestreado un perfil vertical de llanura de inundación, dividiéndolo en tramos que generalmente corresponden a diferentes episodios de inundación. Estos sedimentos se depositan durante los episodios de crecidas una vez que la corriente sobrepasa los límites del canal. Se caracterizan normalmente por tener un tamaño de grano fino y una estructura en capas horizontales que corresponden a los sucesivos episodios de inundación. Muestran dos ventajas principales con respecto a otros medios de muestreo en geoquímica como suelos o sedimentos de corriente: • Pueden almacenar sedimento antiguo así como actual, con lo que se puede estudiar la historia geoquímica de una zona específica en un mismo punto de muestreo (perfil vertical). • Los sedimentos de llanuras de inundación son capaces de caracterizar grandes áreas de drenaje. El origen de los sedimentos es más diverso que en el sedimento de corriente, debido a las mayores áreas de donde proceden las aguas de las avenidas. Las cuencas han sido seleccionadas según las actividades antropogénicas que en ellas se llevan a cabo, en concreto, actividades urbanas e industriales, minería y agricultura. Así mismo, se han estudiado, como referencia, dos cuencas donde no se espera encontrar ningún tipo de actividad contaminante. Una vez hecha la selección, los sedimentos aluviales de las cuencas se han estudiado cuidadosamente para asegurar que no existen depósitos de acreción lateral en el punto seleccionado. Posteriormente se ha procedido al muestreo del perfil vertical. Las muestras han sido analizadas mediante ICP-MS (ataque total) e INAA para conocer los contenidos totales de los elementos traza y mayoritarios. Los análisis de la fracción extraíble se han llevado a cabo mediante ICP-MS (ataque con agua regia). Así mismo, algunas muestras seleccionadas han sido sometidas a una extracción secuencial para un estudio más detallado. La presencia de materia orgánica ha sido estimada mediante el análisis de Carbono Orgánico Total (TOC). Finalmente, se ha llevado a cabo un análisis de isótopos de Pb en muestras escogidas en los perfiles, con el objetivo de hacer una evaluación ambiental. Los contenidos metálicos aumentan hacia la superficie en algunos de los perfiles, mientras en otros muestran una distribución muy constante exceptuando algún nivel específico con un aumento de los contenidos de la mayoría de los metales. Ha sido posible determinar la influencia de las actividades antropogénicas en algunos de los perfiles. Aquellos que pertenecen a cuencas mineras, urbanas o industrializadas muestran generalmente altos contenidos en elementos metálicos. Es el caso de los perfiles muestreados en los ríos Odiel y Tinto, Besaya, Besós y Manzanares. Algunos de estos perfiles pueden incluso correlacionarse con periodos de tiempo en los que ha tenido lugar una actividad antropogénica más intensa. Los perfiles que mejor se correlacionan con la actividad antropogénica de la cuenca son el perfil de Rivas en el río Manzanares (Madrid), que refleja un crecimiento de la contaminación producida por las actividades urbana e industrial en las últimas décadas, y el río Tinto, que muestra un crecimiento llamativo de los contenidos en su mayoría metálicos que puede estar relacionado con el incremento de la actividad minera que tuvo lugar hace aproximadamente 125 años. El análisis de los isótopos de Pb ha resultado ser una herramienta útil en la evaluación ambiental de estos sedimentos. Con este estudio y mediante la comparación con fuentes naturales y antropogénicas, ha sido posible diferenciar las muestras afectadas por diferentes fuentes de plomo, así como detectar las más afectadas antropogénicamente. ABSTRACT The geochemical composition of overbank sediments of some selected river basins is studied in this thesis in order to contribute to a better knowledge of the environmental conditions surrounding them. In each basin a vertical overbank profile has been sampled, dividing it into stretches that usually correspond to different flood events. The overbank sediments are those deposited during a flood event once the flow spills over the channel banks. They are usually characterized by a very fine grain size and a structure of horizontal layers, which correspond to successive flood events. These sediments show two main advantages regarding other sampling media in geochemistry, like soils or stream sediments: • They can store sediment deposited in the past as well as in current times, so that the history of a specific location can be studied at the very same point (vertical profile). • The overbank sediments are able to characterize a large drainage area. The origin of the sediment is wider than in the stream sediments due to the larger areas where the flood water comes from. The basins have been selected depending on the anthropogenic activities developed in them, namely, urban and industrial activities, mining activities and agricultural activities. As well, two pristine basins have been studied as a reference. Afterwards, the alluvial sediments in the basins have been carefully studied in order to sample a vertical profile and make sure that lateral accretion materials are not present in the profile. The samples have been analysed by ICP-MS (total digestion) and INAA to know the total contents of trace and major elements. Analysis of the mobile fraction has been carried out by ICP-MS (aqua regia); as well some of the samples have been subjected to sequential extraction for a more detailed study. The presence of organic matter has been estimated by the analysis of the Total Organic Carbon (TOC). Finally, a lead isotope analysis of some of the samples in the profiles was carried out in order to make an environmental assessment. Metal contents grow towards the surface in some of the profiles, while others show a very steady distribution, except for some of them with a growth of most of the metals in a specific level. XI It has been possible to determine the influence of the anthropogenic activities in some of the profiles. The ones that belong to mining and urban or industrialized basins show generally high contents of metal elements. This is the case of the profiles sampled in the Odiel and Tinto Rivers, the Besaya River, the Besós River and the Manzanares River. Some of these profiles can even correlate with the periods of time when a more intense activity in their respective basins has taken place. The profiles which best correlate with the anthropogenic activity are the Rivas profile in the Manzanares River, which reflects a growth of the pollution produced by urban and industrial activities in the city of Madrid in the last decades and the Tinto profile, which shows a very dramatic growth of the elemental contents (mostly metals) which can be related to the increase of the mining activities that took place in the last 125 years. The analysis of lead isotopes has turned out to be a powerful tool in the environmental assessment in this type of sediments. With this study and through the comparison with natural and anthropogenic sources, it has been possible to determine samples affected by different sources of lead and to detect the most anthropogenicaly affected ones.

Origen y distribución de los hidrocarburos policíclicos aromáticos en sedimentos actuales de la Laguna de El Hito (España central)

Se llevó a cabo la evaluación ambiental y el estudio del estado actual de la cuenca de la Laguna de El Hito referido a 18 hidrocarburos policíclicos aromáticos (PAHs) de 2 a 6 anillos bencénicos. Se determinó su origen a partir de diversos índices (%naftaleno, Fen/Ant y Flu/Pir), interpretándose tanto fuentes no antropogénicas (petrogénica) como antropogénicas (pirogénica). Se obtuvieron los mapas de distribución de las concentraciones de PAHs y de sus índices para localizar los puntos de concentraciones más elevadas. Ningún PAH superó las concentraciones marcadas por los Niveles Genéricos de Referencia (NGR) para la salud humana en los distintos usos del suelo del R.D.09/2005. Los PAHs con las mayores concentraciones fueron el naftaleno y el fenantreno.Environmental evaluation and analysis of the current state of El Hito Lake Basin referred to 18 polycyclic aromatic hydrocarbons (PAHs) with 2 to 6 benzene rings was carried out. Different indexes were used to determine the source of PAHs (% naftalene, Phe/Ant and Flu/Pyr). Both non anthropogenic (petrogenic) and anthropogenic (pyrogenic) sources were interpreted. Distribution maps for PAHs and indexes were plotted to locate the position of the higher concentrations and, therefore, their possible sources. None of these compounds showed concentrations above the Soil Screening Levels (SSL) for human health in the different uses of soil as is established in R.D.09/2005. The ones that reached the highest concentrations were naphthalene and phenanthrene.

Origen y distribución de los pesticidas organoclorados (OCPs) en sedimentos actuales de la Laguna de El Hito (Cuenca, España Central)

Se llevó a cabo el estudio del estado ambiental actual de la Laguna de El Hito basado en el análisis y cuantificación de 24 pesticidas organoclorados (OCPs): hexaclorobenceno, hexaclorociclohexanos, diclorodifenil-tricloroetano sus homólogos y metabolitos (DDTs), y ciclodienos (aldrín, dieldrín, endrín, endrín aldehido, endrín cetona, ?-clordano, ?-clordano, endosulfán I, endosulfán II, endosulfán sulfato, heptacloro, heptacloro epóxido B y metoxicloro). Algunos compuestos superaron los Niveles Genéricos de Referencia (NGR) para la salud humana y los ecosistemas del R.D.09/2005. Se determinó su origen a partir de diversos índices interpretándose tanto aplicaciones históricas como más recientes. Se obtuvieron los mapas de distribución de los índices y de las concentraciones de OCPs.Environmental evaluation of the current state of El Hito Lake referred to 24 organochlorine pesticides (OCPs) was carried out: hexachlorobenzene, hexachlorocyclohexanes, dichlorodiphenyltrichloroethane their homologs and metabolites (DDTs), and cyclodienes (aldrin, dieldrin, endrin, endrin aldehyde, endrin ketone, α- chlordane, γ- chlordane, endosulfan I, endosulfan II, endosulfan sulphate, heptachlor, heptachlor epóxide B and metoxichlor). Some of these compounds showed concentrations above the Soil Screening Levels (SSLs) for human health and ecosystems as is established in R.D.09/2005. Different indexes were used to determine their origin. Both historical and recent applications sources were interpreted. Distribution maps for concentrations and indexes were plotted

Origen y distribución de tierras raras (RREs) en el agua y sedimentos de la Laguna de El Hito (Cuenca, España Central).

Determination of Rare Earth Elements (REEs) in the waters and sediments of El Hito Lake and its drainage basin was carried out. Distribution maps for REEs concentrations show increasing values in the lake sediments from the center to the eastern edge reaching 166.5 mg/kg. In the drainage basin, higher values were observed with a maximum of 350.9 mg/kg in the southern part. Concentrations in the water were considerably lower with a maximum value of 1.3 µg/L. Ce, La and Nd were the most abundant elements. When normalized REE concentrations against NASC (North American Shale Composite), a positive anomaly of Eu in the water and of Nd in the sediments (lake and basin) was observed. The (La/Gd)NASC and (La/Yb)NASC ratios determined the predominance of light REE (LREE) over medium (MREE) and heavy REE (HREE). Key-words: Rare Earth Elements, lacustrine sediments and water, Eu anomaly, sulfates RESUMEN: Se llevó a cabo el estudio de Tierras Raras (REEs) del agua y los sedimentos de la Laguna de El Hito así como de su cuenca de recepción. Los mapas de distribución de las concentraciones de REE en los sedimentos de la laguna mostraron valores crecientes desde el centro hacia el margen este alcanzando 166.5 mg/kg. En la cuenca se midieron valores más elevados llegando hasta 350.9 mg/kg en la mitad sur. Los valores en el agua fueron notablemente inferiores, con un máximo de 1.3 µg/L. El elemento más abundante fue el Ce seguido del Nd. Se normalizaron las concentraciones de REEs frente a las del NASC (North American Shale Composite), observándose una anomalía positiva de Eu en el agua, así como de Nd en los sedimentos (laguna y cuenca). Mediante los índices (La/Gd)NASC y (La/Yb)NASC se determinó el predominio de las REEs ligeras (LREEs) frente a las medias (MREE) y pesadas (HREE).

Empleo de biomarcadores para la reconstrucción paleoclimática. El registro pleistoceno del Mar de Fuentillejo (Ciudad Real)

El registro sedimentario del lago volcánico conocido como maar de Fuentillejo (Centro de España) ofrece la oportunidad de determinar la evolución paleoambiental durante el Pleistoceno medio y superior de la parte Centro-Sur de la Península Ibérica, objetivo principal de esta Tesis Doctoral. Se ha recopilado información sobre otras reconstrucciones paleoambientales en maares de todo el mundo y estudios paleoclimáticos realizados en el registro de la turbera de Padul, relativamente próxima al maar de Fuentillejo. Con ello, se pudo realizar el estado del arte.A lo largo del sondeo continuo (142,25 m) perforado en la parte central del maar de Fuentillejo en el año 2002, se han tomando un total de 613 muestras a intervalos de 20 cm para el análisis de biomarcadores. Los biomarcadores se obtuvieron mediante un equipo de extracción acelerada de disolventes (Dionex ASE 200) y fueron analizados con un cromatógrafo de gases (HP 6890) con detector selectivo de masas (HP 5973). La datación de sedimentos del maar de Fuentillejo ha sido complicada debido a la ausencia de restos orgánicos terrestres. Se seleccionaron seis muestras de sedimento para intentar datar por el método de racemización de aminoácidos: los resultados obtenidos no fueron buenos. Finalmente, el modelo de edad-profundidad se ha llevado a cabo mediante paleomagnetismo, dataciones 14C y U/Th. En base a una interpolación lineal de las edades medias calibradas se ha datado la base del registro en aproximadamente 355 ka. Se han identificado y cuantificado los n-alcanos, n-metil cetonas, ácidos n- alcanoicos y azufre orgánico del sedimento a lo largo del registro. La interpretación de los diferentes índices establecidos a partir de éstos permite independizar tres Unidades Geoquímicas generales (UG) y nueve Subunidades (S) que definen condiciones paleoambientales distintas (predominio en el aporte de vegetación terrestre/de macrofitas/acuática, lámina de agua más alta vs. más baja, periodos más húmedos vs. más secos).A partir de análisis estadísticos de los índices de n-alcanos (matriz de correlación y análisis de agrupaciones) se han definido tres asociaciones: Grupo 1 (condiciones ambientales más húmedas y fase con lámina de agua elevada), Grupo 2 (con condiciones intermedias) y Grupo 3 (las fases más secas y con lámina de agua más baja). Además, se establece un Indicador Paleoambiental (I.P.) que muestra variaciones notables correlacionables con escenarios paleoambientales distintos (más húmedos vs. más secos) distinguiendo veintiséis Episodios (E): cinco Secos, nueve Húmedos y doce Intermedios.En general, se encuentra una buena correspondencia entre estos Episodios (E) y los eventos paleoclimáticos cortos descritos por Vegas et al. (2010) en este registro, en base a la palinología y geoquímica, así como los Episodios del Oxígeno Marino (MIS) y Eventos Heinrich, que permiten suponer que las oscilaciones paleoclimatológicas registradas en el maar tienen un origen global. Esta correspondencia se confirma por el análisis espectral, que revela que los ciclos interpretados (103, 41, 26 y 23 ka) evidencian influencia orbital. Además, se ha realizado una correlación tentativa entre los episodios paleoambientales de larga duración obtenidos a partir del perfil suavizado del porcentaje relativo de C27 de la Cuenca de Padul (Granada, Sur de la Península Ibérica) y el perfil I.P. suavizado del maar de Fuentillejo, aunque hay diferencias que podrían ser explicadas por la vecindad de Padul a la zona de nieves perpetuas de Sierra Nevada lo que origina marcadas diferencias hidrológicas entre ambos sitios. El análisis por cadenas de Markov ha permitido determinar que los cambios paleoambientales reflejados por los biomarcadores no fueron, en la mayoría de los casos, bruscos pasando, normalmente, por estados transicionales.

Optimización de la succión de sólidos sumergidos utilizando bombas centrífugas con velocidad variable

La acumulación de material sólido en embalses, cauces fluviales y en zonas marítimas hace que la extracción mecánica de estos materiales por medio de succión sea cada vez mas frecuente, por ello resulta importante estudiar el rendimiento de la succión de estos materiales analizando la forma de las boquillas y los parámetros del flujo incluyendo la bomba. Esta tesis estudia, mediante equipos experimentales, la eficacia de distintos dispositivos de extracción de sólidos (utilizando boquillas de diversas formas y bombas de velocidad variable). El dispositivo experimental ha sido desarrollado en el Laboratorio de Hidráulica de la E.T.S.I. de Caminos, C. y P. de la Universidad Politécnica de Madrid. Dicho dispositivo experimental incluye un lecho sumergido de distintos tipos de sedimentos, boquillas de extracción de sólidos y bomba de velocidad variable, así como un elemento de separación del agua y los sólidos extraídos. Los parámetros básicos analizados son el caudal líquido total bombeado, el caudal sólido extraído, diámetro de la tubería de succión, forma y sección de la boquilla extractora, así como los parámetros de velocidad y rendimiento en la bomba de velocidad variable. Los resultados de las medidas obtenidas en el dispositivo experimental han sido estudiados por medio del análisis dimensional y con métodos estadísticos. A partir de este estudio se ha desarrollado una nueva formulación, que relaciona el caudal sólido extraído con los diámetros de tubería y boquilla, caudal líquido bombeado y velocidad de giro de la bomba. Así mismo, desde el punto de vista práctico, se han analizado la influencia de la forma de la boquilla con la capacidad de extracción de sólidos a igualdad del resto de los parámetros, de forma que se puedan recomendar que forma de la boquilla es la más apropiada. The accumulation of solid material in reservoirs, river channels and sea areas causes the mechanical extraction of these materials by suction is becoming much more common, so it is important to study the performance of the suction of these materials analyzing the shape of the nozzles and flow parameters, including the pump. This thesis studies, using experimental equipment, the effectiveness of different solids removal devices (using nozzles of different shapes and variable speed pumps). The experimental device was developed at the Hydraulics Laboratory of the Civil University of the Polytechnic University of Madrid. The device included a submerged bed with different types of sediment solids, different removal nozzles and variable speed pump. It also includes a water separation element and the solids extracted. The key parameters analyzed are the total liquid volume pumped, the solid volume extracted, diameter of the suction pipe, a section of the nozzle and hood, and the parameters of speed and efficiency of the variable speed pump. The basic parameters analyzed are the total liquid volume pumped, the removed solid volume, the diameter of the suction pipe, the shape and cross-section of the nozzle, and the parameters of speed, efficiency and energy consumed by the variable speed pump. The measurements obtained on the experimental device have been studied with dimensional analysis and statistical methods. The outcome of this study is a new formulation, which relates the solid volume extracted with the pipe and nozzle diameters, the pumped liquid flow and the speed of the pump. Also, from a practical point of view, the influence of the shape of the nozzle was compared with the solid extraction capacity, keeping equal the rest of the parameters. So, a recommendation of the best shape of the nozzle can be given.

Afectación de los agrosistemas de la zona del Caño San Miguel por la salinización de los suelos

La degradación por salinización de los suelos regados con aguas salobres viene aumentando a escala mundial. El problema de la concentración de sales más solubles que el yeso depende principalmente del agua de riego, la aridez climática y la ausencia de drenaje. Estas condiciones se dan en el aluvium del río Limón, que es un tributario del lago Maracaibo, sito en el estado de Zulia de Venezuela. La regulación del río Limón mediante el cierre de los embalses de Manuelote y Tulé ha disminuido los aportes de aguas y sedimentos de las avenidas de inundación, que tienen carácter diluyente. Por otro lado, el balance de sales solubles en el suelo ha registrado una acumulación neta en los años de extrema aridez anteriores al año 2006, dado que la mayor dilución de las aguas ombrogénicas embalsadas procedentes de las lluvias no ha sido suficiente para compensar la concentración por evapotranspiración “in situ” de las aguas retenidas en la cuenca baja, sobre todo en ausencia de desagüe superficial y drenaje profundo. Las inundaciones posteriores a 2006 fueron suficientes para disminuir la salinidad superficial hasta los valores encontrados en 2010. El estudio experimental de esta problemática en el sector del caño San Miguel ha sido abordado mediante el establecimiento del perfil de salinidad acoplado con el perfil hipotético de humedad usado en la taxonomía de suelos. Este perfil define la disponibilidad del agua del suelo para la vegetación en función de tres potenciales: 1) el potencial físico-químico o matricial, que depende de la energía de adsorción a la superficie de las partículas; 2) el potencial gravitatorio, que depende de la profundidad; y 3) el potencial osmótico, que depende de la concentración de la solución del suelo; lo que supone un avance respecto a tener en cuenta sólo el perfil de humedad, que solamente considera el potencial gravi-químico integrado por el matricial y el gravitatorio. El perfil normalizado de 200 mm de de agua útil, retenida entre 33 y 1500 kPa de succión, incluye ocho fases gravi-químicas de 25 mm. La presente investigación incluye el potencial osmótico estimado por la conductividad eléctrica del extracto de pasta saturada. Los experimentos de lavado de sales en columnas de suelo, simulando la distribución de las lluvias en cinco años representativos de los cuartiles estadísticos de la serie disponible de 38 años completos, han determinado el comportamiento de las sales solubles en un suelo sometido a drenaje. Los resultados han evidenciado que el balance de sales unido al balance de agua controla la degradación de los agrosistemas por salinización. La alternativa frutícola puede ser aumentada en estas condiciones, porque el balance de sales favorece el establecimiento de cultivos permanentes a costa de otros usos del suelo de menor interés económico, como el cultivo de forrajes en regadío y el aprovechamiento de los pastizales en secano durante el barbecho de desalinización, cuya caracterización se ha completado con el estudio de la vegetación indicadora del grado de salinidad. ABSTRACT Saline degradation of soils irrigated with brackish water is increasing worldwide. The problem of salts concentration more soluble than gypsum depends on irrigation water quality, climatic aridity, and drainage limitations. These conditions meet in Limón River alluvium, which is tributary to Maracaibo´s Lake in Zulia State, Venezuela. Limón River regulation by closing Manuelote and Tulé reservoirs has diminished the input of water and sediments from inundations, which exerted dilutive effects. On the other hand, the soil balance of soluble salts has registered a net accumulation during those extremely dry years before 2006 because the greater dilution of ombrogenic dammed water coming from rain has not been enough to compensate salt concentration by “in situ” evapotranspiration in middle basin water, mainly in the absence of superficial runoff and deep drainage. Floods after 2006 were enough to reduce the high superficial salinity figures to those addressed in 2010. The experimental study of this trouble in San Miguel´s pipe area has been addressed through of the establishment of its salinity profile together to the hypothetic moisture profile typically used in soil taxonomy. This salinity profile describes soil water availability for vegetation according to three potentials: 1) physico-chemical or matrix potential, which depends on the adsorption energy of the soil solution to the surface of soil particles; 2) gravitational potential, which depends on soil depth; and 3) osmotic potential, which depends on the concentration of the soil solution. This represents an advance from just using moisture regime, which only considers the matrix and gravitational components of a gravi-chemical potential. The standardized moisture profile of 200 mm useful water being retained between 33 and 1500 kPa includes eight gravi-chemical stages of 25 mm. This research also includes the osmotic component, which is estimated by the electric conductivity of the saturated paste extract. Salts leaching trials in soil columns simulating rain distribution along five model years, representing the statistical quartiles of the available series of 38 complete years, have determined the behaviour of soluble salts in a soil being subjected to drainage. Results have evidenced that salt and water balances considered together are able to control the agrosystem’s degradation by salinization. The fruit production alternative could be improved under these conditions because the salts balance favours the establishment of permanent crops to the detriment of other soil uses of lower economical interest such as irrigated forage and non-irrigated pasture during desalinization fallow, which characterization has been completed through assessing the presence of salinity-indicator vegetation.

Interactions between roughness and topography in hydraulic models

Analysis of river flow using hydraulic modelling and its implications in derived environ-mental applications are inextricably connected with the way in which the river boundary shape is represented. This relationship is scale-dependent upon the modelling resolution which in turn determines the importance of a subscale performance of the model and the way subscale (surface and flow) processes are parameterised. Commonly, the subscale behaviour of the model relies upon a roughness parameterisation whose meaning depends on the dimensionality of the hydraulic model and the resolution of the topographic represen¬tation scale. This latter is, in turn, dependent on the resolution of the computational mesh as well as on the detail of measured topographic data. Flow results are affected by this interactions between scale and subscale parameterisation according to the dimensionality approach. The aim of this dissertation is the evaluation of these interactions upon hy¬draulic modelling results. Current high resolution topographic source availability induce this research which is tackled using a suitable roughness approach according to each di¬mensionality with the purpose of the interaction assessment. A 1D HEC-RAS model, a 2D raster-based diffusion-wave model with a scale-dependent distributed roughness parame-terisation and a 3D finite volume scheme with a porosity algorithm approach to incorporate complex topography have been used. Different topographic sources are assessed using a 1D scheme. LiDAR data are used to isolate the mesh resolution from the topographic content of the DEM effects upon 2D and 3D flow results. A distributed roughness parameterisation, using a roughness height approach dependent upon both mesh resolution and topographic content is developed and evaluated for the 2D scheme. Grain-size data and fractal methods are used for the reconstruction of topography with microscale information, required for some applications but not easily available. Sensitivity of hydraulic parameters to this topographic parameterisation is evaluated in a 3D scheme at different mesh resolu¬tions. Finally, the structural variability of simulated flow is analysed and related to scale interactions. Model simulations demonstrate (i) the importance of the topographic source in a 1D models; (ii) the mesh resolution approach is dominant in 2D and 3D simulations whereas in a 1D model the topographic source and even the roughness parameterisation impacts are more critical; (iii) the increment of the sensitivity to roughness parameterisa-tion in 1D and 2D schemes with detailed topographic sources and finer mesh resolutions; and (iv) the topographic content and microtopography impact throughout the vertical profile of computed 3D velocity in a depth-dependent way, whereas 2D results are not affected by topographic content variations. Finally, the spatial analysis shows that the mesh resolution controls high resolution model scale results, roughness parameterisation control 2D simulation results for a constant mesh resolution; and topographic content and micro-topography variations impacts upon the organisation of flow results depth-dependently in a 3D scheme. Resumen La topografía juega un papel fundamental en la distribución del agua y la energía en los paisajes naturales (Beven and Kirkby 1979; Wood et al. 1997). La simulación hidráulica combinada con métodos de medición del terreno por teledetección constituyen una poderosa herramienta de investigación en la comprensión del comportamiento de los flujos de agua debido a la variabilidad de la superficie sobre la que fluye. La representación e incorporación de la topografía en el esquema hidráulico tiene una importancia crucial en los resultados y determinan el desarrollo de sus aplicaciones al campo medioambiental. Cualquier simulación es una simplificación de un proceso del mundo real, y por tanto el grado de simplificación determinará el significado de los resultados simulados. Este razonamiento es particularmente difícil de trasladar a la simulación hidráulica donde aspectos de la escala tan diferentes como la escala de los procesos de flujo y de representación del contorno son considerados conjuntamente incluso en fases de parametrización (e.g. parametrización de la rugosidad). Por una parte, esto es debido a que las decisiones de escala vienen condicionadas entre ellas (e.g. la dimensionalidad del modelo condiciona la escala de representación del contorno) y por tanto interaccionan en sus resultados estrechamente. Y por otra parte, debido a los altos requerimientos numéricos y computacionales de una representación explícita de alta resolución de los procesos de flujo y discretización de la malla. Además, previo a la modelización hidráulica, la superficie del terreno sobre la que el agua fluye debe ser modelizada y por tanto presenta su propia escala de representación, que a su vez dependerá de la escala de los datos topográficos medidos con que se elabora el modelo. En última instancia, esta topografía es la que determina el comportamiento espacial del flujo. Por tanto, la escala de la topografía en sus fases de medición y modelización (resolución de los datos y representación topográfica) previas a su incorporación en el modelo hidráulico producirá a su vez un impacto que se acumulará al impacto global resultante debido a la escala computacional del modelo hidráulico y su dimensión. La comprensión de las interacciones entre las complejas geometrías del contorno y la estructura del flujo utilizando la modelización hidráulica depende de las escalas consideradas en la simplificación de los procesos hidráulicos y del terreno (dimensión del modelo, tamaño de escala computacional y escala de los datos topográficos). La naturaleza de la aplicación del modelo hidráulico (e.g. habitat físico, análisis de riesgo de inundaciones, transporte de sedimentos) determina en primer lugar la escala del estudio y por tanto el detalle de los procesos a simular en el modelo (i.e. la dimensionalidad) y, en consecuencia, la escala computacional a la que se realizarán los cálculos (i.e. resolución computacional). Esta última a su vez determina, el detalle geográfico con que deberá representarse el contorno acorde con la resolución de la malla computacional. La parametrización persigue incorporar en el modelo hidráulico la cuantificación de los procesos y condiciones físicas del sistema natural y por tanto debe incluir no solo aquellos procesos que tienen lugar a la escala de modelización, sino también aquellos que tienen lugar a un nivel subescalar y que deben ser definidos mediante relaciones de escalado con las variables modeladas explícitamente. Dicha parametrización se implementa en la práctica mediante la provisión de datos al modelo, por tanto la escala de los datos geográficos utilizados para parametrizar el modelo no sólo influirá en los resultados, sino también determinará la importancia del comportamiento subescalar del modelo y el modo en que estos procesos deban ser parametrizados (e.g. la variabilidad natural del terreno dentro de la celda de discretización o el flujo en las direcciones laterales y verticales en un modelo unidimensional). En esta tesis, se han utilizado el modelo unidimensional HEC-RAS, (HEC 1998b), un modelo ráster bidimensional de propagación de onda, (Yu 2005) y un esquema tridimensional de volúmenes finitos con un algoritmo de porosidad para incorporar la topografía, (Lane et al. 2004; Hardy et al. 2005). La geometría del contorno viene definida por la escala de representación topográfica (resolución de malla y contenido topográfico), la cual a su vez depende de la escala de la fuente cartográfica. Todos estos factores de escala interaccionan en la respuesta del modelo hidráulico a la topografía. En los últimos años, métodos como el análisis fractal y las técnicas geoestadísticas utilizadas para representar y analizar elementos geográficos (e.g. en la caracterización de superficies (Herzfeld and Overbeck 1999; Butler et al. 2001)), están promoviendo nuevos enfoques en la cuantificación de los efectos de escala (Lam et al. 2004; Atkinson and Tate 2000; Lam et al. 2006) por medio del análisis de la estructura espacial de la variable (e.g. Bishop et al. 2006; Ju et al. 2005; Myint et al. 2004; Weng 2002; Bian and Xie 2004; Southworth et al. 2006; Pozd-nyakova et al. 2005; Kyriakidis and Goodchild 2006). Estos métodos cuantifican tanto el rango de valores de la variable presentes a diferentes escalas como la homogeneidad o heterogeneidad de la variable espacialmente distribuida (Lam et al. 2004). En esta tesis, estas técnicas se han utilizado para analizar el impacto de la topografía sobre la estructura de los resultados hidráulicos simulados. Los datos de teledetección de alta resolución y técnicas GIS también están siendo utilizados para la mejor compresión de los efectos de escala en modelos medioambientales (Marceau 1999; Skidmore 2002; Goodchild 2003) y se utilizan en esta tesis. Esta tesis como corpus de investigación aborda las interacciones de esas escalas en la modelización hidráulica desde un punto de vista global e interrelacionado. Sin embargo, la estructura y el foco principal de los experimentos están relacionados con las nociones espaciales de la escala de representación en relación con una visión global de las interacciones entre escalas. En teoría, la representación topográfica debe caracterizar la superficie sobre la que corre el agua a una adecuada (conforme a la finalidad y dimensión del modelo) escala de discretización, de modo que refleje los procesos de interés. La parametrización de la rugosidad debe de reflejar los efectos de la variabilidad de la superficie a escalas de más detalle que aquellas representadas explícitamente en la malla topográfica (i.e. escala de discretización). Claramente, ambos conceptos están físicamente relacionados por un

Identificación de episodios climáticos fríos mediante el registro de susceptibilidad magnética en la secuencia lacustre del maar de Fuentillejo (Ciudad Real).

Resumen: El maar de Fuentillejo está localizado en la Región Volcánica Central de Campo de Calatrava (Ciudad Real). Desde su inicio se ha comportado como un sistema cerrado y presenta una potencia total de 142 m de sedimentos lacustres, organizados en 23 unidades sedimentarias. Se ha realizado un estudio de la susceptibilidad magnética y densidad aparente de los sedimentos, mediante un equipo de testificación multisensor Geotek. Para obtener la edad del registro se han efectuado dataciones absolutas mediante radiocarbono, U-Th y con el estudio de polaridad magnética, obteniéndose un modelo de edad que permite datar el sondeo FUENT-1 en torno a los 350 ka. Los datos del registro de susceptibilidad magnética, junto con las dataciones efectuadas, permiten identificar una secuencia de eventos erosivos relacionados con los últimos estadios glaciares (LGM: Ultimo Máximo Glaciar y los estadios isotópicos MIS 6 y 8). El registro de susceptibilidad magnética se encuentra atenuado por procesos de disolución de los óxidos de Fe-Ti. Abstract: The Fuentillejo maar is located in the Central Spanish Volcanic Field of Campo de Calatrava (Ciudad Real). Fuentillejo maar-lake was a closed system where up to 142 m depth of lacustrine sediments were deposited. Magnetic susceptibility and bulk density were measured by a GEOTEK multisensor core logger. The chronological framework was constructed based on radiocarbon and U-Th methods, as well as a detailed study of magnetic polarity, yielding an age model that covers last 350 ka. Intervals with terrigenous sediments correspond to high magnetic susceptibility values. These intervals were correlated with erosive events during the Last Glacial Maximum and the MIS 6 and 8 isotopic stages. The record of magnetic susceptibility is attenuated by disolution processes of Fe-Ti oxides.

Reparación de los desagües de fondo de la presa de la Fuensanta utilizando las válvulas de seguridad como escudos frontales

La presa de la Fuensanta se construyó hace más de 75 años. Los desagües de fondo, alojados en el tapón de cierre del túnel de desvío del río, son dos conductos rectangulares de sección abocinada, maniobrados con válvulas de compuerta que evacuan directamente el vertido a la galería. Dada la antigüedad de los equipos mecánicos, y su estado de deterioro, se ha procedido a su renovación. Los problemas habituales de estas actuaciones (existencia de sedimentos y desprendimientos en el área de la embocadura, falta de conocimiento de las formas geométricas de las zonas ocultas por el agua, defectos de construcción no documentados, deterioro de los blindajes, dificultad de acceso, etc) se han resuelto en esta obra acondicionando las antiguas válvulas de seguridad de los conductos para utilizarlas como escudos frontales, lo que ha permitido sustituir y ampliar los blindajes de acero, e instalar las nuevas compuertas de aguas abajo.

Reordenación ambiental y recuperación del entorno urbanístico de la Marina Baja en Benidorm, Alicante

El objeto de estudio de este proyecto es la Cala Finestrat, situada en la Marina Baja, en la frontera sur del municipio de Benidorm. La zona se encuentra bastante deteriorada, especialmente por el volumen turístico y de desarrollo de las inmediaciones de Benidorm. Esto ha originado en la cala una disminución de la calidad de las aguas, que bajan principalmente de la zona norte. Por otra parte, al no pertenecer administrativamente a Benidorm, presenta un abandono urbanístico, ya que no es incluida en los desarrollos y planes de ordenación urbana de éste. Además, las características fisiográficas de la cala, hacen que se encuentre en un continuo proceso de transporte de sedimentos del extremo Nor-Este al Sur-Este. Por lo tanto, las necesidades primarias del emplazamiento serán la regeneración ambiental del entorno de la cala, para garantizar el desarrollo de la fauna marina y, en especial, de la flora mediterránea, la cual es única en el Mundo y se encuentra en situación de protección, ya que han sido gravemente dañadas durante los años de desarrollo urbanístico boyante en el entorno. Se debe proceder a una regeneración litoral que reduzca el transporte de sedimentos, mediante elementos de contención de finos que nos garanticen una playa estable. Dado que la cala tiene una alta ocupación, habrá que aumentar el ancho de playa, especialmente en la zona de erosión. A sí mismo, para paliar los efectos del desarrollo humano en la zona seca, será necesaria la reordenación de la fachada urbana, mediante la construcción de un paseo marítimo que haga de conexión entre la playa y la ciudad, un linde muy abrupto en la actualidad. Todas estas medidas deberán realizarse en un régimen de interacción con el paisaje, para provocar el menor impacto ambiental y visual posible.

Consideraciones sobre el cálculo sísmico de presas-bóveda

En esta comunicación se ponen de manifiesto las características que debe reunir un modelo para el cálculo sísmico de presas-bóveda. Se analizan los factores locales que influyen sobre la solicitación sísmica que actúa en una ubicación concreta. Se destaca la importancia que en el comportamiento sísmico de estas construcciones tiene la interacción entre presa, suelo y embalse los cuales constituyen un único sistema dinámico que se ve afectado por factores tales como: la geometría de la presa y el vaso, las propiedades del suelo, la posible falta de homogeneidad subterránea y la presencia de sedimentos de fondo. A la luz de los anteriores factores y en general de las características propias de las presas-bóveda, el suelo y los embalses, se discuten las posibilidades que ofrecen y los inconvenientes que tienen algunos de los modelos de Elementos Finitos existentes. Se presenta un modelo de Elementos de Contorno que, por las características propias de este método, resulta muy adecuado.

Alteración del comportamiento hidráulico natural de la vertiente noreste de la marisma de El Rocío por efectos antrópicos

El geo-dinamismo torrencial es resultado de la incidencia de precipitaciones torrenciales en la cuenca hidrográfica, y se manifiesta por una intensificación del ciclo del agua (con la presencia de mayores volúmenes de escorrentía y la consiguiente aparición de repentinos y elevados caudales de avenida) e intensificación del ciclo de los sedimentos (mayor erosión del suelo, incremento del transporte de sedimentos en los cauces de drenaje y la formación de grades depósitos de sedimentos en el área dominada de la cuenca). Estos procesos naturales se pueden agravar si se realizan en la cuenca actuaciones que alteren sustancialmente su morfología; es el caso del arroyo del Partido, tributario a la marisma de El Rocío en el Parque Nacional de Doñana, cuyo encauzamiento del tramo final sobre su antiguo cono de sedimentación fue el causante de la formación de un nuevo cono de sedimentación sobre la marisma. Pero las intervenciones antrópicas también pueden dirigir el sistema hacia la estabilidad sedimentaria y el tramo final del arroyo del Partido ha experimentado asimismo actuaciones de carácter restaurador. El antiguo cono de sedimentación del arroyo del Partido es un lugar complejo en el que drenan también otros tres cursos que intercambian sus flujos con dicho arroyo durante las avenidas. Esta Memoria analiza cuatro situaciones del paraje objeto de estudio: la primera responde al estado original del lugar en cuestión (situación A); la segunda a la creada tras el encauzamiento del arroyo del Partido en 1981 (situación B); la tercera a la formada pasado un tiempo desde el encauzamiento, cuando el geo-dinamismo del lugar hizo sus efectos sobre el paraje en cuestión (situación C); y la cuarta la generada tras la restauración en el verano de 2006 del régimen hidráulico-sedimentario del arroyo del Partido con la Actuación núm. 3 del Proyecto Doñana 2005. Se han investigado los efectos que la alteración de la morfología del sistema de drenaje del antiguo cono de sedimentación del arroyo del Partido causa en los parámetros hidráulicos que desencadenan el geo-dinamismo torrencial. Para ello se han reproducido las condiciones morfológicas de las cuatro situaciones más representativas ya comentadas que ha experimentado el paraje en los últimos 30 años y que han condicionado el comportamiento de las avenidas torrenciales ocurridas en dicho periodo. En este contexto se ha elaborado un modelo hidrológico con el programa HECHMS a partir de la información disponible y un modelo hidrodinámico en 2D con el programa IBER; con este último se han comparado ocho escenarios característicos derivados de la incidencia de dos avenidas representativas de los eventos torrenciales de la zona sobre las morfologías contempladas. También se ha estudiado la tensión de arrastre de las avenidas dada su incidencia en el terreno sobre el que transitan, analizando de manera separada los efectos en los cuatro cauces de drenaje del antiguo cono de sedimentación y del resto de la superficie del mismo. La investigación realizada ha permitido describir minuciosamente el comportamiento de los drenajes locales del antiguo y complejo cono de sedimentación que bordea el noroeste de la marisma del Parque Nacional de Doñana, desde la marisma de El Rocío hasta la marisma gallega, lo que redundará en beneficio de la gestión del Parque Nacional de Doñana y de su entorno, y especialmente en la seguridad del perímetro oriental de la aldea de El Rocío. ABSTRACT The torrential geo-dynamism results from the effect of torrential precipitations on the hydrological watershed, and it shows ups with the intensification of the water cycle (with presence of higher volumes of run-off and the consequent appearance of sudden and high food flows) and intensification of the sediment cycle (higher soil erosion, increase of sediment transport in the drainage courses and the formation of big sediment deposits in the area downstream the watershed). These natural processes can worsen if actions that change the morphology of the watershed are carried out. It is the case of the Partido stream, tributary of the marsh of El Rocío in the Doñana National Park, whose channelization of the final stretch over its old alluvial fan was the cause of the formation of a new alluvial fan over the marsh. But human actions can also go toward the sediment stabilization, and the final stretch of the Partido stream has experienced restoration actions as well. The old alluvial fan of the Partido stream is a complex terrain where also other three courses exchange their flows during floods. This document analyses four situations of the area studied: the first corresponds to the original state of the place (situation A); the second is the state created after the channel of the Partido stream was built in 1981 (situation B); the third is the state formed after some time passed and the geo-dynamism caused effect in the channelization (situation C); and the forth situation is the one created after the restoration in summer of 2006 of the hydraulic-sedimentary regime of the Partido stream with de Action number 3 of the Doñana 2005 Project. It has been researched the effects that the alteration of the morphology of the drainage system of the old alluvial fan of the Partido stream causes to the hydraulic parameters that trigger the torrential geo-dynamism. For that, the morphologic conditions of the four more representative situations were reproduced. Those situations experienced in the area in the last 30 years have conditioned the behaviour of torrential floods in that period. In this context, a hydrologic model with HEC-HMS and a 2D hydro-dynamic model with IBER have been set up. With IBER, eight scenes have been compared resulting from the impact of two representative floods from torrential events in the region against the morphologies considered. Also, the shear stress of the floods on the terrain has been studied, analysing separately the effect in the four drainage courses of the old alluvial fan and in the rest of its surfaces. The research carried out has allowed describing in minute detail the behavior of the local drainages of the old and complex alluvial fan that border the northwest part of the marsh of Doñana National Park, from the marsh of El Rocío to the gallega marsh. This will result in benefit for the management of Doñana National Park and its surroundings, and especially for the safety of the eastern limit of El Rocío village.

Residuos orgánicos como fuentes de fósforo

En este trabajo, se han llevado a cabo distintos experimentos en laboratorio, con el objetivo de estudiar el efecto de la aplicación de residuos orgánicos como fuentes de P en las pérdidas de este elemento, que se producen en suelo, tanto por escorrentía superficial como por lixiviación. El interés por evaluar las pérdidas de P se debe a la necesidad de conocer mejor los factores que influyen en los procesos de pérdidas de este elemento y así, poder reducir los problemas de eutrofización de aguas, tanto superficiales como subterráneas, provocadas por un exceso de este nutriente, junto con otros como el nitrógeno. Los trabajos experimentales que se han llevado a cabo se detallan a continuación: Se ha realizado el estudio de las formas de P contenidas en una serie de 14 residuos orgánicos, de distinto origen y tratamiento (compost, lodos, purines y digestato), comparando la información aportada por tres protocolos de fraccionamientos de P diferentes, seleccionados entre los principales métodos utilizados: protocolo de fraccionamiento de Ruttemberg (1992), protocolo de Normas, medidas y ensayos (Ruban et al., 2001a) y protocolo de Huang et al. (2008). Todos los métodos de fraccionamiento empleados aportaron información útil para conocer las formas de P de los residuos, a pesar de que alguno de ellos fue descrito para sedimentos o suelos. Sin embargo, resulta difícil comparar los resultados entre unos y otros, ya que cada uno emplea extractantes y tiempos de ensayos diferentes. Las cantidades de P total determinadas por cada método mantienen una relación lineal, aunque el método SMT, por ser más directo, obtiene las cantidades más elevadas para todos los residuos. Los métodos secuenciales (métodos de Huang y Ruttemberg), a pesar de ser más largos y tediosos, aportan información más detallada de la disponibilidad de las formas de P, y con ello, permiten obtener una mejor estimación de las pérdidas potenciales de este elemento tras su aplicación en suelo. Se han encontrado relaciones positivas entre las fracciones determinadas por estos dos métodos. Así mismo, se encuentra una correlación entre las formas solubles de P y la concentración de los iones [Ca + Fe + Al] de los residuos, útiles como indicadores de la disponibilidad de este elemento. Sin embargo, el protocolo SMT, no obtiene información de la solubilidad o disponibilidad de las formas de P contenidas, ni las formas determinadas mantienen relaciones directas con las de los otros métodos, con excepción del P total. Para el estudio del comportamiento de los residuos aplicados en suelos, se pusieron a punto sistemas de lluvia simulada, con el objetivo de caracterizar las pérdidas de P en la escorrentía superficial generada. Por otra parte, se emplearon columnas de suelos enmendados con residuos orgánicos, para el estudio de las pérdidas de P por lixiviación. Los ensayos de simulación de lluvia se llevaron a cabo de acuerdo al “National Phosphorus Research proyect“ (2001), que consigue simular eventos sucesivos de lluvia en unas condiciones semejantes a la realidad, empleando cajas llenas de suelo del horizonte superficial, con residuos aplicados tanto superficialmente como mediante mezcla con el propio suelo. Los ensayos se realizaron con seis residuos de diferente naturaleza y sometidos a distintos tratamientos. Se encontraron diferencias significativas en las pérdidas de las formas de P analizadas, tanto disueltas como particuladas, en las aguas de escorrentía generadas. En general, las pérdidas en el primer evento de lluvia tras la aplicación de los residuos fueron mayores a las generadas en el segundo evento, predominando las formas de P particuladas respecto a las disueltas en ambos. Se encontró una relación positiva entre las pérdidas de P en las aguas de escorrentía generadas en cada ensayo, con los contenidos de P soluble en agua y fácilmente disponible de los residuos empleados, determinados por los protocolos de fraccionamientos secuenciales. Además, se emplearon los modelos matemáticos desarrollados por Vadas et al. (2005, 2007), de evaluación de las pérdidas de P por escorrentía para fertilizantes y estiércoles. La predicción de estos modelos no se cumple en el caso de todos los residuos. Las distintas propiedades físicas de los residuos pueden afectar a las diferencias entre las pérdidas experimentales y las esperadas. Los ensayos de simulación del proceso de lixiviación se llevaron a cabo en columnas de percolación, con suelos enmendados con residuos orgánicos, de acuerdo a la norma “CEN/TS 14405–2004: Caracterización de los residuos – Test de comportamiento de lixiviación – Test de flujo ascendente”. Las pérdidas de P por procesos de lixiviación de agua, han sido despreciadas durante mucho tiempo respecto a las pérdidas por escorrentía. Sin embargo, se ha demostrado que deben tenerse en consideración, principalmente en algunos tipos de suelos o zonas cercanas a acuíferos. Se utilizaron tres suelos de distinta procedencia para los ensayos, de manera que se pudo estudiar la influencia del tipo de suelo en las pérdidas de P para cada tipo de residuo (purín, compost, digestato y lodo de EDAR). Los índices de adsorción de P determinados para cada suelo permiten evaluar aquellos que presentarán más riesgo de producir pérdidas de este elemento al aplicarse fuentes externas de P, encontrando una relación positiva entre ambos. Las pérdidas de P en los lixiviados varían en función tanto del residuo como del suelo empleado. Para el compost, el purín y el lodo, se encontró una relación entre las pérdidas generadas en el agua lixiviada de las columnas y las formas de P soluble contenidas en los residuos. Sin embargo, en el caso del digestato, no existía esta correlación. Las pérdidas para este residuo fueron en todos los casos menores a las estimadas, considerando las formas de P contenido. El estudio de la mojabilidad, propiedad física del residuo que evalúa la capacidad de interacción residuo-agua, permitió explicar el comportamiento anómalo de este residuo, con una mayor resistencia a que el agua entrara en su estructura y por tanto, una mayor dificultad de solubilizar el P contenido en el propio residuo, que en el caso de otros residuos. En general, podemos considerar que el estudio de las formas de P más disponibles o solubles en agua, aporta información útil de las pérdidas potenciales de P. Sin embargo, es necesario estudiar las propiedades físicas de los residuos orgánicos aplicados y la capacidad de adsorción de P de los suelos, para estimar las pérdidas de P y con ello, colaborar a controlar los procesos de eutrofización en aguas. ABSTRACT This dissertation explores the effect of organic wastes application as sources of P in losses of this element that occur by both surface runoff and leaching in soil. To do so, diverse laboratory experiments are conducted and presented here. Evaluating P losses is necessary to better understand the factors that influence the processes behind the loss of this element. Reducing P losses reduces eutrophication problems of both surface water and groundwater caused by an excess of this nutrient, along with other as nitrogen. Details of the experiments are presented below: The first experiment studies the forms of P contained in a series of 14 organic wastes of different origin and treatment (compost, sludge, slurry and digestate), comparing the information provided by three methods of P fractionation. The methods selected were: Ruttemberg protocol (1992); Standards, Measurements and Testing protocol (Ruban et al., 2001a); and Huang protocol (Huang et al., 2008). All fractionation methods employed successfully contribute to our knowledge of P forms in wastes, even though one of them was originally described for sediments or soils information. However, it is difficult to compare results among each other, as each protocol employs different extractants and time in the trials. Total amounts of P obtained by each method show a linear relationship, although the SMT method, which is more direct, obtains the highest amounts for all residues. Sequential methods (Huang and Ruttemberg’s protocols), despite being longer and more tedious, provide more detailed information on the availability of the forms of P. Therefore, allow the estimation of the potential losses of P after application in soil. Furthermore, positive relationships have been found among fractions obtained by these methods. Positive relationship has been found also among soluble forms of P and the concentration of ions Fe + Ca + Al, which is useful as an indicator of the availability of this element. However, the SMT protocol does not collect information about solubility or availability of forms of P contained; neither do certain forms maintain direct relations with the forms from other methods, with the exception of total P methods. To study the behavior of wastes applied to soils two experiments were conducted. Simulated rain systems were prepared to characterize P losses in the surface runoff generated. In addition, columns of soils amended with organic waste were developed for the study of P leaching losses. Simulated rain systems were carried out according to the ’National Phosphorus Research Project’ (2001), which manages to simulate successive rainfall events in conditions resembling reality. The experiment uses boxes filled with soil from the surface horizon amended with residues, both superficially and by mixing with the soil. Tests were conducted with six residues of different type and subjected to diverse treatments. Findings show significant differences in losses of the P forms analyzed in the generated runoff water, in both solution and particulate forms. In general, losses in the first rainfall event after application of waste were higher than the losses generated in the second event, predominating particulate forms of P over dissolved forms in both events. In all trials, a positive relationship was found between various P forms determined by sequential fractionation protocols (water soluble P and readily available P forms) and P losses in runoff. Furthermore, results from Vadas´s mathematical models (Vadas et al., 2005; 2007) to assess P losses by runoff fertilizers and manures indicate that the prediction of this model is not fulfilled in the case of all residues. The diverse physical properties of wastes may affect the differences between experimental and expected losses. Finally, leaching simulation processes were carried out in percolation columns, filled with soils amended with organic wastes, following the ‘CEN/TS 14405-2004 standard: Characterization of waste - Leaching behavior test - Test Flow ascending ’. P losses by leaching have been neglected for a long time with respect to runoff losses. However, findings corroborate previous studies showing that these P losses have to be taken into account, especially in certain types of soils and in zones near aquifers. To study the influence of soil type on P losses, experiments were carried out with three different soils and for each type of waste (manure, compost, digestate and sludge WWTP). Each soil’s P adsorption rates allow assessing which soils imply a higher risk of P losses when external sources of P are applied. P losses in leachate vary according to the type of soil employed and according to the specific residue. In the case of compost, manure and sludge, there is a relationship between leaching losses and residues’ soluble forms of P. The exception being the digestate, where there was no such correlation. Digestate P losses by leaching were lower than expected in all cases considering the forms of P contained. Moreover, examining digestate wettability -- that is, the physical property of the residue that assesses the capacity of waste-water interaction -- allowed explaining the anomalous behavior of this residue. Digestate has a high resistance to water entering its structure and thus higher difficulty to solubilize the P contained. Overall, studying the more available or soluble P forms provides useful information about the potential loss of P. However, this dissertation shows that it is necessary to examine the physical properties of organic residues applied as well as the P adsorption capacity of soils to estimate P losses, and thus to control eutrophication in water.

Los efectos del bosque en las cuencas de montaña ante la incidencia en ellas de eventos torrenciales: Torrente de Arratiecho en el Pirineo Aragonés

Este documento describe sintéticamente el fenómeno del geo-dinamismo torrencial que se desencadena en las cuencas hidrográficas, cuando sobre ellas inciden precipitaciones torrenciales extraordinarias, que intensifican en las mismas los ciclos del agua y de los sedimentos; especialmente cuando se trata de cuencas de montaña. Se comenta a continuación la importancia del bosque como amortiguador de dicho geo-dinamismo y su repercusión en los trabajos de corrección de torrentes de montaña, a través de la rehabilitación forestal de sus cuencas vertientes. Finalmente, se ilustra con un ejemplo la conveniencia de llevar a cabo la restauración hidrológico-forestal de las cuencas de montaña que lo necesitan. Se trata de la restauración efectuada entre 1901-05 al torrente de Arratiecho en el Pirineo de Huesca.