Mineralogy and stable oxygen isotope ratios of the Cretaceous-Tertiary boundary

Results of detailed mineralogical, chemical, and oxygen isotope analyses of the clay minerals and zeolites from two Cretaceous-Tertiary (K/T) boundary regions, Stevns Klint, Denmark, and Deep Sea Drilling Project (DSDP) Hole 465A in the north central Pacific Ocean, are presented. In the central part of the Stevns Klint K/T boundary layer, the only clay mineral detected by x-ray diffraction is a pure smectite with > 95 percent expandable layers. No detrital clay minerals or quartz were observed in the clay size fraction in these beds, whereas the clay minerals above and below the boundary layer are illite and mixed-layer smectite-illite of detrital origin as well as quartz. The mineralogical purity of the clay fraction, the presence of smectite only at the boundary, and the d18O value of the smectite (27.2 ± 0.2 per mil) suggest that it formed in situ by alteration of glass. Formation from impact rather than from volcanic glass is supported by its major element chemistry. The high content of iridium and other siderophile elements is not due to the cessation of calcium carbonate deposition and resulting slow sedimentation rates. At DSDP Hole 465A, the principal clay mineral in the boundary zone (80 to 143 centimeters) is a mixed-layer smectite-illite with >=90 percent expandable layers, accompanied by some detrital quartz and small amounts of a euhedral authigenic zeolite (clinoptilolite). The mixed-layer smectite-illite from the interval 118 to 120 centimeters in the zone of high iridium abundance has a very low rare earth element content; the negative cerium anomaly indicates formation in the marine environment. This conclusion is corroborated by the d18O value of this clay mineral (27.1 ± 0.2 per mil). Thus, this mixed-layer smectite-illite formed possibly from the same glass as the K/T boundary smectite at Stevns Klint, Denmark.

Marine osmium isotope record across the Cretaceous-Tertiary boundary

A composite late Maastrichtian (65.5 to 68.5 Ma) marine osmium (Os) isotope record, based on samples from the Southern Ocean (ODP Site 690), the Tropical Pacific Ocean (DSDP Site 577), the South Atlantic (DSDP Site 525) and the paleo-Tethys Ocean demonstrates that subaerially exposed pelagic carbonates can record seawater Os isotope variations with a fidelity comparable to sediments recovered from the seafloor. New results provide robust evidence of a 20% decline in seawater 187Os/188Os over a period of about 200 kyr early in magnetochron C29r well below the Cretaceous-Paleogene Boundary (KPB), confirming previously reported low-resolution data from the South Atlantic Ocean. New results also confirm a second more rapid decline in 187Os/188Os associated with the KPB that is accompanied by a significant increase in Os concentrations. Complementary platinum (Pt) and iridium (Ir) concentration data indicate that the length scale of diagenetic remobilization of platinum group elements from the KPB is less than 1 m and does not obscure the pre-KPB decline in 187Os/188Os. Increases in bulk sediment Ir concentrations and decreases in bulk carbonate content that coincide with the Os isotope shift suggest that carbonate burial flux may have been lower during the initial decline in 187Os/188Os. We speculate that diminished carbonate burial rate may have been the result of ocean acidification caused by Deccan volcanism.

Geochemistry and isotopes at DSDP Leg 80 Holes

Thirty-eight samples from DSDP Sites 549 to 551 were analyzed for major and minor components and trace element abundances. Multivariate statistical analysis of geochemical data groups the samples into two major classes: an organic-carbon- rich group (> 1% TOC) containing high levels of marine organic matter and certain trace elements (Cu, Zn, V, Ni, Co, Ba, and Cr) and an organic-carbon-lean group depleted in these components. The greatest organic and trace metal enrichments occur in the uppermost Albian to Turanian sections of Sites 549 to 551. Carbon-isotopic values of bulk carbonate for the middle Cenomanian section of Site 550 (2.35 to 2.70 per mil) and the upper Cenomanian-Turonian sections of Sites 549 (3.35 to 4.47 per mil) and 551 (3.13 to 3.72 per mil) are similar to coeval values reported elsewhere in the region. The relatively heavy d13C values from Sites 549 and 551 indicate that this interval was deposited during the global "oceanic anoxic event" that occurred at the Cenomanian/Turonian boundary. Variation in the d18O of bulk carbonate for Section 550B-18-1 of middle Cenomanian age suggests that paleosalinity and/or paleotemperature variations may have occurred concurrently with periodic anoxia at this site. Climatically controlled increases in surface-water runoff may have caused surface waters to periodically freshen, resulting in stable salinity stratification

Paleomagnetism of DSDP Leg 80 holes

