Spurious transients of projection methods in microflow simulations

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Utilização da modelagem sistêmica para analisar o efeito da variação climática no parâmetro produtividade da água na laranja Natal (Citrus sinensis L. Osbeck) no Estado de São Paulo

Pós-graduação em Agronomia (Irrigação e Drenagem) - FCA

A new predictive solar radiation numerical model

[EN]A predictive solar radiation numerical model is presented. Starting from the works of, a solar radiation numerical model is developed considering the terrain surface through 2-D adaptive meshes of triangles which are constructed using a refinement/derefinement procedure in accordance with the variations of terrain surface and albedo. The effect of shadows is considered in each time step. Solar radiation is first computed for clear-sky (CS) conditions and then, real-sky values are computed daily in terms of the CS index computed using all the observational data which are available for each day at several points of the studied zone…

Analisi numerica della combustione di gas a basso potere calorifico per applicazioni a microturbine a gas

Il presente lavoro si occupa dell’analisi numerica di combustione di gas a basso potere calorifico (gas di sintesi derivanti da pirolisi di biomasse). L’analisi è stata condotta su due principali geometrie di camera di combustione. La prima è un bruciatore sperimentale da laboratorio adatto allo studio delle proprietà di combustione del singas. Esso è introdotto in camera separatamente rispetto ad una corrente d’aria comburente al fine di realizzare una combustione non-premiscelata diffusiva in presenza di swirl. La seconda geometria presa in considerazione è la camera di combustione anulare installata sulla microturbina a gas Elliott TA 80 per la quale si dispone di un modello installato al banco al fine dell’esecuzione di prove sperimentali. I principali obbiettivi conseguiti nello studio sono stati la determinazione numerica del campo di moto a freddo su entrambe le geometrie per poi realizzare simulazioni in combustione mediante l’utilizzo di diversi modelli di combustione. In particolare è stato approfondito lo studio dei modelli steady laminar flamelet ed unsteady flamelet con cui sono state esaminate le distribuzioni di temperatura e delle grandezze tipiche di combustione in camera, confrontando i risultati numerici ottenuti con altri modelli di combustione (Eddy Dissipation ed ED-FR) e con i dati sperimentali a disposizione. Di importanza fondamentale è stata l’analisi delle emissioni inquinanti, realizzata per entrambe le geometrie, che mostra l’entità di tali emissioni e la loro tipologia. Relativamente a questo punto, il maggior interesse si sposta sui risultati ottenuti numericamente nel caso della microturbina, per la quale sono a disposizione misure di emissione ottenute sperimentalmente. Sempre per questa geometria è stato inoltre eseguito il confronto fra microturbina alimentata con singas a confronto con le prestazioni emissive ottenute con il gas naturale. Nel corso dei tre anni, l’esecuzione delle simulazioni e l’analisi critica dei risultati ha suggerito alcuni limiti e semplificazioni eseguite sulle griglie di calcolo realizzate per lo studio numerico. Al fine di eliminare o limitare le semplificazioni o le inesattezze, le geometrie dei combustori e le griglie di calcolo sono state migliorate ed ottimizzate. In merito alle simulazioni realizzate sulla geometria del combustore della microturbina Elliott TA 80 è stata condotta dapprima l’analisi numerica di combustione a pieno carico per poi analizzare le prestazioni ai carichi parziali. Il tutto appoggiandosi a tecniche di simulazione RANS ed ipotizzando alimentazioni a gas naturale e singas derivato da biomasse. Nell’ultimo anno di dottorato è stato dedicato tempo all’approfondimento e allo studio della tecnica Large Eddy Simulation per testarne una applicazione alla geometria del bruciatore sperimentale di laboratorio. In tale simulazione è stato implementato l’SGS model di Smagorinsky-Lilly completo di combustione con modelli flamelet. Dai risultati sono stati estrapolati i profili di temperatura a confronto con i risultati sperimentali e con i risultati RANS. Il tutto in diverse simulazioni a diverso valore del time-step imposto. L’analisi LES, per quanto migliorabile, ha fornito risultati sufficientemente precisi lasciando per il futuro la possibilità di approfondire nuovi modelli adatti all’applicazione diretta sulla MTG.

