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.

Improved geoid solution for the Weddell Sea region

We present a geoid solution for the Weddell Sea and adjacent continental Antarctic regions. There, a refined geoid is of interest, especially for oceanographic and glaciological applications. For example, to investigate the Weddell Gyre as a part of the Antarctic Circumpolar Current and, thus, of the global ocean circulation, the mean dynamic topography (MDT) is needed. These days, the marine gravity field can be inferred with high and homogeneous resolution from altimetric height profiles of the mean sea surface. However, in areas permanently covered by sea ice as well as in coastal regions, satellite altimetry features deficiencies. Focussing on the Weddell Sea, these aspects are investigated in detail. In these areas, ground-based data that have not been used for geoid computation so far provide additional information in comparison with the existing high-resolution global gravity field models such as EGM2008. The geoid computation is based on the remove-compute-restore approach making use of least-squares collocation. The residual geoid with respect to a release 4 GOCE model adds up to two meters and more in the near-coastal and continental areas of the Weddell Sea region, also in comparison with EGM2008. Consequently, the thus refined geoid serves to compute new estimates of the regional MDT and geostrophic currents.

(Table 2) Atmospheric electric measurements during the Atlantic Expedition 1965 of METEOR

During the Atlantic expedition potential gradient, small ion density and space charge density have been recorded. Laborious efforts have been taken for receiving an exact estimation of the reduction factor for the field measurements. The mean value of the potential gradient on the free Atlantic Ocean was 105 V/m. The mean daily course is in very good agreement with the results of the Carnegie Institution. Even records taken on individual days near the quator show this course. For the first time it has been attempted to correlate the potential gradient at sea and the voltage between ionosphere and earth measured over land. A narrow relation has been found in 10 cases of balloon ascents with radiosondes. A further remarkable result is, that the short periodical fluctuations of the air electric field at sea with periods of 2 to 20 minutes have amplitudes of the magnitude of the mean field strength and exist all over the oceans. Recordings of the space charge density show, that positively charged air parcels drift in the first hectometer of the air near the sea surface and produce the fluctuation of the potential gradient. A period analysis did not indicate a recognizable relation to the wind velocity up to now, although an effect of air turbulence must be involved. The concentration of small ions also has been measured occasionally. With this and mean values of the potential gradient the air earth curent density has been computed. With n+ = 310 cm**-3, n- = 220 cm**-3 the air conductivity would be Lambda = 1,14 * 10**-14 Ohm**-1 m**-1. These values are smaller than values of other authors by a factor of 2 or 3. Therefore the computed air earth current density is also smaller. The discrepancy could not be explained yet.

Radiocarbon ages, grain sizes, paleomagentic properties and XRF data of Lake Nam Co sediments, Tibetan Plateau

In September 2008 several cores (68 cm-115 cm length) (water depth: 93 m) were retrieved from Lake Nam Co (southern-central Tibetan Plateau; 4718 m a.s.l.). This study focuses on the interpretation of high-resolution (partly 0.2 cm) data from three gravity cores and the upper part of a 10.4 m long piston core, i.e., the past 4000 cal BP in terms of lake level changes, hydrological variations in the catchment area and consequently variations in monsoon strength. A wide spectrum of sedimentological, geochemical and mineralogical investigations was carried out. Results are presented for XRF core-scans, grain size distribution, XRD-measurements and SEM-image analyses. These data are complemented by an age-depth model using 210Pb and 137Cs analyses as well as eleven AMS-14C-ages. This model is supported by excellent agreement between secular variations determined on one of the gravity cores to geomagnetic field models. This is a significant improvement of the chronology as most catchments of lacustrine systems on the Tibetan Plateau contain carbonates resulting in an unknown reservoir effect for radiocarbon dates. The good correlation of our record to the geomagnetic field models confirms our age-depth model and indicates only insignificant changes in the reservoir effect throughout the last 4 ka. High (summer-) monsoonal activity, i.e. moist environmental conditions, was detected in our record between approximately 4000 and 1950 cal BP as well as between 1480 and 1200 cal BP. Accordingly, lower monsoon activity prevails in periods between the two intervals and thereafter. This pattern shows a good correlation to the variability of the Indian Ocean Summer Monsoon (IOSM) as recorded in a peat bog ~1000 km in NE direction from Lake Nam Co. This is the first time that such a supra regional homogenous monsoon activity is shown on the Tibetan Plateau and beyond. Finally our data show a significant lake level rise after the Little Ice Age (LIA) in Lake Nam Co which is suggested to be linked to glacier melting in consequence of rising temperatures occurring on the whole Tibetan Plateau during this time.

(Table 2) Radiocarbon dates of sediment cores obtained during Bio Hespérides cruise SVAIS and OGS Explora cruise EGLACOM between 2007-2008

