Distribution of benthic foraminifera in lower Cretaceous sediments of the Cuvuer Abyssal Plain in the Indian Ocean

Analysis of 66 samples from DSDP Site 263 (Cores 263-4R-4 to 263-29R-4) reveals a unique faunal composition with a predominance of agglutinated taxa, many of them previously unrecorded from any other DSDP and ODP Indian Ocean sites. A total of 66 agglutinated and 31 calcareous taxa are documented and five new species are described: Hippocrepina gracilis n.sp., "Textuluriopsis" elegans n.sp., Aaptotoichus challengeri n.sp., "Gaudryinopsis" pseudobettenstaedti n.sp. and "Gaudryina" cuvierensis n.sp. Three assemblages are recognized based on changes in the composition of dominant taxa and occurrences of stratigraphically important species: (1) a high-diversity Valanginian to Barremian Bulbobaculites-Recurvoides Assemblage (Cores 263-29R to - 18R), comprised of numerous elongate agglutinated forms, rare nodosariids, and variable numbers of tubes and ammodiscids; (2) a moderately diverse Aptian to Albian Rhizammina-Ammodiscus-Glomospira Assemblage (Cores 263-18R to -7R) with highly fluctuating numbers of the nominate taxa and Haplophragmoides, Trochammina, Verneuilinoides spp., and Verneuilina howchini; (3) a very low diversity Albian or younger Assemblage (Cores 263-6R to -4R) containing sparse agglutinated foraminifera, rare nodosariids and rotaliids. We interpret the assemblages as shelf to lower slope and consider them to reflect a deepening palaeobathymetry as the Cuvier margin subsided after the initial breakup of East Gondwana during the Valanginian. Our interpretation is in sharp contrast with initial palaeodepth estimates of less than 100 m, as well as with original chronostratigraphic interpretations based on foraminifera and nannofossils which correlated the base of the recovered interval with the Aptian. The absence of many cosmoplitan forms, despite high diversity suggests strong faunal differentiation in the Austral realm or endemisn within the Cuvier Basin during the Early Cretaceous.

Age determination and faunal composition of sediments off northern Greenland

We present a record encompassing marine isotope stages 7-1 from a hitherto unexplored and heavily ice-covered area of the Arctic Ocean, the Lomonosov Ridge off the northern Greenland-Canada continental margin, using nannofossil and benthic foraminifera stratigraphy. Planktic foraminifera assemblages are used as a key paleoceanographic proxy, and a surprisingly large variability is found for an interior Arctic Ocean site. Abundant small (63-125 µm) subpolar Turborotalita quinqueloba occur in two sections, possibly representing substages 5e (last interglacial) and 5a (warm interstadial). However, the present-day circulation pattern and the very distant location of high productive regions cannot explain such high abundances of subpolar specimens in the interior, perennially sea ice-covered Arctic Ocean. Hence our proxy record indicates that last interglacial sea ice concentrations were reduced off some areas of northern Greenland-Canada. Whether this was part of a larger regional pattern or it represents the influence of polynya areas with locally increased productivity remains to be solved. With respect to glacial conditions, increased ice-rafted debris (IRD) deposition in the area appears to be associated with glacial stages 6, 4, and late 3. Stage 2 sediments (including the Last Glacial Maximum) are condensed with a sparse IRD content only.

Silicoflagellates in ODP Leg 108 holes

The occurrence of Cenozoic silicoflagellates at three Ocean Drilling Program (ODP) Holes (660A, 662A, and 667A) was investigated to determine biostratigraphic and relative paleotemperature relations in the tropical Atlantic Ocean. This report presents the data obtained from a study of 37 samples and some preliminary comments on the data. The age of the single sparse assemblage at Hole 660A is late middle Eocene or late Eocene (Dictyocha hexacantha Zone); the sparse to common assemblages of Hole 667A are Oligocene and early Miocene and the common to abundant assemblages of Hole 662A are early Pliocene to Quaternary. Dissolution thinning of silicoflagellates is noted in most samples, even in Hole 662A, which is under the present productive Benguela Current.

Neogene palaeoceanographic and palaeoclimatic events inferred from palynological data: Cape Basin off South Africa, ODP Leg 175