Techniques for Lagrangian modelling of dispersion in geophysical flows

Basic concepts and definitions relative to Lagrangian Particle Dispersion Models (LPDMs)for the description of turbulent dispersion are introduced. The study focusses on LPDMs that use as input, for the large scale motion, fields produced by Eulerian models, with the small scale motions described by Lagrangian Stochastic Models (LSMs). The data of two different dynamical model have been used: a Large Eddy Simulation (LES) and a General Circulation Model (GCM). After reviewing the small scale closure adopted by the Eulerian model, the development and implementation of appropriate LSMs is outlined. The basic requirement of every LPDM used in this work is its fullfillment of the Well Mixed Condition (WMC). For the dispersion description in the GCM domain, a stochastic model of Markov order 0, consistent with the eddy-viscosity closure of the dynamical model, is implemented. A LSM of Markov order 1, more suitable for shorter timescales, has been implemented for the description of the unresolved motion of the LES fields. Different assumptions on the small scale correlation time are made. Tests of the LSM on GCM fields suggest that the use of an interpolation algorithm able to maintain an analytical consistency between the diffusion coefficient and its derivative is mandatory if the model has to satisfy the WMC. Also a dynamical time step selection scheme based on the diffusion coefficient shape is introduced, and the criteria for the integration step selection are discussed. Absolute and relative dispersion experiments are made with various unresolved motion settings for the LSM on LES data, and the results are compared with laboratory data. The study shows that the unresolved turbulence parameterization has a negligible influence on the absolute dispersion, while it affects the contribution of the relative dispersion and meandering to absolute dispersion, as well as the Lagrangian correlation.

Stima dei parametri di modelli idrologici mediante ottimizzazione dell’utilità

Tradizionalmente, l'obiettivo della calibrazione di un modello afflussi-deflussi è sempre stato quello di ottenere un set di parametri (o una distribuzione di probabilità dei parametri) che massimizzasse l'adattamento dei dati simulati alla realtà osservata, trattando parzialmente le finalità applicative del modello. Nel lavoro di tesi viene proposta una metodologia di calibrazione che trae spunto dell'evidenza che non sempre la corrispondenza tra dati osservati e simulati rappresenti il criterio più appropriato per calibrare un modello idrologico. Ai fini applicativi infatti, può risultare maggiormente utile una miglior rappresentazione di un determinato aspetto dell'idrogramma piuttosto che un altro. Il metodo di calibrazione che viene proposto mira a valutare le prestazioni del modello stimandone l'utilità nell'applicazione prevista. Tramite l'utilizzo di opportune funzioni, ad ogni passo temporale viene valutata l'utilità della simulazione ottenuta. La calibrazione viene quindi eseguita attraverso la massimizzazione di una funzione obiettivo costituita dalla somma delle utilità stimate nei singoli passi temporali. Le analisi mostrano come attraverso l'impiego di tali funzioni obiettivo sia possibile migliorare le prestazioni del modello laddove ritenute di maggior interesse per per le finalità applicative previste.

Diskontinuierliche Galerkin-Verfahren für die operationelle Wettervorhersage

