Preliminary results and documentation of sediment cores CRP2 and CRP-2A from the Ross Sea off Cape Roberts, Antarctica

The site for CRP-2, 14 km east of Cape Roberts (77.006°S; 163.719°E), was selected to overlap the early Miocene strata cored in nearby CRP-1, and to sample deeper into the east-dipping strata near the western margin ofe he Victoria Land Basin to investigate Palaeogene climatic and tectonic history. CRP-2 was cored from 5 to 57 mbsf (metres below the sea floor) (core recovery 91 %), with a deviation resulting in CRP-2A being cored at the same site. CRP-2A reached down to 624mbsf (recovery 95%), and to strata with an age of c. 33-35 Ma. Drilling took place from 16 October to 25 November 1998, on 2.0-2.2 m of sea ice and through 178 m of water. Core fractures and other physical properties, such as sonic velocity, density and magnetic susceptibility, were measured throughout the core. Down-hole logs for these and other properties were run from 63 to 167 mbsf and subsequently from 200 to 623 mbsf, although density and velocity data could be obtained only to 440 mbsf because of hole collapse. Sonic velocity averages c. 2.0 km S-1 for the upper part of the hole, but there is an sharp increase to c. 3.0 km s-1 and also a slight angular unconformity, at 306 mbsf, corresponding most likely to the early/late Oligocene boundary (c. 28-30 Ma). Velocity then increases irregularly to around 3.6 km s-1 at the bottom of the hole, which is estimated to lie 120 m above the V4/V5 boundary. The higher velocities below 306 mbsf probably reflect more extensive carbonate and common pyrite cementation, in patches, nodules, bedding-parallel masses and as vein infills. Dip of the strata also increases down-hole from 3° in the upper 300 in to over 10° at the bottom. Temperature gradient is 21° k-1. Over 2 000 fractures were logged through the hole. Borehole televiewer imagery was obtained for the interval from 200 to 440 mbsf to orient the fractures for stress field analysis. Lithostratigraphical descriptions on a scale of 1:20 are presented for the full length of the core, along with core box images, as a 200 page supplement to this issue. The hole initially passed through a layer of muddy gravel to 5.5 mbsf (Lithological Sub-Unit or LSU 1.1), and then into a Quaternary diatom-bearing clast-rich diamicton to 21 mbsf (LSU 2. l), with an interval of alternating compact diamicton and loose sand, and containing a rich Pliocene foraminiferal fauna, to 27 mbsf (LSU 2.2). The unit beneath this (LSU 3.1) has similar physical properties (sonic velocity, porosity, magnetic susceptibility) and includes diamictites of similar character to those of LSU 2.1 and 2.2, but an early Miocene (c. 19 Ma) diatom assemblage at 28 mbsf (top of LSU 3.1) shows that this sub-unit is part of the older section. The strata beneath 27 mbsf, primary target for the project, extend from early Miocene to perhaps latest Eocene age, and are largely cyclic glacimarine nearshore to offshore sediments. They are described as 41 lithological sub-units and interpreted in terms of 12 recurrent lithofacies. These are 1) mudstone, 2) inter-stratified mudstone and sandstone, 3) muddy very fine to coarse sandstone, 4) well-sorted stratified fine sandstone, 5) moderately to well-sorted, medium-grained sandstone, 6) stratified diamictite, 7) massive diamictite, 8) rhythmically inter-stratified sandstone and mudstone, 9) clast-supported conglomerate, 10) matrix-supported conglomerate, 11) mudstone breccia and 12) volcaniclastic sediment. Sequence stratigraphical analysis has identified 22 unconformity-bounded depositional sequences in pre- Pliocene strata. They typically comprise a four-part architecture involving, in ascending order, 1) a sharp-based coarse-grained unit (Facies 6,7,9 or 10), 2) a fining-upward succession of sandstones (Facies 3 and 4), 3) a mudstone interval (Facies l), in some cases coarsening upward to muddy sandstones (Facies 3), and 4) a sharp-based sandstone dominated succession (mainly Facies 4). The cyclicity recorded by the strata is interpreted in terms of a glacier ice margin retreating and advancing from land to the west, and of rises and falls in sea level. Analysis of sequence periodicity awaits afirmer chronology. However, apreliminary spectral analysis of magnetic susceptibility for a deepwater mudstone within one of the sequences (from 339 to 347 mbsf) reveals ratios between hierarchical levels that are similar to those of the three Milankovitch orbital forcing periodicities. The strata contain a wide range of fossils, the most abundant being marine diatoms. These commonly form up to 5% of the sediment, though in places the core is barren (notably between 300 and 412 mbsf). Fifty samples out of 250 reviewed were studied in detail. The assemblages define ten biostratigraphical zones, some of them based on local or as yet undescribed forms. The assemblages are neritic, and largely planktonic, suggesting that the sea floor was mostly below the photic zone throughout deposition of the corcd sequence. Calcareous nannofossils, representing incursions of ocean surface waters, are much less common (72 out of 183 samples examined) and restricted to mudstone intervals a few tens of metres thick, but are important for dating. Foraminifera are also sparse (73 out of 135 samples) and represented only by calcareous benthic species. Changing assemblages indicate a shift from inshore environments in the early Oligocenc to outer shelf in the late Oligocenc, returning to inshore in the early Miocene. Marine palynomorplis yielded large numbers of well-preserved forms from most of the 116 samples examined. The new in situ assemblagc found last year in CRP-1 is extended down into the late Oligocene and a further new assemblage is found in the early Oligoccnc. Many taxa are new, and cannot us yet contribute to an improved understanding of chronology or ecology. Marine invertebrate macrofossils, mostly molluscs and serpulid tubes, are scattered throughout the core. Preservation is good in mudstones but poor in other lithologies. Climate on land is reflected in the content of terrestrial palynomorphs, which are extremely scarce down to c. 300 mbsf. Some forms are reworked, and others represent a low growing sparse tundra with at least one species of Nothofagus. Beneath this level, a significantly greater diversity and abundance suggests a milder climate and a low diversity woody vegetation in the early Oligocene, but still far short of the richness found in known Eocene strata of the region. Sedimentary facies in the oldest strata also suggest a milder climate in the oldest strata cored, with indications of substantial glacial melt-water discharges, but are typical of a coldcr climate in late Oligocene and early Miocene times. Clast analyses from diamictites reveal weak to random fabrics, suggesting either lack of ice-contact deposition or post-depositional modification, but periods when ice grounded at the drill site are inferred from thin zones of in-situ brecciated rock and soft-sediment folding. These are more common above c. 300 mbsf, perhaps reflecting more extensive glacial advances during deposition of those strata. Erosion of the adjacent Transantarctic Mountains through Jurassic basalt and dolerite-intruded Beacon strata into basement rocks beneath is recorded by petrographical studies of clast and sand grain assemblages. Core below 310 mbsf contains a dominance of fine-grained Jurassic dolerite and basalt fragments along with Beacon-derived coal debris and rounded quartz grains, whereas the strata above this level have a much higher proportion of basement derived granitoids, implying that the large areas of the adjacent mountains had been eroded to basement by the end of the early Oligocene. There is little indication of rift-related volcanism below 310 mbsf. Above this, however, basaltic and trachytic tephras are common, especially from 280 to 200 mbsf, from 150 to 46 mbsf, and in Pliocene LSU 2.2 from 21 to 27 mbsf. The largest volcanic eruptions generated layers of coarse (up to 1 cm) trachytic pumice lapilli between 97 and 114 mbsf. The thickest of these (1.2 m at 112 mbsf) may have produced an eruptive column extending tens of km into the stratosphere. A source within a few tens of km of the drill site is considered most likely. Present age estimates for the pre-Pliocene sequence are based mainly on biostratigraphy (using mainly marine diatoms and to a lesser extent calcareous nannofossils), with the age of the tephra from 112 to 114 mbsf (21.44k0.05 Ma from 84 crystals by Ar-Ar) as a key reference point. Although there are varied and well-preserved microfossil assemblages through most of the sequence (notably of diatoms and marine palynomorphs), they comprise largely taxa either known only locally or as yet undescribed. In addition, sequence stratigraphical analysis and features in the core itself indicate numerous disconformities. The present estimate from diatom assemblages is that the interval from 27 to 130 mbsf is early Miocene in age (c. 19 to 23.5 Ma), consistent with the Ar-Ar age from 112 to 114 mbsf. Diatom assemblages also indicate that the late Oligocene epoch extends from c. 130 to 307 mbsf, which is supported by late Oligocene nannofossils from 130 to 185 mbsf. Strata from 307 to 412 mbsf have no age-diagnostic assemblages, but below this early Oligocene diatoms and nannofossils have been recovered. A nannoflora at the bottom of the hole is consistent with an earliest Oligocene or latest Eocene age. Magnetostratigraphical studies based on about 1000 samples, 700 of which have so far undergone demagnetisation treatment, have provided a polarity stratigraphy of 12 pre-Pliocene magnetozones. Samples above 270 mbsf are of consistently high quality. Below this, magnetic behaviour is more variable. A preliminary age-depth plot using the Magnetic Polarity Time Scale (MPTS) and constrained by biostratigraphical data suggests that episodes of relatively rapid sedimentation took place at CRP-2 during Oligocene times (c. 100 m/My), but that more than half of the record was lost in a few major and many minor disconformities. Age estimates from Sr isotopes in shell debris and further tephra dating are expected to lead to a better comparison with the MPTS. CRP-2/2A has recorded a history of subsidence of the Victoria Land Basin margin that is similar to that found in CIROS-170 km to the south, reflecting stability in both basin and the adjacent mountains in late Cenozoic times, but with slow net accumulation in the middle Cenozoic. The climatic indicators from both drill holes show a similar correspondence, indicating polar conditions for the Quaternary but with sub-polar conditions in the early Miocene-late Oligocene and indications of warmer conditions still in the early Oligocene. Correlation between the CRP-2A core and seismic records shows that seismic units V3 and V4, both widespread in the Victoria Land Basin, represent a period of fluctuating ice margins and glacimarine sedimentation. The next drill hole, CRP-3, is expected to core deep into V5 and extend this record of climate and tectonics still further back in time.

