Biogeochemical indicators of elevated nitrogen deposition in semiarid Mediterranean ecosystems

Nitrogen (N) deposition has doubled the natural N inputs received by ecosystems through biological N fixation and is currently a global problem that is affecting the Mediterranean regions. We evaluated the existing relationships between increased atmospheric N deposition and biogeochemical indicators related to soil chemical factors and cryptogam species across semiarid central, southern, and eastern Spain. The cryptogam species studied were the biocrust-forming species Pleurochaete squarrosa (moss) and Cladonia foliacea (lichen). Sampling sites were chosen in Quercus coccifera (kermes oak) shrublands and Pinus halepensis (Aleppo pine) forests to cover a range of inorganic N deposition representative of the levels found in the Iberian Peninsula (between 4.4 and 8.1 kg N ha(-1) year(-1)). We extended the ambient N deposition gradient by including experimental plots to which N had been added for 3 years at rates of 10, 20, and 50 kg N ha(-1) year(-1). Overall, N deposition (extant plus simulated) increased soil inorganic N availability and caused soil acidification. Nitrogen deposition increased phosphomonoesterase (PME) enzyme activity and PME/nitrate reductase (NR) ratio in both species, whereas the NR activity was reduced only in the moss. Responses of PME and NR activities were attributed to an induced N to phosphorus imbalance and to N saturation, respectively. When only considering the ambient N deposition, soil organic C and N contents were positively related to N deposition, a response driven by pine forests. The PME/NR ratios of the moss were better predictors of N deposition rates than PME or NR activities alone in shrublands, whereas no correlation between N deposition and the lichen physiology was observed. We conclude that integrative physiological measurements, such as PME/NR ratios, measured on sensitive species such as P. squarrosa, can provide useful data for national-scale biomonitoring programs, whereas soil acidification and soil C and N storage could be useful as additional corroborating ecosystem indicators of chronic N pollution.

Guidelines for the design of an environmental flow regime program at the regional scale: Application of the ELOHA framework in the Ebro River basin

