Plant species, biomass and environmental characteristics of relevés along the North America Arctic bioclimate gradient

Question: How do interactions between the physical environment and biotic properties of vegetation influence the formation of small patterned-ground features along the Arctic bioclimate gradient? Location: At 68° to 78°N: six locations along the Dalton Highway in arctic Alaska and three in Canada (Banks Island, Prince Patrick Island and Ellef Ringnes Island). Methods: We analysed floristic and structural vegetation, biomass and abiotic data (soil chemical and physical parameters, the n-factor [a soil thermal index] and spectral information [NDVI, LAI]) on 147 microhabitat releves of zonalpatterned-ground features. Using mapping, table analysis (JUICE) and ordination techniques (NMDS). Results: Table analysis using JUICE and the phi-coefficient to identify diagnostic species revealed clear groups of diagnostic plant taxa in four of the five zonal vegetation complexes. Plant communities and zonal complexes were generally well separated in the NMDS ordination. The Alaska and Canada communities were spatially separated in the ordination because of different glacial histories and location in separate floristic provinces, but there was no single controlling environmental gradient. Vegetation structure, particularly that of bryophytes and total biomass, strongly affected thermal properties of the soils. Patterned-ground complexes with the largest thermal differential between the patterned-ground features and the surrounding vegetation exhibited the clearest patterned-ground morphologies.

Peatland depths and radiocarbon dates, recent decades, North America

Peatland ecosystems store about 500-600 Pg of organic carbon, largely accumulated since the last glaciation. Whether they continue to sequester carbon or release it as greenhouse gases, perhaps in large amounts, is important in Earth's temperature dynamics. Given both ages and depths of numerous dated sample peatlands, their rate of carbon sequestration can be estimated throughout the Holocene. Here we use average values for carbon content per unit volume, the geographical extent of peatlands, and ecological models of peatland establishment and growth, to reconstruct the time-trajectory of peatland carbon sequestration in North America and project it into the future. Peatlands there contain ~163 Pg of carbon. Ignoring effects of climate change and other major anthropogenic disturbances, the rate of carbon accumulation is projected to decline slowly over millennia as reduced net carbon accumulation in existing peatlands is largely balanced by new peatland establishment. Peatlands are one of few long-term terrestrial carbon sinks, probably important for global carbon regulation in future generations. This study contributes to a better understanding of these ecosystems that will assist their inclusion in earth-system models, and therefore their management to maintain carbon storage during climate change.

Monthly air temperatures, snow density and snow and ice thickness at Malcolm Ramsay Lake, North America

In northern regions where observational data is sparse, lake ice models are ideal tools as they can provide valuable information on ice cover regimes. The Canadian Lake Ice Model was used to simulate ice cover for a lake near Churchill, Manitoba, Canada throughout the 2008/2009 and 2009/2010 ice covered seasons. To validate and improve the model results, in situ measurements of the ice cover through both seasons were obtained using an upward-looking sonar device Shallow Water Ice Profiler (SWIP) installed on the bottom of the lake. The SWIP identified the ice-on/off dates as well as collected ice thickness measurements. In addition, a digital camera was installed on shore to capture images of the ice cover through the seasons and field measurements were obtained of snow depth on the ice, and both the thickness of snow ice (if present) and total ice cover. Altering the amounts of snow cover on the ice surface to represent potential snow redistribution affected simulated freeze-up dates by a maximum of 22 days and break-up dates by a maximum of 12 days, highlighting the importance of accurately representing the snowpack for lake ice modelling. The late season ice thickness tended to be under estimated by the simulations with break-up occurring too early, however, the evolution of the ice cover was simulated to fall between the range of the full snow and no snow scenario, with the thickness being dependant on the amount of snow cover on the ice surface.

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).

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.

(Table 1) n-alkane distribution, and stable carbon and deuterium isotope ratios of DSDP Hole 10-94

Hydrology, source region, and timing of precipitation are important controls on the climate of the Great Plains of North America and the composition of terrestrial ecosystems. Moisture delivered to the Great Plains varies seasonally and predominately derives from the Gulf of Mexico/Atlantic Ocean with minor contributions from the Pacific Ocean and Arctic region. For this work, we evaluate long-term relationships for the past ~ 35 million years between North American hydrology, climate, and floral change, using isotopic records and average carbon chain lengths of higher plant n-alkanes from Gulf of Mexico sediments (DSDP Site 94). We find that carbon isotope values (d13C) of n-alkanes, corrected for variations in the d13C value of atmospheric CO2, provide minor evidence for contributions of C4 plants prior to the Middle Miocene. A sharp spike in C4 input is identified during the Middle Miocene Climatic Optimum, and the influence of C4 plants steadily increased during the Late Miocene into the Pleistocene - consistent with other North American records. Chain-length distributions of n-alkanes, indicative of the composition of higher plant communities, remained remarkably constant from 33 to 4 Ma. However, a trend toward longer chain lengths occurred during the past 4 million years, concurrent with an increase in d13C values, indicating increased C4 plant influence and potentially aridity. The hydrogen isotope values (dD) of n-alkanes are relatively invariant between 33 and 9 Ma, and then become substantially more negative (75 per mil) from 9 to 2 Ma. Changes in the plant community and temperature of precipitation can solely account for the observed variations in dD from 33 to 5 Ma, but cannot account for Plio-Pleistocene dD variations and imply substantial changes in the source region of precipitation and seasonality of moisture delivery. We posit that hydrological changes were linked to tectonic and oceanographic processes including the shoaling and closure of the Panamanian Seaway, amplification of North Atlantic Deep Water Production and an associated increase of meridional winds. The southerly movement of the Intertropical Convergence Zone near 4 Ma allowed for the development of a near-modern pressure/storm track system, driving increased aridity and changes in seasonality within the North American interior.

