Geochemistry of terrigenous material in deep-sea sediments of the equatorial Pacific

Biological productivity in the modern equatorial Pacific Ocean, a region with high nutrients and low chlorophyll, is currently limited by the micronutrient Fe. In order to test whether Fe was limiting in the past and to identify potential pathways of Fe delivery that could drive Fe fertilization (i.e., dust delivery from eolian inputs vs. Fe supplied by the Equatorial Undercurrent), we chemically isolated the terrigenous material from sediment along a cross-equatorial transect in the central equatorial Pacific at 140°W and at Ocean Drilling Program Site 850 in the eastern equatorial Pacific. We quantified the contribution from each potential Fe-bearing terrigenous source using a suite of chemical- and isotopic discrimination strategies as well as multivariate statistical techniques. We find that the distribution of the terrigenous sources (i.e., Asian loess, South American ash, Papua New Guinea, and ocean island basalt) varies through time, latitude, and climate. Regardless of which method is used to determine accumulation rate, there also is no relationship between flux of any particular Fe source and climate. Moreover, there is no connection between a particular Fe source or pathway (eolian vs. Undercurrent) to total productivity during the Last Glacial Maximum, Pleistocene glacial episodes, and the Miocene "Biogenic Bloom". This would suggest an alternative process, such as an interoceanic reorganization of nutrient inventories, may be responsible for past changes in total export in the open ocean, rather than simply Fe supply from dust and/or Equatorial Undercurrent processes. Additionally, perhaps a change in Fe source or flux is related to a change in a particular component of the total productivity (e.g., the production of organic matter, calcium carbonate, or biogenic opal).

Biostratigraphy and Sr isotopic composition of sediments from DSDP Legs 12, 14, 15, 17, 20-22, and 30

We present the data used to construct the Cenozoic and Cretaceous portion of the Phanerozoic curve of seawater 87Sr/86Sr that had been given in summary form by W.H. Burke co-workers. All Cenozoic samples (128) and 22 Cretaceous samples are foram-nannofossil oozes and limestones from DSDP cores distributed among 13 sites in the Atlantic, Pacific and Indian Oceans, and the Caribbean Sea. Non-DSDP Cretaceous samples (126) include limestone, anhydrite and phosphate samples from North America, Europe and Asia. Determination of the 87Sr/86Sr value of seawater at particular times in the past is based on comparison of ratios derived from coeval marine samples from widely separated geographic areas. These samples are characterized by a wide variety of diagenetic and burial histories. The large size and cosmopolitan nature of the data set decreases the likelihood that, among coeval data, systematic error has been introduced by a similar pattern of diagenetic alteration of the ratios. There is good clustering of data points throughout the Cenozoic and Cretaceous curve. The consistency of data is illustrated by Cenozoic and Cretaceous data plots that include a separate symbol for each DSDP site and non-DSDP sample location. More than 98% of the data points are enclosed by upper and lower lines that define a narrow band. For any given time, the correct seawater ratio probably lies within this band. A line drawn within the band represents our estimate of the actual seawater ratio as a function of time. The general configuration of the Cenozoic and Cretaceous curve appears to be strongly influenced by the history of plate interactions and sea-floor spreading. Specific rises and falls in the 87Sr/86Sr of seawater, however, may be caused by a variety of factors such as variation in lithologic composition of the crust exposed to weathering, configuration and topographic relief of continents, volcanic activity, rate of sea-floor spreading, extent of continental inundation by epeiric seas, and variations in both climate and paleooceanographic conditions. Many or all of these factors are probably related to global tectonic processes, yet their combined effect on the temporal variation of seawater 87Sr/86Sr can complicate a direct platetectonic interpretation for portions of the seawater curve.

Initial biomass and biomass change of tundra and forest plants in three herbivore treatments in Norway and Sweden

Recent Pan-Arctic shrub expansion has been interpreted as a response to a warmer climate. However, herbivores can also influence the abundance of shrubs in arctic ecosystems. We addressed these alternative explanations by following the changes in plant community composition during the last 10 years in permanent plots inside and outside exclosures with different mesh sizes that exclude either only reindeer or all mammalian herbivores including voles and lemmings. The exclosures were replicated at three forest and tundra sites at four different locations along a climatic gradient (oceanic to continental) in northern Fennoscandia. Since the last 10 years have been exceptionally warm, we could study how warming has influenced the vegetation in different grazing treatments. Our results show that the abundance of the dominant shrub, Betula nana, has increased during the last decade, but that the increase was more pronounced when herbivores were excluded. Reindeer have the largest effect on shrubs in tundra, while voles and lemmings have a larger effect in the forest. The positive relationship between annual mean temperature and shrub growth in the absence of herbivores and the lack of relationships in grazed controls is another indication that shrub abundance is controlled by an interaction between herbivores and climate. In addition to their effects on taller shrubs (> 0.3 m), reindeer reduced the abundance of lichens, whereas microtine rodents reduced the abundance of dwarf shrubs (< 0.3 m) and mosses. In contrast to short-term responses, competitive interactions between dwarf shrubs and lichens were evident in the long term. These results show that herbivores have to be considered in order to understand how a changing climate will influence tundra ecosystems.