Sites 1085, 1086 and 1087 were drilled off South Africa during Ocean Drilling Program (ODP) Leg 175 to investigate the Benguela Current System. While previous studies have focused on reconstructing the Neogene palaeoceanographic and palaeoclimatic history of these sites, palynology has been largely ignored, except for the Late Pliocene and Quaternary. This study presents palynological data from the upper Middle Miocene to lower Upper Pliocene sediments in Holes 1085A, 1086A and 1087C that provide complementary information about the history of the area. Abundant and diverse marine palynomorphs (mainly dinoflagellate cysts), rare spores and pollen, and dispersed organic matter have been recovered. Multivariate statistical analysis of dispersed organic matter identified three palynofacies assemblages (A, B, C) in the most continuous hole (1085A), and they were defined primarily by amorphous organic matter (AOM), and to a lesser extent black debris, structured phytoclasts, degraded phytoclasts, and marine palynomorphs. Ecostratigraphic interpretation based on dinoflagellate cyst, spore-pollen and palynofacies data allowed us to identify several palaeoceanographic and palaeoclimatic signals. First, the late Middle Miocene was subtropical, and sediments contained the highest percentages of land-derived organic matter, even though they are rich in AOM (palynofacies assemblage A). Second, the Late Miocene was cool-temperate and characterized by periods of intensified upwelling, increase in productivity, abundant and diverse oceanic dinoflagellate cysts, and the highest percentages of AOM (palynofacies assemblage C). Third, the Early to early Late Pliocene was warm-temperate with some dry intervals (increase in grass pollen) and intensified upwelling. Fourth, the Neogene "carbonate crash" identified in other southern oceans was recognized in two palynofacies A samples in Hole 1085A that are nearly barren of dinoflagellate cysts: one Middle Miocene sample (590 mbsf, 13.62 Ma) and one Upper Miocene sample (355 mbsf, 6.5 Ma). Finally, the extremely low percentages of pollen suggest sparse vegetation on the adjacent landmass, and Namib desert conditions were already in existence during the late Middle Miocene.

10Be data for sediments from DSDP Leg 56 and ODP Leg 170

Sediment accretion and subduction at convergent margins play an important role in the nature of hazardous interplate seismicity (the seismogenic zone) and the subduction recycling of volatiles and continentally derived materials to the Earth's mantle. Identifying and quantifying sediment accretion, essential for a complete mass balance across the margin, can be difficult. Seismic images do not define the processes by which a prism was built, and cored sediments may show disturbed magnetostratigraphy and sparse biostratigraphy. This contribution reports the first use of cosmogenic 10Be depth profiles to define the origin and structural evolution of forearc sedimentary prisms. Biostratigraphy and 10Be model ages generally are in good agreement for sediments drilled at Deep Sea Drilling Project Site 434 in the Japan forearc, and support an origin by imbricate thrusting for the upper section. Forearc sediments from Ocean Drilling Program Site 1040 in Costa Rica lack good fossil or paleomagnetic age control above the decollement. Low and homogeneous 10Be concentrations show that the prism sediments are older than 3-4 Ma, and that the prism is either a paleoaccretionary prism or it formed largely from slump deposits of apron sediments. Low 10Be in Costa Rican lavas and the absence of frontal accretion imply deeper sediment underplating or subduction erosion.

Late Quaternary siliciclastic/carbonate sedimentation in the Cabricorn Channel

