Regional geoid model for the area of the subglacial Lake Vostok

We present a regional geoid model for the area of Lake Vostok, Antarctica, from a combination of local airborne gravity, ice-surface and ice-thickness data and a lake bathymetry model. The topography data are used for residual terrain modelling (RTM) in a remove-compute-restore approach together with the GOCE satellite model GOCO03S. The disturbing potential at the Earth's surface, i.e. the quasigeoid, is predicted by least-squares collocation (LSC) and subsequently converted to geoid heights. Compared to GOCO03S our regional solution provides an additional short-wavelength signal of up to 1.48 m, or 0.56 m standard deviation, respectively. More details can be found in Schwabe et. al (2014).

Pollen record and calculated holocene vegetation and climate dynamics from Sihailongwan Maar Lake, NE China

High-resolution palynological analysis on annually laminated sediments of Sihailongwan Maar Lake (SHL) provides new insights into the Holocene vegetation and climate dynamics of NE China. The robust chronology of the presented record is based on varve counting and AMS radiocarbon dates from terrestrial plant macro-remains. In addition to the qualitative interpretation of the pollen data, we provide quantitative reconstructions of vegetation and climate based on the method of biomization and weighted averaging partial least squares regression (WA-PLS) technique, respectively. Power spectra were computed to investigate the frequency domain distribution of proxy signals and potential natural periodicities. Pollen assemblages, pollen-derived biome scores and climate variables as well as the cyclicity pattern indicate that NE China experienced significant changes in temperature and moisture conditions during the Holocene. Within the earliest phase of the Holocene, a large-scale reorganization of vegetation occurred, reflecting the reconstructed shift towards higher temperatures and precipitation values and the initial Holocene strengthening and northward expansion of the East Asian summer monsoon (EASM). Afterwards, summer temperatures remain at a high level, whereas the reconstructed precipitation shows an increasing trend until approximately 4000 cal. yr BP. Since 3500 cal. yr BP, temperature and precipitation values decline, indicating moderate cooling and weakening of the EASM. A distinct periodicity of 550-600 years and evidence of a Mid-Holocene transition from a temperature-triggered to a predominantly moisture-triggered climate regime are derived from the power spectra analysis. The results obtained from SHL are largely consistent with other palaeoenvironmental records from NE China, substantiating the regional nature of the reconstructed vegetation and climate patterns. However, the reconstructed climate changes contrast with the moisture evolution recorded in S China and the mid-latitude (semi-)arid regions of N China. Whereas a clear insolation-related trend of monsoon intensity over the Holocene is lacking from the SHL record, variations in the coupled atmosphere-Pacific Ocean system can largely explain the reconstructed changes in NE China.

Paleomagnetic of lower Cretaceous sediments from the Argo and Gascoyne Abyssal Plain

Lower Cretaceous sediments were sampled for magnetostratigraphy at three sites. ODP Site 765 and DSDP Site 261, in the Argo Abyssal Plain, consist primarily of brownish-red to gray claystone having hematite and magnetite carriers of characteristic magnetization. ODP Site 766, in the Gascoyne Abyssal Plain, consists mainly of dark greenish-gray volcaniclastic turbidites with magnetite as the carrier of characteristic magnetization. Progressive thermal demagnetization (Sites 765 and 261) or alternating field demagnetization (Site 766) yielded well-defined polarity zones and a set of reliable paleolatitudes. Magnetic polarity chrons were assigned to polarity zones using biostratigraphic correlations. Late Aptian chron M"-1"r, a brief reversed-polarity chron younger than MOr, is a narrow, 40-cm feature delimited in Hole 765C. Early Aptian reversed-polarity chron MOr is also present in Hole 765C. Polarity chrons Mir through M3r were observed in the Barremian of all three sites. Valanginian and Hauterivian polarity chrons can be tentatively assigned to polarity zones only in Hole 766A. The paleolatitude of this region remained at 35° to 37°S from the Berriasian through the Aptian. During this interval, there was approximately 16° of clockwise rotation, with the oriented sample suites of Site 765 displaying a Berriasian declination of 307° to an Aptian declination of 323°. These results are consistent with the interpolated Early Cretaceous apparent polar wander for Australia, but indicate that this region was approximately 5? farther north than predicted.

