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.

Model simulation output data (CCSM3) related to uplift of African topography and the intensification of upwelling in the Benguela region

The Benguela Current, located off the west coast of southern Africa, is tied to a highly productive upwelling system**1. Over the past 12 million years, the current has cooled, and upwelling has intensified**2, 3, 4. These changes have been variously linked to atmospheric and oceanic changes associated with the glaciation of Antarctica and global cooling**5, the closure of the Central American Seaway**1, 6 or the further restriction of the Indonesian Seaway**3. The upwelling intensification also occurred during a period of substantial uplift of the African continent**7, 8. Here we use a coupled ocean-atmosphere general circulation model to test the effect of African uplift on Benguela upwelling. In our simulations, uplift in the East African Rift system and in southern and southwestern Africa induces an intensification of coastal low-level winds, which leads to increased oceanic upwelling of cool subsurface waters. We compare the effect of African uplift with the simulated impact of the Central American Seaway closure9, Indonesian Throughflow restriction10 and Antarctic glaciation**11, and find that African uplift has at least an equally strong influence as each of the three other factors. We therefore conclude that African uplift was an important factor in driving the cooling and strengthening of the Benguela Current and coastal upwelling during the late Miocene and Pliocene epochs.

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

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

Depth profiles of pore-water and solid-phase constituents, respiration rates and grain size distribution in sediments of the Clarion-Clipperton Fracture Zone

Manganese nodules of the Clarion-Clipperton Fracture Zone (CCFZ) in the NE Pacific Ocean are highly enriched in Ni, Cu, Co, Mo and rare-earth elements, and thus may be the subject of future mining operations. Elucidating the depositional and biogeochemical processes that contribute to nodule formation, as well as the respective redox environment in both, water column and sediment, supports our ability to locate future nodule deposits and evaluates the potential ecological and environmental effects of future deep-sea mining. For these purposes we evaluated the local hydrodynamics and pore-water geochemistry with respect to the nodule coverage at four sites in the eastern CCFZ. Furthermore, we carried out selective leaching experiments at these sites in order to assess the potential mobility of Mn in the solid phase, and compared them with the spatial variations in sedimentation rates. We found that the oxygen penetration depth is 180 - 300 cm at all four sites, while reduction of Mn and NO3- is only significant below the oxygen penetration depth at sites with small or no nodules on the sediment surface. At the site without nodules, potential microbial respiration rates, determined by incubation experiments using 14C-labelled acetate, are slightly higher than at sites with nodules. Leaching experiments showed that surface sediments covered with big or medium-sized nodules are enriched in mobilizable Mn. Our deep oxygen measurements and pore-water data suggest that hydrogenetic and oxic-diagenetic processes control the present-day nodule growth at these sites, since free manganese from deeper sediments is unable to reach the sediment surface. We propose that the observed strong lateral contrasts in nodule size and abundance are sensitive to sedimentation rates, which in turn, are controlled by small-scale variations in seafloor topography and bottom-water current intensity.

Antarctic ice sheet topography, cavity geometry, and global bathymetry (RTopo 1.0.5-beta)

Sub-ice shelf circulation and freezing/melting rates in ocean general circulation models depend critically on an accurate and consistent representation of cavity geometry. Existing global or pan-Antarctic data sets have turned out to contain various inconsistencies and inaccuracies. The goal of this work is to compile independent regional fields into a global data set. We use the S-2004 global 1-minute bathymetry as the backbone and add an improved version of the BEDMAP topography for an area that roughly coincides with the Antarctic continental shelf. Locations of the merging line have been carefully adjusted in order to get the best out of each data set. High-resolution gridded data for upper and lower ice surface topography and cavity geometry of the Amery, Fimbul, Filchner-Ronne, Larsen C and George VI Ice Shelves, and for Pine Island Glacier have been carefully merged into the ambient ice and ocean topographies. Multibeam survey data for bathymetry in the former Larsen B cavity and the southeastern Bellingshausen Sea have been obtained from the data centers of Alfred Wegener Institute (AWI), British Antarctic Survey (BAS) and Lamont-Doherty Earth Observatory (LDEO), gridded, and again carefully merged into the existing bathymetry map. The global 1-minute dataset (RTopo-1 Version 1.0.5) has been split into two netCDF files. The first contains digital maps for global bedrock topography, ice bottom topography, and surface elevation. The second contains the auxiliary maps for data sources and the surface type mask. A regional subset that covers all variables for the region south of 50 deg S is also available in netCDF format. Datasets for the locations of grounding and coast lines are provided in ASCII format.

