Age model and grain size distribution of the fine fraction of ODP Site 114-704

The physical properties of sediments beneath an upwelling area in the southern part of the Atlantic Ocean (ODP Hole 704A) were investigated. Highly significant correlations characterize the relationship of carbonate content to bulk density (R = 0.85), carbonate content to porosity (R = 0.84), and carbonate content to impedance (R = 0.84). No relationship exists between carbonate content and compressional-wave velocity (R = 0.24), indicating that amplitude variations in impedance are primarily controlled by variations in bulk density, which, in turn, are controlled by climatically driven biogenic opal and carbonate deposition. In general, maxima in impedance correspond to maxima in carbonate content (minima in opal content). The impedance record exhibits its most drastic change at about 2.4 Ma, marking dramatic increases in the average content of biogenic opal and the beginning of large-amplitude fluctuations. Between 0.7 and 0.4 Ma carbonate content, bulk density, and grain density decrease while opal content drastically increases. Similar changes have been observed in sediments beneath an upwelling cell off northwest Africa, indicating an oceanwide enhancement in upwelling or in the calcite corrosiveness of bottom water that appears to be isochronous.

Glacial-interglacial record of ODP Site 181-1119 in the Canterbury Bright, Southwest Pacific Ocean

Ocean Drilling Program (ODP) Site 1119 is located at water depth 395 m near the subtropical front (STF; here represented by the Southland Front), just downslope from the shelf edge of eastern South Island, New Zealand. The upper 86.19 metres composite depth (mcd) of Site 1119 sediment was deposited at an average sedimentation rate of 34 cm/kyr during Marine Isotope Stages (MIS) 1-8 (0-252 ka), and is underlain across a ~25 kyr intra-MIS 8 unconformity by MIS 8.5-11 (277-367 ka) and older sediment deposited at ~14 cm/kyr. A time scale is assigned to Site 1119 using radiocarbon dates for the period back to ~39 ka, and, prior to then, by matching its climatic record with that of the Vostok ice core, which it closely resembles. Four palaeoceanographic proxy measures for surface water masses vary together with the sandy-muddy, glacial-interglacial (G/I) cyclicity at the site. Interglacial intervals are characterised by heavy delta13C, high colour reflectance (a proxy for carbonate content), low Q-ray (a proxy for clay content) and light delta18O; conversely, glacial intervals exhibit light delta13C, low reflectance, high Q-ray and heavy delta18O signatures. Early interglacial intervals are represented by silty clays with 10-105-cm-thick beds of sharp-based (Chondrites-burrowed), shelly, graded, fine sand. The sands are rich in foraminifera, and were deposited distant from the shoreline under the influence of longitudinal flow in relatively deep water. Glacial intervals comprise mostly micaceous silty clay, though with some thin (2-10 cm thick) sands present also at peak cold periods, and contain the cold-water scallop Zygochlamys delicatula. Interglacial sandy intervals are characterised by relatively low sedimentation rates of 5-32 cm/kyr; cold climate intervals MIS 10, 6 and 2 have successively higher sedimentation rates of 45, 69 and 140 cm/kyr. Counter-intuitively,and forced by the bathymetric control of a laterally-moving shoreline during G/I and I/G transitions, the 1119 core records a southeasterly (seaward) movement of the STF during early glacial periods, accompanied by the incursion of subtropical water (STW) above the site, and northwesterly (landward) movement during late glacial and interglacial times, resulting in a dominant influence then of subantarctic surface water (SAW). The history of passage of these different water masses at the site is clearly delineated by their characteristic delta13C values. The intervals of thin, graded sands-muds which occur within MIS 2-3, 6, 7.4 and 10 indicate the onset at times of peak cold of intermittent bottom currents caused by strengthened and expanded frontal flows along the STF, which at such times lay near Site 1119 in close proximity to seaward-encroaching subantarctic waters within the Bounty gyre. In common with other nearby Southern Hemisphere records, the cold period which represents the last glacial maximum lasted between ~23-18 ka at Site 1119, during which time the STF and Subantarctic Front (SAF) probably merged into a single intense frontal zone around the head of the adjacent Bounty Trough.

Core alignment and composite depth scale of DSDP Site 86-577

