Stable oxygen isotope record and Calcium concentrations of sediment core GeoB3910-2

There is increasing evidence that the preceding Holocene climate was as unstable as the last glacial period, although variations occurred at much lower amplitudes. However, low-latitude climate records that confirm this variability are sparse. Here we present a radiocarbon-dated Holocene marine record from the tropical western Atlantic. Aragonite dissolution derived from the degree of preservation of the pteropod Limacina inflata records changes in the corrosiveness of the bottom water at the core site due to the changing influence of northern versus southern water masses. The delta18O difference between the shallow-living planktonic foraminifera Globigerinoides sacculifer and the deep-living Globorotalia tumida is used as proxy for changes in the vertical stratification of the surface water, hence the trade wind strength at this latitude. We compared our data to high-latitude records of the North Atlantic region. A good agreement is found between the aragonite dissolution and the strength in the Island-Scotland Overflow Water, which contributes significantly to the North Atlantic Deep Water. This suggests that large-scale variations in the Atlantic thermohaline circulation occurred throughout the Holocene. Concurrently, the comparison of our Delta delta18O with the GISP2 glaciochemical records points to global Holocene atmospheric reorganizations seen in both the tropics and high northern latitudes.

Stable isotope ratios and organic carbon accumulation rates of late Pleistocene sediments of ODP Hole 117-724C

Stable isotope records of coexisting benthic foraminifers Uvigerina spp. and Cibicidoides spp. and planktonic G. ruber (white variety) from Site 724 are used to study the late Pleistocene evolution of surface and intermediate water hydrography (593 m water depth) at the Oman Margin. Glacial-interglacial d18O amplitudes recorded by the benthic foraminifers are reduced when compared to the estimated mean ocean changes of d18Oseawater . Epibenthic d13C remains at its modern level or is increased during glacial times. This implies that Red Sea outflow waters which are enriched in d18Oseawater and d13C (Sum CO2) have been replaced during glacial periods by intermediate waters still positive in d13C (Sum CO2) but more negative in d18Oseawater. Glacial-interglacial amplitudes of the planktonic d18O record exceed those of the mean ocean d18Oseawater variation and imply decreased surface water temperatures (SST) during glacial times. Throughout most of the records these cooling events correlate with enhanced rates of carbon accumulation. However, both negative (colder) SST and positive Corg accumulation rate anomalies do not correlate with potential physical upwelling maxima as inferred from the orbital monsoon index. This is in conflict with the established hypothesis that upwelling in the estern Arabia Sea should be strongest during maxima of the southwest monsoon.

Oxygen and Carbon stable isotope ratios for Site U1334

The Oligocene-Miocene transition (OMT) (~23 Ma) is interpreted as a transient global cooling event, associated with a large-scale Antarctic ice sheet expansion. Here we present a 2.23 Myr long high-resolution (~3 kyr) benthic foraminiferal oxygen and carbon isotope (d18O and d13C) record from Integrated Ocean Drilling Program Site U1334 (eastern equatorial Pacific Ocean), covering the interval from 21.91 to 24.14 Ma. To date, five other high-resolution benthic foraminiferal stable isotope stratigraphies across this time interval have been published, showing a ~1 per mil increase in benthic foraminiferal d18O across the OMT. However, these records are still few and spatially limited and no clear understanding exists of the global versus local imprints. We show that trends and the amplitudes of change are similar at Site U1334 as in other high-resolution stable isotope records, suggesting that these represent global deep water signals. We create a benthic foraminiferal stable isotope stack across the OMT by combining Site U1334 with records from ODP Sites 926, 929, 1090, 1264, and 1218 to best approximate the global signal. We find that isotopic gradients between sites indicate interbasinal and intrabasinal variabilities in deep water masses and, in particular, note an offset between the equatorial Atlantic and the equatorial Pacific, suggesting that a distinct temperature gradient was present during the OMT between these deep water masses at low latitudes. A convergence in the d18O values between infaunal and epifaunal species occurs between 22.8 and 23.2 Ma, associated with the maximum d18O excursion at the OMT, suggesting climatic changes associated with the OMT had an effect on interspecies offsets of benthic foraminifera. Our data indicate a maximum glacioeustatic sea level change of ~50 m across the OMT.

Seismic Data from MSM21/4 to Knipovich Ridge on the western Svalbard margin

