The performance of single- and multi-proxy transfer functions (testate amoebae, bryophytes, vascular plants) for reconstructing mire surface wetness and pH

Peatlands are widely exploited archives of paleoenvironmental change. We developed and compared multiple transfer functions to infer peatland depth to the water table (DWT) and pH based on testate amoeba (percentages, or presence/absence), bryophyte presence/absence, and vascular plant presence/absence data from sub-alpine peatlands in the SE Swiss Alps in order to 1) compare the performance of single-proxy vs. multi-proxy models and 2) assess the performance of presence/absence models. Bootstrapping cross-validation showing the best performing single-proxy transfer functions for both DWT and pH were those based on bryophytes. The best performing transfer functions overall for DWT were those based on combined testate amoebae percentages, bryophytes and vascular plants; and, for pH, those based on testate amoebae and bryophytes. The comparison of DWT and pH inferred from testate amoeba percentages and presence/absence data showed similar general patterns but differences in the magnitude and timing of some shifts. These results show new directions for paleoenvironmental research, 1) suggesting that it is possible to build good-performing transfer functions using presence/absence data, although with some loss of accuracy, and 2) supporting the idea that multi-proxy inference models may improve paleoecological reconstruction. The performance of multi-proxy and single-proxy transfer functions should be further compared in paleoecological data.

Integrated Segmentation of Brain Tumor Images for Radiotherapy and Neurosurgery

Image-based modeling of tumor growth combines methods from cancer simulation and medical imaging. In this context, we present a novel approach to adapt a healthy brain atlas to MR images of tumor patients. In order to establish correspondence between a healthy atlas and a pathologic patient image, tumor growth modeling in combination with registration algorithms is employed. In a first step, the tumor is grown in the atlas based on a new multi-scale, multi-physics model including growth simulation from the cellular level up to the biomechanical level, accounting for cell proliferation and tissue deformations. Large-scale deformations are handled with an Eulerian approach for finite element computations, which can operate directly on the image voxel mesh. Subsequently, dense correspondence between the modified atlas and patient image is established using nonrigid registration. The method offers opportunities in atlasbased segmentation of tumor-bearing brain images as well as for improved patient-specific simulation and prognosis of tumor progression.

The Aeolian Flux of Calcium, Chloride and Nitrate to the McMurdo Dry Valleys Landscape: Evidence from Snow Pit Analysis

We have determined the flux of calcium, chloride and nitrate to the McMurdo Dry Valleys region by analysing snow pits for their chemical composition and their snow accumulation using multiple records spanning up to 48 years. The fluxes demonstrate patterns related to elevation and proximity to the ocean. In general, there is a strong relationship between the nitrate flux and snow accumulation, indicating that precipitation rates may have a great influence over the nitrogen concentrations in the soils of the valleys. Aeolian dust transport plays an important role in the deposition of some elements (e.g. C(2+)) into the McMurdo Dry Valleys' soils. Because of the antiquity of some of the soil surfaces in the McMurdo Dry Valleys regions, the accumulated atmospheric flux of salts to the soils has important ecological consequences. Although precipitation may be an important mechanism of salt deposition to the McMurdo Dry Valley surfaces, it is poorly understood because of difficulties in measurement and high losses from sublimation.

Stable-Isotope Time Series and Precipitation Origin from Firn-Core and Snow Samples, Altai Glaciers, Siberia