Paleomagnetic analysis of ODP Site 133-820

Paleomagnetic analysis of sediment samples from Ocean Drilling Program (ODP) Leg 133, Site 820, 10 km from the outer edge of the Great Barrier Reef, is undertaken to investigate the mineral magnetic response to environmental (sea level) changes. Viscous remanent magnetization (VRM) of both multidomain and near-superparamagnetic origin is prevalent and largely obscures the primary remanence, except in isolated high-magnetization zones. The Brunhes/Matuyama boundary cannot be identified, but is expected to be below 120 mbsf. The only evidence that exists for a geomagnetic excursion occurs at about 33 mbsf (-135 k.y.). Only one-half the cores were oriented, and many suffered from internal rotation about the core axis, caused by coring and/or slicing. The decay of magnetic remanence below the surface layer (0-2 mbsf) is attributed to sulfate reduction processes. The magnetic susceptibility (K) record is central for describing and understanding the magnetic properties of the sediments, and their relationship to glacio-eustatic fluctuations in sea level. Three prominent magnetic susceptibility peaks, at about 7, 32, and 64 mbsf, are superimposed on a background of smaller susceptibility oscillations. Fluctuations in susceptibility and remanence in the ôbackgroundö zone are controlled predominantly by variations in the concentration, rather than the composition of ferrimagnetics, with carbonate dilution playing an important role (type-A properties). The sharp susceptibility maxima occur at the start of the marine transgressions following low stands in sea level (high d18O, glacial maxima), and are characterized by a stable single-domain remanence, with a significant contribution from ultra-fine, superparamagnetic grains (type-C properties). During the later marine transgression, the susceptibility gradually returns to low values and the remanence is carried by stable single-domain magnetite (type-B properties). The A, B, and C types of sediment have distinctive ARM/K ratios. Throughout most of the sequence a strong inverse correlation exists between magnetic susceptibility and both CaCO3 and d18O variations. However, in the sharp susceptibility peaks (early transgression), more complex phase relationships are apparent among these parameters. In particular, the K-d18O correlation switches to positive, then reverts to negative during the course of the late transgression, indicating that two distinct mechanisms are responsible for the K-d18O correlation. Lower in the sequence, where sea-level-controlled cycles of upward-coarsening sediments, we find that the initial, mud phase of each cycle has been enriched in high-coercivity magnetic material, which is indicative of more oxic conditions. The main magnetic characteristics of the sediments are thought to reflect sea-level-controlled variations in the sediment source regions and related run-off conditions. Some preliminary evidence is seen that biogenic magnetite may play a significant role in the magnetization of these sediments.

