Metodología de simulación numérica de la transferencia de calor por ebullición en película en torno a superficies esféricas refrigeradas por líquido saturado en convección libre: aplicación a reactores nucleares de lecho de bolas (PBR)

En la presente Tesis se realizó un análisis numérico, usando el código comercial Ansys-Fluent, de la refrigeración de una bola de combustible de un reactor de lecho de bolas (PBR, por sus siglas en inglés), ante un escenario de emergencia en el cual el núcleo sea desensamblado y las bolas se dejen caer en una piscina de agua, donde el mecanismo de transferencia de calor inicialmente sería la ebullición en película, implicando la convección y la radiación al fluido. Previamente se realizaron pruebas de validación, comparando los resultados numéricos con datos experimentales disponibles en la literatura para tres geometrías diferentes, lo cual permitió seleccionar los esquemas y modelos numéricos con mejor precisión y menor costo computacional. Una vez identificada la metodología numérica, todas las pruebas de validación fueron ampliamente satisfactorias, encontrándose muy buena concordancia en el flujo de calor promedio con los datos experimentales. Durante estas pruebas de validación se lograron caracterizar numéricamente algunos parámetros importantes en la ebullición en película con los cuales existen ciertos niveles de incertidumbre, como son el factor de acoplamiento entre convección y radiación, y el factor de corrección del calor latente de vaporización. El análisis térmico de la refrigeración de la bola del reactor por ebullición en película mostró que la misma se enfría, a pesar del calor de decaimiento, con una temperatura superficial de la bola que desciende de forma oscilatoria, debido al comportamiento inestable de la película de vapor. Sin embargo, la temperatura de esta superficie tiene una buena uniformidad, notándose que las áreas mejor y peor refrigeradas están localizadas en la parte superior de la bola. Se observó la formación de múltiples domos de vapor en diferentes posiciones circunferenciales, lo cual causa que el área más caliente de la superficie se localice donde se forman los domos más grandes. La separación entre los domos de vapor fue consistente con la teoría hidrodinámica, con la adición de que la separación entre domos se reduce a medida que evolucionan y crecen, debido a la curvatura de la superficie. ABSTRACT A numerical cooling analysis of a PBR fuel pebble, after an emergency scenario in which the nucleus disassembly is made and the pebbles are dropped into a water pool, transmitting heat by film boiling, involving convection and radiation to the fluid, is carried out in this Thesis. First, were performed validation tests comparing the numerical results with experimental works available for three different geometries, which allowed the selection of numerical models and schemes with better precision and lower computational cost. Once identified the numerical methodology, all validation tests were widely satisfactory, finding very good agreement with experimental works in average heat flux. During these validation tests were achieved numerically characterize some important parameters in film boiling with which there are certain levels of uncertainty, such as the coupling factor between convection and radiation, and the correction factor of the latent heat of vaporization. The thermal analysis of pebble cooling by film boiling shows that despite its decay heat, cooling occurs, with pebble surface temperature descending from an oscillatory manner, due to the instability of the vapor film. However, the temperature of this surface has a good uniformity, noting that the best and worst refrigerated area is located at the top of the pebble. The formation of multiple vapor domes at different circumferential positions is observed, which cause that the hottest area of the surface was located where biggest vapor domes were formed. The separation between vapor domes was consistent with the hydrodynamic theory, with the addition that the separation is reduced as the vapor dome evolves and grows, due to the surface curvature.

Some aspects of lubrication in heavy regimes: thermal effects, stability and turbulence

In this work, a combination of numerical methods applied to thermohydrodynamic lubrication problems with cavitation is presented. It should be emphasized the difficulty of the nonlinear mathematical coupled model involving a free boundary problem, but also the simplicity of the algorithms employed to solve it. So, finite element discretizations for the hydrodynamic and thermal equations combined with upwind techniques for the convection terms and duality methods for nonlinear features are proposed. Additionally, a model describing the movement of the shaft is provided. Considering the shaft as a rigid body this model will consist of an ODE system relating acceleration of the center of gravity and external and pressure loads. The numerical experiments of mechanical stability try to clarify the position of the neutral stability curve. Finally, a rotating machine for ship propulsion involving both axial and radial bearings operating with nonconventional lubricants (seawater to avoid environmental pollution) is analyzed by using laminar and turbulent inertial flows.

