Daily polynya area and ice production in the North Water polynya region

The North Water (NOW) Polynya is a regularly-forming area of open-water and thin-ice, located between northwestern Greenland and Ellesmere Island (Canada) at the northern tip of Baffin Bay. Due to its large spatial extent, it is of high importance for a variety of physical and biological processes, especially in wintertime. Here, we present a long-term remote sensing study for the winter seasons 1978/1979 to 2014/2015. Polynya characteristics are inferred from (1) sea ice concentrations and brightness temperatures from passive microwave satellite sensors (Advanced Microwave Scanning Radiometer (AMSR-E and AMSR2), Scanning Multichannel Microwave Radiometer (SMMR), Special Sensor Microwave Imager/Sounder (SSM/I-SSMIS)) and (2) thin-ice thickness distributions, which are calculated using MODIS ice-surface temperatures and European Center for Medium-Range Weather Forecasts (ECMWF) atmospheric reanalysis data in a 1D thermodynamic energy-balance model. Daily ice production rates are retrieved for each winter season from 2002/2003 to 2014/2015, assuming that all heat loss at the ice surface is balanced by ice growth. Two different cloud-cover correction schemes are applied on daily polynya area and ice production values to account for cloud gaps in the MODIS composites. Our results indicate that the NOW polynya experienced significant seasonal changes over the last three decades considering the overall frequency of polynya occurrences, as well as their spatial extent. In the 1980s, there were prolonged periods of a more or less closed ice cover in northern Baffin Bay in winter. This changed towards an average opening on more than 85% of the days between November and March during the last decade. Noticeably, the sea ice cover in the NOW polynya region shows signs of a later-appearing fall freeze-up, starting in the late 1990s. Different methods to obtain daily polynya area using passive microwave AMSR-E/AMSR2 data and SSM/I-SSMIS data were applied. A comparison with MODIS data (thin-ice thickness < 20 cm) shows that the wintertime polynya area estimates derived by MODIS are about 30 to 40% higher than those derived using the polynya signature simulation method (PSSM) with AMSR-E data. In turn, the difference in polynya area between PSSM and a sea ice concentration (SIC) threshold of 70% is fairly low (approximately 10%) when applied to AMSR-E data. For the coarse-resolution SSM/I-SSMIS data, this difference is much larger, particularly in November and December. Instead of a sea ice concentration threshold, the PSSM method should be used for SSM/I-SSMIS data. Depending on the type of cloud-cover correction, the calculated ice production based on MODIS data reaches an average value of 264.4 ± 65.1 km**3 to 275.7 ± 67.4 km**3 (2002/2003 to 2014/2015) and shows a high interannual variability. Our achieved long-term results underline the major importance of the NOW polynya considering its influence on Arctic ice production and associated atmosphere/ocean processes.

XCO2 and XCH4 total column measurement during POLARSTERN cruise PS83 (ANT-XXIX/10), north-south gradient

A portable Fourier transform spectrometer (FTS), model EM27/SUN, was deployed onboard the research vessel Polarstern to measure the column-average dry air mole fractions of carbon dioxide (XCO2) and methane (XCH4) by means of direct sunlight absorption spectrometry. We report on technical developments as well as data calibration and reduction measures required to achieve the targeted accuracy of fractions of a percent in retrieved XCO2 and XCH4 while operating the instrument under field conditions onboard the moving platform during a 6-week cruise on the Atlantic from Cape Town (South Africa, 34° S, 18° E; 5 March 2014) to Bremerhaven (Germany, 54° N, 19° E; 14 April 2014). We demonstrate that our solar tracker typically achieved a tracking precision of better than 0.05° toward the center of the sun throughout the ship cruise which facilitates accurate XCO2 and XCH4 retrievals even under harsh ambient wind conditions. We define several quality filters that screen spectra, e.g., when the field of view was partially obstructed by ship structures or when the lines-of-sight crossed the ship exhaust plume. The measurements in clean oceanic air, can be used to characterize a spurious air-mass dependency. After the campaign, deployment of the spectrometer alongside the TCCON (Total Carbon Column Observing Network) instrument at Karlsruhe, Germany, allowed for determining a calibration factor that makes the entire campaign record traceable to World Meteorological Organization (WMO) standards. Comparisons to observations of the GOSAT satellite and concentration fields modeled by the European Centre for Medium-Range Weather Forecasts (ECMWF) Copernicus Atmosphere Monitoring Service (CAMS) demonstrate that the observational setup is well suited to provide validation opportunities above the ocean and along interhemispheric transects.

