(Appendix) Recent surface areas and historical lengths of glaciers from Patagonia and Tierra del Fuego, South America

Fluctuations in the length of 72 glaciers in the Northern and Southern Patagonia Icefield (NPI and SPI, respectively) and the Cordillera Darwin Icefield (CDI) were estimated between 1945 and 2005. The information obtained from historical maps based on 1945 aerial photographs was compared to ASTER and Landsat satellite images and to information found in the literature. The majority of glaciers have retreated considerably, with maximum values of 12.2 km for Marinelli Glacier in the CDI, 11.6 km for O'Higgins Glacier in the SPI and 5.7 km for San Rafael Glacier in the NPI. Among the 20 glaciers that have retreated the most relative to their size, small (less than 50 km**2) and medium (between 50 and 200 km**2) glaciers are the most affected. However, no direct relation between glacier retreat and size was found for the 72 glaciers studied. The highest percentage retreat in the CDI was by the CDI-03 Glacier (37.9%) and Marinelli Glacier (37.6%). In the SPI, relative retreats were heterogeneous and fluctuated between 27.2% (Amelia Glacier) and 0.4% (Viedma Glacier). In the NPI, relative retreat was very high for Strindberg and Cachet glaciers (35.9% and 27.6%, respectively) but for the remaining glaciers in this icefield it ranged between 11.8% (Piscis Glacier) and 3.6% (San Quintin Glacier). In addition to surface area, the surface slope (calculated on the basis of the DEM SRTM) was also related to the relative retreat and no straightforward relation was found. From a global point of view, we suggest that glacier retreat in the region is controlled firstly by atmospheric warming, as it has been reported in this area. Besides the general increase in temperature observed, no signal of a geographical pattern for the fluctuations in glacier length was found. Consequently, glaciers appear to initially react to local conditions most probably induced by their exposition, geometry and hypsometry. The heterogeneity of rates of retreat suggests that differences in basin geometry, glacier dynamics and response time are key features to explain fluctuations of each glacier.

Clay mineral and stratigraphic record of two sediment cores from the Aegean Sea

Different source areas, oceanography and climate regimes influenced the clay mineral assemblages and grain size distribution of two sediment cores from the North and South Aegean Sea during the last glacial and the Holocene. In the North Aegean Sea, clay mineral composition is mainly controlled by sea level evolution, melting of southeastern European glaciers, and establishment of the connection between the Black Sea and Aegean Sea. The long-term development of clay mineral assemblages in the South Aegean Sea reflects changes in the Nile discharge and African dust input. At this site, the establishment of pluvial conditions in the Nile catchment during the early to middle Holocene resulted in a substantial rise in smectite/illite ratios. In the late Holocene, stepwise aridification of the southern borderlands caused an increase in windblown sediment material and a decrease in Nile suspended material. The clay mineral records exhibit periodic millennial-scale fluctuations. In the North Aegean Sea, the changes are centred at a period of 1.3-1.8 ka and can be attributed to short-term climate and weathering changes in the northern borderlands. The changes in the South Aegean Sea are centred at periods of 3.2-4.3, 1.9-2.4 and 1.3-1.7 ka reflecting short-term changes in wind strength and Northeast African hydrology.

Paleogene planktonic foraminifera from the southern Kerguelen Plateau

An essentially complete Paleogene record was recovered on the Central and Southern Kerguelen plateaus (55°-59°S) in a calcareous biofacies. Recovery deteriorated in the middle Eocene and down to the upper Paleocene because of the presence of interbedded cherts and chalks. The stratigraphic distribution of about 70 taxa of planktonic foraminifers recovered at Sites 747-749 is reported in this paper. Faunas exhibited fairly high diversity (approximately 20-25 species) in the early Eocene, followed by a gradual reduction in diversity in the middle Eocene. A brief incursion of tropical keeled morozovellids occurred near the Paleocene/Eocene boundary, similar to that recorded on the Maud Rise (ODP Sites 689 and 690). The high-latitude Paleogene zonal scheme developed for ODP Leg 113 sites has been adopted (with minor modifications) for the lower Eocene-Oligocene part of the Kerguelen Plateau record. A representative Oligocene (polarity chronozones 7-13) and late Eocene-late middle Eocene (questionably polarity chronozones 16-18) magnetostratigraphic record has allowed the calibration of several biostratigraphic datum levels to the standard Global Polarity Time Scale (GPTS) and established their essential synchrony between low and high latitudes.