Deep Sea Drilling Project Site 577 on Shatsky Rise (North Pacific Ocean) recovered a series of cores at three holes that contain calcareous nannofossil ooze of latest Cretaceous (late Maastrichtian) through early Eocene age. Several important records have been generated using samples from these cores, but the stratigraphy has remained outdated and confusing. Here we revise the stratigraphy at Site 577. This includes refining several age datums, realigning cores in the depth domain, and placing all stratigraphic markers on a current time scale. The work provides a template for appropriately bringing latest Cretaceous and Paleogene data sets at old drill sites into current paleoceanographic literature for this time interval. While the Paleocene Eocene Thermal Maximum (PETM) lies within core gaps at Holes 577* and 577A, the sedimentary record at the site holds other important events and remains crucially relevant to understanding changes in oceanographic conditions from the latest Cretaceous through early Paleogene.

Modelling calving front dynamics using a level-set method: application to Jakobshavn Isbræ, West Greenland, links model results in MATLAB format

Calving is a major mechanism of ice discharge of the Antarctic and Greenland ice sheets, and a change in calving front position affects the entire stress regime of marine terminating glaciers. The representation of calving front dynamics in a 2-D or 3-D ice sheet model remains non-trivial. Here, we present the theoretical and technical framework for a level-set method, an implicit boundary tracking scheme, which we implement into the Ice Sheet System Model (ISSM). This scheme allows us to study the dynamic response of a drainage basin to user-defined calving rates. We apply the method to Jakobshavn Isbræ, a major marine terminating outlet glacier of the West Greenland Ice Sheet. The model robustly reproduces the high sensitivity of the glacier to calving, and we find that enhanced calving triggers significant acceleration of the ice stream. Upstream acceleration is sustained through a combination of mechanisms. However, both lateral stress and ice influx stabilize the ice stream. This study provides new insights into the ongoing changes occurring at Jakobshavn Isbræ and emphasizes that the incorporation of moving boundaries and dynamic lateral effects, not captured in flow-line models, is key for realistic model projections of sea level rise on centennial timescales.

Age model for ODP Leg 175 sites

Upwelling along the western coast of Africa south of the equator may be partitioned into three major areas, each having its own dynamics and history: (1) the eastern equatorial region, comprising the Congo Fan and the area of Mid-Angola; (2) the Namibia upwelling system, extending from the Walvis Ridge to Lüderitz; and (3) the Cape Province region, where upwelling is subdued. The highest nutrient contents in thermocline waters are in the northern region, the lowest in the southern one. Wind effects are at a maximum near the southern end of the Namibia upwelling system, and maximum productivity occurs near Walvis Bay, where the product between upwelling rate and nutrient content of upwelled waters is at a maximum. In the Congo/Angola region, opal tends to follow organic carbon quite closely in the Quaternary record. However, organic carbon has a strong precessional component, while opal does not. Despite relatively low opal content, sediments off Angola show the same patterns as those off the Congo; thus, they are part of the same regime. The spectrum shows nonlinear interference patterns between high- and low-latitude forcing, presumably tied to thermocline fertility and wind. On Walvis Ridge, as in the Congo-Angola region, the organic matter record behaves normally; that is, supply is high during glacial periods. In contrast, interglacial periods are favorable for opal deposition. The pattern suggests reduction in silicate content of the thermocline during glacial periods. The reversed phase (opal abundant during interglacials) persists during the entire Pleistocene and can be demonstrated deep into the Pliocene, not just on Walvis Ridge but all the way to the Oranje River and off the Cape Province. From comparison with other regions, it appears that silicate is diminished in the global thermocline, on average, whenever winds become strong enough to substantially shorten the residence time of silicate in upper waters (Walvis Hypothesis, solving the Walvis Paradox of reversed phase in opal deposition). The central discovery during Leg 175 was the documentation of a late Pliocene opal maximum for the entire Namibia upwelling system (early Matuyama Diatom Maximum [MDM]). The maximum is centered on the period between the end of the Gauss Chron and the beginning of the Olduvai Chron. A rather sharp increase in both organic matter deposition and opal deposition occurs near 3 Ma in the middle of the Gauss Chron, in association with a series of major cooling steps. As concerns organic matter, high production persists at least to 1 Ma, when there are large changes in variability, heralding subsequent pulsed production periods. From 3 to 2 Ma, organic matter and opal deposition run more or less parallel, but after 2 Ma opal goes out of phase with organic matter. Apparently, this is the point when silicate becomes limiting to opal production. Thus, the MDM conundrum is solved by linking planetary cooling to increased mixing and upwelling (ramping up to the MDM) and a general removal of silicate from the upper ocean through excess precipitation over global supply (ramping down from the MDM). The hypothesis concerning the origin of the Namibia opal acme or MDM is fundamentally the same as the Walvis Hypothesis, stating that glacial conditions result in removal of silicate from the thermocline (and quite likely from the ocean as a whole, given enough time). The Namibia opal acme, and other opal maxima in the latest Neogene in other regions of the ocean, marks the interval when a cooling ocean selectively removes the abundant silicate inherited from a warm ocean. When the excess silicate is removed, the process ceases. According to the data gathered during Leg 175, major upwelling started in the late part of the late Miocene. Presumably, this process contributed to the drawing down of carbon dioxide from the atmosphere, helping to prepare the way for Northern Hemisphere glaciation.