The seismic data were acquired north of the Knipovich Ridge on the western Svalbard margin during cruise MSM21/4. They were recorded using a Geometrics GeoEel streamer of either 120 channels (profiles p100-p208) or 88 channels (profiles p300-p805) with a group spacing of 1.56 m and a sampling rate of 2 kHz. A GI-Gun (2×1.7 l) with a main frequency of ~150 Hz was used as a source and operated at a shot interval of 6-8 s. Processing of profiles p100-p208 and p600-p805: Positions for each channel were calculated by backtracking along the profiles from the GI-Gun GPS positions. The shot gathers were analyzed for abnormal amplitudes below the seafloor reflection by comparing neighboring traces in different frequency bands within sliding time windows. To suppress surface-generated water noise, a tau-p filter was applied in the shot gather domain. Common mid-point (CMP) profiles were then generated through crooked-line binning with a CMP spacing of 1.5625 m. A zero-phase band-pass filter with corner frequencies of 60 Hz and 360 Hz was applied to the data. Based on regional velocity information from MCS data [Sarkar, 2012], an interpolated and extrapolated 3D interval velocity model was created below the digitized seafloor reflection of the high-resolution streamer data. This velocity model was used to apply a CMP stack and an amplitude-preserving Kirchhoff post-stack time migration. Processing of profiles p400-p500: Data were sampled at 0.5 ms and sorted into common midpoint (CMP) domain with a bin spacing of 5 m. Normal move out correction was carried out with a velocity of 1500 m s-1 and an Ormsby bandpass filter with corner frequencies at 40, 80, 600 and 1000 Hz was applied. The data were time migrated using the water velocity.

Ecosystem regimes and responses in a coupled ancient lake system from MIS 5b to present: the diatom record of lakes Ohrid and Prespa

We reconstruct the aquatic ecosystem interactions since the last interglacial period in the oldest, most diverse, hydrologically connected European lake system, by using palaeolimnological diatom and selected geochemistry data from Lake Ohrid “DEEP site” core and equivalent data from Lake Prespa core, Co1215. Driven by climate forcing, the lakes experienced two adaptive cycles during the last 92 ka: "interglacial and interstadial" and "glacial" cycle. The short-term ecosystems reorganizations, e.g. regime shifts within these cycles substantially differ between the lakes, as evident from the inferred amplitudes of variation. The deeper Lake Ohrid shifted between ultra oligo- and oligotrophic regimes in contrast to the much shallower Lake Prespa, which shifted from a deeper, (oligo-) mesotrophic to a shallower, eutrophic lake and vice versa. Due to the high level of ecosystem stability (e.g. trophic state, lake level), Lake Ohrid appears relatively resistant to external forcing, such as climate and environmental change. Recovering in a relatively short time from major climate change, Lake Prespa is a resilient ecosystem. At the DEEP site, the decoupling between the lakes' response to climate change is marked in the prolonged and gradual changes during the MIS 5/4 and 2/1 transitions. These response differences and the lakes' different physical and chemical properties may limit the influence of Lake Prespa on Lake Ohrid. Regime shifts of Lake Ohrid due to potential hydrological change in Lake Prespa are not evident in the data presented here. Moreover, a complete collapse of the ecosystems functionality and loss of their diatom communities did not happen in either lake for the period presented in the study.

Idiopathic multifocal choroiditis/punctate inner choroidopathy with acute photoreceptor loss or dysfunction out of proportion to clinically visible lesions.

PURPOSE To report acute/subacute vision loss and paracentral scotomata in patients with idiopathic multifocal choroiditis/punctate inner choroidopathy due to large zones of acute photoreceptor attenuation surrounding the chorioretinal lesions. METHODS Multimodal imaging case series. RESULTS Six women and 2 men were included (mean age, 31.5 ± 5.8 years). Vision ranged from 20/20-1 to hand motion (mean, 20/364). Spectral domain optical coherence tomography demonstrated extensive attenuation of the external limiting membrane, ellipsoid and interdigitation zones, adjacent to the visible multifocal choroiditis/punctate inner choroidopathy lesions. The corresponding areas were hyperautofluorescent on fundus autofluorescence and were associated with corresponding visual field defects. Full-field electroretinogram (available in three cases) showed markedly decreased cone/rod response, and multifocal electroretinogram revealed reduced amplitudes and increased implicit times in two cases. Three patients received no treatment, the remaining were treated with oral corticosteroids (n = 4), oral acyclovir/valacyclovir (n = 2), intravitreal/posterior subtenon triamcinolone acetate (n = 3), and anti-vascular endothelial growth factor (n = 2). Visual recovery occurred in only three cases of whom two were treated. Varying morphological recovery was found in six cases, associated with decrease in hyperautofluorescence on fundus autofluorescence. CONCLUSION Multifocal choroiditis/punctate inner choroidopathy can present with transient or permanent central photoreceptor attenuation/loss. This presentation is likely a variant of multifocal choroiditis/punctate inner choroidopathy with chorioretinal atrophy. Associated changes are best evaluated using multimodal imaging.

