Zircons, isotope ratios and element analyses from tonalite and oxide gabbro of the mid-ocean ridge from ODP Hole 176-735B

A limiting factor in the accuracy and precision of U/Pb zircon dates is accurate correction for initial disequilibrium in the 238U and 235U decay chains. The longest-lived-and therefore most abundant-intermediate daughter product in the 235U isotopic decay chain is 231Pa (T1/2 = 32.71 ka), and the partitioning behavior of Pa in zircon is not well constrained. Here we report high-precision thermal ionization mass spectrometry (TIMS) U-Pb zircon data from two samples from Ocean Drilling Program (ODP) Hole 735B, which show evidence for incorporation of excess 231Pa during zircon crystallization. The most precise analyses from the two samples have consistent Th-corrected 206Pb/238U dates with weighted means of 11.9325 ± 0.0039 Ma (n = 9) and 11.920 ± 0.011 Ma (n = 4), but distinctly older 207Pb/235U dates that vary from 12.330 ± 0.048 Ma to 12.140 ± 0.044 Ma and 12.03 ± 0.24 to 12.40 ± 0.27 Ma, respectively. If the excess 207Pb is due to variable initial excess 231Pa, calculated initial (231Pa)/(235U) activity ratios for the two samples range from 5.6 ± 1.0 to 9.6 ± 1.1 and 3.5 ± 5.2 to 11.4 ± 5.8. The data from the more precisely dated sample yields estimated DPazircon/DUzircon from 2.2-3.8 and 5.6-9.6, assuming (231Pa)/(235U) of the melt equal to the global average of recently erupted mid-ocean ridge basaltic glasses or secular equilibrium, respectively. High precision ID-TIMS analyses from nine additional samples from Hole 735B and nearby Hole 1105A suggest similar partitioning. The lower range of DPazircon/DUzircon is consistent with ion microprobe measurements of 231Pa in zircons from Holocene and Pleistocene rhyolitic eruptions (Schmitt (2007; doi:10.2138/am.2007.2449) and Schmitt (2011; doi:10.1146/annurev-earth-040610-133330)). The data suggest that 231Pa is preferentially incorporated during zircon crystallization over a range of magmatic compositions, and excess initial 231Pa may be more common in zircons than acknowledged. The degree of initial disequilibrium in the 235U decay chain suggested by the data from this study, and other recent high precision datasets, leads to resolvable discordance in high precision dates of Cenozoic to Mesozoic zircons. Minor discordance in zircons of this age may therefore reflect initial excess 231Pa and does not require either inheritance or Pb loss.

Global macrozoobenthos production and energy budget data base V150731, with link to database

I developed a new model for estimating annual production-to-biomass ratio P/B and production P of macrobenthic populations in marine and freshwater habitats. Self-learning artificial neural networks (ANN) were used to model the relationships between P/B and twenty easy-to-measure abiotic and biotic parameters in 1252 data sets of population production. Based on log-transformed data, the final predictive model estimates log(P/B) with reasonable accuracy and precision (r2 = 0.801; residual mean square RMS = 0.083). Body mass and water temperature contributed most to the explanatory power of the model. However, as with all least squares models using nonlinearly transformed data, back-transformation to natural scale introduces a bias in the model predictions, i.e., an underestimation of P/B (and P). When estimating production of assemblages of populations by adding up population estimates, accuracy decreases but precision increases with the number of populations in the assemblage.

Rare earth elements determined in Antarctic ice by inductively coupled plasma

Inductively coupled plasma mass spectrometry (ICP-MS) is a suitable tool for multi-element analysis at low concentration levels. Rare earth element (REE) determinations in standard reference materials and small volumes of molten ice core samples from Antarctica have been performed with an ICP-time of flight-MS (ICP-TOF-MS) system. Recovery rates for REE in e.g. SPS-SW1 amounted to not, vert, similar ~103%, and the relative standard deviations were 3.4% for replicate analysis at REE concentrations in the lower ng/l range. Analyses of REE concentrations in Antarctic ice core samples showed that the ICP-TOF-MS technique meets the demands of restricted sample mass. The data obtained are in good agreement with ICP-Quadrupole-MS (ICP-Q-MS) and ICP-Sector Field-MS (ICP-SF-MS) results. The ICP-TOF-MS system determines accurately and precisely REE concentrations exceeding 5 ng/l while between 0.5 and 5 ng/l accuracy and precision are element dependent.

Paleomagnetic, stable oxygen isotope ratios, and U/Th activities of ODP Site 172-1062

