Isotopic composition of terrestrial plant waxes (dD and d13C of n-alkanes) of sediment core SO188_342 from the Northern Bay of Bengal for the last 18 ka

The Indian Summer Monsoon (ISM) is a major global climatic phenomenon. Long-term precipitation proxy records of the ISM, however, are often fragmented and discontinuous, impeding an estimation of the magnitude of precipitation variability from the Last Glacial to the present. To improve our understanding of past ISM variability, we provide a continuous reconstructed record of precipitation and continental vegetation changes from the lower Ganges-Brahmaputra-Meghna catchment and the Indo-Burman ranges over the last 18,000 years (18 ka). The records derive from a marine sediment core from the northern Bay of Bengal (NBoB), and are complemented by numerical model results of spatial moisture transport and precipitation distribution over the Bengal region. The isotopic composition of terrestrial plant waxes (dD and d13C of n-alkanes) are compared to results from an isotope-enabled general atmospheric circulation model (IsoCAM) for selected time slices (pre-industrial, mid-Holocene and Heinrich Stadial 1). Comparison of proxy and model results indicate that past changes in the dD of precipitation and plant waxes were mainly driven by the amount effect, and strongly influenced by ISM rainfall. Maximum precipitation is detected for the Early Holocene Climatic Optimum (EHCO; 10.5-6 ka BP), whereas minimum precipitation occurred during the Heinrich Stadial 1 (HS1; 16.9-15.4 ka BP). The IsoCAM model results support the hypothesis of a constant moisture source (i.e. the NBoB) throughout the study period. Relative to the pre-industrial period the model reconstructions show 20% more rain during the mid-Holocene (6 ka BP) and 20% less rain during the Heinrich Stadial 1 (HS1), respectively. A shift from C4-plant dominated ecosystems during the glacial to subsequent C3/C4-mixed ones during the interglacial took place. Vegetation changes were predominantly driven by precipitation variability, as evidenced by the significant correlation between the dD and d13C alkane records. When compared to other records across the ISM domain, precipitation and vegetation changes inferred from our records and the numerical model results provide evidence for a coherent regional variability of the ISM from the Last Glacial to the present.

Geochemistry and Mg-isotopes of carbonate sediments and porefluids of ODP Site 130-807

Deep-sea pore fluids are potential archives of ancient seawater chemistry. However, the primary signal recorded in pore fluids is often overprinted by diagenetic processes. Recent studies have suggested that depth profiles of Mg concentration in deep-sea carbonate pore fluids are best explained by a rapid rise in seawater Mg over the last 10-20 Myr. To explore this possibility we measured the Mg isotopic composition of pore fluids and carbonate sediments from Ocean Drilling Program (ODP) site 807. Whereas the concentration of Mg in the pore fluid declines with depth, the isotopic composition of Mg in the pore fluid increases from -0.78 per mil near the sediment-water interface to -0.15 per mil at 778 mbsf. The Mg isotopic composition of the sediment, with few important exceptions, does not change with depth and has an average d26Mg value of -4.72 per mil. We reproduce the observed changes in sediment and pore-fluid Mg isotope values using a numerical model that incorporates Mg, Ca and Sr cycling and satisfies existing pore-fluid Ca isotope and Sr data. Our model shows that the observed trends in magnesium concentrations and isotopes are best explained as a combination of two processes: a secular rise in the seawater Mg over the Neogene and the recrystallization of low-Mg biogenic carbonate to a higher-Mg diagenetic calcite. These results indicate that burial recrystallization will add Mg to pelagic carbonate sediments, leading to an overestimation of paleo-temperatures from measured Mg/Ca ratios. The Mg isotopic composition of foraminiferal calcite appears to be only slightly altered by recrystallization making it possible to reconstruct the Mg isotopic composition of seawater through time.

Decade Wind and Wave Data for NW Iberia, Cabo Sillerio