In dieser Arbeit wird ein neuer Dynamikkern entwickelt und in das bestehendernnumerische Wettervorhersagesystem COSMO integriert. Für die räumlichernDiskretisierung werden diskontinuierliche Galerkin-Verfahren (DG-Verfahren)rnverwendet, für die zeitliche Runge-Kutta-Verfahren. Hierdurch ist ein Verfahrenrnhoher Ordnung einfach zu realisieren und es sind lokale Erhaltungseigenschaftenrnder prognostischen Variablen gegeben. Der hier entwickelte Dynamikkern verwendetrngeländefolgende Koordinaten in Erhaltungsform für die Orographiemodellierung undrnkoppelt das DG-Verfahren mit einem Kessler-Schema für warmen Niederschlag. Dabeirnwird die Fallgeschwindigkeit des Regens, nicht wie üblich implizit imrnKessler-Schema diskretisiert, sondern explizit im Dynamikkern. Hierdurch sindrndie Zeitschritte der Parametrisierung für die Phasenumwandlung des Wassers undrnfür die Dynamik vollständig entkoppelt, wodurch auch sehr große Zeitschritte fürrndie Parametrisierung verwendet werden können. Die Kopplung ist sowohl fürrnOperatoraufteilung, als auch für Prozessaufteilung realisiert.rnrnAnhand idealisierter Testfälle werden die Konvergenz und die globalenrnErhaltungseigenschaften des neu entwickelten Dynamikkerns validiert. Die Massernwird bis auf Maschinengenauigkeit global erhalten. Mittels Bergüberströmungenrnwird die Orographiemodellierung validiert. Die verwendete Kombination ausrnDG-Verfahren und geländefolgenden Koordinaten ermöglicht die Behandlung vonrnsteileren Bergen, als dies mit dem auf Finite-Differenzenverfahren-basierendenrnDynamikkern von COSMO möglich ist. Es wird gezeigt, wann die vollernTensorproduktbasis und wann die Minimalbasis vorteilhaft ist. Die Größe desrnEinflusses auf das Simulationsergebnis der Verfahrensordnung, desrnParametrisierungszeitschritts und der Aufteilungsstrategie wirdrnuntersucht. Zuletzt wird gezeigt dass bei gleichem Zeitschritt die DG-Verfahrenrnaufgrund der besseren Skalierbarkeit in der Laufzeit konkurrenzfähig zurnFinite-Differenzenverfahren sind.

IMEX finite volume methods for the shallow water equations

Die Flachwassergleichungen (SWE) sind ein hyperbolisches System von Bilanzgleichungen, die adäquate Approximationen an groß-skalige Strömungen der Ozeane, Flüsse und der Atmosphäre liefern. Dabei werden Masse und Impuls erhalten. Wir unterscheiden zwei charakteristische Geschwindigkeiten: die Advektionsgeschwindigkeit, d.h. die Geschwindigkeit des Massentransports, und die Geschwindigkeit von Schwerewellen, d.h. die Geschwindigkeit der Oberflächenwellen, die Energie und Impuls tragen. Die Froude-Zahl ist eine Kennzahl und ist durch das Verhältnis der Referenzadvektionsgeschwindigkeit zu der Referenzgeschwindigkeit der Schwerewellen gegeben. Für die oben genannten Anwendungen ist sie typischerweise sehr klein, z.B. 0.01. Zeit-explizite Finite-Volume-Verfahren werden am öftersten zur numerischen Berechnung hyperbolischer Bilanzgleichungen benutzt. Daher muss die CFL-Stabilitätsbedingung eingehalten werden und das Zeitinkrement ist ungefähr proportional zu der Froude-Zahl. Deswegen entsteht bei kleinen Froude-Zahlen, etwa kleiner als 0.2, ein hoher Rechenaufwand. Ferner sind die numerischen Lösungen dissipativ. Es ist allgemein bekannt, dass die Lösungen der SWE gegen die Lösungen der Seegleichungen/ Froude-Zahl Null SWE für Froude-Zahl gegen Null konvergieren, falls adäquate Bedingungen erfüllt sind. In diesem Grenzwertprozess ändern die Gleichungen ihren Typ von hyperbolisch zu hyperbolisch.-elliptisch. Ferner kann bei kleinen Froude-Zahlen die Konvergenzordnung sinken oder das numerische Verfahren zusammenbrechen. Insbesondere wurde bei zeit-expliziten Verfahren falsches asymptotisches Verhalten (bzgl. der Froude-Zahl) beobachtet, das diese Effekte verursachen könnte.Ozeanographische und atmosphärische Strömungen sind typischerweise kleine Störungen eines unterliegenden Equilibriumzustandes. Wir möchten, dass numerische Verfahren für Bilanzgleichungen gewisse Equilibriumzustände exakt erhalten, sonst können künstliche Strömungen vom Verfahren erzeugt werden. Daher ist die Quelltermapproximation essentiell. Numerische Verfahren die Equilibriumzustände erhalten heißen ausbalanciert.rnrnIn der vorliegenden Arbeit spalten wir die SWE in einen steifen, linearen und einen nicht-steifen Teil, um die starke Einschränkung der Zeitschritte durch die CFL-Bedingung zu umgehen. Der steife Teil wird implizit und der nicht-steife explizit approximiert. Dazu verwenden wir IMEX (implicit-explicit) Runge-Kutta und IMEX Mehrschritt-Zeitdiskretisierungen. Die Raumdiskretisierung erfolgt mittels der Finite-Volumen-Methode. Der steife Teil wird mit Hilfe von finiter Differenzen oder au eine acht mehrdimensional Art und Weise approximniert. Zur mehrdimensionalen Approximation verwenden wir approximative Evolutionsoperatoren, die alle unendlich viele Informationsausbreitungsrichtungen berücksichtigen. Die expliziten Terme werden mit gewöhnlichen numerischen Flüssen approximiert. Daher erhalten wir eine Stabilitätsbedingung analog zu einer rein advektiven Strömung, d.h. das Zeitinkrement vergrößert um den Faktor Kehrwert der Froude-Zahl. Die in dieser Arbeit hergeleiteten Verfahren sind asymptotisch erhaltend und ausbalanciert. Die asymptotischer Erhaltung stellt sicher, dass numerische Lösung das "korrekte" asymptotische Verhalten bezüglich kleiner Froude-Zahlen besitzt. Wir präsentieren Verfahren erster und zweiter Ordnung. Numerische Resultate bestätigen die Konvergenzordnung, so wie Stabilität, Ausbalanciertheit und die asymptotische Erhaltung. Insbesondere beobachten wir bei machen Verfahren, dass die Konvergenzordnung fast unabhängig von der Froude-Zahl ist.