La gestión de los recursos hídricos se convierte en un reto del presente y del futuro frente a un panorama de continuo incremento de la demanda de agua debido al crecimiento de la población, el crecimiento del desarrollo económico y los posibles efectos del calentamiento global. La política hidráulica desde los años 60 en España se ha centrado en la construcción de infraestructuras que han producido graves alteraciones en el régimen natural de los ríos. Estas alteraciones han provocado y acrecentado los impactos sobre los ecosistemas fluviales y ribereños. Desde los años 90, sin embargo, ha aumentado el interés de la sociedad para conservar estos ecosistemas. El concepto de caudales ambientales consiste en un régimen de caudales que simula las características principales del régimen natural. Los caudales ambientales están diseñados para conservar la estructura y funcionalidad de los ecosistemas asociados al régimen fluvial, bajo la hipótesis de que los elementos que conforman estos ecosistemas están profundamente adaptados al régimen natural de caudales, y que cualquier alteración del régimen natural puede provocar graves daños a todo el sistema. El método ELOHA (Ecological Limits of Hydrological Alteration) tiene como finalidad identificar las componentes del régimen natural de caudales que son clave para mantener el equilibrio de los ecosistemas asociados, y estimar los límites máximos de alteración de estas componentes para garantizar su buen estado. Esta tesis presenta la aplicación del método ELOHA en la cuenca del Ebro. La cuenca del Ebro está profundamente regulada e intervenida por el hombre, y sólo las cabeceras de los principales afluentes del Ebro gozan todavía de un régimen total o cuasi natural. La tesis se estructura en seis capítulos que desarrollan las diferentes partes del método. El primer capítulo explica cómo se originó el concepto “caudales ambientales” y en qué consiste el método ELOHA. El segundo capítulo describe el área de estudio. El tercer capítulo realiza una clasificación de los regímenes naturales de la cuenca (RNC) del Ebro, basada en series de datos de caudal mínimamente alterado y usando exclusivamente parámetros hidrológicos. Se identificaron seis tipos diferentes de régimen natural: pluvial mediterráneo, nivo-pluvial, pluvial mediterréaneo con una fuerte componente del caudal base, pluvial oceánico, pluvio-nival oceánico y Mediterráneo. En el cuarto capítulo se realiza una regionalización a toda la cuenca del Ebro de los seis RNC encontrados en la cueca. Mediante parámetros climáticos y fisiográficos se extrapola la información del tipo de RNC a puntos donde no existen datos de caudal inalterado. El patrón geográfico de los tipos de régimen fluvial obtenido con la regionalización resultó ser coincidente con el patrón obtenido a través de la clasificación hidrológica. El quinto capítulo presenta la validación biológica de los procesos de clasificación anteriores: clasificación hidrológica y regionalización. La validación biológica de los tipos de regímenes fluviales es imprescindible, puesto que los diferentes tipos de régimen fluvial van a servir de unidades de gestión para favorecer el mantenimiento de los ecosistemas fluviales. Se encontraron diferencias significativas entre comunidades biológicas en cinco de los seis tipos de RNC encontrados en la cuenca. Finalmente, en el sexto capítulo se estudian las relaciones hidro-ecológicas existentes en tres de los seis tipos de régimen fluvial encontrados en la cuenca del Ebro. Mediante la construcción de curvas hidro-ecológicas a lo largo de un gradiente de alteración hidrológica, se pueden sugerir los límites de alteración hidrológica (ELOHAs) para garantizar el buen estado ecológico en cada uno de los tipos fluviales estudiados. Se establecieron ELOHAs en tres de los seis tipos de RNC de la cuenca del Ebro Esta tesis, además, pone en evidencia la falta de datos biológicos asociados a registros de caudal. Para llevar a cabo la implantación de un régimen de caudales ambientales en la cuenca, la ubicación de los puntos de muestreo biológico cercanos a estaciones de aforo es imprescindible para poder extraer relaciones causa-efecto de la gestión hidrológica sobre los ecosistemas dependientes. ABSTRACT In view of a growing freshwater demand because of population raising, improvement of economies and the potential effects of climate change, water resources management has become a challenge for present and future societies. Water policies in Spain have been focused from the 60’s on constructing hydraulic infrastructures, in order to dampen flow variability and granting water availability along the year. Consequently, natural flow regimes have been deeply altered and so the depending habitats and its ecosystems. However, an increasing acknowledgment of societies for preserving healthy freshwater ecosystems started in the 90’s and agreed that to maintain healthy freshwater ecosystems, it was necessary to set environmental flow regimes based on the natural flow variability. The Natural Flow Regime paradigm (Richter et al. 1996, Poff et al. 1997) bases on the hypothesis that freshwater ecosystems are made up by elements adapted to natural flow conditions, and any change on these conditions can provoke deep impacts on the whole system. Environmental flow regime concept consists in designing a flow regime that emulates natural flow characteristics, so that ecosystem structure, functions and services are maintained. ELOHA framework (Ecological Limits of Hydrological Alteration) aims to identify key features of the natural flow regime (NFR) that are needed to maintain and preserve healthy freshwater and riparian ecosystems. Moreover, ELOHA framework aims to quantify thresholds of alteration of these flow features according to ecological impacts. This thesis describes the application of the ELOHA framework in the Ebro River Basin. The Ebro River basin is the second largest basin in Spain and it is highly regulated for human demands. Only the Ebro headwaters tributaries still have completely unimpaired flow regime. The thesis has six chapters and the process is described step by step. The first chapter makes an introduction to the origin of the environmental flow concept and the necessity to come up. The second chapter shows a description of the study area. The third chapter develops a classification of NFRs in the basin based on natural flow data and using exclusively hydrological parameters. Six NFRs were found in the basin: continental Mediterranean-pluvial, nivo-pluvial, continental Mediterranean pluvial (with groundwater-dominated flow pattern), pluvio-oceanic, pluvio-nival-oceanic and Mediterranean. The fourth chapter develops a regionalization of the six NFR types across the basin by using climatic and physiographic variables. The geographical pattern obtained from the regionalization process was consistent with the pattern obtained with the hydrologic classification. The fifth chapter performs a biological validation of both classifications, obtained from the hydrologic classification and the posterior extrapolation. When the aim of flow classification is managing water resources according to ecosystem requirements, a validation based on biological data is compulsory. We found significant differences in reference macroinvertebrate communities between five over the six NFR types identified in the Ebro River basin. Finally, in the sixth chapter we explored the existence of significant and explicative flow alteration-ecological response relationships (FA-E curves) within NFR types in the Ebro River basin. The aim of these curves is to find out thresholds of hydrological alteration (ELOHAs), in order to preserve healthy freshwater ecosystem. We set ELOHA values in three NFR types identified in the Ebro River basin. During the development of this thesis, an inadequate biological monitoring in the Ebro River basin was identified. The design and establishment of appropriate monitoring arrangements is a critical final step in the assessment and implementation of environmental flows. Cause-effect relationships between hydrology and macroinvertebrate community condition are the principal data that sustain FA-E curves. Therefore, both data sites must be closely located, so that the effects of external factors are minimized. The scarce hydro-biological pairs of data available in the basin prevented us to apply the ELOHA method at all NFR types.

