Annotated record of the detailed examination of Mn nodules recovered from the VA04 expedition (R/V Valdivia) in the Pacific Ocean

Twenty-four manganese nodules from the surface of the sea floor and fifteen buried nodules were studied. With three exceptions, the nodules were collected from the area covered by Valdivia Cruise VA 04 some 1200 nautical miles southeast of Hawaii. Age determinations were made using the ionium method. In order to get a true reproduction of the activity distribution in the nodules, they were cut in half and placed for one month on nuclear emulsion plates to determine the alpha-activity of the ionium and its daughter products. Special methods of counting the alpha-tracks resolution to depth intervals of 0.125 mm. For the first time it was possible to resolve zones of rapid growth (impulse growth) with growth rates, s > 50 mm/106 yr and interruptions in growth. With few exceptions the average rate of growth of all nodules was surprisingly uniform at 4-9 mm/10 yr. No growth could be recognized radioactively in the buried nodules. One exceptional nodule has had recent impulse growth and, in the material formed, the ionium is not yet in equilibrium with its daughter products. Individual layers in one nodule from the Indian Ocean could be dated and an average time interval of t = 2600±400 yr was necessary to form one layer. The alternation between iron and manganese-rich parts of the nodules was made visible by colour differences resulting from special treatment of cut surfaces with HCl vapour. The zones of slow growth of one nodule are relatively enriched in iron. Earlier attempts to find paleomagnetic reversals in manganese nodules have been continued. Despite considerable improvement in areal resolution, reversals were not detected in the nodules studied. Comparisons of the surface structure, microstructure in section and the radiometric dating show that there are erosion surfaces and growth surfaces on the outer surfaces of the manganese nodules. The formation of cracks in the nodules was studied in particular. The model of age-dependent nodule shrinkage and cracking surprisingly indicates that the nodules break after exceeding a certain age and/or size. Consequently, the breaking apart of manganese nodules is a continuous process not of catastrophic or discontinuous origin. The microstructure of the nodules exhibits differences in the mechanism of accretion and accretion rate of material, shortly referred to as accretion form. Thus non-directional growth inside the nodules as well as a directional growth may be observed. Those nodules with large accretion forms have grown faster than smaller ones. Consequently, parallel layers indicate slow growth. The upper surfaces of the nodules, protruding into the bottom water appear to be more prone to growth disturbances than the lower surfaces, immersed in the sediment. Features of some nodules show, that as they develop, they neither turned nor rolled. Yet unknown is the mechanism that keeps the nodules at the surface during continuous sedimentation. All in all, the nodules remain the objects of their own distinctive problems. The hope of using them as a kind of history book still seems to be very remote.

(Table 1) Water column biogenic silica concentrations and Si isotopic composition during Marion Dufresne cruise MD166 off South Africa

Southern Ocean biogeochemical processes have an impact on global marine primary production and global elemental cycling, e.g. by likely controlling glacial-interglacial pCO2 variation. In this context, the natural silicon isotopic composition (d30Si) of sedimentary biogenic silica has been used to reconstruct past Si-consumption:supply ratios in the surface waters. We present a new dataset in the Southern Ocean from a IPY-GEOTRACES transect (Bonus-GoodHope) which includes for the first time summer d30Si signatures of suspended biogenic silica (i) for the whole water column at three stations and (ii) in the mixed layer at seven stations from the subtropical zone up to the Weddell Gyre. In general, the isotopic composition of biogenic opal exported to depth was comparable to the opal leaving the mixed layer and did not seem to be affected by any diagenetic processes during settling, even if an effect of biogenic silica dissolution cannot be ruled out in the northern part of the Weddell Gyre. We develop a mechanistic understanding of the processes involved in the modern Si-isotopic balance, by implementing a mixed layer model. We observe that the accumulated biogenic silica (sensu Rayleigh distillation) should satisfactorily describe the d30Si composition of biogenic silica exported out of the mixed layer, within the limit of the current analytical precision on the d30Si. The failures of previous models (Rayleigh and steady state) become apparent especially at the end of the productive period in the mixed layer, when biogenic silica production and export are low. This results from (1) a higher biogenic silica dissolution:production ratio imposing a lower net fractionation factor and (2) a higher Si-supply:Si-uptake ratio supplying light Si-isotopes into the mixed layer. The latter effect is especially expressed when the summer mixed layer becomes strongly Si-depleted, together with a large vertical silicic acid gradient, e.g. in the Polar Front Zone and at the Polar Front.