Seismic modelling and ava analysis of hidrocarbon traps

In this study, forward seismic modelling of four geological models with Hydrocarbon (HC) traps were performed by ray tracing method to produce synthetic seismogram of each model. The idea is to identify the Hydrocarbon Indicators (HCI‟s) such as bright spot, flat spot, dim spot and Bottom Simulating Reflector (BSR) in the synthethic seismogram. The modelling was performed in DISCO/FOCUS 5.0 seismic data processing programme. Strong positive and negative reflection amplitudes and some artifact reflection horizons were observed on produced seismograms due to rapid changes in subsurface velocity and geometry respectively Additionally, Amplitude-versus-angle (AVA) curves of each HCIs was calculated by the Crewes Zoeppritz Explorer programme. AVA curves show that how the reflection coefficients change with the density and the P and S wave velocities of each layer such as oil, gas, gas hydrate or water saturated sediments. Due to AVA curves, an increase in reflection amplitude with incident angle of seismic waves corresponds to an indicator of a hydrocarbon reservoir

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.

Estudio de la variación del contenido de humedad en el Pinus pinaster Ait. por medio de la técnica no destructiva del georradar / Assessment of the variation of the moisture content in the Pinus pinaster Ait using the non destructive gpr technique

La variación en el contenido de humedad (CH) tiene una influencia significativa tanto en las propiedades físico- químicas de la madera, como en sus propiedades electromagnéticas y por tanto afecta a las características de la propagación de las ondas. En concreto, en este trabajo se estudia la capacidad del georradar (GR) empleando una antena de 1.6GHz de frecuencia central para registrar las variaciones que se producen en la velocidad y en la amplitud de las ondas electromagnéticas cuando se propagan en unas viguetas de madera de Pinus pinaster Ait de uso estructural cuyo CH va disminuyendo. Se ha comprobado como cuando el CH descendía la velocidad de propagación y las amplitudes, tanto de la onda directa como de la reflejada aumentaba. Los altos factores de correlación encontrados demuestran que el GR es una técnica capaz de evaluar, de forma no destructiva, el CH de la madera de uso estructural. The moisture content variations in wood have a significant influence in wood?s physicochemical properties, as well as in its electromagnetic properties and to specific effects upon waves? characteristics. In particular, this paper focuses on the analysis of the Ground-penetrating Radar?s (GPR) using an antenna of 1.6 GHz central frequency capacity to register the velocity and the amplitude of the electromagnetic waves? variation during the drying process of Pinus pinaster Ait timber joists. The results showed that when timber MC descended, the propagation velocity and amplitude of both the direct and the reflected wave increased. The high correlation found between the variables studied demonstrates GPR efficiency and the innovative application of this technique as a non-destructive evaluation tool for timber structures, particularly when studying its moisture content.

ALIFE: Caracterización Térmica De La Fermentación Del Café

La fermentación del café es una de las operaciones críticas en el proceso de beneficiado por su impacto en la calidad organoléptica del producto final. De forma general la mayor parte de las explotaciones carecen de tecnificación alguna del fermentador, no teniendo el operario ninguna información objetiva sobre la evolución de los parámetros físico-químicos de control del proceso que le ayuden a la toma de decisión. El objetivo de este trabajo es la integración multi-distribuida de sensores de bajo coste y alta fiabilidad para supervisión de la operación de fermentación del café en origen. Concretamente se ha centrado en la caracterización térmica del proceso de fermentación, mediante la instalación multidistribuida de sensores inalámbricos de temperatura (TurboTag ®) en fermentadores de café en Colombia. El análisis: temporal de la evolución de la temperatura a lo largo del proceso de fermentación y complejo mediante diagrama de fases o de cinética de la evolución de las temperaturas, permiten en este trabajo caracterizar térmicamente fermentadores industriales y prototipos plásticos. Se ha constatado la alta heterogeneidad de las fermentaciones tanto intra-lote (amplitudes térmicas de hasta 8.5 ºC y temperaturas en algunos puntos de los fermentadores de 17 ºC) como inter-lote en fermentaciones sujetas a condiciones exógenas y endógenas muy variables

Reduced order adaptive models for systems of PDEs using POD