Understanding and improving the chemical vapor deposition process for solar grade silicon production

Esta Tesis Doctoral se centra en la investigación del proceso de producción de polisilicio para aplicaciones fotovoltaicas (FV) por la vía química; mediante procesos de depósito en fase vapor (CVD). El polisilicio para la industria FV recibe el nombre de silicio de grado solar (SoG Si). Por un lado, el proceso que domina hoy en día la producción de SoG Si está basado en la síntesis, destilación y descomposición de triclorosilano (TCS) en un reactor CVD -denominado reactor Siemens-. El material obtenido mediante este proceso es de muy alta pureza, pero a costa de un elevado consumo energético. Así, para alcanzar los dos principales objetivos de la industria FV basada en silicio, bajos costes de producción y bajo tiempo de retorno de la energía invertida en su fabricación, es esencial disminuir el consumo energético de los reactores Siemens. Por otro lado, una alternativa al proceso Siemens considera la descomposición de monosilano (MS) en un reactor de lecho fluidizado (FBR). Este proceso alternativo tiene un consumo energético mucho menor que el de un reactor Siemens, si bien la calidad del material resultante es también menor; pero ésta puede ser suficiente para la industria FV. A día de hoy los FBR deben aún abordar una serie de retos para que su menor consumo energético sea una ventaja suficiente comparada con otras desventajas de estos reactores. En resumen, la investigación desarrollada se centra en el proceso de depósito de polysilicio por CVD a partir de TCS -reactor Siemens-; pero también se investiga el proceso de producción de SoG Si en los FBR exponiendo las fortalezas y debilidades de esta alternativa. Para poder profundizar en el conocimiento del proceso CVD para la producción de polisilicio es clave el conocimiento de las reacciones químicas fundamentales y cómo éstas influencian la calidad del producto resultante, al mismo tiempo que comprender los fenómenos responsables del consumo energético. Por medio de un reactor Siemens de laboratorio en el que se llevan a cabo un elevado número de experimentos de depósito de polisilicio de forma satisfactoria se adquiere el conocimiento previamente descrito. Se pone de manifiesto la complejidad de los reactores CVD y de los problemas asociados a la pérdidas de calor de estos procesos. Se identifican las contribuciones a las pérdidas de calor de los reactores CVD, éstas pérdidas de calor son debidas principalmente a los fenómenos de radiación y, conducción y convección vía gases. En el caso de los reactores Siemens el fenómeno que contribuye en mayor medida al alto consumo energético son las pérdidas de calor por radiación, mientras que en los FBRs tanto la radiación como el calor transferido por transporte másico contribuyen de forma importante. Se desarrolla un modelo teórico integral para el cálculo de las pérdidas de calor en reactores Siemens. Este modelo está formado a su vez por un modelo para la evaluación de las pérdidas de calor por radiación y modelos para la evaluación de las pérdidas de calor por conducción y convección vía gases. Se ponen de manifiesto una serie de limitaciones del modelo de pérdidas de calor por radiación, y se desarrollan una serie de modificaciones que mejoran el modelo previo. El modelo integral se valida por medio un reactor Siemens de laboratorio, y una vez validado se presenta su extrapolación a la escala industrial. El proceso de conversión de TCS y MS a polisilicio se investiga mediante modelos de fluidodinámica computacional (CFD). Se desarrollan modelados CFD para un reactor Siemens de laboratorio y para un prototipo FBR. Los resultados obtenidos mediante simulación son comparados, en ambos casos, con resultados experimentales. Los modelos desarrollados se convierten en herramientas para la identificación de aquellos parámetros que tienen mayor influencia en los procesos CVD. En el caso del reactor Siemens, ambos modelos -el modelo integral y el modelado CFD permiten el estudio de los parámetros que afectan en mayor medida al elevado consumo energético, y mediante su análisis se sugieren modificaciones para este tipo de reactores que se traducirían en un menor número de kilovatios-hora consumidos por kilogramo de silicio producido. Para el caso del FBR, el modelado CFD permite analizar el efecto de una serie de parámetros sobre la distribución de temperaturas en el lecho fluidizado; y dicha distribución de temperaturas está directamente relacionada con los principales retos de este tipo de reactores. Por último, existen nuevos conceptos de depósito de polisilicio; éstos se aprovechan de la ventaja teórica de un mayor volumen depositado por unidad de tiempo -cuando una mayor superficie de depósito está disponible- con el objetivo de reducir la energía consumida por los reactores Siemens. Estos conceptos se exploran mediante cálculos teóricos y pruebas en el reactor Siemens de laboratorio. ABSTRACT This Doctoral Thesis comprises research on polysilicon production for photovoltaic (PV) applications through the chemical route: chemical vapor deposition (CVD) process. PV polysilicon is named solar grade silicon (SoG Si). On the one hand, the besetting CVD process for SoG Si production is based on the synthesis, distillation, and decomposition of thriclorosilane (TCS) in the so called Siemens reactor; high purity silicon is obtained at the expense of high energy consumption. Thus, lowering the energy consumption of the Siemens process is essential to achieve the two wider objectives for silicon-based PV technology: low production cost and low energy payback time. On the other hand, a valuable variation of this process considers the use of monosilane (MS) in a fluidized bed reactor (FBR); lower output material quality is obtained but it may fulfil the requirements for the PV industry. FBRs demand lower energy consumption than Siemens reactors but further research is necessary to address the actual challenges of these reactors. In short, this work is centered in polysilicon CVD process from TCS -Siemens reactor-; but it also offers insights on the strengths and weaknesses of the FBR for SoG Si production. In order to aid further development in polysilicon CVD is key the understanding of the fundamental reactions and how they influence the product quality, at the same time as to comprehend the phenomena responsible for the energy consumption. Experiments conducted in a laboratory Siemens reactor prove the satisfactory operation of the prototype reactor, and allow to acquire the knowledge that has been described. Complexity of the CVD reactors is stated and the heat loss problem associated with polysilicon CVD is addressed. All contributions to the energy consumption of Siemens reactors and FBRs are put forward; these phenomena are radiation and, conduction and convection via gases heat loss. In a Siemens reactor the major contributor to the energy consumption is radiation heat loss; in case of FBRs radiation and heat transfer due to mass transport are both important contributors. Theoretical models for radiation, conduction and convection heat loss in a Siemens reactor are developed; shaping a comprehensive theoretical model for heat loss in Siemens reactors. Limitations of the radiation heat loss model are put forward, and a novel contribution to the existing model is developed. The comprehensive model for heat loss is validated through a laboratory Siemens reactor, and results are scaled to industrial reactors. The process of conversion of TCS and MS gases to solid polysilicon is investigated by means of computational fluid-dynamics models. CFD models for a laboratory Siemens reactor and a FBR prototype are developed. Simulated results for both CVD prototypes are compared with experimental data. The developed models are used as a tool to investigate the parameters that more strongly influence both processes. For the Siemens reactors, both, the comprehensive theoretical model and the CFD model allow to identify the parameters responsible for the great power consumption, and thus, suggest some modifications that could decrease the ratio kilowatts-hour per kilogram of silicon produced. For the FBR, the CFD model allows to explore the effect of a number of parameters on the thermal distribution of the fluidized bed; that is the main actual challenge of these type of reactors. Finally, there exist new deposition surface concepts that take advantage of higher volume deposited per time unit -when higher deposition area is available- trying to reduce the high energy consumption of the Siemens reactors. These novel concepts are explored by means of theoretical calculations and tests in the laboratory Siemens prototype.

