Geochemistry and isotopes at DSDP Site 76-533

Concentrations and d34S and d13C values were determined on SO4, HCO3, CO2, and CH4 in interstitial water and gas samples from the uppermost 400 m of sediment on the Blake Outer Ridge. These measurements provide the basis for detailed interpretation of diagenetic processes associated with anaerobic respiration of electrons generated by organic- matter decomposition. The sediments are anaerobic at very shallow depths (<1 m) below the seafloor. Sulfate reduction is confined to the uppermost 15 m of sediment and results in a significant outflux of oxidized carbon from the sediments. At the base of the sulfate reduction zone, upward-diffusing CH4 is being oxidized, apparently in conjunction with SO4 reduction. CH4 generation by CO2 reduction is the most important metabolic process below the 15-m depth. CO2 removal is more rapid than CO2 input over the depth interval from 15 to 100 m, and results in a slight decrease in HCO3 concentration accompanied by a 40 per mil positive shift in d13C. The differences among coexisting CH4, CO2, and HCO3 are consistent with kinetic fractionation between CH4 and dissolved CO2, and equilibrium fractionation between CO2 and HCO3. At depths greater than 100 m, the rate of input of CO2 (d13C = -25 per mil) exceeds by 2 times the rate of removal of CO2 by conversion to CH4 (d13C of -60 to -65 per mil). This results in an increase of dissolved HCO3 concentration while maintaining d13C of HCO3 relatively constant at +10 per mil. Non-steady-state deposition has resulted in significantly higher organic carbon contents and unusually high (70 meq/l) pore-water alkalinities below 150 m. These high alkalinities are believed to be related more to spontaneous decarboxylation reactions than to biological processes. The general decrease in HCO3 concentration with constant d13C over the depth interval of 200 to 400 m probably reflects increased precipitation of authigenic carbonate. Input-output carbon isotope-mass balance calculations, and carbonate system equilibria in conjunction with observed CO2-CH4 ratios in the gas phase, independently suggest that CH4 concentrations on the order of 100 mmol/kg are present in the pore waters of Blake Outer Ridge sediments. This quantity of CH4 is believed to be insufficient to saturate pore waters and stabilize the CH4*6H2O gas hydrate. Results of these calculations are in conflict with the physical recovery of gas hydrate from 238 m, and with the indirect evidence (seismic reflectors, sediment frothing, slightly decreasing salinity and chlorinity with depth, and pressure core barrel observations) of gas-hydrate occurrence in these sediments. Resolution of this apparent conflict would be possible if CH4 generation were restricted to relatively thin (1-10 m) depth intervals, and did not occur uniformly at all depths throughout the sediment column, or if another methanogenic process (e.g., acetate fermentation) were a major contributor of gas.

(Table 1) Snow and sea ice thickness, solar heat input and solar heat used for melting at seven ice mass-balance buoys in the Arctic in 2008

There has been a marked decline in the summer extent of Arctic sea ice over the past few decades. Data from autonomous ice mass-balance buoys can enhance our understanding of this decline. These buoys monitor changes in snow deposition and ablation, ice growth, and ice surface and bottom melt. Results from the summer of 2008 showed considerable large-scale spatial variability in the amount of surface and bottom melt. Small amounts of melting were observed north of Greenland, while melting in the southern Beaufort Sea was quite large. Comparison of net solar heat input to the ice and heat required for surface ablation showed only modest correlation. However, there was a strong correlation between solar heat input to the ocean and bottom melting. As the ice concentration in the Beaufort Sea region decreased, there was an increase in solar heat to the ocean and an increase in bottom melting.

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.

Chert intervals in DSDP and ODP sites (Table 1)

Radiolarian cherts in the Tethyan realm of Jurassic age were recently interpreted as resulting from high biosiliceous productivity along upwelling zones in subequatorial paleolatitudes the locations of which were confirmed by revised paleomagnetic estimates. However, the widespread occurrence of cherts in the Eocene suggests that cherts may not always be reliable proxies of latitude and upwelling zones. In a new survey of the global spatio-temporal distribution of Cenozoic cherts in Deep Sea Drilling Project (DSDP) and Ocean Drilling Program (ODP) sediment cores, we found that cherts occur most frequently in the Paleocene and early Eocene, with a peak in occurrences at ~50 Ma that is coincident with the time of highest bottom water temperatures of the early Eocene climatic optimum (EECO) when the global ocean was presumably characterized by reduced upwelling efficiency and biosiliceous productivity. Cherts occur less commonly during the subsequent Eocene global cooling trend. Primary paleoclimatic factors rather than secondary diagenetic processes seem therefore to control chert formation. This timing of peak Eocene chert occurrence, which is supported by detailed stratigraphic correlations, contradicts currently accepted models that involve an initial loading of large amounts of dissolved silica from enhanced weathering and/or volcanism in a supposedly sluggish ocean of the EECO, followed during the subsequent middle Eocene global cooling by more vigorous oceanic circulation and consequent upwelling that made this silica reservoir available for enhanced biosilicification, with the formation of chert as a result of biosilica transformation during diagenesis. Instead, we suggest that basin-basin fractionation by deep-sea circulation could have raised the concentration of EECO dissolved silica especially in the North Atlantic, where an alternative mode of silica burial involving widespread direct precipitation and/or absorption of silica by clay minerals could have been operative in order to maintain balance between silica input and output during the upwelling-deficient conditions of the EECO. Cherts may therefore not always be proxies of biosiliceous productivity associated with latitudinally focused upwelling zones.