Paleomagnetic results from sediments acquired from the continental margin at DSDP Sites 548, 549, 550, and 551 are described. Where possible, the results were used to construct a polarity reversal stratigraphy for the sections sampled, thus enabling the biostratigraphic dating of the sediments to be refined. Several sections in this study were found to be suitable for magnetostratigraphic work, in particular the upper Paleocene to middle Eocene sediments from Site 549, which contained rich faunal assemblages. These sediments are underlain by a thick sequence of Cretaceous sediments that formed during the Long Cretaceous normal polarity interval. Sediments that formed during the later part of this magnetically quiet interval were also recovered at Site 550. Three short reverse polarity intervals were also recovered at this site; they lie directly over basement and are thought to represent a mixed-polarity interval of late Albian age. They may therefore provide important evidence concerning the age of the earliest sediments at this site. In addition, measurements of the magnetic susceptibility and intensity of remanent magnetism proved to be of interest. A significant decrease in the susceptibility and intensity values close to the early/middle Eocene boundary was noted at Sites 548 and 549. This decrease may be correlated with the results from Holes 400A and 401, which were drilled on DSDP Leg 48 in the northeast Bay of Biscay. The decrease may represent an abrupt reduction in the supply of terrigenous material at the end of the early Eocene, reflecting, perhaps, a change in sediment transport processes at that time

Primary composition, alteration, and origin of Cretaceous volcaniclastic rocks at DSDP Site 89-585

An upper Aptian to middle Albian series of volcaniclastic rocks more than 300 m thick was drilled at Site 585 in the East Mariana Basin. On the basis of textural and compositional (bulk-rock chemistry, primary and secondary mineral phases) evidence, the volcaniclastic unit is subdivided into a lower (below 830 m sub-bottom) and an upper (about 670-760 m) sequence; the boundary in the interval between is uncertain owing to lack of samples. The rocks are dominantly former vitric basaltic tuffs and minor lapillistones with lesser amounts of crystals and basaltic lithic clasts. They are mixed with shallow-water carbonate debris (ooids, skeletal debris), and were transported by mass flows to their site of deposition. The lower sequence is mostly plagioclase- and olivine-phyric with lesser amounts of Ti-poor clinopyroxene. Mineralogical and bulk-rock chemical data indicate a tholeiitic composition slightly more enriched than N-MORB (normal mid-ocean ridge basalt). Transport was by debris flows from shallow-water sites, as indicated by admixed ooids. Volcanogenic particles are chiefly moderately vesicular to nonvesicular blocky shards (former sideromelane) and less angular tachylite with quench plagioclase and pyroxene, indicating generation of volcanic clasts predominantly by spalling and breakage of submarine pillow and/or sheet-flow lavas. The upper sequence is mainly clinopyroxene- and olivine-phyric with minor plagioclase. The more Ti-rich clinopyroxene and the bulk-rock analyses show that the moderately alkali basaltic composition throughout is more mafic than the basal tholeiitic sequence. Transport was by turbidity currents. Rounded epiclasts of crystalline basalts are more common than in the lower sequence, and, together with the occurrence of oxidized olivine pseudomorphs and vesicular tachylite, are taken as evidence of derivation from eroded subaerially exposed volcanics. Former sideromelane shards are more vesicular than in the lower sequence; vesicularity exceeds 60 vol.% in some clasts. The dominant clastic process is interpreted to be by shallow-water explosive eruptions. All rocks have undergone low-temperature alteration; the dominant secondary phases are "palagonite," chlorite/smectite mixed minerals, analcite, and chabazite. Smectite, chlorite, and natrolite occur in minor amounts. Phillipsite is recognized as an early alteration product, now replaced by other zeolites. During alteration, the rocks have lost up to 50% of their Ca, compared with a fresh shard and fresh glass inclusions in primary minerals, but have gained much less K, Rb, and Ba than expected, indicating rapid deposition prior to significant seafloor weathering.

Tracing CO2 emissions in global value chains

This paper integrates two lines of research into a unified conceptual framework: trade in global value chains and embodied emissions. This allows both value added and emissions to be systematically traced at the country, sector, and bilateral levels through various production network routes. By combining value-added and emissions accounting in a consistent way, the potential environmental cost (amount of emissions per unit of value added) along global value chains can be estimated. Using this unified accounting method, we trace CO2 emissions in the global production and trade network among 41 economies in 35 sectors from 1995 to 2009, basing our calculations on the World Input–Output Database, and show how they help us to better understand the impact of cross-country production sharing on the environment.

Effects of RoHs and REACH regulations on firm-level production and export, and the role of global value chains : the cases of Malaysia and Vietnam