Transient thermal analysis during the ascent phase of a balloon-borne payload. Comparison with SUNRISE test flight measurements

The thermal design of stratospheric balloon payloads usually focuses on the cruise phase of the missions, that is, the floating altitude conditions. The ascent phase usually takes between 2 and 4 h, a very small period compared to the duration of the whole mission, which can last up to 4 weeks. However, during this phase payloads are subjected to very harsh conditions due mainly to the convective cooling that occurs as the balloon passes through the cold atmosphere, with minimum temperatures in the tropopause. The aim of this work is to study the thermal behaviour of a payload carried by a long duration balloon during the ascent phase. Its temperature has been calculated as a function of the altitude from sea level to floating conditions. To perform this analysis it has been assumed that the thermal interactions (convection and radiation) depend on the altitude, on the environmental conditions (which in turn depend also on the altitude) and on the temperature of the system itself. The results have been compared with the measurements taken during the SUNRISE test flight, launched in October 2007 by CSBF from Fort Sumner (New Mexico).

A numerical strategy for efficient modeling of the equatorial wave guide

Convection in the tropics is observed to involve a wide-ranging hierarchy of scales from a few kilometers to the planetary scales and also has a profound impact on short-term climate. The mechanisms responsible for this behavior present a major unsolved problem. A promising emerging approach to address these issues is cloud-resolving modeling. Here a family of numerical models is introduced specifically to model the feedback of small-scale deep convection on tropical planetary waves and tropical circulation in a highly efficient manner compatible with the approach through cloud-resolving modeling. Such a procedure is also useful for theoretical purposes. The basic idea in the approach is to use low-order truncation in the meriodonal direction through Gauss–Hermite quadrature projected onto a simple discrete radiation condition. In this fashion, the cloud-resolving modeling of equatorially trapped planetary waves reduces to the solution of a small number of purely zonal two-dimensional wave systems along a few judiciously chosen meriodonal layers that are coupled only by some additional source terms. The approach is analyzed in detail with full mathematical rigor for linearized equatorial primitive equations with source terms.