A model is presented for hemipelagic siliciclastic and carbonate sedimentation during the last glacial-interglacial cycle in the Capricorn Channel, southern Great Barrier Reef (GBR). Stable isotope ratios, grainsize, carbonate content and mineralogy were analysed for seven cores in a depth transect from 166 to 2892 m below sea level (mbsl). Results show variations in the flux of terrigenous, neritic and pelagic sediments to the continental slope over the last sea level cycle. During the glacial lowstand terrigenous sediment influenced all the cores down to 2000 mbsl. The percentages of quartz and feldspar in the cores decreased with water depth, while the percentage of clay increased. X-ray diffraction analysis of the glacial lowstand clay mineralogy suggests that the siliciclastic sediment was primarily sourced from the Fitzroy River, which debouched directly into the northwest sector of the Capricorn Channel at this time. The cores also show a decrease in pelagic calcite and an increase in aragonite and high magnesium calcite (HMC) during the glacial. The influx of HMC and aragonite is most likely from reworking of coral reefs exposed on the continental shelf during the glacial, and also from HMC ooids precipitated at the head of the Capricorn Channel at this time. Mass accumulation rates (MARs) are high (13.5 g/cm**/kyr) during the glacial and peak at ~20 g/cm** 3/kyr in the early transgression (16-14 ka BP). MARs then decline with further sea level rise as the Fitzroy River mouth retreats from the edge of the continental shelf after 13.5 ka BP. MARs remain low (4 g/cm**3/kyr) throughout the Holocene highstand. Data for the Holocene highstand indicate there is a reduction in siliciclastic influx to the Capricorn Channel with little quartz and feldspar below 350 mbsl. However, fine-grained fluvial sediments, presumably from the Fitzroy River, were still accumulating on the mid slope down to 2000 mbsl. The proportion of pelagic calcite in the core tops increases with water depth, while HMC decreases, and is present only in trace amounts in cores below 1500 mbsl. The difference in the percentage of HMC in the deeper cores between the glacial and Holocene may reflect differences in supply or deepening of the HMC lysocline during the glacial. Sediment accumulation rates also vary between cores in the Capricorn Channel and do not show the expected exponential decrease with depth. This may be due to intermediate or deep water currents reworking the sediments. It is also possible that present bathymetry data are too sparse to detect the potential role that submarine channels may play in the distribution and accumulation of sediments. Comparison of the Capricorn Channel MARs with those for other mixed carbonate/siliciclastic provinces from the northeast margin of Australia indicates that peak MARs in the early transgression in the Capricorn Channel precede those from the central GBR and south of Fraser Island. The difference in the timing of the carbonate and siliciclastic MAR peaks along the northeast margin is primarily related to differences in the physiography and climate of the provinces. The only common trend in the MARs from the northeast margin of Australia is the near synchronicity of the carbonate and siliciclastic MAR peaks in individual sediment cores, which supports a coeval sedimentation model.

(Table 1) Calcite, siderite, and stable isotopes d13C and d18O at DSDP Hole 93-603B

Diagenesis of the fine-grained, feldspathic sandstones in the Lower Cretaceous submarine fan complex cored in DSDP Hole 603B can be considered to have occurred in three stages: (1) replacement of matrix and framework grains by pyrite, siderite, phillipsite (?), and particularly by ferroan calcite; (2) dissolution of ferroan calcite and feldspars to produce secondary macroporosity; and (3) development of sparse feldspar and quartz overgrowths, and authigenic modification of remnant matrix. Only ferroan calcite is a volumetrically important diagenetic mineral phase (up to 50 vol.%). Matrix in thin sandstone turbidite deposits has been extensively replaced by ferroan calcite. Carbon stable isotope data suggest that organic diagenesis had only a minor influence on calcite precipitation. Oxygen stable isotope data indicate that the minimum average calcite precipitation temperature was 40° C. Preliminary calculations show that steadystate diffusion of Ca+ + from the dissolution of nannoplankton skeletal material in the interbedded pelagic marls to the associated sandstones is a feasible transport mechanism. A thick sandstone unit from 1234-1263 m sub-bottom is extensively replaced by calcite near the upper and lower contacts. Farther into the sand body away from the contacts, the sandstone has good secondary porosity resulting from the dissolution of ferroan calcite that partially replaced matrix and framework grains. The central portion of the thick sand appears to be a channel with high-energy clean sand. We believe that the channel provided a conduit for focused flow of diagenetic compactional fluids responsible for dissolution. Focused flow may be the result of the earlier lithification of the pelagic limestones and thin-bedded sandstones which, then formed vertical permeability barriers. Calcite dissolution has occurred and may still be occurring at temperatures less than 65°C.

On nonlinearity in neural encoding models applied to the primary visual cortex

