Processes controlling surface, bottom and lateral melt of Arctic sea ice in a state of the art sea ice model

We present a modelling study of processes controlling the summer melt of the Arctic sea ice cover. We perform a sensitivity study and focus our interest on the thermodynamics at the ice–atmosphere and ice–ocean interfaces. We use the Los Alamos community sea ice model CICE, and additionally implement and test three new parametrization schemes: (i) a prognostic mixed layer; (ii) a three equation boundary condition for the salt and heat flux at the ice–ocean interface; and (iii) a new lateral melt parametrization. Recent additions to the CICE model are also tested, including explicit melt ponds, a form drag parametrization and a halodynamic brine drainage scheme. The various sea ice parametrizations tested in this sensitivity study introduce a wide spread in the simulated sea ice characteristics. For each simulation, the total melt is decomposed into its surface, bottom and lateral melt components to assess the processes driving melt and how this varies regionally and temporally. Because this study quantifies the relative importance of several processes in driving the summer melt of sea ice, this work can serve as a guide for future research priorities.

Quantification of the probability of containment failure caused by an in-vessel steam explosion for the Sizewell B PWR

AEA Technology has provided an assessment of the probability of α-mode containment failure for the Sizewell B PWR. After a preliminary review of the methodologies available it was decided to use the probabilistic approach described in the paper, based on an extension of the methodology developed by Theofanous et al. (Nucl. Sci. Eng. 97 (1987) 259–325). The input to the assessment is 12 probability distributions; the bases for the quantification of these distributions are discussed. The α-mode assessment performed for the Sizewell B PWR has demonstrated the practicality of the event-tree method with input data represented by probability distributions. The assessment itself has drawn attention to a number of topics, which may be plant and sequence dependent, and has indicated the importance of melt relocation scenarios. The α-mode failure probability following an accident that leads to core melt relocation to the lower head for the Sizewell B PWR has been assessed as a few parts in 10 000, on the basis of current information. This assessment has been the first to consider elevated pressures (6 MPa and 15 MPa) besides atmospheric pressure, but the results suggest only a modest sensitivity to system pressure.

The influence of viscous and latent heating on crystal-rich magma flow in a conduit

The flow dynamics of crystal-rich high-viscosity magma is likely to be strongly influenced by viscous and latent heat release. Viscous heating is observed to play an important role in the dynamics of fluids with temperature-dependent viscosities. The growth of microlite crystals and the accompanying release of latent heat should play a similar role in raising fluid temperatures. Earlier models of viscous heating in magmas have shown the potential for unstable (thermal runaway) flow as described by a Gruntfest number, using an Arrhenius temperature dependence for the viscosity, but have not considered crystal growth or latent heating. We present a theoretical model for magma flow in an axisymmetric conduit and consider both heating effects using Finite Element Method techniques. We consider a constant mass flux in a 1-D infinitesimal conduit segment with isothermal and adiabatic boundary conditions and Newtonian and non-Newtonian magma flow properties. We find that the growth of crystals acts to stabilize the flow field and make the magma less likely to experience a thermal runaway. The additional heating influences crystal growth and can counteract supercooling from degassing-induced crystallization and drive the residual melt composition back towards the liquidus temperature. We illustrate the models with results generated using parameters appropriate for the andesite lava dome-forming eruption at Soufriere Hills Volcano, Montserrat. These results emphasize the radial variability of the magma. Both viscous and latent heating effects are shown to be capable of playing a significant role in the eruption dynamics of Soufriere Hills Volcano. Latent heating is a factor in the top two kilometres of the conduit and may be responsible for relatively short-term (days) transients. Viscous heating is less restricted spatially, but because thermal runaway requires periods of hundreds of days to be achieved, the process is likely to be interrupted. Our models show that thermal evolution of the conduit walls could lead to an increase in the effective diameter of flow and an increase in flux at constant magma pressure.

Influence of molecular composition on the development of microstructure from sheared polyethylene melts: Molecular and lamellar templating