Seafloor sediment mobilization on the inner Northwest Iberian continental shelf is caused largely by ocean surface waves. The temporal and spatial variability in the wave height, wave period, and wave direction has a profound effect on local sediment mobilization, leading to distinct sediment mobilization scenarios. Six grain-size specific sediment mobilization scenarios, representing seasonal average and storm conditions, were simulated with a physics-based numerical model. Model inputs included meteorological and oceanographic data in conjunction with seafloor grain-size and the shelf bathymetric data. The results show distinct seasonal variations, most importantly in wave height, leading to sediment mobilization, specifically on the inner shelf shallower than 30 m water depth where up to 49% of the shelf area is mobilized. Medium to severe storm events are modeled to mobilize up to 89% of the shelf area above 150 m water depth. The frequency of each of these seasonal and storm-related sediment mobilization scenarios is addressed using a decade of meteorological and oceanographic data. The temporal and spatial patterns of the modeled sediment mobilization scenarios are discussed in the context of existing geological and environmental processes and conditions to assist scientific, industrial and environmental efforts that are directly affected by sediment mobilization. Examples, where sediment mobilization plays a vital role, include seafloor nutrient advection, recurrent arrival of oil from oil-spill-laden seafloor sediment, and bottom trawling impacts.

DISPRO datasets for validation of coastal hydrodynamic models

Appropriate field data are required to check the reliability of hydrodynamic models simulating the dispersion of soluble substances in the marine environment. This study deals with the collection of physical measurements and soluble tracer data intended specifically for this kind of validation. The intensity of currents as well as the complexity of topography and tides around the Cap de La Hague in the center of the English Channel makes it one of the most difficult areas to represent in terms of hydrodynamics and dispersion. Controlled releases of tritium - in the form of HTO - are carried out in this area by the AREVA-NC plant, providing an excellent soluble tracer. A total of 14 493 measurements were acquired to track dispersion in the hours and days following a release. These data, supplementing previously gathered data and physical measurements (bathymetry, water-surface levels, Eulerian and Lagrangian current studies) allow us to test dispersion models from the hour following release to periods of several years which are not accessible with dye experiments. The dispersion characteristics are described and methods are proposed for comparing models against measurements. An application is proposed for a 2 dimensions high-resolution numerical model. It shows how an extensive dataset can be used to build, calibrate and validate several aspects of the model in a highly dynamic and macrotidal area: tidal cycle timing, tidal amplitude, fixed-point current data, hodographs. This study presents results concerning the model's ability to reproduce residual Lagrangian currents, along with a comparison between simulation and high-frequency measurements of tracer dispersion. Physical and tracer data are available from the SISMER database of IFREMER (www.ifremer.fr/sismer/catal). This tool for validation of models in macro-tidal seas is intended to be an open and evolving resource, which could provide a benchmark for dispersion model validation.

Measured ice surface velocities on the King George Island ice cap with DGPS

King George Island is located at the northern tip of the Antarctic Peninsula, which is influenced by maritime climate conditions. The observed mean annual air temperature at sea level is -2.4°C. Thus, the ice cap is regarded as sensitive to changing climatic conditions. Ground-penetrating radar surveys indicate a partly temperate ice cap with an extended water layer at the firn/ice transition of the up to 700 m high ice cap. Measured firn temperatures are close to 0°C at the higher elevations, and they differ considerably from the measured mean annual air temperature. The aim of this paper is to present ice-flow dynamics by means of observations and simulations of the flow velocities. During several field campaigns from 1997/98 to 2008/09, ice surface velocities were derived with repeated differential GPS measurements. Ice velocities vary from 0.7 m/a at the dome to 112.1 m/a along steep slopes. For the western part of the ice cap a three-dimensional diagnostic full-Stokes model was applied to calculate ice flow. Parameters of the numerical model were identified with respect to measured ice surface velocities. The simulations indicate cold ice at higher elevations, while temperate ice at lower elevations is consistent with the observations.

Bedform measurements from the San Francisco Bay Coastal System