In the summers of 2001 and 2002, glacio-climatological research was performed at 4110-4120 m a.s.l. on the Belukha snow/firn plateau, Siberian Altai. Hundreds of samples from snow pits and a 21 m snow/firn core were collected to establish the annual/seasonal/monthly depth-accumulation scale, based on stable-isotope records, stratigraphic analyses and meteorological and synoptic data. The fluctuations of water stable-isotope records show well-preserved seasonal variations. The delta(18)O and delta D relationships in precipitation, snow pits and the snow/firn core have the same slope to the covariance as that of the global meteoric water line. The origins of precipitation nourishing the Belukha plateau were determined based on clustering analysis of delta(18)O and d-excess records and examination of synoptic atmospheric patterns. Calibration and validation of the developed clusters occurred at event and monthly timescales with about 15% uncertainty. Two distinct moisture sources were shown: oceanic sources with d-excess < 12 parts per thousand, and the Aral-Caspian closed drainage basin sources with d-excess > 12 parts per thousand. Two-thirds of the annual accumulation was from oceanic precipitation, of which more than half had isotopic ratios corresponding to moisture evaporated over the Atlantic Ocean. Precipitation from the Arctic/Pacific Ocean had the lowest deuterium excess, contributing one-tenth to annual accumulation.

A review of air–ice chemical and physical interactions (AICI): liquids, quasi-liquids, and solids in snow

Snow in the environment acts as a host to rich chemistry and provides a matrix for physical exchange of contaminants within the ecosystem. The goal of this review is to summarise the current state of knowledge of physical processes and chemical reactivity in surface snow with relevance to polar regions. It focuses on a description of impurities in distinct compartments present in surface snow, such as snow crystals, grain boundaries, crystal surfaces, and liquid parts. It emphasises the microscopic description of the ice surface and its link with the environment. Distinct differences between the disordered air–ice interface, often termed quasi-liquid layer, and a liquid phase are highlighted. The reactivity in these different compartments of surface snow is discussed using many experimental studies, simulations, and selected snow models from the molecular to the macro-scale. Although new experimental techniques have extended our knowledge of the surface properties of ice and their impact on some single reactions and processes, others occurring on, at or within snow grains remain unquantified. The presence of liquid or liquid-like compartments either due to the formation of brine or disorder at surfaces of snow crystals below the freezing point may strongly modify reaction rates. Therefore, future experiments should include a detailed characterisation of the surface properties of the ice matrices. A further point that remains largely unresolved is the distribution of impurities between the different domains of the condensed phase inside the snowpack, i.e. in the bulk solid, in liquid at the surface or trapped in confined pockets within or between grains, or at the surface. While surface-sensitive laboratory techniques may in the future help to resolve this point for equilibrium conditions, additional uncertainty for the environmental snowpack may be caused by the highly dynamic nature of the snowpack due to the fast metamorphism occurring under certain environmental conditions. Due to these gaps in knowledge the first snow chemistry models have attempted to reproduce certain processes like the long-term incorporation of volatile compounds in snow and firn or the release of reactive species from the snowpack. Although so far none of the models offers a coupled approach of physical and chemical processes or a detailed representation of the different compartments, they have successfully been used to reproduce some field experiments. A fully coupled snow chemistry and physics model remains to be developed.

Genetic diversity and differentiation follow secondary succession in a multi-species study on woody plants from subtropical China

Aims: Species diversity and genetic diversity may be affected in parallel by similar environmental drivers. However, genetic diversity may also be affected independently by habitat characteristics. We aim at disentangling relationships between genetic diversity, species diversity and habitat characteristics of woody species in subtropical forest. Methods: We studied 11 dominant tree and shrub species in 27 plots in Gutianshan, China, and assessed their genetic diversity (Ar) and population differentiation (F’ST) with microsatellite markers. We tested if Ar and population specific F’ST were correlated to local species diversity and plot characteristics. Multi-model inference and model averaging were used to determine the relative importance of each predictor. Additionally we tested for isolation-by-distance and isolation-by-elevation by regressing pairwise F’ST against pairwise spatial and elevational distances. Important findings: Genetic diversity was not related to species diversity for any of the study species. Thus, our results do not support joint effects of habitat characteristics on these two levels of biodiversity. Instead, genetic diversity in two understory shrubs, Rhododendron simsii and Vaccinium carlesii, was affected by plot age with decreasing genetic diversity in successionally older plots. Population differentiation increased with plot age in Rhododendron simsii and Lithocarpus glaber. This shows that succession can reduce genetic diversity within, and increase genetic diversity between populations. Furthermore, we found four cases of isolation-by-distance and two cases of isolation-by-elevation. The former indicates inefficient pollen and seed dispersal by animals whereas the latter might be due to phenological asynchronies. These patterns indicate that succession can affect genetic diversity without parallel effects on species diversity and that gene flow in a continuous subtropical forest can be restricted even at a local scale.