Dissolved Si isotope composition measured in water samples during METEOR cruises M77/3 and M77/4

Silicon isotopes are a powerful tool to investigate the cycling of dissolved silicon (Si). In this study the distribution of the Si isotope composition of dissolved silicic acid (d30Si(OH)4) was analyzed in the water column of the Eastern Equatorial Pacific (EEP) where one of the globally largest Oxygen Minimum Zones (OMZs) is located. Samples were collected at 7 stations along two meridional transects from the equator to 14°S at 85°50'W and 82°00'W off the Ecuadorian and Peruvian coast. Surface waters show a large range in isotope compositions d30Si(OH)4 (+2.2 per mil to +4.4 per mil) with the highest values found at the southernmost station at 14°S. This station also revealed the most depleted silicic acid concentrations (0.2 µmol/kg), which is a function of the high degree of Si utilization by diatoms and admixture with waters from highly productive areas. Samples within the upper water column and the OMZ at oxygen concentrations below 10 µmol/kg are characterized by a large range in d30Si(OH)4, which mainly reflects advection and mixing of different water masses, even though the highly dynamic hydrographic system of the upwelling area off Peru does not allow the identification of clear Si isotope signals for distinct water masses. Therefore we cannot rule out that also dissolution processes have an influence on the d30Si(OH)4 signature in the subsurface water column. Deep water masses (>2000 m) in the study area show a mean d30Si(OH)4 of +1.2±0.2 per mil, which is in agreement with previous studies from the eastern and central Pacific. Comparison of the new deep water data of this study and previously published data from the central Pacific and Southern Ocean reveal substantially higher d30Si(OH)4 values than deep water signatures from the North Pacific. As there is no clear correlation between d30Si(OH)4 and silicic acid concentrations in the entire data set the distribution of d30Si(OH)4 signatures in deep waters of the Pacific is considered to be mainly a consequence of the mixing of several end member water masses with distinct Si isotope signatures including Lower Circumpolar Deep Water (LCDW) and North Pacific Deep Water (NPDW).

Planktonic foraminifera and sea surface temperature reconstruction of sediments from the Mediterranean Sea

