Reconstruction of sea-surface conditions at middle to high latitudes of the Northern Hemisphere

A new calibration database of census counts of organic-walled dinoflagellate cyst (dinocyst) assemblages has been developed from the analyses of surface sediment samples collected at middle to high latitudes of the Northern Hemisphere after standardisation of taxonomy and laboratory procedures. The database comprises 940 reference data points from the North Atlantic, Arctic and North Pacific oceans and their adjacent seas, including the Mediterranean Sea, as well as epicontinental environments such as the Estuary and Gulf of St. Lawrence, the Bering Sea and the Hudson Bay. The relative abundance of taxa was analysed to describe the distribution of assemblages. The best analogue technique was used for the reconstruction of Last Glacial Maximum (LGM) sea-surface temperature and salinity during summer and winter, in addition to sea-ice cover extent, at sites from the North Atlantic (n=63), Mediterranean Sea (n=1) and eastern North Pacific (n=1). Three of the North Atlantic cores, from the continental margin of eastern Canada, revealed a barren LGM interval, probably because of quasi-permanent sea ice. Six other cores from the Greenland and Norwegian seas were excluded from the compilation because of too sparse assemblages and poor analogue situation. At the remaining sites (n= 54), relatively close modern analogues were found for most LGM samples, which allowed reconstructions. The new LGM results are consistent with previous reconstructions based on dinocyst data, which show much cooler conditions than at present along the continental margins of Canada and Europe, but sharp gradients of increasing temperature offshore. The results also suggest low salinity and larger than present contrasts in seasonal temperatures with colder winters and more extensive sea-ice cover, whereas relatively warm conditions may have prevailed offshore in summer. From these data, we hypothesise low thermal inertia in a shallow and low-density surface water layer.

Photochemical responses of the diatom Skeletonema costatum grown under elevated CO2 concentrations to short-term changes in pH

Variability in pH is a common occurrence in many aquatic environments, due to physical, chemical and biological processes. In coastal waters, lagoons, estuaries and inland waters, pH can change very rapidly (within seconds or hours) in addition to daily and seasonal changes. At the same time, progressive ocean acidification caused by anthropogenic CO2 emissions is superimposed on these spatial and temporal pH changes. Photosynthetic organisms are therefore unavoidably subject to significant pH variations at the cell surface. Whether this will affect their response to long-term ocean acidification is still unknown, nor is it known whether the short-term sensitivity to pH change is affected by the pCO2 to which the cells are acclimated. We posed the latter open question as our experimental hypothesis: Does acclimation to seawater acidification affect the response of phytoplankton to acute pH variations? The diatom Skeletonema costatum, commonly found in coastal and estuarine waters where short-term acute changes in pH frequently occur, was selected to test the hypothesis. Diatoms were grown at both 390 (pH 8.2, low CO2; LC) and 1000 (pH 7.9, high CO2; HC) µatm CO2 for at least 20 generations, and photosynthetic responses to short-term and acute changes in pH (between 8.2 and 7.6) were investigated. The effective quantum yield of LC-grown cells decreased by ca. 70% only when exposed to pH 7.6; this was not observed when exposed to pH 7.9 or 8.2. HC-grown cells did not show significant responses in any pH treatment. Non-photochemical quenching showed opposite trends. In general, our results indicate that while LC-grown cells are rather sensitive to acidification, HC-grown cells are relatively unresponsive in terms of photochemical performance.

Underway measurements of halocarbons (air) during POSEIDON cruise POS399 in June 2010

