Magnetic properties, grain sizes and minerals at DSDP Site 83-504B

Twenty-seven samples from the Leg 83 section of Hole 504B have been investigated using magnetic, optical, and electron optical methods. The primary magnetic mineral to crystallize was titanomagnetite of approximate composition Fe2.4Ti0.6O4 (TM60), but none survives, nor is there evidence of titanomaghemite produced by oxidation of TM60. The average measured magnetic properties can be interpreted in terms of magnetite, Fe3O4, having average grain size of <1 µm and present in average volume concentration of - 0.5%. The intensity of the natural remanent magnetization (NRM) of the rocks could also be accounted for as being a thermoremanence carried by this mineral. Although the heterogeneity of the titanomagnetite grains could be detected optically, the texture of the intergrown phases is poorly developed. In some samples from the massive units of the lower part of the section, trellis patterns were visible. The Fe3O4 present in the intergrowths is too intimately mixed with the other intergrown phases to be revealed by electron microprobe analysis that simply returns the bulk composition of the intergrowth (oxidized TM60). The path by which the mineral assemblage evolved from TM60 to an Fe304-containing intergrowth, under the temperature and pressure conditions obtaining in the Leg 83 section, makes interesting speculation. Deuteric oxidation, maghemitization/inversion, or some hypothetical low-temperature/high-pressure oxidation by a leaching-of-iron process may all play roles.

Stratigraphic correlation of ODP Site 199-1218 and 199-1219

During Ocean Drilling Program Leg 199 a high-resolution (~1-2 cm/k.y.) biogenic sediment record from the late Paleocene to the early Miocene was recovered, containing an uninterrupted set of geomagnetic chrons as well as a detailed record of calcareous and siliceous biostratigraphic datum events. Shipboard lithologic proxy measurements and shore-based determinations of CaCO3 revealed regular cycles that can be attributed to climatic forcing. Discovering drill sites with well defined magneto- and biostratigraphic records that also show clear lithologic cycles is rare and valuable and creates the opportunity to develop a detailed stratigraphic intersite correlation, providing the basis to study paleoceanographic processes and mass accumulation rates at high resolution. Here we present extensive postcruise work that extends the shipboard composite depth stratigraphy by providing a high-resolution revised meters composite depth (rmcd) scale to compensate for depth distortion within individual cores. The depth-aligned data were then used to generate stacked records of lithologic proxy measurements. Making use of the increased signal-to-noise ratio in the stacked records, we then proceeded to generate a detailed site-to-site correlation between Sites 1218 and 1219 in order to decrease the depth uncertainty for magneto- and biostratigraphic datums. Stacked lithologic proxy records in combination with discrete measurements of CaCO3 were then exploited to calculate high-resolution carbonate concentration curves by regression of the multisensor track data with discrete measurements. By matching correlative features between the cores and wireline logging data, we also rescaled our core rmcd back to in situ depths. Our study identifies lithology-dependent core expansion due to unloading as the mechanism of varying stratigraphic thicknesses between cores.

Radiolarian chronology and fauna across the middle/late Eocene boundary at ODP Hole 171-1052A

Quantitative radiolarian assemblage analysis has been conducted on middle and upper Eocene sediments (Zones RP16 to RP18) from Ocean Drilling Program Site 1052 in order to establish the radiolarian magnetobiochronology and determine the nature of the faunal turnover across the middle/late Eocene boundary in the western North Atlantic Ocean. We recognize and calibrate forty-five radiolarian bioevents to the magneto- and cyclo-stratigraphy from Site 1052 to enhance the biochronologic resolution for the middle and late Eocene. Our data is compared to sites in the equatorial Pacific (Leg 199) to access the diachrony of biostratigraphic events. Eleven bioevents are good biostratigraphic markers for tropical/subtropical locations (south of 30°N). The primary markers (lowest occurrences of Cryptocarpium azyx and Calocyclas bandyca) which are tropical zonal boundary markers for Zones RP17 and RP18 provide robust biohorizons for correlation and age determination from the low to middle latitudes and between the Atlantic and Pacific Oceans. Some other radiolarian bioevents are highly diachronous (<1 million years) between oceanic basins. A significant faunal turnover of radiolarians is recognized within Chron C17n.3n (37.7 Ma) where 13 radiolarian species disappear rapidly in less than 100 kyr and 4 new species originate. The radiolarian faunal turnover coincides with a major extinction in planktonic foraminifera. We name the turnover phase, the Middle/Late Eocene Turnover (MLET). Assemblage analysis reveals the MLET to be associated with a decrease in low-mid latitude taxa and increase in cosmopolitan taxa and radiolarian accumulation rates. The MLET might be related to increased biological productivity rather than to surface-water cooling.