Mineralogy and major element oxide composition of sediments from ODP Hole 172-1063D

Dansgaard-Oeschger (D-O) cycles in sediment at Site 1063 are characterized by distinct fluctuations in physical properties. Stadials are marked by low bulk density and interstadials by high bulk density. Compressional (P-)wave velocity is in phase with bulk density over some but not all depth intervals. Four of the D-O cycles straddling the oxygen isotope Stage 4/5 boundary have been studied in detail to understand the origin of the physical properties changes. Sediment on the Bermuda Rise is comprised of three main components: calcite, aluminosilicate minerals, and biogenic silica. Calcite concentrations vary from 1% to 43% of bulk sediment and are highest during interstadials. Aluminosilicate concentrations vary from 52% to 92% of bulk sediment and are highest during stadials. The major element ratios Al2O3/TiO2 and K2O/Al2O3 show increases across bulk density cycles, suggesting a change in the composition of aluminosilicates. This interpretation is supported by mineralogical analyses, which show a subtle change in clay composition. Biogenic silica concentrations vary from 0% to 23% of bulk sediment and are also highest during stadials. However, the abundance of silica varies significantly from one D-O cycle to another. Silt and fine sand abundance also increase during the first of the four stadials. This coarsening of sediment coincides with the increase in biogenic silica. The low grain density and high porosity associated with biogenic silica result in intervals of low bulk-sediment density. The abundance of biogenic silica closely matches P-wave velocity, suggesting that silica imparts a greater rigidity to the sediment.

