Geochemical reconstruction of late Holocene drainage and mixing in Kluane Lake,Yukon Territory

The level of Kluane Lake in southwest Yukon Territory, Canada, has fluctuated tens of metres during the late Holocene. Contributions of sediment from different watersheds in the basin over the past 5,000 years were inferred from the elemental geochemistry of Kluane Lake sediment cores. Elements associated with organic material and oxyhydroxides were used to reconstruct redox fluctuations in the hypolimnion of the lake. The data reveal complex relationships between climate and river discharge during the late Holocene. A period of influx of Duke River sediment coincides with a relatively warm climate around 1,300 years BP. Discharge of Slims River into Kluane Lake occurred when Kaskawulsh Glacier advanced to the present drainage divide separating flow to the Pacific Ocean via Kaskawulsh and Alsek rivers from flow to Bering Sea via tributaries of Yukon River. During periods when neither Duke nor Slims river discharged into Kluane Lake, the level of the lake was low and stable thermal stratification developed, with anoxic and eventually euxinic conditions in the hypolimnion.     

An inter-comparison of tidal solutions computed with a range of unstructured grid models of the Irish and Celtic Sea Regions

Three finite element codes, namely TELEMAC, ADCIRC and QUODDY, are used to compute the spatial distributions of the M₂, M₄ and M₆ components of the tide in the sea region off the west coast of Britain. This region is chosen because there is an accurate topographic dataset in the area and detailed open boundary M₂ tidal forcing for driving the model. In addition, accurate solutions (based upon comparisons with extensive observations) using uniform grid finite difference models forced with these open boundary data exist for comparison purposes. By using boundary forcing, bottom topography and bottom drag coefficients identical to those used in an earlier finite difference model, there is no danger of comparing finite element solutions for “untuned unoptimised solutions” with those from a “tuned optimised solution”. In addition, by placing the open boundary in all finite element calculations at the same location as that used in a previous finite difference model and using the same M₂ tidal boundary forcing and water depths, a like with like comparison of solutions derived with the various finite element models was possible. In addition, this open boundary was well removed from the shallow water region, namely the eastern Irish Sea where the higher harmonics were generated. Since these are not included in the open boundary, forcing their generation was determined by physical processes within the models. Consequently, an inter-comparison of these higher harmonics generated by the various finite element codes gives some indication of the degree of variability in the solution particularly in coastal regions from one finite element model to another. Initial calculations using high-resolution near-shore topography in the eastern Irish Sea and including “wetting and drying” showed that M₂ tidal amplitudes and phases in the region computed with TELEMAC were in good agreement with observations. The ADCIRC code gave amplitudes about 30 cm lower and phases about 8° higher. For the M₄ tide, in the eastern Irish Sea amplitudes computed with TELEMAC were about 4 cm higher than ADCIRC on average, with phase differences of order 5°. For the M₆ component, amplitudes and phases showed significant small-scale variability in the eastern Irish Sea, and no clear bias between the models could be found. Although setting a minimum water depth of 5 m in the near-shore region, hence removing wetting and drying, reduced the small-scale variability in the models, the differences in M₂ and M₄ tide between models remained. For M₆, a significant reduction in variability occurred in the eastern Irish Sea when a minimum 5-m water depth was specified. In this case, TELEMAC gave amplitudes that were 1 cm higher and phases 30° lower than ADCIRC on average. For QUODDY in the eastern Irish Sea, average M₂ tidal amplitudes were about 10 cm higher and phase 8° higher than those computed with TELEMAC. For M₄, amplitudes were approximately 2 cm higher with phases of order 15° higher in the northern part of the region and 15° lower in the southern part. For M₆ in the north of the region, amplitudes were 2 cm higher and about 2 cm lower in the south. Very rapid M₆ tidal-phase changes occurred in the near-shore regions. The lessons learned from this model inter-comparison study are summarised in the final section of the paper. In addition, the problems of performing a detailed model–model inter-comparison are discussed, as are the enormous difficulties of conducting a true model skill assessment that would require detailed measurements of tidal boundary forcing, near-shore topography and precise knowledge of bed types and bed forms. Such data are at present not available.     

Patterns of pH change inferred from chrysophyte microfossils in Adirondack and northern New England lakes

The silica-scaled chrysophytes are being used increasingly in paleolimnological studies of lake acidification. This study compares the distribution of scaled chrysophytes in the surface sediment of a calibration set of lakes from the Adirondack Park and from northern New England. With the exception of two taxa, species distributions with respect to pH are similar in the two regions. The stratigraphic distribution of chrysophytes is then compared in five Adirondack and three New England lakes. All the presently acid lakes indicate recent lake acidification, with more acid tolerant chrysophytes increasing in relative frequency in the recent sediments. Because the timing of these species changes is unrelated to any watershed disturbances that may have occurred in the lakes' catchments, atmospheric deposition of acids is believed responsible for these recent pH changes.This is the seventh of a series of papers to be published by this journal which is a contribution of the Paleoecological Investigation of Recent Lake Acidification (PIRLA) project. Drs. D. F. Charles and D. R. Whitehead are guest editors for this series.     

A robust algorithm for ellipse-based discrete element modelling of granular materials

With recent research indicating the importance of the rolling mechanism of deformation in granular systems consisting of perfectly round particles, it has become popular to use ellipse-shaped particles in the Discrete Element Method (DEM) numerical model. Inherent in this technique is the need for accurately computing ellipse to ellipse intersection, in order to properly detect contact formation and compute relative contact velocities. However, the commonly used algorithms for computing ellipse-ellipse intersection are generally poorly conditioned and can be inaccurate. An alternate method for computing ellipse-ellipse intersection is developed and presented which results in a well-conditioned, stable and accurate contact detection method. These modification are incorporated into the general DEM algorithm.     

Enigmatic Efflorescence in Kärkevagge, Swedish Lapland: The Key to Chemical Weathering?

The second marvel to catch the eye of the visitor to Kärkevagge, after the impressive boulder deposit on the floor of the valley, is the series of prominent white stripes running down the valley's dark cliffs. Streams and springs descending the eastern flank of Kärkevagge are marked by the presence of whitish coatings on the black rock surfaces and on cobbles lining ephemeral waterways. These were referred to as 'lime crusts' by early investigators, but they are not reactive to HCl. We believe that they are a precipitate resulting from acid attack on the local rocks. Pyrite is common in many of the rocks in the valley and its oxidation produces sulfuric acid. As the dissolved mineral elements are carried in the drainage water, efflorescence forms on the surfaces where the water flows due to evaporation or to changes in temperature. The exact mineralogy of the white crusts is unknown, but the crusts are dominated by Al, S, and O, and in some cases by Ca, depending on the substrate and local conditions. Gypsum, illite, and chlorite have been identified by X–ray diffraction of some scrapings of white–coated rocks. However, we believe that some unidentified oxy–hydroxy aluminum sulfates make up the bulk of the precipitates.