Spatially explicit estimates of length and biomass of Clupea harengus (Norwegian spring spawning herring) from acoustic and trawl surveys in the North-East Atlantic in May 2010

Acoustic estimates of herring and blue whiting abundance were obtained during the surveys using the Simrad ER60 scientific echosounder. The allocation of NASC-values to herring, blue whiting and other acoustic targets were based on the composition of the trawl catches and the appearance of echo recordings. To estimate the abundance, the allocated NASC -values were averaged for ICES-squares (0.5° latitude by 1° longitude). For each statistical square, the unit area density of fish (rA) in number per square nautical mile (N*nm-2) was calculated using standard equations (Foote et al., 1987; Toresen et al., 1998). To estimate the total abundance of fish, the unit area abundance for each statistical square was multiplied by the number of square nautical miles in each statistical square and then summed for all the statistical squares within defined subareas and over the total area. Biomass estimation was calculated by multiplying abundance in numbers by the average weight of the fish in each statistical square then summing all squares within defined subareas and over the total area. The Norwegian BEAM soft-ware (Totland and Godø 2001) was used to make estimates of total biomass.

Spatially explicit estimates of length and biomass of Clupea harengus (Norwegian spring spawning herring) from acoustic and trawl surveys in the North-East Atlantic in May 2012

Acoustic estimates of herring and blue whiting abundance were obtained during the surveys using the Simrad ER60 scientific echosounder. The allocation of NASC-values to herring, blue whiting and other acoustic targets were based on the composition of the trawl catches and the appearance of echo recordings. To estimate the abundance, the allocated NASC -values were averaged for ICES-squares (0.5° latitude by 1° longitude). For each statistical square, the unit area density of fish (rA) in number per square nautical mile (N*nm-2) was calculated using standard equations (Foote et al., 1987; Toresen et al., 1998). To estimate the total abundance of fish, the unit area abundance for each statistical square was multiplied by the number of square nautical miles in each statistical square and then summed for all the statistical squares within defined subareas and over the total area. Biomass estimation was calculated by multiplying abundance in numbers by the average weight of the fish in each statistical square then summing all squares within defined subareas and over the total area. The Norwegian BEAM soft-ware (Totland and Godø 2001) was used to make estimates of total biomass.

Spatially explicit estimates of length and biomass of Micromesistius poutassou (Blue Whiting) from acoustic and trawl surveys in the North-East Atlantic in May 2004

Acoustic estimates of herring and blue whiting abundance were obtained during the surveys using the Simrad ER60 scientific echosounder. The allocation of NASC-values to herring, blue whiting and other acoustic targets were based on the composition of the trawl catches and the appearance of echo recordings. To estimate the abundance, the allocated NASC -values were averaged for ICES-squares (0.5° latitude by 1° longitude). For each statistical square, the unit area density of fish (rA) in number per square nautical mile (N*nm-2) was calculated using standard equations (Foote et al., 1987; Toresen et al., 1998). To estimate the total abundance of fish, the unit area abundance for each statistical square was multiplied by the number of square nautical miles in each statistical square and then summed for all the statistical squares within defined subareas and over the total area. Biomass estimation was calculated by multiplying abundance in numbers by the average weight of the fish in each statistical square then summing all squares within defined subareas and over the total area. The Norwegian BEAM soft-ware (Totland and Godø 2001) was used to make estimates of total biomass.

(Table 1) Circulating levels of TT3 and TT4, and age, length, girth, and estimated body mass in polar bears (Ursus maritimus) from East Greenland 1999-2001

We investigated the multivariate relationships between adipose tissue residue levels of 48 individual organohalogen contaminants (OHCs) and circulating thyroid hormone (TH) levels in polar bears (Ursus maritimus) from East Greenland (1999-2001, n = 62), using projection to latent structure (PLS) regression for four groupings of polar bears; subadults (SubA), adult females with cubs (AdF_N), adult females without cubs (AdF_S) and adult males (AdM). In the resulting significant PLS models for SubA, AdF_N and AdF_S, some OHCs were especially important in explaining variations in circulating TH levels: polybrominated diphenylether (PBDE)-99, PBDE-100, PBDE-153, polychlorinated biphenyl (PCB)-52, PCB-118, cis-nonachlor, trans-nonachlor, trichlorobenzene (TCB) and pentachlorobenzene (QCB), and both negative and positive relationships with THs were found. In addition, the models revealed that DDTs had a positive influence on total 3,5,3'-triiodothyronine (TT3) in AdF_S, and that a group of 17 higher chlorinated ortho-PCBs had a positive influence on total 3,5,3',5'-tetraiodothyronine (thyroxine, TT4) in AdF_N. TH levels in AdM seemed less influenced by OHCs because of non-significant PLS models. TH levels were also influenced by biological factors such as age, sex, body size, lipid content of adipose tissue and sampling date. When controlling for biological variables, the major relationships from the PLS models for SubA, AdF_N and AdF_S were found significant in partial correlations. The most important OHCs that influenced TH levels in the significant PLS models may potentially act through similar mechanisms on the hypothalamic-pituitary-thyroid (HPT) axis, suggesting that both combined effects by dose and response addition and perhaps synergistic potentiation may be a possibility in these polar bears. Statistical associations are not evidence per se of biological cause-effect relationships. Still, the results of the present study indicate that OHCs may affect circulating TH levels in East Greenland polar bears, adding to the "weight of evidence" suggesting that OHCs might interfere with thyroid homeostasis in polar bears.