Late Cenozoic eolian dust accumulation in the eastern equatorial Pacific Ocean

Sediments recovered during Ocean Drilling Program (ODP) Leg 138 in the eastern equatorial Pacific Ocean were analyzed for variations in eolian accumulation rate and mean grain-size. Latitudinal and temporal patterns of these parameters showed important changes in the intensity of atmospheric circulation and eolian flux associated with the intertropical convergence zone (ITCZ) and suggested that eolian input parameters could be used to define its paleoposition through time. Modern atmospheric circulation in the equatorial region is weakest in the intertropical convergence zone and increases as the trade winds are approached to the north and south. Thus, the expected spatial pattern of eolian grain size would have the finest material deposited beneath the ITCZ and a coarsening of material in both directions away from this zone. Sediments from ODP Leg 138 show this pattern for much of the Pleistocene and Pliocene but, prior to about 4 Ma, begin to lose the northern coarse component suggesting that the ITCZ was located north of its present position during the late Miocene. Eolian flux records also show a latitudinal pattern of deposition associated with the position of the ITCZ that, similar to eolian grain-size variability, suggests a more northerly position of the ITCZ during the late Miocene. Overall, the regional input of eolian material to the equatorial Pacific has decreased throughout the late Neogene. This reduction in eolian input reflects climatic changes to relatively wetter conditions in the continental eolian source regions beginning during the late Pliocene.

Sediment texture, carbonate content, and carbonate mineralogy of ODP Holes 182-1130A and 182-1132B

At Sites 1130 and 1132 of Ocean Drilling Program Leg 182 in the Great Australian Bight, we recovered an expanded Pleistocene section dominated by packstone and wackestone, deposited at unusually high rates of >20 cm/k.y. Shipboard observations detected an intermittent meter-scale alternation of light gray intervals with olive-gray intervals. Meter-scale samples were collected from the upper 250 m at both sites and decimeter-scale samples from four selected 2.5- to 4.0-m intervals in order to determine the texture and composition of sediments deposited along the upper slope throughout the Quaternary. Detailed textural and compositional data are presented from a total of 540 samples collected from both sites. Results indicate a general coarsening upward at both sites, with an accompanying upcore increase in high-Mg calcite (HMC) and aragonite and a decrease in low-Mg calcite (LMC). Samples collected at decimeter-scale intervals substantiate that the alternating light gray and olive-gray units detected on board ship are lithologically distinct. Light gray units consist of an LMC-rich silt, whereas olive-gray units consist of an aragonite and HMC-rich sand and silt. Sediment sources as well as timing and controls of this cyclic depositional pattern will be the subject of further investigations.

Silicic tephras in ODP Legs 151, 152, and 162

A 6-m.y.-long composite marine record of explosive silicic volcanism from five Ocean Drilling Program sites in the subpolar North Atlantic was compared with several marine records of global and local paleoclimate proxies (benthic d18O and ice-rafted debris records). Coarsening and high frequency of occurrence of Icelandic tephras were recorded in 3.6-3 Ma sediments, suggesting that these tephras were dispersed farther from the source by enhanced westerly winds over the subpolar North Atlantic. The 40Ar/39Ar ages were determined by laser probe on K-feldspar and biotite phenocrysts of tephras that were erupted from the Jan Mayen volcanic system. Compared to the tuned paleomagnetic age model, the 40Ar/39Ar dating (0.618+/-0.007 Ma to 4.90+/-0.05 Ma) yields a new age model that postdates by 155 k.y. the inception of ice rafting on the Iceland Plateau during the cold marine isotope stage M2 (i.e., 3.3-3.14 Ma).

(Table 1) Heavy minerals in silt and sand samples of IODP Site 301-U1301

The northwestern Cascadia Basin of western North America accumulated high-sedimentation-rate sequences during the Pleistocene sea-level low-stands. The continental shelf was largely exposed at that time, and rivers and estuaries delivered large sediment fluxes directly to the deep ocean. The IODP EXP1301 core, which was taken from the middle portion of the Cascadia Basin, is well preserved and exhibits the deeper and - more distal sedimentary facies. The lithology in this location is composed of two units, 1) hemipelagic mud with a thin sand layer and 2) thick, coarsening upward silt-sand turbidites with a small proportion of granules at the top. We will focus on the detailed sand-grain proportions in order to understand the origin of these sediments. We determined the modal proportions of the heavy minerals, and the chemical composition of olivine and orthopyroxene in fourteen samples. These are characterized by an abundance of amphibole, pyroxenes and epidote, and the presence of minerals derived from peridotite. There is no drastic change in the modal and mineral compositions of the sands and silts between the turbidite and hemipelagic sequences. There were two probable drainage systems on the continent, the Frazer and Columbia rivers, which shed turbidites into the Cascadia Basin after 1.6 Ma, especially at 0.46-0.76 Ma. Based on a comparison of the modal and mineral compositions, the Northern Cascadia Basin has been supplied with sediments, mainly from the Frazer River, through the Straits of Juan de Fuca, by Pleistocene to Holocene turbidites.