Calcium isotopes in planktonic foraminifera - core tops and sediment trap samples

For paleoceanographic studies, it is important to understand the processes that influence the calcium (Ca) isotopic composition of foraminiferal calcite tests preserved in the sediment record. Seven species of planktonic foraminifera from coretop sediments collectively exhibited a Ca temperature dependent fractionation of 0.013 per mil per °C. This is in agreement with previously published estimates for most species of planktonic foraminifera as well as biogenic and inorganic calcite and aragonite. Four species of planktonic foraminifera collected from a sediment trap showed a considerable amount of scatter and no consistent temperature dependent fractionation. Analyzed size fractions of coretop samples show no significant relationship with d44/40Ca. However, preliminary results suggest that the symbiotic and spinose foraminifera G. sacculifer might exhibit a relationship between test size and d44/40Ca. A one-box model in which Ca isotopes are allowed to fractionate by Rayleigh distillation from a biomineralization reservoir (internal pool) was used to constrain the isotopic composition of the original biomineralization Ca reservoir, assuming around 85% of the Ca reservoir is precipitated and the fractionation factor during precipitation is 0.9985 + 0.00002(T ºC). To explain the foraminiferal Ca isotope data, this model indicates that the Ca isotopic composition of the biomineralization reservoir is offset from seawater (approximately -0.8per mil).

Sr/Ca and Cd/Ca ratios of benthic foraminifera in surface deep-sea sediments

Measurements of Sr/Ca of benthic foraminifera show a linear decrease with water depth which is superimposed upon significant variability identified by analyses of individual foraminifera. New data for Cd/Ca support previous work in defining a contrast between waters shallower and deeper than ~2500 m. Measured element partition coefficients in foraminiferal calcium carbonate relative to sea water (D) have been described by means of a one-box model in which elements are extracted by Rayleigh distillation from a biomineralization reservoir that serves for calcification with a constant fractionation factor (alpha), such that D = (1 - f**alpha)/(l - f), where f is the proportion of Ca remaining after precipitation. A modification to the model recognises differences in element speciation. The model is consistent with differences between D[Sr], D[Ba], and D[Cd] in benthic but not planktonic foraminifera. Depth variations in D for Sr and Ba are consistent with the model, as are differences in depth variation of D[Cd] in calcitic and aragonitic benthic foraminifera. The shallower depth variations may reflect increasing calcification rates with increasing water depth to an optimum of about 2500 m. Observations of unusually lower DCd for some deep waters, not accompanied by similar [Sr], or D[Ba] may be because of dissolution or a calcification response to a lower carbonate saturation state.

Mean dynamic topography and oceanographic parameters estimated from an inverse model and satellite geodesy, with link to model result in one single NetCDF file (392 MB), including the inverse data error covariance