Planimetric ice-flow convergence and flow-orthonormal strain rate across the Antarctic Ice Sheet, with links to model result files

Fast-flowing ice streams discharge most of the ice from the interior of the Antarctic Ice Sheet coastward. Understanding how their tributary organisation is governed and evolves is essential for developing reliable models of the ice sheet's response to climate change. Despite much research on ice-stream mechanics, this problem is unsolved, because the complexity of flow within and across the tributary networks has hardly been interrogated. Here I present the first map of planimetric flow convergence across the ice sheet, calculated from satellite measurements of ice surface velocity, and use it to explore this complexity. The convergence map of Antarctica elucidates how ice-stream tributaries draw ice from the interior. It also reveals curvilinear zones of convergence along lateral shear margins of streaming, and abundant convergence ripples associated with nonlinear ice rheology and changes in bed topography and friction. Flow convergence on ice-stream tributaries and their feeding zones is markedly uneven, and interspersed with divergence at distances of the order of kilometres. For individual drainage basins as well as the ice sheet as a whole, the range of convergence and divergence decreases systematically with flow speed, implying that fast flow cannot converge or diverge as much as slow flow. I therefore deduce that flow in ice-stream networks is subject to mechanical regulation that limits flow-orthonormal strain rates. These properties and the gridded data of convergence and flow-orthonormal strain rate in this archive provide targets for ice- sheet simulations and motivate more research into the origin and dynamics of tributarization.

Chemical and Sr isotopic compositions of pore fluids and marine carbonates from ODP Hole 130-807A

The 87Sr/86Sr ratios and Sr concentrations in sediment and pore fluids are used to evaluate the rates of calcite recrystallization at ODP Site 807A on the Ontong Java Plateau, an 800-meter thick section of carbonate ooze and chalk. A numerical model is used to evaluate the pore fluid chemistry and Sr isotopes in an accumulating section. The deduced calcite recrystallization rate is 2% per million years (%/Myr) near the top of the section and decreases systematically in older parts of the section such that the rate is close to 0.1/age (in years). The deduced recrystallization rates have important implications for the interpretation of Ca and Mg concentration profiles in the pore fluids. The effect of calcite recrystallization on pore fluid chemistry is described by the reaction length, L, which varies by element, and depends on the concentration in pore fluid and solid. When L is small compared to the thickness of the sedimentary section, the pore fluid concentration is controlled by equilibrium or steady-state exchange with the solid phase, except within a distance L of the sediment-water interface. When L is large relative to the thickness of sediment, the pore fluid concentration is mostly controlled by the boundary conditions and diffusion. The values of L for Ca, Sr, and Mg are of order 15, 150, and 1500 meters, respectively. L_Sr is derived from isotopic data and modeling, and allows us to infer the values of L_Ca and L_Mg. The small value for L_Ca indicates that pore fluid Ca concentrations, which gradually increase down section, must be equilibrium values that are maintained by solution-precipitation exchange with calcite and do not reflect Ca sources within or below the sediment column. The pore fluid Ca measurements and measured alkalinity allow us to calculate the in situ pH in the pore fluids, which decreases from 7.6 near the sediment-water interface to 7.1+/-0.1 at 400-800 mbsf. While the calculated pH values are in agreement with some of the values measured during ODP Leg 130, most of the measurements are artifacts. The large value for L_Mg indicates that the pore fluid Mg concentrations at 807A are not controlled by calcite-fluid equilibrium but instead are determined by the changing Mg concentration of seawater during deposition, modified by aqueous diffusion in the pore fluids. We use the pore fluid Mg concentration profile at Site 807A to retrieve a global record for seawater Mg over the past 35 Myr, which shows that seawater Mg has increased rapidly over the past 10 Myr, rather than gradually over the past 60 Myr. This observation suggests that the Cenozoic rise in seawater Mg is controlled by continental weathering inputs rather than by exchange with oceanic crust. The relationship determined between reaction rate and age in silicates and carbonates is strikingly similar, which suggests that reaction affinity is not the primary determinant of silicate dissolution rates in nature.

