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.

Age model and color measurements of sediment cores from the eastern equatorial Pacific

High-resolution records of glacial-interglacial variations in biogenic carbonate, opal, and detritus (derived from non-destructive core log measurements of density, P-wave velocity and color; r >= 0.9) from 15 sediment sites in the eastern equatorial (sampling resolution is ~1 kyr) clear response to eccentricity and precession forcing. For the Peru Basin, we generate a high-resolution (21 kyr increment) orbitally-based chronology for the last 1.3 Ma. Spectral analysis indicates that the 100 kyr cycle became dominant at roughly 1.2 Ma, 200-300 kyr earlier than reported for other paleoclimatic records. The response to orbital forcing is weaker since the Mid-Brunhes Dissolution Event (at 400 ka). A west-east reconstruction of biogenic sedimentation in the Peru Basin (four cores; 91-85°W) distinguishes equatorial and coastal upwelling systems in the western and eastern sites, respectively. A north-south reconstruction perpendicular to the equatorial upwelling system (11 cores, 11°N-°3S) shows high carbonate contents (>= 50%) between 6°N and 4°S and highly variable opal contents between 2°N and 4°S. Carbonate cycles B-6, B-8, B-10, B-12, B-14, M-2, and M-6 are well developed with B-10 (430 ka) as the most prominent cycle. Carbonate highs during glacials and glacial-interglacial transitions extended up to 400 km north and south compared to interglacial or interglacial^glacial carbonate lows. Our reconstruction thus favors glacial-interglacial expansion and contraction of the equatorial upwelling system rather than shifting north or south. Elevated accumulation rates are documented near the equator from 6°N to 4°S and from 2°N to 4°S for carbonate and opal, respectively. Accumulation rates are higher during glacials and glacial-interglacial transitions in all cores, whereas increased dissolution is concentrated on Peru Basin sediments close to the carbonate compensation depth and occurred during interglacials or interglacial-glacial transitions.

ODP Site 162-980 age model and color reflectance record of sediments from ODP Leg 162 sites

Sediment spectral reflectance measurements were generated aboard the JOIDES Resolution during Ocean Drilling Program Leg 162 shipboard operations. The large size of the raw data set (over 1.3 gigabytes) and limited computer hard disk storage space precluded detailed analysis of the data at sea, although broad band averages were used as aids in developing splices and determining lithologic boundaries. This data report describes the methods used to collect these data and their shipboard and postcruise processing. These initial results provide the basis for further postcruise research.

Age model and Plio-Pleistocene benthic stable oxygen isotope ratios of ODP Site 184-1143 in the Sourh China Sea

Based on benthic foraminiferal delta18O from ODP Site 1143, a 5-Myr astronomical timescale for the West Pacific Plio-Pleistocene was established using an automatic orbital tuning method. The tuned Brunhes/Matuyama paleomagnetic polarity reversal age agrees well with the previously published age of 0.78 Ma. The tuned ages for several planktonic foraminifer bio-events also agree well with published dates, and new ages for some other bio-events in the South China Sea were also estimated. The benthic delta18O from Site 1143 is highly coherent with the Earth's orbit (ETP) both at the obliquity and precession bands for the last 5 Myr, and at the eccentricity band for the last 2 Myr. In general, the 41-kyr cycle was dominant through the Plio-Pleistocene although the 23-kyr cycle was also very strong. The 100-kyr cycle became dominant only during the last 1 Myr. A comparison of the benthic delta18O between the Atlantic (ODP 659) and the East and West Pacific (846 and 1143) reveals that the Atlantic-Pacific benthic oxygen isotope difference ratio (Delta delta18OAtl-Pac) displays an increasing trend in three time intervals: 3.6-2.7 Ma, 2.7-2.1 Ma and 1.5-0.25 Ma. Each of the intervals begins with a rapid negative shift in Delta delta18OAtl-Pac, followed by a long period with an increasing trend, corresponding to the growth of the Northern Hemisphere ice sheet. This means that all three intervals of ice sheet growth in the Northern Hemisphere were accompanied at the beginning by a rapid relative warming of deep water in the Atlantic as compared to that of the Pacific, followed by its gradual relative cooling. This general trend, superimposed on the frequent fluctuations with glacial cycles, should yield insights into the processes leading to the boreal glaciation. Cross-spectral analyses of the Delta delta18OAtl-Pac with the Earth's orbit suggests that after the initiation of Northern Hemisphere glaciation at about 2.5 Ma, obliquity rather than precession had become the dominant force controlling the vertical structure or thermohaline circulation in the paleo-ocean.

Stable isotope record and age model for sapropel stratigraphy of ODP Site 160-967