Geostrophic surface velocities can be derived from the gradients of the mean dynamic topography-the difference between the mean sea surface and the geoid. Therefore, independently observed mean dynamic topography data are valuable input parameters and constraints for ocean circulation models. For a successful fit to observational dynamic topography data, not only the mean dynamic topography on the particular ocean model grid is required, but also information about its inverse covariance matrix. The calculation of the mean dynamic topography from satellite-based gravity field models and altimetric sea surface height measurements, however, is not straightforward. For this purpose, we previously developed an integrated approach to combining these two different observation groups in a consistent way without using the common filter approaches (Becker et al. in J Geodyn 59(60):99-110, 2012, doi:10.1016/j.jog.2011.07.0069; Becker in Konsistente Kombination von Schwerefeld, Altimetrie und hydrographischen Daten zur Modellierung der dynamischen Ozeantopographie, 2012, http://nbn-resolving.de/nbn:de:hbz:5n-29199). Within this combination method, the full spectral range of the observations is considered. Further, it allows the direct determination of the normal equations (i.e., the inverse of the error covariance matrix) of the mean dynamic topography on arbitrary grids, which is one of the requirements for ocean data assimilation. In this paper, we report progress through selection and improved processing of altimetric data sets. We focus on the preprocessing steps of along-track altimetry data from Jason-1 and Envisat to obtain a mean sea surface profile. During this procedure, a rigorous variance propagation is accomplished, so that, for the first time, the full covariance matrix of the mean sea surface is available. The combination of the mean profile and a combined GRACE/GOCE gravity field model yields a mean dynamic topography model for the North Atlantic Ocean that is characterized by a defined set of assumptions. We show that including the geodetically derived mean dynamic topography with the full error structure in a 3D stationary inverse ocean model improves modeled oceanographic features over previous estimates.

ETOPO1 Global Relief Model converted to PanMap layer format

ETOPO1 is a 1 arc-minute global relief model of Earth's surface that integrates land topography and ocean bathymetry. It was built from numerous global and regional data sets. Data were converted to the PanMap layer format in 14 contour lines from 500 to 7000 meter in steps of 500 m. The link provides a zip-archive (1.1 MB) with *.lay files. The PanMap Mini-GIS software is published at doi:10.1594/PANGAEA.104840.

Static GPS measurements and elevation model in the vicinity of the EDML deep-drilling site

Interpretation of ice-core records requires accurate knowledge of the past and present surface topography and stress-strain fields. The European Project for Ice Coring in Antarctica (EPICA) drilling site (0.0684° E and 75.0025° S, 2891.7 m) in Dronning Maud Land, Antarctica, is located in the immediate vicinity of a transient and splitting ice divide. A digital elevation model is determined from the combination of kinematic GPS measurements with the GLAS12 data sets from the ICESat satellite. Based on a network of stakes, surveyed with static GPS, the velocity field around the EDML drilling site is calculated. The annual mean velocity magnitude of 12 survey points amounts to 0.74 m/a. Flow directions mainly vary according to their distance from the ice divide. Surface strain rates are determined from a pentagon-shaped stake network with one center point, close to the drilling site. The strain field is characterised by along flow compression, lateral dilatation, and vertical layer thinning.

Simulated Bromoform emission and concentration using the ocean biogeochemistry model MPIOM/HAMOCC forced by 6-hourly NCEP data

Bromoform (CHBr3) is one important precursor of atmospheric reactive bromine species that are involved in ozone depletion in the troposphere and stratosphere. In the open ocean bromoform production is linked to phytoplankton that contains the enzyme bromoperoxidase. Coastal sources of bromoform are higher than open ocean sources. However, open ocean emissions are important because the transfer of tracers into higher altitude in the air, i.e. into the ozone layer, strongly depends on the location of emissions. For example, emissions in the tropics are more rapidly transported into the upper atmosphere than emissions from higher latitudes. Global spatio-temporal features of bromoform emissions are poorly constrained. Here, a global three-dimensional ocean biogeochemistry model (MPIOM-HAMOCC) is used to simulate bromoform cycling in the ocean and emissions into the atmosphere using recently published data of global atmospheric concentrations (Ziska et al., 2013) as upper boundary conditions. Our simulated surface concentrations of CHBr3 match the observations well. Simulated global annual emissions based on monthly mean model output are lower than previous estimates, including the estimate by Ziska et al. (2013), because the gas exchange reverses when less bromoform is produced in non-blooming seasons. This is the case for higher latitudes, i.e. the polar regions and northern North Atlantic. Further model experiments show that future model studies may need to distinguish different bromoform-producing phytoplankton species and reveal that the transport of CHBr3 from the coast considerably alters open ocean bromoform concentrations, in particular in the northern sub-polar and polar regions.

Using ground-penetrating radar, topography and classification of vegetation to model the sediment and active layer thickness in a periglacial lake catchment, Western Greenland, link to shapefiles

The geometries of a catchment constitute the basis for distributed physically based numerical modeling of different geoscientific disciplines. In this paper results from ground-penetrating radar (GPR) measurements, in terms of a 3D model of total sediment thickness and active layer thickness in a periglacial catchment in western Greenland, is presented. Using the topography, thickness and distribution of sediments is calculated. Vegetation classification and GPR measurements are used to scale active layer thickness from local measurements to catchment scale models. Annual maximum active layer thickness varies from 0.3 m in wetlands to 2.0 m in barren areas and areas of exposed bedrock. Maximum sediment thickness is estimated to be 12.3 m in the major valleys of the catchment. A method to correlate surface vegetation with active layer thickness is also presented. By using relatively simple methods, such as probing and vegetation classification, it is possible to upscale local point measurements to catchment scale models, in areas where the upper subsurface is relatively homogenous. The resulting spatial model of active layer thickness can be used in combination with the sediment model as a geometrical input to further studies of subsurface mass-transport and hydrological flow paths in the periglacial catchment through numerical modelling.