The metabolic implications of intracellular circulation

Two views currently dominate research into cell function and regulation. Model I assumes that cell behavior is quite similar to that expected for a watery bag of enzymes and ligands. Model II assumes that three-dimensional order and structure constrain and determine metabolite behavior. A major problem in cell metabolism is determining why essentially all metabolite concentrations are remarkably stable (are homeostatic) over large changes in pathway fluxes—for convenience, this is termed the [s] stability paradox. For muscle cells, ATP and O2 are the most perfectly homeostatic, even though O2 delivery and metabolic rate correlate in a 1:1 fashion. In total, more than 60 metabolites are known to be remarkably homeostatic in differing metabolic states. Several explanations of [s] stability are usually given by traditional model I studies—none of which apply to all enzymes in a pathway, and all of which require diffusion as the means for changing enzyme–substrate encounter rates. In contrast, recent developments in our understanding of intracellular myosin, kinesin, and dyenin motors running on actin and tubulin tracks or cables supply a mechanistic basis for regulated intracellular circulation systems with cytoplasmic streaming rates varying over an approximately 80-fold range (from 1 to >80 μm × sec−1). These new studies raise a model II hypothesis of intracellular perfusion or convection as a primary means for bringing enzymes and substrates together under variable metabolic conditions. In this view, change in intracellular perfusion rates cause change in enzyme–substrate encounter rates and thus change in pathway fluxes with no requirement for large simultaneous changes in substrate concentrations. The ease with which this hypothesis explains the [s] stability paradox is one of its most compelling features.

Order and disorder in fluid motion.

The development of complex states of fluid motion is illustrated by reviewing a series of experiments, emphasizing film flows, surface waves, and thermal convection. In one dimension, cellular patterns bifurcate to states of spatiotemporal chaos. In two dimensions, even ordered patterns can be surprisingly intricate when quasiperiodic patterns are included. Spatiotemporal chaos is best characterized statistically, and methods for doing so are evolving. Transport and mixing phenomena can also lead to spatial complexity, but the degree depends on the significance of molecular or thermal diffusion.