A model for turbulent combustion simulation of large scale hydrogen explosions

El agotamiento, la ausencia o, simplemente, la incertidumbre sobre la cantidad de las reservas de combustibles fósiles se añaden a la variabilidad de los precios y a la creciente inestabilidad en la cadena de aprovisionamiento para crear fuertes incentivos para el desarrollo de fuentes y vectores energéticos alternativos. El atractivo de hidrógeno como vector energético es muy alto en un contexto que abarca, además, fuertes inquietudes por parte de la población sobre la contaminación y las emisiones de gases de efecto invernadero. Debido a su excelente impacto ambiental, la aceptación pública del nuevo vector energético dependería, a priori, del control de los riesgos asociados su manipulación y almacenamiento. Entre estos, la existencia de un innegable riesgo de explosión aparece como el principal inconveniente de este combustible alternativo. Esta tesis investiga la modelización numérica de explosiones en grandes volúmenes, centrándose en la simulación de la combustión turbulenta en grandes dominios de cálculo en los que la resolución que es alcanzable está fuertemente limitada. En la introducción, se aborda una descripción general de los procesos de explosión. Se concluye que las restricciones en la resolución de los cálculos hacen necesario el modelado de los procesos de turbulencia y de combustión. Posteriormente, se realiza una revisión crítica de las metodologías disponibles tanto para turbulencia como para combustión, que se lleva a cabo señalando las fortalezas, deficiencias e idoneidad de cada una de las metodologías. Como conclusión de esta investigación, se obtiene que la única estrategia viable para el modelado de la combustión, teniendo en cuenta las limitaciones existentes, es la utilización de una expresión que describa la velocidad de combustión turbulenta en función de distintos parámetros. Este tipo de modelos se denominan Modelos de velocidad de llama turbulenta y permiten cerrar una ecuación de balance para la variable de progreso de combustión. Como conclusión también se ha obtenido, que la solución más adecuada para la simulación de la turbulencia es la utilización de diferentes metodologías para la simulación de la turbulencia, LES o RANS, en función de la geometría y de las restricciones en la resolución de cada problema particular. Sobre la base de estos hallazgos, el crea de un modelo de combustión en el marco de los modelos de velocidad de la llama turbulenta. La metodología propuesta es capaz de superar las deficiencias existentes en los modelos disponibles para aquellos problemas en los que se precisa realizar cálculos con una resolución moderada o baja. Particularmente, el modelo utiliza un algoritmo heurístico para impedir el crecimiento del espesor de la llama, una deficiencia que lastraba el célebre modelo de Zimont. Bajo este enfoque, el énfasis del análisis se centra en la determinación de la velocidad de combustión, tanto laminar como turbulenta. La velocidad de combustión laminar se determina a través de una nueva formulación capaz de tener en cuenta la influencia simultánea en la velocidad de combustión laminar de la relación de equivalencia, la temperatura, la presión y la dilución con vapor de agua. La formulación obtenida es válida para un dominio de temperaturas, presiones y dilución con vapor de agua más extenso de cualquiera de las formulaciones previamente disponibles. Por otra parte, el cálculo de la velocidad de combustión turbulenta puede ser abordado mediante el uso de correlaciones que permiten el la determinación de esta magnitud en función de distintos parámetros. Con el objetivo de seleccionar la formulación más adecuada, se ha realizado una comparación entre los resultados obtenidos con diversas expresiones y los resultados obtenidos en los experimentos. Se concluye que la ecuación debida a Schmidt es la más adecuada teniendo en cuenta las condiciones del estudio. A continuación, se analiza la importancia de las inestabilidades de la llama en la propagación de los frentes de combustión. Su relevancia resulta significativa para mezclas pobres en combustible en las que la intensidad de la turbulencia permanece moderada. Estas condiciones son importantes dado que son habituales en los accidentes que ocurren en las centrales nucleares. Por ello, se lleva a cabo la creación de un modelo que permita estimar el efecto de las inestabilidades, y en concreto de la inestabilidad acústica-paramétrica, en la velocidad de propagación de llama. El modelado incluye la derivación matemática de la formulación heurística de Bauwebs et al. para el cálculo de la incremento de la velocidad de combustión debido a las inestabilidades de la llama, así como el análisis de la estabilidad de las llamas con respecto a una perturbación cíclica. Por último, los resultados se combinan para concluir el modelado de la inestabilidad acústica-paramétrica. Tras finalizar esta fase, la investigación se centro en la aplicación del modelo desarrollado en varios problemas de importancia para la seguridad industrial y el posterior análisis de los resultados y la comparación de los mismos con los datos experimentales correspondientes. Concretamente, se abordo la simulación de explosiones en túneles y en contenedores, con y sin gradiente de concentración y ventilación. Como resultados generales, se logra validar el modelo confirmando su idoneidad para estos problemas. Como última tarea, se ha realizado un analisis en profundidad de la catástrofe de Fukushima-Daiichi. El objetivo del análisis es determinar la cantidad de hidrógeno que explotó en el reactor número uno, en contraste con los otros estudios sobre el tema que se han centrado en la determinación de la cantidad de hidrógeno generado durante el accidente. Como resultado de la investigación, se determinó que la cantidad más probable de hidrogeno que fue consumida durante la explosión fue de 130 kg. Es un hecho notable el que la combustión de una relativamente pequeña cantidad de hidrogeno pueda causar un daño tan significativo. Esta es una muestra de la importancia de este tipo de investigaciones. Las ramas de la industria para las que el modelo desarrollado será de interés abarca la totalidad de la futura economía de hidrógeno (pilas de combustible, vehículos, almacenamiento energético, etc) con un impacto especial en los sectores del transporte y la energía nuclear, tanto para las tecnologías de fisión y fusión. ABSTRACT The exhaustion, absolute absence or simply the uncertainty on the amount of the reserves of fossil fuels sources added to the variability of their prices and the increasing instability and difficulties on the supply chain are strong incentives for the development of alternative energy sources and carriers. The attractiveness of hydrogen in a context that additionally comprehends concerns on pollution and emissions is very high. Due to its excellent environmental impact, the public acceptance of the new energetic vector will depend on the risk associated to its handling and storage. Fromthese, the danger of a severe explosion appears as the major drawback of this alternative fuel. This thesis investigates the numerical modeling of large scale explosions, focusing on the simulation of turbulent combustion in large domains where the resolution achievable is forcefully limited. In the introduction, a general description of explosion process is undertaken. It is concluded that the restrictions of resolution makes necessary the modeling of the turbulence and combustion processes. Subsequently, a critical review of the available methodologies for both turbulence and combustion is carried out pointing out their strengths and deficiencies. As a conclusion of this investigation, it appears clear that the only viable methodology for combustion modeling is the utilization of an expression for the turbulent burning velocity to close a balance equation for the combustion progress variable, a model of the Turbulent flame velocity kind. Also, that depending on the particular resolution restriction of each problem and on its geometry the utilization of different simulation methodologies, LES or RANS, is the most adequate solution for modeling the turbulence. Based on these findings, the candidate undertakes the creation of a combustion model in the framework of turbulent flame speed methodology which is able to overcome the deficiencies of the available ones for low resolution problems. Particularly, the model utilizes a heuristic algorithm to maintain the thickness of the flame brush under control, a serious deficiency of the Zimont model. Under the approach utilized by the candidate, the emphasis of the analysis lays on the accurate determination of the burning velocity, both laminar and turbulent. On one side, the laminar burning velocity is determined through a newly developed correlation which is able to describe the simultaneous influence of the equivalence ratio, temperature, steam dilution and pressure on the laminar burning velocity. The formulation obtained is valid for a larger domain of temperature, steam dilution and pressure than any of the previously available formulations. On the other side, a certain number of turbulent burning velocity correlations are available in the literature. For the selection of the most suitable, they have been compared with experiments and ranked, with the outcome that the formulation due to Schmidt was the most adequate for the conditions studied. Subsequently, the role of the flame instabilities on the development of explosions is assessed. Their significance appears to be of importance for lean mixtures in which the turbulence intensity remains moderate. These are important conditions which are typical for accidents on Nuclear Power Plants. Therefore, the creation of a model to account for the instabilities, and concretely, the acoustic parametric instability is undertaken. This encloses the mathematical derivation of the heuristic formulation of Bauwebs et al. for the calculation of the burning velocity enhancement due to flame instabilities as well as the analysis of the stability of flames with respect to a cyclic velocity perturbation. The results are combined to build a model of the acoustic-parametric instability. The following task in this research has been to apply the model developed to several problems significant for the industrial safety and the subsequent analysis of the results and comparison with the corresponding experimental data was performed. As a part of such task simulations of explosions in a tunnel and explosions in large containers, with and without gradient of concentration and venting have been carried out. As a general outcome, the validation of the model is achieved, confirming its suitability for the problems addressed. As a last and final undertaking, a thorough study of the Fukushima-Daiichi catastrophe has been carried out. The analysis performed aims at the determination of the amount of hydrogen participating on the explosion that happened in the reactor one, in contrast with other analysis centered on the amount of hydrogen generated during the accident. As an outcome of the research, it was determined that the most probable amount of hydrogen exploding during the catastrophe was 130 kg. It is remarkable that the combustion of such a small quantity of material can cause tremendous damage. This is an indication of the importance of these types of investigations. The industrial branches that can benefit from the applications of the model developed in this thesis include the whole future hydrogen economy, as well as nuclear safety both in fusion and fission technology.

