Thermochronometry and palaeomagnetism of the Archaean Nelshoogte Pluton, South Africa

²⁰⁷Pb/²⁰⁶Pb single-grain zircon, ⁴⁰Ar/³⁹Ar single-grain hornblende and biotite, and ⁴⁰Ar/³⁹Ar bulk-sample muscovite and biotite ages from the Nelshoogte trondhjemite pluton located in eastern Transvaal, South Africa, show that this granitoid had a protracted thermal history spanning 3213±4  Ma to about 3000  Ma. Whole-rock ⁴⁰Ar/³⁹Ar ages from cross-cutting dolerite dykes indicate that these were intruded at about 1900  Ma. There is no evidence of this or other, later events significantly affecting the argon systematics of the minerals from the pluton dated by the ⁴⁰Ar/³⁹Ar method.
  The pluton has a well-defined palaeomagnetic pole which is dated at 3179±18 (2 σ ) Ma by ⁴⁰Ar/³⁹Ar dating of hornblende. This pole (18°N, 310°E, A ₉₅=9°) yields a palaeolatitude of 0°, significantly different from other Archaean poles from the Kaapvaal Craton. The palaeolatitude difference implies that there was significant apparent polar wander during the Archaean. A second, overprinting magnetization seen in the pluton is also seen in the lower-Proterozoic dolerite dykes, and is consistent with other lower-Proterozoic (2150–1950  Ma) poles for southern Africa.     

222Rn and 226Ra: indicators of sea-ice effects on air-sea gas exchange

²²² Rn and ²²⁶Ra distributions beneath the sea ice of the Barents Sea revealed that ice cover has varied effects on air-sea gas exchange. Twice, once in late summer and once in late winter, seawater samples from the top meter below drill holes had ²²²Rn activities that were not lower than their ^226,Ra activities, indicating the existence of secular equilibrium and a negligible net exchange of ²²²Rn and other gases with the atmosphere. However, seawater in the upper 20-85 m usually exhibited at least some ²²²Rn depletion; ²²²Rn-to-²²⁶Ra activity ratios tended to have 'ice-free' values (0.3-0.9) in the summer and values between 0.9 and 1.0 in the winter. Integrated ²²²Rn depletions and piston velocities in both seasons typically fell in the lower 25% of the ranges for ice-free seawater, suggesting that a moderate but far from total reduction in gas exchange is normally caused by ice cover and/or meltwater. The results demonstrate that sea-ice interference with the oceanic uptake of atmospheric gases such as CO, is not well understood and needs further investigation.     

Deglaciation history and post-glacial mass movements in Balsfjord, northern Norway

The deglaciation history of Balsfjord, northern Norway, and post-glacial mass movement events were investigated. Radiocarbon dates indicate that the Balsfjord glacier retreated from the Tromsø–Lyngen moraines about 10.4 ¹⁴C Ky BP. Between ca. 10.3 ¹⁴C Ky BP and 9.9 ¹⁴C Ky BP, deposition of a distinct end moraine–the Skjevelnes moraine–in the central part of Balsfjord occurred. The transition from glacimarine to open marine sedimentary environment took place before 9.6 ¹⁴C Ky BP. Between ca. 9.5 ¹⁴C Ky BP and 8.4 ¹⁴C Ky BP, at least one local and three regional mass movement events occurred. After this period, no gravity flow activity is preserved in the cores. The high frequency of mass movements in the early post-glacial period is presumed to be due to fast sea level changes and/or tectonic activity induced by rapid isostatic uplift.     

40Ar—39Ar ages of guyots in the western Pacific and discussion of their evolution