Variability of the ocean heat content during the last millennium - an assessment with the ECHO-g Model

Studies addressing climate variability during the last millennium generally focus on variables with a direct influence on climate variability, like the fast thermal response to varying radiative forcing, or the large-scale changes in atmospheric dynamics (e. g. North Atlantic Oscillation). The ocean responds to these variations by slowly integrating in depth the upper heat flux changes, thus producing a delayed influence on ocean heat content (OHC) that can later impact low frequency SST (sea surface temperature) variability through reemergence processes. In this study, both the externally and internally driven variations of the OHC during the last millennium are investigated using a set of fully coupled simulations with the ECHO-G (coupled climate model ECHAMA4 and ocean model HOPE-G) atmosphere-ocean general circulation model (AOGCM). When compared to observations for the last 55 yr, the model tends to overestimate the global trends and underestimate the decadal OHC variability. Extending the analysis back to the last one thousand years, the main impact of the radiative forcing is an OHC increase at high latitudes, explained to some extent by a reduction in cloud cover and the subsequent increase of short-wave radiation at the surface. This OHC response is dominated by the effect of volcanism in the preindustrial era, and by the fast increase of GHGs during the last 150 yr. Likewise, salient impacts from internal climate variability are observed at regional scales. For instance, upper temperature in the equatorial Pacific is controlled by ENSO (El Nino Southern Oscillation) variability from interannual to multidecadal timescales. Also, both the Pacific Decadal Oscillation (PDO) and the Atlantic Multidecadal Oscillation (AMO) modulate intermittently the interdecadal OHC variability in the North Pacific and Mid Atlantic, respectively. The NAO, through its influence on North Atlantic surface heat fluxes and convection, also plays an important role on the OHC at multiple timescales, leading first to a cooling in the Labrador and Irminger seas, and later on to a North Atlantic warming, associated with a delayed impact on the AMO.

Reconstituição da paleoprecipitação no sul do Nordeste Brasileiro durante os dois últimos ciclos glaciais a partir da aplicação de registros isotópicos de oxigênio de estalagmites da Chapada Diamantina, Bahia.