Magnetic properties of two cores from the Iberian Margin

Magnetic properties of late Quaternary sediments on the SW Iberian Margin are dominated by bacterial magnetite, observed by transmission electron microscopy (TEM), with contributions from detrital titanomagnetite and hematite. Reactive hematite from eolian dust, together with low organic matter concentrations and the lack of sulfate reduction, lead to dissimilatory iron reduction and availability of Fe(II) for abundant magnetotactic bacteria. Magnetite grain-size proxies (kARM/k and ARM/IRM) and S-ratios (sensitive to hematite) vary on stadial/interstadial timescales, contain orbital power, and mimic planktic d18O. The detrital/biogenic magnetite ratio and hematite concentration are greater during stadials and glacial isotopic stages, reflecting increased detrital (magnetite) input during times of lowered sea level, coinciding with atmospheric conditions favoring hematitic dust supply. Magnetic susceptibility, on the other hand, has a very different response being sensitive to coarse detrital multidomain (MD) magnetite associated with ice-rafted debris (IRD). High susceptibility and/or magnetic grain size coarsening, mark Heinrich stadials (HS), particularly HS2, HS3, HS4, HS5, HS6 and HS7, as well as older Heinrich-like detrital layers, indicating the sensitivity of this region to fluctuations in the position of the polar front. Relative paleointensity (RPI) records have well-constrained age models based on planktic d18O correlation to ice-core chronologies, however, they differ from reference records (e.g. PISO) particularly in the vicinity of glacial maxima, mainly due to inefficient normalization of RPI records in intervals of enhanced detrital/eolian hematite input.

Grain size and XRD analysis on Site 162-983

The grain size of deep-sea sediments provides an apparently simple proxy for current speed. However, grain size-based proxies may be ambiguous when the size distribution reflects a combination of processes, with current sorting only one of them. In particular, such sediment mixing hinders reconstruction of deep circulation changes associated with ice-rafting events in the glacial North Atlantic because variable ice-rafted detritus (IRD) input may falsely suggest current speed changes. Inverse modeling has been suggested as a way to overcome this problem. However, this approach requires high-precision size measurements that register small changes in the size distribution. Here we show that such data can be obtained using electrosensing and laser diffraction techniques, despite issues previously raised on the low precision of electrosensing methods and potential grain shape effects on laser diffraction. Down-core size patterns obtained from a sediment core from the North Atlantic are similar for both techniques, reinforcing the conclusion that both techniques yield comparable results. However, IRD input leads to a coarsening that spuriously suggests faster current speed. We show that this IRD influence can be accounted for using inverse modeling as long as wide size spectra are taken into account. This yields current speed variations that are in agreement with other proxies. Our experiments thus show that for current speed reconstruction, the choice of instrument is subordinate to a proper recognition of the various processes that determine the size distribution and that by using inverse modeling meaningful current speed reconstructions can be obtained from mixed sediments.

Bulk rock magnetic measurements and thermomagnetic analyses of sediment core GeoB4901-8 from the southeastern flank of the Niger deep-sea fan

Diagenesis has extensively affected the magnetic mineral inventory of organic-rich late Quaternary sediments in the Niger deep-sea fan. Changes in concentration, grain size, and coercivity document modifications of the primary magnetic mineral assemblages at two horizons. The first front, the modern iron redox boundary, is characterized by a drastic decline in magnetic mineral content, coarsening of the grain size spectrum, and reduction in coercivity. Beneath a second front, the transition from the suboxic to the sulfidic anoxic domain, a further but less pronounced decrease in concentration and bulk grain size occurs. Finer grains and higher coercive magnetic constituents substantially increase in the anoxic environment. Low- and high-temperature experiments were performed on bulk sediments and on extracts which have also been examined by X-ray diffraction. Thermomagnetic analyses proved ferrimagnetic titanomagnetites of terrigenous provenance as the principal primary magnetic mineral components. Their broad range of titanium contents reflects the volcanogenic traits of the Niger River drainage areas. Diagenetic alteration is not only a grain size selective process but also critically depends on titanomagnetite composition. Low-titanium compounds are less resistant to diagenetic dissolution. Intermediate titanium content titanomagnetite thus persists as the predominant magnetic mineral fraction in the sulfidic anoxic sediments. At the Fe redox boundary, precipitation of authigenic, possibly bacterial, magnetite is documented. The presence of hydrogen sulfide in the pore water suggests a formation of secondary magnetic iron sulfides in the anoxic domain. Grain size-specific data argue for a gradual development of a superparamagnetic and single-domain iron sulfide phase in this milieu, most likely greigite.