Precipitation over sea: Validation against shipboard in-situ measurements

The satellite derived HOAPS (Hamburg Ocean Atmosphere Parameters and Fluxes from Satellite data) and ECMWF (European Centre for Medium-Range Weather Forecasts) ERA-Interim reanalysis data sets have been validated against in-situ precipitation measurements from ship rain gauges and optical disdrometers over the open-ocean by applying a statistical analysis for binary forecasts. For this purpose collocated pairs of data were merged within a certain temporal and spatial threshold into single events, according to the satellites' overpass, the observation and the forecast times. HOAPS detects the frequency of precipitation well, while ERA-Interim strongly overestimates it, especially in the tropics and subtropics. Although precipitation rates are difficult to compare because along-track point measurements are collocated with areal estimates and the numbers of available data are limited, we find that HOAPS underestimates precipitation rates, while ERA-Interim's Atlantic-wide average precipitation rate is close to measurements. However, regionally averaged over latitudinal belts, there are deviations between the observed mean precipitation rates and ERA-Interim. The most obvious ERA-Interim feature is an overestimation of precipitation in the area of the intertropical convergence zone and the southern sub-tropics over the Atlantic Ocean. For a limited number of snow measurements by optical disdrometers it can be concluded that both HOAPS and ERA-Interim are suitable to detect the occurrence of solid precipitation.

Land Surface Temperature continuous time series from SSM/I passive microwave radiometers, over northern latitudes (Latitudes > 45°N), during summer snow free period (2000-2011), with links to NetCDF files

Remote sensing instruments are key players to map land surface temperature (LST) at large temporal and spatial scales. In this paper, we present how we combine passive microwave and thermal infrared data to estimate LST during summer snow-free periods over northern high latitudes. The methodology is based on the SSM/I-SSMIS 37 GHz measurements at both vertical and horizontal polarizations on a 25 km × 25 km grid size. LST is retrieved from brightness temperatures introducing an empirical linear relationship between emissivities at both polarizations as described in Royer and Poirier (2010). This relationship is calibrated at pixel scale, using cloud-free independent LST data from MODIS instruments. The SSM/I-SSMIS and MODIS data are synchronized by fitting a diurnal cycle model built on skin temperature reanalysis provided by the European Centre for Medium-Range Weather Forecasts (ECMWF). The resulting temperature dataset is provided at 25 km scale and at an hourly time step during the ten-year analysis period (2000-2011). This new product was locally evaluated at five experimental sites of the EU-PAGE21 project against air temperature measurements and meteorological model reanalysis, and compared to the MODIS LST product at both local and circumpolar scale. The results giving a mean RMSE of the order of 2.2 K demonstrate the usefulness of the microwave product, which is unaffected by clouds as opposed to thermal infrared products and offers a better resolution compared to model reanalysis.

Helicopter-borne sea ice thickness measurements in the Storfjord, Spitsbergen from flight HEM_STJ06_01

Results from electromagnetic induction surveys of sea-ice thickness in Storfjorden, Svalbard, reveal large interannual ice-thickness variations in a region which is typically characterized by a reoccurring polynya. The surveys were performed in March 2003, May 2006 and March 2007 with helicopter- and ship-based sensors. The thickness distributions are influenced by sea-ice and atmospheric boundary conditions 2 months prior to the surveys, which are assessed with synthetic aperture radar (SAR) images, regional QuikSCAT backscatter maps and wind information from the European Centre for Medium-Range Weather Forecasts (ECMWF) reanalysis dataset. Locally formed thin ice from the Storfjorden polynya was frequently observed in 2003 and 2007 (mean thickness 0.55 and 0.37 m, respectively) because these years were characterized by prevailing northeasterly winds. In contrast, the entire fjord was covered with thick external sea ice in 2006 (mean thickness 2.21 m), when ice from the Barents Sea was driven into the fjord by predominantly southerly winds. The modal thickness of this external ice in 2006 increased from 1.2 m in the northern fjord to 2.4 m in the southern fjord, indicating stronger deformation in the southern part. This dynamically thickened ice was even thicker than multi-year ice advected from the central Arctic Ocean in 2003 (mean thickness 1.83 m). The thermodynamic ice thickness of fast ice as boundary condition is investigated with a one-dimensional sea-ice growth model (1DICE) forced with meteorological data from the weather station at the island of Hopen, southeast of Storfjorden. The model results are in good agreement with the modal thicknesses of fast-ice measurements in all years.