Grain size distribution in surface sediments samples of the continental shelf and slope of southeast Greenland

During "Meteor" Cruise 6/1966 in the northwest Atlantic a systematic survey of the bottom topography of the southeast Greenland continental margin was undertaken. Eighty-seven profiles transverse to the shelf edge at distances of 3-4 nautical miles and two longitudinal profiles parallel to the coast were carried out with the ELAC Narrow Beam Echo-Sounder giving a reliable record of even steep slopes. On the basis of the echo soundings the topography and morphology of the continental shelf and slope are evaluated. A detailed bathymetric chart and a serial profile chart were designed as working material for the morphological research. These maps along with the original echograms are morphometrically evaluated. The analysis of the sea bottom features is the basis of a subsequent morphogenetical interpretation, verified and extended by means of interpretation of magnetic data and sediment analysis (grain size, roundness, lithology). The results of the research are expressed in a geomorphological map. The primary findings can be summarized as follows: 1) The southeast Greenland shelf by its bottom topography can be clearly designated as a glacially formed area. The glacial features of the shelf can be classified into two zones nearly parallel to the coast: glacial erosion forms on the inner shelf and glacial accumulation forms on the outer shelf. The inner shelf is characterized by the rugged and hummocky topography of ice scoured plains with clear west/east slope asymmetry. On the outer shelf three types of glacial accumulation forms can be recognized: ice margin deposits with clearly expressed terminal moraines, glacial till plains and glaciomarine outwash fans. Both zones of the shelf can be subdivided into two levels of relief. The ice scoured plains, with average depths of 240 meters (m), are dissected to a maximum depth of 1060 m (Gyldenloves Trough) by trough valleys, which are the prolongations of the Greenland fjords. The banks of the outer shelf, with an average depth of 180 m, surround glacial basins with a maximum depth of 670 meters. 2) The sediments of the continental shelf can be classified as glacial due to their grain size distribution and the degree of roundness of the gravel particles. The ice margin deposits on the outer shelf can be recognized by their high percentage of gravels. On the inner shelf a rock surface is suggested, intermittently covered by glacial deposits. In the shelf troughs fine-grained sediments occur mixed with gravels. 3) Topography and sediments show that the southeast Greenland shelf was covered by an ice sheet resting on the sea floor during the Pleistocene ice-age. The large end moraines along the shelf edge probably indicate the maximum extent of the Wurm shelf ice resting on the sea floor. The breakthroughs of the end moraines in front of the glacial basins suggest that the shelf ice has floated further seaward over the increasing depths. 4) Petrographically the shelf sediments consist of gneisses, granites and basalts. While gneisses and granites occire on the nearby coast, basalt is not known to exist here. Either this material has been drifted by icebergs from the basalt province to the north or exists on the southeast Greenland shelf itself. The last interpretation is supported bythe high portion of basalt contained in the sediment samples taken and the strong magnetic anomalies probably caused by basaltic intrusions. 5) A magnetic profile allows the recognition of two magnetically differing areas which approximately coincide with the glacial erosion and accumulation zones. The inner shelf shows a strong and variable magnetic field because the glacially eroded basement forms the sea floor. The outer shelf is characterized by a weak and homogenous magnetic field, as the magnetized basement lies at greater depthy, buried by a thick cover of glacial sediments. The strong magnetic anomalies of the inner shelf are probably caused by dike swarms, similar to those observed further to the north in the Kangerdlugssuaq Fjord region. This interpretation is supported by the high basalt content of the sediment samples and the rough topography of the ice scoured plains which correlates in general with the magnetic fluctuations. The dike structures of the basement have been differentially eroded by the shelf ice. 6) The continental slope, extending from the shelf break at 313 m to a depth of 1270 m with an average slope of 11°, is characterized by delta-shaped projections in front of the shelf basins, by marginal plateaus, ridges and hills, by canyons and slumping features. The projections could be identified as glaciomarine sediment fans. This conclusion is supported by the strong decrease of magnetic field intensity. The deep sea hills and ridges with their greater magnetic intensities have to be regarded as basement outcrops projecting through the glaciomarine sediment cover. The upper continental rise, sloping seaward at about 2°, is composed of wide sediment fans and slump material. A marginal depression on the continental rise running parallel to the shelf edge has been identified. In this depression bottom currents capable of erosion have been recorded. South of Cape Farvel the depression extends to the accumulation zone of the "Eirik" sedimentary ridge. 7) By means of a study of the recent marine processes, postglacial modification of the ice-formed relief can be postulated. The retention effect of the fjord troughs and the high velocity of the East Greenland stream prevents the glacial features from being buried by sediments. Bottom currents capable of active erosion have only been found in the marginal depression on the continental rise. In addition, at the time of the lowest glacio-eustatic sea level, the shelf bottom was not situated in the zone of wave erosion. Only on the continental slope and rise bottom currents, sediment slumps and turbidity currents have led to significant recent modifications. Considering these results, the geomorphological development of the southeast Greenland continental terrace can be suggested as follows: 1. initial formation of a "peneplain", 2. fluvial incision, 3. submergence, and finally 4. glacial modification.