During the DRIVE (Diurnal and Regional Variability of Halogen Emissions) ship campaign we investigated the variability of the halogenated very short-lived substances (VSLS) bromoform (CHBr3), dibromomethane (CH2Br2) and methyl iodide (CH3I) in the marine atmospheric boundary layer in the eastern tropical and subtropical North Atlantic Ocean during May/June 2010. The highest VSLS mixing ratios were found near the Mauritanian coast and close to Lisbon (Portugal). With backward trajectories we identified predominantly air masses from the open North Atlantic with some coastal influence in the Mauritanian upwelling area, due to the prevailing NW winds. The maximum VSLS mixing ratios above the Mauritanian upwelling were 8.92 ppt for bromoform, 3.14 ppt for dibromomethane and 3.29 ppt for methyl iodide, with an observed maximum range of the daily mean up to 50% for bromoform, 26% for dibromomethane and 56% for methyl iodide. The influence of various meteorological parameters - such as wind, surface air pressure, surface air and surface water temperature, humidity and marine atmospheric boundary layer (MABL) height - on VSLS concentrations and fluxes was investigated. The strongest relationship was found between the MABL height and bromoform, dibromomethane and methyl iodide abundances. Lowest MABL heights above the Mauritanian upwelling area coincide with highest VSLS mixing ratios and vice versa above the open ocean. Significant high anti-correlations confirm this relationship for the whole cruise. We conclude that especially above oceanic upwelling systems, in addition to sea-air fluxes, MABL height variations can influence atmospheric VSLS mixing ratios, occasionally leading to elevated atmospheric abundances. This may add to the postulated missing VSLS sources in the Mauritanian upwelling region (Quack et al., 2007).

The photo-physiological response of a model cnidarian-dinoflagellate symbiosis to CO2-induced acidification at the cellular level

We measured the relationship between CO2-induced seawater acidification, photo-physiological performance and intracellular pH (pHi) in a model cnidarian-dinoflagellate symbiosis - the sea anemone Aiptasia sp. -under ambient (289.94 ± 12.54 µatm), intermediate (687.40 ± 25.10 µatm) and high (1459.92 ± 65.51 µatm) CO2 conditions. These treatments represented current CO2 levels, in addition to CO2 stabilisation scenarios IV and VI provided by the Intergovernmental Panel on Climate Change (IPCC). Anemones were exposed to each treatment for two months and sampled at regular intervals. At each time-point we measured a series of physiological responses: maximum dark-adapted fluorescent yield of PSII (Fv/Fm), gross photosynthetic rate, respiration rate, symbiont population density, and light-adapted pHi of both the dinoflagellate symbiont and isolated host anemone cell. We observed increases in all but one photo-physiological parameter (Pgross:R ratio). At the cellular level, increases in light-adapted symbiont pHi were observed under both intermediate and high CO2 treatments, relative to control conditions (pHi 7.35 and 7.46 versus pHi 7.25, respectively). The response of light-adapted host pHi was more complex, however, with no change observed under the intermediate CO2 treatment, but a 0.3 pH-unit increase under the high CO2 treatment (pHi 7.19 and 7.48, respectively). This difference is likely a result of a disproportionate increase in photosynthesis relative to respiration at the higher CO2 concentration. Our results suggest that, rather than causing cellular acidosis, the addition of CO2 will enhance photosynthetic performance, enabling both the symbiont and host cell to withstand predicted ocean acidification scenarios.

Parental exposure to elevated pCO2 influences the reproductive success of copepods

Substantial variations are reported for egg production and hatching rates of copepods exposed to elevated carbon dioxide concentrations (pCO2). One possible explanation, as found in other marine taxa, is that prior parental exposure to elevated pCO2 (and/or decreased pH) affects reproductive performance. Previous studies have adopted two distinct approaches, either (1) expose male and female copepoda to the test pCO2/pH scenarios, or (2) solely expose egg-laying females to the tests. Although the former approach is more realistic, the majority of studies have used the latter approach. Here, we investigated the variation in egg production and hatching success of Acartia tonsa between these two experimental designs, across five different pCO2 concentrations (385-6000 µatm pCO2). In addition, to determine the effect of pCO2 on the hatching success with no prior parental exposure, eggs produced and fertilized under ambient conditions were also exposed to these pCO2 scenarios. Significant variations were found between experimental designs, with approach (1) resulting in higher impacts; here >20% difference was seen in hatching success between experiments at 1000 µatm pCO2 scenarios (2100 year scenario), and >85% at 6000 µatm pCO2. This study highlights the potential to misrepresent the reproductive response of a species to elevated pCO2 dependent on parental exposure.