Theta burst transcranial magnetic stimulation is associated with increased EEG synchronization in the stimulated relative to unstimulated cerebral hemisphere

Theta burst transcranial magnetic stimulation (TBS) may induce behavioural changes that outlast the stimulation period. The neurophysiological basis of these behavioural changes are currently under investigation. Given the evidence that cortical information processing relies on transient synchronization and desynchronization of neuronal assemblies, we set out to test whether TBS is associated with changes of neuronal synchronization as assessed by surface EEG. In four healthy subjects one TBS train of 600 pulses (200 bursts, each burst consisting of 3 pulses at 30 Hz, repeated at intervals of 100 ms) was applied over the right frontal eye field and EEG synchronization was assessed in a time-resolved manner over 60 min by using a non-overlapping moving window. For each time step the linear cross-correlation matrix for six EEG channels of the right and for the six homotopic EEG channels of the left hemisphere were computed and their largest eigenvalues used to assess changes of synchronization. Synchronization was computed for broadband EEG and for the delta, theta, alpha, beta and gamma frequency bands. In all subjects EEG synchronization of the stimulated hemisphere was significantly and persistently increased relative to EEG synchronization of the unstimulated hemisphere. This effect occurred immediately after TBS for the theta, alpha, beta and gamma frequency bands and 10-20 min after TBS for broadband and delta frequency band EEG. Our results demonstrate that TBS is associated with increased neuronal synchronization of the cerebral hemisphere ipsilateral to the stimulation site relative to the unstimulated hemisphere. We speculate that enhanced synchronization interferes with cortical information processing and thus may be a neurophysiological correlate of the impaired behavioural performance detected previously.

Structural response and dynamics of fluid-structure-control interaction in wind turbine blades

