Glacial mass balance changes in the Karakoram and Himalaya based on CMIP5 multi-model climate projections

The impact of future climate change on the glaciers in the Karakoram and Himalaya (KH) is investigated using CMIP5 multi-model temperature and precipitation projections, and a relationship between glacial accumulation-area ratio and mass balance developed for the region based on the last 30 to 40 years of observational data. We estimate that the current glacial mass balance (year 2000) for the entire KH region is -6.6 +/- 1 Gta(-1), which decreases about sixfold to -35 +/- 2 Gta(-1) by the 2080s under the high emission scenario of RCP8.5. However, under the low emission scenario of RCP2.6 the glacial mass loss only doubles to -12 +/- 2 Gta(-1) by the 2080s. We also find that 10.6 and 27 % of the glaciers could face `eventual disappearance' by the end of the century under RCP2.6 and RCP8.5 respectively, underscoring the threat to water resources under high emission scenarios.

Non-Gravitational Effects with Density-Matching in Evaluating the Influence of Sedimentation on Colloidal Coagulation

The method of density matching between the solid and liquid phases is often adopted to effectively eliminate the effect of sedimentation of suspensions in studies on dynamic behaviour of a colloidal system. However, the associated changes in the solvent composition may bring side effects to the properties investigated and therefore might lead to a faulty conclusion if the relevant correction is not made. To illustrate the importance of this side effect, we present an example of the sedimentation influence on the coagulation rate of suspensions of 2 μm (diameter) polystyrene. The liquid mixtures, in the proper proportions of water (H2O), deuterium oxide (D2O) and methanol (MeOH) as the liquid phase, density-matched and unmatched experiments are performed. Besides the influence of viscosity, the presence of methanol in solvent media, used to enhance the sedimentation effect, causes significant changes (reduction) in rapid coagulation rates compared to that in pure water. Without the relevant corrections for those non-gravitational factors it seems that gravitational sedimentation would retard the coagulation. The magnitude of the contribution from the non-gravitational factor is quantitatively determined, making the relevant correction possible. After necessary corrections for all factors, our experiments show that the influence of the sedimentation on coagulation rates at the initial stage of the coagulation is not observable.

Towards an Understanding of the Influence of Sedimentation on Colloidal Aggregation by Peclet Number

The Peclet number is a useful index to estimate the importance of sedimentation as compared to the Brownian motion. However, how to choose the characteristic length scale for the Peclet number evaluation is rather critical because the diffusion length increases as the square root of the time whereas the drifting length is linearly related to time. Our Brownian dynamics simulation shows that the degree of sedimentation influence on the coagulation decreases when the dispersion volume fraction increases. Therefore using a fixed length, such as the diameter of particle, as the characteristic length scale for Peclet number evaluation is not a good choice when dealing with the influence of sedimentation on coagulation. The simulations demonstrated that environmental factors in the coagulation process, such as dispersion volume fraction and size distribution, should be taken into account for more reasonable evaluation of the sedimentation influence.

Sedimentation Process of Saturated Sand under Impact Loading

The initial small inhomogeneity of saturated sand could be amplified during the sedimentation process after liquefaction, and cracks could be observed in the sand column. Layers of fine sand could also be found at the exact place where cracks developed and disappeared. The phenomena and the whole process were experimentally shown by X-rays images. To account for the phenomena, a linearized stability analysis of the sedimentation of saturated sand was conducted; however, it did not produce a satisfactory result. A three-phase flow model describing the transportation of fine sand is presented in this paper. It is shown that such a kind of erosion/deposition model was qualitatively in good agreement with the experimental observation.

Single Curved Fiber Sedimentation Under Gravity

Dynamics of single curved fiber sedimentation under gravity are simulated by using the lattice Boltzmann method. The results of migration and rotation of the curved fiber at different Reynolds numbers are reported. The results show that the rotation and migration processes are sensitive to the curvature of the fiber. (c) 2007 Elsevier Ltd. All rights reserved.

Computer simulation of influence of sedimentation on rapid coagulation

A computer simulation was performed to explore the features and effects of sedimentation on rapid coagulation. To estimate the accumulated influence of gravity on coagulation for dispersions, a sedimentation influence ratio is defined. Some factors possibly related to the influence of sedimentation were considered in the simulation and analysed by comparing the size distribution of aggregates, the change in collision number, and coagulation rates at different gravity levels (0 g, 1 g and more with g being the gravitational constant).

A coupled model for multiphase fluid flow and sedimentation deformation in oil reservoir and its numerical simulation

A mathematical model for coupled multiphase fluid flow and sedimentation deformation is developed based on fluid-solid interaction mechanism. A finite difference-finite element numerical approach is presented. The results of an example show that the fluid-solid coupled effect has great influence on multiphase fluid flow and reservoir recovery performances, and the coupled model has practical significance for oilfield development.

Single curved fiber sedimentation under gravity

Dynamics of single curved fiber sedimentation under gravity are simulated by using the lattice Boltzmann method. The results of migration and rotation of the curved fiber at different Reynolds numbers are reported. The results show that the rotation and migration processes are sensitive to the curvature of the fiber. (c) 2007 Elsevier Ltd. All rights reserved.

An algorithm for coarse particle sedimentation simulation by Stokesian dynamics

Coarse Particle sedimentation is studied by using an algorithm with no adjustable parameters based on stokesian dynamics. Only inter-particle interactions of hydrodynamic force and gravity are considered. The sedimentation of a simple cubic array of spheres is used to verify the computational results. The scaling and parallelism with OpenMP of the method are presented. Random suspension sedimentation is investigated with Mont Carlo simulation. The computational results are shown in good agreement with experimental fitting at the lower computational cost of O(N In N).

Movement of a Spherical Brownian Particle Between Infinite Parallel Plates: Hindered Dispersion and Sedimentation

The dispersion of an isolated, spherical, Brownian particle immersed in a Newtonian fluid between infinite parallel plates is investigated. Expressions are developed for both a 'molecular' contribution to dispersion, which arises from random thermal fluctuations, and a 'convective' contribution, arising when a shear flow is applied between the plates. These expressions are evaluated numerically for all sizes of the particle relative to the bounding plates, and the method of matched asymptotic expansions is used to develop analytical expressions for the dispersion coefficients as a function of particle size to plate spacing ratio for small values of this parameter.

It is shown that both the molecular and convective dispersion coefficients decrease as the size of the particle relative to the bounding plates increase. When the particle is small compared to the plate spacing, the coefficients decrease roughly proportional to the particle size to plate spacing ratio. When the particle closely fills the space between the plates, the molecular dispersion coefficient approaches zero slowly as an inverse logarithmic function of the particle size to plate spacing ratio, and the convective dispersion coefficent approaches zero approximately proportional to the width of the gap between the edges of the sphere and the bounding plates.

The short-timescale behavior of glacial ice

Glaciers are often assumed to deform only at slow (i.e., glacial) rates. However, with the advent of high rate geodetic observations of ice motion, many of the intricacies of glacial deformation on hourly and daily timescales have been observed and quantified. This thesis explores two such short timescale processes: the tidal perturbation of ice stream motion and the catastrophic drainage of supraglacial meltwater lakes. Our investigation into the transmission length-scale of a tidal load represents the first study to explore the daily tidal influence on ice stream motion using three-dimensional models. Our results demonstrate both that the implicit assumptions made in the standard two-dimensional flow-line models are inherently incorrect for many ice streams, and that the anomalously large spatial extent of the tidal influence seen on the motion of some glaciers cannot be explained, as previously thought, through the elastic or viscoelastic transmission of tidal loads through the bulk of the ice stream. We then discuss how the phase delay between a tidal forcing and the ice stream’s displacement response can be used to constrain in situ viscoelastic properties of glacial ice. Lastly, for the problem of supraglacial lake drainage, we present a methodology for implementing linear viscoelasticity into an existing model for lake drainage. Our work finds that viscoelasticity is a second-order effect when trying to model the deformation of ice in response to a meltwater lake draining to a glacier’s bed. The research in this thesis demonstrates that the first-order understanding of the short-timescale behavior of naturally occurring ice is incomplete, and works towards improving our fundamental understanding of ice behavior over the range of hours to days.

The deep ocean density structure at the Last Glacial Maximum : What was it and why?

The search for reliable proxies of past deep ocean temperature and salinity has proved difficult, thereby limiting our ability to understand the coupling of ocean circulation and climate over glacial-interglacial timescales. Previous inferences of deep ocean temperature and salinity from sediment pore fluid oxygen isotopes and chlorinity indicate that the deep ocean density structure at the Last Glacial Maximum (LGM, approximately 20,000 years BP) was set by salinity, and that the density contrast between northern and southern sourced deep waters was markedly greater than in the modern ocean. High density stratification could help explain the marked contrast in carbon isotope distribution recorded in the LGM ocean relative to that we observe today, but what made the ocean's density structure so different at the LGM? How did it evolve from one state to another? Further, given the sparsity of the LGM temperature and salinity data set, what else can we learn by increasing the spatial density of proxy records?

We investigate the cause and feasibility of a highly and salinity stratified deep ocean at the LGM and we work to increase the amount of information we can glean about the past ocean from pore fluid profiles of oxygen isotopes and chloride. Using a coupled ocean--sea ice--ice shelf cavity model we test whether the deep ocean density structure at the LGM can be explained by ice--ocean interactions over the Antarctic continental shelves, and show that a large contribution of the LGM salinity stratification can be explained through lower ocean temperature. In order to extract the maximum information from pore fluid profiles of oxygen isotopes and chloride we evaluate several inverse methods for ill-posed problems and their ability to recover bottom water histories from sediment pore fluid profiles. We demonstrate that Bayesian Markov Chain Monte Carlo parameter estimation techniques enable us to robustly recover the full solution space of bottom water histories, not only at the LGM, but through the most recent deglaciation and the Holocene up to the present. Finally, we evaluate a non-destructive pore fluid sampling technique, Rhizon samplers, in comparison to traditional squeezing methods and show that despite their promise, Rhizons are unlikely to be a good sampling tool for pore fluid measurements of oxygen isotopes and chloride.

Long-term natural acidification of upland sites in Cumbria: Evidence from post-glacial lake sediments

Studies by the Freshwater Biological Association over the last 25 years have supplied data relevant to the levels of acidity in local soils and water before the onset of industrial pollution and current interest in acid rain. This article reviews published analysis from cores of lake sediments, in or near the catchment of the River Duddon. Electron spin resonance spectra of humic acids and iodine values confirm evidence from pollen analysis for a history of progressive acidification of the source material of lake sediments since before 5000 radiocarbon years, in upland catchments of the Lake District. Processes involved included: removal of basic ions from soils by rainfall, the effects of which were intensified by removal by man of deciduous forest; acidification of soils and waters by decomposition products of Calluna and further acidification of waters by Sphagnum species which colonized habitats where drainage became impeded by paludification processes.

The significance of sedimentation and sediments to phytoplankton growth in drinking-water reservoirs

In the mesotrophic-eutrophic Saidenbach Reservoir in Saxony, the nanoplankton and cyanobacteria have increased at the expense of diatom dominance, due to a doubling of the external phosphorus load in the last 15 years. However, the phosphorus sedimentation flux is still very high (up to 80% of the input), corresponding to more than 2 g m2 d-1 in terms of dry weight. There is a strong correlation between the abundance of diatoms in the euphotic zone and their sedimentation flux (with a delay of about 2 weeks). Only about 25% of the deposited material could be clearly attributed to plankton biomass; the remainder resulted from flocculation and precipitation processes or directly from the inflow of clay minerals. The ash content of the deposited material was high (73%). Thus the sedimentation flux can be considered to operate as an internal water-treatment/oligotrophication process within the lake. The neighbouring Neunzehnhain Reservoir still has a very clear water with a transparency up to 18 m depth. Though the sediment was not much lower than Saidenbach sediment in total phosphorus and total numbers of bacteria, sulphide was always absent and the ratio of Fe 2+ to Fe 3+ was very low in the upper (0- 5 cm) layer. Thus the external and internal phosphorus loads do not attain the critical level necessary to induce a ”phosphorus - phytoplankton” feedback loop.