(Table 1) Porewater analysis with respect to gypsum from different ODP holes

The stability of gypsum in marine sediments has been investigated through the calculation of its saturation index at the sediment in situ temperature and pressure, using the entire ODP/IODP porewater composition database (14416 samples recovered from sediments collected during 95 ODP and IODP Legs). Saturation is reached in sediment porewaters of 26 boreholes drilled at 23 different sites, during 12 ODP/IODP Legs. As ocean bottom seawater is largely undersaturated with respect to gypsum, the porewater Ca content or its SO4 concentration, or both, must increase in order to reach equilibrium. At several sites equilibrium is reached either through the presence of evaporitic gypsum layers found in the sedimentary sequence, and/or through a salinity increase due to the presence of evaporitic brines with high concentrations of Ca and SO4. Saturation can also be reached in porewaters of seawater-like salinity (~ 35 per mil), provided sulfate reduction is limited. In this case, saturation is due to the alteration of volcanogenic material which releases large amounts of Ca to the porewaters, where the Ca concentration can reach 55 times its seawater value as for example at ODP Leg 134 site 833. At a few sites, saturation is reached in hydrothermal environments, or as a consequence of the alteration of the basaltic basement. In addition to the well known influence of brines on the formation of gypsum, these results indicate that the alteration of sediments rich in volcanogenic material is a major process leading to gypsum saturation in marine sediment porewaters. Therefore, the presence of gypsum in ancient and recent marine sediments should not be systematically interpreted as due to hypersaline waters, especially if volcanogenic material is present.

Acoustic records of the underwater soundscape at PALAOA with links to audio stream files, 2005-2011

Scientific background: Marine mammals use sound for communication, navigation and prey detection. Acoustic sensors therefore allow the detection of marine mammals, even during polar winter months, when restricted visibility prohibits visual sightings. The animals are surrounded by a permanent natural soundscape, which, in polar waters, is mainly dominated by the movement of ice. In addition to the detection of marine mammals, acoustic long-term recordings provide information on intensity and temporal variability of characteristic natural and anthropogenic background sounds, as well as their influence on the vocalization of marine mammals Scientific objectives: The PerenniAL Acoustic Observatory in the Antarctic Ocean (PALAOA, Hawaiian "whale") near Neumayer Station is intended to record the underwater soundscape in the vicinity of the shelf ice edge over the duration of several years. These long-term recordings will allow studying the acoustic repertoire of whales and seals continuously in an environment almost undisturbed by humans. The data will be analyzed to (1) register species specific vocalizations, (2) infer the approximate number of animals inside the measuring range, (3) calculate their movements relative to the observatory, and (4) examine possible effects of the sporadic shipping traffic on the acoustic and locomotive behaviour of marine mammals. The data, which are largely free of anthropogenic noise, provide also a base to set up passive acoustic mitigation systems used on research vessels. Noise-free bioacoustic data thereby represent the foundation for the development of automatic pattern recognition procedures in the presence of interfering sounds, e.g. propeller noise.

Combined use of gliadins and SSRs to analyse the genetic variability of the Spanish collection of cultivated diploid wheat (Triticum monococcum L. ssp. monococcum)

This work studied the combined use of gliadins and SSRs to analyse inter- and intra-accession variability of the Spanish collection of cultivated einkorn (Triticum monococcum L. ssp. monococcum) maintained at the CRF-INIA. In general, gliadin loci presented higher discrimination power than SSRs, reflecting the high variability of the gliadins. The loci on chromosome 6A were the most polymorphic with similar PIC values for both marker systems, showing that these markers are very useful for genetic variability studies in wheat. The gliadin results indicated that the Spanish einkorn collection possessed high genetic diversity, being the differentiation large between varieties and small within them. Some associations between gliadin alleles and geographical and agro-morphological data were found. Agro-morphological relations were also observed in the clusters of the SSRs dendrogram. A high concordance was found between gliadins and SSRs for genotype identification. In addition, both systems provide complementary information to resolve the different cases of intra-accession variability not detected at the agro-morphological level, and to identify separately all the genotypes analysed. The combined use of both genetic markers is an excellent tool for genetic resource evaluation in addition to agro-morphological evaluation.

