Results of calculated surface air temperature anomaly of COSMOS Pliocene experiment 2 in NetCDF format

In this manuscript we describe the experimental procedure employed at the Alfred Wegener Institute in Germany in the preparation of the simulations for the Pliocene Model Intercomparison Project (PlioMIP). We present a description of the utilized Community Earth System Models (COSMOS, version: COSMOS-landveg r2413, 2009) and document the procedures that we applied to transfer the Pliocene Research, Interpretation and Synoptic Mapping (PRISM) Project mid-Pliocene reconstruction into model forcing fields. The model setup and spin-up procedure are described for both the paleo- and preindustrial (PI) time slices of PlioMIP experiments 1 and 2, and general results that depict the performance of our model setup for mid-Pliocene conditions are presented. The mid-Pliocene, as simulated with our COSMOS setup and PRISM boundary conditions, is both warmer and wetter in the global mean than the PI. The globally averaged annual mean surface air temperature in the mid-Pliocene standalone atmosphere (fully coupled atmosphere-ocean) simulation is 17.35 °C (17.82 °C), which implies a warming of 2.23 °C (3.40 °C) relative to the respective PI control simulation.

Continuously measured 222Rn activity concentrations (daily means) in near surface air at Neumayer Station in the period 1995 through 2011

We report on continuously measured 222Rn activity concentrations in near-surface air at Neumayer Station in the period 1995-2011. This 17-year record showed no long-term trend and has overall mean ± standard deviation of (0.019 ± 0.012) Bq/m**3. A distinct and persistent seasonality could be distinguished with maximum values of (0.028 ± 0.013) Bq/m**3 from January to March and minimum values of (0.015 ± 0.009) Bq/m**3 from May to October. Elevated 222Rn activity concentrations were typically associated with air mass transport from the Antarctic Plateau. Our results do not support a relation between enhanced 222Rn activity concentrations at Neumayer and cyclonic activity or long-range transport from South America. The impact of oceanic 222Rn emissions could not be properly assessed but we tentatively identified regional sea ice extent (SIE) variability as a significant driver of the annual 222Rn cycle.

Underway measurements of halocarbons (air) during POSEIDON cruise POS399 in June 2010

During the DRIVE (Diurnal and Regional Variability of Halogen Emissions) ship campaign we investigated the variability of the halogenated very short-lived substances (VSLS) bromoform (CHBr3), dibromomethane (CH2Br2) and methyl iodide (CH3I) in the marine atmospheric boundary layer in the eastern tropical and subtropical North Atlantic Ocean during May/June 2010. The highest VSLS mixing ratios were found near the Mauritanian coast and close to Lisbon (Portugal). With backward trajectories we identified predominantly air masses from the open North Atlantic with some coastal influence in the Mauritanian upwelling area, due to the prevailing NW winds. The maximum VSLS mixing ratios above the Mauritanian upwelling were 8.92 ppt for bromoform, 3.14 ppt for dibromomethane and 3.29 ppt for methyl iodide, with an observed maximum range of the daily mean up to 50% for bromoform, 26% for dibromomethane and 56% for methyl iodide. The influence of various meteorological parameters - such as wind, surface air pressure, surface air and surface water temperature, humidity and marine atmospheric boundary layer (MABL) height - on VSLS concentrations and fluxes was investigated. The strongest relationship was found between the MABL height and bromoform, dibromomethane and methyl iodide abundances. Lowest MABL heights above the Mauritanian upwelling area coincide with highest VSLS mixing ratios and vice versa above the open ocean. Significant high anti-correlations confirm this relationship for the whole cruise. We conclude that especially above oceanic upwelling systems, in addition to sea-air fluxes, MABL height variations can influence atmospheric VSLS mixing ratios, occasionally leading to elevated atmospheric abundances. This may add to the postulated missing VSLS sources in the Mauritanian upwelling region (Quack et al., 2007).

Fallout and concentration of aerosol-borne radioactivity over the Atlantic (Appendix)

Measurements of atmospheric radioactivity attached to aerosols are described. Fallout was collected in a vessel of large area. Emphasis was on separation of "wet" and "dry" samples. For strontium 90 a ratio of "wet" to "dry" fallout of 5:1 has been found independent of latitude. The total fallout was smaller than comparable values from continents because of very small amounts of rainfall in the equatorial zone. In order to achieve consistency in the global balance a better knowledge not only of radioactivity but also of precipitation over the ocean is required. Fallout of Ra-D clearly shows the ITC as a barrier for the latitudinal movement of near sea-surface air masses. The concentration of short-lived emanation daughters shows large variations according to varying geographic conditions. A variation with time could not be explained. The specific activity of long-lived radioactive substances shows the expected effect of the ITC as well as a seasonal diminuation of average concentration, similar to that measured at Heidelberg.

Meteorological measurements from the Fourcade and Polar Club Glacier, Warszawa Icefield, King George Island, West Antarctica