Paleomagneitc of early Cretaceous sediments from the western Central Atlantic

Drilling at Sites 534 and 603 of the Deep Sea Drilling Project recovered thick sections of Berriasian through Aptian white limestones to dark gray marls, interbedded with claystone and clastic turbidites. Progressive thermal demagnetization removed a normal-polarity overprint carried by goethite and/or pyrrhotite. The resulting characteristic magnetization is carried predominantly by magnetite. Directions and reliability of characteristic magnetization of each sample were computed by using least squares line-fits of magnetization vectors. The corrected true mean inclinations of the sites suggest that the western North Atlantic underwent approximately 6° of steady southward motion between the Berriasian and Aptian stages. The patterns of magnetic polarity of the two sites, when plotted on stratigraphic columns of the pelagic sediments without turbidite beds, display a fairly consistent magnetostratigraphy through most of the Hauterivian-Barremian interval, using dinoflagellate and nannofossil events and facies changes in pelagic sediment as controls on the correlations. The composite magnetostratigraphy appears to include most of the features of the M-sequence block model of magnetic anomalies from Ml to Ml ON (Barremian-Hauterivian) and from M16 to M23 (Berriasian-Tithonian). The Valanginian magnetostratigraphy of the sites does not exhibit reversed polarity intervals corresponding to Ml 1 to M13 of the M-sequence model; this may be the result of poor magnetization, of a major unrecognized hiatus in the early to middle Valanginian in the western North Atlantic, or of an error in the standard block model. Based on these tentative polarity-zone correlations, the Hauterivian/Barremian boundary occurs in or near the reversed-polarity Chron M7 or M5, depending upon whether the dinoflagellate or nannofossil zonation, respectively, is used; the Valanginian/Hauterivian boundary, as defined by the dinoflagellate zonation, is near reversed-polarity Chron M10N.

Basalt density, basement age, and intrusive/extrusive relations from DSDP Legs 2 through 7, 9, and 14

Densities of layer 2 basalt recovered during the Deep Sea Drilling Project have been found to decrease steadily with age, a finding ascribed to progressive submarine weathering in the context of sea-floor spreading. The least-squares solution for 52 density measurements gives a rate of decrease in density of (Delta p)/(Delta t) = -0.0046 g per ccm m.y. = -16 percent per 100 m.y., which is in excellent agreement with earlier estimates based on observed chemical depletion rates of dredged oceanic basalt. Weathering of sea-floor basalt, should it penetrate to any considerable depth in layer 2, will decrease layer 2 seismic refraction velocities, act as a source of geothermal heat, and substantially influence the chemistry of sea water and the overlying column of sediment.

Physical properties and sedimentology on core GeoB1510-1

Detailed information about the sediment properties and microstructure can be provided through the analysis of digital ultrasonic P wave seismograms recorded automatically during full waveform core logging. The physical parameter which predominantly affects the elastic wave propagation in water-saturated sediments is the P wave attenuation coefficient. The related sedimentological parameter is the grain size distribution. A set of high-resolution ultrasonic transmission seismograms (ca. 50-500 kHz), which indicate downcore variations in the grain size by their signal shape and frequency content, are presented. Layers of coarse-grained foraminiferal ooze can be identified by highly attenuated P waves, whereas almost unattenuated waves are recorded in fine-grained areas of nannofossil ooze. Color-encoded pixel graphics of the seismograms and instantaneous frequencies present full waveform images of the lithology and attenuation. A modified spectral difference method is introduced to determine the attenuation coefficient and its power law a = kfn. Applied to synthetic seismograms derived using a "constant Q" model, even low attenuation coefficients can be quantified. A downcore analysis gives an attenuation log which ranges from ca. 700 dB/m at 400 kHz and a power of n = 1-2 in coarse-grained sands to few decibels per meter and n ? 0.5 in fine-grained clays. A least squares fit of a second degree polynomial describes the mutual relationship between the mean grain size and the attenuation coefficient. When it is used to predict the mean grain size, an almost perfect coincidence with the values derived from sedimentological measurements is achieved.

Improved geoid solution for the Weddell Sea region