A combination of in situ and ex situ X-ray scattering techniques and transmission electron microscopy has been used to study the crystallization behaviour of polyethylene, following the imposition of melt shear. In the case of a branched material, the imposition of shear flow up to a rate of 30 s(-1) was found to induce no anisotropy. Although shearing the linear material only ever induced a very small degree of anisotropy in the melt, for shear rates > 0.15 s(-1), subsequent crystallization resulted in increasing anisotropy. Blends of the above two polyethylenes were produced, in which the linear material constituted the minority fraction (similar to 10%). Isothermal crystallization at temperatures where extensive crystallization of the branched material does not occur demonstrated that the behaviour of the linear component of the sheared blend mirrored that of the linear polyethylene alone. However, in addition, it was found that when crystallized in the presence of an oriented morphology, the branched polymer also formed anisotropic structures. We have termed the process templating, in which the crystallization behaviour of the bulk of the system (similar to 90% branched material) is completely altered (spherulitic to oriented lamellar) by mapping it onto a pre-existing minority structure (similar to 10% linear polymer). (c) 2006 Elsevier Ltd. All rights reserved.

Linear viscoelasticity from molecular dynamics simulation of entangled polymers

The linear viscoelastic (LVE) spectrum is one of the primary fingerprints of polymer solutions and melts, carrying information about most relaxation processes in the system. Many single chain theories and models start with predicting the LVE spectrum to validate their assumptions. However, until now, no reliable linear stress relaxation data were available from simulations of multichain systems. In this work, we propose a new efficient way to calculate a wide variety of correlation functions and mean-square displacements during simulations without significant additional CPU cost. Using this method, we calculate stress−stress autocorrelation functions for a simple bead−spring model of polymer melt for a wide range of chain lengths, densities, temperatures, and chain stiffnesses. The obtained stress−stress autocorrelation functions were compared with the single chain slip−spring model in order to obtain entanglement related parameters, such as the plateau modulus or the molecular weight between entanglements. Then, the dependence of the plateau modulus on the packing length is discussed. We have also identified three different contributions to the stress relaxation:  bond length relaxation, colloidal and polymeric. Their dependence on the density and the temperature is demonstrated for short unentangled systems without inertia.

Enhanced Templating in the Crystallisation of Poly(E-caprolactone) Using 1,3:2,4-di(4-Chlorobenzylidene) Sorbitol

We show that small quantities of 1,3:2,4-di(4-chlorobenzylidene) sorbitol dispersed in poly(epsilon-caprolactone) provide a very effective self-assembling nanoscale framework which, with a flow field, yields extremely high levels of polymer crystal orientation. During modest shear flow of the polymer melt, the additive forms highly extended nano-particles which adopt a preferred alignment with respect to the flow field. On cooling, polymer crystallisation is directed by these particles. This chloro substituted dibenzylidene sorbitol is considerably more effective at directing the crystal growth of poly(epsilon-caprolactone) than the unsubstituted compound.

An Integrated Assessment of changes in the thermohaline circulation

This paper discusses the risks of a shutdown of the thermohaline circulation (THC) for the climate system, for ecosystems in and around the North Atlantic as well as for fisheries and agriculture by way of an Integrated Assessment. The climate model simulations are based on greenhouse gas scenarios for the 21st century and beyond. A shutdown of the THC, complete by 2150, is triggered if increased freshwater input from inland ice melt or enhanced runoff is assumed. The shutdown retards the greenhouse gas-induced atmospheric warming trend in the Northern Hemisphere, but does not lead to a persistent net cooling. Due to the simulated THC shutdown the sea level at the North Atlantic shores rises by up to 80 cm by 2150, in addition to the global sea level rise. This could potentially be a serious impact that requires expensive coastal protection measures. A reduction of marine net primary productivity is associated with the impacts of warming rather than a THC shutdown. Regional shifts in the currents in the Nordic Seas could strongly deteriorate survival chances for cod larvae and juveniles. This could lead to cod fisheries becoming unprofitable by the end of the 21st century. While regional socioeconomic impacts might be large, damages would be probably small in relation to the respective gross national products. Terrestrial ecosystem productivity is affected much more by the fertilization from the increasing CO2 concentration than by a THC shutdown. In addition, the level of warming in the 22nd to 24th century favours crop production in northern Europe a lot, no matter whether the THC shuts down or not. CO2 emissions corridors aimed at limiting the risk of a THC breakdown to 10% or less are narrow, requiring departure from business-as-usual in the next few decades. The uncertainty about THC risks is still high. This is seen in model analyses as well as in the experts’ views that were elicited. The overview of results presented here is the outcome of the Integrated Assessment project INTEGRATION.