The morphology of ~45,000 bedforms from 13 multibeam bathymetry surveys was used as a proxy for identifying net bedload sediment transport directions and pathways throughout the San Francisco Bay estuary and adjacent outer coast. The spatially-averaged shape asymmetry of the bedforms reveals distinct pathways of ebb and flood transport. Additionally, the region-wide, ebb-oriented asymmetry of 5% suggests net seaward-directed transport within the estuarine-coastal system, with significant seaward asymmetry at the mouth of San Francisco Bay (11%), through the northern reaches of the Bay (7-8%), and among the largest bedforms (21% for lambda > 50 m). This general indication for the net transport of sand to the open coast strongly suggests that anthropogenic removal of sediment from the estuary, particularly along clearly defined seaward transport pathways, will limit the supply of sand to chronically eroding, open-coast beaches. The bedform asymmetry measurements significantly agree (up to ~ 76%) with modeled annual residual transport directions derived from a hydrodynamically-calibrated numerical model, and the orientation of adjacent, flow-sculpted seafloor features such as mega-flute structures, providing a comprehensive validation of the technique. The methods described in this paper to determine well-defined, cross-validated sediment transport pathways can be applied to estuarine-coastal systems globally where bedforms are present. The results can inform and improve regional sediment management practices to more efficiently utilize often limited sediment resources and mitigate current and future sediment supply-related impacts.

(Table S1) Stable oxygen isotope record of barite from DSDP Leg 85 sites

Oxygen isotopes in marine sulfate (d18O SO4) measured in marine barite show variability over the past 10 million years, including a 5per mil decrease during the Plio-Pleistocene, with near-constant values during the Miocene that are slightly enriched over the modern ocean. A numerical model suggests that sea level fluctuations during Plio-Pleistocene glacial cycles affected the sulfur cycle by reducing the area of continental shelves and increasing the oxidative weathering of pyrite. The data also require that sulfate concentrations were 10 to 20% lower in the late Miocene than today.

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.

(Table 2) Foehn characteristics of station in the McMurdo Dry Valleys, Antarctica

Foehn winds resulting from topographic modification of airflow in the lee of mountain barriers are frequently experienced in the McMurdo Dry Valleys (MDVs) of Antarctica. Strong foehn winds in the MDVs cause dramatic warming at onset and have significant effects on landscape forming processes; however, no detailed scientific investigation of foehn in the MDVs has been conducted. As a result, they are often misinterpreted as adiabatically warmed katabatic winds draining from the polar plateau. Herein observations from surface weather stations and numerical model output from the Antarctic Mesoscale Prediction System (AMPS) during foehn events in the MDVs are presented. Results show that foehn winds in the MDVs are caused by topographic modification of south-southwesterly airflow, which is channeled into the valleys from higher levels. Modeling of a winter foehn event identifies mountain wave activity similar to that associated with midlatitude foehn winds. These events are found to be caused by strong pressure gradients over the mountain ranges of the MDVs related to synoptic-scale cyclones positioned off the coast of Marie Byrd Land. Analysis of meteorological records for 2006 and 2007 finds an increase of 10% in the frequency of foehn events in 2007 compared to 2006, which corresponds to stronger pressure gradients in the Ross Sea region. It is postulated that the intra- and interannual frequency and intensity of foehn events in the MDVs may therefore vary in response to the position and frequency of cyclones in the Ross Sea region.

Geochemistry of sediments from the Eastern Cape Basin

Sediment cores retrieved in the Benguela coastal upwelling system off Namibia show very distinct enrichments of solid phase barium at the sulfate/methane transition (SMT). These barium peaks represent diagenetic barite (BaSO4) fronts which form by the reaction of upwardly diffusing barium with interstitial sulfate. Calculated times needed to produce these barium enrichments indicate a formation time of about 14,000 yr. Barium spikes a few meters below the SMT were observed at one of the investigated sites (GeoB 8455). Although this sulfate-depleted zone is undersaturated with respect to barite, the dominant mineral phase of these buried barium enrichments was identified as barite by scanning electron microscopy (SEM). This is the first study which reports the occurrence/preservation of pronounced barite enrichments in sulfate-depleted sediments buried a few meters below the SMT. At site GeoB 8455 high concentrations of dissolved barium in pore water as well as barium in the solid phase were observed. Modeling the measured barium concentrations at site GeoB 8455 applying the numerical model CoTReM reveals that the dissolution rate of barite directly below the SMT is about one order of magnitude higher than at the barium enrichments deeper in the sediment core. This indicates that the dissolution of barite at these deeper buried fronts must be retarded. Thus, the occurrence of the enrichments in solid phase barium at site GeoB 8455 could be explained by decreased dissolution rates of barite due to the changes in the concentration of barite in the sediment, as well as changes in the saturation state of fluids. Furthermore, the alteration of barite into witherite (BaCO3) via the transient phase barium sulfide could lead to the preservation of a former barite front as BaCO3. The calculations and modeling indicate that a relocation of the barite front to a shallower depth occurred between the last glacial maxium (LGM) and the Pleistocene/Holocene transition. We suggest that an upward shift of the SMT occurred at that time, most likely as a result of an increase in the methanogenesis rates due to the burial of high amounts of organic matter below the SMT.

