Seal census raw data of campaigns EMAGE I to V between 1996 and 2000 with links to datasets

The development of models of marine ecosystems in the Southern Ocean is becoming increasingly important as a means of understanding and managing impacts such as exploitation and climate change. Collating data from disparate sources, and understanding biases or uncertainties inherent in those data, are important first steps for improving ecosystem models. This review focuses on seals that breed in ice habitats of the Southern Ocean (i.e. the crabeater seal, Lobodon carcinophaga; Ross seal, Ommatophoca rossii; leopard seal, Hydrurga leptonyx; and Weddell seal, Leptonychotes weddellii). Data on populations (abundance and trends in abundance), distribution and habitat use (movement, key habitat and environmental features) and foraging (diet) are summarised, and potential biases and uncertainties inherent in those data are identified and discussed. Spatial and temporal gaps in knowledge of the populations, habitats and diet of each species are also identified.

X-ray diffraction analysis, binocular component analysis and radiocarbon dating of sediment cores from the NW Iberian shelf system

Based on a well-established stratigraphic framework and 47 AMS-14C dated sediment cores, the distribution of facies types on the NW Iberian margin is analysed in response to the last deglacial sea-level rise, thus providing a case study on the sedimentary evolution of a high-energy, low-accumulation shelf system. Altogether, four main types of sedimentary facies are defined. (1) A gravel-dominated facies occurs mostly as time-transgressive ravinement beds, which initially developed as shoreface and storm deposits in shallow waters on the outer shelf during the last sea-level lowstand; (2) A widespread, time-transgressive mixed siliceous/biogenic-carbonaceous sand facies indicates areas of moderate hydrodynamic regimes, high contribution of reworked shelf material, and fluvial supply to the shelf; (3) A glaucony-containing sand facies in a stationary position on the outer shelf formed mostly during the last-glacial sea-level rise by reworking of older deposits as well as authigenic mineral formation; and (4) A mud facies is mostly restricted to confined Holocene fine-grained depocentres, which are located in mid-shelf position. The observed spatial and temporal distribution of these facies types on the high-energy, low-accumulation NW Iberian shelf was essentially controlled by the local interplay of sediment supply, shelf morphology, and strength of the hydrodynamic system. These patterns are in contrast to high-accumulation systems where extensive sediment supply is the dominant factor on the facies distribution. This study emphasises the importance of large-scale erosion and material recycling on the sedimentary buildup during the deglacial drowning of the shelf. The presence of a homogenous and up to 15-m thick transgressive cover above a lag horizon contradicts the common assumption of sparse and laterally confined sediment accumulation on high-energy shelf systems during deglacial sea-level rise. In contrast to this extensive sand cover, laterally very confined and maximal 4-m thin mud depocentres developed during the Holocene sea-level highstand. This restricted formation of fine-grained depocentres was related to the combination of: (1) frequently occurring high-energy hydrodynamic conditions; (2) low overall terrigenous input by the adjacent rivers; and (3) the large distance of the Galicia Mud Belt to its main sediment supplier.

Mapped cetacean habitat suitability and richness in the Azores, links to ArcGIS files

Marine spatial planning and ecological research call for high-resolution species distribution data. However, those data are still not available for most marine large vertebrates. The dynamic nature of oceanographic processes and the wide-ranging behavior of many marine vertebrates create further difficulties, as distribution data must incorporate both the spatial and temporal dimensions. Cetaceans play an essential role in structuring and maintaining marine ecosystems and face increasing threats from human activities. The Azores holds a high diversity of cetaceans but the information about spatial and temporal patterns of distribution for this marine megafauna group in the region is still very limited. To tackle this issue, we created monthly predictive cetacean distribution maps for spring and summer months, using data collected by the Azores Fisheries Observer Programme between 2004 and 2009. We then combined the individual predictive maps to obtain species richness maps for the same period. Our results reflect a great heterogeneity in distribution among species and within species among different months. This heterogeneity reflects a contrasting influence of oceanographic processes on the distribution of cetacean species. However, some persistent areas of increased species richness could also be identified from our results. We argue that policies aimed at effectively protecting cetaceans and their habitats must include the principle of dynamic ocean management coupled with other area-based management such as marine spatial planning.

Fitting model-based psychometric functions to simultaneity and temporal-order judgment data: MATLAB and R routines