Sources and seasonality of polycyclic aromatic hydrocarbons transported to Alert, Canadian Arctic

A probabilistic function (integrated source contribution function, ISCF) based on backward air mass trajectory calculation was developed to track sources and atmospheric pathways of polycyclic aromatic hydrocarbons (PAHs) to the Canadian High Arctic station of Alert. In addition to the movement of air masses, the emission intensities at the sources and the major processes of partition, indirect photolysis, and deposition occurring on the way to the Arctic were incorporated into the ISCF. The predicted temporal trend of PAHs at Alert was validated by measured PAH concentrations throughout 2004. The PAH levels in the summer are orders of magnitude lower than those in the winter and spring when long-range atmospheric transport events occur more frequently. PAHs observed at Alert are mostly from East Asia (including Russia Far East), North Europe (including European Russia), and North America. These sources account for 25, 45, and 27% of PAHs atmospheric level at Alert, respectively. Source regions and transport pathways contributing to the PAHs contamination in the Canadian High Arctic vary seasonally. In the winter, Russia and Europe are the major sources. PAHs from these sources travel eastward and turn to the north at approximately 120°E before reaching Alert, in conjunction with the well- known Arctic haze events. In the spring, PAHs from Russia and Europe first migrate to the west and then turn to the north at 60°W toward Alert. The majority of PAHs in the summer are from northern Canada where they are carried to Alert via low- level transport pathways. In the fall, 70% of PAHs arriving at Alert are delivered from North American sources.

Long-term chemical characterization of aerosols at CVAO from 2007 to 2011

The first long-term aerosol sampling and chemical characterization results from measurements at the Cape Verde Atmospheric Observatory (CVAO) on the island of São Vicente are presented and are discussed with respect to air mass origin and seasonal trends. In total 671 samples were collected using a high-volume PM10 sampler on quartz fiber filters from January 2007 to December 2011. The samples were analyzed for their aerosol chemical composition, including their ionic and organic constituents. Back trajectory analyses showed that the aerosol at CVAO was strongly influenced by emissions from Europe and Africa, with the latter often responsible for high mineral dust loading. Sea salt and mineral dust dominated the aerosol mass and made up in total about 80% of the aerosol mass. The 5-year PM10 mean was 47.1 ± 55.5 µg/m**2, while the mineral dust and sea salt means were 27.9 ± 48.7 and 11.1 ± 5.5 µg/m**2, respectively. Non-sea-salt (nss) sulfate made up 62% of the total sulfate and originated from both long-range transport from Africa or Europe and marine sources. Strong seasonal variation was observed for the aerosol components. While nitrate showed no clear seasonal variation with an annual mean of 1.1 ± 0.6 µg/m**3, the aerosol mass, OC (organic carbon) and EC (elemental carbon), showed strong winter maxima due to strong influence of African air mass inflow. Additionally during summer, elevated concentrations of OM were observed originating from marine emissions. A summer maximum was observed for non-sea-salt sulfate and was connected to periods when air mass inflow was predominantly of marine origin, indicating that marine biogenic emissions were a significant source. Ammonium showed a distinct maximum in spring and coincided with ocean surface water chlorophyll a concentrations. Good correlations were also observed between nss-sulfate and oxalate during the summer and winter seasons, indicating a likely photochemical in-cloud processing of the marine and anthropogenic precursors of these species. High temporal variability was observed in both chloride and bromide depletion, differing significantly within the seasons, air mass history and Saharan dust concentration. Chloride (bromide) depletion varied from 8.8 ± 8.5% (62 ± 42%) in Saharan-dust-dominated air mass to 30 ± 12% (87 ± 11%) in polluted Europe air masses. During summer, bromide depletion often reached 100% in marine as well as in polluted continental samples. In addition to the influence of the aerosol acidic components, photochemistry was one of the main drivers of halogenide depletion during the summer; while during dust events, displacement reaction with nitric acid was found to be the dominant mechanism. Positive matrix factorization (PMF) analysis identified three major aerosol sources: sea salt, aged sea salt and long-range transport. The ionic budget was dominated by the first two of these factors, while the long-range transport factor could only account for about 14% of the total observed ionic mass.

Particle concentration and size distribution in the NGRIP ice core

