Committed terrestrial ecosystem changes due to climate change

Targets for stabilizing climate change are often based on considerations of the impacts of different levels of global warming, usually assessing the time of reaching a particular level of warming. However, some aspects of the Earth system, such as global mean temperatures1 and sea level rise due to thermal expansion2 or the melting of large ice sheets3, continue to respond long after the stabilization of radiative forcing. Here we use a coupled climate–vegetation model to show that in turn the terrestrial biosphere shows significant inertia in its response to climate change. We demonstrate that the global terrestrial biosphere can continue to change for decades after climate stabilization. We suggest that ecosystems can be committed to long-term change long before any response is observable: for example, we find that the risk of significant loss of forest cover in Amazonia rises rapidly for a global mean temperature rise above 2 °C. We conclude that such committed ecosystem changes must be considered in the definition of dangerous climate change, and subsequent policy development to avoid it.

The effects of rotation and ice shelf topography on frazil-laden ice shelf water plumes

A model of the dynamics and thermodynamics of a plume of meltwater at the base of an ice shelf is presented. Such ice shelf water plumes may become supercooled and deposit marine ice if they rise (because of the pressure decrease in the in situ freezing temperature), so the model incorporates both melting and freezing at the ice shelf base and a multiple-size-class model of frazil ice dynamics and deposition. The plume is considered in two horizontal dimensions, so the influence of Coriolis forces is incorporated for the first time. It is found that rotation is extremely influential, with simulated plumes flowing in near-geostrophy because of the low friction at a smooth ice shelf base. As a result, an ice shelf water plume will only rise and become supercooled (and thus deposit marine ice) if it is constrained to flow upslope by topography. This result agrees with the observed distribution of marine ice under Filchner–Ronne Ice Shelf, Antarctica. In addition, it is found that the model only produces reasonable marine ice formation rates when an accurate ice shelf draft is used, implying that the characteristics of real ice shelf water plumes can only be captured using models with both rotation and a realistic topography.

Modelling the rheology of sea ice as a collection of diamond-shaped floes

In polar oceans, seawater freezes to form a layer of sea ice of several metres thickness that can cover up to 8% of the Earth’s surface. The modelled sea ice cover state is described by thickness and orientational distribution of interlocking, anisotropic diamond-shaped ice floes delineated by slip lines, as supported by observation. The purpose of this study is to develop a set of equations describing the mean-field sea ice stresses that result from interactions between the ice floes and the evolution of the ice floe orientation, which are simple enough to be incorporated into a climate model. The sea ice stress caused by a deformation of the ice cover is determined by employing an existing kinematic model of ice floe motion, which enables us to calculate the forces acting on the ice floes due to crushing into and sliding past each other, and then by averaging over all possible floe orientations. We describe the orientational floe distribution with a structure tensor and propose an evolution equation for this tensor that accounts for rigid body rotation of the floes, their apparent re-orientation due to new slip line formation, and change of shape of the floes due to freezing and melting. The form of the evolution equation proposed is motivated by laboratory observations of sea ice failure under controlled conditions. Finally, we present simulations of the evolution of sea ice stress and floe orientation for several imposed flow types. Although evidence to test the simulations against is lacking, the simulations seem physically reasonable.

Modeling the summertime evolution of sea-ice melt ponds

1] We present a mathematical model describing the summer melting of sea ice. We simulate the evolution of melt ponds and determine area coverage and total surface ablation. The model predictions are tested for sensitivity to the melt rate of unponded ice, enhanced melt rate beneath the melt ponds, vertical seepage, and horizontal permeability. The model is initialized with surface topographies derived from laser altimetry corresponding to first-year sea ice and multiyear sea ice. We predict that there are large differences in the depth of melt ponds and the area of coverage between the two types of ice. We also find that the vertical seepage rate and the melt rate of unponded ice are important in determining the total surface ablation and area covered by melt ponds.

Validity of the ice shelf water plume concept beneath Filchner-Ronne Ice Shelf