Experimental study of the thermal-energy performance of an insulated vegetal facade under summer conditions in a continental mediterranean climate

The use of vegetal systems in facades affects the reduction of the buildings' energy demand, the attenuation of the urban heat island (UHI) and the filtration of pollutants present in the air. Even so, up to now the knowledge about the effect of this type of systems on the thermal performance of insulated facades is limited. This article presents the results of an experimental study carried out in a vegetal facade located in a continental Mediterranean climate zone. The objective is to study the effect of a vegetal finishing, formed by plants and substrate, on the thermal-energy performance of an insulated facade under summer conditions. To this effect, the thermal data obtained from two full-scale experimental mock-ups of the same dimensions and composition of the enclosure and only different in the south facade's enclosure where one incorporates a vegetation layer are compared and analysed. The results show that, in spite of the high thermal resistance of the enclosure, the effect of the vegetation is very positive, particularly in the warmer hours of the day. Therefore, vegetal facades can be used as a passive cooling strategy, reducing the consumption of energy for refrigeration and improving the comfort conditions of the users. (C) 2014 Elsevier Ltd. All rights reserved.

Polychaete distribution pattern on the Valencian Community coast, Spanish Mediterranean

Changes in benthic community structure are strongly related to environmental factors, and we need to determine how these natural changes occur in order to interpret the possible changes associated with anthropogenic impacts. The aim of this survey was to characterize and classify the polychaete assemblages inhabiting unpolluted soft bottoms in the Spanish Mediterranean in relation to environmental factors. Thirteen localities were sampled at depths between 9 and 31 m, from 2004 to 2006. Multivariate techniques showed that the structure of polychaete assemblages detected in 2004 was consistent over time and correlations between polychaetes and environmental factors were detected. The study area comprises four kinds of communities mainly characterized by polychaete assemblages, sediment types, and depth.