This paper uses firm-level data to examine the impact of foreign chemical safety regulations such as RoHS and REACH on the production costs and export performance of firms in Malaysia and Vietnam. This paper also investigates the role of global value chains in enhancing the likelihood that a firm complies with RoHS and REACH. We find that in addition to the initial setup costs for compliance, EU RoHS (REACH) implementation imposes on firms additional variable production costs by requiring additional labor and capital expenditures of around 57% (73%) of variable costs. We also find that compliance with RoHS and REACH significantly increases the probability of export and that compliance with EU RoHS and REACH helps firms enter a greater variety of countries. Furthermore, firms participating in global value chains have higher compliance with RoHS and REACH regulations, regardless of whether the firm is directly exporting, when the firm operates in upstream or downstream industries of the countries' supply chain.

A unified framework of trade in value added : physical, monetary, exchange rates, and GHG emissions

Koopman et al. (2014) developed a method to consistently decompose gross exports in value-added terms that accommodate infinite repercussions of international and inter-sector transactions. This provides a better understanding of trade in value added in global value chains than does the conventional gross exports method, which is affected by double-counting problems. However, the new framework is based on monetary input--output (IO) tables and cannot distinguish prices from quantities; thus, it is unable to consider financial adjustments through the exchange market. In this paper, we propose a framework based on a physical IO system, characterized by its linear programming equivalent that can clarify the various complexities relevant to the existing indicators and is proved to be consistent with Koopman's results when the physical decompositions are evaluated in monetary terms. While international monetary tables are typically described in current U.S. dollars, the physical framework can elucidate the impact of price adjustments through the exchange market. An iterative procedure to calculate the exchange rates is proposed, and we also show that the physical framework is also convenient for considering indicators associated with greenhouse gas (GHG) emissions.

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.

Theoretical Analysis of the No-Slip Boundary Condition Enforcement in SPH Methods

The aim of the present work is to provide an in-depth analysis of the most representative mirroring techniques used in SPH to enforce boundary conditions (BC) along solid profiles. We specifically refer to dummy particles, ghost particles, and Takeda et al. [Prog. Theor. Phys. 92 (1994), 939] boundary integrals. The analysis has been carried out by studying the convergence of the first- and second-order differential operators as the smoothing length (that is, the characteristic length on which relies the SPH interpolation) decreases. These differential operators are of fundamental importance for the computation of the viscous drag and the viscous/diffusive terms in the momentum and energy equations. It has been proved that close to the boundaries some of the mirroring techniques leads to intrinsic inaccuracies in the convergence of the differential operators. A consistent formulation has been derived starting from Takeda et al. boundary integrals (see the above reference). This original formulation allows implementing no-slip boundary conditions consistently in many practical applications as viscous flows and diffusion problems.

A set of canonical problems in sloshing. Part 0: Experimental setup and data processing

In a series of attempts to research and document relevant sloshing type phenomena, a series of experiments have been conducted. The aim of this paper is to describe the setup and data processing of such experiments. A sloshing tank is subjected to angular motion. As a result pressure registers are obtained at several locations, together with the motion data, torque and a collection of image and video information. The experimental rig and the data acquisition systems are described. Useful information for experimental sloshing research practitioners is provided. This information is related to the liquids used in the experiments, the dying techniques, tank building processes, synchronization of acquisition systems, etc. A new procedure for reconstructing experimental data, that takes into account experimental uncertainties, is presented. This procedure is based on a least squares spline approximation of the data. Based on a deterministic approach to the first sloshing wave impact event in a sloshing experiment, an uncertainty analysis procedure of the associated first pressure peak value is described.

Tranversal motion and flow structure of fully nonlinear streaks in a laminar boundary layer

Typical streak computations present in the literature correspond to linear streaks or to small amplitude nonlinear streaks computed using DNS or nonlinear PSE. We use the Reduced Navier-Stokes (RNS) equations to compute the streamwise evolution of fully non-linear streaks with high amplitude in a laminar flat plate boundary layer. The RNS formulation provides Reynolds number independent solutions that are asymptotically exact in the limit $Re \gg 1$, it requires much less computational effort than DNS, and it does not have the consistency and convergence problems of the PSE. We present various streak computations to show that the flow configuration changes substantially when the amplitude of the streaks grows and the nonlinear effects come into play. The transversal motion (in the wall normal-streamwise plane) becomes more important and strongly distorts the streamwise velocity profiles, that end up being quite different from those of the linear case. We analyze in detail the resulting flow patterns for the nonlinearly saturated streaks and compare them with available experimental results.