Cenozoic radiolaria in the western North Atlantic

Eight Cenozoic radiolarian zones were recognized in samples from two holes at Site 603, drilled on the lower continental rise off North America during Leg 93 of the Deep Sea Drilling Project. Paleocene to early Eocene radiolarian zones (Bekoma bidartensis, Buryella clinata, and Phormocyrtis striata striata zones) and early to late Miocene radiolarian zones (Calocycletta costata, Dorcadospyris alata, Diartus petterssoni, and Didymocyrtis antepenultima zones) were recognized in sediments from Holes 603 and 603B. In addition, a new Paleocene Bekoma campechensis radiolarian Zone is defined by the interval between the first morphotypic appearance of B. campechensis and the B. campechensis-B. bidartensis evolutionary transition. This zone is immediately below the B. bidartensis Zone of Foreman (1973), and has previously been discussed as a Paleocene "unnamed zone" by other investigators. A hiatus between Neogene and Paleogene sequences was also recognized in the radiolarian faunas.

Lead isotopes of feldspar in North Atlantic Ocean sediments

The provenance of ice-rafted debris (IRD) deposited in the North Atlantic before, during, and after Heinrich event 2 has been determined through measuring the lead isotopic composition of single feldspar grains and multiple-grain composites from the larger than 150-µm size fraction, from cores from the eastern and western North Atlantic and from the Labrador Sea. Single-grain analyses are used to identify the specific continental sources of the IRD, whereas composite samples are used to assess the relative IRD contributions from different sources. All single grains from Heinrich layer 2 (H 2) as well as H 2 composites plot along a correlation line on a 207Pb/204Pb versus 206Pb/204Pb diagram characteristic of the Churchill province of the Canadian shield. This is yet another strong piece of evidence that this Heinrich event was dominated by a massive iceberg discharge of the Laurentide ice sheet lobe located over Hudson Bay. In contrast, single grains from the ambient glacial sediment (above and below H 2) have multiple sources: many of them also lie along the correlation line with H 2 grains, but many others have Pb signatures consistent with derivation from the Grenville province and the Appalachian range in North America and possibly from Scandinavia and Greenland. Composites from the ambient sediment generally lie well to the right of the H 2 reference line in agreement with the results of the single-grain analyses. The evidence provided by lead isotopes regarding the dominant role played by the Hudson Bay lobe of the Laurentide ice sheet in the development of the Heinrich events lends support to the binge/purge model advanced by MacAyeal [1993a, b] that invokes trapping of geothermal heat by the base of the icecap and subsequent basal melting as the mechanism that triggered the Heinrich events.

Parasound profiles and multibeam surveys in the Arctic Ocean

During the Pleistocene glaciations, Arctic ice sheets on western Eurasia, Greenland and North America terminated at their continental margins. In contrast, the exposed continental shelves in the Beringian region of Siberia are thought to have been covered by a tundra landscape. Evidence of grounded ice on seafloor ridges and plateaux off the coast of the Beringian margin, at depths of up to 1,000 m, have generally been attributed to ice shelves or giant icebergs that spread oceanwards during glacial maxima. Here we identify marine glaciogenic landforms visible in seismic profiles and detailed bathymetric maps along the East Siberian continental margin. We interpret these features, which occur in present water depths of up to 1,200 m, as traces from grounding events of ice sheets and ice shelves. We conclude that the Siberian Shelf edge and parts of the Arctic Ocean were covered by ice sheets of about 1 km in thickness during several Pleistocene glaciations before the most recent glacial period, which must have had a significant influence on albedo and oceanic and atmospheric circulation.

Pb isotopes in individual sand-sized (>150 µm) ice-rafted feldspars deposited at subpolar North Atlantic Ocean Deep Sea Drilling Project Site 611 during the late Pliocene (Marine Isotope Stages, MIS, G6 through 100)

