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.     

Influence of wind upon water quality in Liverpool Bay

Comprehensive distribution patterns of physical and chemical characteristics have been obtained from a series of cruises in Liverpool Bay. The marked feature of these distributions was their temporal variability, suggesting that the surface residual circulation is also temporally variable. The influence of wind stress upon the circulation pattern and hence water quality of this sea area is illustrated.     

Rock salt mass in the Paleoproterozoic sequence of the Onega trough in Karelia (from the Onega parametric well data)

A thick (200 m) rock salt mass covering Late Archean granitoids was exposed for the first time in the Early Proterozoic volcanogenic-sedimentary formations in the Onega trough of the east Baltic Shield by the Onega parametric well. The mineral composition of salts, their geochemical features, and the isotopic composition of carbonate carbon and oxygen have been studied. After fluid inclusions present in salts, their metamorphism temperature and isotopic composition of helium and argon were determined. The obtained results give evidence of the fact that rock salts and magnesites associated with them were formed in an evaporate basin with participation of deep crustal processes. The age of the underlying granitoids (2.716 ± 9 Ma) is determined using the Pb—Pb method.     

Developing a modern pollen–climate calibration data set for Norway

Bjune, A. E., Birks, H. J. B., Peglar, S. M. & Odland, A. 2010: Developing a modern pollen–climate calibration data set for Norway. Boreas, Vol. 39, pp. 674–688. 10.1111/j.1502‐3885.2010.00158.x. ISSN 0300‐9483. Modern pollen–climate data sets consisting of modern pollen assemblages and modern climate data (mean July temperature and mean annual precipitation) have been developed for Norway based on 191 lakes and 321 lakes. The original 191‐lake data set was designed to optimize the distribution of the lakes sampled along the mean July temperature gradient, thereby fulfilling one of the most critical assumptions of weighted‐averaging regression and calibration and its relative, weighted‐averaging partial least‐squares regression. A further 130 surface samples of comparable taphonomy, taxonomic detail and analyst became available as a result of other projects. These 130 samples, all from new lakes, were added to the 191‐lake data set to create the 321‐lake data set. The collection and construction of these data sets are outlined. Numerical analyses involving generalized linear modelling, constrained ordination techniques, weighted‐averaging partial least‐squares regression, and two different cross‐validation procedures are used to asses the effects of increasing the size of the calibration data set from 191 to 321 lakes. The two data sets are used to reconstruct mean July temperature and mean annual precipitation for a Holocene site in northwest Norway and a Lateglacial site in west‐central Norway. Overall, little is to be gained by increasing the modern data set beyond about 200 lakes in terms of modern model performance statistics, but the down‐core reconstructions show less between‐sample variability and are thus potentially more plausible and realistic when based on the 321‐lake data set.     

The Karoo Basin of South Africa: type basin for the coal-bearing deposits of southern Africa

The coal-bearing sediments and coal seams of the Karoo Basin, Southern Africa are described and discussed. The Karoo Basin is bounded on its southern margin by the Cape Fold Belt, onlaps onto the Kaapvaal Craton in the north and is classified as a foreland basin. Coal seams are present within the Early Permian Vryheid Formation and the Triassic Molteno Formation.The peats of the Vryheid Formation accumulated within swamps in a cool temperate climatic regime. Lower and upper delta plain, back-barrier and fluvial environments were associated with peat formation. Thick, laterally extensive coal seams have preferentially accumulated in fluvial environments. The coals are in general inertinite-rich and high in ash. However, increasing vitrinite and decreasing ash contents within seams occur from west to east across the coalfields. The Triassic Molteno coal seams accumulated with aerially restricted swamps in fluvial environments. These Molteno coals are thin, laterally impersistent, vitrinite-rich and shaly, and formed under a warm temperate climatic regime.Palaeoclimate, depositional systems, differential subsidence and basin tectonics influence to varying degrees, the maceral content, thickness and lateral extent of coal seams. However, the geographic position of peat-forming swamps within a foreland basin, coupled with basin tectonics and differential subsidence are envisaged as the primary controls on coal parameters. The Permian coals are situated in proximal positions on the passive margin of the foreland basin. Here, subsidence was limited which enhanced oxidation of organic matter and hence the formation of inertinitic coals. The coals in this tectonic setting are thick and laterally extensive. The Triassci coals are situated within the tectonically active foreland basin margin. Rapid subsidence and sedimentation rates occurred during peat formation which resulted in the preservation of thin, laterally impersistent, high ash, vitrinite-rich, shaly coals.     

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.     

Meteorological excitation of the annual polar motion

Summary. Numerous studies have indicated that the annual term in the polar motion cannot be explained in any detail by meteorological/hydrological excitation and no reasonable alternative excitations have been put forward. Part of the problem has been that the hydrostatic adjustment of the oceans to the atmospheric pressure changes has traditionally been computed using the inverse barometer approach. This approach does not properly model the gravitational interaction between the atmosphere and oceans, and the inverse barometer theory is modified in this paper to account for this properly. The information necessary to compute the ocean tide and polar excitation caused by any change in the atmospheric pressure pattern is presented. The results of the application of this theory to two global atmospheric pressure data sets are examined and compared to results of other workers.
It is concluded that the atmosphere is observed well enough to answer the question of the annual excitation of polar motion and it is argued that the ground water excitation is the component with the largest error and remains the chief obstacle to the successful solution of this problem.