Average fluorescence and dissolved iron and Fe-binding ligand characteristics during POLARSTERN expedition ANT-XXIV/3 in 2008

Organic complexation of dissolved iron (dFe) was investigated in the Atlantic sector of the Southern Ocean in order to understand the distribution of Fe over the whole water column. The total concentration of dissolved organic ligands ([Lt]) measured by voltammetry ranged between 0.54 and 1.84 nEq of M Fe whereas the conditional binding strength (K') ranged between 10**21.4 and 10**22.8. For the first time, trends in Fe-organic complexation were observed in an ocean basin by examining the ratio ([Lt]/[dFe]), defined as the organic ligand concentration divided by the dissolved Fe concentration. The [Lt]/[dFe] ratio indicates the saturation state of the natural ligands with Fe; a ratio near 1 means saturation of the ligands leading to precipitation of Fe. Reversely, high ratios mean Fe depletion and show a high potential for Fe solubilisation. In surface waters where phytoplankton is present low dissolved Fe and high variable ligand concentrations were found. Here the [Lt]/[dFe] ratio was on average 4.4. It was especially high (5.6-26.7) in the HNLC (High Nutrient, Low Chlorophyll) regions, where Fe was depleted. The [Lt]/[dFe] ratio decreased with depth due to increasing dissolved Fe concentrations and became constant below 450 m, indicating a steady state between ligand and Fe. Relatively low [Lt]/[dFe] ratios (between 1.1 and 2.7) existed in deep water north of the Southern Boundary, facilitating Fe precipitation. The [Lt]/[dFe] ratio increased southwards from the Southern Boundary on the Zero Meridian and from east to west in the Weddell Gyre due to changes both in ligand characteristics and in dissolved iron concentration. High [Lt]/[dFe] ratio expresses Fe depletion versus ligand production in the surface. The decrease with depth reflects the increase of [dFe] which favours scavenging and (co-) precipitation, whereas a horizontal increase in the deep waters results from an increasing distance from Fe sources. This increase in the [Lt]/[dFe] ratio at depth shows the very resistant nature of the dissolved organic ligands.

Tide amplitude at Ilulissat, and surface height and flow velocities of Jakobshavn Isbrae in August 2004

During the summer of 2004, the front area of the Jakobshavn Isbræ was monitored using a geodetic-photogrammetric survey with temporarily coincident precise observations of local ocean tides in the Disko Bay close to Ilulissat. The geodetic and photogrammetric observations were conducted at the southern margin of the glacier front. The largest observed horizontal flow velocities are in the central part of the front with values up to 45 m/d. This is a factor of 2 greater than the average velocities at the front area observed in the last century. Our new observations confirm previous estimates of an acceleration of glacier flow during the last decade. The photogrammetric survey provided flow trajectories for 4000 surface points with a time resolution of 30 min. These flow trajectories were used to compare the vertical motion of the glacier with the observed tides. The existence of a free-floating glacier tongue in 2004 was confirmed by these data. However, it occupied only a small belt, of at most a few 100 m width, in the central part of the glacier front. Horizontal motion did not appear to depend on the tidal phase, unlike some of the fast-moving ice streams of West Antarctica.

Model of interactive visualization of diffusion processes from different strains, and source routes and spatial diffusion of capuchin monkeys over the past 6 million years

Source routes and Spatial Diffusion of capuchin monkeys over the past 6 million years, rebuilt in the SPREAD 1.0.6 from the MCC tree. The map shows the 10 different regions to which distinctive samples were associated. The different transmission routes have been calculated from the average rate over time. Only rates with Bayes factor> 3 were considered as significantly different from zero. Significant diffusion pathways are highlighted with color varying from dark brown to red, being the dark brown less significant rates and deep red the most significant rates.

Topographic model of Greenland (elevation lines in PanMap and ASCII format)

In this paper, a new high-resolution elevation model of Greenland, including the ice sheet as well as the ice free regions, is presented. It is the first published full coverage model, computed with an average resolution of 2 km and providing an unprecedented degree of detail. The topography is modeled from a wide selection of data sources, including satellite radar altimetry from Geosat and ERS 1, airborne radar altimetry and airborne laser altimetry over the ice sheet, and photogrammetric and manual map scannings in the ice free region. The ice sheet model accuracy is evaluated by omitting airborne laser data from the analysis and treating them as ground truth observations. The mean accuracy of the ice sheet elevations is estimated to be 12-13 m, and it is found that on surfaces of a slope between 0.2° and 0.8°, corresponding to approximately 50% of the ice sheet, the model presents a 40% improvement over models based on satellite altimetry alone. On coastal bedrock, the model is compared with stereo triangulated reference points, and it is found that the model accuracy is of the order of 25-35 m in areas covered by stereo photogrammetry scannings and between 200 and 250 m elsewhere.

(Table 1) Age model of sediment core MD03-2705

A reconstruction of northwest African summer monsoon strength during the cold marine isotopic stage (MIS) 6 indicates a link to the seasonal migration of the Intertropical Convergence Zone (ITCZ). High-resolution studies of eolian dust supply and sea surface temperature recorded in marine core MD03-2705, on the Mauritanian margin, provide a better understanding about the penultimate glacial history of northwestern African aridity/humidity and upwelling coastal activity. Today, site MD03-2705 experiences increased upwelling and dust flux during the winter months, when the ITCZ is in a southerly position. Analyses of foraminifera isotopic composition suggest that during MIS 6.5 (180-168 ka) the average position of the ITCZ migrated north, marked by an increase in the strength of the summer monsoon, which decreased eolian dust transport and the coastal upwelling activity. The northward migration is in phase with a specific orbital combination of a low precessional index with a high obliquity signal. High-resolution analysis of stable isotopes (d18O and d13C) and microscale resolution geochemical (Ti/Al and quartz grain counts) determinations reveal that the transition between monsoonal humid (MIS 6.5) and dry (MIS 6.4) conditions has occurred in less than 1.3 ka. Such rapid changes suggest a nonlinear link between the African monsoonal rainfall system and environmental changes over the continent. This study provides new insights about the influence of vegetation and oceanic temperature feedbacks on the onset of African summer monsoon and demonstrates that, during the penultimate glacial period, changes in tropical dynamics had regional and global impacts.

Results of analysis and model simulation of the bed profiles

Bedforms such as dunes and ripples are ubiquitous in rivers and coastal seas, and commonly described as triangular shapes from which height and length are calculated to estimate hydrodynamic and sediment dynamic parameters. Natural bedforms, however, present a far more complicated morphology; the difference between natural bedform shape and the often assumed triangular shape is usually neglected, and how this may affect the flow is unknown. This study investigates the shapes of natural bedforms and how they influence flow and shear stress, based on four datasets extracted from earlier studies on two rivers (the Rio Paraná in Argentina, and the Lower Rhine in The Netherlands). The most commonly occurring morphological elements are a sinusoidal stoss side made of one segment and a lee side made of two segments, a gently sloping upper lee side and a relatively steep (6 to 21°) slip face. A non-hydrostatic numerical model, set up using Delft3D, served to simulate the flow over fixed bedforms with various morphologies derived from the identified morphological elements. Both shear stress and turbulence increase with increasing slip face angle and are only marginally affected by the dimensions and positions of the upper and lower lee side. The average slip face angle determined from the bed profiles is 14°, over which there is no permanent flow separation. Shear stress and turbulence above natural bedforms are higher than above a flat bed but much lower than over the often assumed 30° lee side angle.

Sea ice thickness distribution (ITD) Model for the Arctic, links to NetCDF files

Sea ice models contain many different parameterizations of which one of the most commonly used is a subgrid-scale ice thickness distribution (ITD). The effect of this model component and the associated ice strength formulation on the reproduction of observed Arctic sea ice is assessed. To this end the model's performance in reproducing satellite observations of sea ice concentration, thickness and drift is evaluated. For an unbiased comparison, different model configurations with and without an ITD are tuned with an automated parameter optimization. The original combination of ITD and ice strength parameterization does not lead to better results than a simple single category model. Yet changing to a simpler ice strength formulation, which depends linearly on the mean ice thickness across all thickness categories, allows to clearly improve the model-data misfit when using an ITD. In the original formulation, the ice strength depends strongly on the number of thickness categories, so that introducing more categories can lead to thicker albeit weaker ice on average.