TERRESTRIAL HEAT FLOW IN THE SWISS ALPS

i
Terrestrial heat flow has been measured in three Alpine railroad tunnels. The geothermal gradients were calculated from temperatures measured during the construction of the tunnels, and corrections for topographic irregularities were made. The thermal conductivity of 113 rock specimens from the vicinity of the tunnels was measured. The heat flow in the Gotthard tunnel was found to be 1.6 10^-6 cal/cm² sec, in the Simplon 2.2 10^-6 cal/cm² sec, and in the Loetschberg 1.9 10^-6 cal/cm² sec. Most of the flux at the surface can be attributed to radioactive decay in a thickened crust, but a non-uniform distribution of radioactive elements may be required to explain the relatively high heat flow in the Simplon and Loetschberg tunnels.     

THERMAL CONVECTION IN THE INTERIOR OF THE EARTH

The paper deals with thermal convection in the shell of the earth, caused by various assumed zonal temperature perturbations. One temperature perturbation here treated is that due to the difference in the temperature distribution under a continental crust made up of 10 km. of granite on top of 20 km. of basaltic material and a sub-oceanic crust consisting of 25 km. of basalt. The kinematic viscosity v was assumed to be 3 × 1021, as estimated recently by N. A. Haskell from a study of the uplift of Fennoscandia after the ice load. It is found that when g, ν, p and a, the coefficient of volume expansion, are constant, the velocities and their gradients are proportional to the amplitude of the temperature perturbation and to (ag/ν), while the stresses are independent of the viscosity. The velocities are found to be of the order of I cm./year. The shearing stress exerted by the convective substratum on the crust is of the order of 10⁷ dyn./cm.², while the normal stresses are about 10 times larger. The crust is pushed upwards under the warmer (continental) regions and pulled downwards under colder (oceanic) regions. The maximum stress-difference occurs at the bottom of the crust over the centre of the oceans or continents. The surface inequalities are nearly compensated.     

THE INFLUENCE OF TEMPERATURE GRADIENT AND VARIATION OF COMPOSITION IN THE MANTLE ON THE COMPUTATION OF DENSITY VALUES IN EARTH MODEL A

i
The effect of errors in the expression for dp/dr used in constructing the Earth Model A is examined. It is shown that, for a 10 per cent error in dp/dr arising from possible deviations from adiabatic gradients and uniform chemical composition in the regions B and D, the maximum error in the density ρ in the mantle of Model A would be only 0.07 g/cm³. This is less than one-quarter of the error that would be entailed if a 10 per cent correction to dp/dr were applied throughout the mantle. Thus temperature deviations to any likely extent cannot be considered as serious sources of error in the density distribution of Model A.     

Estimating the crust permeability from fluid-injection-induced seismic emission at the KTB site

During the hydraulic-fracturing experiment in the German Continental Deep Drilling Borehole (KTB) in December 1994, microseismic activity was induced. Here we develop a technique for estimating permeability using the spatio-temporal distribution of the fluid-injection-induced seismic emission. The values we have obtained for the KTB experiment (0.25times10^-16 to 1.0times10^-16 1.0times10^-16 m²) are in a very good agreement with the previous hydraulic-type permeability estimates from KTB deep-observatory studies. In addition, our estimates of the hydraulic diffusivity support the previously calculated value for the upper crust, which is of the order of 1 m² s^-1. However, this estimate now relates to the depth range 7.5-9 km.     

The tidal dynamics of the Irish and Celtic Seas

Summary. Current meter data collected over periods of more than 14 day from the Irish and Celtic Seas are harmonically analysed and presented in maps of tidal stream information. Making use of the analysed current data, and by constructing time series of frictional and inertial stresses which are also harmonically analysed, harmonic constituents of the surface tidal slopes at current meter stations are obtained. Using these with data collected from offshore tide gauges, and in conjunction with coastal tide data, cotidal maps are drawn with some confidence for M ₂, S ₂, O ₁ and K ₁, the M ₂ chart resolving the discrepancy which exists between the different charts of the Celtic Sea already produced. Cotidal maps for M ₃ and M ₄ are also presented.
The mean over a tidal cycle of the energy flux for M ₂, S ₂ and O ₁ is also presented in the form of the total energy flux in these constituents which crosses different sectional lines. A flux of 44 × 10⁶ kW is observed to enter the Celtic Sea from the Atlantic and this is compared with previous estimates. An energy budget is also performed for M ₂, including all the effects of astronomical forcing and Earth tides to enable comparison to be made between the true energy inflow and the estimated frictional dissipation. Finally, comparison is made between the mean of the instantaneous energy flux and the sum of the energy fluxes associated with the major harmonics.     

