Field data, numerical simulations and probability analyses to assess lava flow hazards at Mount Etna

Improving lava flow hazard assessment is one of the most important and challenging fields of volcanology, and has an immediate and practical impact on society. Here, we present a methodology for the quantitative assessment of lava flow hazards based on a combination of field data, numerical simulations and probability analyses. With the extensive data available on historic eruptions of Mt. Etna, going back over 2000 years, it has been possible to construct two hazard maps, one for flank and the other for summit eruptions, allowing a quantitative analysis of the most likely future courses of lava flows. The effective use of hazard maps of Etna may help in minimizing the damage from volcanic eruptions through correct land use in densely urbanized area with a population of almost one million people. Although this study was conducted on Mt. Etna, the approach used is designed to be applicable to other volcanic areas.

The SUGAR Toolbox

The SUGAR Toolbox contains scripts coded in MATLAB for calculating various thermodynamic, kinetic, and geologic properties of substances occurring in the marine environment, particularly gas hydrate and seep systems. Brief descriptions of the toolbox scripts and some notes on the underlying basic theory as well as tables of additional property values can be found in the accompanying documentation.

The SUGAR Toolbox - Version 7

The SUGAR Toolbox contains scripts coded in MATLAB for calculating various thermodynamic, kinetic, and geologic properties of substances occurring in the marine environment, particularly gas hydrate and seep systems. Brief descriptions of the toolbox scripts and some notes on the underlying basic theory as well as tables of additional property values can be found in the accompanying documentation.

Stable carbon isotope in bottom water and Cibicidoides wuellerstorfi on the Ontong Java Plateau and Emperor Seamount

We have measured the carbon isotopic composition of dissolved inorganic carbon in bottom waters of the Ontong Java Plateau (western equatorial Pacific) and on the northern Emperor Seamounts (northwest Pacific). Each of these locations is several hundred miles from the nearest Geochemical Ocean Sections Study (GEOSECS) stations, and the observed delta13C values at each site differ substantially from regionally averaged GEOSECS delta13C profiles. We discuss the possible causes of these differences, including horizontal variability, near-bottom effects, and problems with the Pacific GEOSECS delta13C data. We also measured the isotopic composition (C and O) of core top C. wuellerstorfi from a depth transect of cores at each location. The delta18O data are used to verify that our samples are Holocene. Comparison of foraminiferal and bottom water delta13C values shows that this species faithfully records bottom water delta13C at both sites and demonstrates that there is no depth-related artifact in the dissolved inorganic carbon-C. wuellerstorfi delta13C relationship at these sites.

Grain entrainment and transport rates of magnetic minerals of coastal sediment samples from the Bay of Plenty, New Zealand

Heavy (magnetic & non-magnetic) minerals are found concentrated by natural processes in many fluvial, estuarine, coastal and shelf environments with a potential to form economic placer deposits. Understanding the processes of heavy mineral transport and enrichment is prerequisite to interpret sediment magnetic properties in terms of hydro- and sediment dynamics. In this study, we combine rock magnetic and sedimentological laboratory measurements with numerical 3D discrete element models to investigate differential grain entrainment and transport rates of magnetic minerals in a range of coastal environments (riverbed, mouth, estuary, beach and near-shore). We analyzed grain-size distributions of representative bulk samples and their magnetic mineral fractions to relate grain-size modes to respective transport modes (traction, saltation, suspension). Rock magnetic measurements showed that distribution shapes, population sizes and grain-size offsets of bulk and magnetic mineral fractions hold information on the transport conditions and enrichment process in each depositional environment. A downstream decrease in magnetite grain size and an increase in magnetite concentration was observed from riverine source to marine sink environments. Lower flow velocities permit differential settling of light and heavy mineral grains creating heavy mineral enriched zones in estuary settings, while lighter minerals are washed out further into the sea. Numerical model results showed that higher heavy mineral concentrations in the bed increased the erosion rate and enhancing heavy mineral enrichment. In beach environments where sediments contained light and heavy mineral grains of equivalent grain sizes, the bed was found to be more stable with negligible amount of erosion compared to other bed compositions. Heavy mineral transport rates calculated for four different bed compositions showed that increasing heavy mineral content in the bed decreased the transport rate. There is always a lag in transport between light and heavy minerals which increases with higher heavy mineral concentration in all tested bed compositions. The results of laboratory experiments were validated by numerical models and showed good agreement. We demonstrate that the presented approach bears the potential to investigate heavy mineral enrichment processes in a wide range of sedimentary settings.