Summary. Six submarine rocks dredged from guyots in the western Pacific were dated by means of ⁴⁰Ar–³⁹Ar step heating dating technique. All ages show mid- to late-Cretaceous ages in accordance with the generally supposed age of the ocean crust estimated from the magnetic anomaly pattern.
The back-tracked paths were calculated for the guyots with the use of their absolute ages, for both the Pacific plate movements proposed by Morgan and by Clague & Jarrard. In both cases these show birth places for the guyots, which are located near the equator.
Fossil ages are nearly concordant with the ⁴⁰Ar–³⁹Ar ages for the guyots which are older than 95 Myr, whereas the fossil ages are considerably younger than the ⁴⁰Ar–³⁹Ar ages for the guyots younger than 90 Myr. This contrast may be due to the world-wide marine transgression, whose climax is believed to have occurred at about 85 Myr BP.     

Micromonas pusilla (Prasinophyceae) as part of pico-and nanoplankton communities of the Barents Sea

Micromonas pusilla (Butcher) Manton & Parke appears to be a prominent member of the Barents Sea picoplankton community as revealed by the serial dilution culture method. Cell numbers frequently exceeded 10⁷ cells 1⁻¹, though they usually varied between 10³and 10⁶ cells l⁻¹. A number of other identified and unidentified taxa were recorded and quantified. Distribution relative to the marginal ice zone is reported.     

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.     

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.     

Palaeointensity results from Ethiopian basalts: implications for the Oligocene geomagnetic field strength

From a large collection of Ethiopian flood basalts (~30  Myr old) sampled for magnetostratigraphy, ⁴⁰Ar/³⁹Ar geochronology and geochemical investigations, 47 samples were selected in order to test their suitability for Thellier palaeointensity experiments. Only 17 samples from eight individual flows yielded reliable palaeointensity estimates, with flow-mean virtual dipole moments ranging from 3.0 to 10.5 × 10²²  A  m² .
  A critical review of the Oligocene palaeointensity data set, including these new Ethiopian data, indicates an Oligocene mean virtual dipole moment of 5.1 ± 2.5 × 10²²  A  m² for the complete data set. After applying mild selection criteria, the reduced data set yields a mean value of only 4.6 ± 1.9 × 10²²  A  m² . This value is significantly lower than the present-day field strength but is higher than the Mesozoic dipole low mean field. This low Oligocene field might be in agreement with the high palaeosecular variation and rather high non-dipole field invoked around 30  Ma. However, the Oligocene data set is largely dependent on the palaeointensity determinations from Armenia, obtained mainly from baked contacts, which show amazingly low dispersion at both flow and between-flow levels. More data are needed to reduce the weight of these determinations on the mean value and avoid a possible bias.     

Early Holocene environment on Bjørnøya (Svalbard) inferred from multidisciplinary lake sediment studies

Bjørnøya, a small (178 km²) island situated between the mainland of Norway and southern Spitsbergen, provides the opportunity for the reconstruction of early Holocene terrestrial and limnic palaeoenvironments in the southwestern Barents Sea. The AMS ¹⁴C dating technique, geochemical, mineral magnetic, micro and macrofossil analyses were applied to sediments recovered from lake Stevatnet and the results are interpreted in terms of palaeoenvironmental conditions between 9800 and 8300 ¹⁴C bp. After the disappearance of local glaciers before ca 9800¹⁴C BP, the lake productivity increased rapidly at the same time as pioneer plant communities developed on soils which gradually became more stable. Insect data indicates that strong seasonal contrasts with mean July temperatures around 9°C and mean January temperatures around −12°C prevailed between 9500 and 8300 ¹⁴C BP. These high summer temperatures, possibly as much as 4-5°C higher than the present, favoured the development of a flora including Dryas and Angelica cf. archangelica . The enhanced freeze/thaw processes led to an increased erosion of minerogenic and organic material. After 8000 ¹⁴C BP the temperatures may have gradually declined. The environmental reconstruction derived from our data set supports the conceptual insolation model which proposes maximum Holocene seasonality for the Northern Hemisphere at ca 9000 ¹⁴C BP.     

Strain accumulation along an oblique plate boundary: the Reykjanes Peninsula, southwest Iceland