A novel laser microparticle detector used in conjunction with continuous sample melting has provided a more than 1500 m long record of particle concentration and size distribution of the NGRIP ice core, covering continuously the period approximately from 9.5-100 kyr before present; measurements were at 1.65 m depth resolution, corresponding to approximately 35-200 yr. Particle concentration increased by a factor of 100 in the Last Glacial Maximum (LGM) compared to the Preboreal, and sharp variations of concentration occurred synchronously with rapid changes in the delta18O temperature proxy. The lognormal mode µ of the volume distribution shows clear systematic variations with smaller modes during warmer climates and coarser modes during colder periods. We find µ ~ 1.7 µm diameter during LGM and µ ~ 1.3 µm during the Preboreal. On timescales below several 100 years µ and the particle concentration exhibit a certain degree of independence present especially during warm periods, when µ generally is more variable. Using highly simplifying considerations for atmospheric transport and deposition of particles we infer that (1) the observed changes of µ in the ice largely reflect changes in the size of airborne particles above the ice sheet and (2) changes of µ are indicative of changes in long range atmospheric transport time. From the observed size changes we estimate shorter transit times by roughly 25% during LGM compared to the Preboreal. The associated particle concentration increase from more efficient long range transport is estimated to less than one order of magnitude.

Stable isotope ratios and chemistry on corals from the Caribbean Sea

Instrumental climate data are limited in length and only available with low spatial coverage before the middle of the 20th century. This is too short to reliably determine and interpret decadal and longer scale climate variability and to understand the underlying mechanisms with sufficient accuracy. A proper knowledge of past variability of the climate system is needed to assess the anthropogenic impact on climate and ecosystems, and also important with regard to long-range climate forecasting. Highly-resolved records of past climate variations that extend beyond pre-industrial times can significantly help to understand long-term climate changes and trends. Indirect information on past environmental and climatic conditions can be deduced from climate-sensitive proxies. Large colonies of massive growing tropical reef corals have been proven to sensitively monitor changes in ambient seawater. Rapid skeletal growth, typically ranging between several millimeters to centimeters per year, allows the development of proxy records at sub-seasonal resolution. Stable oxygen isotopic composition and trace elemental ratios incorporated in the aragonitic coral skeleton can reveal a detailed history of past environmental conditions, e.g., sea surface temperature (SST). In general, coral-based reconstructions from the tropical Atlantic region have lagged behind the extensive work published using coral records from the Indian and Pacific Oceans. Difficulties in the analysis of previously utilized coral archives from the Atlantic, typically corals of the genera Montastrea and Siderastrea, have so far exacerbated the production of long-term high-resolution proxy records. The objective of this study is the evaluation of massive fast-growing corals of the species Diploria strigosa as a new marine archive for climate reconstructions from the tropical Atlantic region. For this purpose, coral records from two study sites in the eastern Caribbean Sea (Guadeloupe, Lesser Antilles; and Archipelago Los Roques, Venezuela) were examined. At Guadeloupe, a century-long monthly resolved multi-proxy coral record was generated. Results present the first d18O (Sr/Ca)-SST calibration equations for the Atlantic braincoral Diploria strigosa, that are robust and consistent with previously published values using other coral species from different regions. Both proxies reflect local variability of SST on a sub-seasonal scale, which is a precondition for studying seasonally phase-locked climate variations, as well as track variability on a larger spatial scale (i.e., in the Caribbean and tropical North Atlantic). Coral Sr/Ca reliably records local annual to interannual temperature variations and is higher correlated to in-situ air temperature than to grid-SST. The warming calculated from coral Sr/Ca is concurrent with the strong surface temperature increase at the study site during the past decades. Proxy data show a close relationship to major climate signals from the tropical Pacific and North Atlantic (the El Niño Southern Oscillation (ENSO) and the North Atlantic Oscillation (NAO)) affecting the seasonal cycle of SST in the North Tropical Atlantic (NTA). Coral oxygen isotopes are also influenced by seawater d18O (d18Osw) which is linked to the hydrological cycle, and capture large-scale climate variability in the NTA region better than Sr/Ca. Results from a quantitative comparison between extreme events in the two most prominent modes of external forcing, namely the ENSO and NAO, and respective events recorded in seasonal coral d18O imply that SST variability at the study site is highly linked to Pacific and North Atlantic variability, by this means supporting the assumptions of observational- and model-based studies which suggest a strong impact of ENSO and NAO forcings onto the NTA region through a modulation of trade wind strength in winter. Results from different spectral analysis tools suggest that interannual climate variability recorded by the coral proxies is II largely dictated by Pacific ENSO forcing, whereas at decadal and longer timescales the influence of the NAO is dominan. tThe Archipelago Los Roques is situated in the southeastern Caribbean Sea, north of the Venezuelan coast. Year-to-year variations in monthly resolved coral d18O of a nearcentury- long Diploria strigosa record are significantly correlated with SST and show pronounced multidecadal variations. About half of the variance in coral d18O can be explained by variations in seawater d18O, which can be estimated by calculating the d18Oresidual via subtracting the SST component from measured coral d18O. The d18Oresidual and a regional precipitation index are highly correlated at low frequencies, suggesting that d18Osw variations are primarily atmospheric-driven. Warmer SSTs at Los Roques broadly coincide with higher precipitation in the southeastern Caribbean at multidecadal time scales, effectively strengthening the climate signal in the coral d18O record. The Los Roques coral d18O record displays a strong and statistically significant relationship to different indices of hurricane activity during the peak of the Atlantic hurricane season in boreal summer and is a particularly good indicator of decadal-multidecadal swings in the latter indices. In general, the detection of long-term changes and trends in Atlantic hurricane activity is hampered due to the limited length of the reliable instrumental record and the known inhomogeneity in the observational databases which result from changes in observing practice and technology over the years. The results suggest that coral-derived proxy data from Los Roques can be used to infer changes in past hurricane activity on timescales that extend well beyond the reliable record. In addition, the coral record exhibits a clear negative trend superimposed on the decadal to multidecadal cycles, indicating a significant warming and freshening of surface waters in the genesis region of tropical cyclones during the past decades. The presented coral d18O time series provides the first and, so far, longest continuous coral-based record of hurricane activity. It appears that the combination of both signals (SST and d18Osw) in coral d18O leads to an amplification of large-scale climate signals in the record, and makes coral d18O even a better proxy for hurricane activity than SST alone. Atlantic hurricane activity naturally exhibits strong multidecadal variations that are associated with the Atlantic Multidecadal Oscillation (AMO), the major mode of lowfrequency variability in the North Atlantic Ocean. However, the mechanisms underlying this multidecadal variability remain controversial, primarily because of the limited instrumental record. The Los Roques coral d18O displays strong multidecadal variability with a period of approximately 60 years that is closely related to the AMO, making the Archipelago Los Roques a very sensitive location for studying low-frequency climate variability in the Atlantic Ocean. In summary, the coral records presented in this thesis capture different key climate variables in the north tropical Atlantic region very well, indicating that fast-growing Diploria strigosa corals represent a promising marine archive for further proxy-based reconstructions of past climate variability on a range of time scales.