Eh and pH values in sediments and composition of interstitial waters from the Frolikha Bay, Baikal Lake

Redox conditions and compositions of bottom sediments and sedimentary pore waters in the area of the hydrothermal vent in the Frolikha Bay (Baikal Lake) are under discussion. According to obtained results, the submarine vent and its companion spring nearby on the land originate from a common source. The most convincing evidence for their relation comes from proximity of stable oxygen and hydrogen isotope compositions in the pore waters and spring water. The isotope composition indicates meteoric origin of the pore waters, but their major- and minor element compositions have influence of deep water, which may seep through the permeable faulted crust. Although the pore waters near the submarine vent have specific enrichment in major and minor constituents, hydrothermal discharge at the Baikal bottom causes minor influence on water composition of the Baikal Lake, unlike freshwater lakes in rifts of the East Africa and North America.

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.

(Table 3) Air temperatures, and snow and ice characteristics of fresh and brackish water lakes in the Northwest Territories, Canada

The algorithms designed to estimate snow water equivalent (SWE) using passive microwave measurements falter in lake-rich high-latitude environments due to the emission properties of ice covered lakes on low frequency measurements. Microwave emission models have been used to simulate brightness temperatures (Tbs) for snowpack characteristics in terrestrial environments but cannot be applied to snow on lakes because of the differing subsurface emissivities and scattering matrices present in ice. This paper examines the performance of a modified version of the Helsinki University of Technology (HUT) snow emission model that incorporates microwave emission from lake ice and sub-ice water. Inputs to the HUT model include measurements collected over brackish and freshwater lakes north of Inuvik, Northwest Territories, Canada in April 2008, consisting of snowpack (depth, density, and snow water equivalent) and lake ice (thickness and ice type). Coincident airborne radiometer measurements at a resolution of 80x100 m were used as ground-truth to evaluate the simulations. The results indicate that subsurface media are simulated best when utilizing a modeled effective grain size and a 1 mm RMS surface roughness at the ice/water interface compared to using measured grain size and a flat Fresnel reflective surface as input. Simulations at 37 GHz (vertical polarization) produce the best results compared to airborne Tbs, with a Root Mean Square Error (RMSE) of 6.2 K and 7.9 K, as well as Mean Bias Errors (MBEs) of -8.4 K and -8.8 K for brackish and freshwater sites respectively. Freshwater simulations at 6.9 and 19 GHz H exhibited low RMSE (10.53 and 6.15 K respectively) and MBE (-5.37 and 8.36 K respectively) but did not accurately simulate Tb variability (R= -0.15 and 0.01 respectively). Over brackish water, 6.9 GHz simulations had poor agreement with airborne Tbs, while 19 GHz V exhibited a low RMSE (6.15 K), MBE (-4.52 K) and improved relative agreement to airborne measurements (R = 0.47). Salinity considerations reduced 6.9 GHz errors substantially, with a drop in RMSE from 51.48 K and 57.18 K for H and V polarizations respectively, to 26.2 K and 31.6 K, although Tb variability was not well simulated. With best results at 37 GHz, HUT simulations exhibit the potential to track Tb evolution, and therefore SWE through the winter season.

(Table 1) Diagnostic characters of specimens of Marenzelleria bastropi sp. n., M. arctia, M. neglecta, M. viridis and M. wireni

Marenzelleria bastropi, a new species of Spionidae (Polychaeta) from the brackish water Currituck Sound, North Carolina, is described. The new species is characterized by the great number of chaetigers between the first neuro- and notopodial hooded hooks, the extension of the nuchal organ up to the end of chaetiger 2/middle of chaetiger 3 and the presence of about 60-90 branchiate chaetigers. Marenzelleria bastropi sp. nov. is closely related to M. neglecta (Sikorski and Bick, 2004) and Marenzelleria viridis (Verrill, 1873). Marenzelleria wireni Augener, 1913 is described here for the first time from western Spitsbergen. Adult specimens are investigated and compared with specimens from other areas of distribution. A key for subadult and adult specimens of all Marenzelleria species is provided.

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.

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

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