An individual based model spatial distribution, production and biomass for Calanus finmarchicus in the Norwegian Sea

The copepod Calanus finmarchicus is the dominant species of the meso-zooplankton in the Norwegian Sea, and constitutes an important link between the phytoplankton and the higher trophic levels in the Norwegian Sea food chain. An individualbased model for C. finmarchicus, based on super-individuals and evolving traits for behaviour, stages, etc., is two-way coupled to the NORWegian ECOlogical Model system (NORWECOM). One year of modelled C. finmarchicus spatial distribution, production and biomass are found to represent observations reasonably well. High C. finmarchicus abundance is found along the Norwegian shelf-break in the early summer, while the overwintering population is found along the slope and in the deeper Norwegian Sea basins. The timing of the spring bloom is generally later than in the observations. Annual Norwegian Sea production is found to be 29 million tonnes of carbon and a production to biomass (P/B) ratio of 4.3 emerges. Sensitivity tests show that the modelling system is robust to initial values of behavioural traits and with regards to the number of super-individuals simulated given that this is above about 50,000 individuals. Experiments with the model system indicate that it provides a valuable tool for studies of ecosystem responses to causative forces such as prey density or overwintering population size. For example, introducing C. finmarchicus food limitations reduces the stock dramatically, but on the other hand, a reduced stock may rebuild in one year under normal conditions. The NetCDF file contains model grid coordinates and bottom topography.

Treeline dynamics in Siberia under changing climates as inferred from an individual-based model for Larix, supplementary material

Siberian boreal forests are expected to expand northwards in the course of global warming. However, processes of the treeline ecotone transition, as well astiming and related climate feedbacks are still not understood. Here, we present 'Larix Vegetation Simulator' LAVESI, an individual-based spatially-explicit model that can simulate Larix gmelinii (RUPR.) RUPR. stand dynamics in an attempt to improve our understanding about past and future treeline movements under changing climates. The relevant processes (growth, seed production and dispersal, establishment and mortality) are incorporated and adjusted to observation data mainly gained from the literature. Results of a local sensitivity analysis support the robustness of the model's parameterization by giving relatively small sensitivity values. We tested the model by simulating tree stands under modern climate across the whole Taymyr Peninsula, north-central Siberia (c. 64-80° N; 92-119° E). We find tree densities similar to observed forests in the northern to mid-treeline areas, but densities are overestimated in the southern parts of the simulated region. Finally, from a temperature-forcing experiment, we detect that the responses of tree stands lag the hypothetical warming by several decades, until the end of 21st century. With our simulation experiments we demonstrate that the newly-developed model captures the dynamics of the Siberian latitudinal treeline.

Chemistry of minerals from amphibole-chlorite veins of ODP Site 209-1270

We examined small-scale shear zones in drillcore samples of abyssal peridotites from the Mid-Atlantic Ridge. These shear zones are associated with veins consisting of chlorite + actinolite/tremolite assemblages, with accessory phases zircon and apatite, and they are interpreted as altered plagiogranite melt impregnations, which originate from hydrous partial melting of gabbroic intrusion in an oceanic detachment fault. Ti-in-zircon thermometry yields temperatures around 820°C for the crystallization of the evolved melt. Reaction path modeling indicates that the alteration assemblage includes serpentine of the adjacent altered peridotites. Based on the model results, we propose that formation of chlorite occurred at higher temperatures than serpentinization, thus leading to strain localization around former plagiogranites during alteration. The detachment fault represents a major pathway for fluids through the oceanic crust, as evidenced by extremely low d18O of altered plagiogranite veins (+3.0-4.2 per mil) and adjacent serpentinites (+ 2.6-3.7 per mil). The uniform oxygen isotope data indicate that fluid flow in the detachment fault system affected veins and adjacent host serpentinites likewise.

Crustal uplift rates of 4 GPS stations near Jakobshavn Isbrae (2006-2009)