An automatic translation scheme from prolog to the andorra kernel language

The Andorra family of languages (which includes the Andorra Kernel Language -AKL) is aimed, in principie, at simultaneously supporting the programming styles of Prolog and committed choice languages. On the other hand, AKL requires a somewhat detailed specification of control by the user. This could be avoided by programming in Prolog to run on AKL. However, Prolog programs cannot be executed directly on AKL. This is due to a number of factors, from more or less trivial syntactic differences to more involved issues such as the treatment of cut and making the exploitation of certain types of parallelism possible. This paper provides basic guidelines for constructing an automatic compiler of Prolog programs into AKL, which can bridge those differences. In addition to supporting Prolog, our style of translation achieves independent and-parallel execution where possible, which is relevant since this type of parallel execution preserves, through the translation, the user-perceived "complexity" of the original Prolog program.

An automatic translation scheme from CLP to AKL

The Andorra Kernel language scheme was aimed, in principle, at simultaneously supporting the programming styles of Prolog and committed choice languages. Within the constraint programming paradigm, this family of languages could also in principle support the concurrent constraint paradigm. This happens for the Agents Kernel Language (AKL). On the other hand, AKL requires a somewhat detailed specification of control by the user. This could be avoided by programming in CLP to run on AKL. However, CLP programs cannot be executed directly on AKL. This is due to a number of factors, from more or less trivial syntactic differences to more involved issues such as the treatment of cut and making the exploitation of certain types of parallelism possible. This paper provides a translation scheme which is a basis of an automatic compiler of CLP programs into AKL, which can bridge those differences. In addition to supporting CLP, our style of translation achieves independent and-parallel execution where possible, which is relevant since this type of parallel execution preserves, through the translation, the user-perceived "complexity" of the original program.

Nonlinear response to ultrasound of encapsulated microbubbles

The acoustic backscatter of encapsulated gas-filled microbubbles immersed in a weak compressible liquid and irradiated by ultrasound fields of moderate to high pressure amplitudes is investigated theoretically. The problem is formulated by considering, for the viscoelastic shell of finite thickness, an isotropic hyperelastic neo-Hookean model for the elastic contribution in addition to a Newtonian viscous component. First and second harmonic scattering cross-sections have been evaluated and the quantitative influence of the driving pressure amplitude on the harmonic resonance frequencies for different initial equilibrium bubble sizes and for different encapsulating physical properties has been determined. Conditions for optimal second harmonic imaging have been also investigated and some regions in the parameters space where the second harmonic intensity is dominant over the fundamental have been identified. Results have been obtained for albumin, lipid and polymer encapsulating shells, respectively.

Virtual Excursions: a new way to explore science in class