Topological measure locating the effective crossover between segregation and integration in a modular network

We introduce an easily computable topological measure which locates the effective crossover between segregation and integration in a modular network. Segregation corresponds to the degree of network modularity, while integration is expressed in terms of the algebraic connectivity of an associated hypergraph. The rigorous treatment of the simplified case of cliques of equal size that are gradually rewired until they become completely merged, allows us to show that this topological crossover can be made to coincide with a dynamical crossover from cluster to global synchronization of a system of coupled phase oscillators. The dynamical crossover is signaled by a peak in the product of the measures of intracluster and global synchronization, which we propose as a dynamical measure of complexity. This quantity is much easier to compute than the entropy (of the average frequencies of the oscillators), and displays a behavior which closely mimics that of the dynamical complexity index based on the latter. The proposed topological measure simultaneously provides information on the dynamical behavior, sheds light on the interplay between modularity and total integration, and shows how this affects the capability of the network to perform both local and distributed dynamical tasks.

Fast modularisation and aomic decomposition of ontologies using axiom dependency hypergraphs

In this paper we define the notion of an axiom dependency hypergraph, which explicitly represents how axioms are included into a module by the algorithm for computing locality-based modules. A locality-based module of an ontology corresponds to a set of connected nodes in the hypergraph, and atoms of an ontology to strongly connected components. Collapsing the strongly connected components into single nodes yields a condensed hypergraph that comprises a representation of the atomic decomposition of the ontology. To speed up the condensation of the hypergraph, we first reduce its size by collapsing the strongly connected components of its graph fragment employing a linear time graph algorithm. This approach helps to significantly reduce the time needed for computing the atomic decomposition of an ontology. We provide an experimental evaluation for computing the atomic decomposition of large biomedical ontologies. We also demonstrate a significant improvement in the time needed to extract locality-based modules from an axiom dependency hypergraph and its condensed version.

Strong Connectivity in Directed Hypergraphs and its Application to the Atomic Decomposition of Ontologies