The assumptions and effects of conservative procedures in low -dose risk assessment for setting safety standards

Conservative procedures in low-dose risk assessment are used to set safety standards for known or suspected carcinogens. However, the assumptions upon which the methods are based and the effects of these methods are not well understood.^ To minimize the number of false-negatives and to reduce the cost of bioassays, animals are given very high doses of potential carcinogens. Results must then be extrapolated to much smaller doses to set safety standards for risks such as one per million. There are a number of competing methods that add a conservative safety factor into these calculations.^ A method of quantifying the conservatism of these methods was described and tested on eight procedures used in setting low-dose safety standards. The results using these procedures were compared by computer simulation and by the use of data from a large scale animal study.^ The method consisted of determining a "true safe dose" (tsd) according to an assumed underlying model. If one assumed that Y = the probability of cancer = P(d), a known mathematical function of the dose, then by setting Y to some predetermined acceptable risk, one can solve for d, the model's "true safe dose".^ Simulations were generated, assuming a binomial distribution, for an artificial bioassay. The eight procedures were then used to determine a "virtual safe dose" (vsd) that estimates the tsd, assuming a risk of one per million. A ratio R = ((tsd-vsd)/vsd) was calculated for each "experiment" (simulation). The mean R of 500 simulations and the probability R $<$ 0 was used to measure the over and under conservatism of each procedure.^ The eight procedures included Weil's method, Hoel's method, the Mantel-Byran method, the improved Mantel-Byran, Gross's method, fitting a one-hit model, Crump's procedure, and applying Rai and Van Ryzin's method to a Weibull model.^ None of the procedures performed uniformly well for all types of dose-response curves. When the data were linear, the one-hit model, Hoel's method, or the Gross-Mantel method worked reasonably well. However, when the data were non-linear, these same methods were overly conservative. Crump's procedure and the Weibull model performed better in these situations. ^

Mean snow accumulation rate, and mean methane sulfonate, chloride and nitrate concentration in snow pits and firn cores of Dronning Maud Land

We quantified postdepositional losses of methane sulfonate (MSA-), nitrate, and chloride at the European Project for Ice Coring in Antarctica (EPICA) drilling site in Dronning Maud Land (DML) (75°S, 0°E). Analyses of four intermediate deep firn cores and 13 snow pits were considered. We found that about 26 ± 13% of the once deposited nitrate and typically 51 ± 20% of MSA- were lost, while for chloride, no significant depletion could be observed in firn older than one year. Assuming a first order exponential decay rate, the characteristic e-folding time for MSA- is 6.4 ± 3 years and 19 ± 6 years for nitrate. It turns out that for nitrate and MSA- the typical mean concentrations representative for the last 100 years were reached after 5.4 and 6.5 years, respectively, indicating that beneath a depth of around 1.2-1.4 m postdepositional losses can be neglected. In the area of investigation, only MSA- concentrations and postdepositional losses showed a distinct dependence on snow accumulation rate. Consequently, MSA- concentrations archived at this site should be significantly dependent on the variability of annual snow accumulation, and we recommend a corresponding correction. With a simple approach, we estimated the partial pressure of the free acids MSA, HNO3, and HCl on the basis of Henry's law assuming that ionic impurities of the bulk ice matrix are localized in a quasi-brine layer (QBL). In contrast to measurements, this approach predicts a nearly complete loss of MSA-, NO3 - , and Cl-.