Geochemical records of sediment cores from the Sea of Japan

Intervals of organic C- and carbonate-rich laminated sediments occur in the Sea of Japan with roughly the same frequency as temperature changes observed in Greenland ice cores, providing clear evidence of rapid oceanographic change during the past 36 kyr. Planktonic foraminiferal d18O data suggest that only the laminated sediments deposited during the Last Glacial Maximum (LGM), and perhaps one other interval formed during a period of increased water column stratification. Sedimentary Re and Mo data are consistent with bottom waters that were sulfidic during the LGM and suboxic during other laminated intervals. Results of a numerical model of Corg and Re burial are consistent with a mechanism whereby an increased Corg flux to the seafloor drove oxygen concentrations toward depletion during times of deposition of the suboxic laminated intervals. Such a process could have resulted from increased upwelling driven either by increased deep water formation due to colder and/or more saline surface waters or by stronger northeasterly monsoonal winds.

Age and oxygen isotope data for bulk carbonate from DSDP Sites 41-366 and 72-516 and ODP Sites 111-677A and 130-807

A numerical model which describes oxygen isotope exchange during burial and recrystallization of deep-sea carbonate is used to obtain information on how sea surface temperatures have varied in the past by correcting measured d18O values of bulk carbonate for diagenetic overprinting. Comparison of bulk carbonate and planktonic foraminiferal d18O records from ODP site 677A indicates that the oxygen isotopic composition of bulk carbonate does reflect changes in sea surface temperature and d18O. At ODP Site 690, we calculate that diagenetic effects are small, and that both bulk carbonate and planktonic foraminiferal d18O records accurately reflect Paleogene warming of high latitude surface oceans, biased from diagenesis by no more than 1°C. The same is likely to be true for other high latitude sites where sedimentation rates are low. At DSDP sites 516 and 525, the effects of diagenesis are more significant. Measured d18O values of Eocene bulk carbonates are more than 2? lower at deeply buried site 516 than at site 525, consistent with the model prediction that the effects of diagenesis should be proportional to sedimentation rate. Model-corrections reconcile the differences in the data between the two sites; the resulting paleotemperature reconstruction indicates a 4°C cooling of mid-latitude surface oceans since the Eocene. At low latitudes, the contrast in temperature between the ocean surface and bottom makes the carbonate d180 values particularly sensitive to diagenetic effects; most of the observed variations in measured d18O values are accounted for by diagenetic effects rather than by sea surface temperature variations. We show that the data are consistent with constant equatorial sea surface temperatures through most of the Cenozoic, with the possible exception of the early Eocene, when slightly higher temperatures are indicated. We suggest that the lower equatorial sea surface temperatures for the Eocene and Oligocene reported in other oxygen isotope studies are artifacts of diagenetic recrystallization, and that it is impossible to reconstruct accurately equatorial sea surface temperatures without explicitly accounting for diagenetic overprinting.

Calcium isotope ratios of pore fluids and solid phase of organic rich sediments