Terrestrial Heat Flow and Heat Generation in South-west England

Three heat flow values for south-west England are presented. Two of the sites, Geevor and South Crofty, are operating tin mines on the northern contacts of the Land's End and Carnmenellis Granites, respectively, while the third, Wilsey Down, is a stratigraphical borehole 5 km north of the Bodmin Moor Granite. After applying topographic corrections values of 128·6 mW m^-2 (3·07 μ cal cm^-2 s^-1) for Geevor, 128·9 mW m^-2 (3·08 (A cal cm^-2 s^-1) for South Crofty and 67·3 mW cm^-2 (1·61 cal cm^-2 s^-1) for Wilsey Down, were determined. The value at Wilsey Down is shown to be consistent with that for an environment in which the Hercynian orogeny was the last significant thermal event. An additional heat source term must clearly be involved at Geevor and South Crofty to account for the unusually high values at these sites. Radiogenic heat production has been determined on granites from these sites and in spite of the fact that it is high it does not fully account for the measured heat flow. A compilation of underground temperature measurements made in the nineteenth century suggests that high heat flow is a general feature of the mineralized belt. At least part of this can be explained in terms of hot spring activity recorded widely throughout the area but the ultimate cause remains to be evaluated.     

Glaciers in Svalbard: mass balance, runoff and freshwater flux

Gain or loss of the freshwater stored in Svalbard glaciers has both global implications for sea level and, on a more local scale, impacts upon the hydrology of rivers and the freshwater flux to fjords. This paper gives an overview of the potential runoff from the Svalbard glaciers. The freshwater flux from basins of different scales is quantified. In small basins (A < 10 km²), the extra runoff due to the negative mass balance of the glaciers is related to the proportion of glacier cover and can at present yield more than 20% higher runoff than if the glaciers were in equilibrium with the present climate. This does not apply generally to the ice masses of Svalbard, which are mostly much closer to being in balance. The total surface runoff from Svalbard glaciers due to melting of snow and ice is roughly 25 ± 5 km³ a⁻¹, which corresponds to a specific runoff of 680 ± 140 mm a⁻¹, only slightly more than the annual snow accumulation. Calving of icebergs from Svalbard glaciers currently contributes significantly to the freshwater flux and is estimated to be 4 ± 1 km³ a⁻¹ or about 110 mm a⁻¹.     

Oxygen consumption in Macrotrachela musculosa and Trichotria truncata (Rotatoria) from the High Arctic

Oxygen consumption by rotifers Macrotrachela musculosa and Trichotria truncata from Spitsbergen tundra (77°N) was measured using the method of Cartesian divers. The metabolic rate of M. musculosa was: 0.205 10⁻³mm³ 0₂ per g 10⁻⁶ per hour at 2°C, 0.201 10⁻⁶mm³ at 6°C and 0.616 10⁻³mm³ 0₂ per g 10⁻⁶ per hour at 10°C. The metabolic rate of Trichotria truncata at 6° was 0.103 10⁻³mm³ per g 10⁻⁶ per hour. The relation between body weight and oxygen consumption by M. musculosa at 2°C is expressed with the equation R = 0.18W^0.67, with R – oxygen consumption in mm³10⁻³ per individual per hour and W – wet weight of an animal in g 10⁻⁶.     

Deglaciation effects on the rotation of the Earth

Summary. The viscoelastic response of the Earth to the mass displacements caused by late Pleistocene deglaciation and concomitant sea level changes is shown to be capable of producing the secular motion of the Earth's rotation pole as deduced from astronomical observations. The calculations for a viscoelastic Earth yield a secular motion in the direction of 72° W meridian which is in excellent agreement with observed values. The average Newtonian viscosity and the relaxation time obtained from polar motion data are about (1.1 ± 0.6)10²³ poise (P) and 10⁴ (1 ± 0.5) yr. The non-tidal secular acceleration of the Earth can also be attributed to the viscoelastic response to deglaciation and results in an independent viscosity estimate of 1.6 × 10²³ P with upper and lower limits of 1.1 × 10²³ and 2.8 × 10²³ P. These values are in agreement with those based on the polar drift analysis and indicate an average mantle viscosity of 1–2 × 10²³ P.     