Composition and age of volcanic rocks from the Sinii Utes Depression, Southern Primoye

New data are reported on structure of sections, chemical composition, and age of volcano-sedimentary and volcanic rocks from the Sinii Utes Depression in the Southern Primorye region. The Sinii Utes Depression is filled with two sequences: the lower sequence composed of sedimentary-volcanogenic coaliferous rocks (the stratotype of the Sinii Utes Formation) and the upper sequence consisting of tephroid with overlying basalts. This work considers chemical composition and problems of K-Ar dating of basalts. The uppermost basaltic flow has K-Ar age 22.0±1.0 Ma. The dates obtained for the middle and upper parts of lava flows are underestimated. It is explained by their heating due to combustion of brown coals of the Sinii Utes Formation underlying the lava flow. Calculations show that argon could only partly have been removed from the basalts owing to conductive heat transfer and was lost largely due to infiltration of hot gases in heterogeneous fissured medium. Basaltic volcanism on continental margins of the southern Primorye region and the adjacent Korean and Chinese areas at the Oligocene-Miocene boundary preceded Early-Middle Miocene spreading and formation of the Sea of Japan basin. Undifferentiated moderately alkaline basalts of intraplate affinity developed in the Amba Depression and some other structures of the southern Primorye region and intraplate alkali basalts of the Phohang Graben in the Korean Peninsula serve as indicators of incipient spreading regime in the Sea of Japan. Potassic basalt-trachybasalt eruptions occurred locally in riftogenic depressions and shield volcanoes. In some structures this volcanism was terminated by eruptions of intermediate and acid lavas. Such evolution of volcanism is explained by selective contamination of basaltic melts during their interaction with crustal acid material and generation of acid anatectic melts.

Equivalent dose, dose rate and age determination of cobble surfaces in the Antarctic

Most current methods of reconstructing past sea levels within Antarctica rely on radiocarbon dating. However, radiocarbon dating is limited by the availability of material for dating and problems inherent with radiocarbon reservoirs in Antarctic marine systems. Here we report on the success of a new approach to dating raised beach deposits in Antarctica for the purpose of reconstructing past sea levels. This new approach is the use of optically stimulated luminescence (OSL) on quartz-grains obtained from the underside of cobbles within raised beaches and boulder pavements. We obtained eight OSL dates from three sites along the shores of Maxwell Bay in the South Shetland Islands of the Antarctic Peninsula. These dates are internally consistent and fit well with previously published radiocarbon ages obtained from the same deposits. In addition, when the technique was applied to a modern beach, it resulted in an age of zero. Our results suggest that this method will provide a valuable tool in the reconstruction of past sea levels in Antarctica and other coarse-grained beach deposits across the globe.

GPS raw data (control points and ground control points) from the Liguria Region, Borghetto Santo Spirito, Italy

Monitoring the impact of sea storms on coastal areas is fundamental to study beach evolution and the vulnerability of low-lying coasts to erosion and flooding. Modelling wave runup on a beach is possible, but it requires accurate topographic data and model tuning, that can be done comparing observed and modeled runup. In this study we collected aerial photos using an Unmanned Aerial Vehicle after two different swells on the same study area. We merged the point cloud obtained with photogrammetry with multibeam data, in order to obtain a complete beach topography. Then, on each set of rectified and georeferenced UAV orthophotos, we identified the maximum wave runup for both events recognizing the wet area left by the waves. We then used our topography and numerical models to simulate the wave runup and compare the model results to observed values during the two events. Our results highlight the potential of the methodology presented, which integrates UAV platforms, photogrammetry and Geographic Information Systems to provide faster and cheaper information on beach topography and geomorphology compared with traditional techniques without losing in accuracy. We use the results obtained from this technique as a topographic base for a model that calculates runup for the two swells. The observed and modeled runups are consistent, and open new directions for future research.

(Table 1) Ice flux, polynya area and sea-ice drift speed in the Kara Sea during winter (Jan-Apr) 1997-2001