We use annual GPS observations on the Reykjanes Peninsula (RP) from 2000 to 2006 to generate maps of surface velocities and strain rates across the active plate boundary. We find that the surface deformation on the RP is consistent with oblique plate boundary motion on a regional scale, although considerable temporal and spatial strain rate variations are observed within the plate boundary zone. A small, but consistent increase in eastward velocity is observed at several stations on the southern part of the peninsula, compared to the 1993–1998 time period. The 2000–2006 velocities can be modelled by approximating the plate boundary as a series of vertical dislocations with left-lateral motion and opening. For the RP plate boundary we estimate left-lateral motion  18⁺⁴₋₃ mm yr⁻¹  and opening of  7⁺³₋₂ mm yr⁻¹  below a locking depth of  7⁺¹₋₂ km  . The resulting deep motion of  20⁺⁴₋₃ mm yr⁻¹  in the direction of  N(100⁺⁸₋₆)°E  agrees well with the predicted relative North America–Eurasia rate. We calculate the areal and shear strain rates using velocities from two periods: 1993–1998 and 2000–2006. The deep motion along the plate boundary results in left-lateral shear strain rates, which are perturbed by shallow deformation due to the 1994–1998 inflation and elevated seismicity in the Hengill–Hrómundartindur volcanic system, geothermal fluid extraction at the Svartsengi power plant, and possibly earthquake activity on the central part of the peninsula.     

Approaches to Modelling the Surface Albedo of a High Arctic Glacier

Broadband surface albedo measurements, made during the summer melt season at three weather stations on John Evans Glacier (79°40 ' N, 74°00 ' W), varied strongly with the solar zenith angle ( θ _z ). Tests were carried out to assess the impact of diurnal variations in surface albedo on seasonal net shortwave radiation ( K ^* ) totals. Removing the diurnal signal from albedo measurements by daily averaging of hourly measurements, or by applying midday measurements to all hours of the day, changed K ^* by up to 16%. Ignoring measurements made at θ _z & 75°, to account for measurement (cosine) error at high θ _z , decreased K ^* by between 5 and 18%. Given the sensitivity of K ^* to diurnal patterns in surface albedo, experiments were carried out with two albedo models. One model accounted for albedo variations with θ _z and one did not. The model driven by θ _z , when implemented within a surface energy balance model for John Evans Glacier, produced better melt estimates. This suggests that diurnal variations in surface albedo should be accounted for in energy balance models of glacier melt.     

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.     

Seismic sources with observable glut moments of spatial degree two

Let ζ_Λ and r _Λ. be the hypocentral position and time of an extended indigenous seismic source. Backus showed that the force moment tensors of the source, Γ^{( m +1, n )}(ζ_Λ, r _Λ), determine and are determined by the motion which the source produces. For small m + n , only the long-period motion is relevant. The glut moment tensor Λ^{( m,n )} (ζ_Λ, r _Λ.) can be calculated uniquely from γ^{( m +1, n )}(ζ_Λ r _Λ) only if m = 0 or m = 1. The tensor G =Λ^(2,0)(ζ_Λ) gives the spatial variance tensor W_Λ of the source, and W_Λ. roughly describes the size, shape and orientation of the source region. Therefore the failure of the observed F =Γ^(3,0)(ζ_Λ) to determine G uniquely is of seismological interest. In the present paper we show that F determines G uniquely if we assume the source to be a simple straight line source (SSLS) or an ideal fault in an isotropic medium with isotropic prestress (IFIMIP). We give tests on F which determine whether it can come from a SSLS, from an IFIMIP or from a simple plane surface source (SPSS). If we assume the source to be a SPSS then knowing F and the fault plane determines G to within an unknown scalar multiple of a certain tensor tangent to the fault plane. Moreover F determines the fault plane uniquely unless F can come from a SSLS. If it can, then F determines this virtual source line uniquely, and F permits the fault plane to be any plane containing the virtual source line.     

Use of 40Ar/39Ar K-feldspar thermochronology in basin thermal history reconstruction: an example from the Big Lake Suite granites, Warburton Basin, South Australia