Seawater carbonate chemistry and shell length of northern abalone (Haliotis kamtschatkana) during experiments, 2011

Increasing levels of anthropogenic carbon dioxide in the world's oceans are resulting in a decrease in the availability of carbonate ions and a drop in seawater pH. This process, known as ocean acidification, is a potential threat to marine populations via alterations in survival and development. To date, however, little research has examined the effects of ocean acidification on rare or endangered species. To begin to assess the impacts of acidification on endangered northern abalone (Haliotis kamtschatkana) populations, we exposed H. kamtschatkana larvae to various levels of CO2 [400 ppm (ambient), 800 ppm, and 1800 ppm CO2] and measured survival, settlement, shell size, and shell development. Larval survival decreased by ca. 40% in elevated CO2 treatments relative to the 400 ppm control. However, CO2 had no effect on the proportion of surviving larvae that metamorphosed at the end of the experiment. Larval shell abnormalities became apparent in approximately 40% of larvae reared at 800 ppm CO2, and almost all larvae reared at 1800 ppm CO2 either developed an abnormal shell or lacked a shell completely. Of the larvae that did not show shell abnormalities, shell size was reduced by 5% at 800 ppm compared to the control. Overall, larval development of H. kamtschatkana was found to be sensitive to ocean acidification. Near future levels of CO2 will likely pose a significant additional threat to this species, which is already endangered with extinction due in part to limited reproductive output and larval recruitment.

(Table 1) Age, fork length, tissue d15N and d13C and mercury content of arctic char (Salvelinus alpinus) across the Canadian Arctic

Among-lake variation in mercury (Hg) concentrations in landlocked Arctic char was examined in 27 char populations from remote lakes across the Canadian Arctic. A total of 520 landlocked Arctic char were collected from 27 lakes, as well as sediments and surface water from a subset of lakes in 1999, 2002, and 2005 to 2007. Size, length, age, and trophic position (d15N) of individual char were determined and relationships with total Hg (THg) concentrations investigated, to identify a common covariate for adjustment using analysis of covariance (ANCOVA). A subset of 216 char from 24 populations was used for spatial comparison, after length-adjustment. The influence of trophic position and food web length and abiotic characteristics such as location, geomorphology, lake area, catchment area, catchment-to-lake area ratio of the lakes on adjusted THg concentrations in char muscle tissue were then evaluated. Arctic char from Amituk Lake (Cornwallis Island) had the highest Hg concentrations (1.31 µg/g wet wt), while Tessisoak Lake (Labrador, 0.07 µg/g wet wt) had the lowest. Concentrations of THg were positively correlated with size, d15N, and age, respectively, in 88,71, and 58% of 24 char populations. Length and d15N were correlated in 67% of 24 char populations. Food chain length did not explain the differences in length-adjusted THg concentrations in char. No relationships between adjusted THg concentrations in char and latitude or longitude were found, however, THg concentrations in char showed a positive correlation with catchment-to-lake area ratio. Furthermore, we conclude that inputs from the surrounding environment may influence THg concentrations, and will ultimately affect THg concentrations in char as a result of predicted climate-driven changes that may occur in Arctic lake watersheds.

Seawater carbonate chemistry, length and survival of Inland silverside, Menidia beryllina, during experiments, 2012

Absorption of anthropogenic carbon dioxide by the world's oceans is causing mankind's 'other CO2 problem', ocean acidification. Although this process will challenge marine organisms that synthesize calcareous exoskeletons or shells, it is unclear how it will affect internally calcifying organisms, such as marine fish. Adult fish tolerate short-term exposures to CO2 levels that exceed those predicted for the next 300 years (~2,000 ppm), but potential effects of increased CO2 on growth and survival during the early life stages of fish remain poorly understood. Here we show that the exposure of early life stages of a common estuarine fish (Menidia beryllina) to CO2 concentrations expected in the world's oceans later this century caused severely reduced survival and growth rates. When compared with present-day CO2 levels (~400 ppm), exposure of M. beryllina embryos to ~1,000 ppm until one week post-hatch reduced average survival and length by 74% and 18%, respectively. The egg stage was significantly more vulnerable to high CO2-induced mortality than the post-hatch larval stage. These findings challenge the belief that ocean acidification will not affect fish populations, because even small changes in early life survival can generate large fluctuations in adult-fish abundance.

Seawater carbonate chemistry and length, weight, survival rate, metabolic rate of coral reef fish Amphiprion melanopus in a laboratory experiment