Cyclically-arranged, storm-controlled, prograding lithosomes in Messinian terrigenous shelves (Bajo Segura Basin, western Mediterranean)

This work focuses on a Messinian shallow-marine terrigenous unit, termed the La Virgen Formation, which forms part of the sedimentary infill of the Bajo Segura Basin (Betic margin of the western Mediterranean). This formation was deposited during a high sea level phase prior to the onset of the Messinian Salinity Crisis. Stratigraphically, it comprises a prograding stack of sandstone lithosomes alternating with marly intervals (1st-order cyclicity). These lithosomes are characterized by a homoclinal geometry that tapers distally, and interfinger with pelagic sediments rich in planktonic and benthic microfauna (Torremendo Formation). An analysis of sedimentary facies of each lithosome reveals a repetitive succession of sandy storm beds (tempestites), occasionally amalgamated, which are separated by thin marly layers (2nd-order cyclicity). Each storm bed contains internal erosional surfaces (3rd-order cyclicity) that delimit sets of laminae. Two categories of storm beds have been differentiated. The first one includes layers formed below storm wave base (SWB), characterized by traction structures associated to unidirectional flows (scoured base, planar lamination, and parting lineation). The second category consists of layers deposited above the SWB which display typical high regime oscillatory flow structures (swaley and hummocky cross lamination). In both cases, the ichnological record is characterized by an oligotypic association of Ophiomorpha nodosa, which can be interpreted as the result of allochthonous tracemakers (crustaceans) transported during storm events together with the sediment. The benthic microfauna in the marly intervals that separate the sandstone lithosomes (1st-order cyclicity) indicates that the storm ebb surges were deposited at depths ranging from those of inner shelf settings (with Elphidium spp. and Cibicides lobatulus) to those of outer shelf (with Valvulineria complanata and Uvigerina cylindrica). At the distal end of the sandstone lithosomes, the planktonic microfauna is characterized by a high content of taxa indicative of warm-oligotrophic waters (Globigerinoides obliquus and Globigerinoides bulloideus). In contrast, in the marly intervals, the microfauna is dominated by species typical of cold-eutrophic waters (Globigerina and Neogloboquadrina). This alternation of planktic foraminiferal assemblages is interpreted as being the expression of climatic cycles, in which every episode of progradation of tempestite-dominated lithosomes corresponds to maximum insolation and warm waters, whereas episodes of marly deposition correspond to minimal insolation and cold waters. The 1st-order cyclicity recorded in the La Virgen Formation, in a context of terrigenous storm-dominated shelf, corresponds to sapropel/homogeneous marl cycles formed in a pelagic basin (Torremendo Fm). These cycles in pelagic sediments are commonplace throughout the Mediterranean during the Messinian and reflect precession orbital changes: repeated periods of maximum insolation – minimum precession (sapropels) and minimal insolation – maximum precession (homogeneous marls). The fact that the example of terrigenous unit studied herein is coetaneous with the well-developed reef complexes in the Mediterranean basins points out the importance of sediment supply in the formation of large-scale sandy lithosomes. This is a crucial aspect to understanding reservoir genesis as well as lateral stratigraphic relationships with potential seal and/or source rocks.