In winter, brine rejection from sea ice formation and export in the Weddell Sea, offshore of Filchner-Ronne Ice Shelf (FRIS), leads to the formation of High Salinity Shelf Water (HSSW). This dense water mass enters the cavity beneath FRIS by sinking southward down the sloping continental shelf towards the grounding line. Melting occurs when the HSSW encounters the ice shelf, and the meltwater released cools and freshens the HSSW to form a water mass known as Ice Shelf Water (ISW). If this ISW rises, the ‘ice pump’ is initiated (Lewis and Perkin, 1986), whereby the ascending ISW becomes supercooled and deposits marine ice at shallower locations due to the pressure increase in the in-situ freezing temperature. Sandh¨ager et al. (2004) were able to infer the thickness patterns of marine ice deposits at the base of FRIS (figure 1), so the primary aim of this work is to try to understand the ocean flows that determine these patterns. The plume model we use to investigate ISW flow is described fully by Holland and Feltham (accepted) so only a relatively brief outline is presented here. The plume is simulated by combining a parameterisation of ice shelf basal interaction and a multiplesize- class frazil dynamics model with an unsteady, depth-averaged reduced-gravity plume model. In the model an active region of ISW evolves above and within an expanse of stagnant ambient fluid, which is considered to be ice-free and has fixed profiles of temperature and salinity. The two main assumptions of the model are that there is a well-mixed layer underneath the ice shelf and that the ambient fluid outside the plume is stagnant with fixed properties. The topography of the ice shelf that the plume flows beneath is set to the FRIS ice shelf draft calculated by Sandh¨ager et al. (2004) masked with the grounding line from the Antarctic Digital Database (ADD Consortium, 2002). To initiate the plumes, we assume that the intrusion of dense HSSW initially causes melting at the points on the grounding line where the glaciological tributaries feeding FRIS go afloat.

Development of frazil modelling in ice shelf water plumes

Abstract Preliminary results are presented from a modelling study directed at the spatial variation of frazil ice formation and its effects on flow underneath large ice shelves. The chosen plume and frazil models are briefly introduced, and results from two simplified cases are outlined. It is found that growth and melting dominate the frazil model in the short term. Secondary nucleation converts larger crystals into several nuclei due to crystal collisions (microattrition) and fluid shear and therefore governs the ice crystal dynamics after the initial supercooling has been quenched. Frazil formation is found to have a significant depth-dependence in an idealised study of an Ice Shelf Water plume. Finally, plans for more extensive and realistic studies are discussed.

September Arctic sea-ice minimum predicted by spring melt-pond fraction

The area of Arctic September sea ice has diminished from about 7 million km2 in the 1990s to less than 5 million km2 in five of the past seven years, with a record minimum of 3.6 million km2 in 2012 (ref. 1). The strength of this decrease is greater than expected by the scientific community, the reasons for this are not fully understood, and its simulation is an on-going challenge for existing climate models2, 3. With growing Arctic marine activity there is an urgent demand for forecasting Arctic summer sea ice4. Previous attempts at seasonal forecasts of ice extent were of limited skill5, 6, 7, 8, 9. However, here we show that the Arctic sea-ice minimum can be accurately forecasted from melt-pond area in spring. We find a strong correlation between the spring pond fraction and September sea-ice extent. This is explained by a positive feedback mechanism: more ponds reduce the albedo; a lower albedo causes more melting; more melting increases pond fraction. Our results help explain the acceleration of Arctic sea-ice decrease during the past decade. The inclusion of our new melt-pond model10 promises to improve the skill of future forecast and climate models in Arctic regions and beyond.

Solution calorimetry as a tool for investigating drug interaction with intestinal fluid

Solution calorimetry offers a reproducible technique for measuring the enthalpy of solution (ΔsolH) of a solute dissolving into a solvent. The ΔsolH of two solutes, propranolol HCl and mannitol were determined in simulated intestinal fluid (SIF) solutions designed to model the fed and fasted states within the gut, and in Hanks’ balanced salt solution (HBSS) of varying pH. The bile salt and lipid within the SIF solutions formed mixed micelles. Both solutes exhibited endothermic reactions in all solvents. The ΔsolH for propranolol HCl in the SIF solutions differed from those in the HBSS and was lower in the fed state than the fasted state SIF solution, revealing an interaction between propranolol and the micellar phase in both SIF solutions. In contrast, for mannitol the ΔsolH was constant in all solutions indicating minimal interaction between mannitol and the micellar phases of the SIF solutions. In this study, solution calorimetry proved to be a simple method for measuring the enthalpy associated with the dissolution of model drugs in complex biological media such as SIF solutions. In addition, the derived power–time curves allowed the time taken for the powdered solutes to form solutions to be estimated.

Development of the Surface Urban Energy and Water balance Scheme (SUEWS) for cold climate cities

The Surface Urban Energy and Water Balance Scheme (SUEWS) is developed to include snow. The processes addressed include accumulation of snow on the different urban surface types: snow albedo and density aging, snow melting and re-freezing of meltwater. Individual model parameters are assessed and independently evaluated using long-term observations in the two cold climate cities of Helsinki and Montreal. Eddy covariance sensible and latent heat fluxes and snow depth observations are available for two sites in Montreal and one in Helsinki. Surface runoff from two catchments (24 and 45 ha) in Helsinki and snow properties (albedo and density) from two sites in Montreal are also analysed. As multiple observation sites with different land-cover characteristics are available in both cities, model development is conducted independent of evaluation. The developed model simulates snowmelt related runoff well (within 19% and 3% for the two catchments in Helsinki when there is snow on the ground), with the springtime peak estimated correctly. However, the observed runoff peaks tend to be smoother than the simulated ones, likely due to the water holding capacity of the catchments and the missing time lag between the catchment and the observation point in the model. For all three sites the model simulates the timing of the snow accumulation and melt events well, but underestimates the total snow depth by 18–20% in Helsinki and 29–33% in Montreal. The model is able to reproduce the diurnal pattern of net radiation and turbulent fluxes of sensible and latent heat during cold snow, melting snow and snow-free periods. The largest model uncertainties are related to the timing of the melting period and the parameterization of the snowmelt. The results show that the enhanced model can simulate correctly the exchange of energy and water in cold climate cities at sites with varying surface cover.