Within the regression framework, we show how different levels of nonlinearity influence the instantaneous firing rate prediction of single neurons. Nonlinearity can be achieved in several ways. In particular, we can enrich the predictor set with basis expansions of the input variables (enlarging the number of inputs) or train a simple but different model for each area of the data domain. Spline-based models are popular within the first category. Kernel smoothing methods fall into the second category. Whereas the first choice is useful for globally characterizing complex functions, the second is very handy for temporal data and is able to include inner-state subject variations. Also, interactions among stimuli are considered. We compare state-of-the-art firing rate prediction methods with some more sophisticated spline-based nonlinear methods: multivariate adaptive regression splines and sparse additive models. We also study the impact of kernel smoothing. Finally, we explore the combination of various local models in an incremental learning procedure. Our goal is to demonstrate that appropriate nonlinearity treatment can greatly improve the results. We test our hypothesis on both synthetic data and real neuronal recordings in cat primary visual cortex, giving a plausible explanation of the results from a biological perspective.

Learning an L1-regularized Gaussian Bayesian Network in the Equivalence Class Space

Learning the structure of a graphical model from data is a common task in a wide range of practical applications. In this paper, we focus on Gaussian Bayesian networks, i.e., on continuous data and directed acyclic graphs with a joint probability density of all variables given by a Gaussian. We propose to work in an equivalence class search space, specifically using the k-greedy equivalence search algorithm. This, combined with regularization techniques to guide the structure search, can learn sparse networks close to the one that generated the data. We provide results on some synthetic networks and on modeling the gene network of the two biological pathways regulating the biosynthesis of isoprenoids for the Arabidopsis thaliana plant

Lazy Lasso for local regression

Locally weighted regression is a technique that predicts the response for new data items from their neighbors in the training data set, where closer data items are assigned higher weights in the prediction. However, the original method may suffer from overfitting and fail to select the relevant variables. In this paper we propose combining a regularization approach with locally weighted regression to achieve sparse models. Specifically, the lasso is a shrinkage and selection method for linear regression. We present an algorithm that embeds lasso in an iterative procedure that alternatively computes weights and performs lasso-wise regression. The algorithm is tested on three synthetic scenarios and two real data sets. Results show that the proposed method outperforms linear and local models for several kinds of scenarios

Automatic Categorization for Improving Spanish into Spanish Sign Language Machine Translation

This paper describes a preprocessing module for improving the performance of a Spanish into Spanish Sign Language (Lengua de Signos Espanola: LSE) translation system when dealing with sparse training data. This preprocessing module replaces Spanish words with associated tags. The list with Spanish words (vocabulary) and associated tags used by this module is computed automatically considering those signs that show the highest probability of being the translation of every Spanish word. This automatic tag extraction has been compared to a manual strategy achieving almost the same improvement. In this analysis, several alternatives for dealing with non-relevant words have been studied. Non-relevant words are Spanish words not assigned to any sign. The preprocessing module has been incorporated into two well-known statistical translation architectures: a phrase-based system and a Statistical Finite State Transducer (SFST). This system has been developed for a specific application domain: the renewal of Identity Documents and Driver's License. In order to evaluate the system a parallel corpus made up of 4080 Spanish sentences and their LSE translation has been used. The evaluation results revealed a significant performance improvement when including this preprocessing module. In the phrase-based system, the proposed module has given rise to an increase in BLEU (Bilingual Evaluation Understudy) from 73.8% to 81.0% and an increase in the human evaluation score from 0.64 to 0.83. In the case of SFST, BLEU increased from 70.6% to 78.4% and the human evaluation score from 0.65 to 0.82.

Linear instability of orthogonal compressible leading-edge boundary layer flow

Instability analysis of compressible orthogonal swept leading-edge boundary layer flow was performed in the context of BiGlobal linear theory. 1, 2 An algorithm was developed exploiting the sparsity characteristics of the matrix discretizing the PDE-based eigenvalue problem. This allowed use of the MUMPS sparse linear algebra package 3 to obtain a direct solution of the linear systems associated with the Arnoldi iteration. The developed algorithm was then applied to efficiently analyze the effect of compressibility on the stability of the swept leading-edge boundary layer and obtain neutral curves of this flow as a function of the Mach number in the range 0 ≤ Ma ≤ 1. The present numerical results fully confirmed the asymptotic theory results of Theofilis et al. 4 Up to the maximum Mach number value studied, it was found that an increase of this parameter reduces the critical Reynolds number and the range of the unstable spanwise wavenumbers.

The PSE-3D instability analysis methodology for flows depending strongly on two and weakly on the third spatial dimension