Precise measurement of seismic travel times from an air gun source

Summary. Seismic travel times for extrema, zero-crossings, or entire body waves need to be determined precisely to one part in 10³ or better in several varieties of seismic studies employing an impulsive artificial source. Examples are crosshole surveys which delineate rock crack distribution separating the holes and monitoring of crustal seismic travel times in earthquake precursor studies. A timing resolution of one part in 10³ has been achieved previously using digitally recorded seismic data. These methods, however, do not use interpolation between digitized data points as a method to increase the timing resolution. We report travel-time determinations based on interpolation between digitized points which achieve a precision of two parts in 10⁴, a five-fold improvement over the existing methods. In addition, the effects of seismic noise on travel-time measurement have been compared for the extremum location, the unnormalized correlation, and the normalized correlation method. The following conclusions are drawn from this comparison: (1) the normalized correlation method provides an 18–55 per cent improvement in the standard deviation of the mean over the extremum location method, and (2) results as accurate as those by the normalized correlation can be obtained by the unnormalized correlation if a complete up-and-down swing of the waveform is used as the master trace and if the master trace is close to being sinusoidal. The advantage of the unnormalized correlation over the normalized correlation is speed; the unnormalized correlation is faster by a factor of 28 in computing time.     

Distribution of Chemical Constituents in Superimposed Ice from Austre Brøggerbreen, Spitsbergen

10 m and 2.3 m ice cores were obtained on Austre Brøggerbreen, Spitsbergen in Svalbard (78°53 ' N, 11°56 ' E, 450 m a.s.l.) in September 1994 and in March 1995, respectively. Stratigraphy, bulk density, pH, electrical conductivity, and major ions were obtained from the core samples.
The chemical effect of meltwater percolation through snow/ice is examined. Good correlation between Cl⁻ and Na⁺ was obtained. The ratio of Cl⁻ to Na⁺ was 1.14 which was nearly the same value as in bulk sea water. However, the variation of Cl⁻/Na⁺ shows that higher ratio occured in the bubble-free ice. Furthermore the Cl⁻ ions remain in higher concentration than SO ₄ ²⁻ or Na⁺ ions.     

STUDIES ON THE GROWTH OF RHIZOCARPON GEOGRAPHICUM IN NW SCOTLAND, AND SOME IMPLICATIONS FOR LICHENOMETRY

Scotland, a maritime subpolar environment (55–60°N), has seen relatively few applications of lichenometry – even though it offers much potential. Perhaps surprisingly, direct measurements of Rhizocarpon geographicum growth rates in Scotland are so far lacking. This study reports on the growth of this crustose areolate species from two sites in Assynt, NW Scotland, between 2002 and 2009. Repeat photography of 23 non-competing thalli growing under identical environmental conditions on a single vertical surface over 5 years at Inchnadamph showed growth rates to be a function of size – with larger thalli (10–30 mm) growing significantly faster than the smallest thalli (<10 mm). Mean diametral growth rates in thalli >10 mm are 0.67 mm yr⁻¹ (s.d. = 0.16). Studies on a second vertical surface near Lochinver, over 7 years, yielded complex growth data on a more mature population of R. geographicum thalli (<50 mm in diameter). Here, mean diametral growth rates in the larger thalli (>10 mm) are slower (0.29 mm yr⁻¹; s.d. = 0.12) than those at Inchnadamph. However, at this site, competition with other species rules out any meaningful comparison of growth rates between the two sites. Other growth processes were monitored over the five to seven-year study period, including hypothallus growth, areolae development, thallus coalescence, and inter-species competition – all have important implications for the use of Rhizocarpon species in lichenometry.     

Velocity structure of upper-mantle transition zones beneath central eurasia from seismic inversion using genetic algorithms