Carbon dioxide concentrations in the surface ocean are increasing owing to rising CO2 concentrations in the atmosphere. Higher CO2 levels are predicted to affect essential physiological processes of many aquatic organisms, leading to widespread impacts on marine diversity and ecosystem function, especially when combined with the effects of global warming. Yet the ability for marine species to adjust to increasing CO2 levels over many generations is an unresolved issue. Here we show that ocean conditions projected for the end of the century (approximately 1,000 µatm CO2 and a temperature rise of 1.5-3.0 °C) cause an increase in metabolic rate and decreases in length, weight, condition and survival of juvenile fish. However, these effects are absent or reversed when parents also experience high CO2 concentrations. Our results show that non-genetic parental effects can dramatically alter the response of marine organisms to increasing CO2 and demonstrate that some species have more capacity to acclimate to ocean acidification than previously thought.

Freshwater mesocosm experiment 2012 in Vancouver, Canada on effects of increased temperature and food chain length on oxygen fluxes and plankton community structure

Over broad thermal gradients, the effect of temperature on aerobic respiration and photosynthesis rates explains variation in community structure and function. Yet for local communities, temperature dependent trophic interactions may dominate effects of warming. We tested the hypothesis that food chain length modifies the temperature-dependence of ecosystem fluxes and community structure. In a multi-generation aquatic food web experiment, increasing temperature strengthened a trophic cascade, altering the effect of temperature on estimated mass-corrected ecosystem fluxes. Compared to consumer-free and 3-level food chains, grazer-algae (2-level) food chains responded most strongly to the temperature gradient. Temperature altered community structure, shifting species composition and reducing zooplankton density and body size. Still, food chain length did not alter the temperature dependence of net ecosystem fluxes. We conclude that locally, food chain length interacts with temperature to modify community structure, but only temperature, not food chain length influenced net ecosystem fluxes.

Spatially explicit estimates of length and biomass of Micromesistius poutassou (Blue Whiting) from acoustic and trawl surveys in the North-East Atlantic in May 2010

Acoustic estimates of herring and blue whiting abundance were obtained during the surveys using the Simrad ER60 scientific echosounder. The allocation of NASC-values to herring, blue whiting and other acoustic targets were based on the composition of the trawl catches and the appearance of echo recordings. To estimate the abundance, the allocated NASC -values were averaged for ICES-squares (0.5° latitude by 1° longitude). For each statistical square, the unit area density of fish (rA) in number per square nautical mile (N*nm-2) was calculated using standard equations (Foote et al., 1987; Toresen et al., 1998). To estimate the total abundance of fish, the unit area abundance for each statistical square was multiplied by the number of square nautical miles in each statistical square and then summed for all the statistical squares within defined subareas and over the total area. Biomass estimation was calculated by multiplying abundance in numbers by the average weight of the fish in each statistical square then summing all squares within defined subareas and over the total area. The Norwegian BEAM soft-ware (Totland and Godø 2001) was used to make estimates of total biomass.

Spatially explicit estimates of length and biomass of Clupea harengus (Norwegian spring spawning herring) from acoustic and trawl surveys in the North-East Atlantic in May 2011

Acoustic estimates of herring and blue whiting abundance were obtained during the surveys using the Simrad ER60 scientific echosounder. The allocation of NASC-values to herring, blue whiting and other acoustic targets were based on the composition of the trawl catches and the appearance of echo recordings. To estimate the abundance, the allocated NASC -values were averaged for ICES-squares (0.5° latitude by 1° longitude). For each statistical square, the unit area density of fish (rA) in number per square nautical mile (N*nm-2) was calculated using standard equations (Foote et al., 1987; Toresen et al., 1998). To estimate the total abundance of fish, the unit area abundance for each statistical square was multiplied by the number of square nautical miles in each statistical square and then summed for all the statistical squares within defined subareas and over the total area. Biomass estimation was calculated by multiplying abundance in numbers by the average weight of the fish in each statistical square then summing all squares within defined subareas and over the total area. The Norwegian BEAM soft-ware (Totland and Godø 2001) was used to make estimates of total biomass.

Spatially explicit estimates of length and biomass of Clupea harengus (Norwegian spring spawning herring) from acoustic and trawl surveys in the North-East Atlantic in May 2007

Acoustic estimates of herring and blue whiting abundance were obtained during the surveys using the Simrad ER60 scientific echosounder. The allocation of NASC-values to herring, blue whiting and other acoustic targets were based on the composition of the trawl catches and the appearance of echo recordings. To estimate the abundance, the allocated NASC -values were averaged for ICES-squares (0.5° latitude by 1° longitude). For each statistical square, the unit area density of fish (rA) in number per square nautical mile (N*nm-2) was calculated using standard equations (Foote et al., 1987; Toresen et al., 1998). To estimate the total abundance of fish, the unit area abundance for each statistical square was multiplied by the number of square nautical miles in each statistical square and then summed for all the statistical squares within defined subareas and over the total area. Biomass estimation was calculated by multiplying abundance in numbers by the average weight of the fish in each statistical square then summing all squares within defined subareas and over the total area. The Norwegian BEAM soft-ware (Totland and Godø 2001) was used to make estimates of total biomass.

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.