Reducing the uncertainties related to blade dynamics by the improvement of the quality of numerical simulations of the fluid structure interaction process is a key for a breakthrough in wind-turbine technology. A fundamental step in that direction is the implementation of aeroelastic models capable of capturing the complex features of innovative prototype blades, so they can be tested at realistic full-scale conditions with a reasonable computational cost. We make use of a code based on a combination of two advanced numerical models implemented in a parallel HPC supercomputer platform: First, a model of the structural response of heterogeneous composite blades, based on a variation of the dimensional reduction technique proposed by Hodges and Yu. This technique has the capacity of reducing the geometrical complexity of the blade section into a stiffness matrix for an equivalent beam. The reduced 1-D strain energy is equivalent to the actual 3-D strain energy in an asymptotic sense, allowing accurate modeling of the blade structure as a 1-D finite-element problem. This substantially reduces the computational effort required to model the structural dynamics at each time step. Second, a novel aerodynamic model based on an advanced implementation of the BEM(Blade ElementMomentum) Theory; where all velocities and forces are re-projected through orthogonal matrices into the instantaneous deformed configuration to fully include the effects of large displacements and rotation of the airfoil sections into the computation of aerodynamic forces. This allows the aerodynamic model to take into account the effects of the complex flexo-torsional deformation that can be captured by the more sophisticated structural model mentioned above. In this thesis we have successfully developed a powerful computational tool for the aeroelastic analysis of wind-turbine blades. Due to the particular features mentioned above in terms of a full representation of the combined modes of deformation of the blade as a complex structural part and their effects on the aerodynamic loads, it constitutes a substantial advancement ahead the state-of-the-art aeroelastic models currently available, like the FAST-Aerodyn suite. In this thesis, we also include the results of several experiments on the NREL-5MW blade, which is widely accepted today as a benchmark blade, together with some modifications intended to explore the capacities of the new code in terms of capturing features on blade-dynamic behavior, which are normally overlooked by the existing aeroelastic models.

Adaptive sensing for target tracking applications

In a statistical inference scenario, the estimation of target signal or its parameters is done by processing data from informative measurements. The estimation performance can be enhanced if we choose the measurements based on some criteria that help to direct our sensing resources such that the measurements are more informative about the parameter we intend to estimate. While taking multiple measurements, the measurements can be chosen online so that more information could be extracted from the data in each measurement process. This approach fits well in Bayesian inference model often used to produce successive posterior distributions of the associated parameter. We explore the sensor array processing scenario for adaptive sensing of a target parameter. The measurement choice is described by a measurement matrix that multiplies the data vector normally associated with the array signal processing. The adaptive sensing of both static and dynamic system models is done by the online selection of proper measurement matrix over time. For the dynamic system model, the target is assumed to move with some distribution and the prior distribution at each time step is changed. The information gained through adaptive sensing of the moving target is lost due to the relative shift of the target. The adaptive sensing paradigm has many similarities with compressive sensing. We have attempted to reconcile the two approaches by modifying the observation model of adaptive sensing to match the compressive sensing model for the estimation of a sparse vector.

Predicting the Mechanical Properties of Carbon-Based Materials Using Molecular Dynamics

The development of innovative carbon-based materials can be greatly facilitated by molecular modeling techniques. Although molecular modeling has been used extensively to predict elastic properties of materials, modeling of more complex phenomenon such as fracture has only recently been possible with the development of new force fields such as ReaxFF, which is used in this work. It is not fully understood what molecular modeling parameters such as thermostat type, thermostat coupling, time step, system size, and strain rate are required for accurate modeling of fracture. Selection of modeling parameters to model fracture can be difficult and non-intuitive compared to modeling elastic properties using traditional force fields, and the errors generated by incorrect parameters may be non-obvious. These molecular modeling parameters are systematically investigated and their effects on the fracture of well-known carbon materials are analyzed. It is determined that for coupling coefficients of 250 fs and greater do not result in substantial differences in the stress-strain response of the materials using any thermostat type. A time step of 0.5 fs of smaller is required for accurate results. Strain rates greater than 2.2 ns-1 are sufficient to obtain repeatable results with slower strain rates for the materials studied. The results of this study indicate that further refinement of the Chenoweth parameter set is required to accurately predict the mechanical response of carbon-based systems. The ReaxFF has been used extensively to model systems in which bond breaking and formation occur. In particular ReaxFF has been used to model reactions of small molecules. Some elastic and fracture properties have been successfully modeled using ReaxFF in materials such as silicon and some metals. However, it is not clear if current parameterizations for ReaxFF are able to accurately reproduce the elastic and fracture properties of carbon materials. The stress-strain response of a new ReaxFF parameterization is compared to the previous parameterization and density functional theory results for well-known carbon materials. The new ReaxFF parameterization makes xv substantial improvements to the predicted mechanical response of carbon materials, and is found to be suitable for modeling the mechanical response of carbon materials. Finally, a new material composed of carbon nanotubes within an amorphous carbon (AC) matrix is modeled using the ReaxFF. Various parameters that may be experimentally controlled are investigated such as nanotube bundling, comparing multi-walled nanotube with single-walled nanotubes, and degree of functionalization of the nanotubes. Elastic and fracture properties are investigated for the composite systems and compared to results of pure-nanotube and pure-AC models. It is found that the arrangement of the nanotubes and degree of crosslinking may substantially affect the properties of the systems, particularly in the transverse directions.