Pb-210-Ra-226 and Ra-228-Th-232 systematics in young arc lavas: implications for magma degassing and ascent rates

New data show that island arc rocks have (Pb-210/Ra-226)(o) ratios which range from as low as 0.24 up to 2.88. In contrast, (Ra-22S/Th-232) appears always within error of I suggesting that the large Ra-226-excesses observed in arc rocks were generated more than 30 years ago. This places a maximum estimate on melt ascent velocities of around 4000 m/year and provides further confidence that the Ra-226 excesses reflect deep (source) processes rather than shallow level alteration or seawater contamination. Conversely, partial melting must have occurred more than 30 years prior to eruption. The Pb-210 deficits are most readily explained by protracted magma degassing. Using published numerical models, the data suggest that degassing occurred continuously for periods up to several decades just prior to eruption but no link with eruption periodicity was found. Longer periods are required if degassing is discontinuous, less than 100% efficient or if magma is recharged or stored after degassing. The long durations suggest much of this degassing occurs at depth with implications for the formation of hydrothermal and copper-porphyry systems. A suite of lavas erupted in 1985-1986 from Sangeang Api volcano in the Sunda arc are characterised by deficits of Pb-210 relative to Ra-226 from which 6-8 years of continuous Rn-222 degassing would be inferred from recent numerical models. These data also form a linear (Pb-210)/Pb-(Ra-226)/Pb array which might be interpreted as a 71-year isochron. However, the array passes through the origin suggesting displacement downwards from the equiline in response to degassing and so the slope of the array is inferred not to have any age significance. Simple modelling shows that the range of (Ra-226)/Pb ratios requires thousands of years to develop consistent with differentiation occurring in response to cooling at the base of the crust. Thus, degassing post-dated, and was not responsible for magma differentiation. The formation, migration and extraction of gas bubbles must be extremely efficient in mafic magma whereas the higher viscosity of more siliceous magmas retards the process and can lead to Pb-210 excesses. A possible negative correlation between (Pb-210/Ra-226)(o) and SO2 emission rate requires further testing but may have implications for future eruptions. (C) 2004 Elsevier B.V. All rights reserved.

Seasonal trends in the stable isotopic composition of snow and meltwater runoff in a subarctic catchment at Okstindan, Norway

Seasonal variations in the stable isotopic composition of snow and meltwater were investigated in a sub-arctic, mountainous, but non-glacial, catchment at Okstindan in northern Norway based on analyses of delta(18)O and deltaD. Samples were collected during four field periods (August 1998; April 1999; June 1999 and August 1999) at three sites lying on an altitudinal transect (740-970 m a.s.l.). Snowpack data display an increase in the mean values of delta(18)O (increasing from a mean value of - 13.51 to - 11.49% between April and August), as well as a decrease in variability through the melt period. Comparison with a regional meteoric water line indicates that the slope of the delta(18)O - deltaD line for the snowpacks decreases over the same period, dropping from 7.49 to approximately 6.2. This change points to the role of evaporation in snowpack ablation and is confirmed by the vertical profile of deuterium excess. Snowpack seepage data, although limited, also suggest reduced values of deltaD, as might be associated with local evaporation during meltwater generation. In general, meltwaters were depleted in delta(18)O relative to the source snowpack at the peak of the melt (June), but later in the year (August) the difference between the two was not statistically significant. The diurnal pattern of isotopic composition indicates that the most depleted meltwaters coincide with the peak in temperature and, hence, meltwater production.