In this work we propose a method to accelerate time dependent numerical solvers of systems of PDEs that require a high cost in computational time and memory. The method is based on the combined use of such numerical solver with a proper orthogonal decomposition, from which we identify modes, a Galerkin projection (that provides a reduced system of equations) and the integration of the reduced system, studying the evolution of the modal amplitudes. We integrate the reduced model until our a priori error estimator indicates that our approximation in not accurate. At this point we use again our original numerical code in a short time interval to adapt the POD manifold and continue then with the integration of the reduced model. Application will be made to two model problems: the Ginzburg-Landau equation in transient chaos conditions and the two-dimensional pulsating cavity problem, which describes the motion of liquid in a box whose upper wall is moving back and forth in a quasi-periodic fashion. Finally, we will discuss a way of improving the performance of the method using experimental data or information from numerical simulations

Reduced order modeling of three dimensional external aerodynamic flows

A method is presented to construct computationally efficient reduced-order models (ROMs) of three-dimensional aerodynamic flows around commercial aircraft components. The method is based on the proper orthogonal decomposition (POD) of a set of steady snapshots, which are calculated using an industrial solver based on some Reynolds averaged Navier-Stokes (RANS) equations. The POD-mode amplitudes are calculated by minimizing a residual defined from the Euler equations, even though the snapshots themselves are calculated from viscous equations. This makes the ROM independent of the peculiarities of the solver used to calculate the snapshots. Also, both the POD modes and the residual are calculated using points in the computational mesh that are concentrated in a close vicinity of the aircraft, which constitute a much smaller number than the total number of mesh points. Despite these simplifications, the method provides quite good approximations of the flow variables distributions in the whole computational domain, including the boundary layer attached to the aircraft surface and the wake. Thus, the method is both robust and computationally efficient, which is checked considering the aerodynamic flow around a horizontal tail plane, in the transonic range 0.4?Mach number?0.8, ?3°?angle of attack?3°.

Nonlinear response to ultrasound of encapsulated microbubbles

The acoustic backscatter of encapsulated gas-filled microbubbles immersed in a weak compressible liquid and irradiated by ultrasound fields of moderate to high pressure amplitudes is investigated theoretically. The problem is formulated by considering, for the viscoelastic shell of finite thickness, an isotropic hyperelastic neo-Hookean model for the elastic contribution in addition to a Newtonian viscous component. First and second harmonic scattering cross-sections have been evaluated and the quantitative influence of the driving pressure amplitude on the harmonic resonance frequencies for different initial equilibrium bubble sizes and for different encapsulating physical properties has been determined. Conditions for optimal second harmonic imaging have been also investigated and some regions in the parameters space where the second harmonic intensity is dominant over the fundamental have been identified. Results have been obtained for albumin, lipid and polymer encapsulating shells, respectively.

Low-dimensional modeling of streaks in a wedge flow boundary layer

This paper is concerned with the low dimensional structure of optimal streaks in a wedge flow boundary layer, which have been recently shown to consist of a unique (up to a constant factor) three-dimensional streamwise evolving mode, known as the most unstable streaky mode. Optimal streaks exhibit a still unexplored/unexploited approximate self-similarity (not associated with the boundary layer self-similarity), namely the streamwise velocity re-scaled with their maximum remains almost independent of both the spanwise wavenumber and the streamwise coordinate; the remaining two velocity components instead do not satisfy this property. The approximate self-similar behavior is analyzed here and exploited to further simplify the description of optimal streaks. In particular, it is shown that streaks can be approximately described in terms of the streamwise evolution of the scalar amplitudes of just three one-dimensional modes, providing the wall normal profiles of the streamwise velocity and two combinations of the cross flow velocity components; the scalar amplitudes obey a singular system of three ordinary differential equations (involving only two degrees of freedom), which approximates well the streamwise evolution of the general streaks.

Asymptotic Description of Flutter Amplitude Saturation by Nonlinear Friction Forces

The computation of the non-linear vibration dynamics of an aerodynamically unstable bladed-disk is a formidable numerical task, even for the simplified case of aerodynamic forces assumed to be linear. The nonlinear friction forces effectively couple dif- ferent travelling waves modes and, in order to properly elucidate the dynamics of the system, large time simulations are typically required to reach a final, saturated state. Despite of all the above complications, the output of the system (in the friction microslip regime) is basically a superposition of the linear aeroelastic un- stable travelling waves, which exhibit a slow time modulation that is much longer than the elastic oscillation period. This slow time modulation is due to both, the small aerodynamic effects and the small nonlinear friction forces, and it is crucial to deter- mine the final amplitude of the flutter vibration. In this presenta- tion we apply asymptotic techniques to obtain a new simplified model that captures the slow time dynamics of the amplitudes of the travelling waves. The resulting asymptotic model is very re- duced and extremely cheap to simulate, and it has the advantage that it gives precise information about the characteristics of the nonlinear friction models that actually play a role in the satura- tion of the vibration amplitude.