Ice thickness measurements in the Lincoln Sea and adjacent Arctic Ocean

Results of helicopter-borne electromagnetic measurements of total (ice plus snow) sea-ice thickness performed in May 2004 and 2005 in the Lincoln Sea and adjacent Arctic Ocean up to 86° N are presented. Thickness distributions south of 84° N are dominated by multi-year ice with modal thicknesses of 3.9 m in 2004 and 4.2 m in 2005 (mean thicknesses 4.67 and 5.18 m, respectively). Modal and mean snow thickness on multi-year ice amounted to 0.18 and 0.30 m in 2004, and 0.28 and 0.35 m in 2005. There are also considerable amounts of 0.9-2.2 m thick first-year ice (modal thickness), mostly representing ice formed in the recurring, refrozen Lincoln Polynya. Results are in good agreement with ground-based electromagnetic thickness measurements and with ice types demarcated in satellite synthetic aperture radar imagery. Four drifting buoys deployed in 2004 between 86° N and 84.5° N show a similar pattern of a mean southward drift of the ice pack of 83 ± 18 km between May 2004 and April 2005, towards the coast of Ellesmere Island and Nares Strait. The resulting area decrease of 26% between the buoys and the coast is larger than the observed thickness increase south of 84° N. This points to the importance of shear in a narrow band along the coast, and of ice export through Nares Strait in removing ice from the study region.

(Table 1) Characteristics of snow pits on De Geerfonna and Vestfonna during the IPY 2007-2010

Extensive glaciological field measurements were carried out on the ice cap Vestfonna as well as on the minor ice body De Geerfonna (Nordaustlandet, Svalbard) within the framework of IPY Kinnvika. Field campaigns were conducted during the period 2007-2010 in spring (April/May) and summer (August). In this study we compile and present snow cover information obtained from 22 snow pits that were dug on Vestfonna during this period. Locations are along two transects on the northwestern, land terminating slope of the ice cap, on its central summit, Ahlmann Summit, and at a set of several other locations in the eastern and northern part of the ice cap. Snow-cover information acquired from four snow pits on adjacent De Geerfonna is also incorporated in this study. Field data are analysed regarding snow stratigraphy, snow density, snow hardness and snow temperature. Results reveal mean snow densities of around 400 kg/m**3 for the snowpack of Vestfonna with no apparent spatial or interannual variability. A distinctly higher value of more than 450 kg/m**3 was obtained for De Geerfonna. A spatial comparison of snow water equivalents above the previous end-of-summer surface serves for obtaining insights into the spatial distribution of snow accumulation across Vestfonna. Altitude was found to be the only significant spatial parameter for controlling snow accumulation across the ice cap.

Field measurements of the atmosphere, ocean, sea ice and sub-ice platelet layer at Atka Bay in 2013

Ice shelves strongly interact with coastal Antarctic sea ice and the associated ecosystem by creating conditions favourable to the formation of a sub-ice platelet layer. The close investigation of this phenomenon and its seasonal evolution remain a challenge due to logistical constraints and a lack of suitable methodology. In this study, we characterize the seasonal cycle of Antarctic fast ice adjacent to the Ekström Ice Shelf in the eastern Weddell Sea. We used a thermistor chain with the additional ability to record the temperature response induced by cyclic heating of resistors embedded in the chain. Vertical sea-ice temperature and heating profiles obtained daily between November 2012 and February 2014 were analyzed to determine sea-ice and snow evolution, and to calculate the basal energy budget. The residual heat flux translated into an ice-volume fraction in the platelet layer of 0.18 ± 0.09, which we reproduced by a independent model simulation and agrees with earlier results. Manual drillings revealed an average annual platelet-layer thickness increase of at least 4m, and an annual maximum thickness of 10m beneath second-year sea ice. The oceanic contribution dominated the total sea-ice production during the study, effectively accounting for up to 70% of second-year sea-ice growth. In summer, an oceanic heat flux of 21 W/m**2 led to a partial thinning of the platelet layer. Our results further show that the active heating method, in contrast to the acoustic sounding approach, is well suited to derive the fast-ice mass balance in regions influenced by ocean/ice-shelf interaction, as it allows sub-diurnal monitoring of the platelet-layer thickness.