The potential use of ⁴⁰Ar/³⁹Ar thermochronologic data from K-feldspars in reconstructing basin thermal history has been evaluated using the example of the Warburton/Cooper/Eromanga Basin, Australia's largest onshore oil- and gas-producing basin. Results from ⁴⁰Ar/³⁹Ar step-heating experiments reveal details of the evolution of the basin system, including the following: (1) the operation of high geothermal gradient regimes during the earliest basin evolution, suggesting that basin formation was active rather than passive; (2) slow cooling from a Permo-Triassic temperature peak of at least 250–300°C; (3) a rise in thermal gradients to contemporary bottom hole temperatures in the last 5–10 Myr; and (4) spatially variable recrystallization events between 100 and 50 Ma and at around 20 Ma. Initial microstructural observations serve as a useful predictor of the quality and nature of the obtainable age information. Data from 'pristine' K-feldspars may constrain the peak temperature conditions experienced in the basin, the basin's early thermal history and also any recent changes in thermal gradient. Contrasting data from texturally modified K-feldspars may constrain times of thermal transients and/or fluid flow, with the preferred interpretation that K-feldspars recrystallize in response to such events. The Warburton/Cooper/Eromanga Basin example suggests that the ⁴⁰Ar/³⁹Ar technique may serve as a useful adjunct to apatite and zircon fission track analysis and conventional organic maturation indices in basin thermal history analysis.     

40Ar/39Ar dates of eastern Iceland lavas

Summary. Twenty-two plateau ages obtained by the ⁴⁰Ar/³⁹Ar step-heating method span most of the time of formation of eastern Iceland. With one exception the dates are consistent with the stratigraphy, the exception having a very unusual argon release pattern. It is concluded that the ⁴⁰Ar/³⁹Ar step-heating method can give correct dates on hydrothermally altered samples the total fusion ages of which are not correct; however, it is not possible to conclude from the results obtained that the method is inherently more accurate than the conventional method using carefully selected samples.
Though it is known from geological evidence that minor irregularities of spreading have occurred the results show that, within the resolution of the dates, spreading has been steady and symmetrical; in particular, neither a large eastward jump nor repetitive westward jumps can have occurred. The average half spreading rate for the period 2Myr to 12 Myr ago was 0.78 ± 0.16 cm yr-¹.     

Deep-water renewal in a Scottish fjord: temperature, salinity and oxygen isotopes

The sea lochs (fjords) of north-west Scotland are located in a region of Europe particularly well situated to monitor changes in westerly air streams. The moisture transported in these air streams has a profound effect on regional precipitation, freshwater run-off and, in turn, sea loch circulation. The gentle slope of the regional salinity:δ¹⁸O mixing-line, defined as 0.18 ‰ per salinity unit, suggests that the temperature: δ¹⁸O relationship may be readily resolved in these coastal waters. Deep-water renewal events, both observed and predicted from empirical models, in the bottom-waters of Loch Etive provide an opportunity to assess the temperature, salinity and δ¹⁸O relationship. Predicted changes in δ¹⁸O_calcite as a function of changing salinity (ΔS) and changing temperature (ΔT) during deep-water renewal events suggest that >80% fall above analytical detection limits. The theoretical likelihood of recording such renewal events in the "palaeoclimate" record appears to be promising, but temperature and salinity change during renewal events may have either sign. Scottish fjords, because of the relatively small impact which salinity has on δ¹⁸O_water, may provide useful study sites in palaeoclimate research, particularly where palaeotemperature is the primary record of interest.     

A linear approximation to the solution of a one-dimensional Stefan problem and its geophysical implications

Summary. The motion of a phase boundary in the Earth caused by temperature and pressure excitations at the Earth's surface is determined under a linear approximation. The solution is found as a sum of convolutions of pressure and temperature Green's functions with the corresponding excitations. The Green's functions are given under the form of Laplace transforms that can be inverted either by numerical evaluation of a branch cut integral or by inversion of a series expansion. This solution is a generalization of a solution previously derived by Gjevik. This latter solution is the first term in the series expansion. The relaxation times associated with the phase boundary motion are of the order of 10⁵–10⁷yr for the olivine—spinel phase transition and of 10⁶–10⁷yr for the basalt—eclogite transition. The linear approximation remains valid for long times only if the phase boundary moves slowly.     