We present velocity constraints for the upper-mantle transition zones beneath Central Siberia based on observations of the 1982 RIFT Deep Seismic Sounding (DSS) profile. The data consist of seismic recordings of a nuclear explosion in north-western Siberia along a 2600 km long seismic profile extending from the Yamal Peninsula to Lake Baikal. We invert seismic data from the mantle transition zones using a non-linear inversion scheme using a genetic algorithm for optimization and the WKBJ method to compute the synthetic seismograms. A statistical error analysis using a graph-binning technique was performed to provide uncertainty values in the velocity models.
Our best model for the upper-mantle velocity discontinuity near 410 km depth has a two-stage velocity-gradient structure, with velocities increasing from 8.70–9.25 km s⁻¹ over a depth range of 400–415 km, a gradient of 0.0433 s⁻¹, and from 9.25–9.60 km s⁻¹ over a depth range of 415–435 km, a gradient of 0.0175 s⁻¹. This derived model is consistent with other seismological observations and mineral-physics models. The model for the velocity discontinuity near 660 km depth is simple, sharp and includes velocities increasing from 10.15 km s⁻¹ at 655 km depth to 10.70 km s⁻¹ at 660 km depth, a gradient of 0.055 s⁻¹.     

DEGLACIATION OF THE EASTERN FLANK OF THE NORTH PATAGONIAN ICEFIELD AND ASSOCIATED CONTINENTAL-SCALE LAKE DIVERSIONS

ABSTRACT. We examine the deglaciation of the eastern flank of the North Patagonian Icefield between latitudes 46° and 48°S in an attempt to link the chronology of the Last Glacial Maximum moraines and those close to present-day outlet glaciers. The main features of the area are three shorelines created by ice-dammed lakes that drained eastwards to the Atlantic. On the basis of 16 ¹⁴C and exposure age dates we conclude that there was rapid glacier retreat at 15–16 ka (calendar ages) that saw glaciers retreat 90–125 km to within 20 km of their present margins. There followed a phase of glacier and lake stability at 13.6–12.8 ka. The final stage of deglaciation occurred at c. 12.8 ka, a time when the lake suddenly drained, discharging nearly 2000 km³ to the Pacific Ocean. This latter event marks the final separation of the North and South Patagonian Icefields. The timing of the onset of deglaciation and its stepped nature are similar to elsewhere in Patagonia and the northern hemisphere. However, the phase of lake stability, coinciding with the Antarctic Cold Reversal and ending during the Younger Dryas interval, mirrors climatic trends as recorded in Antarctic ice cores. The implication is that late-glacial changes in southern Patagonia were under the influence of the Antarctic realm and out of phase with those of the northern hemisphere.     

Nitrogen uptake rates in phytoplankton and ice algae in the Barents Sea

Uptake rates of NH₄⁺, NO₃⁻ and dissolved organic nitrogen (urea) were measured in phytoplankton and in ice algae in the Barents Sea using a ¹⁵N-technique. NO₃⁻ was the most important nitrogen source for the ice algae (f-ratio = 0.92). The in situ irradiances in the subsurface chlorophyll maximum and in the ice algal communities were low. The in situ NO₃⁻ uptake rate in the ice algal communities was light-limited The in situ NO₃⁻ and NH₄⁻ uptake rates in the subsurface chlorophyll maximum were at times light-limited. It is hypothesised that NH₄⁺ may accumulate in low light in the bottom of the euphotic zone and inhibit the in situ NO₃⁻ uptake rate.     

The traveltime perturbations for seismic body waves in factorized anisotropic inhomogeneous media