Land-atmosphere coupling over North America in CRCM5

Land-atmosphere coupling and its impact on extreme precipitation and temperature events over North America are studied using the fifth generation of the Canadian Regional Climate Model (CRCM5). To this effect, two 30 year long simulations, spanning the 1981–2010 period, with and without land-atmosphere coupling, have been performed with CRCM5, driven by the European Centre for Medium-Range Weather Forecasts reanalysis at the boundaries. In the coupled simulation, the soil moisture interacts freely with the atmosphere at each time step, while in the uncoupled simulation, soil moisture is replaced with its climatological value computed from the coupled simulation, thus suppressing the soil moisture-atmosphere interactions. Analyses of the coupled and uncoupled simulations, for the summer period, show strong soil moisture-temperature coupling over the Great Plains, consistent with previous studies. The maxima of soil moisture-precipitation coupling is more spread out and covers the semiarid regions of the western U.S. and parts of the Great Plains. However, the strength of soil moisture-precipitation coupling is found to be generally weaker than that of soil moisture-temperature coupling. The study clearly indicates that land-atmosphere coupling increases the interannual variability of the seasonal mean daily maximum temperature in the Great Plains. Land-atmosphere coupling is found to significantly modulate selected temperature extremes such as the number of hot days, frequency, and maximum duration of hot spells over the Great Plains. Results also suggest additional hot spots, where soil moisture modulates extreme events. These hot spots are located in the southeast U.S. for the hot days/hot spells and in the semiarid regions of the western U.S. for extreme wet spells. This study thus demonstrates that climatologically wet/dry regions can become hot spots of land-atmosphere coupling when the soil moisture decreases/increases to an intermediate transitional level where evapotranspiration becomes moisture sensitive and large enough to affect the climate.

Oxygen and carbon isotope stratigraphy of ODP Holes 107-653A and 107-654A