Research on temporal-order perception uses temporal-order judgment (TOJ) tasks or synchrony judgment (SJ) tasks in their binary SJ2 or ternary SJ3 variants. In all cases, two stimuli are presented with some temporal delay, and observers judge the order of presentation. Arbitrary psychometric functions are typically fitted to obtain performance measures such as sensitivity or the point of subjective simultaneity, but the parameters of these functions are uninterpretable. We describe routines in MATLAB and R that fit model-based functions whose parameters are interpretable in terms of the processes underlying temporal-order and simultaneity judgments and responses. These functions arise from an independent-channels model assuming arrival latencies with exponential distributions and a trichotomous decision space. Different routines fit data separately for SJ2, SJ3, and TOJ tasks, jointly for any two tasks, or also jointly for the three tasks (for common cases in which two or even the three tasks were used with the same stimuli and participants). Additional routines provide bootstrap p-values and confidence intervals for estimated parameters. A further routine is included that obtains performance measures from the fitted functions. An R package for Windows and source code of the MATLAB and R routines are available as Supplementary Files.

On the discrepant results in synchrony judgment and temporal-order judgment tasks: a quantitative model.

Research on the perception of temporal order uses either temporal-order judgment (TOJ) tasks or synchrony judgment (SJ) tasks, in both of which two stimuli are presented with some temporal delay and observers must judge the order of presentation. Results generally differ across tasks, raising concerns about whether they measure the same processes. We present a model including sensory and decisional parameters that places these tasks in a common framework that allows studying their implications on observed performance. TOJ tasks imply specific decisional components that explain the discrepancy of results obtained with TOJ and SJ tasks. The model is also tested against published data on audiovisual temporal-order judgments, and the fit is satisfactory, although model parameters are more accurately estimated with SJ tasks. Measures of latent point of subjective simultaneity and latent sensitivity are defined that are invariant across tasks by isolating the sensory parameters governing observed performance, whereas decisional parameters vary across tasks and account for observed differences across them. Our analyses concur with other evidence advising against the use of TOJ tasks in research on perception of temporal order.

A geographical perspective on the decline and extermination of the Irish wolf Canis lupus - an initial assessment

Wolves were a component of the Irish landscape until 1786 when the last one was killed. It had taken a concerted effort by Cromwell and his Government in Ireland to bring this about particularly through deforestation and landscape change, legislation, bounties and the efforts of a few professional wolf hunters. This paper estimates the wolf population in Ireland at three lime periods in the 1600s and examines how each of the forces already mentioned led to their eventual extermination. The 87 dated and documented wolf incidents which include wolf attacks on both animals and humans, wolf observations and the hunting and killing of wolves over the period 1560-1789 show both spatial and temporal variations.

Mechanochemistry for Active Materials and Devices

The coupling of mechanical stress fields in polymers to covalent chemistry (polymer mechanochemistry) has provided access to previously unattainable chemical reactions and polymer transformations. In the bulk, mechanochemical activation has been used as the basis for new classes of stress-responsive polymers that demonstrate stress/strain sensing, shear-induced intermolecular reactivity for molecular level remodeling and self-strengthening, and the release of acids and other small molecules that are potentially capable of triggering further chemical response. The potential utility of polymer mechanochemistry in functional materials is limited, however, by the fact that to date, all reported covalent activation in the bulk occurs in concert with plastic yield and deformation, so that the structure of the activated object is vastly different from its nascent form. Mechanochemically activated materials have thus been limited to “single use” demonstrations, rather than as multi-functional materials for structural and/or device applications. Here, we report that filled polydimethylsiloxane (PDMS) elastomers provide a robust elastic substrate into which mechanophores can be embedded and activated under conditions from which the sample regains its original shape and properties. Fabrication is straightforward and easily accessible, providing access for the first time to objects and devices that either release or reversibly activate chemical functionality over hundreds of loading cycles.

While the mechanically accelerated ring-opening reaction of spiropyran to merocyanine and associated color change provides a useful method by which to image the molecular scale stress/strain distribution within a polymer, the magnitude of the forces necessary for activation had yet to be quantified. Here, we report single molecule force spectroscopy studies of two spiropyran isomers. Ring opening on the timescale of tens of milliseconds is found to require forces of ~240 pN, well below that of previously characterized covalent mechanophores. The lower threshold force is a combination of a low force-free activation energy and the fact that the change in rate with force (activation length) of each isomer is greater than that inferred in other systems. Importantly, quantifying the magnitude of forces required to activate individual spiropyran-based force-probes enables the probe behave as a “scout” of molecular forces in materials; the observed behavior of which can be extrapolated to predict the reactivity of potential mechanophores within a given material and deformation.