The educational platform Virtual Science Hub (ViSH) has been developed as part of the GLOBAL excursion European project. ViSH (http://vishub.org/) is a portal where teachers and scientist interact to create virtual excursions to science infrastructures. The main motivation behind the project was to connect teachers - and in consequence their students - to scientific institutions and their wide amount of infrastructures and resources they are working with. Thus the idea of a hub was born that would allow the two worlds of scientists and teachers to connect and to innovate science teaching. The core of the ViSH?s concept design is based on virtual excursions, which allow for a number of pedagogical models to be applied. According to our internal definition a virtual excursion is a tour through some digital context by teachers and pupils on a given topic that is attractive and has an educational purpose. Inquiry-based learning, project-based and problem-based learning are the most prominent approaches that a virtual excursion may serve. The domain specific resources and scientific infrastructures currently available on the ViSH are focusing on life sciences, nano-technology, biotechnology, grid and volunteer computing. The virtual excursion approach allows an easy combination of these resources into interdisciplinary teaching scenarios. In addition, social networking features support the users in collaborating and communicating in relation to these excursions and thus create a community of interest for innovative science teaching. The design and development phases were performed following a participatory design approach. An important aspect in this process was to create design partnerships amongst all actors involved, researchers, developers, infrastructure providers, teachers, social scientists, and pedagogical experts early in the project. A joint sense of ownership was created and important changes during the conceptual phase were implemented in the ViSH due to early user feedback. Technology-wise the ViSH is based on the latest web technologies in order to make it cross-platform compatible so that it works on several operative systems such as Windows, Mac or Linux and multi-device accessible, such as desktop, tablet and mobile devices. The platform has been developed in HTML5, the latest standard for web development, assuring that it can run on any modern browser. In addition to social networking features a core element on the ViSH is the virtual excursions editor. It is a web tool that allows teachers and scientists to create rich mash-ups of learning resources provided by the e-Infrastructures (i.e. remote laboratories and live webcams). These rich mash-ups can be presented in either slides or flashcards format. Taking advantage of the web architecture supported, additional powerful components have been integrated like a recommendation engine to provide personalized suggestions about educational content or interesting users and a videoconference tool to enhance real-time collaboration like MashMeTV (http://www.mashme.tv/).

VOTESCRIPT: Telematic voting system designed to enable final count verification

In this paper we present a global description of a telematic voting system based on advanced cryptography and on the use of smart cards (VOTESCRIPT system) whose most outstanding characteristic is the ability to verify that the tally carried out by the system is correct, meaning that the results published by the system correspond with votes cast. The VOTESCRIPT system provides an individual verification mechanism allowing each Voter to confirm whether his vote has been correctly counted. The innovation with respect to other solutions lies in the fact that the verification process is private so that Voters have no way of proving what they voted in the presence of a non-authorized third party. Vote buying and selling or any other kind of extortion are prevented. The existence of the Intervention Systems allows the whole electoral process to be controlled by groups of citizens or authorized candidatures. In addition to this the system can simply make an audit not only of the final results, but also of the whole process. Global verification provides the Scrutineers with robust cryptographic evidence which enables unequivocal proof if the system has operated in a fraudulent way.

The NACIR Project: New Applications for CPV'S: A fast way to improve reliability and Technology progress

The main goal of this proposal is to join together the owners of the most advanced CPV technology, with respect to the state of the art, in order to research from its leading position new applications for CPV systems. In addition to opening up new markets, it will unveil possible sources of failure in new environments outside Europe, in order to assure component reliability. The proposed project will also try to improve the current technology of the industrial partners (ISOFOTON and CONCENTRIX) by accelerating the learning curve that CPV must follow in order to reach the competitive market, and lowering the cost under the current flat panel PV significantly within 3-4 years. The use of CPV systems in remote areas, together with harsher radiation, ambient and infrastructure conditions will help to increase the rate of progress of this technology. In addition, the ISFOC s contribution, which brings together seven power plants from seven CPV technologies up to 3 MWpeak, will allow creating the most complete database of components and systems performance to be generated as well as the effects of radiation and meteorology on systems operations. Finally, regarding the new applications for CPV subject, the project will use a CPV system sized 25 kWp in a stand-alone station in Egypt (NWRC) for the first time for water pumping and irrigation purposes. In a similar way ISOFOTON will connect up to 25 kWp CPV to the Moroccan ONE utility grid. From the research content point of view of this project, which is directly addressed by the scope of the call, the cooperative research between UPM, FhG-ISE and the two companies will be favoured by the fact that all are progressing in similar directions: developing two-stage optics CPV systems. In addition to these technology improvements the UPM is very interested in developing a new concept of module, recently patented, which will fulfil all required characteristics of a good CPV with less components and reducing cost.

Estimate of the total ice volume of Svalbard glaciers and their potential contribution to sea-level rise = Estimación del volumen de hielo de Svalbard y su contribución potencial a la subida del nivel del mar

El objetivo final de las investigaciones recogidas en esta tesis doctoral es la estimación del volumen de hielo total de los ms de 1600 glaciares de Svalbard, en el Ártico, y, con ello, su contribución potencial a la subida del nivel medio del mar en un escenario de calentamiento global. Los cálculos más exactos del volumen de un glaciar se efectúan a partir de medidas del espesor de hielo obtenidas con georradar. Sin embargo, estas medidas no son viables para conjuntos grandes de glaciares, debido al coste, dificultades logísticas y tiempo requerido por ellas, especialmente en las regiones polares o de montaña. Frente a ello, la determinación de áreas de glaciares a partir de imágenes de satélite sí es viable a escalas global y regional, por lo que las relaciones de escala volumen-área constituyen el mecanismo más adecuado para las estimaciones de volúmenes globales y regionales, como las realizadas para Svalbard en esta tesis. Como parte del trabajo de tesis, hemos elaborado un inventario de los glaciares de Svalbard en los que se han efectuado radioecosondeos, y hemos realizado los cálculos del volumen de hielo de más de 80 cuencas glaciares de Svalbard a partir de datos de georradar. Estos volúmenes han sido utilizados para calibrar las relaciones volumen-área desarrolladas en la tesis. Los datos de georradar han sido obtenidos en diversas campañas llevadas a cabo por grupos de investigación internacionales, gran parte de ellas lideradas por el Grupo de Simulación Numérica en Ciencias e Ingeniería de la Universidad Politécnica de Madrid, del que forman parte la doctoranda y los directores de tesis. Además, se ha desarrollado una metodología para la estimación del error en el cálculo de volumen, que aporta una novedosa técnica de cálculo del error de interpolación para conjuntos de datos del tipo de los obtenidos con perfiles de georradar, que presentan distribuciones espaciales con unos patrones muy característicos pero con una densidad de datos muy irregular. Hemos obtenido en este trabajo de tesis relaciones de escala específicas para los glaciares de Svalbard, explorando la sensibilidad de los parámetros a diferentes morfologías glaciares, e incorporando nuevas variables. En particular, hemos efectuado experimentos orientados a verificar si las relaciones de escala obtenidas caracterizando los glaciares individuales por su tamaño, pendiente o forma implican diferencias significativas en el volumen total estimado para los glaciares de Svalbard, y si esta partición implica algún patrón significativo en los parámetros de las relaciones de escala. Nuestros resultados indican que, para un valor constante del factor multiplicativo de la relacin de escala, el exponente que afecta al área en la relación volumen-área decrece según aumentan la pendiente y el factor de forma, mientras que las clasificaciones basadas en tamaño no muestran un patrón significativo. Esto significa que los glaciares con mayores pendientes y de tipo circo son menos sensibles a los cambios de área. Además, los volúmenes de la población total de los glaciares de Svalbard calculados con fraccionamiento en grupos por tamaño y pendiente son un 1-4% menores que los obtenidas usando la totalidad de glaciares sin fraccionamiento en grupos, mientras que los volúmenes calculados fraccionando por forma son un 3-5% mayores. También realizamos experimentos multivariable para obtener estimaciones óptimas del volumen total mediante una combinación de distintos predictores. Nuestros resultados muestran que un modelo potencial simple volumen-área explica el 98.6% de la varianza. Sólo el predictor longitud del glaciar proporciona significación estadística cuando se usa además del área del glaciar, aunque el coeficiente de determinación disminuye en comparación con el modelo más simple V-A. El predictor intervalo de altitud no proporciona información adicional cuando se usa además del área del glaciar. Nuestras estimaciones del volumen de la totalidad de glaciares de Svalbard usando las diferentes relaciones de escala obtenidas en esta tesis oscilan entre 6890 y 8106 km3, con errores relativos del orden de 6.6-8.1%. El valor medio de nuestras estimaciones, que puede ser considerado como nuestra mejor estimación del volumen, es de 7.504 km3. En términos de equivalente en nivel del mar (SLE), nuestras estimaciones corresponden a una subida potencial del nivel del mar de 17-20 mm SLE, promediando 19_2 mm SLE, donde el error corresponde al error en volumen antes indicado. En comparación, las estimaciones usando las relaciones V-A de otros autores son de 13-26 mm SLE, promediando 20 _ 2 mm SLE, donde el error representa la desviación estándar de las distintas estimaciones. ABSTRACT The final aim of the research involved in this doctoral thesis is the estimation of the total ice volume of the more than 1600 glaciers of Svalbard, in the Arctic region, and thus their potential contribution to sea-level rise under a global warming scenario. The most accurate calculations of glacier volumes are those based on ice-thicknesses measured by groundpenetrating radar (GPR). However, such measurements are not viable for very large sets of glaciers, due to their cost, logistic difficulties and time requirements, especially in polar or mountain regions. On the contrary, the calculation of glacier areas from satellite images is perfectly viable at global and regional scales, so the volume-area scaling relationships are the most useful tool to determine glacier volumes at global and regional scales, as done for Svalbard in this PhD thesis. As part of the PhD work, we have compiled an inventory of the radio-echo sounded glaciers in Svalbard, and we have performed the volume calculations for more than 80 glacier basins in Svalbard from GPR data. These volumes have been used to calibrate the volume-area relationships derived in this dissertation. Such GPR data have been obtained during fieldwork campaigns carried out by international teams, often lead by the Group of Numerical Simulation in Science and Engineering of the Technical University of Madrid, to which the PhD candidate and her supervisors belong. Furthermore, we have developed a methodology to estimate the error in the volume calculation, which includes a novel technique to calculate the interpolation error for data sets of the type produced by GPR profiling, which show very characteristic data distribution patterns but with very irregular data density. We have derived in this dissertation scaling relationships specific for Svalbard glaciers, exploring the sensitivity of the scaling parameters to different glacier morphologies and adding new variables. In particular, we did experiments aimed to verify whether scaling relationships obtained through characterization of individual glacier shape, slope and size imply significant differences in the estimated volume of the total population of Svalbard glaciers, and whether this partitioning implies any noticeable pattern in the scaling relationship parameters. Our results indicate that, for a fixed value of the factor in the scaling relationship, the exponent of the area in the volume-area relationship decreases as slope and shape increase, whereas size-based classifications do not reveal any clear trend. This means that steep slopes and cirque-type glaciers are less sensitive to changes in glacier area. Moreover, the volumes of the total population of Svalbard glaciers calculated according to partitioning in subgroups by size and slope are smaller (by 1-4%) than that obtained considering all glaciers without partitioning into subgroups, whereas the volumes calculated according to partitioning in subgroups by shape are 3-5% larger. We also did multivariate experiments attempting to optimally predict the volume of Svalbard glaciers from a combination of different predictors. Our results show that a simple power-type V-A model explains 98.6% of the variance. Only the predictor glacier length provides statistical significance when used in addition to the predictor glacier area, though the coefficient of determination decreases as compared with the simpler V-A model. The predictor elevation range did not provide any additional information when used in addition to glacier area. Our estimates of the volume of the entire population of Svalbard glaciers using the different scaling relationships that we have derived along this thesis range within 6890-8106 km3, with estimated relative errors in total volume of the order of 6.6-8.1% The average value of all of our estimates, which could be used as a best estimate for the volume, is 7,504 km3. In terms of sea-level equivalent (SLE), our volume estimates correspond to a potential contribution to sea-level rise within 17-20 mm SLE, averaging 19 _ 2 mm SLE, where the quoted error corresponds to our estimated relative error in volume. For comparison, the estimates using the V-A scaling relations found in the literature range within 13-26 mm SLE, averaging 20 _ 2 mm SLE, where the quoted error represents the standard deviation of the different estimates.