Cobalt incorporation in calcite: Thermochemistry of (Ca,Co)CO3 solid solutions from density functional theory simulations

The incorporation of cobalt in mixed metal carbonates is a possible route to the immobilization of this toxic element in the environment. However, the thermodynamics of (Ca,Co)CO3 solid solutions are still unclear due to conflicting data from experiment and from the observation of natural ocurrences. We report here the results of a computer simulation study of the mixing of calcite (CaCO3) and spherocobaltite (CoCO3), using density functional theory calculations. Our simulations suggest that previously proposed thermodynamic models, based only on the range of observed compositions, significantly overestimate the solubility between the two solids and therefore underestimate the extension of the miscibility gap under ambient conditions. The enthalpy of mixing of the disordered solid solution is strongly positive and moderately asymmetric: calcium incorporation in spherocobaltite is more endothermic than cobalt incorporation in calcite. Ordering of the impurities in (0001) layers is energetically favourable with respect to the disordered solid solution at low temperatures and intermediate compositions, but the ordered phase is still unstable to demixing. We calculate the solvus and spinodal lines in the phase diagram using a sub-regular solution model, and conclude that many Ca1-xCoxCO3 mineral solid solutions (with observed compositions of up to x=0.027, and above x=0.93) are metastable with respect to phase separation. We also calculate solid/aqueous distribution coefficients to evaluate the effect of the strong non-ideality of mixing on the equilibrium with aqueous solution, showing that the thermodynamically-driven incorporation of cobalt in calcite (and of calcium in spherocobaltite) is always very low, regardless of the Co/Ca ratio of the aqueous environment.

Enhancing the crystallization and orientation of electrospinning poly (lactic acid) (PLLA) by combining with additives

PLLA is a thermoplastic biopolymer and can be used in industrial applications for medical and filtration applications. The brittleness of PLLA is attributed to slow crystallization rates and its glass transition temperature (Tg) is high (60 °C); for this reason, its applications are limited. The orientation, morphology, and crystal structure of the electrospun fibers was investigated by SEM, POM, DSC, FTIR, XRD, and SAXS. Combining with additives leads to a large decrease of fiber diameter, viscosity, and changes of fiber morphology and crystal structure compared to pure PLLA. DSC showed that the Tg of PLLA decreased about 15 °C and there was no change in relaxation enthalpy by the addition of plasticizer. FT-IR indicate a strong interaction between PLLA and additives; a new band appears in the PLLA blend at 1,756 cm−1 at room temperature as a crystalline band without any annealing. In addition, WAXD indicated that the intensities of the two peaks at (200/110) and (203) increased for the blend at room temperature without any annealing in comparison with PLLA; this means that PHB crystallizes in the amorphous region of PLLA. The POM experiments agree with the results from DSC, FTIR, and WAXS measurements, confirming that adding PHB results in an increase in the number of nuclei with much smaller spherulites and enhances the crystallization behavior of this material, thereby improving its potential for applications.

An AeroCom assessment of black carbon in Arctic snow and sea ice