We subsequently translated the design platform to existing dynamic soft technologies to fabricate the first mechanochemically responsive devices; first, by remotely inducing dielectric patterning of an elastic substrate to produce assorted fluorescent patterns in concert with topological changes; and second, by adopting a soft robotic platform to produce a color change from the strains inherent to pneumatically actuated robotic motion. Shown herein, covalent polymer mechanochemistry provides a viable mechanism to convert the same mechanical potential energy used for actuation into value-added, constructive covalent chemical responses. The color change associated with actuation suggests opportunities for not only new color changing or camouflaging strategies, but also the possibility for simultaneous activation of latent chemistry (e.g., release of small molecules, change in mechanical properties, activation of catalysts, etc.) in soft robots. In addition, mechanochromic stress mapping in a functional actuating device might provide a useful design and optimization tool, revealing spatial and temporal force evolution within the actuator in a way that might also be coupled to feedback loops that allow autonomous, self-regulation of activity.

In the future, both the specific material and the general approach should be useful in enriching the responsive functionality of soft elastomeric materials and devices. We anticipate the development of new mechanophores that, like the materials, are reversibly and repeatedly activated, expanding the capabilities of soft, active devices and further permitting dynamic control over chemical reactivity that is otherwise inaccessible, each in response to a single remote signal.

Retrieval of grassland biophysical parameters using multitemporal optical and radar satellite data

The amount and quality of available biomass is a key factor for the sustainable livestock industry and agricultural management related decision making. Globally 31.5% of land cover is grassland while 80% of Ireland’s agricultural land is grassland. In Ireland, grasslands are intensively managed and provide the cheapest feed source for animals. This dissertation presents a detailed state of the art review of satellite remote sensing of grasslands, and the potential application of optical (Moderate–resolution Imaging Spectroradiometer (MODIS)) and radar (TerraSAR-X) time series imagery to estimate the grassland biomass at two study sites (Moorepark and Grange) in the Republic of Ireland using both statistical and state of the art machine learning algorithms. High quality weather data available from the on-site weather station was also used to calculate the Growing Degree Days (GDD) for Grange to determine the impact of ancillary data on biomass estimation. In situ and satellite data covering 12 years for the Moorepark and 6 years for the Grange study sites were used to predict grassland biomass using multiple linear regression, Neuro Fuzzy Inference Systems (ANFIS) models. The results demonstrate that a dense (8-day composite) MODIS image time series, along with high quality in situ data, can be used to retrieve grassland biomass with high performance (R2 = 0:86; p < 0:05, RMSE = 11.07 for Moorepark). The model for Grange was modified to evaluate the synergistic use of vegetation indices derived from remote sensing time series and accumulated GDD information. As GDD is strongly linked to the plant development, or phonological stage, an improvement in biomass estimation would be expected. It was observed that using the ANFIS model the biomass estimation accuracy increased from R2 = 0:76 (p < 0:05) to R2 = 0:81 (p < 0:05) and the root mean square error was reduced by 2.72%. The work on the application of optical remote sensing was further developed using a TerraSAR-X Staring Spotlight mode time series over the Moorepark study site to explore the extent to which very high resolution Synthetic Aperture Radar (SAR) data of interferometrically coherent paddocks can be exploited to retrieve grassland biophysical parameters. After filtering out the non-coherent plots it is demonstrated that interferometric coherence can be used to retrieve grassland biophysical parameters (i. e., height, biomass), and that it is possible to detect changes due to the grass growth, and grazing and mowing events, when the temporal baseline is short (11 days). However, it not possible to automatically uniquely identify the cause of these changes based only on the SAR backscatter and coherence, due to the ambiguity caused by tall grass laid down due to the wind. Overall, the work presented in this dissertation has demonstrated the potential of dense remote sensing and weather data time series to predict grassland biomass using machine-learning algorithms, where high quality ground data were used for training. At present a major limitation for national scale biomass retrieval is the lack of spatial and temporal ground samples, which can be partially resolved by minor modifications in the existing PastureBaseIreland database by adding the location and extent ofeach grassland paddock in the database. As far as remote sensing data requirements are concerned, MODIS is useful for large scale evaluation but due to its coarse resolution it is not possible to detect the variations within the fields and between the fields at the farm scale. However, this issue will be resolved in terms of spatial resolution by the Sentinel-2 mission, and when both satellites (Sentinel-2A and Sentinel-2B) are operational the revisit time will reduce to 5 days, which together with Landsat-8, should enable sufficient cloud-free data for operational biomass estimation at a national scale. The Synthetic Aperture Radar Interferometry (InSAR) approach is feasible if there are enough coherent interferometric pairs available, however this is difficult to achieve due to the temporal decorrelation of the signal. For repeat-pass InSAR over a vegetated area even an 11 days temporal baseline is too large. In order to achieve better coherence a very high resolution is required at the cost of spatial coverage, which limits its scope for use in an operational context at a national scale. Future InSAR missions with pair acquisition in Tandem mode will minimize the temporal decorrelation over vegetation areas for more focused studies. The proposed approach complements the current paradigm of Big Data in Earth Observation, and illustrates the feasibility of integrating data from multiple sources. In future, this framework can be used to build an operational decision support system for retrieval of grassland biophysical parameters based on data from long term planned optical missions (e. g., Landsat, Sentinel) that will ensure the continuity of data acquisition. Similarly, Spanish X-band PAZ and TerraSAR-X2 missions will ensure the continuity of TerraSAR-X and COSMO-SkyMed.

