(Table 1) Age determination of basalts from DSDP Hole 30-289

The Pacific plate has undergone a substantial northward displacement during the late Mesozoic and the Cainozoic. Here we give additional documentation for such motion based on palaeomagnetic measurements of a sequence of sedimentary and basalt samples collected from middle Oligocene to Aptian sections of Deep Sea Drilling Project (DSDP) site 289 (Andrews, 1975; 00° 29.92'S, 158° 30.69'E) drilled on the Ontong Java Plateau.

EPANET Input Files of New York tunnels and Pacific City used in a metamodel-based optimization study

Metamodels have proven be very useful when it comes to reducing the computational requirements of Evolutionary Algorithm-based optimization by acting as quick-solving surrogates for slow-solving fitness functions. The relationship between metamodel scope and objective function varies between applications, that is, in some cases the metamodel acts as a surrogate for the whole fitness function, whereas in other cases it replaces only a component of the fitness function. This paper presents a formalized qualitative process to evaluate a fitness function to determine the most suitable metamodel scope so as to increase the likelihood of calibrating a high-fidelity metamodel and hence obtain good optimization results in a reasonable amount of time. The process is applied to the risk-based optimization of water distribution systems; a very computationally-intensive problem for real-world systems. The process is validated with a simple case study (modified New York Tunnels) and the power of metamodelling is demonstrated on a real-world case study (Pacific City) with a computational speed-up of several orders of magnitude.

Atmospheric, oceanic and hydrological polar motion excitation mechanisms based on a combination of geometric and gravimetric space observations

The goal of our study is to determine accurate time series of geophysical Earth rotation excitations to learn more about global dynamic processes in the Earth system. For this purpose, we developed an adjustment model which allows to combine precise observations from space geodetic observation systems, such as Satellite Laser Ranging (SLR), Global Navigation Satellite Systems (GNSS), Very Long Baseline Interferometry (VLBI), Doppler Orbit determination and Radiopositioning Integrated on Satellite (DORIS), satellite altimetry and satellite gravimetry in order to separate geophysical excitation mechanisms of Earth rotation. Three polar motion time series are applied to derive the polar motion excitation functions (integral effect). Furthermore we use five time variable gravity field solutions from Gravity Recovery and Climate Experiment (GRACE) to determine not only the integral mass effect but also the oceanic and hydrological mass effects by applying suitable filter techniques and a land-ocean mask. For comparison the integral mass effect is also derived from degree 2 potential coefficients that are estimated from SLR observations. The oceanic mass effect is also determined from sea level anomalies observed by satellite altimetry by reducing the steric sea level anomalies derived from temperature and salinity fields of the oceans. Due to the combination of all geodetic estimated excitations the weaknesses of the individual processing strategies can be reduced and the technique-specific strengths can be accounted for. The formal errors of the adjusted geodetic solutions are smaller than the RMS differences of the geophysical model solutions. The improved excitation time series can be used to improve the geophysical modeling.

Terrigenous input and marine productivity for the Pacific sector of the Southern Ocean based on lipid biomarkers

In this study, we present a new multiproxy data set of terrigenous input, marine productivity and sea surface temperature (SST) from 52 surface sediment samples collected along E-W transects in the Pacific sector of the Southern Ocean. Allochtonous terrigenous input was characterized by the distribution of plant wax n-alkanes and soil-derived branched glycerol dialkyl glycerol tetraethers (brGDGTs). 230Th-normalized burial rates of both compound groups were highest close to the potential sources in Australia and New Zealand and are strongly related to lithogenic contents, indicating common sources and transport. Detection of both long-chain n-alkanes and brGDGTs at the most remote sites in the open ocean strongly suggests a primarily eolian transport mechanism to at least 110°W, i.e. by prevailing westerly winds. Two independent organic SST proxies were used, the UK'37 based on long-chain alkenones, and the TEX86 based on isoprenoid GDGTs. Both, UK'37 and TEX86 indices show robust relationships with temperature over a temperature range between 0.5 and 20°C, likely implying different seasonal and regional imprints on the temperature signal. While alkenone-based temperature estimates reliably reflect modern SST even at the low temperature end, large temperature residuals are observed for the polar ocean using the TEX86 index. 230Th-normalized burial rates of alkenones are highest close to the Subtropical Front and are positively related to lithogenic fluxes throughout the study area. In contrast, highest isoGDGT burial south of the Antarctic Polar Front is not related with dust flux but may be largely controlled by diatom blooms, and thus high opal fluxes during austral summer.