Ice and snow characteristics in Prydz Bay in 2006

Under the Chinese National Antarctic Research Expedition program in 2006, the annual thermal mass balance of landfast ice in the vicinity of Zhongshan Station, Prydz Bay, east Antarctica, was investigated. Sea ice formed from mid-February onward, and maximum ice thickness occurred in late November. Snow cover remained thin, and blowing snow caused frequent redistribution of the snow. The vertical ice salinity showed a 'question-mark-shaped' profile for most of the ice growth season, which only turned into an 'I-shaped' profile after the onset of ice melt. The oceanic heat flux as estimated from a flux balance at ice-ocean interface using internal ice temperatures decreased from 11.8 (±3.5) W/m**2 in April to an annual minimum of 1.9 (±2.4) W/m**2 in September. It remained low through late November, in mid-December it increased sharply to about 20.0 W/m**2. Simulations applying the modified versions of Stefan's law, taking account the oceanic heat flux and ice-atmosphere coupling, compare well with observed ice growth. There was no obvious seasonal cycle for the thermal conductivity of snow cover, which was also derived from internal ice temperatures. Its annual mean was 0.20 (±0.04) W/m/°C.

(Table 1) Summary statistics of in situ sea-ice characteristics of Ice Station Belgica sites Fabra, Patria/Liège and Brussels in October 2007

Three ice type regimes at Ice Station Belgica (ISB), during the 2007 International Polar Year SIMBA (Sea Ice Mass Balance in Antarctica) expedition, were characterized and assessed for elevation, snow depth, ice freeboard and thickness. Analyses of the probability distribution functions showed great potential for satellite-based altimetry for estimating ice thickness. In question is the required altimeter sampling density for reasonably accurate estimation of snow surface elevation given inherent spatial averaging. This study assesses an effort to determine the number of laser altimeter 'hits' of the ISB floe, as a representative Antarctic floe of mixed first- and multi-year ice types, for the purpose of statistically recreating the in situ-determined ice-thickness and snow depth distribution based on the fractional coverage of each ice type. Estimates of the fractional coverage and spatial distribution of the ice types, referred to as ice 'towns', for the 5 km**2 floe were assessed by in situ mapping and photo-visual documentation. Simulated ICESat altimeter tracks, with spot size ~70 m and spacing ~170 m, sampled the floe's towns, generating a buoyancy-derived ice thickness distribution. 115 altimeter hits were required to statistically recreate the regional thickness mean and distribution for a three-town assemblage of mixed first- and multi-year ice, and 85 hits for a two-town assemblage of first-year ice only: equivalent to 19.5 and 14.5 km respectively of continuous altimeter track over a floe region of similar structure. Results have significant implications toward model development of sea-ice sampling performance of the ICESat laser altimeter record as well as maximizing sampling characteristics of satellite/airborne laser and radar altimetry missions for sea-ice thickness.

Particulate matter and organic compounds in snow, ice, and water of the Lazarev and Cooperation Seas

Integrated studies of particulate matter and organic compounds in surface waters and the snow-ice cover by means of geochemical (concentrations of the particulate matter, Corg, hydrocarbons, lipids, and chlorophyll a) and optical techniques were performed in the Southern Ocean and in the East Atlantic Ocean along the vessel's route: Africa - Antarctica - Africa - St. Petersburg. Correlations between studied compounds were found. It was shown that supply of pollutants affects not only concentrations but also proportions of the considered compounds. New data were obtained on the processes of accumulation of particulate matter and organic compounds under ice formation.