Selection of a network of large lakes and reservoirs suitable for global environmental change analysis using Earth Observation

The GloboLakes project, a global observatory of lake responses to environmental change, aims to exploit current satellite missions and long remote-sensing archives to synoptically study multiple lake ecosystems, assess their current condition, reconstruct past trends to system trajectories, and assess lake sensitivity to multiple drivers of change. Here we describe the selection protocol for including lakes in the global observatory based upon remote-sensing techniques and an initial pool of the largest 3721 lakes and reservoirs in the world, as listed in the Global Lakes and Wetlands Database. An 18-year-long archive of satellite data was used to create spatial and temporal filters for the identification of waterbodies that are appropriate for remote-sensing methods. Further criteria were applied and tested to ensure the candidate sites span a wide range of ecological settings and characteristics; a total 960 lakes, lagoons, and reservoirs were selected. The methodology proposed here is applicable to new generation satellites, such as the European Space Agency Sentinel-series.

Sources, Fate and Transformation of Organic Matter in Wetlands and Estuaries

Dissolved organic matter (DOM) is a complex mixture of organic compounds and represents the largest reservoirs of carbon (C) on earth. Particulate organic matter (POM) is another important carbon component in C cycling and controls a variety of biogeochemical processes. Estuaries, as important interfaces between land and ocean, play important roles in retaining and transforming such organic matter (OM) and serve as both sources and sinks of DOM and POM. There is a diverse array of both autochthonous and allochthonous OM sources in wetland/estuarine ecosystems. A comprehensive study on the sources, transformation and fate of OM in such ecosystems is essential in advancing our understanding of C cycling and better constraining the global C budget. In this work, DOM characteristics were investigated in different estuaries. Dissolved organic matter source strengths and dynamics were assessed in a seagrass-dominated subtropical estuarine lagoon. DOM dynamics controlled by hydrology and seagrass primary productivity were confirmed, and the primary source of DOM was quantified using the combination of excitation emission matrix fluorescence with parallel factor analysis (EEM-PARAFAC) and stable C isotope analysis. Seagrass can contribute up to 72% of the DOM in the study area. The spatial and temporal variation of DOM dynamics was also studied in a freshwated dominated estuary fringed with extensive salt marshes. The data showed that DOM was primarily derived from freshwater marshes and controlled by hydrology while salt marsh plants play a significant role in structuring the distribution patterns of DOM quality and quantity. The OM dynamics was also investigated in a mangrove-dominate estuary and a comparative study was conducted between the DOM and POM pools. The results revealed both similarity and dissimilarity in DOM and POM composition. The dynamics of both OM pools are largely uncoupled as a result of source differences. Fringe mangrove swamps are suggested to export similar amounts of DOM and POM and should be considered as an important source in coastal C budgets. Lastly, chemical characterizations were conducted on the featured fluorescence component in OM in an attempt to better understand the composition and origins of the specific PARAFAC component. The traditionally defined ‘protein-like’ fluorescence was found to contain both proteinaceous and phenolic compounds, suggesting that the application of this parameter as a proxy for amino acid content and bioavailability may be limited.