Stable carbon and hydrogen isotope record of n-Alkanoic acids of sediment core SO189/2_144KL

The Indo-Pacific Warm Pool (IPWP) is a key site for the global hydrologic cycle, and modern observations indicate that both the Indian Ocean Zonal Mode (IOZM) and the El Niño Southern Oscillation exert strong influence on its regional hydrologic characteristics. Detailed insight into the natural range of IPWP dynamics and underlying climate mechanisms is, however, limited by the spatial and temporal coverage of climate data. In particular, long-term (multimillennial) precipitation patterns of the western IPWP, a key location for IOZM dynamics, are poorly understood. To help rectify this, we have reconstructed rainfall changes over Northwest Sumatra (western IPWP, Indian Ocean) throughout the past 24,000 y based on the stable hydrogen and carbon isotopic compositions (dD and d13C, respectively) of terrestrial plant waxes. As a general feature of western IPWP hydrology, our data suggest similar rainfall amounts during the Last Glacial Maximum and the Holocene, contradicting previous claims that precipitation increased across the IPWP in response to deglacial changes in sea level and/or the position of the Intertropical Convergence Zone. We attribute this discrepancy to regional differences in topography and different responses to glacioeustatically forced changes in coastline position within the continental IPWP. During the Holocene, our data indicate considerable variations in rainfall amount. Comparison of our isotope time series to paleoclimate records from the Indian Ocean realm reveals previously unrecognized fluctuations of the Indian Ocean precipitation dipole during the Holocene, indicating that oscillations of the IOZM mean state have been a constituent of western IPWP rainfall over the past ten thousand years.

Stable isotopes of foraminifers, mineralogy and ages of sediment core GeoB1523-1

Today the western tropical Atlantic is the most important passage for cross-equatorial transfer of heat in the form of warm surface water flowing from the South into the North Atlantic. Circulation changes north of South America may thus have influenced the global thermohaline circulation system and high northern latitude climate. Here we reconstruct late Quaternary variations of western equatorial Atlantic surface circulation and Amazon lowland climate obtained from a multiproxy sediment record from Ceará Rise. Variations in the illite/smectite ratio suggest drier climatic conditions in the Amazon Basin during glacials relative to interglacials. The 230Thex-normalized fluxes and the 13C/12C record of organic carbon indicate that sea level fluctuations, shelf topography, and changes of the surface circulation pattern controlled variations and amplitude of terrigenous sediment supply to the Ceará Rise. We attribute variations in thermocline depth, reconstructed from vertical planktic foraminiferal oxygen isotope gradients and abundances of the phytoplankton species Florisphaera profunda, to changes in southeast trade wind intensity. Strong trade winds during ice volume maxima are associated with a deep western tropical Atlantic thermocline, strengthening of the North Brazil Current retroflection, and more vigorous eastward flow of surface waters.

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.