The onset of abundant ice-rafted debris (IRD) deposition in the Nordic Seas and subpolar North Atlantic Ocean 2.72 millions of years ago (Ma) is thought to record the Pliocene onset of major northern hemisphere glaciation (NHG) due to a synchronous advance of North American Laurentide, Scandinavian and Greenland ice-sheets to their marine calving margins during marine isotope stage (MIS) G6. Numerous marine and terrestrial records from the Nordic Seas region indicate that extensive ice sheets on Greenland and Scandinavia increased IRD inputs to these seas from 2.72 Ma. The timing of ice-sheet expansion on North America as tracked by IRD deposition in the subpolar North Atlantic Ocean, however, is less clear because both Europe and North America are potential sources for icebergs in this region. Moreover, cosmogenic-dating of terrestrial tills on North America indicate that the Laurentide Ice Sheet did not extend to ~39°N until 2.4 ±0.14 Ma, at least 180 ka after the onset of major IRD deposition at 2.72 Ma. To address this problem,we present the first detailed analysis of the geochemical provenance of individual sand-sized IRD deposited in the subpolar North Atlantic Ocean between MIS G6 and 100 (~2.72-2.52 Ma). IRD provenance is assessed using laser ablation lead (Pb) isotope analyses of single ice-rafted (>150 mm) feldspar grains. To track when an ice-rafting setting consistent with major NHG first occurred in the North Atlantic Ocean during the Pliocene intensification of NHG (iNHG), we investigate when the Pb-isotope composition (206Pb/204Pb, 207Pb/204Pb, 208Pb/204Pb) of feldspars deposited at DSDP Site 611 first resembles that determined for IRD deposited at this site during MIS 100, the oldest glacial for which there exists convincing evidence for widespread glaciation of North America. Whilst Quaternary-magnitude IRD fluxes exist at Site 611 during glacials from 2.72 Ma, we find that the provenance of this IRD is not constant. Instead, we find that the Pb isotope composition of IRD at our study site is not consistent with major NHG until MIS G2 (2.64 Ma). We hypothesise that IRD deposition in the North Atlantic Ocean prior to MIS G2 was dominated by iceberg calving from Greenland and Scandinavia. We further suggest that the grounding line of continental ice on Northeast America may not have extended onto the continental shelf and calved significant numbers of icebergs to the North Atlantic Ocean during glacials until 2.64 Ma.

Database of Ice-Rich Yedoma Permafrost (IRYP)