We analyze 2006-2009 data from four continuous Global Positioning System (GPS) receivers located between 5 and 150 km from the glacier Jakobshavn Isbrae, West Greenland. The GPS stations were established on bedrock to determine the vertical crustal motion due to the unloading of ice from Jakobshavn Isbrae. All stations experienced uplift, but the uplift rate at Kangia North, only 5 km from the glacier front, was about 10 mm/yr larger than the rate at Ilulissat, located only ~45 km further away. This suggests that most of the uplift is due to the unloading of the Earth's surface as Jakobshavn thins and loses mass. Our estimate of Jakobshavn's contribution to uplift rates at Kangia North and Ilulissat are 14.6 ± 1.7 mm/yr and 4.9 ± 1.1 mm/yr, respectively. The observed rates are consistent with a glacier thinning model based on repeat altimeter surveys from NASA's Airborne Topographic Mapper (ATM), which shows that Jakobshavn lost mass at an average rate of 22 ± 2 km**3/yr between 2006 and 2009. At Kangia North and Ilulissat, the predicted uplift rates computed using thinning estimates from the ATM laser altimetry are 12.1 ± 0.9 mm/yr and 3.2 ± 0.3 mm/yr, respectively. The observed rates are slightly larger than the predicted rates. The fact that the GPS uplift rates are much larger closer to Jakobshavn than further away, and are consistent with rates inferred using the ATM-based glacier thinning model, shows that GPS measurements of crustal motion are a potentially useful method for assessing ice-mass change models.

(Table 1) Bioclimatic subzones and soil characteristics along the Yamal Arctic Transect, northwestern Siberia

Sustainability of tundra vegetation under changing climate on the Yamal Peninsula, northwestern Siberia, home to the world's largest area of reindeer husbandry, is of crucial importance to the local native community. An integrated investigation is needed for better understanding of the effects of soils, climate change and grazing on tundra vegetation in the Yamal region. In this study we applied a nutrient-based plant community model - ArcVeg - to evaluate how two factors (soil organic nitrogen (SON) levels and grazing) interact to affect tundra responses to climate warming across a latitudinal climatic gradient on the Yamal Peninsula. Model simulations were driven by field-collected soil data and expected grazing patterns along the Yamal Arctic Transect (YAT), within bioclimate subzones C (high arctic), D (northern low arctic) and E (southern low arctic). Plant biomass and NPP (net primary productivity) were significantly increased with warmer bioclimate subzones, greater soil nutrient levels and temporal climate warming, while they declined with higher grazing frequency. Temporal climate warming of 2 °C caused an increase of 665 g/m**2 in total biomass at the high SON site in subzone E, but only 298 g/m**2 at the low SON site. When grazing frequency was also increased, total biomass increased by only 369 g/m**2 at the high SON site in contrast to 184 g/m**2 at the low SON site in subzone E. Our results suggest that high SON can support greater plant biomass and plant responses to climate warming, while low SON and grazing may limit plant response to climate change. In addition to the first order factors (SON, bioclimate subzones, grazing and temporal climate warming), interactions among these significantly affect plant biomass and productivity in the arctic tundra and should not be ignored in regional scale studies.

Geochemistry of ODP Site 135-834 basalts

New major, trace element, and isotope data (Pb, Sr, and Nd) reveal an impressive compositional variation in the basalts recovered from Site 834. Major element compositions span almost the entire range observed in basalts from the modern axial systems of the Lau Basin, and variations are consistent with low-pressure fractionation of a mid-ocean-ridge-basalt (MORB)-like parent, in which plagioclase crystallization has been somewhat suppressed. Trace element compositions deviate from MORB in all but one unit (Unit 7) and show enrichments in large-ion-lithophile elements (LILEs) relative to high-field-strength elements (HFSEs) more typically associated with island-arc magmas. The Pb-isotope ratios define linear trends that extend from the field of Pacific MORB to highly radiogenic values similar to those observed in rocks from the northernmost islands of the Tofua Arc. The Sr-isotope compositions also show significant variation, and these too project from radiogenic values back into the field for Pacific MORB. The variations in key trace element and isotopic features are consistent with magma mixing between two relatively mafic melts: one represented by Pacific MORB, and the other by a magma similar to those erupted on 'Eua when it was part of the original Tongan arc, or perhaps members of the Lau Volcanic Group (LVG). Based on our model, the most radiogenic compositions (Units 2 and 8) represent approximately 50:50 mixtures of these MORB and arc end-members. Magma mixing requires that both components are simultaneously available, and implies that melts have not shown a compositional progression from arc-like to MORB-like with extension at this locality. Rather, it is apparent that essentially pristine MORB can erupt as one of the earliest products of backarc initiation. Indeed, repetition of isotopic and trace element signatures with depth suggests that eruptions have been triggered by periodic injections of fresh MORB melts into the source regions of these magmas. The slow and almost amagmatic extension of the original arc complex envisaged to explain the observed chemistry is also consistent with the horst-and-graben topography of the western side of the Lau Basin. Given the similarities between basalts erupted at the modern Lau Basin spreading centers and MORB from the Indian Ocean, the overwhelming evidence for involvement of mantle similar to Pacific MORB in the petrogenesis of basalts from Site 834 is a new and important observation. It indicates that the original arc was underlain by asthenospheric material derived from the Pacific mantle convection cell, and that this has somehow been replaced by Indian Ocean MORB during the last ~5.5 Ma.