The South Shetland Islands are located at the northern tip of the AP which is among the fastest warming regions on Earth. The islands are especially vulnerable to climate change due to their exposure to transient low-pressure systems and their maritime climate. Surface air temperature increases (2.5K in 50 years) are concurrent with retreating glacier fronts, an increase in melt areas, ice surface lowering and rapid break-up and disintegration of ice shelves. We have compiled a unique meteorological data set for the King George Island (KGI)/Isla 25 de Mayo, the largest of the South Shetland Islands. It comprises high-temporal resolution and spatially distributed observations of surface air temperature, wind directions and wind velocities, as well as glacier ice temperatures in profile with a fully equipped automatic weather station on the Warszawa Icefield, from November 2010 and ongoing. In combination with two long-term synoptic datasets (40 and 10 years, respectively) and NCEP/NCAR reanalysis data, we have looked at changes in the climatological drivers of the glacial melt processes, and the sensitivity of the inland ice cap with regard to winter melting periods and pressure anomalies. The analysis has revealed, a positive trend of 5K over four decades in minimum surface air temperatures for winter months, clearly exceeding the published annual mean statistics, associated to a decrease in mean monthly winter sea level pressure. This concurs with a positive trend in the Southern Annular Mode (SAM) index, which gives a measure for the strength and extension of the Antarctic vortex. We connect this trend with a higher frequency of low-pressure systems hitting the South Shetland Islands during austral winter, bringing warm and moist air masses from lower latitudes. Due to its exposure to the impact of transient synoptic weather systems, the ice cap of KGI is especially vulnerable to changes during winter glacial mass accumulation period. A revision of seasonal changes in adiabatic air temperature lapse rates and their dependency on exposure and elevation has shown a clear decoupling of atmospheric surface layers between coastal areas and the higher-elevation ice cap, showing the higher sensitivity to free atmospheric flow and synoptic changes. Observed surface air temperature lapse rates show a high variability during winter months (standard deviations up to ±1.0K/100 m), and a distinct spatial variability reflecting the impact of synoptic weather patterns. The observed advective conditions bringing warm, moist air with high temperatures and rain, lead to melt conditions on the ice cap, fixating surface air temperatures to the melting point. This paper assesses the impact of large-scale atmospheric circulation variability and climatic changes on the atmospheric surface layer and glacier mass accumulation of the upper ice cap during winter season for the Warszawa Icefield on KGI.

Annual-layer thickness, d18O and sodium values on Akademii Nauk ice core (AD 900-1998) based on core chronology AN 2012

Understanding recent Arctic climate change requires detailed information on past changes, in particular on a regional scale. The extension of the depth-age relation of the Akademii Nauk (AN) ice core from Severnaya Zemlya (SZ) to the last 1100 yr provides new perspectives on past climate fluctuations in the Barents and Kara seas region. Here, we present the easternmost high-resolution ice-core climate proxy records (d18O and sodium) from the Arctic. Multi-annual AN d18O data as near-surface air-temperature proxies reveal major temperature changes over the last millennium, including the absolute minimum around 1800 and the unprecedented warming to a double-peak maximum in the early 20th century. The long-term cooling trend in d18O is related to a decline in summer insolation but also to the growth of the AN ice cap as indicated by decreasing sodium concentrations. Neither a pronounced Medieval Climate Anomaly nor a Little Ice Age are detectable in the AN d18O record. In contrast, there is evidence of several abrupt warming and cooling events, such as in the 15th and 16th centuries, partly accompanied by corresponding changes in sodium concentrations. These abrupt changes are assumed to be related to sea-ice cover variability in the Barents and Kara seas region, which might be caused by shifts in atmospheric circulation patterns. Our results indicate a significant impact of internal climate variability on Arctic climate change in the last millennium.

Antarctic dataset in NetCDF format

The dataset described in this document has been put together for the purposes of numerical ice sheet modelling of the Antarctic Ice Sheet (AIS), containing data on the ice sheet configuration (e.g. ice surface and ice thickness) and boundary conditions, such as the surface air temperature and accumulation. It is now possible to download a community ice sheet model (e.g. Glimmer-CISM, Rutt et al., 2009 doi:10.1029/2008JF001015), but without adequate data it is difficult to utilise such models. More specifically, ice sheet models that are initialised and run forward from the present day ice sheet configuration, need input data to represent the present-day ice sheet configuration as closely as possible (unlike those spun-up from ice free conditions, which only require the bed/bathymetry). Whilst the BEDMAP dataset (Lythe et al., 2001) was a step forward when it was made, there are a number of inconsistencies within the dataset (see Section 3), and since its release, more data has become available. The dataset described here incorporates some major new datasets (e.g. AGASEA/BBAS ice thickness, Nitsche et al. (2006) bathymetry doi:10.1029/2007GC001694), but by no means incorporates all the new data available. This considerable task is left for a 'BEDMAP2', (an updated version of BEDMAP), however, the processing carried out in this document illustrates the requirements of a dataset for the purpose of high resolution ice sheet modelling, and bridges the gap until a BEDMAP2 is published. It is envisaged, however, that updated versions of the data set will be made available periodically when new regional data sets become available and can be readily incorporated.