Lagrangian back-tracing results for Kohnen station; output of time dependent numerical experiment with a 3D nested Antarctic ice sheet model

A nested ice flow model was developed for eastern Dronning Maud Land to assist with the dating and interpretation of the EDML deep ice core. The model consists of a high-resolution higher-order ice dynamic flow model that was nested into a comprehensive 3-D thermomechanical model of the whole Antarctic ice sheet. As the drill site is on a flank position the calculations specifically take into account the effects of horizontal advection as deeper ice in the core originated from higher inland. First the regional velocity field and ice sheet geometry is obtained from a forward experiment over the last 8 glacial cycles. The result is subsequently employed in a Lagrangian backtracing algorithm to provide particle paths back to their time and place of deposition. The procedure directly yields the depth-age distribution, surface conditions at particle origin, and a suite of relevant parameters such as initial annual layer thickness. This paper discusses the method and the main results of the experiment, including the ice core chronology, the non-climatic corrections needed to extract the climatic part of the signal, and the thinning function. The focus is on the upper 89% of the ice core (appr. 170 kyears) as the dating below that is increasingly less robust owing to the unknown value of the geothermal heat flux. It is found that the temperature biases resulting from variations of surface elevation are up to half of the magnitude of the climatic changes themselves.

Phosphorus components, CaCO3 and biogenic opal in sediments from ODP Leg 167 sites

The ocean history of reactive phosphorus (P) (i.e., dissolved P available to fuel oceanic primary productivity) is of interest because of the role of P as a biolimiting nutrient, and knowledge of P burial in marine sediments is key to testing hypotheses about temporal changes in P input or output fluxes. Our understanding of the history of the P cycle over the Cenozoic has increased substantially with temporal records of reactive P mass accumulation rates from open-ocean Pacific and Atlantic equatorial sites. However, questions about the relative importance of nutrient burial in ocean-margin sediments relative to burial in open-ocean sediments and about the extent of P remobilization in organic-rich, reducing environments characteristic of margin sediments remain unresolved. Nutrient burial in oceanic boundary current systems has been suggested to have a controlling role in oceanic nutrient budgets in certain time intervals (Vincent and Berger, 1985, doi:10.1029/GM032p0455), with higher sediment accumulation rates balancing the limited spatial extent of these sediments. Some investigators suggest that remobilization of P from reducing sediments in margin settings is a significant positive feedback to primary productivity (e.g., Van Cappellan and Ingall, 1994, doi:10.1029/94PA01455), whereas other results indicate that both P uptake and P release may occur in these settings depending on the balance of organic carbon and iron supply to the sediments and on the oxygenation of bottom waters (McManus et al., 1997, doi:10.1016/S0016-7037(97)00138-5). It is important to quantitatively understand the geochemistry of reactive P in margin sediments, where productivity and delivery of organic-rich material to the sediments in relatively shallow-water settings is often sufficient to promote anoxia in interstitial waters. To address these questions, we determined the P concentrations and geochemistry in sediment samples from eight sites drilled during Ocean Drilling Program (ODP) Leg 167, California margin (Sites 1010-1012, 1014, 1016-1017, and 1021-1022). These results are the first records of reactive P concentrations on long time scales-required for the calculation of P accumulation rates-for sediments from a highly productive eastern boundary current setting. In addition, we determined calcium carbonate contents and biogenic silica concentrations to define the environments of sedimentary production, burial, and diagenesis.

Ocean acidification decreases the light-use efficiency in an Antarctic diatom under dynamic but not constant light, link to supplementary data