Middle Eocene to Pleistocene calcareous nannofossils of ODP Leg 113 holes

Long sequences of Upper Cretaceous through Quaternary sediments rich in calcareous and siliceous microfossils were recovered at Ocean Drilling Program Sites 689 and 690 on Maud Rise off East Antarctica. These sites have become the southernmost anchor in the Atlantic Basin for bio-, magneto-, chemostratigraphic, and paleobiogeographic studies. ODP Sites 692 and 693 on the Weddell Sea margin of East Antarctica and Site 696 on the South Orkney microcontinent of West Antarctica yielded calcareous nannofossils within some stratigraphic intervals. Sites 691, 692, 694, 695, and 697 did not recover Cenozoic calcareous nannofossils. Calcareous nannofossil biostratigraphy suggests a major hiatus across the Paleogene/Neogene boundary at Sites 689 and 690, and two additional hiatuses in the middle Eocene-lower Oligocene section at Site 690. Correlation with magnetostratigraphy reveals: the last occurrence (LO) of Reticulofenestra umbilica at Maud Rise is over 1 m.y. younger than that at the middle-latitude sites; the LO of Isthmolithus recurvus is synchronous in the middle-latitude and high-latitude areas (about 34.8 Ma); Reticulofenestra oamaruensis ranges from 38.0 to 36.0 Ma at Maud Rise; Reticulofenestra reticulata has a shorter range at Maud Rise (42.1 to 38.9 Ma) than at the middle-latitude DSDP Site 516; the range of Chiasmolithus oamaruensis is diachronous over different latitudes; and the LO of Chiasmolithus solitus is a good datum at 41.3 Ma from 30°S to 65°S in the South Atlantic Ocean. Comparison of calcareous nannofossil abundances in a latitudinal transect shows: Reticulofenestra bisecta is a temperate-water species and its LO, which crosses below that of Chiasmolithus altus at Maud Rise, is not applicable for the Paleogene/Neogene boundary in high southern latitude areas; Clausicoccus fenestratus is rare or absent at Maud Rise and can not be used as a marker; Coccolithus formosus is a warm-water species which disappeared earlier toward higher latitudes. Calcareous nannofossil assemblages indicate that by at least the middle Eocene, surface water temperatures became considerably lower in the high southern latitudes than in the middle-latitude areas and that there have been more extreme cold events in the high latitudes during the Neogene. Bicolumnus ovatus n. gen., n. sp. is proposed in this paper.

X-ray fluorescence data (calcium and iron), bulk organic d13C, and dinoflagellate cysts of ODP Hole 189-1172D

The 'Paleocene/Eocene Thermal Maximum' or PETM (~55 Ma) was associated with dramatic warming of the oceans and atmosphere, pronounced changes in ocean circulation and chemistry, and upheaval of the global carbon cycle. Many relatively complete PETM sequences have by now been reported from around the world, but most are from ancient low- to midlatitude sites. ODP Leg 189 in the Tasman Sea recovered sediments from this critical phase in Earth history at Sites 1171 and 1172, potentially representing the southernmost PETM successions ever encountered (at ~70° to 65° S paleolatitude). Downhole and core logging data, in combination with dinoflagellate cyst biostratigraphy, magneto-stratigraphy, and stable isotope geochemistry indicate that the sequences at both sites were deposited in a high accumulation-rate, organic rich, marginal marine setting. Furthermore, Site 1172 indeed contains a fairly complete P-E transition, whereas at Site 1171, only the lowermost Eocene is recovered. However, at Site 1172, the typical PETM-indicative acme of the dinocyst Apectodinium was not recorded. We conclude that unfortunately, the critical latest Paleocene and PETM intervals are missing at Site 1172. We relate the missing section to a sea level driven hiatus and/or condensed section and recovery problems. Nevertheless, our integrated records provide a first-ever portrait of the trend toward, and aftermath of, the PETM in a marginal marine, southern high-latitude setting.

(Table 1) Viscosity coefficients at various temperatures of samples from ODP Hole 176-735B