We present a geoid solution for the Weddell Sea and adjacent continental Antarctic regions. There, a refined geoid is of interest, especially for oceanographic and glaciological applications. For example, to investigate the Weddell Gyre as a part of the Antarctic Circumpolar Current and, thus, of the global ocean circulation, the mean dynamic topography (MDT) is needed. These days, the marine gravity field can be inferred with high and homogeneous resolution from altimetric height profiles of the mean sea surface. However, in areas permanently covered by sea ice as well as in coastal regions, satellite altimetry features deficiencies. Focussing on the Weddell Sea, these aspects are investigated in detail. In these areas, ground-based data that have not been used for geoid computation so far provide additional information in comparison with the existing high-resolution global gravity field models such as EGM2008. The geoid computation is based on the remove-compute-restore approach making use of least-squares collocation. The residual geoid with respect to a release 4 GOCE model adds up to two meters and more in the near-coastal and continental areas of the Weddell Sea region, also in comparison with EGM2008. Consequently, the thus refined geoid serves to compute new estimates of the regional MDT and geostrophic currents.

Iberian Margin 420 kyr geochemical and color variations record

Here we report 420 kyr long records of sediment geochemical and color variations from the southwestern Iberian Margin. We synchronized the Iberian Margin sediment record to Antarctic ice cores and speleothem records on millennial time scales and investigated the phase responses relative to orbital forcing of multiple proxy records available from these cores. Iberian Margin sediments contain strong precession power. Sediment "redness" (a* and 570-560 nm) and the ratio of long-chain alcohols to n-alkanes (C26OH/(C26OH + C29)) are highly coherent and in-phase with precession. Redder layers and more oxidizing conditions (low alcohol ratio) occur near precession minima (summer insolation maxima). We suggest these proxies respond rapidly to low-latitude insolation forcing by wind-driven processes (e.g., dust transport, upwelling, precipitation). Most Iberian Margin sediment parameters lag obliquity maxima by 7-8 ka, indicating a consistent linear response to insolation forcing at obliquity frequencies driven mainly by high-latitude processes. Although the lengths of the time series are short (420 ka) for detecting 100 kyr eccentricity cycles, the phase relationships support those obtained by Shackleton []. Antarctic temperature and the Iberian Margin alcohol ratios (C26OH/(C26OH + C29)) lead eccentricity maxima by 6 kyr, with lower ratios (increased oxygenation) occurring at eccentricity maxima. CO2, CH4, and Iberian SST are nearly in phase with eccentricity, and minimum ice volume (as inferred from Pacific d18Oseawater) lags eccentricity maxima by 10 kyr. The phase relationships derived in this study continue to support a potential role of the Earth's carbon cycle in contributing to the 100 kyr cycle.

Integrated Paleocene calcareous plankton magnetobiochronology of DSDP Hole 43-384 in the Northwest Atlantic Ocean

At Deep Sea Drilling Site 384 (J-Anomaly Ridge, Grand Banks Continental Rise, NW Atlantic Ocean) Paleocene nannofossil chalks and oozes (~70 m thick) are unconformably/disconformably underlain (~168 m; upper Maastrichtian) and overlain (~98.7 m; upper lower Eocene) by sediments of comparable lithologies. The chalks are more indurated in stratigraphically higher levels of the Paleocene reflecting increasing amounts of biosiliceous (radiolarians and diatoms) components. This site serves as an excellent location for an integrated calcareous and siliceous microfossil zonal stratigraphy and stable isotope stratigraphy. We report the results of a magnetostratigraphic study which, when incorporated with published magnetostratigraphic results, reveals an essentially complete magnetostratigraphic record spanning the interval from Magnetochron C31n (late Maastrichtian) to C25n (partim) (late Paleocene, Thanetian). Integrated magnetobiochronology and stable isotope stratigraphy support the interpretation of, and constrain the estimated duration of, a short hiatus (~0.9 my) within the younger part of Chron C29r (including the K/P boundary) and an ~6 my hiatus separating upper Paleocene (Magnetozone C25n) and upper lower Eocene (Magnetozone C22r) sediments. Some 30 planktonic foraminiferal datum levels [including the criteria used to denote the Paleocene planktonic foraminiferal (sub)tropical zonal scheme of Berggren and Miller, Micropaleontology 34 (4) (1988) 362-380 and Berggren et al., SEPM Spec. Publ. 54 (1995) 129-212, Geol. Soc. Am. Bull. 107 (11) (1995) 1272-1287], and nearly two dozen calcareous nannoplankton datum levels have been recognized and calibrated to the magnetochronology. Planktonic foraminiferal Subzones P4a and P4b of (upper Paleocene) Zone P4 are emended/redefined based on the discovery of a longer stratigraphic extension of Acarinina subsphaerica (into at last Magnetozone C25n). Stable isotope stratigraphies from benthic foraminifera and fine fraction (<38 µm) carbonate have been calibrated to the biochronology and magnetostratigraphy. A minimum in benthic foraminifer delta13C was reached near the Danian/Selandian boundary (within Chron C26r, planktonic foraminiferal Zone P3a and calcareous nannoplankton Zone NP4) and is followed by the rise to maximum delta13C values in the late Thanetian (near the base of C25n, in Zone P4c and NP9a, respectively) that can be used for global correlation in the Paleocene.

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.