Registros isotópicos de oxigênio obtidos em alta resolução das estalagmites CL2 e MAG das cavernas Calixto e Marota, região da Chapada Diamantina (CD) (12ºS), Estado da Bahia, sul do Nordeste brasileiro (sNEB), permitiram reconstituir as mudanças passadas da precipitação entre 165-128 e 59-39 mil anos A.P. Para a reconstituição paleoclimática considerou-se resultados de um estudo de calibração realizado em duas cavernas da CD o qual demonstrou uma relação entre composição isotópica da água meteórica e de gotejamento e sugeriu um ambiente adequado para a deposição do espeleotema em condições equilíbrio e/ou próximas com a água de gotejamento. A interpretação da paleoprecipitação através dos registros isotópicos \'\'delta\' POT.18\'O das estalagmites também foi baseada na relação entre composição isotópica da água da precipitação e a quantidade de chuva obtidos em estações da IAEA-GNIP no Brasil e de simulações das variações do \'\'delta\' POT.18\'O da chuva através do modelo climático ECHAM-4. Esses dados indicaram o efeito quantidade (amount effect) como fator preponderante de controle isotópico da água da chuva que formam os espeleotemas na CD, significando que a diminuição dos valores de \'\'delta\' POT.18\'O está associada ao aumento do volume de chuvas e vice-versa. Os registros de \'\'delta\' POT.18\'O dos espeleotemas permitiram reconstituir a variação da paleoprecipitação na escala orbital e milenar durante o penúltimo glacial bem como correlacionar mudanças na paleoprecipitação no sNEB com eventos milenares registrados na Groelândia no último glacial. Os registros da CD indicaram um aumento (diminuição) da paleoprecipitação na Bahia relacionado a diminuição (aumento) da insolação austral de verão a 10ºS durante o penúltimo glacial, similar ao observado no último ciclo precessional. Na escala orbital os registros da CD estiveram em antifase com os paleoindicadores isotópicos do Sudeste brasileiro e em fase com os valores de\'\'delta\' POT.18\'O dos espeleotemas do leste da China. Esse padrão de precipitação é similar ao observado na última glaciação e sugere que a variação na insolação de verão afetou as monções sul-americanas (MSA) promovendo mudanças na precipitação no sNEB no penúltimo glacial. Condições áridas no sNEB durante o aumento da insolação de verão estariam provavelmente associadas ao aprofundamento da subsidência de ar provocado pelo fortalecimento da circulação leste-oeste da MSA devido ao aumento das atividades convectivas na Amazônia o que teria, favorecido um posicionamento mais a sul da Zona de Convergência do Atlântico Sul (ZCAS). O oposto também ocorreria durante as fases de baixa insolação de verão quando a MSA estaria provavelmente mais desintensificada. Durante o penúltimo glacial (Terminação Glacial II) abruptas oscilações nos registros da CD para valores mais baixos de \'\'delta\' POT.18\'O indicaram um profundo aumento da precipitação coincidente com o evento Heinrich (H11). Nesse período a paleoprecipitação no sNEB esteve correlacionada negativamente com as mudanças climáticas na China e no oeste amazônico (Peru) e positivamente com o Sudeste brasileiro. Interpretou-se que as anomalias positivas da precipitação no sNEB podem ter estado relacionadas ao deslocamento para sul da Zona de Convergência Intertropical (ZCIT) bem como com a intensificação da MSA e ZCAS nesse período. Finalmente, oscilações isotópicas abruptas para valores mais altos observadas durante o estágio marinho isotópico 3 coincidentes com os eventos quentes registrados na Groelândia, denominados de eventos Dansgaard-Oeschger (DO), foram interpretados como a ocorrência de eventos muito secos no sNEB. Essas variações da precipitação na escala milenar, que estão em fase com os registros no Peru, podem ter estado relacionadas ao deslocamento para norte da ZCIT o que teria promovido uma profunda desintensificação da MSA.

Stable isotope record, Mg/Ca ratios and sea surface reconstruction for sediment core MD06-3067

We present sea surface, upper thermocline, and benthic d18O data, as well as temperature and paleoproductivity proxy data, from the International Marine Global Change Study Program (IMAGES) Core MD06-3067 (6°31'N, 126°30'E, 1575 m water depth), located in the western equatorial Pacific Ocean within the flow path of the Mindanao Current. Our records reveal considerable glacial-interglacial and suborbital variability in the Mindanao Dome upwelling over the last 160 kyr. Dome activity generally intensified during glacial intervals resulting in cooler thermocline waters, whereas it substantially declined during interglacials, in particular in the early Holocene and early marine oxygen isotope stage (MIS) 5e, when upwelling waters did not reach the thermocline. During MIS 3 and MIS 2, enhanced surface productivity together with remarkably low SST and low upper ocean thermal contrast provide evidence for episodic glacial upwelling to the surface, whereas transient surface warming marks periodic collapses of the Mindanao Dome upwelling during Heinrich events. We attribute the high variability during MIS 3 and MIS 2 to changes in the El Niño Southern Oscillation state that affected boreal winter monsoonal winds and upper ocean circulation. Glacial upwelling intensified when a strong cyclonic gyre became established, whereas El Niño-like conditions during Heinrich events tended to suppress the cyclonic circulation, reducing Ekman transport. Thus, our findings demonstrate that variations in the Mindanao Dome upwelling are closely linked to the position and intensity of the tropical convection and also reflect far-field influences from the high latitudes.

Data compilation of ESOP