Water exchange between the Black Sea and the Mediterranean Sea has been a major focus of the paleohydrography of the eastern Mediterranean. Glacial melt water released from the Black Sea is a potential factor in the formation of sapropel S1, an organic-rich sediment layer that accumulated during the Early Holocene. A high-resolution study done on sediments from the Marmara Sea, the gateway between the Mediterranean and the Black Sea, sheds light on the Holocene exchange processes. Past sea surface temperature and sea surface salinity (SSS) were derived from stable oxygen isotope ratios (delta18O) of foraminiferal calcite and alkenone unsaturation ratios (Uk'37). Heavy delta18O values and high SSS in the Marmara Sea suggest absence of low salinity water from the Black Sea during S1. The comparison with data from the Levantine Basin and southern Aegean Sea outlines gradients of freshening in the eastern Mediterranean Sea, whereby the major sources of freshwater were closer to the Levantine Basin. It is thus concluded that the Black Sea was not a major freshwater source contributing to formation of S1. Given the absence of a low salinity layer, the deposition of organic-rich sediments corresponding to S1 in the Marmara Sea is likely the result of the global transgression and the concomitant re-organization of biogeochemical cycles, leading to enhanced productivity as shown by Globigerina bulloides.

Abundance and size composition of ctenophore Mnemiopsis leidyi population in the Caspian Sea on 2001-06-24 to 2001-07-01

Abundance and size distribution of ctenophore Mnemiopsis leidyi in different parts of the Caspian Sea were studied in summer 2001 in relation to environmental conditions. In general, principal differences were found in M. leidyi abundance and population reproduction activity in northern-, middle- and southern Caspian waters. Ctenophore was practically absent in the northern Caspian. In the west of the middle Caspian Sea it penetrated far to the north demonstrating low reproduction activity. In the east the first single comb jellies were pointed out only in the most south of the region. In the warmest and most productive southern part of the Caspian Sea several zones of M. leidyi active breeding were found with total abundance exceeding 6000 #/m**2. Breeding activity and abundance of ctenophores increased here from the east to the west exceeding maximum values along the western coast of the southern Caspian Sea in regions of intensive sprat catching. Dependence of M. leidyi population development on temperature conditions was mentioned. On the base of remote sensed surface temperature, chlorophyll, and suspended mater distribution analysis possible ctenophore settling mechanisms by mesoscale dynamic structures were examined. Practical applications of obtained results are discussed for using effective biological methods to prevent catastrophic consequences of M. leidyi invasion to the Caspian Sea.

Global compilation of last glacial maximum oxygen isotopic ratios for planktonic and benthic foraminifera

We use the fully coupled atmosphere-ocean three-dimensional model of intermediate complexity iLOVECLIM to simulate the climate and oxygen stable isotopic signal during the Last Glacial Maximum (LGM, 21 000 yr). By using a model that is able to explicitly simulate the sensor (d18O), results can be directly compared with data from climatic archives in the different realms. Our results indicate that iLOVECLIM reproduces well the main feature of the LGM climate in the atmospheric and oceanic components. The annual mean d18O in precipitation shows more depleted values in the northern and southern high latitudes during the LGM. The model reproduces very well the spatial gradient observed in ice core records over the Greenland ice-sheet. We observe a general pattern toward more enriched values for continental calcite d18O in the model at the LGM, in agreement with speleothem data. This can be explained by both a general atmospheric cooling in the tropical and subtropical regions and a reduction in precipitation as confirmed by reconstruction derived from pollens and plant macrofossils. Data-model comparison for sea surface temperature indicates that iLOVECLIM is capable to satisfyingly simulate the change in oceanic surface conditions between the LGM and present. Our data-model comparison for calcite d18O allows investigating the large discrepancies with respect to glacial temperatures recorded by different microfossil proxies in the North Atlantic region. The results argue for a trong mean annual cooling between the LGM and present (>6°C), supporting the foraminifera transfer function reconstruction but in disagreement with alkenones and dinocyst reconstructions. The data-model comparison also reveals that large positive calcite d18O anomaly in the Southern Ocean may be explained by an important cooling, although the driver of this pattern is unclear. We deduce a large positive d18Osw anomaly for the north Indian Ocean that contrasts with a large negative d18Osw anomaly in the China Sea between the LGM and present. This pattern may be linked to changes in the hydrological cycle over these regions. Our simulation of the deep ocean suggests that changes in d18Osw between the LGM and present are not spatially homogenous. This is supported by reconstructions derived from pore fluids in deep-sea sediments. The model underestimates the deep ocean cooling thus biasing the comparison with benthic calcite d18O data. Nonetheless, our data-model comparison support a heterogeneous cooling of few degrees (2-4°C) in the LGM Ocean.

Chemical and isotope compositions of sediments from DSDP Hole 86-576A and Core LL44-GPC-3

A record of changes in Pb and Sr isotopic composition of two cores (DSDP 86-576A and LL44- GPC-3) from the red clay region of the central North Pacific has been determined for the past 60-65 million years. The isotope records of the eolian silicate fraction of the red clays reflect the change in source area as the core sites migrated under different wind systems. The Sr isotope compositions of eolian silicate material are consistent with Asian loess and North American arc volcanism that has been recognized from mineralogical studies. The silicate-bound eolian Pb isotopic compositions similarly reflect Asian loess and arc volcanism. The isotope records of three ferromanganese crusts from similar locations in the central Pacific are similar to the eolian component of red clays, but offset to less radiogenic values. This may be due to two mechanisms: (1) Pb that can be removed from eolian material by seawater is much less radiogenic, or less likely (2) hydrothermal Pb can be transported further away from venting sites through particle exchange with seawater, despite hydrothermal venting acting as a net sink of oceanic Pb. The temporal changes in Pb isotopes in the ferromanganese crusts, bulk red clays and eolian silicates are similar although offset from each other suggesting that eolian deposition is an important source of Pb to seawater and to ferromanganese crusts. This contrasts with the Atlantic and Southern Ocean where more intense deep water flow leads to isotopic gradients in FeMn crusts that do not reflect surface water conditions immediately above the crust. A mechanism is proposed which accounts for Pacific deepwater Pb being isotopically influenced by eolian deposition.

Oxygen and hydrogen isotopes measured in the Northeast Atlantic

Only a few studies have examined the variation of oxygen and hydrogen isotopes of seawater in NE Atlantic water masses, and data are especially sparse for intermediate and deep-water masses. The current study greatly expands this record with 527 d18O values from 47 stations located throughout the mid- to low-latitude NE Atlantic. In addition, dD was analyzed in the 192 samples collected along the GEOTRACES North Atlantic Transect GA03 (GA03_e=KN199-4) and the 115 Iberia-Forams cruise samples from the western and southern Iberian margin. An intercomparison study between the two stable isotope measurement techniques (cavity ring-down laser spectroscopy and magnetic-sector isotope ratio mass spectrometry) used to analyze GA03_e samples reveals relatively good agreement for both hydrogen and oxygen isotope ratios. The surface (0-100 m) and central (100-500 m) water isotope data show the typical, evaporation related trend of increasing values equatorward with the exception for the zonal transect off Cape Blanc, NW Africa. Off Cape Blanc, surface water isotope signatures are modified by the upwelling of fresher Antarctic Intermediate Water (AAIW) that generally has isotopic values of 0.0 to 0.5 per mil for d18O and 0 to 2 per mil for dD. Along the Iberian margin the Mediterranean Outflow Water (MOW) is clearly distinguished by its high d18O (0.5-1.1 per mil) and dD (3-6 per mil) values that can be traced into the open Atlantic. Isotopic values in the NE Atlantic Deep Water (NEADW) are relatively low (d18O: -0.1 to 0.5 per mil; dD: -1 to 4 per mil) and show a broader range than observed previously in the northern and southern convection areas. The NEADW is best observed at GA03_e Stations 5 and 7 in the central NE Atlantic basin. Antarctic Bottom Water isotope values are relatively high indicating modification of the original Antarctic source water along the flow path. The reconstructed d18O-salinity relationship for the complete data set has a slope of 0.51, i.e., slightly steeper than the 0.46 described previously by Pierre et al. (1994, J. Mar. Syst. 5 (2), 159-170.) for the tropical to subtropical Northeast Atlantic. This slope decreases to 0.46 for the subtropical North Atlantic Central Water (NACW) and the MOW and to 0.32 for the surface waters of the upper 50 m. The dD-salinity mixing lines have estimated slopes of 3.01 for the complete data, 1.26 for the MOW, 3.47 for the NACW, and 2.63 for the surface waters. The slopes of the d18O-dD relationship are significantly lower than the one for the Global Meteoric Water Line with 5.6 for the complete data set, 2.30 for the MOW, 4.79 for the NACW, and 3.99 for the surface waters. The lower slopes in all the relationships clearly reflect the impact of the evaporation surplus in the subtropics.

Chemical composition of DSDP and ODP cherts

Rare earth element (REE), major, and trace element abundances and relative fractionations in forty nodular cherts sampled by the Deep Sea Drilling Project (DSDP) and Ocean Drilling Program (ODP) indicate that the REE composition of chert records the interplay between terrigenous sources and scavenging from the local seawater. Major and (non-REE) trace element ratios indicate that the aluminosilicate fraction within the chert is similar to NASC (North American Shale Composite), with average Pacific chert including ~7% NASC-like particles, Indian chert ~11% NASC, Atlantic chert ~17% NASC, and southern high latitude (SHL) chert 53% NASC. Using La as a proxy for sum REE, approximations of excessive La (the amount of La in excess of that supplied by the detrital aluminosilicate fraction) indicate that Pacific chert contains the greatest excessive La (85% of total La) and SHL chert the least (38% of total La). As shown by interelement associations, this excessive La is most likely an adsorbed component onto aluminosilicate and phosphatic phases. Accordingly, chert from the large Pacific Ocean, where deposition occurs relatively removed from significant terrigenous input, records a depositional REE signal dominated by adsorption of dissolved REEs from seawater. Pacific chert Ce/Ce* <<1 and normative La/Yb ~ 0.8-1, resulting from adsorption of local Ce-depleted seawater and preferential adsorption of LREEs from seawater (e.g., normative La/Yb ~0.4), which increases the normative La/Yb ratio recorded in chert. Chert from the Atlantic basin, a moderately sized ocean basin lined by passive margins and with more terrigenous input than the Pacific, records a mix of adsorptive and terrigenous REE signals, with moderately negative Ce anomalies and normative La/Yb ratios intermediate to those of the Pacific and those of terrigenous input. Chert from the SHL region is dominated by the large terrigenous input on the Antarctic passive margin, with inherited Ce/Ce* ~1 and inherited normative La/Yb values of ~1.2-1.4. Ce/Ce* does not vary with age, either throughout the entire data base or within a particular basin. Overall, Ce/Ce* does not correlate with P2O5 concentrations, even though phosphatic phases may be an important REE carrier.

Age determination and Argon geochemistry of Hawaiian Emperor seamounts

The Hawaiian-Emperor bend has played a prominent yet controversial role in deciphering past Pacific plate motions and the tempo of plate motion change. New ages for volcanoes of the central and southern Emperor chain define large changes in volcanic migration rate with little associated change in the chain's trend, which suggests that the bend did not form by slowing of the Hawaiian hot spot. Initiation of the bend near Kimmei seamount about 50 million years ago (MA) was coincident with realignment of Pacific spreading centers and early magmatism in western Pacific arcs, consistent with formation of the bend by changed Pacific plate motion.

Petrography and geochemical composition of the Atlantis II Fracture Zone

SeaBeam echo sounding, seismic reflection, magnetics, and gravity profiles were run along closely spaced tracks (5 km) parallel to the Atlantis II Fracture Zone on the Southwest Indian Ridge, giving 80% bathymetric coverage of a 30- * 170-nmi strip centered over the fracture zone. The southern and northern rift valleys of the ridge were clearly defined and offset north-south by 199 km. The rift valleys are typical of those found elsewhere on the Southwest Indian Ridge, with relief of more than 2200 m and widths from 22 to 38 km. The ridge-transform intersections are marked by deep nodal basins lying on the transform side of the neovolcanic zone that defines the present-day spreading axis. The walls of the transform generally are steep (25°-40°), although locally, they can be more subdued. The deepest point in the transform is 6480 m in the southern nodal basin, and the shallowest is an uplifted wave-cut terrace that exposes plutonic rocks from the deepest layer of the ocean crust at 700 m. The transform valley is bisected by a 1.5-km-high median tectonic ridge that extends from the northern ridge-transform intersection to the midpoint of the active transform. The seismic survey showed that the floor of the transform contains up to 0.5 km of sediment. Piston-coring at two locations on the transform floor recovered more than 1 m of sand and gravel, which appears to be turbidites shed from the walls of the fracture zone. Extensive dredging showed that more than two-thirds of the crust exposed in the transform valley and its walls were plutonic rocks, principally gabbros and residual mantle peridotites. In contrast, based on dredging and seafloor morphology, only relatively undisrupted pillow basalt flows have been exposed on crust of the same age spreading away from the transform. Magnetic anomalies are well defined out to 11 m.y. over the flanking transverse ridges and transform valley, even where layer 2 appears to be absent. The total opening rate is 1.6 cm/yr, but the arrangement of the anomalies indicates that the spreading for each ridge is asymmetric, with the ridge flanks facing the transform spreading at a rate of 1.0 cm/yr. Such an asymmetric spreading pattern requires that both the northern and southern ridges migrate away from each other at 0.2 cm/yr, thus lengthening the transform at 0.4 cm/yr for the last 11 m.y. To the north, the fracture zone valley is oriented differently from the present-day transform, indicating a paleospreading direction change at 17 m.y. from N10°E to due north-south. This change placed the transform into extension for the 11-m.y. period required for simple orthogonal ridge-transform geometry to be reestablished and produced a large transtensional basin within the transform valley. This basin was split by continued transform slip after 11 m.y., with the larger half moving to the north with the African Plate.