There is increasing evidence that different light intensities strongly modulate the effects of ocean acidification (OA) on marine phytoplankton. The aim of the present study was to investigate interactive effects of OA and dynamic light, mimicking natural mixing regimes. The Antarctic diatom Chaetoceros debilis was grown under two pCO2 (390 and 1000 latm) and light conditions (constant and dynamic), the latter yielding the same integrated irradiance over the day. To characterize interactive effects between treatments, growth, elemental composition, primary production and photophysiology were investigated. Dynamic light reduced growth and strongly altered the effects of OA on primary production, being unaffected by elevated pCO2 under constant light, yet significantly reduced under dynamic light. Interactive effects between OA and light were also observed for Chl production and particulate organic carbon (POC) quotas. Response patterns can be explained by changes in the cellular energetic balance. While the energy transfer efficiency from photochemistry to biomass production (Phi_e,C) was not affected by OA under constant light, it was drastically reduced under dynamic light. Contrasting responses under different light conditions need to be considered when making predictions regarding a more stratified and acidified future ocean.

Mass balance and acceleration of the Greenland ice sheet

Within the last decade, the Greenland ice sheet (GrIS) and its surroundings have experienced record high surface temperatures (Mote, 2007, doi:10.1029/2007GL031976; Box et al., 2010), ice sheet melt extent (Fettweis et al., 2011, doi:10.5194/tc-5-359-2011) and record-low summer sea-ice extent (Nghiem et al., 2007, doi:10.1029/2007GL031138). Using three independent data sets, we derive, for the first time, consistent ice-mass trends and temporal variations within seven major drainage basins from gravity fields from the Gravity Recovery and Climate Experiment (GRACE; Tapley et al., 2004, doi:10.1029/2004GL019920), surface-ice velocities from Inteferometric Synthetic Aperture Radar (InSAR; Rignot and Kanagaratnam, 2006, doi:10.1126/science.1121381) together with output of the regional atmospheric climate modelling (RACMO2/ GR; Ettema et al., 2009, doi:10.1029/2009GL038110), and surface-elevation changes from the Ice, cloud and land elevation satellite (ICESat; Sorensen et al., 2011, doi:10.5194/tc-5-173-2011). We show that changing ice discharge (D), surface melting and subsequent run-off (M/R) and precipitation (P) all contribute, in a complex and regionally variable interplay, to the increasingly negative mass balance of the GrIS observed within the last decade. Interannual variability in P along the northwest and west coasts of the GrIS largely explains the apparent regional mass loss increase during 2002-2010, and obscures increasing M/R and D since the 1990s. In winter 2002/2003 and 2008/2009, accumulation anomalies in the east and southeast temporarily outweighed the losses by M/R and D that prevailed during 2003-2008, and after summer 2010. Overall, for all basins of the GrIS, the decadal variability of anomalies in P, M/R and D between 1958 and 2010 (w.r.t. 1961-1990) was significantly exceeded by the regional trends observed during the GRACE period (2002-2011).

Grain size analysis of ODP Site 162-984

Past changes in the freshwater balance of the surface North Atlantic Ocean are thought to have influenced the rate of deep-water formation, and consequently climate (Broecker and Denton, 1989, doi:10.1016/0016-7037(89)90123-3; Manabe and Stouffer, 1996; doi:10.1038/378165a0). Although water-mass proxies are generally consistent with an impact of freshwater input on meridional overturning circulation (Boyle and Keigwin, 1987, doi:10.1038/330035a0), there has been little dynamic evidence to support this linkage. Here we present a 25,000 year record of variations in sediment grain size from south of Iceland, which indicates vigorous bottom-water currents during both the last glacial maximum and the Holocene period. Together with reconstructions of North Atlantic water-mass distribution, vigorous bottom currents suggest a shorter residence time of northern-source waters during the last glacial maximum, relative to the Holocene period. The most significant reductions in flow strength occur during periods that have been associated with freshening of the surface North Atlantic. The short-term deglacial oscillations in bottom current strength are closely coupled to changes in Greenland air temperature, with a minimum during the Younger Dryas cold reversal and a maximum at the time of rapid warming at the onset of the Holocene. Our results support a strong connection between ocean circulation and rapid climate change.

Chemical composition of ODP Leg 138 sediments

A major oceanographic event preserved in the Cocos plate sedimentary column survived subduction and is recorded in the changing composition of Nicaraguan magmas. A uranium increase in these magmas since the latest Miocene (after 7 Ma) resulted from the 'carbonate crash' at 10 Ma and the ensuing high organic carbon burial in the sediments. The response of the arc to this paleoceanographic event requires near steady-state sediment recycling at this margin since 20 Ma. This relative stability in sediment subduction invites one of the first attempts to balance sedimentary input and arc output across a subduction zone. Calculations based on Th indicate that as much as 75% of the sedimentary column was subducted beneath the arc. The Nicaraguan margin is one of the few places to observe such strong links between the oceans and the solid earth.