Though many global aerosols models prognose surface deposition, only a few models have been used to directly simulate the radiative effect from black carbon (BC) deposition to snow and sea ice. Here, we apply aerosol deposition fields from 25 models contributing to two phases of the Aerosol Comparisons between Observations and Models (AeroCom) project to simulate and evaluate within-snow BC concentrations and radiative effect in the Arctic. We accomplish this by driving the offline land and sea ice components of the Community Earth System Model with different deposition fields and meteorological conditions from 2004 to 2009, during which an extensive field campaign of BC measurements in Arctic snow occurred. We find that models generally underestimate BC concentrations in snow in northern Russia and Norway, while overestimating BC amounts elsewhere in the Arctic. Although simulated BC distributions in snow are poorly correlated with measurements, mean values are reasonable. The multi-model mean (range) bias in BC concentrations, sampled over the same grid cells, snow depths, and months of measurements, are −4.4 (−13.2 to +10.7) ng g−1 for an earlier phase of AeroCom models (phase I), and +4.1 (−13.0 to +21.4) ng g−1 for a more recent phase of AeroCom models (phase II), compared to the observational mean of 19.2 ng g−1. Factors determining model BC concentrations in Arctic snow include Arctic BC emissions, transport of extra-Arctic aerosols, precipitation, deposition efficiency of aerosols within the Arctic, and meltwater removal of particles in snow. Sensitivity studies show that the model–measurement evaluation is only weakly affected by meltwater scavenging efficiency because most measurements were conducted in non-melting snow. The Arctic (60–90° N) atmospheric residence time for BC in phase II models ranges from 3.7 to 23.2 days, implying large inter-model variation in local BC deposition efficiency. Combined with the fact that most Arctic BC deposition originates from extra-Arctic emissions, these results suggest that aerosol removal processes are a leading source of variation in model performance. The multi-model mean (full range) of Arctic radiative effect from BC in snow is 0.15 (0.07–0.25) W m−2 and 0.18 (0.06–0.28) W m−2 in phase I and phase II models, respectively. After correcting for model biases relative to observed BC concentrations in different regions of the Arctic, we obtain a multi-model mean Arctic radiative effect of 0.17 W m−2 for the combined AeroCom ensembles. Finally, there is a high correlation between modeled BC concentrations sampled over the observational sites and the Arctic as a whole, indicating that the field campaign provided a reasonable sample of the Arctic.

High-temperature order–disorder transitions in the skutterudites CoGe1.5Q1.5 (Q=S, Te)

The temperature dependence of anion ordering in the skutterudites CoGe1.5Q1.5 (Q=S, Te) has been investigated by powder neutron diffraction. Both materials adopt a rhombohedral structure at room temperature (space group R-3 ) in which the anions are ordered trans to each other within Ge2Q2 rings. In CoGe1.5S1.5, anion ordering is preserved up to the melting point of 950 °C. However, rhombohedral CoGe1.5Te1.5 undergoes a phase transition at 610 °C involving a change to cubic symmetry (space group Im-3). In the high-temperature modification, there is a statistical distribution of anions over the available sites within the Ge2Te2 rings. The structural transition involves a reduction in the degree of distortion of the Ge2Te2 rings which progressively transform from a rhombus to a rectangular shape. The effect of this transition on the thermoelectric properties has been investigated.

Stabilization of the i-motif structure by the intra-strand cross-link formation

The i-motif structures are formed by oligonucleotides containing cytosine tracts under acidic conditions. The folding of the i-motif under physiological conditions is of great interest because of its biological role. In this study, we investigated the effect of the intra-strand cross-link on the stability of the i-motif structure. The 4-vinyl-substituted analog of thymidine (T-vinyl) was incorporated into the 5′-end of the human telomere complementary strand, which formed the intra-strand cross-link with the internal adenine. The intra-strand cross-linked i-motif displayed CD spectra similar to that of the natural i-motif at acidic pH, which was transformed into a random coil with the increasing pH. The pH midpoint for the transition from the i-motif to random coil increased from pH 6.1 for the natural one to pH 6.8 for the cross-linked one. The thermodynamic parameters were obtained by measuring the thermal melting behaviors by CD and UV, and it was determined that the intra-strand cross-linked i-motif is stabilized due to a favorable entropy effect. Thus, this study has clearly indicated the validity of the intra-strand cross-linking for stabilization of the i-motif structure.

Observed anomalous atmospheric patterns in summers of unusual Arctic sea ice melt

The Arctic sea ice retreat has accelerated over the last decade. The negative trend is largest in summer, but substantial interannual variability still remains. Here we explore observed atmospheric conditions and feedback mechanisms during summer months of anomalous sea ice melt in the Arctic. Compositing months of anomalous low and high sea ice melt over 1979–2013, we find distinct patterns in atmospheric circulation, precipitation, radiation, and temperature. Compared to summer months of anomalous low sea ice melt, high melt months are characterized by anomalous high sea level pressure in the Arctic (up to 7 hPa), with a corresponding tendency of storms to track on a more zonal path. As a result, the Arctic receives less precipitation overall and 39% less snowfall. This lowers the albedo of the region and reduces the negative feedback the snowfall provides for the sea ice. With an anticyclonic tendency, 12 W/m2 more incoming shortwave radiation reaches the surface in the start of the season. The melting sea ice in turn promotes cloud development in the marginal ice zones and enhances downwelling longwave radiation at the surface toward the end of the season. A positive cloud feedback emerges. In midlatitudes, the more zonally tracking cyclones give stormier, cloudier, wetter, and cooler summers in most of northern Europe and around the Sea of Okhotsk. Farther south, the region from the Mediterranean Sea to East Asia experiences significant surface warming (up to 2.4◦C), possibly linked to changes in the jet stream.