Pore fluid calcium isotope, calcium concentration and strontium concentration data are used to measure the rates of diagenetic dissolution and precipitation of calcite in deep-sea sediments containing abundant clay and organic material. This type of study of deep-sea sediment diagenesis provides unique information about the ultra-slow chemical reactions that occur in natural marine sediments that affect global geochemical cycles and the preservation of paleo-environmental information in carbonate fossils. For this study, calcium isotope ratios (d44/40Ca) of pore fluid calcium from Ocean Drilling Program (ODP) Sites 984 (North Atlantic) and 1082 (off the coast of West Africa) were measured to augment available pore fluid measurements of calcium and strontium concentration. Both study sites have high sedimentation rates and support quantitative sulfate reduction, methanogenesis and anaerobic methane oxidation. The pattern of change of d44/40Ca of pore fluid calcium versus depth at Sites 984 and 1082 differs markedly from that of previously studied deep-sea Sites like 590B and 807, which are composed of nearly pure carbonate sediment. In the 984 and 1082 pore fluids, d44/40Ca remains elevated near seawater values deep in the sediments, rather than shifting rapidly toward the d44/40Ca of carbonate solids. This observation indicates that the rate of calcite dissolution is far lower than at previously studied carbonate-rich sites. The data are fit using a numerical model, as well as more approximate analytical models, to estimate the rates of carbonate dissolution and precipitation and the relationship of these rates to the abundance of clay and organic material. Our models give mutually consistent results and indicate that calcite dissolution rates at Sites 984 and 1082 are roughly two orders of magnitude lower than at previously studied carbonate-rich sites, and the rate correlates with the abundance of clay. Our calculated rates are conservative for these sites (the actual rates could be significantly slower) because other processes that impact the calcium isotope composition of sedimentary pore fluid have not been included. The results provide direct geochemical evidence for the anecdotal observation that the best-preserved carbonate fossils are often found in clay or organic-rich sedimentary horizons. The results also suggest that the presence of clay minerals has a strong passivating effect on the surfaces of biogenic carbonate minerals, slowing dissolution dramatically even in relation to the already-slow rates typical of carbonate-rich sediments.

Predicting bed form roughness: the influence of lee side angle

Flow transverse bedforms (ripples and dunes) are ubiquitous in rivers and coastal seas. Local hydrodynamics and transport conditions depend on the size and geometry of these bedforms, as they constitute roughness elements at the bed. Bedform influence on flow energy must be considered for the understanding of flow dynamics, and in the development and application of numerical models. Common estimations or predictors of form roughness (friction factors) are based mostly on data of steep bedforms (with angle-of-repose lee slopes), and described by highly simplified bedform dimensions (heights and lengths). However, natural bedforms often are not steep, and differ in form and hydraulic effect relative to idealised bedforms. Based on systematic numerical model experiments, this study shows how the hydraulic effect of bedforms depends on the flow structure behind bedforms, which is determined by the bedform lee side angle, aspect ratio and relative height. Simulations reveal that flow separation behind bedform crests and, thus, a hydraulic effect is induced at lee side angles steeper than 11 to 18° depending on relative height, and that a fully developed flow separation zone exists only over bedforms with a lee side angle steeper than 24°. Furthermore, the hydraulic effect of bedforms with varying lee side angle is evaluated and a reduction function to common friction factors is proposed. A function is also developed for the Nikuradse roughness (k s), and a new equation is proposed which directly relates k s to bedform relative height, aspect ratio and lee side angle.

(Table 1) Physical oceanography from several POLARSTERN cruises to the Arctic Ocean, 1993 - 2012

Freshwater in the Arctic Ocean plays an important role in the regional ocean circulation, sea ice, and global climate. From salinity observed by a variety of platforms, we are able, for the first time, to estimate a statistically reliable liquid freshwater trend from monthly gridded fields over all upper Arctic Ocean basins. From 1992 to 2012 this trend was 600±300 km**3/yr. A numerical model agrees very well with the observed freshwater changes. A decrease in salinity made up about two thirds of the freshwater trend and a thickening of the upper layer up to one third. The Arctic Ocean Oscillation index, a measure for the regional wind stress curl, correlated well with our freshwater time series. No clear relation to Arctic Oscillation or Arctic Dipole indices could be found. Following other observational studies, an increased Bering Strait freshwater import to the Arctic Ocean, a decreased Davis Strait export, and enhanced net sea ice melt could have played an important role in the freshwater trend we observed.