Temperature measured in permafrost borehole at ICDP site 5011-3 in El'gygytgyn Crater, 2008-2011

This study focuses on the temperature field observed in boreholes drilled as part of interdisciplinary scientific campaign targeting the El'gygytgyn Crater Lake in NE Russia. Temperature data are available from two sites: the lake borehole 5011-1 located near the center of the lake reaching 400 m depth, and the land borehole 5011-3 at the rim of the lake, with a depth of 140 m. Constraints on permafrost depth and past climate changes are derived from numerical simulation of the thermal regime associated with the lake-related talik structure. The thermal properties of the subsurface needed for these simulations are based on laboratory measurements of representative cores from the quaternary sediments and the underlying impact-affected rock, complemented by further information from geophysical logs and data from published literature. The temperature observations in the lake borehole 5011-1 are dominated by thermal perturbations related to the drilling process, and thus only give reliable values for the lowermost value in the borehole. Undisturbed temperature data recorded over more than two years are available in the 140 m deep land-based borehole 5011-3. The analysis of these observations allows determination of not only the recent mean annual ground surface temperature, but also the ground surface temperature history, though with large uncertainties. Although the depth of this borehole is by far too insufficient for a complete reconstruction of past temperatures back to the Last Glacial Maximum, it still affects the thermal regime, and thus permafrost depth. This effect is constrained by numerical modeling: assuming that the lake borehole observations are hardly influenced by the past changes in surface air temperature, an estimate of steady-state conditions is possible, leading to a meaningful value of 14 ± 5 K for the post-glacial warming. The strong curvature of the temperature data in shallower depths around 60 m can be explained by a comparatively large amplitude of the Little Ice Age (up to 4 K), with low temperatures prevailing far into the 20th century. Other mechanisms, like varying porosity, may also have an influence on the temperature profile, however, our modeling studies imply a major contribution from recent climate changes.

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.

Cenozoic global ice volume and temperature simulations over the past 40 million years

Variations in global ice volume and temperature over the Cenozoic era have been investigated with a set of one-dimensional (1-D) ice-sheet models. Simulations include three ice sheets representing glaciation in the Northern Hemisphere, i.e. in Eurasia, North America and Greenland, and two separate ice sheets for Antarctic glaciation. The continental mean Northern Hemisphere surface-air temperature has been derived through an inverse procedure from observed benthic d18O records. These data have yielded a mutually consistent and continuous record of temperature, global ice volume and benthic d18O over the past 35 Ma. The simple 1-D model shows good agreement with a comprehensive 3-D ice-sheet model for the past 3 Ma. On average, differences are only 1.0°C for temperature and 6.2 m for sea level. Most notably, over the 35 Ma period, the reconstructed ice volume-temperature sensitivity shows a transition from a climate controlled by Southern Hemisphere ice sheets to one controlled by Northern Hemisphere ice sheets. Although the transient behaviour is important, equilibrium experiments show that the relationship between temperature and sea level is linear and symmetric, providing limited evidence for hysteresis. Furthermore, the results show a good comparison with other simulations of Antarctic ice volume and observed sea level.

(Table 1) Characteristics of air temperature, sea surface temperature, sea surface salinity, nutrients, and total particulate matter for Admiralty Bay and Potter Cove, King George Island, Antarctica

Since the early 1990s, phytoplankton has been studied and monitored in Potter Cove (PC) and Admiralty Bay (AB), King George/25 de Mayo Island (KGI), South Shetlands. Phytoplankton biomass is typically low compared to other Antarctic shelf environments, with average spring - summer values below 1 mg chlorophyll a (Chl a)/m**3. The physical conditions in the area (reduced irradiance induced by particles originated from the land, intense winds) limit the coastal productivity at KGI, as a result of shallow Sverdrup's critical depths (Zc) and large turbulent mixing depths (Zt). In January 2010 a large phytoplankton bloom with a maximum of around 20 mg Chl a/m**3, and monthly averages of 4 (PC) and 6 (AB) mg Chl a/m**3, was observed in the area, making it by far the largest recorded bloom over the last 20 yr. Dominant phytoplankton species were the typical bloom-forming diatoms that are usually found in the western Antarctic Peninsula area. Anomalously cold air temperature and dominant winds from the eastern sector seem to explain adequate light : mixing environment. Local physical conditions were analyzed by means of the relationship between Zc and Zt, and conditions were found adequate for allowing phytoplankton development. However, a multiyear analysis indicates that these conditions may be necessary but not sufficient to guarantee phytoplankton accumulation. The relation between maximum Chl a values and air temperature suggests that bottom-up control would render such large blooms even less frequent in KGI under the warmer climate expected in the area during the second half of the present century.