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.

Code and readme of a MATLAB-based graphical user interface program for computing functionals of the geopotential

We present a novel graphical user interface program GrafLab (GRAvity Field LABoratory) for spherical harmonic synthesis (SHS) created in MATLAB®. This program allows to comfortably compute 38 various functionals of the geopotential up to ultra-high degrees and orders of spherical harmonic expansion. For the most difficult part of the SHS, namely the evaluation of the fully normalized associated Legendre functions (fnALFs), we used three different approaches according to required maximum degree: (i) the standard forward column method (up to maximum degree 1800, in some cases up to degree 2190); (ii) the modified forward column method combined with Horner's scheme (up to maximum degree 2700); (iii) the extended-range arithmetic (up to an arbitrary maximum degree). For the maximum degree 2190, the SHS with fnALFs evaluated using the extended-range arithmetic approach takes only approximately 2-3 times longer than its standard arithmetic counterpart, i.e. the standard forward column method. In the GrafLab, the functionals of the geopotential can be evaluated on a regular grid or point-wise, while the input coordinates can either be read from a data file or entered manually. For the computation on a regular grid we decided to apply the lumped coefficients approach due to significant time-efficiency of this method. Furthermore, if a full variance-covariance matrix of spherical harmonic coefficients is available, it is possible to compute the commission errors of the functionals. When computing on a regular grid, the output functionals or their commission errors may be depicted on a map using automatically selected cartographic projection.

Medem Channel profiles based on stand-alone model simulation and corrections with nudging

State-of-the-art process-based models have shown to be applicable to the simulation and prediction of coastal morphodynamics. On annual to decadal temporal scales, these models may show limitations in reproducing complex natural morphological evolution patterns, such as the movement of bars and tidal channels, e.g. the observed decadal migration of the Medem Channel in the Elbe Estuary, German Bight. Here a morphodynamic model is shown to simulate the hydrodynamics and sediment budgets of the domain to some extent, but fails to adequately reproduce the pronounced channel migration, due to the insufficient implementation of bank erosion processes. In order to allow for long-term simulations of the domain, a nudging method has been introduced to update the model-predicted bathymetries with observations. The model-predicted bathymetry is nudged towards true states in annual time steps. Sensitivity analysis of a user-defined correlation length scale, for the definition of the background error covariance matrix during the nudging procedure, suggests that the optimal error correlation length is similar to the grid cell size, here 80-90 m. Additionally, spatially heterogeneous correlation lengths produce more realistic channel depths than do spatially homogeneous correlation lengths. Consecutive application of the nudging method compensates for the (stand-alone) model prediction errors and corrects the channel migration pattern, with a Brier skill score of 0.78. The proposed nudging method in this study serves as an analytical approach to update model predictions towards a predefined 'true' state for the spatiotemporal interpolation of incomplete morphological data in long-term simulations.

Presence-absence species/sample matrix for the measured sections at the area of Maghara, North Sinai, Egypt

As age-diagnostic fossils are rare in the Middle to Upper Jurassic sedimentary succession of Gebel Maghara, North Sinai, Egypt, and in order to ensure maximal stratigraphic resolution, chronostratigraphic boundaries were determined based on quantitative biostratigraphy. A data matrix comprising 231 macrofaunal taxa in 93 samples from four sections has been processed with the Unitary Association (UA) Method. This led to construction of a sequence of 29 UAs (maximal sets of actually or virtually coexisting taxa), which have been grouped into 14 laterally reproducible association zones. The UA method allowed an in-depth analysis of the stratigraphically conflicting taxa, enabled the biostratigraphic subdivision of the studied interval, and also provided stratigraphic correlation among the measured sections and with the Tethyan ammonite zones.

(Table 2) Trace element composition of seleceted clast and matrix materials from ODP Site 125-786