Stable isotope analysis of two species (or groups of species) of planktonic foraminifers: Globigerinoides ruber (or G. obliquus and G. obliquus extremus) and Globigerina bulloides (or G. falconensis and G. obesa) from ODP Hole 653A and Site 654 in the Tyrrhenian basin, records the Pliocene-Pleistocene glacial history of the Northern Hemisphere. The overall increase in mean d18O values through the interval 4.6-0.08 Ma is 1.7 per mil for G. bulloides and 1.5 per mil for G. ruber. The time interval 3.1-2.5 Ma corresponds to an important phase of 18O enrichment for planktonic foraminifers. In this interval, glacial d18O values of both species G. bulloides and G. ruber increase by about l per mil, this increase being more progressive for G. ruber than for G. bulloides. The increase of interglacial d18O values is higher for G. bulloides (1.5 per mil) than for the Gruber group (1 per mil). These data suggest a more pronounced seasonal stratification of the water masses during interglacial phases. Large positive d18O fluctuations of increasing magnitude are also recorded at 2.25 and 2.15 Ma by G bulloides and appear to be diachronous with those of Site 606 in the Atlantic Ocean. Other events of increasing d18O values are recorded between 1.55 and 1.3 Ma, at 0.9 Ma, 0.8 Ma, and near 0.34 Ma. In the early Pliocene the d18O variability recorded by the planktonic species G. bulloides was higher in the Mediterranean than in the Atlantic at the same latitude. This suggests that important cyclic variations in the water budget of the Mediterranean occurred since that time. Step increases in the d18O variability are synchronous with those of the open ocean at 0.9 and 0.34 Ma. The higher variability as well as the higher amplitude of the peaks of 18O enrichment may be partly accounted for by increase of dryness over the Mediterranean area. In particular the high amplitude d18O fluctuations recorded between 3.1 and 2.1 Ma are correlated with the onset of a marked seasonal contrast and a summer dryness, revealed by pollen analyses. Strong fluctuations towards d13C values higher than modern ones are recorded by the G. ruber group species before 1.7 Ma and suggest a high production of phytoplankton. When such episodes of high primary production are correlated with episodes of decreasing 13C content of G. bulloides, they are interpreted as the consequence of a higher stratification of the upper water masses resulting itself from a marked seasonality. Such episodes occur between 4.6 and 4.05 Ma, 3.9 and 3.6 Ma, and 3.25 and 2.66 Ma. The interval 2.66-1.65 Ma corresponds to a weakening of the stratification of the upper water layers. This may be related to episodes of cooling and increasing dryness induced by the Northern Hemisphere Glaciations. The Pleistocene may have been a less productive period. The transition from highly productive to less productive surface waters also coincides with a new step increase in dryness and cooling, between 1.5 and 1.3 Ma. The comparison of the 13C records of G ruber and G. bulloides in fact suggests that a high vertical convection became a dominant feature after 2.6 Ma. Increases in the nutrient input and the stratification of the upper water masses may be suspected, however, during short episodes near 0.86 Ma (isotopic stage 25), 0.57-0.59 Ma (isotopic stage 16), 0.49 Ma (isotopic stage 13), 0.4-0.43 Ma (isotopic stage 11), and 0.22 and 0.26 Ma (part of isotopic stage 7 and transition 7/8). In fact, changes in the C02 balance within the different water masses of the Tyrrhenian basin as well as in the local primary production did not follow the general patterns of the open ocean.

Ages and geochemical and geophysical characteristics of sediment core MD99-2294, Lofoten Contourite Drift

A sediment core from the Lofoten Contourite Drift on the continental slope off Northern Norway, proximal to the former Vestfjorden-Trsnadjupet Ice Stream, details the development, variability and decline of marine margins of the northwestern Fennoscandian Ice Sheet during the time interval 25.3-14 cal ka BP, including the Last Glacial Maximum and onset of the deglaciation based on high-resolution IRD records. From the core interval between 25.3 and 17.7 cal ka BP we report data points with a mean time step of 10 years, between 17.7 cal ka BP and the Holocene time steps are typically 50 years. The core is divided into 7 informal ice-rafted debris (IRD) zones based on the variations in IRD including 7 major IRD maxima (A-G), inferred to represent periods of high iceberg production. Petrological identification reveals dominance of crystalline IRD (monocrystalline, plutonic and metamorphic rock fragments) accounting for 75-80% of total IRD assemblages, while sedimentary fragments generally account for 15-20%. The crystalline fragments (including eclogite and mangerite from a nearby terrestrial source) increase across the IRD peaks while the sedimentary fragments remain constant. This points to the importance of erosional products from icebergs originating from fast-flowing paleo-ice streams including the Vestfjorden-Trsnadjupet Ice Stream draining from the Fennoscandian mainland during the IRD maxima periods. Increased temperature of the adjacent surface water masses was probably an important external forcing factor on the Fennoscandian Ice Sheet behavior because some IRD maxima and plumite deposition from meltwater plumes post-date periods of increased sea surface temperatures. The peak IRD depositions occur in centennial and millennial time cycles (~200, 1030 and 3900 year) indicating some external forcing by solar variation. Both mechanisms could explain the observed synchronous instability of the northwestern Fennoscandian Ice Sheet to other European Ice Sheets.