Geodetic mass balance of Azarova glacier, Kodar Mountains, eastern Siberia, and its links to observed and projected climatic change

The Kodar Mountains in eastern Siberia accommodate 30 small, cold-based glaciers with a combined surface area of about 19 km2. Very little is known about these glaciers, with the first survey conducted in the late 1950s. In this paper, we use terrestrial photogrammetry to calculate changes in surface area, elevation, volume and geodetic mass balance of the Azarova Glacier between 1979 and 2007 and relate these to meteorological data from nearby Chara weather station (1938-2007). The glacier surface area declined by 20±6.9% and surface lowered on average by 20±1.8 m (mean thinning: 0.71 m a-1) resulting in a strongly negative cumulative and average mass balance of -18±1.6 m w.e. and -640±60 mm w.e.a-1 respectively. The July-August air temperature increased at a rate of 0.036oC a-1 between 1979 and 2007 and the 1980-2007 period was, on average, around 1oC warmer than 1938-1979. The regional climate projections for A2 and B2 CO2 emission scenarios developed using PRECIS regional climate model indicate that summer temperatures will increase in 2071–2100 by 2.6-4.7°C and 4.9-6.2°C respectively in comparison with 1961–1990. The annual total of solid precipitation will increase by 20% under B2 scenario but decline by 3% under A2 scenario. The length of the ablation season will extend from July–August to June-September. The Azarova Glacier exhibits high sensitivity to climatic warming due to its low elevation, exposure to comparatively high summer temperatures, and the absence of a compensating impact of cold season precipitation. Further summer warming and decline of solid precipitation projected under the A2 scenario will force Azarova to retreat further while impacts of an increase in solid precipitation projected under the B2 scenario require further investigation.

The advection of high-resolution tracers by low-resolution winds

The usefulness of any simulation of atmospheric tracers using low-resolution winds relies on both the dominance of large spatial scales in the strain and time dependence that results in a cascade in tracer scales. Here, a quantitative study on the accuracy of such tracer studies is made using the contour advection technique. It is shown that, although contour stretching rates are very insensitive to the spatial truncation of the wind field, the displacement errors in filament position are sensitive. A knowledge of displacement characteristics is essential if Lagrangian simulations are to be used for the inference of airmass origin. A quantitative lower estimate is obtained for the tracer scale factor (TSF): the ratio of the smallest resolved scale in the advecting wind field to the smallest “trustworthy” scale in the tracer field. For a baroclinic wave life cycle the TSF = 6.1 ± 0.3 while for the Northern Hemisphere wintertime lower stratosphere the TSF = 5.5 ± 0.5, when using the most stringent definition of the trustworthy scale. The similarity in the TSF for the two flows is striking and an explanation is discussed in terms of the activity of potential vorticity (PV) filaments. Uncertainty in contour initialization is investigated for the stratospheric case. The effect of smoothing initial contours is to introduce a spinup time, after which wind field truncation errors take over from initialization errors (2–3 days). It is also shown that false detail from the proliferation of finescale filaments limits the useful lifetime of such contour advection simulations to 3σ−1 days, where σ is the filament thinning rate, unless filaments narrower than the trustworthy scale are removed by contour surgery. In addition, PV analysis error and diabatic effects are so strong that only PV filaments wider than 50 km are at all believable, even for very high-resolution winds. The minimum wind field resolution required to accurately simulate filaments down to the erosion scale in the stratosphere (given an initial contour) is estimated and the implications for the modeling of atmospheric chemistry are briefly discussed.

Rheology of magmas with bimodal crystal size and shape distributions: insights from analog experiments

Magmas in volcanic conduits commonly contain microlites in association with preexisting phenocrysts, as often indicated by volcanic rock textures. In this study, we present two different experiments that inves- tigate the flow behavior of these bidisperse systems. In the first experiments, rotational rheometric methods are used to determine the rheology of monodisperse and polydisperse suspensions consisting of smaller, prolate particles (microlites) and larger, equant particles (phenocrysts) in a bubble‐free Newtonian liquid (silicate melt). Our data show that increasing the relative proportion of prolate microlites to equant pheno- crysts in a magma at constant total particle content can increase the relative viscosity by up to three orders of magnitude. Consequently, the rheological effect of particles in magmas cannot be modeled by assuming a monodisperse population of particles. We propose a new model that uses interpolated parameters based on the relative proportions of small and large particles and produces a considerably improved fit to the data than earlier models. In a second series of experiments we investigate the textures produced by shearing bimodal suspensions in gradually solidifying epoxy resin in a concentric cylinder setup. The resulting textures show the prolate particles are aligned with the flow lines and spherical particles are found in well‐organized strings, with sphere‐depleted shear bands in high‐shear regions. These observations may explain the measured variation in the shear thinning and yield stress behavior with increasing solid fraction and particle aspect ratio. The implications for magma flow are discussed, and rheological results and tex- tural observations are compared with observations on natural samples.