Consolidation properties and permeability values of sediments from three ODP Leg 205 sites

Vertical permeability and sediment consolidation measurements were taken on seven whole-round drill cores from Sites 1253 (three samples), 1254 (one sample), and 1255 (three samples) drilled during Ocean Drilling Program Leg 205 in the Middle America Trench off of Costa Rica's Pacific Coast. Consolidation behavior including slopes of elastic rebound and virgin compression curves (Cc) was measured by constant rate of strain tests. Permeabilities were determined from flow-through experiments during stepped-load tests and by using coefficient of consolidation (Cv) values continuously while loading. Consolidation curves and the Casagrande method were used to determine maximum preconsolidation stress. Elastic slopes of consolidation curves ranged from 0.097 to 0.158 in pelagic sediments and 0.0075 to 0.018 in hemipelagic sediments. Cc values ranged from 1.225 to 1.427 for pelagic carbonates and 0.504 to 0.826 for hemipelagic clay-rich sediments. In samples consolidated to an axial stress of ~20 MPa, permeabilities determined by flow-through experiments ranged from a low value of 7.66 x 10**-20 m**2 in hemipelagic sediments to a maximum value of 1.03 x 10**-16 m**2 in pelagic sediments. Permeabilities calculated from Cv values in the hemipelagic sediments ranged from 4.81 x 10**-16 to 7.66 x 10**-20 m**2 for porosities 49.9%-26.1%.

Global and Regional Sea level budget components from GRACE and radar altimetry (2002-2014)

This dataset contains the result of a joint least squares inversion of GRACE and altimetry data. The results are evaluated in terms of sea level change for the global ocean as well as dedicated areas. In addition, some auxiliary data is provided to enable reproducibility of the results in Rietbroek et al. 2016, and a google Earth kmz file is provided which visualizes the trends derived from the inversion results.

Grain size analysis of sediment core GeoB1510-2

Detailed information about the sediment properties and microstructure can be provided through the analysis of digital ultrasonic P wave seismograms recorded automatically during full waveform core logging. The physical parameter which predominantly affects the elastic wave propagation in water-saturated sediments is the P wave attenuation coefficient. The related sedimentological parameter is the grain size distribution. A set of high-resolution ultrasonic transmission seismograms (-50-500 kHz), which indicate downcore variations in the grain size by their signal shape and frequency content, are presented. Layers of coarse-grained foraminiferal ooze can be identified by highly attenuated P waves, whereas almost unattenuated waves are recorded in fine-grained areas of nannofossil ooze. Color -encoded pixel graphics of the seismograms and instantaneous frequencies present full waveform images of the lithology and attenuation. A modified spectral difference method is introduced to determine the attenuation coefficient and its power law a = kF. Applied to synthetic seismograms derived using a "constant Q" model, even low attenuation coefficients can be quantified. A downcore analysis gives an attenuation log which ranges from -700 dB/m at 400 kHz and a power of n=1-2 in coarse-grained sands to few decibels per meter and n :s; 0.5 in fine-grained clays. A least squares fit of a second degree polynomial describes the mutual relationship between the mean grain size and the attenuation coefficient. When it is used to predict the mean grain size, an almost perfect coincidence with the values derived from sedimentological measurements is achieved.