Geomagnetic excursions are recognized as intrinsic features of the Earth's magnetic field. High-resolution records of field behaviour, captured in marine sedimentary cores, present an opportunity to determine the temporal and geometric character of the field during geomagnetic excursions and provide constraints on the mechanisms producing field variability. We present here the highest resolution record yet published of the Blake geomagnetic excursion (~125 ka) measured in three cores from Ocean Drilling Program (ODP) Site 1062 on the Blake-Bahama Outer Ridge. The Blake excursion has a controversial structure and timing but these cores have a sufficiently high sedimentation rate (~10cm/ka) to allow detailed reconstruction of the field behaviour at this site during the excursion. Palaeomagnetic measurements of the cores reveal rapid transitions (<500 yr) between the contemporary stable normal polarity and a completely reversed state of long duration which spans a stratigraphic interval of 0.7 m. We determine the duration of the reversed state during the Blake excursion using oxygen isotope stratigraphy, combined with 230Th excess measurements to assess variations in the sedimentation rates through the sections of interest. This provides an age and duration for the Blake excursion with greater accuracy and with constrained uncertainty. We date the directional excursion as falling between 129 and 122 ka with a duration for the deviation of 6.5±1.3 kyr. The long duration of this interval and the fully reversed field suggest the existence of a pseudo-stable, reversed dipole field component during the excursion and challenge the idea that excursions are always of short duration.

(Table 2) Temporal changes in phenological observations between 1978 and 2007, northern Sweden

A 30-year series (1978-2007) of photographic records were analysed to determine changes in lake ice cover, local (low elevation) and montane (high elevation) snow cover and phenological stages of mountain birch (Betula pubescens ssp. czerepanovii) at the Abisko Scientific Research Station, Sweden. In most cases, the photographic-derived data showed no significant difference in phenophase score from manually observed field records from the same period, demonstrating the accuracy and potential of using weekly repeat photography as a quicker, cheaper and more adaptable tool to remotely study phenology in both biological and physical systems. Overall, increases in ambient temperatures coupled with decreases in winter ice and snow cover, and earlier occurrence of birch foliage, signal a reduction in the length of winter, a shift towards earlier springs and an increase in the length of available growing season in the Swedish sub-arctic.

Multi-temporal mapping of seagrass cover, species and biomass of the Eastern Banks, Moreton Bay, Australia, with links to shapefiles

The spatial and temporal dynamics of seagrasses have been studied from the leaf to patch (100 m**2) scales. However, landscape scale (> 100 km**2) seagrass population dynamics are unresolved in seagrass ecology. Previous remote sensing approaches have lacked the temporal or spatial resolution, or ecologically appropriate mapping, to fully address this issue. This paper presents a robust, semi-automated object-based image analysis approach for mapping dominant seagrass species, percentage cover and above ground biomass using a time series of field data and coincident high spatial resolution satellite imagery. The study area was a 142 km**2 shallow, clear water seagrass habitat (the Eastern Banks, Moreton Bay, Australia). Nine data sets acquired between 2004 and 2013 were used to create seagrass species and percentage cover maps through the integration of seagrass photo transect field data, and atmospherically and geometrically corrected high spatial resolution satellite image data (WorldView-2, IKONOS and Quickbird-2) using an object based image analysis approach. Biomass maps were derived using empirical models trained with in-situ above ground biomass data per seagrass species. Maps and summary plots identified inter- and intra-annual variation of seagrass species composition, percentage cover level and above ground biomass. The methods provide a rigorous approach for field and image data collection and pre-processing, a semi-automated approach to extract seagrass species and cover maps and assess accuracy, and the subsequent empirical modelling of seagrass biomass. The resultant maps provide a fundamental data set for understanding landscape scale seagrass dynamics in a shallow water environment. Our findings provide proof of concept for the use of time-series analysis of remotely sensed seagrass products for use in seagrass ecology and management.

Benthic foraminifera in Northwest Indian Ocean surface sediments

The relationship between the distribution of benthic foraminifera and sediment type and depositional environment in the Arabian Sea is discussed. The benthic foraminiferal fauna were sampled in nineteen Recent surface sediment samples, and geochemical variables of the sediment of the same samples were measured. The water depths for the box core samples varies from 440 to 4040 m. A total of 103 species and six species-complexes were identified. The geochemical properties were found to correspond well to the sediment type and depositional environment and six different sediment/depositional environment types could be distinguished. Analysis of the benthic foraminiferal fauna reveals specific faunal assemblages that are closely related to these sediment/depositional environment types.

Water temperature and salinity data from southern elephant seals from Marion Island

Over the last decade, several hundred seals have been equipped with conductivity-temperature-depth sensors in the Southern Ocean for both biological and physical oceanographic studies. A calibrated collection of seal-derived hydrographic data is now available, consisting of more than 165,000 profiles. The value of these hydrographic data within the existing Southern Ocean observing system is demonstrated herein by conducting two state estimation experiments, differing only in the use or not of seal data to constrain the system. Including seal-derived data substantially modifies the estimated surface mixedlayer properties and circulation patterns within and south of the Antarctic Circumpolar Current. Agreement with independent satellite observations of sea ice concentration is improved, especially along the East Antarctic shelf. Instrumented animals efficiently reduce a critical observational gap, and their contribution to monitoring polar climate variability will continue to grow as data accuracy and spatial coverage increase.

(Table 1) Carbon and isotopic composition of organic carbon in sediments from DSDP Leg 56 Holes