Environmental variables and multivariate seabed classification maps via k means clustering of Potter Cove, Antarctica, link to input and result files in GeoTIFF format

This study subdivides the Potter Cove, King George Island, Antarctica, into seafloor regions using multivariate statistical methods. These regions are categories used for comparing, contrasting and quantifying biogeochemical processes and biodiversity between ocean regions geographically but also regions under development within the scope of global change. The division obtained is characterized by the dominating components and interpreted in terms of ruling environmental conditions. The analysis includes in total 42 different environmental variables, interpolated based on samples taken during Australian summer seasons 2010/2011 and 2011/2012. The statistical errors of several interpolation methods (e.g. IDW, Indicator, Ordinary and Co-Kriging) with changing settings have been compared and the most reasonable method has been applied. The multivariate mathematical procedures used are regionalized classification via k means cluster analysis, canonical-correlation analysis and multidimensional scaling. Canonical-correlation analysis identifies the influencing factors in the different parts of the cove. Several methods for the identification of the optimum number of clusters have been tested and 4, 7, 10 as well as 12 were identified as reasonable numbers for clustering the Potter Cove. Especially the results of 10 and 12 clusters identify marine-influenced regions which can be clearly separated from those determined by the geological catchment area and the ones dominated by river discharge.

Code and readme of a MATLAB-based graphical user interface program for computing functionals of the geopotential

We present a novel graphical user interface program GrafLab (GRAvity Field LABoratory) for spherical harmonic synthesis (SHS) created in MATLAB®. This program allows to comfortably compute 38 various functionals of the geopotential up to ultra-high degrees and orders of spherical harmonic expansion. For the most difficult part of the SHS, namely the evaluation of the fully normalized associated Legendre functions (fnALFs), we used three different approaches according to required maximum degree: (i) the standard forward column method (up to maximum degree 1800, in some cases up to degree 2190); (ii) the modified forward column method combined with Horner's scheme (up to maximum degree 2700); (iii) the extended-range arithmetic (up to an arbitrary maximum degree). For the maximum degree 2190, the SHS with fnALFs evaluated using the extended-range arithmetic approach takes only approximately 2-3 times longer than its standard arithmetic counterpart, i.e. the standard forward column method. In the GrafLab, the functionals of the geopotential can be evaluated on a regular grid or point-wise, while the input coordinates can either be read from a data file or entered manually. For the computation on a regular grid we decided to apply the lumped coefficients approach due to significant time-efficiency of this method. Furthermore, if a full variance-covariance matrix of spherical harmonic coefficients is available, it is possible to compute the commission errors of the functionals. When computing on a regular grid, the output functionals or their commission errors may be depicted on a map using automatically selected cartographic projection.

Determination of millennial-scale climate variability based on alkenones in the western Atlantic Ocean

High-resolution records of alkenone-derived sea surface temperatures and elemental Ti/Ca ratios from a sediment core retrieved off northeastern Brazil (4° S) reveal short-term climate variability throughout the past 63,000 a. Large pulses of terrigenous sediment discharge, caused by increased precipitation in the Brazilian hinterland, coincide with Heinrich events and the Younger Dryas period. Terrigenous input maxima related to Heinrich events H6-H2 are characterized by rapid cooling of surface water ranging between 0.5 and 2° C. This signature is consistent with a climate model experiment where a reduction of the Atlantic meridional overturning circulation (AMOC) and related North Atlantic cooling causes intensification of NE trade winds and a southward movement of the Intertropical Convergence Zone, resulting in enhanced precipitation off northeastern Brazil. During deglaciation the surface temperature evolution at the core site predominantly followed the Antarctic warming trend, including a cooling, prior to the Younger Dryas period. An abrupt temperature rise preceding the onset of the Bølling/Allerød transition agrees with model experiments suggesting a Southern Hemisphere origin for the abrupt resumption of the AMOC during deglaciation caused by Southern Ocean warming and associated with northward flow anomalies of the South Atlantic western boundary current.