Los hipergrafos dirigidos se han empleado en problemas relacionados con lógica proposicional, bases de datos relacionales, linguística computacional y aprendizaje automático. Los hipergrafos dirigidos han sido también utilizados como alternativa a los grafos (bipartitos) dirigidos para facilitar el estudio de las interacciones entre componentes de sistemas complejos que no pueden ser fácilmente modelados usando exclusivamente relaciones binarias. En este contexto, este tipo de representación es conocida como hiper-redes. Un hipergrafo dirigido es una generalización de un grafo dirigido especialmente adecuado para la representación de relaciones de muchos a muchos. Mientras que una arista en un grafo dirigido define una relación entre dos de sus nodos, una hiperarista en un hipergrafo dirigido define una relación entre dos conjuntos de sus nodos. La conexión fuerte es una relación de equivalencia que divide el conjunto de nodos de un hipergrafo dirigido en particiones y cada partición define una clase de equivalencia conocida como componente fuertemente conexo. El estudio de los componentes fuertemente conexos de un hipergrafo dirigido puede ayudar a conseguir una mejor comprensión de la estructura de este tipo de hipergrafos cuando su tamaño es considerable. En el caso de grafo dirigidos, existen algoritmos muy eficientes para el cálculo de los componentes fuertemente conexos en grafos de gran tamaño. Gracias a estos algoritmos, se ha podido averiguar que la estructura de la WWW tiene forma de “pajarita”, donde más del 70% del los nodos están distribuidos en tres grandes conjuntos y uno de ellos es un componente fuertemente conexo. Este tipo de estructura ha sido también observada en redes complejas en otras áreas como la biología. Estudios de naturaleza similar no han podido ser realizados en hipergrafos dirigidos porque no existe algoritmos capaces de calcular los componentes fuertemente conexos de este tipo de hipergrafos. En esta tesis doctoral, hemos investigado como calcular los componentes fuertemente conexos de un hipergrafo dirigido. En concreto, hemos desarrollado dos algoritmos para este problema y hemos determinado que son correctos y cuál es su complejidad computacional. Ambos algoritmos han sido evaluados empíricamente para comparar sus tiempos de ejecución. Para la evaluación, hemos producido una selección de hipergrafos dirigidos generados de forma aleatoria inspirados en modelos muy conocidos de grafos aleatorios como Erdos-Renyi, Newman-Watts-Strogatz and Barabasi-Albert. Varias optimizaciones para ambos algoritmos han sido implementadas y analizadas en la tesis. En concreto, colapsar los componentes fuertemente conexos del grafo dirigido que se puede construir eliminando ciertas hiperaristas complejas del hipergrafo dirigido original, mejora notablemente los tiempos de ejecucion de los algoritmos para varios de los hipergrafos utilizados en la evaluación. Aparte de los ejemplos de aplicación mencionados anteriormente, los hipergrafos dirigidos han sido también empleados en el área de representación de conocimiento. En concreto, este tipo de hipergrafos se han usado para el cálculo de módulos de ontologías. Una ontología puede ser definida como un conjunto de axiomas que especifican formalmente un conjunto de símbolos y sus relaciones, mientras que un modulo puede ser entendido como un subconjunto de axiomas de la ontología que recoge todo el conocimiento que almacena la ontología sobre un conjunto especifico de símbolos y sus relaciones. En la tesis nos hemos centrado solamente en módulos que han sido calculados usando la técnica de localidad sintáctica. Debido a que las ontologías pueden ser muy grandes, el cálculo de módulos puede facilitar las tareas de re-utilización y mantenimiento de dichas ontologías. Sin embargo, analizar todos los posibles módulos de una ontología es, en general, muy costoso porque el numero de módulos crece de forma exponencial con respecto al número de símbolos y de axiomas de la ontología. Afortunadamente, los axiomas de una ontología pueden ser divididos en particiones conocidas como átomos. Cada átomo representa un conjunto máximo de axiomas que siempre aparecen juntos en un modulo. La decomposición atómica de una ontología es definida como un grafo dirigido de tal forma que cada nodo del grafo corresponde con un átomo y cada arista define una dependencia entre una pareja de átomos. En esta tesis introducimos el concepto de“axiom dependency hypergraph” que generaliza el concepto de descomposición atómica de una ontología. Un modulo en una ontología correspondería con un componente conexo en este tipo de hipergrafos y un átomo de una ontología con un componente fuertemente conexo. Hemos adaptado la implementación de nuestros algoritmos para que funcionen también con axiom dependency hypergraphs y poder de esa forma calcular los átomos de una ontología. Para demostrar la viabilidad de esta idea, hemos incorporado nuestros algoritmos en una aplicación que hemos desarrollado para la extracción de módulos y la descomposición atómica de ontologías. A la aplicación la hemos llamado HyS y hemos estudiado sus tiempos de ejecución usando una selección de ontologías muy conocidas del área biomédica, la mayoría disponibles en el portal de Internet NCBO. Los resultados de la evaluación muestran que los tiempos de ejecución de HyS son mucho mejores que las aplicaciones más rápidas conocidas. ABSTRACT Directed hypergraphs are an intuitive modelling formalism that have been used in problems related to propositional logic, relational databases, computational linguistic and machine learning. Directed hypergraphs are also presented as an alternative to directed (bipartite) graphs to facilitate the study of the interactions between components of complex systems that cannot naturally be modelled as binary relations. In this context, they are known as hyper-networks. A directed hypergraph is a generalization of a directed graph suitable for representing many-to-many relationships. While an edge in a directed graph defines a relation between two nodes of the graph, a hyperedge in a directed hypergraph defines a relation between two sets of nodes. Strong-connectivity is an equivalence relation that induces a partition of the set of nodes of a directed hypergraph into strongly-connected components. These components can be collapsed into single nodes. As result, the size of the original hypergraph can significantly be reduced if the strongly-connected components have many nodes. This approach might contribute to better understand how the nodes of a hypergraph are connected, in particular when the hypergraphs are large. In the case of directed graphs, there are efficient algorithms that can be used to compute the strongly-connected components of large graphs. For instance, it has been shown that the macroscopic structure of the World Wide Web can be represented as a “bow-tie” diagram where more than 70% of the nodes are distributed into three large sets and one of these sets is a large strongly-connected component. This particular structure has been also observed in complex networks in other fields such as, e.g., biology. Similar studies cannot be conducted in a directed hypergraph because there does not exist any algorithm for computing the strongly-connected components of the hypergraph. In this thesis, we investigate ways to compute the strongly-connected components of directed hypergraphs. We present two new algorithms and we show their correctness and computational complexity. One of these algorithms is inspired by Tarjan’s algorithm for directed graphs. The second algorithm follows a simple approach to compute the stronglyconnected components. This approach is based on the fact that two nodes of a graph that are strongly-connected can also reach the same nodes. In other words, the connected component of each node is the same. Both algorithms are empirically evaluated to compare their performances. To this end, we have produced a selection of random directed hypergraphs inspired by existent and well-known random graphs models like Erd˝os-Renyi and Newman-Watts-Strogatz. Besides the application examples that we mentioned earlier, directed hypergraphs have also been employed in the field of knowledge representation. In particular, they have been used to compute the modules of an ontology. An ontology is defined as a collection of axioms that provides a formal specification of a set of terms and their relationships; and a module is a subset of an ontology that completely captures the meaning of certain terms as defined in the ontology. In particular, we focus on the modules computed using the notion of syntactic locality. As ontologies can be very large, the computation of modules facilitates the reuse and maintenance of these ontologies. Analysing all modules of an ontology, however, is in general not feasible as the number of modules grows exponentially in the number of terms and axioms of the ontology. Nevertheless, the modules can succinctly be represented using the Atomic Decomposition of an ontology. Using this representation, an ontology can be partitioned into atoms, which are maximal sets of axioms that co-occur in every module. The Atomic Decomposition is then defined as a directed graph such that each node correspond to an atom and each edge represents a dependency relation between two atoms. In this thesis, we introduce the notion of an axiom dependency hypergraph which is a generalization of the atomic decomposition of an ontology. A module in the ontology corresponds to a connected component in the hypergraph, and the atoms of the ontology to the strongly-connected components. We apply our algorithms for directed hypergraphs to axiom dependency hypergraphs and in this manner, we compute the atoms of an ontology. To demonstrate the viability of this approach, we have implemented the algorithms in the application HyS which computes the modules of ontologies and calculate their atomic decomposition. In the thesis, we provide an experimental evaluation of HyS with a selection of large and prominent biomedical ontologies, most of which are available in the NCBO Bioportal. HyS outperforms state-of-the-art implementations in the tasks of extracting modules and computing the atomic decomposition of these ontologies.