Sulfur isotope composition of pore fluids at ODP Sites 128-798 and 128-799

Micro-crystalline barites recovered by deep-sea drilling from Site 684 on the Peru margin and Site 799 in the Japan Sea are highly enriched in the heavy sulfur isotope relative to seawater ( d34S up to +84?). This isotopic composition is consistent with remobilization of biogenic barite triggered by sulfate reduction, and subsequent reprecipitation as a diagenetic barite front. The high levels of barium sulfate in these deposits (10-50%) cannot be explained by a diffusive transport model in sediments experiencing a constant rate of sedimentation. When sedimentation rates change radically, the barite front will remain at a given depth interval leading to large accumulations of barium sulfate. Such conditions may have generated the barite deposits at Site 799. At Site 684, on the other hand, there is evidence that the barite deposits are a result of the tectonically-driven advection of sulfate-bearing fluids through the sediment column.

Pleistocene stable isotope record of planktonic foraminifera of ODP Site 138-847

We present Pleistocene oxygen and carbon isotope records from two planktonic foraminifer species (Globigerinoides sacculifer and Neogloboquadrina dutertrei) from Ocean Drilling Program Site 847 (0°16'N, 95°19'W; 3334 m water depth). An average sample resolution of 4500 yr was obtained by sampling at an interval of 15 cm through a continuous 35-m section from 0 to 1.15 Ma. Our d18O-based chronology is similar to that derived independently by astronomically tuning the gamma-ray attenuation porosity evaluator (GRAPE) record (Shackleton et al., 1995), though offsets as large as ± 30 k.y. occur on occasion. The surface waters at eastern equatorial Pacific Site 847, 380 km west of the Galapagos, are characterized by strong and constant upwelling, elevated nutrient concentrations, and high productivity. The isotopic composition of G. sacculifer (300-355 µm) reflects conditions in the thin-surface mixed layer, and the composition of N. dutertrei (355-425 µm) monitors the subsurface waters of the permanent shallow (10-40 m) thermocline. The Pleistocene d18O difference (N. dutertrei minus G. sacculifer, Dd18Od-s) averages 0.9 per mil and ranges from 0 per mil to 1.7 per mil. Neglecting species effects and shell size, the average Pleistocene d13C difference (G. sacculifer minus N. dutertrei, Dd13Cs-d) is 0.0 per mil and ranges from -0.5 per mil to 0.5 per mil. The Dd18Od-s and Dd13Cs-d records are used to infer vertical contrasts in upper ocean water temperature and nutrient concentration, though d13C may also be influenced by other factors, such as CO2 gas exchange. Variations in the isotopic differences are often synchronous with glacial/interglacial climate change. Glacial periods are characterized by smaller vertical contrasts in both temperature and nutrient concentration, and by notably greater accumulation rates of N. dutertrei and CaCO3. We attribute these responses to greater upwelling at the equatorial divergence. Superimposed on the glacial/interglacial Dd18Od-s pattern is a long-term trend possibly associated with the advection of Peru Current waters. The temporal fluctuations in the isotopic contrasts are strikingly similar to those observed at Site 851 (Ravelo and Shackleton, this volume), suggesting that the inferred changes in thermal and chemical profiles occurred over a broad region in the equatorial Pacific.

Chemical and isotopic compositions of deep-sea carbonate sediments and interstitial waters from DSDP Sites 30-288 and 30-289

Interstitial waters and sediments from DSDP sites 288 and 289 contain information on the chemistry and diagenesis of carbonate in deep-sea sediments and on the role of volcanic matter alteration processes. Sr/Ca ratios are species dependent in unaltered foraminifera from site 289 and atom ratios (0.0012-0.0016) exceed those predicted by distribution coefficent data (~0.0004). During diagenesis Sr/Ca ratios of carbonates decrease and reach the theoretical distribution at a depth which is identical to the depth of Sr isotopic equilibration, where 87Sr/86Sr ratios of interstitial waters and carbonates converge. Mg/Ca ratios in the carbonates do not increase with depth as found in some other DSDP sites, possibly because of diagenetic re-equilibration with interstitial waters showing decreasing Mg(2+)/Ca(2+) ratios with depth due to Ca input and Mg removal by alteration of volcanic matter. Interstitial 18O/16O ratios increase with depth at site 289 to d18O = 0.67? (SMOW), reflecting carbonate recrystallization at elevated temperatures (>/= 20°C), the first recorded evidence of this effect in interstitial waters. Interstitial Sr2+ concentrations reach high levels, up to 1 mM, chiefly because of carbonate recrystallization. However, 87Sr/86Sr ratios decrease from 0.7092 to less than 0.7078, lower than for contemporaneous sea water, showing that there is a volcanic input of strontium at depth. This volcanic component is recorded in the Sr isotopic composition of recrystallized calcites. Isotopic compositions of the unrecrystallized calcites suggests that the rate of increase of the 87Sr/86Sr ratio of sea water with time has been faster since 3 my ago than in the preceding 13 my.