(Figures S2a-c) Phase diagrams of the upper crust, the lower crust and the upper mantle based on averaged compositions of the oceanic crust and mantle

The age of the subducting Nazca Plate off Chile increases northwards from 0 Ma at the Chile Triple Junction (46°S) to 37 Ma at the latitude of Valparaíso (32°S). Age-related variations in the thermal state of the subducting plate impact on (a) the water influx to the subduction zone, as well as on (b) the volumes of water that are released under the continental forearc or, alternatively, carried beyond the arc. Southern Central Chile is an ideal setting to study this effect, because other factors for the subduction zone water budget appear constant. We determine the water influx by calculating the crustal water uptake and by modeling the upper mantle serpentinization at the outer rise of the Chile Trench. The water release under forearc and arc is determined by coupling FEM thermal models of the subducting plate with stability fields of water-releasing mineral reactions for upper and lower crust and hydrated mantle. Results show that both the influx of water stored in, and the outflux of water released from upper crust, lower crust and mantle vary drastically over segment boundaries. In particular, the oldest and coldest segments carry roughly twice as much water into the subduction zone as the youngest and hottest segments, but their release flux to the forearc is only about one fourth of the latter. This high variability over a subduction zone of < 1500 km length shows that it is insufficient to consider subduction zones as uniform entities in global estimates of subduction zone fluxes. This article is protected by copyright. All rights reserved.

Abundance and opal input of radiolarian species of two sediment cores from the Atlantic sector of the Southern Ocean

Radiolarian-based paleoceanographic reconstructions generally use the abundance of selected radiolarian species. However, the recent focus on the opal flux and the development of isotope measurements in biogenic opal and the organic matter embedded in it demands a better knowledge of the origin of the opal. We present here an estimation of the opal content of the skeleton of 63 radiolarian species from two sites in the Southern Ocean. The skeletons are modelled as associations of simple geometrical shapes, and the volume thus obtained is combined with opal density to obtain the amount of opal. These data are, thus, used to determine the most important opal carriers in the radiolarian assemblage in both cores.

Modes of sapropel formation in the eastern Mediterranean: some constraints based on pyrite properties

Pyrite formation within and directly below sapropels in the eastern Mediterranean was governed by the relative rates of sulphide production and Fe liberation and supply to the organic-rich layers. At times of relatively high [SO4]2- reduction, sulphide could diffuse downward from the sapropel and formed pyrite in underlying sediments. The sources of Fe for pyrite formation comprised detrital Fe and diagenetically liberated Fe(II) from sapropel-underlying sediments. In organic-rich sapropels, input of Fe from the water column via Fe sulphide formation in the water may have been important as well. Rapid pyrite formation at high saturation levels resulted in the formation of framboidal pyrite within the sapropels, whereas below the sapropels slow euhedral pyrite formation at low saturation levels occurred. d34S values of pyrite are -33 per mil to -50 per mil. Below the sapropels d34S is lower than within the sapropels, as a result of increased sulphide re-oxidation at times of relatively high sulphide production and concentration when sulphide could escape from the sediment. The percentage of initially formed sulphide that was re-oxidized was estimated from organic carbon fluxes and burial efficiencies in the sediment. It ranges from 34% to 80%, varying significantly between sapropels. Increased palaeoproductivity as well as enhanced preservation contributed to magnified accumulation of organic matter in sapropels.