The traveltime perturbation equations for the quasi-compressional and the two quasi-shear waves propagating in a factorized anisotropic inhomogeneous (FAI) media are derived. The concept of FAI media simplifies considerably these equations. In the FAI medium, the density normalized elastic parameters a _ijkl ( X _i ) can be described by the relation a _ijkl ( X _i) = f ²( x _i ) A _ijkl, where A _ijkl are constants, independent of coordinates x _i and f ²( x _i) is a continuous smooth function of x _i . The types of anisotropy ( A _ijkl ) and inhomogeneity [ f ( x _i)] are not restricted. The traveltime perturbations of individual seismic body waves ( q ^P , qS 1 and qS 2) propagating in the FAI medium depend, of course, both on the structural pertubations [δ f ²( x _i)] and on the anisotropy perturbations (δ A _ijkl ), but both these effects are fully separated. The perturbation equations for the time delay between the two qS -waves propagating in the FAI medium are simplified even more. If the unperturbed (background) medium is isotropic, the perturbation of the time delay does not depend on the structural perturbations (δ f ²( x _i) at all. This striking result, valid of course only in the framework of first-order perturbation theory, will simplify considerably the interpretation of the time delay between the two split qS -waves in inhomogeneous anisotropic media. Numerical examples are presented.     

Crustal and upper mantle structure of the northwestern North Island, New Zealand, from seismic refraction data

The crustal and upper mantle structure of the northwestern North Island of New Zealand is derived from the results of a seismic refraction experiment; shots were fired at the ends and middle of a 575 km-long line extending from Lake Taupo to Cape Reinga. The principal finding from the experiment is that the crust is 25 ± 2 km thick, and is underlain by what is interpreted to be an upper mantle of seismic velocity 7.6 ± 0.1 km s⁻¹, that increases to 7.9 km s⁻¹ at a depth of about 45 km. Crustal seismic velocities vary between 5.3 and 6.36 km s⁻¹ with an average value of 6.04 km s⁻¹. There are close geophysical and geological similarities between the north-western North Island of New Zealand and the Basin and Range province of the western United States. In particular, the conditions of low upper-mantle seismic velocities, thin crust with respect to surface elevation, and high heat-flow (70–100 mW m⁻²) observed in these two areas can be ascribed to their respective positions behind an active convergent margin for about the past 20 Myr.     

Palaeoclimatic record from a loess–soil profile in northeastern Bulgaria—II. Correlation with global climatic events during the Pleistocene

Magnetic susceptibility ( χ ) variations and the behaviour of the ratio of susceptibility to saturation magnetization ( χ/J _s ) along the loess/palaeosol section at Koriten (NE Bulgaria) are used to deduce climatic changes during the Pleistocene in southeastern Europe. A good correlation of susceptibility variations with the astronomically tuned oxygen isotope record from ODP site 677 enables us to propose a more precise dating of the upper part of the Bulgarian loess complex. Close correspondence between susceptibility and δ ^{1 8} O records demonstrates the global significance of the palaeoclimatic signal recorded, although differences in relative amplitudes of χ and χ/J _s and δ ^{1 8} O create difficulties in making quantitative estimates of the climatic humidity in the past. The role of local factors affecting the palaeoclimatic mineral magnetic record deduced from the profile studied in Bulgaria is discussed.     

Dissolved and total organic and inorganic carbon in some British rivers

Rivers transport both organic and inorganic carbon from their sources to the sea. Results of ~800 organic and inorganic analyses from various British rivers of contrasting size and land use are presented here: (1) the headwater River Tern, a rural river of 852 km² catchment; (2) the Ouseburn, a small urban 55 km² catchment; (3) the River Tyne, a larger river system of ~3000 km² catchment; (4) a spatial survey from 205 sample sites on ~60 rivers from SW England. We found that, with the exception of peat-rich headwaters, DIC concentration is always greater than DOC. DIC is primarily in the form HCO₃ ⁻ , with DIC concentrations highest in highly urbanised catchments, typically greater than those observed in catchments with carbonate bedrock, demonstrating a significant and previously unrecognised anthropogenic inorganic carbon input to urban rivers.     

Long-period geomagnetic variations and mantle conductivity: an inversion using Bailey's method

Summary. We have implemented an algorithm which is based on Bailey's solution of the inverse problem of electromagnetic induction in the Earth. The study was motivated by recent determinations of very long period data and also benefited from recent redeterminations of high frequency data. The algorithm has been successfully tested to provide reliable estimates of conductivity down to a depth of 2000 km, using synthetic data in the period range from 4 days to 11 years. Smooth data sets, which are required for the inversion, were constructed from various sources. At a given depth, the range of inverted models is less than one order of magnitude. Due to the lack of high frequency data, the conductivity of the upper 600 km of the mantle, which is found to be of the order of 10⁻¹Ω⁻¹ m⁻¹, may be overestimated. The algorithm performs well in the middle mantle, where conductivity rises steadily from 1 to 50 Ω⁻¹ m⁻¹. The lack of very low frequency data and limitations of the algorithm prevent one from obtaining meaningful estimates in the lower mantle. However, the study of the propagation of the late 1960s secular variation acceleration provides an estimate of the mean conductivity of the whole mantle. Thus, a complete mantle profile can be constructed. It is found that deep mantle conductivity probability does not exceed a few hundred Ω⁻¹ m⁻¹.