Constraints on the magnitude and patterns of ocean cooling at the Last Glacial Maximum - Supplementary material

Observation-based reconstructions of sea surface temperature from relatively stable periods in the past, such as the Last Glacial Maximum, represent an important means of constraining climate sensitivity and evaluating model simulations. The first quantitative global reconstruction of sea surface temperatures during the Last Glacial Maximum was developed by the Climate Long-Range Investigation, Mapping and Prediction (CLIMAP) project in the 1970s and 1980s. Since that time, several shortcomings of that earlier effort have become apparent. Here we present an updated synthesis of sea surface temperatures during the Last Glacial Maximum, rigorously defined as the period between 23 and 19 thousand years before present, from the Multiproxy Approach for the Reconstruction of the Glacial Ocean Surface (MARGO) project. We integrate microfossil and geochemical reconstructions of surface temperatures and include assessments of the reliability of individual records. Our reconstruction reveals the presence of large longitudinal gradients in sea surface temperature in all of the ocean basins, in contrast to the simulations of the Last Glacial Maximum climate available at present.

Crete fault scarp point cloud data

This paper assesses the along strike variation of active bedrock fault scarps using long range terrestrial laser scanning (t-LiDAR) data in order to determine the distribution behaviour of scarp height and the subsequently calculate long term throw-rates. Five faults on Cretewhich display spectacular limestone fault scarps have been studied using high resolution digital elevation model (HRDEM) data. We scanned several hundred square metres of the fault system including the footwall, fault scarp and hanging wall of the investigated fault segment. The vertical displacement and the dip of the scarp were extracted every metre along the strike of the detected fault segment based on the processed HRDEM. The scarp variability was analysed by using statistical and morphological methods. The analysis was done in a geographical information system (GIS) environment. Results show a normal distribution for the scanned fault scarp's vertical displacement. Based on these facts, the mean value of height was chosen to define the authentic vertical displacement. Consequently the scarp can be divided into above, below and within the range of mean (within one standard deviation) and quantify the modifications of vertical displacement. Therefore, the fault segment can be subdivided into areas which are influenced by external modification like erosion and sedimentation processes. Moreover, to describe and measure the variability of vertical displacement along strike the fault, the semi-variance was calculated with the variogram method. This method is used to determine how much influence the external processes have had on the vertical displacement. By combining of morphological and statistical results, the fault can be subdivided into areas with high external influences and areas with authentic fault scarps, which have little or no external influences. This subdivision is necessary for long term throw-rate calculations, because without this differentiation the calculated rates would be misleading and the activity of a fault would be incorrectly assessed with significant implications for seismic hazard assessment since fault slip rate data govern the earthquake recurrence. Furthermore, by using this workflow areas with minimal external influences can be determined, not only for throw-rate calculations, but also for determining samples sites for absolute dating techniques such as cosmogenic nuclide dating. The main outcomes of this study include: i) there is no direct correlation between the fault's mean vertical displacement and dip (R² less than 0.31); ii) without subdividing the scanned scarp into areas with differing amounts of external influences, the along strike variability of vertical displacement is ±35%; iii) when the scanned scarp is subdivided the variation of the vertical displacement of the authentic scarp (exposed by earthquakes only) is in a range of ±6% (the varies depending on the fault from 7 to 12%); iv) the calculation of the long term throw-rate (since 13 ka) for four scarps in Crete using the authentic vertical displacement is 0.35 ± 0.04 mm/yr at Kastelli 1, 0.31 ± 0.01 mm/yr at Kastelli 2, 0.85 ± 0.06 mm/yr at the Asomatos fault (Sellia) and 0.55 ± 0.05 mm/yr at the Lastros fault.