Isotopic composition of rocks from the Lau/Tonga back-arc basin, data of DSDP Leg 21 and ODP Legs 91 an 135

New trace element, Sr-, Nd-, Pb- and Hf isotope data provide insights into the evolution of the Tonga-Lau Basin subduction system. The involvement of two separate mantle domains, namely Pacific MORB mantle in the pre-rift and early stages of back-arc basin formation, and Indian MORB mantle in the later stages, is confirmed by these results. Contrary to models proposed in recent studies on the basis of Pb isotope and other compositional data, this change in mantle wedge character best explains the shift in the isotopic composition, particularly 143Nd/144Nd ratios, of modern Tofua Arc magmas relative to all other arc products from this region. Nevertheless, significant changes in the slab-derived flux during the evolution of the arc system are also required to explain second order variations in magma chemistry. In this region, the slab-derived flux is dominated by fluid; however, these fluids carry Pb with sediment-influenced isotopic signatures, indicating that their source is not restricted to the subducting altered mafic oceanic crust. This has been the case from the earliest magmatic activity in the arc (Eocene) until the present time, with the exception of two periods of magmatic activity recorded in samples from the Lau Islands. Both the Lau Volcanic Group, and Korobasaga Volcanic Group lavas preserve trace element and isotope evidence for a contribution from subducted sediment that was not transported as a fluid, but possibly in the form of a melt. This component shares similarities with that influencing the chemistry of the northern Tofua Arc magmas, suggesting some caution may be required in the adoption of constraints for the latter dependent upon the involvement of sediments from the Louisville Ridge. A key outcome of this study is to demonstrate that the models proposed to explain subduction zone magmatism cannot afford to ignore the small but important contributions made by the mantle wedge to the incompatible trace element inventory of arc magmas.

(Table 1) Lake elevation change derived from 4 or more years of measurements

In this study, ICESat altimetry data are used to provide precise lake elevations of the Tibetan Plateau (TP) during the period of 2003-2009. Among the 261 lakes examined ICESat data are available on 111 lakes: 74 lakes with ICESat footprints for 4-7 years and 37 lakes with footprints for 1 -3 years. This is the first time that precise lake elevation data are provided for the 111 lakes. Those ICESat elevation data can be used as baselines for future changes in lake levels as well as for changes during the 2003-2009 period. It is found that in the 74 lakes (56 salt lakes) examined, 62 (i.e. 84%) of all lakes and 50 (i.e. 89%) of the salt lakes show tendency of lake level increase. The mean lake water level increase rate is 0.23 m/year for the 56 salt lakes and 0.27 m/year for the 50 salt lakes of water level increase. The largest lake level increase rate (0.80 m/year) found in this study is the lake Cedo Caka. The 74 lakes are grouped into four subareas based on geographical locations and change tendencies in lake levels. Three of the four subareas show increased lake levels. The mean lake level change rates for subareas I, II, III, IV, and the entire TP are 0.12, 0.26, 0.19, -0.11, and 0.2 m/year, respectively. These recent increases in lake level, particularly for a high percentage of salt lakes, supports accelerated glacier melting due to global warming as the most likely cause.

Siliceous microfossil assemblage and accumulation rate variations across the Middle Eocene Climatic Optimum at ODP Site 171-1051

The Middle Eocene Climatic Optimum (MECO; ~ 40 million years ago [Ma]) is one of the most prominent transient global warming events in the Paleogene. Although the event is well documented in geochemical and isotopic proxy records at many locations, the marine biotic response to the MECO remains poorly constrained. We present new high-resolution, quantitative records of siliceous microplankton assemblages from the MECO interval of Ocean Drilling Program (ODP) Site 1051 in the subtropical western North Atlantic Ocean, which are interpreted in the context of published foraminiferal and bulk carbonate stable isotope (d18O and d13C) records. High diatom, radiolarian and silicoflagellate accumulation rates between 40.5 and 40.0 Ma are interpreted to reflect an ~ 500 thousand year (kyr) interval of increased nutrient supply and resultant surface-water eutrophication that was associated with elevated sea-surface temperatures during the prolonged onset of the MECO. Relatively low pelagic siliceous phytoplankton sedimentation accompanied the peak MECO warming interval and the termination of the MECO during an ~ 70 kyr interval centered at ~ 40.0 Ma. Following the termination of the MECO, an ~ 200-kyr episode of increased siliceous plankton abundance indicates enhanced nutrient levels between ~ 39.9 and 39.7 Ma. Throughout the Site 1051 record, abundance and accumulation rate fluctuations in neritic diatom taxa are similar to the trends observed in pelagic taxa, implying either similar controls on diatom production in the neritic and pelagic zones of the western North Atlantic or fluctuations in sea level and/or shelf accommodation on the North American continental margin to the west of Site 1051. These results, combined with published records based on multiple proxies, indicate a geographically diverse pattern of surface ocean primary production changes across the MECO. Notably, however, increased biosiliceous accumulation is recorded at both ODP Sites 1051 and 748 (Southern Ocean) in response to MECO warming. This may suggest that increased biosiliceous sediment accumulation, if indeed a widespread phenomenon, resulted from higher continental silicate weathering rates and an increase in silicic acid supply to the oceans over several 100 kyr during the MECO.