Climate change, glacier retreat, and water availability in the Caucasus region

The paper discusses the observed and projected warming in the Caucasus region and its implications for glacier melt, water availability and potential hazards. A strong positive trend in summer air temperatures of 0.05 degrees C year(-1) is observed in the high-altitude areas (above 2000 m) providing for a strong glacier melt. A widespread glacier retreat has also been reported between 1985 and 2000, with an average rate of 8 m year(-1). A warming of 5-7 degrees C is projected for the Sum mer months in the 2071-2100 period under the A2 emission group of scenarios, Suggesting that enhanced glacier melt and a changing water balance can be expected.

An Assessment of the Recent Past and Future Climate Change, Glacier Retreat, and Runoff in the Caucasus Region Using Dynamical and Statistical Downscaling and HBV-ETH Hydrological Model

The paper discusses the observed and projected warming in the Caucasus region and its implications for glacier melt and runoff. A strong positive trend in summer air temperatures of 0.05 degrees C a(-1) is observed in the high-altitude areas providing for a strong glacier melt and continuous decline in glacier mass balance. A warming of 4-7 degrees C and 3-5 degrees C is projected for the summer months in 2071-2100 under the A2 and B2 emission scenarios respectively, suggesting that enhanced glacier melt can be expected. The expected changes in winter precipitation will not compensate for the summer melt and glacier retreat is likely to continue. However, a projected small increase in both winter and summer precipitation combined with the enhanced glacier melt will result in increased summer runoff in the currently glaciated region of the Caucasus (independent of whether the region is glaciated at the end of the twenty-first century) by more than 50% compared with the baseline period.

Long-term change, interannual and intra-seasonal variability in climate and glacier mass balance in the Central Greater Caucasus, Russia

Long-term trends, interannual and intra-seasonal variability in the mass-balance record from Djankuat glacier, central Greater Caucasus, Russia, are related to local climate change, synoptic and large-scale anomalies in atmospheric circulation. A clear warming signal emerged in the central Greater Caucasus in the early 1990s, leading to a strong increase in ablation. In the absence of a compensating change in winter accumulation, the net mass balance of Djankuat has declined. The highest value of seasonal ablation on record was registered in the summer of 2000. At the beginning of the 21st century these trends reversed. Ablation was below average even in the summer of 2003, which was unusually warm in western Europe. Precipitation and winter accumulation were high, allowing for a partial recovery of net mass balance. The interannual variability in the components of mass balance is weakly related to the North Atlantic Oscillation (NAO) and the Scandinavian teleconnection patterns, but there is a clear link with the large-scale circulation anomalies represented by the Rossby pattern. Five synoptic categories have been identified for the ablation season of 2005, revealing a strong separation between components of radiation budget, air temperature and daily melt. Air temperature is the main control over melt. The highest values of daily ablation are related to the strongly positive NAO which forces high net radiation, and to the warm and moist advection from the Black Sea.

Interactions between mass balance, atmospheric circulation, and recent climate change on the Djankuat Glacier, Caucasus Mountains, Russia