Feldspar size measurements of tephra layers from ODP Leg 165 sites

Crystal size measurements have been carried out on tephra fall layers of Miocene to recent age from Sites 998, 999, and 1000 in the western Caribbean Sea. Maximum crystal size is used as a proxy for the grain size characteristics of the layers and an index of atmospheric dispersal from source eruptions. Crystal sizes range from 50 to 650 µm with the majority falling between 200 and 300 µm. All three sites exhibit a coarsening in the grain size of tephra layers with increasing age to the early Miocene that broadly correlates with an increase in the frequency of layers. Analysis of the present lower and upper level atmospheric circulation in the western Caribbean suggests that the layers were derived from source eruptions to the west of the sites somewhere in the Central American region. Minimum distances to these sources are of the order of 700 km. Crystal sizes in tephra layers at these distances are consistent with their derivation from energetic pyroclastic flow-forming eruptions that injected tephra to stratospheric levels by large-scale co-ignimbrite and plinian-style plumes. Coarsening of the layers during the Miocene peak of explosive volcanism cannot be attributed to any major change in paleowind intensity and is taken to represent the occurrence of more energetic eruptions that were able to disperse tephra over larger areas.

(Appendix) Grain-size distribution, calcium carbonate content and proportion of cemented aggregates in ODP Site 133-820

The textural and compositional characteristics of the 400 m sequence of Pleistocene wackestones and packstones intersected at Ocean Drilling Program (ODP) Site 820 reflect deposition controlled by fluctuations in sea-level, and by variations in the rate of sediment supply. The development of an effective reefal barrier adjacent to Site 820, between 760 k.y. and 1.01 Ma, resulted in a marked reduction in sediment accumulation rates on the central Great Barrier Reef outermost shelf and upper slope. This marked change corresponds with the transition from sigmoidal prograding seismic geometry in the lower 254 m of the sequence, to aggradational geometry in the top 146 m. The reduction in the rate of sediment accumulation that followed development of the reefal barrier also caused a fundamental change in the way in which fluctuations in sea-level controlled sediment deposition. In the lower, progradational portion of the sequence, sea-level cyclicity is represented by superimposed coarsening-upward cycles. Although moderately calcareous throughout (mostly 35%-75% CaCO3), the depositional system acted in a similar manner to siliciclastic shelf depositional systems. Relative sea-level rises resulted in deposition of more condensed, less calcareous, fine, muddy wackestones at the base of each cycle. Sea-level highstands resulted in increased sedimentation rates and greater influx of coarse bioclastic material. Continued high rates of sedimentation of both coarse bioclastic material and mixed carbonate and terrigenous mud marked falling and low sea-levels. This lower part of the sequence therefore is dominated by coarse packstones, with only thin wackestone intervals representing transgressions. In contrast, sea-level fluctuations following formation of an effective reefal barrier produced a markedly different sedimentary record. The more slowly deposited aggradational sequence is characterized by discrete thin interbeds of relatively coarse packstone within a predominantly fine wackestone sequence. These thin packstone beds resulted from relatively low sedimentation rates during falling and low sea-levels, with much higher rates of muddy sediment accumulation during rising and high sea-levels. The transition from progradational to aggradational sequence geometry therefore corresponds to a transition from a "siliciclastic-type" to a "carbonate-type" depositional system.