A high-resolution paleomagnetic and rock magnetic study has been carried out on sediment cores collected in glaciomarine silty-clay sequences from the continental shelf and slope of the southern Storfjorden trough-mouth fan, on the northwestern Barents Sea continental margin. The Storfjorden sedimentary system was investigated during the SVAIS and EGLACOM cruises, when 10 gravity cores, with a variable length from 1.03 m to 6.41 m, were retrieved. Accelerator mass spectrometry (AMS) 14C analyses on 24 samples indicate that the cores span a time interval that includes the Holocene, the last deglaciation phase and in some cores the last glacial maximum. The sediments carry a well-defined characteristic remanent magnetization and have a valuable potential to reconstruct the paleosecular variation (PSV) of the geomagnetic field, including relative paleointensity (RPI) variations. The paleomagnetic data allow reconstruction of past dynamics and amplitude of the geomagnetic field variations at high northern latitudes (75°-76° N). At the same time, the rock magnetic and paleomagnetic data allow a high-resolution correlation of the sedimentary sequences and a refinement of their preliminary age models. The Holocene PSV and RPI records appear particularly sound, since they are consistent between cores and they can be correlated to the closest regional stacking curves (UK PSV, FENNOSTACK and FENNORPIS) and global geomagnetic model for the last 7 ka (CALS7k.2). The computed amplitude of secular variation is lower than that outlined by some geomagnetic field models, suggesting that it has been almost independent from latitude during the Holocene.

Age models, sedimentology and geochemistry from sediment cores in Drake Passage throughflow

The Drake Passage (DP) is the major geographic constriction for the Antarctic Circumpolar Current (ACC) and exerts a strong control on the exchange of physical, chemical, and biological properties between the Atlantic, Pacific, and Indian Ocean basins. Resolving changes in the flow of circumpolar water masses through this gateway is, therefore, crucial for advancing our understanding of the Southern Ocean's role in global ocean and climate variability. Here, we reconstruct changes in DP throughflow dynamics over the past 65,000 y based on grain size and geochemical properties of sediment records from the southernmost continental margin of South America. Combined with published sediment records from the Scotia Sea, we argue for a considerable total reduction of DP transport and reveal an up to ~40% decrease in flow speed along the northernmost ACC pathway entering the DP during glacial times. Superimposed on this long-term decrease are high-amplitude, millennial-scale variations, which parallel Southern Ocean and Antarctic temperature patterns. The glacial intervals of strong weakening of the ACC entering the DP imply an enhanced export of northern ACC surface and intermediate waters into the South Pacific Gyre and reduced Pacific-Atlantic exchange through the DP ("cold water route"). We conclude that changes in DP throughflow play a critical role for the global meridional overturning circulation and interbasin exchange in the Southern Ocean, most likely regulated by variations in the westerly wind field and changes in Antarctic sea ice extent.

Age models and summer sea surface temperature and winter sea ice concentration for the EPILOG-LGM time slice in the Pacific Southern Ocean

Sea surface temperatures and sea-ice extent are the most critical variables to evaluate the Southern Ocean paleoceanographic evolution in relation to the development of the global carbon cycle, atmospheric CO2 variability and ocean-atmosphere circulation. In contrast to the Atlantic and the Indian sectors, the Pacific sector of the Southern Ocean has been insufficiently investigated so far. To cover this gap of information we present diatom-based estimates of summer sea surface temperature (SSST) and winter sea-ice concentration (WSI) from 17 sites in the polar South Pacific to study the Last Glacial Maximum (LGM) at the EPILOG time slice (19,000-23,000 cal. years BP). Applied statistical methods are the Imbrie and Kipp Method (IKM) and the Modern Analog Technique (MAT) to estimate temperature and sea-ice concentration, respectively. Our data display a distinct LGM east-west differentiation in SSST and WSI with steeper latitudinal temperature gradients and a winter sea-ice edge located consistently north of the Pacific-Antarctic Ridge in the Ross sea sector. In the eastern sector of our study area, which is governed by the Amundsen Abyssal Plain, the estimates yield weaker latitudinal SSST gradients together with a variable extended winter sea-ice field. In this sector, sea-ice extent may have reached sporadically the area of the present Subantarctic Front at its maximum LGM expansion. This pattern points to topographic forcing as major controller of the frontal system location and sea-ice extent in the western Pacific sector whereas atmospheric conditions like the Southern Annular Mode and the ENSO affected the oceanographic conditions in the eastern Pacific sector. Although it is difficult to depict the location and the physical nature of frontal systems separating the glacial Southern Ocean water masses into different zones, we found a distinct temperature gradient in latitudes straddled by the modern Southern Subtropical Front. Considering that the glacial temperatures north of this zone are similar to the modern, we suggest that this represents the Glacial Southern Subtropical Front (GSSTF), which delimits the zone of strongest glacial SSST cooling (>4K) to its North. The southern boundary of the zone of maximum cooling is close to the glacial 4°C isotherm. This isotherm, which is in the range of SSST at the modern Antarctic Polar Front (APF), represents a circum-Antarctic feature and marks the northern edge of the glacial Antarctic Circumpolar Current (ACC). We also assume that a glacial front was established at the northern average winter sea ice edge, comparable with the modern Southern Antarctic Circumpolar Current Front (SACCF). During the glacial, this front would be located in the area of the modern APF. The northward deflection of colder than modern surface waters along the South American continent leads to a significant cooling of the glacial Humboldt Current surface waters (4-8K), which affects the temperature regimes as far north as into tropical latitudes. The glacial reduction of ACC temperatures may also result in the significant cooling in the Atlantic and Indian Southern Ocean, thus may enhance thermal differentiation of the Southern Ocean and Antarctic continental cooling. Comparison with temperature and sea ice simulations for the last glacial based on numerical simulations show that the majority of modern models overestimate summer and winter sea ice cover and that there exists few models that reproduce our temperature data rather well.