Shipboard examination of volcanic and sedimentary strata at Site 786 suggested that at least four types of breccias are present: flow-top breccias, associated with cooling and breakup on the upper surface of lava flows; autobreccias, formed by in-situ alteration at the base of flows; fault-gouge breccias; and true sedimentary breccias derived from weathering and erosion of underlying flows. It is virtually impossible to assess the origin of breccia matrix by textural and mineralogical analyses alone. However, it is fundamental for our understanding of breccia provenance to determine the source component of the matrix material. Whether the matrix is uniquely clastderived can be determined by geochemical fingerprinting. Trace elements that are immobile during weathering and alteration do not change their relative abundances. A contribution to the matrix from any source with an immobile trace element signature different from that of the clasts would appear as a perturbation of the trace element signature of the matrix. Trace element analysis of bulk samples from clasts and matrix material in individual breccia units was undertaken in a fashion similar to that used by Brimhall and Dietrich (1987, doi:10.1016/0016-7037(87)90070-6) in analyzing soil provenance: (1) to help distinguish between sedimentary and volcanic breccias, (2) to determine the degree of mixing and depth of erosion in sedimentary breccias, and (3) to analyze the local provenance of the individual breccia components (matrix and clasts). The following elements were analyzed by X-ray fluorescence (XRF): Rb, Sr, Ba, U, Zr, Cu, Zn, Ti, Cr, and V. Of these elements, Zr and Ti probably exhibit truly immobile behavior (Humphris and Thompson, 1978, doi:10.1016/0016-7037(78)90222-3 ). The remaining elements are useful as a reference for the extent of compositional change during the formation of matrix material (Brimhall and Dietrich, 1987, doi:10.1016/0016-7037(87)90070-6).

Composition of matrix glasses and melt inclusions of fallout tephra from ODP Hole 125-782A

We studied the systematics of Cl, F and H2O in Izu arc front volcanic rocks using basaltic through rhyolitic glass shards and melt inclusions (Izu glasses) from Oligocene to Quaternary distal fallout tephra. These glasses are low-K basalts to rhyolites that are equivalent to the Quaternary lavas of the Izu arc front (Izu VF). Most of the Izu glasses have Cl ~400-4000 ppm and F ~70-400 ppm (normal-group glasses). Rare andesitic melt inclusions (halogen-rich andesites; HRA) have very high abundances of Cl (~6600-8600 ppm) and F (~780-910 ppm), but their contents of incompatible large ion lithophile elements (LILE) are similar to the normal-group glasses. The preeruptive H2O of basalt to andesite melt inclusions in plagioclase is estimated to range from ~2 to ~10 wt% H2O. The Izu magmas should be undersaturated in H2O and the halogens at their preferred levels of crystallization in the middle to lower crust (~3 to ~11 kbar, ~820° to ~1200°C). A substantial portion of the original H2O is lost due to degassing during the final ascent to surface. By contrast, halogen loss is minor, except for loss of Cl from siliceous dacitic and rhyolitic compositions. The behavior of Cl, F and H2O in undegassed melts resembles the fluid mobile LILE (e.g.; K, Rb, Cs, Ba, U, Pb, Li). Most of the Cl (>99%), H2O (>95%) and F (>53%) in the Izu VF melts appear to originate from the subducting slab. At arc front depths, the slab fluid contains Cl = 0.94+/-0.25 wt%, F = 990+/-270 ppm and H2O = 25+/-7 wt%. If the subducting sediment and the altered basaltic crust were the only slab sources, then the subducted Cl appears to be almost entirely recycled at the Izu arc (~77-129%). Conversely, H2O (~13-22% recycled at arc) and F (~4-6% recycled) must be either lost during shallow subduction or retained in the slab to greater depths. If a seawater-impregnated serpentinite layer below the basaltic crust were an additional source of Cl and H2O, the calculated percentage of Cl and H2O recycled at arc would be lower. Extrapolating the Izu data to the total length of global arcs (~37000 km), the global arc outflux of fluid-recycled Cl and H2O at subduction zones amounts to Cl ~2.9-3.8 mln ton/yr and H2O ~70-100 mln ton/yr, respectively - comparable to previous estimates. Further, we obtain a first estimate of global arc outflux of fluid-recycled F of ~0.3-0.4 mln ton/yr. Despite the inherent uncertainties, our results support models suggesting that the slab becomes strongly depleted in Cl and H2O in subduction zones. In contrast, much of the subducted F appears to be returned to the deep mantle, implying efficient fractionation of Cl and H2O from F during the subduction process. However, if slab devolatilization produces slab fluids with high Cl/F (~9.5), slab melting will still produce components with low Cl/F ratios (~0.9), similar to those characteristic of the upper continental crust (Cl/F ~0.3-0.9).