On the penetration of a hot diapir through a strongly temperature-dependent viscosity medium

Summary. The ascent of a hot spherical body through a fluid with a strongly temperature-dependent viscosity has been studied using an axisymmetric finite element method. Numerical solutions range over Peclet numbers of 10⁻¹– 10³ from constant viscosity up to viscosity variations of 10⁵. Both rigid and stress-free boundary conditions were applied at the surface of the sphere. The dependence of drag on viscosity variation was shown to have no dependence on the stress boundary condition except for a Stokes flow scaling factor. A Nusselt number parameterization based on the stress-free constant viscosity functional dependence on the Peclet number scaled by a parameter depending on the viscosity structure fits both stress-free and rigid boundary condition data above viscosity variations of 100. The temperature scale height was determined as a function of sphere radius. For the simple physical model studied in this paper pre-heating is required to reduce the ambient viscosity of the country rock to less than 10²² cm² s⁻¹ in order for a 10 km diapir to penetrate a distance of several radii.     

Ice-free conditions in Novaya Zemlya 35 000–30 000 cal years B.P., as indicated by radiocarbon ages and amino acid racemization evidence from marine molluscs

Novaya Zemlya was covered by the eastern part of the Barents–Kara ice sheet during the glacial maximum of marine isotope stage 2 (MIS 2). We obtained ¹⁴C ages on 37 samples of mollusc shells from various sites on the islands. Most samples yielded ages in the range of 48–26 ¹⁴C Ky. Such old samples are sensitive to contamination by young ¹⁴C, and therefore their reliability was assessed using replicate analyses and amino acid geochronology. The extent of aspartic acid racemization (Asp D/L) indicates that many of the ¹⁴C ages are correct, whereas some are minimum ages only. The results indicate that a substantial part of Novaya Zemlya was ice-free about 35–27 ¹⁴C Kya, and probably even earlier. Corresponding shorelines up to >140 m a.s.l. indicate a large Barents–Kara ice sheet during early MIS 3. These results are consistent with findings from Svalbard and northern Russia: in both places a large MIS 4/3 Barents–Kara ice sheet is postulated to have retreated about 50 Kya, followed by an ice-free interstadial that lasted until up to ca. 25 Kya. The duration of the MIS 2 glaciation in Novaya Zemlya was calculated by applying the D/L values to a kinetic equation for Asp racemization. This indicates that the islands were ice covered for less than 3000 years if the basal temperature was 0^oC, and for less than 10 000 years if it was −5^oC.     

Effect of viscosity structure on fault potential and stress orientations in eastern Canada

Previous investigations of the causal relationship between postglacial rebound and earthquakes in eastern Canada have focused on the mode of failure and the observed timing of the pulse of earthquake/faulting activity following deglaciation. In this study, the observational database has been extended to include observed orientations of the contemporary stress field and the rotation of stress since deglacial times. It is shown that many of these observations can be explained by a realistic ice history and a viscoelastic earth with a uniform 10²¹ Pa s mantle.
The effects of viscosity structure on the above predictions are also examined. It is shown that, since most of the above observations are found within the ice margin, they are not very sensitive to lithospheric thickness. Also, the inclusion of a 25 or 50 km ductile layer within the lithosphere will not decouple the seismogenic upper crust. High viscosity (10²² Pa s) in the lower mantle is rejected by the stress orientation and rotation observations. A low-viscosity (6 times 10²⁰Pa s) upper mantle with 1.6 times 10²¹ Pa s in the upper part of the lower mantle and 3 times 10²¹ Pa s in the lower part of the lower mantle below 1200 km depth has been found to give predictions that are in general agreement with the observations.