Sulfur-isotope analyses of gypsium samples from the Southeast Greenland margin, southern Alaska and Gulf of Alaska (Table 1)

Gypsum grains were identified in Miocene-Pleistocene sediment cores from two deep-water ODP sites, Site 918 off the SE Greenland margin and Site 887 in the Gulf of Alaska, and in Holocene sediment cores from shallow-water localities in Disenchantment Bay and Muir Inlet in southern Alaska. Although initial morphologic and textural observations suggested a complex system in which the gypsum may have had more than one origin, quantitative sulfur isotope analyses of the gypsum provide evidence of its detrital nature. d34S values in gypsum from southern Alaska range between +0.0 and +7.1 per mil. Gypsum has d34S values between -27.1 and -27.5 per mil in the Gulf of Alaska and values between -28.5 and +0.2 per mil off the SE Greenland margin. All of these isotopic signatures are too highly depleted in d34S to have precipitated from seawater, present or past. In addition there is no significant change in d34S values for gypsum crystals with differing physical characteristics (abraded vs. unabraded) from the same stratigraphic horizon, suggesting all the gypsum is detrital regardless of the degree of abrasion. The isotopic and physical evidence, in combination with the onshore geology the environmental setting, and site characteristics of the gypsum-bearing marine localities, lead us to propose that the ultimate source of the gypsum is precipitation from freeze-induced terrestrial sediment or soil brines. Furthermore the combined evidence suggests that the subsequent occurrence of gypsum in glacimarine sediments results from ice-rafting (by icebergs or sea ice) of the frozen regolith and/or, in the proximal glacimarine setting of southern Alaska, very rapid burial via turbidity currents.

Geochemistry, isotopic ratios and aldehyde concentrations of frost flowers, snow and ice near Barrow, Alaska

Frost flowers, intricate featherlike crystals that grow on refreezing sea ice leads, have been implicated in lower atmospheric chemical reactions. Few studies have presented chemical composition information for frost flowers over time and many of the chemical species commonly associated with Polar tropospheric reactions have never been reported for frost flowers. We undertook this study on the sea ice north of Barrow, Alaska to quantify the major ion, stable oxygen and hydrogen isotope, alkalinity, light absorbance by soluble species, organochlorine, and aldehyde composition of seawater, brine, and frost flowers. For many of these chemical species we present the first measurements from brine or frost flowers. Results show that major ion and alkalinity concentrations, stable isotope values, and major chromophore (NO3- and H2O2) concentrations are controlled by fractionation from seawater and brine. The presence of these chemical species in present and future sea ice scenarios is somewhat predictable. However, aldehydes, organochlorine compounds, light absorbing species, and mercury (part 2 of this research and Sherman et al. (2012, doi:10.1029/2011JD016186)) are deposited to frost flowers through less predictable processes that probably involve the atmosphere as a source. The present and future concentrations of these constituents in frost flowers may not be easily incorporated into future sea ice or lower atmospheric chemistry scenarios. Thinning of Arctic sea ice will likely present more open sea ice leads where young ice, brine, and frost flowers form. How these changing ice conditions will affect the interactions between ice, brine, frost flowers and the lower atmosphere is unknown.

(Table S1) Age determination of sediment cores of the Gulf of Alaska

Oxygen isotope data from planktonic and benthic foraminifera, on a high-resolution age model (44 14C dates spanning 17,400 years), document deglacial environmental change on the southeast Alaska margin (59°33.32'N, 144°9.21'W, 682 m water depth). Surface freshening (i.e., d18O reduction of 0.8 per mil) began at 16,650 ± 170 cal years B.P. during an interval of ice proximal sedimentation, likely due to freshwater input from melting glaciers. A sharp transition to laminated hemipelagic sediments constrains retreat of regional outlet glaciers onto land circa 14,790 ± 380 cal years B.P. Abrupt warming and/or freshening of the surface ocean (i.e., additional d18O reduction of 0.9 per mil) coincides with the Bølling Interstade of northern Europe and Greenland. Cooling and/or higher salinities returned during the Allerød interval, coincident with the Antarctic Cold Reversal, and continue until 11,740 ± 200 cal years B.P., when onset of warming coincides with the end of the Younger Dryas. An abrupt 1 per mil reduction in benthic d18O at 14,250 ± 290 cal years B.P. likely reflects a decrease in bottom water salinity driven by deep mixing of glacial meltwater, a regional megaflood event, or brine formation associated with sea ice. Two laminated opal-rich intervals record discrete episodes of high productivity during the last deglaciation. These events, precisely dated here at 14,790 ± 380 to 12,990 ± 190 cal years B.P. and 11,160 ± 130 to 10,750 ± 220 cal years B.P., likely correlate to similar features observed elsewhere on the margins of the North Pacific and are coeval with episodes of rapid sea level rise. Remobilization of iron from newly inundated continental shelves may have helped to fuel these episodes of elevated primary productivity and sedimentary anoxia.

High-resolution dynamical downscaling products for the North Slope of Alaska and surrounding areas, links to model result files

Climatic changes are most pronounced in northern high latitude regions. Yet, there is a paucity of observational data, both spatially and temporally, such that regional-scale dynamics are not fully captured, limiting our ability to make reliable projections. In this study, a group of dynamical downscaling products were created for the period 1950 to 2100 to better understand climate change and its impacts on hydrology, permafrost, and ecosystems at a resolution suitable for northern Alaska. An ERA-interim reanalysis dataset and the Community Earth System Model (CESM) served as the forcing mechanisms in this dynamical downscaling framework, and the Weather Research & Forecast (WRF) model, embedded with an optimization for the Arctic (Polar WRF), served as the Regional Climate Model (RCM). This downscaled output consists of multiple climatic variables (precipitation, temperature, wind speed, dew point temperature, and surface air pressure) for a 10 km grid spacing at three-hour intervals. The modeling products were evaluated and calibrated using a bias-correction approach. The ERA-interim forced WRF (ERA-WRF) produced reasonable climatic variables as a result, yielding a more closely correlated temperature field than precipitation field when long-term monthly climatology was compared with its forcing and observational data. A linear scaling method then further corrected the bias, based on ERA-interim monthly climatology, and bias-corrected ERA-WRF fields were applied as a reference for calibration of both the historical and the projected CESM forced WRF (CESM-WRF) products. Biases, such as, a cold temperature bias during summer and a warm temperature bias during winter as well as a wet bias for annual precipitation that CESM holds over northern Alaska persisted in CESM-WRF runs. The linear scaling of CESM-WRF eventually produced high-resolution downscaling products for the Alaskan North Slope for hydrological and ecological research, together with the calibrated ERA-WRF run, and its capability extends far beyond that. Other climatic research has been proposed, including exploration of historical and projected climatic extreme events and their possible connections to low-frequency sea-atmospheric oscillations, as well as near-surface permafrost degradation and ice regime shifts of lakes. These dynamically downscaled, bias corrected climatic datasets provide improved spatial and temporal resolution data necessary for ongoing modeling efforts in northern Alaska focused on reconstructing and projecting hydrologic changes, ecosystem processes and responses, and permafrost thermal regimes. The dynamical downscaling methods presented in this study can also be used to create more suitable model input datasets for other sub-regions of the Arctic.