Temporal variation in saproxylic beetle assemblages in a Mediterranean ecosystem

One of the main challenges in biological conservation has been to understand species distribution across space and time. Over the last decades, many diversity and conservation surveys have been conducted that have revealed that habitat heterogeneity acts as a major factor that determines saproxylic assemblages. However, temporal dynamics have been poorly studied, especially in Mediterranean forests. We analyzed saproxylic beetle distribution at inter and intra-annual scales in a “dehesa” ecosystem, which is a traditional Iberian agrosilvopastoral ecosystem that is characterized by the presence of old and scattered trees that dominate the landscape. Significant differences in effective numbers of families/species and species richness were found at the inter-annual scale, but this was not the case for composition. Temperature and relative humidity did not explain these changes which were mainly due to the presence of rare species. At the intra-annual scale, significant differences in the effective numbers of families/species, species richness and composition between seasons were found, and diversity partitioning revealed that season contributed significantly to gamma-diversity. Saproxylic beetle assemblages exhibited a marked seasonality in richness but not in abundance, with two peaks of activity, the highest between May and June, and the second between September and October. This pattern is mainly driven by the seasonality of the climate in the Mediterranean region, which influences ecosystem dynamics and imposes a marked seasonality on insect assemblages. An extended sampling period over different seasons allowed an overview of saproxylic dynamics, and revealed which families/species were restricted to particular seasons. Recognizing that seasons act as a driver in modelling saproxylic beetle assemblages might be a valuable tool in monitoring and for conservation strategies in Mediterranean forests.

Association Patterns in Saproxylic Insect Networks in Three Iberian Mediterranean Woodlands and Their Resistance to Microhabitat Loss