The present contribution discusses the development of a PSE-3D instability analysis algorithm, in which a matrix forming and storing approach is followed. Alternatively to the typically used in stability calculations spectral methods, new stable high-order finitedifference-based numerical schemes for spatial discretization 1 are employed. Attention is paid to the issue of efficiency, which is critical for the success of the overall algorithm. To this end, use is made of a parallelizable sparse matrix linear algebra package which takes advantage of the sparsity offered by the finite-difference scheme and, as expected, is shown to perform substantially more efficiently than when spectral collocation methods are used. The building blocks of the algorithm have been implemented and extensively validated, focusing on classic PSE analysis of instability on the flow-plate boundary layer, temporal and spatial BiGlobal EVP solutions (the latter necessary for the initialization of the PSE-3D), as well as standard PSE in a cylindrical coordinates using the nonparallel Batchelor vortex basic flow model, such that comparisons between PSE and PSE-3D be possible; excellent agreement is shown in all aforementioned comparisons. Finally, the linear PSE-3D instability analysis is applied to a fully three-dimensional flow composed of a counter-rotating pair of nonparallel Batchelor vortices.

Global Linear Instability at the Dawn of its 4th Decade: A List of Challenges (A Practical Guide on how to Contain the Euphoria and Avoid the Oversell)

Global linear instability theory is concerned with the temporal or spatial development of small-amplitude perturbations superposed upon laminar steady or time-periodic threedimensional flows, which are inhomogeneous in two (and periodic in one) or all three spatial directions.1 The theory addresses flows developing in complex geometries, in which the parallel or weakly nonparallel basic flow approximation invoked by classic linear stability theory does not hold. As such, global linear theory is called to fill the gap in research into stability and transition in flows over or through complex geometries. Historically, global linear instability has been (and still is) concerned with solution of multi-dimensional eigenvalue problems; the maturing of non-modal linear instability ideas in simple parallel flows during the last decade of last century2–4 has given rise to investigation of transient growth scenarios in an ever increasing variety of complex flows. After a brief exposition of the theory, connections are sought with established approaches for structure identification in flows, such as the proper orthogonal decomposition and topology theory in the laminar regime and the open areas for future research, mainly concerning turbulent and three-dimensional flows, are highlighted. Recent results obtained in our group are reported in both the time-stepping and the matrix-forming approaches to global linear theory. In the first context, progress has been made in implementing a Jacobian-Free Newton Krylov method into a standard finite-volume aerodynamic code, such that global linear instability results may now be obtained in compressible flows of aeronautical interest. In the second context a new stable very high-order finite difference method is implemented for the spatial discretization of the operators describing the spatial BiGlobal EVP, PSE-3D and the TriGlobal EVP; combined with sparse matrix treatment, all these problems may now be solved on standard desktop computers.

Efficient methodologies for system matrix modelling in iterative image reconstruction for rotating high-resolution PET

A fully 3D iterative image reconstruction algorithm has been developed for high-resolution PET cameras composed of pixelated scintillator crystal arrays and rotating planar detectors, based on the ordered subsets approach. The associated system matrix is precalculated with Monte Carlo methods that incorporate physical effects not included in analytical models, such as positron range effects and interaction of the incident gammas with the scintillator material. Custom Monte Carlo methodologies have been developed and optimized for modelling of system matrices for fast iterative image reconstruction adapted to specific scanner geometries, without redundant calculations. According to the methodology proposed here, only one-eighth of the voxels within two central transaxial slices need to be modelled in detail. The rest of the system matrix elements can be obtained with the aid of axial symmetries and redundancies, as well as in-plane symmetries within transaxial slices. Sparse matrix techniques for the non-zero system matrix elements are employed, allowing for fast execution of the image reconstruction process. This 3D image reconstruction scheme has been compared in terms of image quality to a 2D fast implementation of the OSEM algorithm combined with Fourier rebinning approaches. This work confirms the superiority of fully 3D OSEM in terms of spatial resolution, contrast recovery and noise reduction as compared to conventional 2D approaches based on rebinning schemes. At the same time it demonstrates that fully 3D methodologies can be efficiently applied to the image reconstruction problem for high-resolution rotational PET cameras by applying accurate pre-calculated system models and taking advantage of the system's symmetries.