40Ar-39Ar dating and isotopic composition of groundmass and felspar samples from core AND-2A

The AND-2A drillcore (Antarctic Drilling Program-ANDRILL) was successfully completed in late 2007 on the Antarctic continental margin (Southern McMurdo Sound, Ross Sea) with the aim of tracking ice proximal to shallow marine environmental fluctuations and to document the 20-Ma evolution of the Erebus Volcanic Province. Lava clasts and tephra layers from the AND-2A drillcore were investigated from a petrographic and stratigraphic point of view and analyzed by the 40Ar-39Ar laser technique in order to constrain the age model of the core and to gain information on the style and nature of sediment deposition in the Victoria Land Basin since Early Miocene. Ten out of 17 samples yielded statistically robust 40Ar-39Ar ages, indicating that the AND-2A drillcore recovered <230 m of Middle Miocene (~128-358 m below sea floor, ~11.5-16.0 Ma) and >780 m of Early Miocene (~358-1093 m below sea floor, ~16.0-20.1 Ma). Results also highlight a nearly continuous stratigraphic record from at least 358 m below sea floor down hole, characterized by a mean sedimentation rate of ~19 cm/ka, possible oscillations of no more than a few hundreds of ka and a break within ~17.5-18.1 Ma. Comparison with available data from volcanic deposits on land, suggests that volcanic rocks within the AND-2A core were supplied from the south, possibly with source areas closer to the drill site for the upper core levels, and from 358 m below sea floor down hole, with the 'proto-Mount Morning' as the main source.

Compositional and isotopic data on core MD84641

The organic carbon isotopic record of the sapropels(S1 and S3-S10) and intercalated marl oozes has been determined in a 12-m piston core from the eastern Mediterranean. The d13C_organic values are systematically lighter (mean=-21.0±0.82 per mil) in all sapropels and heavier (mean=-18.8±1.07 per mil) in the marl oozes. These differences are not due to variable marine and terrestrial organic matter mixtures because all values are heavier than modern plankton in the Mediterranean, there is no relationship between the C_organic/N ratios and the isotopic values, and published information on the abundance and distribution of organic biomarkers shows that terrestrial material constitutes a minor fraction of the total organic matter. Temperature effects on isotope fractionation are also discounted because the change in d13C_organic values between glacial and interglacial horizons is in the opposite sense. Diagenesis, which can produce relatively small changes in the carbon isotopic composition of sedimentary organic matter under certain circumstances, is unlikely to have caused the observed differences because this mechanism would cause an enrichmet in 12C, implying that all values were even heavier originally, and there is no secular trend in the d13C_organic record. The observed differences in d13C_organic between the two lithologies are probably produced by changes in the isotopic composition and the concentration of dissolved CO2. First, freshwater flooding during the formation of the sapropels caused the isotopic composition of the dissolved inorganic carbon in the surface waters of the Mediterranean to become lighter because of the 13C deficiency in fresh waters. Hence photosynthesis would have produced isotopically lighter organic material. Second, changes in atmospheric pCO2 between glacial and interglacial periods, as shown by the Vostok ice core, caused marked changes in the concentration of free dissolved CO2 in the mixed layer; lower values during glacial maxima caused a smaller fractionation of the carbon isotopes by phytoplankton, whereas levels were less limiting during the interglacials. Concentrations of dissolved CO2 could also have been much higher during the deposition of the sapropels because of the supply of regenerated CO2 to the mixed layer by upwelling, and this could have further lightened the d13C_organic values in the sapropels themselves. Carbon isotope records may provide an alternative method for estimating atmospheric pCO2 levels over longer time periods than can be obtained from ice cores.