The assessment of the relationship between species diversity, species interactions and environmental characteristics is indispensable for understanding network architecture and ecological distribution in complex networks. Saproxylic insect communities inhabiting tree hollow microhabitats within Mediterranean woodlands are highly dependent on woodland configuration and on microhabitat supply they harbor, so can be studied under the network analysis perspective. We assessed the differences in interacting patterns according to woodland site, and analysed the importance of functional species in modelling network architecture. We then evaluated their implications for saproxylic assemblages’ persistence, through simulations of three possible scenarios of loss of tree hollow microhabitat. Tree hollow-saproxylic insect networks per woodland site presented a significant nested pattern. Those woodlands with higher complexity of tree individuals and tree hollow microhabitats also housed higher species/interactions diversity and complexity of saproxylic networks, and exhibited a higher degree of nestedness, suggesting that a higher woodland complexity positively influences saproxylic diversity and interaction complexity, thus determining higher degree of nestedness. Moreover, the number of insects acting as key interconnectors (nodes falling into the core region, using core/periphery tests) was similar among woodland sites, but the species identity varied on each. Such differences in insect core composition among woodland sites suggest the functional role they depict at woodland scale. Tree hollows acting as core corresponded with large tree hollows near the ground and simultaneously housing various breeding microsites, whereas core insects were species mediating relevant ecological interactions within saproxylic communities, e.g. predation, competitive or facilitation interactions. Differences in network patterns and tree hollow characteristics among woodland sites clearly defined different sensitivity to microhabitat loss, and higher saproxylic diversity and woodland complexity showed positive relation with robustness. These results highlight that woodland complexity goes hand in hand with biotic and ecological complexity of saproxylic networks, and together exhibited positive effects on network robustness.

Coasts, Mediterranean Sea Region, East, 1680-1689 (Raster Image)

This layer is a georeferenced raster image of the historic paper map entitled: Orientalior districtus Maris Mediterranei = T' Ooster gedeelte van de Middelandse Zee. It was published by Fred. de Wit, between 1680 and 1689. Scale [ca. 1:4,500,000]. Covers Mediterranean Sea, Black Sea, and coasts of Europe, North Africa, and the Middle East. Map in Latin and Dutch.The image inside the map neatline is georeferenced to the surface of the earth and fit to the World Miller Cylindrical projected coordinate system. All map collar and inset information is also available as part of the raster image, including any inset maps, profiles, statistical tables, directories, text, illustrations, index maps, legends, or other information associated with the principal map. This map shows features such as drainage, cities and other human settlements, territorial boundaries, shoreline features, and more. Includes index.This layer is part of a selection of digitally scanned and georeferenced historic maps from the Harvard Map Collection. These maps typically portray both natural and manmade features. The selection represents a range of originators, ground condition dates, scales, and map purposes.

Coasts, Mediterranean Sea Region, West, 1680-1689 (Raster Image)

This layer is a georeferenced raster image of the historic paper map entitled: Occidentalior Tractus Maris Mediterranei = Wester gedeelte van de Middelandse Zee. It was published by Fred. de Wit, between 1680 and 1689. Scale [ca. 1:4,500,000]. Scale [ca. 1:4,500,000]. Covers Mediterranean Sea and coasts of Europe and North Africa. Map in Latin and Dutch. The image inside the map neatline is georeferenced to the surface of the earth and fit to the World Miller Cylindrical projected coordinate system. All map collar and inset information is also available as part of the raster image, including any inset maps, profiles, statistical tables, directories, text, illustrations, index maps, legends, or other information associated with the principal map. This map shows features such as drainage, cities and other human settlements, territorial boundaries, shoreline features, and more. Includes index.This layer is part of a selection of digitally scanned and georeferenced historic maps from the Harvard Map Collection. These maps typically portray both natural and manmade features. The selection represents a range of originators, ground condition dates, scales, and map purposes.

Mediterranean Sea Region, West, 1690-1699 (Raster Image)

This layer is a georeferenced raster image of the historic paper map entitled: Accuratissima occidentalioris districtus maris Mediterranei tabula, authore Iusto Danckerts. It was published by Iusto Danckerts, between 1690 and 1699. Scale [ca. 1:5,300,000]. Covers the western Mediterranean Sea region. Map in Latin. The image inside the map neatline is georeferenced to the surface of the earth and fit to the World Miller Cylindrical projected coordinate system. All map collar and inset information is also available as part of the raster image, including any inset maps, profiles, statistical tables, directories, text, illustrations, index maps, legends, or other information associated with the principal map. This map shows features such as drainage, cities and other human settlements, roads, shoreline features, and more. Relief shown pictorially.This layer is part of a selection of digitally scanned and georeferenced historic maps from the Harvard Map Collection. These maps typically portray both natural and manmade features. The selection represents a range of originators, ground condition dates, scales, and map purposes.