This paper reports recent changes in the mass balance record from the Djankuat Glacier, central greater Caucasus, Russia, and investigates possible relationships between the components of mass balance, local climate, and distant atmospheric forcing. The results clearly show that a strong warming signal has emerged in the central greater Caucasus, particularly since the 1993/1994 mass balance year, and this has led to a significant increase in the summer ablation of Djankuat. At the same time, there has been no compensating consistent increase in winter precipitation and accumulation leading to the strong net loss of mass and increase in glacier runoff. Interannual variability in ablation and accumulation is partly associated with certain major patterns of Northern Hemisphere climatic variability. The positive phase of the North Pacific (NP) teleconnection pattern forces negative geopotential height and temperature anomalies over the Caucasus in summer and results in reduced summer melt, such as in the early 1990s, when positive NP extremes resulted in a temporary decline in ablation rates. The positive phase of the NP is related to El Nino-Southern Oscillation, and it is possible that a teleconnection between the tropical Pacific sea surface temperatures and summer air temperatures in the Caucasus is bridged through the NP pattern. More recently, the NP pattern was predominantly negative, and this distant moderating forcing on summer ablation in the Caucasus was absent. Statistically significant correlations are observed between accumulation and the Scandinavian (SCA) teleconnection pattern. The frequent occurrence of the positive SCA phase at the beginning of accumulation season results in lower than average snowfall and reduced accumulation. The relationship between the North Atlantic Oscillation (NAO), Arctic Oscillation, and accumulation is weak, although positive precipitation anomalies in the winter months are associated with the negative phase of the NAO. A stronger positive correlation is observed between accumulation on Djankuat and geopotential height over the Bay of Biscay unrelated to the established modes of the Northern Hemisphere climatic variability. These results imply that the mass balance of Djankuat is sensitive to the natural variability in the climate system. Distant forcing, however, explains only 16% of the variance in the ablation record and cannot fully explain the recent increase in ablation and negative mass balance.

Pb-210-Ra-226 and Ra-228-Th-232 systematics in young arc lavas: implications for magma degassing and ascent rates

New data show that island arc rocks have (Pb-210/Ra-226)(o) ratios which range from as low as 0.24 up to 2.88. In contrast, (Ra-22S/Th-232) appears always within error of I suggesting that the large Ra-226-excesses observed in arc rocks were generated more than 30 years ago. This places a maximum estimate on melt ascent velocities of around 4000 m/year and provides further confidence that the Ra-226 excesses reflect deep (source) processes rather than shallow level alteration or seawater contamination. Conversely, partial melting must have occurred more than 30 years prior to eruption. The Pb-210 deficits are most readily explained by protracted magma degassing. Using published numerical models, the data suggest that degassing occurred continuously for periods up to several decades just prior to eruption but no link with eruption periodicity was found. Longer periods are required if degassing is discontinuous, less than 100% efficient or if magma is recharged or stored after degassing. The long durations suggest much of this degassing occurs at depth with implications for the formation of hydrothermal and copper-porphyry systems. A suite of lavas erupted in 1985-1986 from Sangeang Api volcano in the Sunda arc are characterised by deficits of Pb-210 relative to Ra-226 from which 6-8 years of continuous Rn-222 degassing would be inferred from recent numerical models. These data also form a linear (Pb-210)/Pb-(Ra-226)/Pb array which might be interpreted as a 71-year isochron. However, the array passes through the origin suggesting displacement downwards from the equiline in response to degassing and so the slope of the array is inferred not to have any age significance. Simple modelling shows that the range of (Ra-226)/Pb ratios requires thousands of years to develop consistent with differentiation occurring in response to cooling at the base of the crust. Thus, degassing post-dated, and was not responsible for magma differentiation. The formation, migration and extraction of gas bubbles must be extremely efficient in mafic magma whereas the higher viscosity of more siliceous magmas retards the process and can lead to Pb-210 excesses. A possible negative correlation between (Pb-210/Ra-226)(o) and SO2 emission rate requires further testing but may have implications for future eruptions. (C) 2004 Elsevier B.V. All rights reserved.

Lava production at Soufriere Hills Volcano, Montserrat: 1995-2009

[1] We estimate that about 1 km3 of andesitic lava has been produced at Soufrière Hills Volcano, Montserrat from 1995 to 2009. There were three major episodes of extrusion, each lasting about 2 to 3.5 years and producing about 280 to 340 M m3 of lava, and one minor episode. Our estimates account for the dense rock equivalent volumetric contributions from the core and talus components of the lava dome, pyroclastic flow deposits and air-fall deposits. By 2005 at least two thirds of the erupted mass has already entered the sea. The average lava flux across the major extrusion episodes has been 3–5 m3s−1, with short-period (10–15 days) pulses up to 10–20 m3s−1. The first and third episodes of extrusion show similar flux histories suggesting similar behaviour of the system ten years apart. Waning flux towards the end of each episode may be caused by declining overpressure in the magma reservoir.