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.     

Amplitudes of long-period PcP and the core-mantle boundary

The amplitudes of the core reflection PcP are sensitive to the wave velocities and densities in the neighborhood of the core-mantle boundary (CMB). We study the amplitude ratio of the long-period phases PcP and P from two South American deep-focus earthquakes with favorable fault-plane solution, depth and magnitude, as recorded by WWNSS and CSN stations in North America.Comparison is made with long-period PcP/P amplitude ratios, derived from theoretical seismograms for a variety of CMB models. Models from previous studies, which were mainly derived from short-period PcP observations and which are characterized by discrete layers above the CMB, are almost all inconsistent with the long-period data. The data also discriminate against low nonzero S velocities below the CMB. Simple first-order-discontinuity models of the CMB, for instance according to the Jeffreys-Bullen earth model or according to recent models based mainly on free oscillations, explain the data reasonably well.Model improvements are attempted by varying the P-velocity gradient above the CMB. The best amplitude fit is obtained for a rather strong decrease in P velocity with depth in this zone which, however, gives no acceptable traveltime fit for PcP. The scatter in body-wave amplitudes is considerable even for long-period waves and may prevent the correct assessment of that part of the amplitude variation of a phase with distance that is due to the variation of velocities and densities with depth alone.     

N and O isotope effects during nitrate assimilation by unicellular prokaryotic and eukaryotic plankton cultures

In order to provide biological systematics from which to interpret nitrogen (N) and oxygen (O) isotope ratios of nitrate (¹⁵N/¹⁴N, ¹⁸O/¹⁶O, respectively) in the environment, we previously investigated the isotopic fractionation of nitrate during its assimilation by mono-cultures of eukaryotic algae (Granger et al., 2004). In this study, we extended our analysis to investigate nitrate assimilation by strains of prokaryotic plankton. We measured the N and O isotope effects, ¹⁵ε and ¹⁸ε, during nitrate consumption by cultures of prokaryotic strains and by additional eukaryotic phytoplankton strains (where ε is the ratio of reaction rate constants of the light vs. heavy isotopologues, ^lightk and ^heavyk; ε = ^lightk/^heavyk − 1 × 1000, expressed in per mil). The observed ¹⁵ε ranged from 5‰ to 8‰ among eukaryotes, whereas it did not exceed 5‰ for three cyanobacterial strains, and was as low as 0.4‰ for a heterotrophic α-protoeobacterium. Eukaryotic phytoplankton fractionated the N and O isotopes of nitrate to the same extent (i.e., ¹⁸ε ∼ ¹⁵ε). The ¹⁸ε:¹⁵ε among the cyanobacteria was also ∼1, whereas the heterotrophic α-proteobacterial strain, which showed the lowest ¹⁵ε, between 0.4‰ and 1‰, had a distinct ¹⁸ε:¹⁵ε of ∼2, unlike any plankton strain observed previously. Equivalent N vs. O isotope discrimination is thought to occur during internal nitrate reduction by nitrate reductase, such that the cellular efflux of the fractionated nitrate into the medium drives the typically observed ¹⁸ε:¹⁵ε of ∼1. We hypothesize that the higher in the ¹⁸ε:¹⁵ε of the α-proteobacterium may result from isotope discrimination by nitrate transport, which is evident only at low amplitude of ε. These observations warrant investigating whether heterotrophic bacterial assimilation of nitrate decreases the community isotope effects at the surface ocean.     

Facies Structure and Quantitative Parameters of Pleistocene Sediments of the Bering Sea

Lithofacies zoning is described for the first time for the Neo- and Eopleistocene of the Bering Sea. Four lithofacies sedimentation zones are distinguished: (I) terrigenous; (II) siliceous–terrigenous; (III) siliceous, and (IV) volcanoterrigenous ones. Corresponding maps were treated using Ronov volumetric method to quantify sedimentation parameters for distinguished lithofacies zones (subzones) and types of Pleistocene sediments. It was revealed that terrigenous sediments predominate over other sediments. Accumulation of the terrigenous sediments was more intense (by 1.4 times) in the Neopleistocene than in the Eopleistocene. The sedimentation rate of siliceous sediments of the Bowers Ridge in the Eopleistocene was two times higher than in the Neopleistocene.     

Rockburst research at Falconbridge's Strathcona Mine,Sudbury, Canada

Microseismic systems at five mines in the Sudbury Basin provide the basic data for Falconbridge Limited's rockburst research. Daily and long-term analysis of this data as well as underground observations have confirmed the fault-slip mechanism at three mines. A detailed analysis of the complete history of Falconbridge Mine is being conducted and Distinct Element numerical models are being used to simulate both the stick-slip behaviour of faults and the dynamic effects of the induced vibrations on rock and backfill. University research includes acoustic tomographic imaging of the rock mass based on seismic wave propagation and collection of full microseismic waveforms to allow application of advanced seismic and statistical analysis techniques.Formerly with Mines Technical Services, Falconbridge Limited, Sudbury Operations, Falconbridge , Ont., P0M 1S0.Presented at the Fred Leighton Memorial Workshop on Mining Induced Seismicity, Montreal, Canada, August 30, 1987.     

Multidimensional statistical technique for detection of low contrast geochemical anomalies

In this article we propose an advanced technique for detecting low contrast geochemical anomalies using a set of features. There are three principal elements in this technique: (1) a statistical measure of the contrast of the anomaly, denoted as τ; (2) selection of a background population; and (3) reduction of the dimensionality of the feature space. In the frame of the model, which describes the statistical distribution of geochemical background as a multidimensional normal distribution of logarithms of concentrations, the index, τ, is a powerful test statistic for the hypothesis of abnormality of an observation. Maps of τ anomalies can be rigorously interpreted on the basis of statistical inferences. Under all equal conditions this technique allows the detection of geochemical anomalies with at least the same contrast (if the chemical elements in a background population are correlated, then even the better) as using selective extractions of metals from soil or other techniques for data processing. The advantages of the proposed technique are demonstrated both theoretically and on examples of rare-metal and copper–nickel mineral deposits.     

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.