Rheology of milk foams produced by steam injection

Rheology of milk foams generated by steam injection was studied during the transient destabilization process using steady flow and dynamic oscillatory techniques: yield stress (τ_y) values were obtained from a stress ramp (0.2 to 25 Pa) and from strain amplitude sweep (0.001 to 3 at 1 Hz of frequency); elastic (G') and viscous (G") moduli were measured by frequency sweep (0.1 to 150 Hz at 0.05 of strain); and the apparent viscosity (η_a) was obtained from the flow curves generated from the stress ramp. The effect of plate roughness and the sweep time on τ_y was also assessed. Yield stress was found to increase with plate roughness whereas it decreased with the sweep time. The values of yield stress and moduli—G' and G"—increased during foam destabilization as a consequence of the changes in foam properties, especially the gas volume fraction, φ, and bubble size, R_32 (Sauter mean bubble radius). Thus, a relationship between τ_y, φ, R_32, and σ (surface tension) was established. The changes in the apparent viscosity, η, showed that the foams behaved like a shear thinning fluid beyond the yield point, fitting the modified Cross model with the relaxation time parameter (λ) also depending on the gas volume fraction. Overall, it was concluded that the viscoelastic behavior of the foam below the yield point and liquid-like behavior thereafter both vary during destabilization due to changes in the foam characteristics.

Self-assembly of a peptide amphiphile containing l-carnosine and its mixtures with a multilamellar vesicle forming lipid

The self-assembly of the peptide amphiphile (PA) hexadecyl-(β-alaninehistidine) is examined in aqueous solution, along with its mixtures with multilamellar vesicles formed by DPPC (dipalmitoyl phosphatidylcholine). This PA, denoted C16-βAH, contains a dipeptide headgroup corresponding to the bioactive molecule L-carnosine. It is found to selfassemble into nanotapes based on stacked layers of molecules. Bilayers are found to coexist with monolayers in which the PA molecules pack with alternating up−down arrangement so that the headgroups decorate both surfaces. The bilayers become dehydrated as PA concentration increases and the number of layers in the stack decreases to produce ultrathin nanotapes comprised of 2−3 bilayers. Addition of the PA to DPPC multilamellar vesicles leads to a transition to well-defined unilamellar vesicles. The unique ability to modulate the stacking of this PA as a function of concentration, combined with its ability to induce a multilamellar to unilamellar thinning of DPPC vesicles, may be useful in biomaterials applications where the presentation of the peptide function at the surface of self-assembled nanostructures is crucial.

Estimating the volume of tephra deposits: a new simple strategy

Volume determination of tephra deposits is necessary for the assessment of the dynamics and hazards of explosive volcanoes. Several methods have been proposed during the past 40 years that include the analysis of crystal concentration of large pumices, integrations of various thinning relationships, and the inversion of field observations using analytical and computational models. Regardless of their strong dependence on tephra-deposit exposure and distribution of isomass/isopach contours, empirical integrations of deposit thinning trends still represent the most widely adopted strategy due to their practical and fast application. The most recent methods involve the best fitting of thinning data using various exponential seg- ments or a power-law curve on semilog plots of thickness (or mass/area) versus square root of isopach area. The exponential method is mainly sensitive to the number and the choice of straight segments, whereas the power-law method can better reproduce the natural thinning of tephra deposits but is strongly sensitive to the proximal or distal extreme of integration. We analyze a large data set of tephra deposits and propose a new empirical method for the deter- mination of tephra-deposit volumes that is based on the integration of the Weibull function. The new method shows a better agreement with observed data, reconciling the debate on the use of the exponential versus power-law method. In fact, the Weibull best fitting only depends on three free parameters, can well reproduce the gradual thinning of tephra deposits, and does not depend on the choice of arbitrary segments or of arbitrary extremes of integration.

Data assimilation with correlated observation errors: experiments with a 1-D shallow water model

Remote sensing observations often have correlated errors, but the correlations are typically ignored in data assimilation for numerical weather prediction. The assumption of zero correlations is often used with data thinning methods, resulting in a loss of information. As operational centres move towards higher-resolution forecasting, there is a requirement to retain data providing detail on appropriate scales. Thus an alternative approach to dealing with observation error correlations is needed. In this article, we consider several approaches to approximating observation error correlation matrices: diagonal approximations, eigendecomposition approximations and Markov matrices. These approximations are applied in incremental variational assimilation experiments with a 1-D shallow water model using synthetic observations. Our experiments quantify analysis accuracy in comparison with a reference or ‘truth’ trajectory, as well as with analyses using the ‘true’ observation error covariance matrix. We show that it is often better to include an approximate correlation structure in the observation error covariance matrix than to incorrectly assume error independence. Furthermore, by choosing a suitable matrix approximation, it is feasible and computationally cheap to include error correlation structure in a variational data assimilation algorithm.

A Lagrangian model of air-mass photochemistry and mixing using a trajectory ensemble: the Cambridge Tropospheric Trajectory model of Chemistry And Transport (CiTTyCAT) version 4.2

A Lagrangian model of photochemistry and mixing is described (CiTTyCAT, stemming from the Cambridge Tropospheric Trajectory model of Chemistry And Transport), which is suitable for transport and chemistry studies throughout the troposphere. Over the last five years, the model has been developed in parallel at several different institutions and here those developments have been incorporated into one "community" model and documented for the first time. The key photochemical developments include a new scheme for biogenic volatile organic compounds and updated emissions schemes. The key physical development is to evolve composition following an ensemble of trajectories within neighbouring air-masses, including a simple scheme for mixing between them via an evolving "background profile", both within the boundary layer and free troposphere. The model runs along trajectories pre-calculated using winds and temperature from meteorological analyses. In addition, boundary layer height and precipitation rates, output from the analysis model, are interpolated to trajectory points and used as inputs to the mixing and wet deposition schemes. The model is most suitable in regimes when the effects of small-scale turbulent mixing are slow relative to advection by the resolved winds so that coherent air-masses form with distinct composition and strong gradients between them. Such air-masses can persist for many days while stretching, folding and thinning. Lagrangian models offer a useful framework for picking apart the processes of air-mass evolution over inter-continental distances, without being hindered by the numerical diffusion inherent to global Eulerian models. The model, including different box and trajectory modes, is described and some output for each of the modes is presented for evaluation. The model is available for download from a Subversion-controlled repository by contacting the corresponding authors.

Cortical thickness changes in the non-lesioned hemisphere associated with non-paretic arm immobilization in modified CI therapy

Recent evidence suggests that immobilization of the upper limb for 2–3 weeks induces changes in cortical thickness as well as motor performance. In constraint induced (CI) therapy, one of the most effective interventions for hemiplegia, the non-paretic arm is constrained to enforce the use of the paretic arm in the home setting. With the present study we aimed to explore whether non-paretic arm immobilization in CI therapy induces structural changes in the non-lesioned hemisphere, and how these changes are related to treatment benefit. 31 patients with chronic hemiparesis participated in CI therapy with (N = 14) and without (N = 17) constraint. Motor ability scores were acquired before and after treatment. Diffusion tensor imaging (DTI) data was obtained prior to treatment. Cortical thickness was measured with the Freesurfer software. In both groups cortical thickness in the contralesional primary somatosensory cortex increased and motor function improved with the intervention. However the cortical thickness change was not associated with the magnitude of motor function improvement. Moreover, the treatment effect and the cortical thickness change were not significantly different between the constraint and the non-constraint groups. There was no correlation between fractional anisotropy changes in the non-lesioned hemisphere and treatment outcome. CI therapy induced cortical thickness changes in contralesional sensorimotor regions, but this effect does not appear to be driven by the immobilization of the non-paretic arm, as indicated by the absence of differences between the constraint and the non-constraint groups. Our data does not suggest that the arm immobilization used in CI therapy is associated with noticeable cortical thinning.

Sources and spread of thermophilic Campylobacter spp. during partial depopulation of broiler chicken flocks

The practice of partial depopulation or ‘thinning’, i.e. early removal of a proportion of birds from a commercial broiler flock, is a reported risk factor for Campylobacter colonization of residual birds because of the difficulty in maintaining biosecurity during the process. Therefore, the effect of this practice was studied in detail for 51 target flocks, each at a different growing farm belonging to one of seven major poultry companies throughout the United Kingdom. On 21 of these farms, the target flock was already colonized by Campylobacter and at slaughter all cecal samples examined were positive, with a mean of log10 8 cfu / g. A further 27 flocks became positive within 2 – 6 days of the start of thinning and had similarly high levels of cecal carriage at slaughter. Just prior to the thinning process, Campylobacter could be isolated frequently from the farm driveways, transport vehicles, equipment and personnel. Strains from seven such farms on which flocks became colonized after thinning were examined by PFGE typing. The study demonstrated an association between strains occurring at specific sampling sites and those isolated subsequently from the thinned flocks. There were also indications that particular strains had spread from one farm to another, when the farms were jointly company-owned and served by the same bird-catching teams and / or vehicles. The results highlighted the need for better hygiene control in relation to catching equipment and personnel, and more effective cleaning and disinfection of vehicles, and bird-transport crates.

Antarctic Peninsula ice sheet evolution during the Cenozoic Era

The Antarctic Peninsula region is currently undergoing rapid environmental change, resulting in the thinning, acceleration and recession of glaciers and the sequential collapse of ice shelves. It is important to view these changes in the context of long-term palaeoenvironmental complexity and to understand the key processes controlling ice sheet growth and recession. In addition, numerical ice sheet models require detailed geological data for tuning and testing. Therefore, this paper systematically and holistically reviews published geological evidence for Antarctic Peninsula Ice Sheet variability for each key locality throughout the Cenozoic, and brings together the prevailing consensus of the extent, character and behaviour of the glaciations of the Antarctic Peninsula region. Major contributions include a downloadable database of 186 terrestrial and marine calibrated dates; an original reconstruction of the LGM ice sheet; and a new series of isochrones detailing ice sheet retreat following the LGM. Glaciation of Antarctica was initiated around the Eocene/Oligocene transition in East Antarctica. Palaeogene records of Antarctic Peninsula glaciation are primarily restricted to King George Island, where glacigenic sediments provide a record of early East Antarctic glaciations, but with modification of far-travelled erratics by local South Shetland Island ice caps. Evidence for Neogene glaciation is derived primarily from King George Island and James Ross Island, where glaciovolcanic strata indicate that ice thicknesses reached 500–850 m during glacials. This suggests that the Antarctic Peninsula Ice Sheet draped, rather than drowned, the topography. Marine geophysical investigations indicate multiple ice sheet advances during this time. Seismic profiling of continental shelf-slope deposits indicates up to ten large advances of the Antarctic Peninsula Ice Sheet during the Early Pleistocene, when the ice sheet was dominated by 40 kyr cycles. Glacials became more pronounced, reaching the continental shelf edge, and of longer duration during the Middle Pleistocene. During the Late Pleistocene, repeated glacials reached the shelf edge, but ice shelves inhibited iceberg rafting. The Last Glacial Maximum (LGM) occurred at 18 ka BP, after which transitional glaciomarine sediments on the continental shelf indicate ice-sheet retreat. The continental shelf contains large bathymetric troughs, which were repeatedly occupied by large ice streams during Pleistocene glaciations. Retreat after the LGM was episodic in the Weddell Sea, with multiple readvances and changes in ice-flow direction, but rapid in the Bellingshausen Sea. The late Holocene Epoch was characterised by repeated fluctuations in palaeoenvironmental conditions, with associated glacial readvances. However, this has been subsumed by rapid warming and ice-shelf collapse during the twentieth century.

Late-Holocene changes in character and behaviour of land-terminating glaciers on James Ross Island, Antarctica

Virtually no information is available on the response of land-terminating Antarctic Peninsula glaciers to climate change on a centennial timescale. This paper analyses the topography, geomorphology and sedimentology of prominent moraines on James Ross Island, Antarctica, to determine geometric changes and to interpret glacier behaviour. The moraines are very likely due to a late-Holocene phase of advance and featured (1) shearing and thrusting within the snout, (2) shearing and deformation of basal sediment, (3) more supraglacial debris than at present and (4) short distances of sediment transport. Retreat of ∼100 m and thinning of 15–20 m has produced a loss of 0.1 km3 of ice. The pattern of surface lowering is asymmetric. These geometrical changes are suggested most simply to be due to a net negative mass balance caused by a drier climate. Comparisons of the moraines with the current glaciological surface structure of the glaciers permits speculation of a transition from a polythermal to a cold-based thermal regime. Small land-terminating glaciers in the northern Antarctic Peninsula region could be cooling despite a warming climate.

The effects of maize planting density and weeding regimes on light and water use

Light and water are among essential resources required for production of photosynthates in plants. A study on the effects of weeding regimes and maize planting density on light and water use was conducted during the 2001/2 short and 2002 long rain seasons at Muguga in - the central highlands of Kenya. Weeding regimes were: weed free (W1), weedy (W2), herbicide (W3) and hand weeding twice (W4). Maize planting densities were 9 (D1) and 18 plants m-2 (D2) intercropped with Phaseolus vulgaris (beans). The experiment was laid as randomized complete block design replicated four times and repeated twice. All plots were thinned to 4 plants m-2 at tasseling stage (96 DAE) and thinnings quantified as forage. Soil moisture content (SMC), photosynthetically active radiation (PAR) interception, evapo-transpiration (ET crop), water use efficiency (WUE), and harvest index (HI), were determined. Percent PAR was higher in D2 than in D1 before thinning but higher in D1 than in D2 after thinning in both seasons. PAR interception was highest in W2 but similar in W1, W3 and W4 in both seasons. SMC was significantly lower in W2 but similar in W1, W3 and W4. D2 had lower SMC than D1 in season two. Weeding regime significantly influenced ET crop, while planting density and weeding regime significantly influenced WUE and HI. D2 maximizes water and light use for forage production but results to increased intra-specific plant competition for water and light severely before thinning (96 DAE) that reduce grain yield in dual purpose maize, relative to D1.

Diagnosing observation error correlations for Doppler radar radial winds in the Met Office UKV model using observation-minus-background and observation-minus-analysis statistics

With the development of convection-permitting numerical weather prediction the efficient use of high resolution observations in data assimilation is becoming increasingly important. The operational assimilation of these observations, such as Dopplerradar radial winds, is now common, though to avoid violating the assumption of un- correlated observation errors the observation density is severely reduced. To improve the quantity of observations used and the impact that they have on the forecast will require the introduction of the full, potentially correlated, error statistics. In this work, observation error statistics are calculated for the Doppler radar radial winds that are assimilated into the Met Office high resolution UK model using a diagnostic that makes use of statistical averages of observation-minus-background and observation-minus-analysis residuals. This is the first in-depth study using the diagnostic to estimate both horizontal and along-beam correlated observation errors. By considering the new results obtained it is found that the Doppler radar radial wind error standard deviations are similar to those used operationally and increase as the observation height increases. Surprisingly the estimated observation error correlation length scales are longer than the operational thinning distance. They are dependent on both the height of the observation and on the distance of the observation away from the radar. Further tests show that the long correlations cannot be attributed to the use of superobservations or the background error covariance matrix used in the assimilation. The large horizontal correlation length scales are, however, in part, a result of using a simplified observation operator.

The refreezing of melt ponds on Arctic sea ice

The presence of melt ponds on the surface of Arctic sea ice significantly reduces its albedo, inducing a positive feedback leading to sea ice thinning. While the role of melt ponds in enhancing the summer melt of sea ice is well known, their impact on suppressing winter freezing of sea ice has, hitherto, received less attention. Melt ponds freeze by forming an ice lid at the upper surface, which insulates them from the atmosphere and traps pond water between the underlying sea ice and the ice lid. The pond water is a store of latent heat, which is released during refreezing. Until a pond freezes completely, there can be minimal ice growth at the base of the underlying sea ice. In this work, we present a model of the refreezing of a melt pond that includes the heat and salt balances in the ice lid, trapped pond, and underlying sea ice. The model uses a two-stream radiation model to account for radiative scattering at phase boundaries. Simulations and related sensitivity studies suggest that trapped pond water may survive for over a month. We focus on the role that pond salinity has on delaying the refreezing process and retarding basal sea ice growth. We estimate that for a typical sea ice pond coverage in autumn, excluding the impact of trapped ponds in models overestimates ice growth by up to 265 million km3, an overestimate of 26%.