Ocean Drilling Program (ODP) Hole 735B was drilled to a depth of 1.5 km in a tectonic window of gabbroic lower oceanic crust created at the Southwest Indian Ridge. The gabbros have a very stable natural remanent magnetization (NRM) of reversed polarity with most unblocking temperatures slightly below the Curie temperature of magnetite. The NRM includes a drilling-induced overprint but its intensity decays strongly towards the interior of the drill core. The demagnetization data yield no or only a very small secondary magnetization component acquired during the present Brunhes chron or an earlier normal chron, suggesting cooling through most of the blocking temperature range during chron C5r and a strong resistance against the acquisition of thermoviscous magnetization. A novel furnace has been designed to measure magnetizations and their time dependences at high temperatures (up to 580 deg C) inside a commercial SQUID magnetometer. Magnetic viscosity experiments have been conducted on the gabbros at temperatures up to 550 deg C to determine the time and temperature stability of remanent magnetization. Viscosities are generally small and increase little with temperature below the main blocking temperature, where the increase becomes almost an order of magnitude. Extrapolations to geological times infer viscous acquisitions that would be 5-25% of a thermoremanence in 100 kyr and at temperatures of 200-500 deg C. At ocean bottom temperature the predicted magnetization of one sample acquired in the present Brunhes chron should be 10% of the NRM. However, this is not recognized during NRM demagnetization and partial thermoremanent magnetization (pTRM) acquisitions at 250 deg C are also much smaller than predicted. It thus appears that the NRMs are generally magnetically harder than magnetizations acquired after heating to 570 deg C in the laboratory. Susceptibility changes during heating are small (<5%) indicating a seemingly stable magneto-mineralogy, but conspicuous minima occur after heating to 520 deg C. Also, quasi paleointensity experiments reveal characteristic patterns in the NRM/pTRM ratios and also large increases in pTRM capacity after heating to 570 deg C. Moreover, anhysteretic remanent magnetization acquisition in the low field range (<=10 mT) is strongly enhanced after heating by factors up to three. The alteration of the magneto-mineralogy is interpreted to result from the annealing of defects in magnetite that originate from tectonically induced strain. The oceanic gabbros of Hole 735B are thus ideal source layer material for marine magnetic anomalies, and secondary thermoviscous acquisition, as a possible cause for anomalous skewness, is essentially absent.

Cretaceous calcareous nannofossils from Exmouth and Wombat Plateau, Northwest Australia

Three sites drilled during Leg 122, Site 761 on the Wombat Plateau and Sites 762 and 763 on the Exmouth Plateau, provide a composite Cretaceous section ranging in age from Berriasian to Maestrichtian. Together, these sites contain an apparently complete, expanded Aptian-Maestrichtian record. Consistently occurring and moderately well-preserved nannofossil assemblages allow reasonably high biostratigraphic resolution. Our data indicate that traditional middle and Upper Cretaceous nannofossil biozonations are not entirely applicable in this region. In this investigation, we compare in detail the relative ranges of key Cretaceous nannofossil markers in the eastern Indian Ocean and in sections from Europe and North Africa. We have determined which previously used events are applicable, and which additional markers have biostratigraphic utility in this region. Significant differences in Campanian-Maestrichtian assemblages exist between the more northern Site 761 and Sites 762 and 763. Such differences are surprising, considering that these sites are only separated by 3° of latitude. We interpret them as marking a strong thermal gradient over the Exmouth Plateau region. Other results include the recovery of an expanded Albian-Cenomanian sequence containing a mixture of Austral and Tethyan floras, which will enable correlation of biozonations established for these two realms; the recovery of two condensed but apparently complete Cenomanian-Turonian boundary sections; correlation of Upper Cretaceous calcareous nannofossil biostratigraphy with magneto- and foraminifer stratigraphy; and correlation of portions of the Barrow Group equivalents to the Berriasian and Valanginian stages.

Paleomagnetic and rock magnetic data from IODP Site U1308

Integrated Ocean Drilling Program (IODP) Site U1308 (central North Atlantic) records paleomagnetic directional and relative paleointensity (RPI) variations for the last 1.5 Myr, in 110 m of the sediment sequence at a mean sedimentation rate of 7.3 cm/kyr. A detailed benthic oxygen isotope record was combined with RPI to produce an integrated, high-resolution magneto-isotopic stratigraphy for Site U1308. Apart from the well-known polarity reversals in this interval, the Punaruu excursion is recorded at 1092 ka and the Cobb Mountain Subchron in the 1182-1208 ka interval. The paleointensity proxies are determined as slopes of NRM versus ARM and NRM versus ARMAQ (ARM acquisition) with linear correlation coefficients to monitor the quality of the linear fit. The RPI record for Site U1308 is compared with the three other paleointensity records (one from the Western Equatorial Pacific and two from the North Atlantic) that cover the same time interval and have accompanying oxygen isotope records. The Match protocol of Lisiecki and Lisiecki (2002) is used to optimize the correlation of paleointensity records. Beginning with the original (published) age models for each record, the Match routine is used to optimize the RPI correlations to Site U1308, with checks to ensure compatibility with oxygen isotope records. Squared wavelet coherence (WTC) indicates significant improvement in RPI (and oxygen isotope) correlations after matching each RPI record to Site U1308, particularly for periods > 10 kyr. The level of coherence for the Atlantic RPI records and the lower resolution Pacific record implies synchronous global variability (at scales > 10 kyr) that can be attributed to the axial dipole geomagnetic field.