Vast portions of Arctic and sub-Arctic Siberia, Alaska and the Yukon Territory are covered by ice-rich silty to sandy deposits that are containing large ice wedges, resulting from syngenetic sedimentation and freezing. Accompanied by wedge-ice growth in polygonal landscapes, the sedimentation process was driven by cold continental climatic and environmental conditions in unglaciated regions during the late Pleistocene, inducing the accumulation of the unique Yedoma deposits up to >50 meters thick. Because of fast incorporation of organic material into syngenetic permafrost during its formation, Yedoma deposits include well-preserved organic matter. Ice-rich deposits like Yedoma are especially prone to degradation triggered by climate changes or human activity. When Yedoma deposits degrade, large amounts of sequestered organic carbon as well as other nutrients are released and become part of active biogeochemical cycling. This could be of global significance for future climate warming as increased permafrost thaw is likely to lead to a positive feedback through enhanced greenhouse gas fluxes. Therefore, a detailed assessment of the current Yedoma deposit coverage and its volume is of importance to estimate its potential response to future climate changes. We synthesized the map of the coverage and thickness estimation, which will provide critical data needed for further research. In particular, this preliminary Yedoma map is a great step forward to understand the spatial heterogeneity of Yedoma deposits and its regional coverage. There will be further applications in the context of reconstructing paleo-environmental dynamics and past ecosystems like the mammoth-steppe-tundra, or ground ice distribution including future thermokarst vulnerability. Moreover, the map will be a crucial improvement of the data basis needed to refine the present-day Yedoma permafrost organic carbon inventory, which is assumed to be between 83±12 (Strauss et al., 2013, doi:10.1002/2013GL058088) and 129±30 (Walter Anthony et al., 2014, doi:10.1038/nature13560) gigatonnes (Gt) of organic carbon in perennially-frozen archives. Hence, here we synthesize data on the circum-Arctic and sub-Arctic distribution and thickness of Yedoma for compiling a preliminary circum-polar Yedoma map. For compiling this map, we used (1) maps of the previous Yedoma coverage estimates, (2) included the digitized areas from Grosse et al. (2013) as well as extracted areas of potential Yedoma distribution from additional surface geological and Quaternary geological maps (1.: 1:500,000: Q-51-V,G; P-51-A,B; P-52-A,B; Q-52-V,G; P-52-V,G; Q-51-A,B; R-51-V,G; R-52-V,G; R-52-A,B; 2.: 1:1,000,000: P-50-51; P-52-53; P-58-59; Q-42-43; Q-44-45; Q-50-51; Q-52-53; Q-54-55; Q-56-57; Q-58-59; Q-60-1; R-(40)-42; R-43-(45); R-(45)-47; R-48-(50); R-51; R-53-(55); R-(55)-57; R-58-(60); S-44-46; S-47-49; S-50-52; S-53-55; 3.: 1:2,500,000: Quaternary map of the territory of Russian Federation, 4.: Alaska Permafrost Map). The digitalization was done using GIS techniques (ArcGIS) and vectorization of raster Images (Adobe Photoshop and Illustrator). Data on Yedoma thickness are obtained from boreholes and exposures reported in the scientific literature. The map and database are still preliminary and will have to undergo a technical and scientific vetting and review process. In their current form, we included a range of attributes for Yedoma area polygons based on lithological and stratigraphical information from the original source maps as well as a confidence level for our classification of an area as Yedoma (3 stages: confirmed, likely, or uncertain). In its current version, our database includes more than 365 boreholes and exposures and more than 2000 digitized Yedoma areas. We expect that the database will continue to grow. In this preliminary stage, we estimate the Northern Hemisphere Yedoma deposit area to cover approximately 625,000 km². We estimate that 53% of the total Yedoma area today is located in the tundra zone, 47% in the taiga zone. Separated from west to east, 29% of the Yedoma area is found in North America and 71 % in North Asia. The latter include 9% in West Siberia, 11% in Central Siberia, 44% in East Siberia and 7% in Far East Russia. Adding the recent maximum Yedoma region (including all Yedoma uplands, thermokarst lakes and basins, and river valleys) of 1.4 million km² (Strauss et al., 2013, doi:10.1002/2013GL058088) and postulating that Yedoma occupied up to 80% of the adjacent formerly exposed and now flooded Beringia shelves (1.9 million km², down to 125 m below modern sea level, between 105°E - 128°W and >68°N), we assume that the Last Glacial Maximum Yedoma region likely covered more than 3 million km² of Beringia. Acknowledgements: This project is part of the Action Group "The Yedoma Region: A Synthesis of Circum-Arctic Distribution and Thickness" (funded by the International Permafrost Association (IPA) to J. Strauss) and is embedded into the Permafrost Carbon Network (working group Yedoma Carbon Stocks). We acknowledge the support by the European Research Council (Starting Grant #338335), the German Federal Ministry of Education and Research (Grant 01DM12011 and "CarboPerm" (03G0836A)), the Initiative and Networking Fund of the Helmholtz Association (#ERC-0013) and the German Federal Environment Agency (UBA, project UFOPLAN FKZ 3712 41 106).

Estimates of %NCW bathing deep North Atlantic (IODP Site 306-U1313) over the past ~3.3 Ma

The circulation and internal structure of the oceans exert a strong influence on Earth's climate because they control latitudinal heat transport and the segregation of carbon between the atmosphere and the abyss (Sigman et al., 2010, doi:10.1038/nature09149). Circulation change, particularly in the Atlantic Ocean, is widely suggested (Bartoli et al., 2005, doi:10.1016/j.epsl.2005.06.020; Haug and Tiedemann, 1998, doi:10.1038/31447; Woodard et al., 2014, doi:10.1126/science.1255586; McKay et al., 2012, doi:10.1073/pnas.1112248109) to have been instrumental in the intensification of Northern Hemisphere glaciation when large ice sheets first developed on North America and Eurasia during the late Pliocene, approximately 2.7 million years ago (Bailey et al., 2013, doi:10.1016/j.quascirev.2013.06.004). Yet the mechanistic link and cause/effect relationship between ocean circulation and glaciation are debated. Here we present new records of North Atlantic Ocean structure using the carbon and neodymium isotopic composition of marine sediments recording deep water for both the Last Glacial to Holocene (35-5 thousand years ago) and the late Pliocene to earliest Pleistocene (3.3-2.4 million years ago). Our data show no secular change. Instead we document major southern-sourced water incursions into the deep North Atlantic during prominent glacials from 2.7 million years ago. Our results suggest that Atlantic circulation acts as a positive feedback rather than as an underlying cause of late Pliocene Northern Hemisphere glaciation. We propose that, once surface Southern Ocean stratification (Sigman, et al., 2004, doi:10.1038/nature02357) and/or extensive sea-ice cover (McKay et al., 2012, doi:10.1073/pnas.1112248109) was established, cold-stage expansions of southern-sourced water such as those documented here enhanced carbon dioxide storage in the deep ocean, helping to increase the amplitude of glacial cycles.

Non-indigenous species in the North and Baltic Seas and the Great Lakes-St. Lawrence River region

List of non-indigenous species (NIS) established in the Great Lakes-St. Lawrence River region and the North and Baltic Seas region, their geographic origin, and taxonomic assignment. Asterisks mark the NIS that occur in both the North and Baltic Seas and the Great Lakes-St. Lawrence River regions. GL, SL, NW, NE, SW and SE denote the Great Lakes, St. Lawrence River, north-west, north-east, south-west, and south-east, respectively. Eurasia represents inland freshwaters except Yangtze River, Indo-Pacific represents Indian Ocean and the archipelago of Indonesia, Malaysia, and Pilipinas, North America (N America) represents inland freshwaters except the Laurentian Great Lakes, St. Lawrence and Mississippi Rivers, while Australia, New Zealand, Africa and South America (S America) cover all inland freshwaters in these areas.

Carbonate sedimentation beneath the Carlifornia Current system

A north-south transect of 17 cores was constructed along the eastern boundary of the California Current system from 33° to 42°N to investigate the changes in biogenic sedimentation over the past 30 kyr. Percentages and mass accumulation rates of CaCO3, Corg, and biogenic opal were assembled at 500 to 1000 years/sample to provide relatively high resolution. Time-space maps reveal a complex pattern of changes that do not follow a simple glacial-interglacial two-mode model. Biogenic sedimentation shows responses that are sometimes time-transgressive and sometimes coeval, and most of the responses show more consistency within a limited geographic area than any temporal consistency. Reconstructed conditions during late oxygen isotope stage 3 were more like early Holocene conditions than any other time during the last 30 kyr. Coastal upwelling and productivity during oxygen isotope stage 3 were relatively strong along the central California margin but were weak along the northern California margin. Precipitation increased during the last glacial interval in the central California region, and the waters of the southern California margin had relatively low productivity. Productivity on the southern Oregon margin was relatively low at the beginning of the last glacial interval, but by about 20 ka, productivity in this area significantly increased. This change suggests that the center of the divergence of the West Wind Drift shifted south at this time. The end of the last glacial interval was characterized by increased productivity in the southern California margin and increased upwelling along the central California margin but upwelling remained weak along the northern California margin. A sudden (<300 years) decrease in CaCO3, Corg, and biogenic opal occurred at 13 ka. The changes suggest a major reorientation of the atmospheric circulation in the North Pacific and western North America and the establishment of a strong seasonality in the central California region. A carbonate preservation event occurred at 10 ka that appears to reflect the uptake of CO2 by the terrestrial biosphere as the northern latitudes were reforested following retreat of the glaciers. The Holocene has been a period of relatively high productivity in the southern California margin, relatively strong coastal upwelling along the central California margin, relatively weak upwelling along the northern California margin, and the northward migration of the divergence zone of the West Wind Drift.

Lithic grain and mineral composition of sand and sandstones from the North Pacific Ocean and the Bering Sea (Table 3)

Sand and sandstone compositions from different types of basins reflect provenance terranes governed by plate tectonics. One hundred and one thin sections of Upper Miocene to Holocene sand-sized material were examined from DSDP/IPOD Sites in the North Pacific Ocean and the Bering Sea. The Gazzi-Dickinson point-counting method was used to establish compositional characteristics of sands from different tectonic settings. Continental margin forearc sands from the western North America continental margin arc system are clearly different from backarc/marginal-sea sands from the Aleutian intraoceanic arc system. The forearc sands have average QFL percentages of 29-42-29, LmLvLst percentages of 32-34-34, 3 Fmwk%M and 0.82 P/F. Aleutian backarc sands have average QFL percentages of 8-22-69. LmLvLst percentages of 9-85-6, 0.5 Fmwk%M and 0.96 P/F. A trend of increasing QFL%Q and decreasing LmLvLst%Lv westward in the backarc region of the Aleutian Ridge reflects the influence of the Asiatic continental margin. Aleutian backarc sands without continental influence have average QFL percentages of 1-20-79, LmLvLst percentages of 1-98-1, 0 Fmwk%M and 0.99 P/F. Of the continental margin forearc samples, sands on the Astoria Fan (west of the Oregon-Washington trench) contain the highest LmLvLst%Lv and lowest P/F; sands from mixed transform-fault and trench settings (Delgada Fan and Gulf of Alaska samples) have slightly higher Qp/Q (0.03); and sands from the Pacific-Juan de Fuca-North America triple junction have the highest Fmwk%M. Delgada Fan and Gulf of Alaska sands have average QFL percentages of 27-38-35, LmLvLst percentages of 37-26-37, 2 Fmwk%M and 0.86 P/F. Astoria Fan sands have average QFL percentages of 35-41-24, LmLvLst percentages of 30-47-23, 3 Fmwk%M and 0.74 P/F. The triple-junction sands have average QFL percentages of 28-59-13, LmLvLst percentages of 25-26-49, 9 Fmwk%M and 0.87 P/F. The petrologic data from the modern ocean basins examined in this study can provide useful analogs for interpretation of ancient oceanic sequences. Our data suggest some refinements of, but generally substantiate, existing petrologic models relating sandstone composition to tectonic setting.