The work in this sub-project of ESOP focuses on the advective and convective transforma-tion of water masses in the Greenland Sea and its neighbouring areas. It includes observational work on the sub-mesoscale and analysis of hydrographic data up to the gyre-scale. Observations of active convective plumes were made with a towed chain equipped with up to 80 CTD sensors, giving a horizontal and vertical resolution of the hydrographic fields of a few metres. The observed scales of the penetrative convective plumes compare well with those given by theory. On the mesoscale the structure of homogeneous eddies formed as a result of deep convection was observed and the associated mixing and renewal of the intermediate layers quantified. The relative importance and efficiency of thermal and haline penetrative convection in relation to the surface boundary conditions (heat and salt fluxes and ice cover) and the ambient stratification are studied using the multi year time series of hydro-graphic data in the central Greenland Sea. The modification of the water column of the Greenland Sea gyre through advection from and mixing with water at its rim is assessed on longer time scales. The relative contributions are quantified using modern water mass analysis methods based on inverse techniques. Likewise the convective renewal and the spreading of the Arctic Intermediate Water from its formation area is quantified. The aim is to budget the heat and salt content of the water column, in particular of the low salinity surface layer, and to relate its seasonal and interannual variability to the lateral fluxes and the fluxes at the air-sea-ice interface. This will allow to estimate residence times for the different layers of the Greenland Sea gyre, a quantity important for the description of the Polar Ocean carbon cycle.

Nd-isotope ratios of Cretaceous sediment samples

The role that meridional overturning circulation (MOC) patterns played in poleward heat transport during the extreme warmth of the Early to Late Cretaceous is a fundamental and unresolved question in climate dynamics. In order to address this question we must determine where deep waters formed, and how they may have circulated during periods of extreme warmth. Here we present late Albian through Maastrichtian (105 to 65 Ma) Nd isotope records from Deep Sea Drilling Project (DSDP) and Ocean Drilling Program (ODP) sites in the proto-Indian Ocean and the tropical Pacific. Comparison of these data with previously published records indicates deep-water formation in the Indian sector of the Southern Ocean began at least ?105 Ma, extending the record of high-latitude convection back into the Early Cretaceous prior to the peak warmth of the mid-Cretaceous. The growing body of data supports a mode of MOC in part characterized by high-latitude downwelling during the peak of greenhouse warmth of the Mesozoic and Cenozoic. However, this mode of MOC likely was characterized by numerous locations of deep convection that were regionally important, but not significant in terms of a globally overturning circulation due to paleogeographic and bathymetric barriers.

(Table 1) Annual melt rates of the South Branch of Jacobshavn Isbrae between July 1984-1985

Photogrammetric reanalysis of 1985 aerial photos has revealed substantial submarine melting of the floating ice tongue of Jakobshavn Isbrae, west Greenland. The thickness of the floating tongue determined from hydrostatic equilibrium tapers from ~940 m near the grounding zone to ~600 m near the terminus. Feature tracking on orthophotos shows speeds on the July 1985 ice tongue to be nearly constant (~18.5 m/d), indicating negligible dynamic thinning. The thinning of the ice tongue is mostly due to submarine melting with average rates of 228 ± 49 m/yr (0.62 ± 0.13 m/d) between the summers of 1984 and 1985. The cause of the high melt rate is the circulation of warm seawater (thermal forcing of up to 4.2°C) beneath the tongue with convection driven by the substantial discharge of subglacial freshwater from the grounding zone. We believe that this buoyancy-driven convection is responsible for a deep channel incised into the sole of the floating tongue. A dramatic thinning, retreat, and speedup began in 1998 and continues today. The timing of the change is coincident with a 1.1°C warming of deep ocean waters entering the fjord after 1997. Assuming a linear relationship between thermal forcing and submarine melt rate, average melt rates should have increased by ~25% (~57 m/yr), sufficient to destabilize the ice tongue and initiate the ice thinning and the retreat that followed.

Zooplankton biomass from METEOR cruise M87/1 in April 2012

Zooplankton samples were taken in five depth strata using a Multinet type Midi, with 50 µm nets. The samples were taken during the second leg only, three times at station 1, two times at station 2 and once at station 3. Zooplankton were identified to species / genus and life-stage, and at least 300 individuals were counted per sample. 10 individuals of each stage / species were measured and the numbers of eggs counted.