Stable oxygen isotope data from four holes drilled at the Ocean Drilling Program Site 967, which is located on the lower northern slope of the Eratosthenes Seamount, provide a continuous record of Eastern Mediterranean surface-water conditions during the last 3.2 Ma. A high-resolution stratigraphy for the Pliocene-Pleistocene sequence was established by using a combination of astronomical calibration of sedimentary cycles, nannofossil stratigraphy, and stable oxygen isotope fluctuations. Sapropels and color cycles are present throughout the last 3.2 Ma at Site 967, and their ages, as determined by calibration against the precessional component of the astronomical record, are consistent with those estimated for the sapropels of the classical land-based marine sequences of the Punta Piccola, San Nicola, Singa, and Vrica sections (southern Italy). The Site 967 oxygen isotope record shows large amplitude fluctuations mainly caused by variations in surface water salinity throughout the entire period. Spectral analysis shows that fluctuations in the d18O record were predominantly influenced by orbital obliquity and precessional forcing from 3.2 to 1 Ma, and all main orbital frequencies characterize the d18O record for the last million years. The start of sapropel formation at 3.2 Ma indicates a possible link between sapropel formation and the build up of northern hemisphere ice sheets. The dominance of the obliquity cycle in the interval from 3.2-1 Ma further points to the sensitivity of Eastern Mediterranean climate to the fluctuations in the volume of Arctic ice sheets. An intensification of negative isotope anomalies at Site 967, relative to the open ocean, supports a link between high run-off (during warm periods) and sapropel formation. freshwater input would have inhibited deep-water formation, which led to stagnation of deeper waters. Comparison with the land sections also confirms that differential preservation and diagenesis play a key role in sapropel occurrence.

Benthic cover map of Heron Reef derived from a high-spatial-resolution multi-spectral satellite image using object based image analysis

Providing accurate maps of coral reefs where the spatial scale and labels of the mapped features correspond to map units appropriate for examining biological and geomorphic structures and processes is a major challenge for remote sensing. The objective of this work is to assess the accuracy and relevance of the process used to derive geomorphic zone and benthic community zone maps for three western Pacific coral reefs produced from multi-scale, object-based image analysis (OBIA) of high-spatial-resolution multi-spectral images, guided by field survey data. Three Quickbird-2 multi-spectral data sets from reefs in Australia, Palau and Fiji and georeferenced field photographs were used in a multi-scale segmentation and object-based image classification to map geomorphic zones and benthic community zones. A per-pixel approach was also tested for mapping benthic community zones. Validation of the maps and comparison to past approaches indicated the multi-scale OBIA process enabled field data, operator field experience and a conceptual hierarchical model of the coral reef environment to be linked to provide output maps at geomorphic zone and benthic community scales on coral reefs. The OBIA mapping accuracies were comparable with previously published work using other methods; however, the classes mapped were matched to a predetermined set of features on the reef.

Sea-suface microlayer (SML) and underlying water (ULW) samples focusing on CDOM spectral characteristics during METEOR cruise M91

The coastal upwelling system off the coast of Peru is characterized by high biological activity and a pronounced subsurface oxygen minimum zone, as well as associated emissions of atmospheric trace gases such as N2O, CH4 and CO2. From 3 to 23 December 2012, R/V Meteor (M91) cruise took place in the Peruvian upwelling system between 4.59 and 15.4°S, and 82.0 to 77.5°W. During M91 we investigated the composition of the sea-surface microlayer (SML), the oceanic uppermost boundary directly subject to high solar radiation, often enriched in specific organic compounds of biological origin like chromophoric dissolved organic matter (CDOM) and marine gels. In the SML, the continuous photochemical and microbial recycling of organic matter may strongly influence gas exchange between marine systems and the atmosphere. We analyzed SML and underlying water (ULW) samples at 38 stations focusing on CDOM spectral characteristics as indicator of photochemical and microbial alteration processes. CDOM composition was characterized by spectral slope (S) values and excitation-emission matrix fluorescence (EEMs), which allow us to track changes in molecular weight (MW) of DOM, and to determine potential DOM sources and sinks. Spectral slope S varied between 0.012 to 0.043 1 nm-1 and was quite similar between SML and ULW, with no significant differences between the two compartments. Higher S values were observed in the ULW of the southern stations below 15°S. By EEMs, we identified five fluorescent components (F1-5) of the CDOM pool, of which two had excitation/emission characteristics of amino-acid-like fluorophores (F1, F4) and were highly enriched in the SML, with a median ratio SML : ULW of 1.5 for both fluorophores. In the study region, values for CDOM absorption ranged from 0.07 to 1.47 m-1. CDOM was generally highly concentrated in the SML, with a median enrichment with respect to the ULW of 1.2. CDOM composition and changes in spectral slope properties suggested a local microbial release of DOM directly in the SML as a response to light exposure in this extreme environment. In a conceptual model of the sources and modifications of optically active DOM in the SML and underlying seawater (ULW), we describe processes we think may take place (Fig. 1); the production of CDOM of higher MW by microbial release through growth, exudation and lysis in the euphotic zone, includes the identified fluorophores (F1, F2, F3, F4, F5). Specific amino-acid-like fluorophores (F1, F4) accumulate in the SML with respect to the ULW, as photochemistry may enhance microbial CDOM release by (a) photoprotection mechanisms and (b) cell-lysis processes. Microbial and photochemical degradation are potential sinks of the amino-acid-like fluorophores (F1, F4), and potential sources of reworked and more refractory humic-like components (F2, F3, F5). In the highly productive upwelling region along the Peruvian coast, the interplay of microbial and photochemical processes controls the enrichment of amino-acid-like CDOM in the SML. We discuss potential implications for air-sea gas exchange in this area.