Carbon and Nitrogen Stable Isotopic Patterns in South Florida Coastal Ecosystems: Modern and Paleoceanographic Perspectives

Long term management plans for restoration of natural flow conditions through the Everglades increase the importance of understanding potential nutrient impacts of increased freshwater delivery on coastal biogeochemistry. The present study sought to increase understanding of the coastal marine system of South Florida under modern conditions and through the anthropogenic changes in the last century, on scales ranging from individual nutrient cycle processes to seasonal patterns in organic material (OM) under varying hydrodynamic regime, to century scale analysis of sedimentary records. In all applications, carbon and nitrogen stable isotopic compositions of OM were examined as natural recorders of change and nutrient cycling in the coastal system. High spatial and temporal variability in stable isotopic compositions were observed on all time scales. During a transient phytoplankton bloom, ä15N values suggested nitrogen fixation as a nutrient source supporting enhanced productivity. Seasonally, particulate organic material (POM) from ten sites along the Florida Reef Tract and in Florida Bay demonstrated variable fluctuations dependent on hydrodynamic setting. Three separate intra-annual patterns were observed, yet statistical differences were observed between groupings of Florida Bay and Atlantic Ocean sites. The POM ä15N values ranged on a quarterly basis by 7‰, while ä13C varied by 22‰. From a sediment history perspective, four cores collected from Florida Bay further demonstrated the spatial and temporal variability of the system in isotopic composition of bulk OM over time. Source inputs of OM varied with location, with terrestrial inputs dominating proximal to Everglades freshwater discharge, seagrasses dominating in open estuary cores, and a marine mixture of phytoplankton and seagrass in a core from the boundary zone between Florida Bay and the Gulf of Mexico. Significant shifts in OM geochemistry were observed coincident with anthropogenic events of the 20th century, including railroad and road construction in the Florida Keys and Everglades, and also the extensive drainage changes in Everglades hydrology. The sediment record also preserved evidence of the major hurricanes of the last century, with excursions in geochemical composition coincident with Category 4-5 storms.

Floral and faunal characteristics and content of photosynthetic pigment in soils of northern Victoria Land, Antarctica

Although soil algae are among the main primary producers in most terrestrial ecosystems of continental Antarctica, there are very few quantitative studies on their relative proportion in the main algal groups and on how their distribution is affected by biotic and abiotic factors. Such knowledge is essential for understanding the functioning of Antarctic terrestrial ecosystems. We therefore analyzed biological soil crusts from northern Victoria Land to determine their pH, electrical conductivity (EC), water content (W), total and organic C (TC and TOC) and total N (TN) contents, and the presence and abundance of photosynthetic pigments. In particular, the latter were tested as proxies for biomass and coarse-resolution community structure. Soil samples were collected from five sites with known soil algal communities and the distribution of pigments was shown to reflect differences in the relative proportions of Chlorophyta, Cyanophyta and Bacillariophyta in these sites. Multivariate and univariate models strongly indicated that almost all soil variables (EC, W, TOC and TN) were important environmental correlates of pigment distribution. However, a significant amount of variation is independent of these soil variables and may be ascribed to local variability such as changes in microclimate at varying spatial and temporal scales. There are at least five possible sources of local variation: pigment preservation, temporal variations in water availability, temporal and spatial interactions among environmental and biological components, the local-scale patchiness of organism distribution, and biotic interactions.

Biomarker distributions in surface sediments from the Central Arctic Ocean and adjacent marginal seas