Radiocarbon ages of sediment core PG1414-2 from Beall Lake, Windmill Islands

A mid-Holocene climate optimum is inferred from a palaeosalinity reconstruction of a closed saline lake (Beall Lake) from the Windmill Islands, East Antarctica using an expanded diatom salinity weighted averaging (WA) regression and calibration model. The addition of 14 lakes and ponds from the Windmill Islands, East Antarctica, to an existing weighted averaging regression and calibration palaeosalinity model of 33 lakes from the Vestfold Hills, East Antarctica expands the number of taxa and lakes and the range of salinity in the existing model and improves the model's predictive ability. This improved model was used to infer Holocene changes in lake water salinity in Beall Lake, Windmill Islands. Six changes in diatom-inferred salinity in Beall Lake are put into broad chronological context based on three radiocarbon dates: as the East Antarctic Ice Sheet (EAIS) retreated from the Windmill Islands during the early Holocene (~9000-8130 corr. yr BP), Beall Lake formed as a melt water-fed freshwater lake, which gradually became more saline as marine influence increased from ~8000 corr. yr BP. Between ~8000 and 4800 corr. yr BP, the diatom assemblage included planktonic marine taxa such as Chaetoceros spp. and cryophilic taxa such as Fragilariopsis cylindrus, which indicate favourable summer growth conditions. A mid-Holocene warm period produced a climate that was warmer and more humid with increased precipitation and snow accumulation. This is reflected in the Beall Lake core as a reduction in the salinity of the lake diatom assemblage from ~4800-4600 corr. yr BP. Holocene isostatic uplift rates in the Windmill Islands vary from 5-6 m/1000 yr. By applying this uplift rate, it is calculated that the bedrock would have risen above sea level by ~4000 yr BP. The Beall Lake core diatom assemblage from ~4600-2900 corr. yr BP includes both marine cryophilic and planktonic taxa together with freshwater benthic and planktonic lacustrine taxa. This mix of species indicates the emergence of the lake from the sea around ~4600 corr. yr BP. From ~2800 corr. yr BP, retreat of the ice margin led to increasing melt water inputs and associated freshening of the lake basin until ~1900 corr. yr BP. The lake basin had no oceanic influence by this time, allowing a terrestrial freshwater flora to establish and thrive for the next ~1000 yr. At ~1850 corr. yr BP, a sudden and rapid salinity change is evident in Beall Lake. A late Holocene warm period between 2000 and 1000 yr BP has been observed in ice core records from Law Dome (an ice cap abutting the Windmill Islands to the east and north). It is therefore inferred that, at ~1850 corr. yr BP, summer temperatures within the Beall Lake catchment area were much higher than present summer temperatures. The climate optimum identified in the Beall Lake core ~4800 yr BP confirms mid-Holocene warming of the Windmill Islands and suggests a synchronous mid-Holocene climate optimum occurred across coastal East Antarctica. In addition, the abrupt climate change inferred at ~1850 yr BP suggests that higher resolution sampling of sediment cores from coastal East Antarctic limnological oases will provide more evidence of rapid climate change events over coastal East Antarctica in future.

Particulate organic matter (POM) of samples from the Lena River delta