Reduced order models for industrial thermo-fluid problems

A new method is presented to generate reduced order models (ROMs) in Fluid Dynamics problems of industrial interest. The method is based on the expansion of the flow variables in a Proper Orthogonal Decomposition (POD) basis, calculated from a limited number of snapshots, which are obtained via Computational Fluid Dynamics (CFD). Then, the POD-mode amplitudes are calculated as minimizers of a properly defined overall residual of the equations and boundary conditions. The method includes various ingredients that are new in this field. The residual can be calculated using only a limited number of points in the flow field, which can be scattered either all over the whole computational domain or over a smaller projection window. The resulting ROM is both computationally efficient(reconstructed flow fields require, in cases that do not present shock waves, less than 1 % of the time needed to compute a full CFD solution) and flexible(the projection window can avoid regions of large localized CFD errors).Also, for problems related with aerodynamics, POD modes are obtained from a set of snapshots calculated by a CFD method based on the compressible Navier Stokes equations and a turbulence model (which further more includes some unphysical stabilizing terms that are included for purely numerical reasons), but projection onto the POD manifold is made using the inviscid Euler equations, which makes the method independent of the CFD scheme. In addition, shock waves are treated specifically in the POD description, to avoid the need of using a too large number of snapshots. Various definitions of the residual are also discussed, along with the number and distribution of snapshots, the number of retained modes, and the effect of CFD errors. The method is checked and discussed on several test problems that describe (i) heat transfer in the recirculation region downstream of a backwards facing step, (ii) the flow past a two-dimensional airfoil in both the subsonic and transonic regimes, and (iii) the flow past a three-dimensional horizontal tail plane. The method is both efficient and numerically robust in the sense that the computational effort is quite small compared to CFD and results are both reasonably accurate and largely insensitive to the definition of the residual, to CFD errors, and to the CFD method itself, which may contain artificial stabilizing terms. Thus, the method is amenable for practical engineering applications. Resumen Se presenta un nuevo método para generar modelos de orden reducido (ROMs) aplicado a problemas fluidodinámicos de interés industrial. El nuevo método se basa en la expansión de las variables fluidas en una base POD, calculada a partir de un cierto número de snapshots, los cuales se han obtenido gracias a simulaciones numéricas (CFD). A continuación, las amplitudes de los modos POD se calculan minimizando un residual global adecuadamente definido que combina las ecuaciones y las condiciones de contorno. El método incluye varios ingredientes que son nuevos en este campo de estudio. El residual puede calcularse utilizando únicamente un número limitado de puntos del campo fluido. Estos puntos puede encontrarse dispersos a lo largo del dominio computacional completo o sobre una ventana de proyección. El modelo ROM obtenido es tanto computacionalmente eficiente (en aquellos casos que no presentan ondas de choque reconstruir los campos fluidos requiere menos del 1% del tiempo necesario para calcular una solución CFD) como flexible (la ventana de proyección puede escogerse de forma que evite contener regiones con errores en la solución CFD localizados y grandes). Además, en problemas aerodinámicos, los modos POD se obtienen de un conjunto de snapshots calculados utilizando un código CFD basado en la versión compresible de las ecuaciones de Navier Stokes y un modelo de turbulencia (el cual puede incluir algunos términos estabilizadores sin sentido físico que se añaden por razones puramente numéricas), aunque la proyección en la variedad POD se hace utilizando las ecuaciones de Euler, lo que hace al método independiente del esquema utilizado en el código CFD. Además, las ondas de choque se tratan específicamente en la descripción POD para evitar la necesidad de utilizar un número demasiado grande de snapshots. Varias definiciones del residual se discuten, así como el número y distribución de los snapshots,el número de modos retenidos y el efecto de los errores debidos al CFD. El método se comprueba y discute para varios problemas de evaluación que describen (i) la transferencia de calor en la región de recirculación aguas abajo de un escalón, (ii) el flujo alrededor de un perfil bidimensional en regímenes subsónico y transónico y (iii) el flujo alrededor de un estabilizador horizontal tridimensional. El método es tanto eficiente como numéricamente robusto en el sentido de que el esfuerzo computacional es muy pequeño comparado con el requerido por el CFD y los resultados son razonablemente precisos y muy insensibles a la definición del residual, los errores debidos al CFD y al método CFD en sí mismo, el cual puede contener términos estabilizadores artificiales. Por lo tanto, el método puede utilizarse en aplicaciones prácticas de ingeniería.