We have analyzed inorganic and organic carbons and determined the isotopic composition of both sedimentary organic carbon and inorganic carbon in carbonates contained in sediments recovered from Holes 434, 434A, 434B, 435, and 435A in the landward slope of Japan and from Hole 436 in the oceanic slope of the Japan Trench. Both inorganic and organic carbons were assayed at the P. P. Shirshov Institute of Oceanology, in the same sample, using the Knopp technique and measuring evolved CO2 gravimetrically. Each sample was analyzed twice in parallel. Measurements were of a ±0.05 per cent accuracy and a probability level of 0.95. Carbon isotopic analysis was carried out on a MI-1305 mass spectrometer at the I. M. Gubkin Institute of Petrochemical and Gas Industry and the results presented as dC13 values related to the PDB standard. The procedure for preparing samples for organic carbon isotopic analysis involved (1) drying damp sediments at 60°C; (2) treating samples, while heating, with 10 N HCl to remove carbonate carbon; and (3) evaporating surplus HCl at 60°C. The organic substance was turned to CO2 by oxidizing it in an oxygen atmosphere. To prepare samples for inorganic carbon isotopic analysis we decomposed the carbonates with orthophosphoric acid and refined the gas evolved. The dC13 measurements, including a full cycle of sample preparation, were of a ±0.5 per cent accuracy and a probability level of 0.95.

Circumarctic wetland dataset based on ENVISAT ASAR Global Monitoring Mode with links to geotiff image and NetCDF files (2005-01 to 2011-12)

A circumpolar representative and consistent wetland map is required for a range of applications ranging from upscaling of carbon fluxes and pools to climate modelling and wildlife habitat assessment. Currently available data sets lack sufficient accuracy and/or thematic detail in many regions of the Arctic. Synthetic aperture radar (SAR) data from satellites have already been shown to be suitable for wetland mapping. Envisat Advanced SAR (ASAR) provides global medium-resolution data which are examined with particular focus on spatial wetness patterns in this study. It was found that winter minimum backscatter values as well as their differences to summer minimum values reflect vegetation physiognomy units of certain wetness regimes. Low winter backscatter values are mostly found in areas vegetated by plant communities typically for wet regions in the tundra biome, due to low roughness and low volume scattering caused by the predominant vegetation. Summer to winter difference backscatter values, which in contrast to the winter values depend almost solely on soil moisture content, show expected higher values for wet regions. While the approach using difference values would seem more reasonable in order to delineate wetness patterns considering its direct link to soil moisture, it was found that a classification of winter minimum backscatter values is more applicable in tundra regions due to its better separability into wetness classes. Previous approaches for wetland detection have investigated the impact of liquid water in the soil on backscatter conditions. In this study the absence of liquid water is utilized. Owing to a lack of comparable regional to circumpolar data with respect to thematic detail, a potential wetland map cannot directly be validated; however, one might claim the validity of such a product by comparison with vegetation maps, which hold some information on the wetness status of certain classes. It was shown that the Envisat ASAR-derived classes are related to wetland classes of conventional vegetation maps, indicating its applicability; 30% of the land area north of the treeline was identified as wetland while conventional maps recorded 1-7%.

Towards convection-resolving, global atmospheric simulations with the Model for Prediction Across Scales (MPAS) v3.1: an extreme scaling experiment, links to sourcecode of model and grides

The Model for Prediction Across Scales (MPAS) is a novel set of Earth system simulation components and consists of an atmospheric model, an ocean model and a land-ice model. Its distinct features are the use of unstructured Voronoi meshes and C-grid discretisation to address shortcomings of global models on regular grids and the use of limited area models nested in a forcing data set, with respect to parallel scalability, numerical accuracy and physical consistency. This concept allows one to include the feedback of regional land use information on weather and climate at local and global scales in a consistent way, which is impossible to achieve with traditional limited area modelling approaches. Here, we present an in-depth evaluation of MPAS with regards to technical aspects of performing model runs and scalability for three medium-size meshes on four different high-performance computing (HPC) sites with different architectures and compilers. We uncover model limitations and identify new aspects for the model optimisation that are introduced by the use of unstructured Voronoi meshes. We further demonstrate the model performance of MPAS in terms of its capability to reproduce the dynamics of the West African monsoon (WAM) and its associated precipitation in a pilot study. Constrained by available computational resources, we compare 11-month runs for two meshes with observations and a reference simulation from the Weather Research and Forecasting (WRF) model. We show that MPAS can reproduce the atmospheric dynamics on global and local scales in this experiment, but identify a precipitation excess for the West African region. Finally, we conduct extreme scaling tests on a global 3?km mesh with more than 65 million horizontal grid cells on up to half a million cores. We discuss necessary modifications of the model code to improve its parallel performance in general and specific to the HPC environment. We confirm good scaling (70?% parallel efficiency or better) of the MPAS model and provide numbers on the computational requirements for experiments with the 3?km mesh. In doing so, we show that global, convection-resolving atmospheric simulations with MPAS are within reach of current and next generations of high-end computing facilities.