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.

Mineralogy and geochemistry of latest Miocene deposits in the Tyrrhenian Sea

Late Miocene sediments from ODP Sites 652 and 654, drilled on the Sardinian margin in the Western Tyrrhenian Sea, are investigated through mineralogical, micromorphological, geochemical, and microgeochemical analyses. Clay associations appear to be little controlled by conditions of deposition, and largely depend on pre- and post-depositional conditions. The sedimentary series from Central Mediterranean gives very different geodynamic information, according to the sector considered. While relatively stable conditions, like those encountered in Caltanissetta Basin, Sicily, favor the mineralogical expression of warm-temperate and subarid Messinian climate, the Eastern Sardinia margin (Site 654) clay suites mainly reflect the transition from tectonically active to relaxed conditions. The series deposited at the foot of the same margin above a thinner crust (Site 652) experienced the effects of burial diagenesis, enhanced by strong geothermal gradient.

Texture analyses of ODP Site 131-808 sediments

The microfabric of 11 mudrock specimens from ODP Site 808 (Nankai accretionary prism) was quantitatively analyzed using X-ray texture goniometry and optical petrography. The objectives of the study were to learn about rock strain and to detect a component of bulk lateral shortening in the deformation of the mudstones. Strain evaluation is based on the predictions of March theory, and on distortions of initially homogeneous marker particle distributions (the Fry technique). The main results are as follows. The specimens underwent a strain path of progressive flattening, which is closely related to loss of pore space by vertical loading. A component of bulk lateral shortening is detectable in the top 550 mbsf at Site 808, but compared with the amount of uniaxial vertical shortening, its relative magnitude is probably small. Moreover, it cannot be said with confidence whether this is caused by toe contraction of the accretionary wedge or by gravitationally induced downslope movement of the sediment pile. The mudstones examined were deposited in a marine environment with an oxic bottom water column. Micropore collapse is an important fabric building mechanism, but below 400 mbsf its effects are at least partly overridden by recrystallization of smectite. We conclude that mud microfabrics are not very precise deformation gauges, but can be used for rough estimations of strain.