Mediterranean Sea Region, East, 1690-1699 (Raster Image)

This layer is a georeferenced raster image of the historic paper map entitled: Accuratissima orientalioris districtus maris Mediterranei tabula, authore Iusto Danckerts. It was published by Iusto Danckerts, between 1690 and 1699. Scale [ca. 1:5,300,000]. Covers the eastern Mediterranean Sea region. Map in Latin. The image inside the map neatline is georeferenced to the surface of the earth and fit to the World Miller Cylindrical projected coordinate system. All map collar and inset information is also available as part of the raster image, including any inset maps, profiles, statistical tables, directories, text, illustrations, index maps, legends, or other information associated with the principal map. This map shows features such as drainage, cities and other human settlements, roads, shoreline features, and more. Relief shown pictorially.This layer is part of a selection of digitally scanned and georeferenced historic maps from the Harvard Map Collection. These maps typically portray both natural and manmade features. The selection represents a range of originators, ground condition dates, scales, and map purposes.

Mediterranean Sea Region, 1690-1699 (Raster Image)

This layer is a georeferenced raster image of the untitled historic paper map: [Mediterranean sea], authore Iusto Danckerts. It was published by Iusto Danckerts, between 1690 and 1699. Scale [ca. 1:5,300,000]. Covers the Mediterranean Sea region. Map in Latin. Source map issued on two sheets, now pasted together. The image inside the map neatline is georeferenced to the surface of the earth and fit to the World Miller Cylindrical projected coordinate system. All map collar and inset information is also available as part of the raster image, including any inset maps, profiles, statistical tables, directories, text, illustrations, index maps, legends, or other information associated with the principal map. This map shows features such as drainage, cities and other human settlements, roads, territorial features, shoreline features, and more. Relief shown pictorially.This layer is part of a selection of digitally scanned and georeferenced historic maps from the Harvard Map Collection. These maps typically portray both natural and manmade features. The selection represents a range of originators, ground condition dates, scales, and map purposes.

Mediterranean Sea Region, 1700-1710 (Raster Image)

This layer is a georeferenced raster image of the historic paper map entitled: Carte nouvelle de la mer Mediterranée : divisée en mer de Levant et de Ponant, subdivisés en leurs principales parties ou mers : avec les observations des Mrs. de l'Académie, dressée par ... Sanson. It was published by chez Pierre Mortier, between 1700 and 1710. Scale [ca. 1:4,375,000]. Covers the Mediterranean Sea and Black Sea regions. Map in French. The image inside the map neatline is georeferenced to the surface of the earth and fit to the World Miller Cylindrical projected coordinate system. All map collar and inset information is also available as part of the raster image, including any inset maps, profiles, statistical tables, directories, text, illustrations, index maps, legends, or other information associated with the principal map. This map shows features such as drainage, cities and other human settlements, roads, territorial boundaries, shoreline features, and more. Relief shown pictorially.This layer is part of a selection of digitally scanned and georeferenced historic maps from the Harvard Map Collection. These maps typically portray both natural and manmade features. The selection represents a range of originators, ground condition dates, scales, and map purposes.

Coasts, Mediterranean Sea Region, 1709 (Raster Image)

This layer is a georeferenced raster image of the historic paper map entitled: Carte de la mer Mediterranée : presentée à Monseigneur le Bailli de Noailles, Lieutenant General des Galeres de France, par Nicolas Therin et H. Michelot, pilotes entretenus sur les Galeres du Roy. It was published by Chez Pierre Mortier, ca. 1709. Scale [ca. 1:3,600,000]. Covers Mediterranean Sea and coasts of Europe, North Africa, and the Middle East. Map in French. The image inside the map neatline is georeferenced to the surface of the earth and fit to the World Miller Cylindrical projected coordinate system. All map collar and inset information is also available as part of the raster image, including any inset maps, profiles, statistical tables, directories, text, illustrations, index maps, legends, or other information associated with the principal map. This map shows features such as drainage, cities and other human settlements, territorial boundaries, shoreline features, and more. Includes index.This layer is part of a selection of digitally scanned and georeferenced historic maps from the Harvard Map Collection. These maps typically portray both natural and manmade features. The selection represents a range of originators, ground condition dates, scales, and map purposes.