Relative abundance of pteropod species of the plankton pump samples in the surface water of the eastern North Atlantic Ocean (Table 4)

Plankton pump samples and plankton tows (size fractions between 0.04 mm and 1.01 mm) from the eastern North Atlantic Ocean contain the following shell- and skeleton-producing planktonic and nektonic organisms, which can be fossilized in the sediments: diatoms, radiolarians, foraminifers, pteropods, heteropods, larvae of benthic gastropods and bivalves, ostracods, and fish. The abundance of these components has been mapped quantitatively in the eastern North Atlantic surface waters in October - December 1971. More ash (after ignition of the organic matter, consisting mostly of these components) per cubic meter of water is found close to land masses (continents and islands) and above shallow submarine elevations than in the open ocean. Preferred biotops of planktonic diatoms in the region described are temperate shallow water and tropical coastal upwelling areas. Radiolarians rarely occur close to the continent, but are abundant in pelagic warm water masses, even near islands. Foraminifers are similar to the radiolarians, rarer in the coastal water mass of the continent than in the open ocean or off oceanic islands. Their abundance is highest outside the upwelling area off NW Africa. Molluscs generally outnumber planktonic foraminifers, implying that the carbonate cycle of the ocean might be influenced considerably by these animals. The molluscs include heteropods, pteropods, and larvae of benthic bivalves and gastropods. Larvae of benthic molluscs occur more frequently close to continental and island margins and above submarine shoals (in this case mostly guyots) than in the open ocean. Their size increases, but they decrease in number with increasing distance from their area of origin. Ostracods and fish have only been found in small numbers concentrated off NW Africa. All of the above-mentioned components occur in higher abundances in the surface water than in subsurface waters. They are closely related to the hydrography of the sampled water masses (here defined through temperature measurements). Relatively warm water masses of the southeastern branches of the Gulf Stream system transport subtropical and southern temperate species to the Bay of Biscay, relatively cool water masses of the Portugal and Canary Currents carry transitional faunal elements along the NW African coast southwards to tropical regions. These mix in the northwest African upwelling area with tropical faunal elements which are generally assumed to live in the subsurface water masses and which probably have been transported northwards to this area by a subsurface counter current. The faunas typical for tropical surface water masses are not only reduced due to the tongue of cool water extending southwards along the coast, but they are also removed from the coastal zone by the upwelling subsurface water masses carrying their own shell and skeleton assemblages. Tropical water masses contain much more shelland skeleton-producing plankters than subtropical and temperate ones. The climatic conditions found at different latitudes control the development and intensity of a separate continental coastal water mass with its own plankton assemblages. Extent of this water mass and steepness of gradients between the pelagic and coastal environment limit the occurrence of pelagic plankton close to the continental coast. A similar water mass in only weakly developed off oceanic islands.