The Shelf Seas of the Arctic are known for their large sea-ice production. This paper presents a comprehensive view of the Kara Sea sea-ice cover from high-resolution numerical modeling and space-borne microwave radiometry. As given by the latter the average polynya area in the Kara Sea takes a value of 21.2 × 10**3 km**2 ± 9.1 × 10**3 km**2 for winters (Jan.-Apr.) 1996/97 to 2000/01, being as high as 32.0 × 10**3 km**2 in 1999/2000 and below 12 × 10**3 km**2 in 1998/99. Day-to-day variations of the Kara Sea polynya area can be as high as 50 × 10**3 km**2. For the seasons 1996/97 to 2000/01 the modeled cumulative winter ice-volume flux out of the Kara Sea varied between 100 km**3/a and 350 km**3/a. Modeled high (low) ice export coincides with a high (low) average and cumulative polynya area, and with a low (high) sea-ice compactness in the Kara Sea from remote sensing data, and with a high (low) sea-ice drift speed across its northern boundary derived from independent model data for the winters 1996/97 to 2000/01.

The CoastColour Round Robin datasets: a database to evaluate the performance of algorithms for the retrieval of water quality parameters in coastal waters

The CoastColour project Round Robin (CCRR) project (http://www.coastcolour.org) funded by the European Space Agency (ESA) was designed to bring together a variety of reference datasets and to use these to test algorithms and assess their accuracy for retrieving water quality parameters. This information was then developed to help end-users of remote sensing products to select the most accurate algorithms for their coastal region. To facilitate this, an inter-comparison of the performance of algorithms for the retrieval of in-water properties over coastal waters was carried out. The comparison used three types of datasets on which ocean colour algorithms were tested. The description and comparison of the three datasets are the focus of this paper, and include the Medium Resolution Imaging Spectrometer (MERIS) Level 2 match-ups, in situ reflectance measurements and data generated by a radiative transfer model (HydroLight). The datasets mainly consisted of 6,484 marine reflectance associated with various geometrical (sensor viewing and solar angles) and sky conditions and water constituents: Total Suspended Matter (TSM) and Chlorophyll-a (CHL) concentrations, and the absorption of Coloured Dissolved Organic Matter (CDOM). Inherent optical properties were also provided in the simulated datasets (5,000 simulations) and from 3,054 match-up locations. The distributions of reflectance at selected MERIS bands and band ratios, CHL and TSM as a function of reflectance, from the three datasets are compared. Match-up and in situ sites where deviations occur are identified. The distribution of the three reflectance datasets are also compared to the simulated and in situ reflectances used previously by the International Ocean Colour Coordinating Group (IOCCG, 2006) for algorithm testing, showing a clear extension of the CCRR data which covers more turbid waters.

Geochemistry and minerals at DSDP Legs 68-70 Holes

The compositions of natural glasses and phenocrysts in basalts from Deep Sea Drilling Project Sites 501, 504, and 505, near the Costa Rica Rift, constitute evidence for the existence of a periodically replenished axial magma chamber that repeatedly erupted lavas of remarkably uniform composition. Magma compositions were affected by three general components: (1) injected magmas carrying (in decreasing order of abundance) Plagioclase, olivine, and chrome-spinel phenocrysts (spinel assemblage); (2) injected magmas carrying Plagioclase, clinopyroxene, and olivine phenocrysts, but no spinel (clinopyroxene assemblage); and (3) moderately evolved hybrids in the magma chamber itself. The compositions of the injected phenocrysts and minerals in glomerocrysts are as follows: Plagioclase - An85-94; olivine - Fo87-89; clinopyroxene - high Cr2O3 (0.7-1.1%), endiopside (Wo42En51Fs7), and aluminous chromian spinel (Cr/Cr + Al = 0.3). These minerals resemble those thought to occur in upper mantle sources (9 kbars and less) of ocean-ridge basalts and to crystallize in magmas near those sources. In the magma chamber, more sodic Plagioclase (An79-85), less magnesian olivine (Fo81-86) and low-Cr2O3 (0.1-0.4%) clinopyroxene formed rims on these crystals, grew as other phenocrysts, and formed cumulus segregations on the walls and floors of the magma chamber. In the spinel-assemblage magmas, magnesiochromite (Cr/Cr + Al = 0.4-0.5) also formed. Some cumulus segregations were later entrained in lavas as xenoliths. The glass compositions define 16 internally homogeneous eruptive units, 13 of which are in stratigraphic order in a single hole, Hole 504B, which was drilled 561.5 meters into the ocean crust. These units are defined as differing from each other by more than analytical uncertainty in one or more oxides. However, many of the glass groups in Hole 504B show virtually no differences in TiO2 contents, Mg/Mg + Fe2+, or normative An/An + Ab, all of which are sensitive indicators of crystallization differentiation. The differences are so small that they are only apparent in the glass compositions; they are almost completely obscured in whole-rock samples by the presence of phenocrysts and the effects of alteration. Moreover, several of the glass units at different depths in Hole 504B are compositionally identical, with all oxides falling within the range of analytical uncertainty, with only small variations in the rest of the suite. The repetition of identical chemical types requires (1) very regular injection of magmas into the magma chamber, (2) extreme similarity of injected magmas, and (3) displacement of very nearly the same proportion of the magmas in the chamber at each injection. Numerical modeling and thermal considerations have led some workers to propose the existence of such conditions at certain types of spreading centers, but the lava and glass compositions at Hole 504B represent the first direct evidence revealed by drilling of the existence of a compositionally nearly steady-state magma chamber, and this chapter examines the processes acting in it in some detail. The glass groups that are most similar are from clinopyroxene-assemblage lavas, which have a range of Mg/Mg + Fe2"1" of 0.59 to 0.65. Spinel-assemblage basalts are less evolved, with Mg/Mg + Fe2+ of 0.65 to 0.69, but both types have nearly identical normative An/An + Ab (0.65-0.66). However, the two lava types contain megacrysts (olivine, Plagioclase, clinopyroxene) that crystallized from melts with Mg/Mg + Fe2+ values of 0.70 to 0.72. Projection of glass compositions into ternary normative systems suggests that spinel-assemblage magmas originated deeper in the mantle than clinopyroxene-assemblage magmas, and mineral data indicate that the two types followed different fractionation paths before reaching the magma chamber. The two magma types therefore represent neither a low- nor a high-pressure fractionation sequence. Some of the spinel-assemblage magmas may have had picritic parents, but were coprecipitating all of the spinel-assemblage phenocrysts before reaching the magma chamber. Clinopyroxene-assemblage magmas did not have picritic parents, but the compositions of phenocrysts suggest that they originated at about 9 kbars, near the transition between plagioclase peridotite and spinel peridotite in the mantle. Two glass groups have higher contents of alkalis, TiO2, and P2O5 than the others, evidently as a result of the compositions of mantle sources. Eruption of these lavas implies that conduits and chambers containing magmas from dissimilar sources were not completely interconnected on the Costa Rica Rift. The data are used to draw comparisons with the East Pacific Rise and to consider the mechanisms that may have prevented the eruption of ferrobasalts at these sites.

Bulk density, ash content, humification, peat accumulation and selected lithogenic data of Store Mosse peat bog

This datafile presents chemical and physical as well as age dating information from the Store Mosse peat bog in southern Sweden. This record dates back to 8900 cal yr BP. The aim of the research was to reconstruct mineral dust deposition over time. As such we have only presented the lithogenic element data (Al, Ga, Rb, Sc, Ti, Y, Zr, Th and the REE) as the sample preparation method was tailored to these. This data is supported by parameters describing the deposit including bulk density, humification, ash content and net peat accumulation rates.

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.

Biogenic opal, biogenic barium, calcium carbonate and U/Th isotope data from core-top sediments of SONNE cruise SO202 (INOPEX)

Sedimentary proxies used to reconstruct marine productivity suffer from variable preservation and are sensitive to factors other than productivity. Therefore, proxy calibration is warranted. Here we map the spatial patterns of two paleoproductivity proxies, biogenic opal and barium fluxes, from a set of core-top sediments recovered in the Subarctic North Pacific. Comparisons of the proxy data with independent estimates of primary and export production, surface water macronutrient concentrations and biological pCO2 drawdown indicate that neither proxy shows a significant correlation with primary or export productivity for the entire region. Biogenic opal fluxes, when corrected for preservation using 230Th-normalized accumulation rates, show a good correlation with primary productivity along the volcanic arcs (tau = 0.71, p = 0.0024) and with export productivity throughout the western Subarctic North Pacific (tau = 0.71, p = 0.0107). Moderate and good correlations of biogenic barium flux with export production (tau = 0.57, p = 0.0022) and with surface water silicate concentrations (tau = 0.70, p = 0.0002) are observed for the central and eastern Subarctic North Pacific. For reasons unknown, however, no correlation is found in the western Subarctic North Pacific between biogenic barium flux and the reference data. Nonetheless, we show that barite saturation, uncertainty in the lithogenic barium corrections and problems with the reference datasets are not responsible for the lack of a significant correlation between biogenic barium flux and the reference data. Further studies evaluating the factors controlling the variability of the biogenic constituents in the sediments are desirable in this region.