Particulate organic matter (POM) derived from permafrost soils and transported by the Lena River represents a quantitatively important terrestrial carbon pool exported to Laptev Sea sediments (next to POM derived from coastal erosion). Its fate in a future warming Arctic, i.e., its remobilization and remineralization after permafrost thawing as well as its transport pathways to and sequestration in marine sediments, is currently under debate. We present one of the first radiocarbon (14C) data sets for surface water POM within the Lena Delta sampled in the summers of 2009 - 2010 and spring 2011 (n = 30 samples). The bulk D14C values varied from -55 to -391 per mil translating into 14C ages of 395 to 3920 years BP. We further estimated the fraction of soil-derived POM to our samples based on (1) particulate organic carbon to particulate nitrogen ratios (POC : PN) and (2) on the stable carbon isotope (d13C) composition of our samples. Assuming that this phytoplankton POM has a modern 14C concentration, we inferred the 14C concentrations of the soil-derived POM fractions. The results ranged from -322 to -884 per mil (i.e., 3060 to 17 250 14C years BP) for the POC : PN-based scenario and from -261 to -944 per mil (i.e., 2370 to 23 100 14C years BP) for the d13C-based scenario. Despite the limitations of our approach, the estimated D14C values of the soil-derived POM fractions seem to reflect the heterogeneous 14C concentrations of the Lena River catchment soils covering a range from Holocene to Pleistocene ages better than the bulk POM D14C values. We further used a dual-carbon-isotope three-end-member mixing model to distinguish between POM contributions from Holocene soils and Pleistocene Ice Complex (IC) deposits to our soil-derived POM fraction. IC contributions are comparatively low (mean of 0.14) compared to Holocene soils (mean of 0.32) and riverine phytoplankton (mean of 0.55), which could be explained with the restricted spatial distribution of IC deposits within the Lena catchment. Based on our newly calculated soil-derived POM D14C values, we propose an isotopic range for the riverine soil-derived POM end member with D14C of -495 ± 153 per mil deduced from our d13C-based binary mixing model and d13C of -26.6 ± 1 per mil deduced from our data of Lena Delta soils and literature values. These estimates can help to improve the dual-carbon-isotope simulations used to quantify contributions from riverine soil POM, Pleistocene IC POM from coastal erosion, and marine POM in Siberian shelf sediments.

Helium isotope concentrations and ratios in sediments from ODP Hole 130-806B

We have determined the helium abundance and isotopic composition of seafloor carbonate sediments from the flanks of the Ontong Java Plateau, western equatorial Pacific Ocean (ODP Site 806). These results provide a two million year record of the burial flux of extraterrestrial 3He, which we believe is a proxy for the terrestrial accretion rate of interplanetary dust particles. The 3He burial flux prior to ~700 ka was relatively low, ~0.5 pcc/cm**2/kyr, but from 700 ka to the present, the burial flux gradually increased to a value of ~1.0 pcc/cm**2/kyr. 100 kyr periodicity in the 3He burial flux is apparent over the last 700 kyr and correlates with the oxygen isotope record of global climate, with high 3He burial fluxes associated with interglacial periods. This periodicity and phase are consistent with previous 3He measurements in North Atlantic sediments. Although 100 kyr periodicity in 3He burial flux is in agreement with recent predictions of the accretion rate of interplanetary dust based on a model of the orbital evolution of asteroidal debris, the measurements and predictions differ by one half cycle in phase. Nevertheless, our observations suggest the terrestrial accretion rate of interplanetary dust is controlled by orbital eccentricity and/or inclination relative to the solar-system invariable plane. Such control is a necessary but not sufficient condition for the hypothesis of that variations in extraterrestrial dust accretion modulates terrestrial climate with a 100 kyr period. We also identify several brief (<25 kyr) intervals of strongly enhanced 3He burial, possibly related to random and transient fluctuations in the accretion rate of asteroidal or cometary dust particles.

Hyperspectral downwelling irradiance during POLARSTERN cruise ANT-XXVI/1

A Monte Carlo based radiative transfer model has been developed for calculating the availability of solar radiation within the top 100 m of the ocean. The model is optimized for simulations of spatial high resolution downwelling irradiance Ed fluctuations that arise from the lensing effect of waves at the water surface. In a first step the accuracy of simulation results has been verified by measurements of the oceanic underwater light field and through intercomparison with an established radiative transfer model. Secondly the potential depth-impact of nonlinear shaped single waves, from capillary to swell waves, is assessed by considering the most favorable conditions for light focusing, i.e. monochromatic light at 490 nm, very clear oceanic water with a low chlorophyll a content of 0.1 mg/m**3 and high sun elevation. Finally light fields below irregular wave profiles accounting for realistic sea states were simulated. Our simulation results suggest that under open ocean conditions light flashes with 50% irradiance enhancements can appear down to 35 m depth, and light variability in the range of ±10% compared to the mean Ed is still possible in 100 m depth.