Records of the spatial and temporal variability of Arctic Ocean sea ice are of significance for understanding the causes of the dramatic decrease in Arctic sea-ice cover of recent years. In this context, the newly developed sea-ice proxy IP25, a mono-unsaturated highly branched isoprenoid alkene with 25 carbon atoms biosynthesized specifically by sea-ice associated diatoms and only found in Arctic and sub-Arctic marine sediments, has been used to reconstruct the recent spatial sea-ice distribution. The phytoplankton biomarkers 24S-brassicasterol and dinosterol were determined alongside IP25 to distinguish ice-free or permanent ice conditions, and to estimate the sea-ice conditions semi-quantitatively by means of the phytoplankton-IP25 index (PIP25). Within our study, for the first time a comprehensive data set of these biomarkers was produced using fresh and deep-frozen surface sediment samples from the Central Arctic Ocean proper (>80°N latitude) characterised by a permanent ice cover today and recently obtained surface sediment samples from the Chukchi Plateau/Basin partly covered by perennial sea ice. In addition, published and new data from other Arctic and sub-Arctic regions were added to generate overview distribution maps of IP25 and phytoplankton biomarkers across major parts of the modern Arctic Ocean. These comprehensive biomarker data indicate perennial sea-ice cover in the Central Arctic, ice-free conditions in the Barents Sea and variable sea-ice situations in other marginal seas. The low but more than zero values of biomarkers in the Central Arctic supported the low in-situ productivity there. The PIP25 index values reflect modern sea-ice conditions better than IP25 alone and show a positive correlation with spring/summer sea ice. When calculating and interpreting PIP25 index as a (semi-quantitative) proxy for reconstructions of present and past Arctic sea-ice conditions from different Arctic/sub-Arctic areas, information of the source of phytoplankton biomarkers and the possible presence of allochthonous biomarkers is needed, and the records of the individual biomarkers always should be considered as well.

(Table 1) Age control points and dating methods for sediment cores of the Sea of Okhotsk

Marine- and terrestrial-derived biomarkers (alkenones, brassicasterol, dinosterol, and long-chain n-alkanes), as well as carbonate, biogenic opal, and ice-rafted debris (IRD), were measured in two sediment cores in the Sea of Okhotsk, which is located in the northwestern Pacific rim and characterized by high primary productivity. Down-core profiles of phytoplankton markers suggest that primary productivity abruptly increased during the global Meltwater Pulse events 1A (about 14 ka) and 1B (about 11 ka) and stayed high in the Holocene. Spatial and temporal distributions of the phytoplankton productivity were found to be consistent with changes in the reconstructed sea ice distribution on the basis of the IRD. This demonstrates that the progress and retreat of sea ice regulated primary productivity in the Sea of Okhotsk with minimum productivity during the glacial period. The mass accumulation rates of alkenones, CaCO3, and biogenic opal indicate that the dominant phytoplankton species during deglaciation was the coccolithophorid, Emiliania huxleyi, which was replaced by diatoms in the late Holocene. Such a phytoplankton succession was probably caused by an increase in silicate supply to the euphotic layer, possibly associated with a change in surface hydrography and/or linked to enhanced upwelling of North Pacific Deep Water.

(Table 1) Age control and temporal resolution of paleoceanographic records of the North Atlantic

Short-term changes in sea surface conditions controlling the thermohaline circulation in the northern North Atlantic are expected to be especially efficient in perturbing global climate stability. Here we assess past variability of sea surface temperature (SST) in the northeast Atlantic and Norwegian Sea during Marine Isotope Stage (MIS) 2 and, in particular, during the Last Glacial Maximum (LGM). Five high-resolution SST records were established on a meridional transect (53°N-72°N) to trace centennial-scale oscillations in SST and sea-ice cover. We used three independent computational techniques (SIMMAX modern analogue technique, Artificial Neural Networks (ANN), and Revised Analog Method (RAM)) to reconstruct SST from planktonic foraminifer census counts. SIMMAX and ANN reproduced short-term SST oscillations of similar magnitude and absolute levels, while RAM, owing to a restrictive analog selection, appears less suitable for reconstructing "cold end" SST. The SIMMAX and ANN SST reconstructions support the existence of a weak paleo-Norwegian Current during Dansgaard-Oeschger (DO) interstadials number 4, 3, 2, and 1. During the LGM, two warm incursions of 7°C water to occurred in the northern North Atlantic but ended north of the Iceland Faroe Ridge. A rough numerical estimate shows that the near-surface poleward heat transfer from 53° across the Iceland-Faroe Ridge up to to 72° N dropped to less than 60% of the modern value during DO interstadials and to almost zero during DO stadials. Summer sea ice was generally confined to the area north of 70°N and only rarely expanded southward along the margins of continental ice sheets. Internal LGM variability of North Atlantic (>40°N) SST in the GLAMAP 2000 compilation (Sarnthein et al., 2003, doi:10.1029/2002PA000771; Pflaumann et al., 2003, doi:10.1029/2002PA000774) indicates maximum instability in the glacial subpolar gyre and at the Iberian Margin, while in the Nordic Seas, SST was continuously low.