Preliminary palaeomagnetic results of an Archaean dolerite dyke of west Greenland: geomagnetic field intensity at 2.8 Ga

The geomagnetic field intensity during Archaean times is evaluated from a palaeomagnetic and chronological study of a dolerite dyke intruded into the 3000 Ma Nuuk Gneisses at Nuuk (64.2°N, 51.7°W), west Greenland. Plagioclase from the dolerite dyke yields a mean K-Ar age of 2752 Ma. Palaeomagnetic directions after thermal demagnetization of the dyke and the gneiss reveal a positive baked-contact test, indicating that the high-temperature-component magnetization of the dyke is primary. Thellier experiments on 12 dyke specimens yield a palaeointensity value of 13.5±4.4 μT. The virtual dipole moment at ca. 2.8 Ga is 1.9±0.6 × 10²² Am², which is about one-quarter of the present value. The present study and other available data imply that the Earth's magnetic field at 2.7 ∼ 2.8 Ga was characterized by a weak dipole moment and that a fairly strong geomagnetic field similar to the present intensity followed the weak field after ca. 2.6 Ga.     

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.     

A new method of palaeofield magnitude correction for thermally altered samples and its application to Lower Carboniferous lavas

Summary. Some lavas of Early Carboniferous age from Scotland and Derbyshire have been analysed for palaeofield magnitude by the application of a new correction to data obtained by the Shaw method. Its application has yielded field magnitude values from data which had previously been rejected. A modern pottery specimen and recent lavas from Sicily and Westman Island have also been analysed.
An average virtual dipole moment (VDM, Smith) of 2.5 ± 0.4 × 10²² A m² was determined for the Early Carboniferous. This is only 30 per cent of the present-day VDM.     

Centennial to millennial geomagnetic secular variation

A time-varying spherical harmonic model of the palaeomagnetic field for 0–7 ka is used to investigate large-scale global geomagnetic secular variation on centennial to millennial scales. We study dipole moment evolution over the past 7 kyr, and estimate its rate of change using the Gauss coefficients of degree 1 (dipole coefficients) from the CALS7K.2 field model and by two alternative methods that confirm the robustness of the predicted variations. All methods show substantial dipole moment variation on timescales ranging from centennial to millennial. The dipole moment from CALS7K.2 has the best resolution and is able to resolve the general decrease in dipole moment seen in historical observations since about 1830. The currently observed rate of dipole decay is underestimated by CALS7K.2, but is still not extraordinarily strong in comparison to the rates of change shown by the model over the whole 7 kyr interval. Truly continuous phases of dipole decrease or increase are decadal to centennial in length rather than longer-term features. The general large-scale secular variation shows substantial changes in power in higher spherical harmonic degrees on similar timescales to the dipole. Comparisons are made between statistical variations calculated directly from CALS7K.2 and longer-term palaeosecular variation models: CALS7K.2 has lower overall variance in the dipole and quadrupole terms, but exhibits an imbalance between dispersion in   g ¹₂  and   h ¹₂  , suggestive of long-term non-zonal structure in the secular variations.     

Reappraisal of surface wave magnitudes in the Eastern Mediterranean region and the Middle East

We report on a detailed palaeomagnetic study of the Miocene Farellones volcanic formation in the Chilean Andes near Santiago (two sections, 37 sites, about 400 orientated cores). Petrological observations show evidence of low-grade metamorphism increasing downwards through the volcanic sequence. Optical observations of opaque minerals and magnetic experiments suggest that in many cases maghemitization is associated with hydrothermal alteration. However, thermal demagnetization data indicate that the low-grade metamorphism did not significantly modify the direction of the primary remanent magnetization recorded at the time of emplacement of the volcanic lava flows. Four intervals of polarity with two intermediate palaeodirections were observed in the ~650-m-thick composite section. According to the dispersion of flow average directions, palaeosecular variation was slightly larger than that observed in general during the Upper Cenozoic. The site mean directions obtained in this study differ significantly from the expected Miocene direction. Clockwise rotations of up to 20° of small blocks are probably associated with the deformation of the Andean Cordillera since middle Miocene times. Geomagnetic palaeointensity data were obtained, using the Thellier method, on 24 samples from eight distinct lava flows. The flow mean VDM varies from 1.4 to 4.0 × 10²² A m⁻². Altogether, our data seem to suggest the existence of a relatively low geomagnetic field undergoing large fluctuations. Although a linear relationship was observed between the natural remanent magnetization and the thermal remanent magnetization acquired during the Thellier–Thellier experiments, undetected chemical alteration of the magnetic minerals during hydrothermalism may also explain the unusually low palaeointensity obtained.     

Observation of electric currents in the Alaska oil pipeline resulting from auroral electrojet current sources

Summary. Currents in the 1.28 × 10³ km Alaska oil pipeline, induced from the ionospheric, auroral electrojet, were measured at three sites, near Fairbanks, Paxson and Valdez, Alaska, using a gradient configuration of two fluxgate magnetometers. The observed pipeline current magnitudes, which reached 50 A during times of mild geomagnetic activity, displayed a linear relationship with the electric earth potential. Using the induction relationship between the electric and magnetic fields and the typical spectral composition of the geomagnetic field at high latitudes, I obtained a spectral appearance of the current that shows a maximum in the range of 4.5- to 10-min period. Near Fairbanks the pipeline current amplitudes, I (Amperes), could be represented, approximately, by I = 0.65 B _x T ^−0.5, where B _x(nT) is the north—south geomagnetic field variation amplitude and T (min) is its apparent period. There is much less pipelines current at the sites south of Fairbanks. A previously established relationship between the local electric field and the planetary geomagnetic activity index, Ap , permitted a prediction of the pipeline current surge amplitudes in the Fairbanks region as approximately I = 5.0 Ap . Current surges larger than 500 A may be expected rather often in the Alaska pipeline during large geomagnetic storms.     

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.     

The secular acceleration of the Earth's spin

Summary. The power spectrum of the Earth's spin has important components with periods ranging from a few days to at least a few thousand years, and probably to the age of the Earth. The secular acceleration, as the term is used here, refers to the components with periods longer than three centuries. In the year 600, the secular acceleration was —19.9 ± 0.8 parts in 10⁹ per century, while the value at the present time is less than half this size. The spin acceleration has important contributions from tidal friction and from an effect that is proportional to the square of the magnetic dipole moment. When these contributions are subtracted from the observed acceleration, we are left with a contribution that amounts to +41 parts in 10⁹ per century. This amount probably results from an unknown combination of changes in the size of the core, in the amount of glaciation, and in the size of the gravitational constant.     

Three-dimensional Rayleigh hysteresis of oriented core samples from the German Continental Deep Drilling Program: susceptibility tensor, Rayleigh tensor, three-dimensional Rayleigh law

Rayleigh hysteresis, as defined by the well-known Rayleigh relations, has been observed not only when magnetization of pyrrhotite-bearing KTB-samples is measured in parallel to a weak dc magnetic field, but also in experiments where field and measuring directions have been adjusted strictly perpendicularly to each other. Nine-tupels of independent Rayleigh hysteresis loops could thus be compiled. Their characteristic coefficients X ^ijk of initial susceptibility together with the Rayleigh loss coefficients α^jk have been proved to determine completely the samples' weak-field magnetic anisotropy. Interpreting the coefficient matrices ( X ^ijk) and (α^jk) as the tensor of initial susceptibility and the Rayleigh tensor, respectively, generalization of the isotropic Rayleigh relations in terms of corresponding tensor relationships has been suggested for the anisotropic case. Application to the KTB samples showed 3-D Rayleigh hysteresis measurements to be an excellent tool for rock magnetic analysis in terms of ore content and crystalline texture. In particular, a magnetocrystalline double texture of the basal planes of pyrrhotite precipitates and their [1120] directions of easy magnetization have been clearly detected. Surprisingly, the welt-known theorem α= const. X ²_I, formulated by Néel (1942) for the isotropic case, has been found to hold true even in tensor generalization (α_jk) = const ( X ²_jk). To reach sufficient resolution for the measurements performed, a sensitive vibrating coil magnetometer (VCM) has been developed.     

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 moment release during slab rupture beneath the Banda Sea

The highest intermediate depth moment release rates in Indonesia occur in the slab beneath the largely submerged segment of the Banda arc in the Banda Sea to the east of Roma, termed the Damar Zone. The most active, western-part of this zone is characterized by downdip extension, with moment release rates (∼10¹⁸ Nm yr^–1 per 50 km strike length) implying the slab is stretching at ∼10⁻¹⁴ s⁻¹ consistent with near complete slab decoupling across the 100–200 km depth range. Differential vertical stretching along the length of the Damar Zone is consistent with a slab rupture front at ∼100–200 km depth beneath Roma propagating eastwards at ∼100 km Myr^–1. Complexities in the slab deformation field are revealed by a narrow zone of anomalous in-plane P -axis trends beneath Damar, where subhorizontal constriction suggests extreme stress concentrations ∼100 km ahead of the slab rupture front. Such stress concentrations may explain the anomalously deep ocean gateways in this region, in which case ongoing slab rupture may have played a key role in modulating the Indonesian throughflow in the Banda Sea over the last few million years.     

The Influence of Surface Inhomogeneities On Deep Electromagnetic Soundings of the Earth

Summary. Asymptotic expressions for components of the electromagnetic field of a grounded electric dipole are considered for the model consisting of a thin surface-layer overlapping a stratified medium with a highly resistive screen on the roof. It is shown that the method of spatial derivatives makes it possible to obtain proper estimates of the impedance at distances of r ≥|λ₀| from the nearest edge of the surface anomaly (|λ₀| being the effective depth of the field penetration in the underlying section). the magnetotelluric methods allow one to obtain the true values of impedance, provided r ≥ max {|λ₀|, |/( S ⁻¹+ Z ⁻₀|^1/2} where S is the integrated conductivity of the surface layer, is the transverse resistance of the screen, and Z ⁻₀ is the Tikhonov—Cagniard impedance for the medium underlying the surface layer.     

Permo-Carboniferous (Kiaman) palaeointensity results from the western Bohemian Massif, Germany

Thellier palaeointensity analyses were carried out on samples from three quartz porphyry intrusions located within the western part of the Bohemian Massif close to the German-Czech Republic border. These intrusives have been dated at 280-275 Ma. which places them in the main part of the Permo-Carboniferous (Kiaman) Reversed Superchron (P-CRS). This age is consistent with the mean characteristic remanence direction calculated at D/ I = 204/- 1, giving a VGP at 36N, 341E. Palaeointensity values range from 6.6 μT to 22.1 μT, with a mean value of 12.0 ± 1.3 μT, which has a corresponding VDM value of (3.0 ± 1.2) × 10²² A m². This suggests that the geomagnetic field strength recorded by the quartz porphyries is only 37 per cent of the current value, supporting recent studies that report a field strength of between 25 and 40 per cent of the present-day value during the P-CRS.     

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.     

Applicability of time-to-failure analysis to accelerated strain before earthquakes and volcanic eruptions

We examine quantitatively the ranges of applicability of the equation Ω= A+B [1− t/t _f ] ^m for predicting 'system-sized' failure times t _f in the Earth. In applications Ω is a proxy measure for strain or crack length, and A , B and the index m are model parameters determined by curve fitting. We consider constitutive rules derived from (a) Charles' law for subcritical crack growth; (b) Voight's equation; and (c) a simple percolation model, and show in each case that this equation holds only when m < 0. When m > 0, the general solution takes the form Ω = A + B [1 + t / T  ] ^m , where T   is a positive time constant, and no failure time can be defined. Reported values for volcanic precursors based on rate data are found to be within the range of applicability of time-to-failure analysis ( m < 0). The same applies to seismic moment release before earthquakes, at the expense of poor retrospective predictability of the time of the a posteriori -defined main shock. In contrast, reported values based on increasing cumulative Benioff strain occur in the region where a system-sized failure time cannot be defined ( m > 0; commonly m ≈ 0.3). We conclude on physical grounds that cumulative seismic moment is preferred as the most direct measure of seismic strain. If cumulative Benioff strain is to be retained on empirical grounds, then it is important that these data either be re-examined with the independent constraint m < 0, or that for the case 0 < m + 1 < 1, a specific correction for the time-integration of cumulative data be applied, of the form ΣΩ = At + B '{1 − [1 − t/t _f ] ^m+1 }.     

Magnetostratigraphy of Palaeocene basalts from the Vaigat Formation of West Greenland

A palaeomagnetic study comprising the directional results from 289 individual lava flows, sampled along eight sections in the Palaeocene basalts of West Greenland, is reported. The eight individual sections are correlated using lithostratigraphical marker horizons to form a single composite profile. Generally, the lithological correlation is in good agreement with the record of geomagnetic secular variation.
  The total composite palaeomagnetic profile represents a stratigraphic thickness of 1.6  km through the Vaigat Formation, which is the lowermost of the two volcanic formations formed during the main stage of plateau volcanism. Only two polarity zones are found in the composite profile, suggesting a very short duration for the West Greenland main plateau-building volcanism. ⁴⁰Ar/³⁹Ar dates support a high extrusion rate and also indicate that the lower normal polarity zone is Chron C27n and that the upper reverse polarity zone is Chron C26r.
  The C27n–C26r transition is fully recorded along one of the sections (Nuusap Qaqqarsua), with intermediate directions covering a 200  m thick succession of lavas. A combined palaeomagnetic, field and geochemical study along this profile showed good agreement; that is, geochemically and geologically derived single magmatic events show groupings of the palaeomagnetic directions. Supposing a duration for the geomagnetic transition of 5000 years, the eruption frequency during this period was as high as one flow every 80 years.     

Palaeomagnetism of alkaline lava flows from El Salto-Almafuerte, Córdoba Province, Argentina

Summary. Stable components of magnetization have been isolated in 15 lava flows (mean K-Ar age 123 ± 4 Myr) from the alkaline sequence outcropping at El Salto-Almafuerte, Province of Cordoba, Argentina. Magnetic and geologic stratigraphy, as well as K-Ar ages indicate that this sequence was probably extruded in the Lower Cretaceous during the first volcanic cycle of the Sierra de los Cóndores Group (Vulcanitas Cerro Colorado Formation).
The palaeomagnetic pole-position for El Salto-Almafuerte lava flows, computed from the mean of 15 virtual geomagnetic poles and denoted SAK₇, is: 25° E, 72° S ( k = 35, α₉₅= 6.5°); it is fairly close to other Lower Cretaceous palaeomagnetic poles for South America. The elongated distribution of Cretaceous palaeomagnetic poles suggest recurrent drift for South America in early Cretaceous time.
The palaeomagnetic and radiometric data for the igneous rocks from El Salto-Almafuerte support the magnetic reversal time-scale for the early Cretaceous suggested by oceanic magnetic lineations.     

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.     

Depositional and tectonic evolution of a supradetachment basin: 40Ar/39Ar geochronology of the Nova Formation, Panamint Range, California

The Nova Basin contains an upper Miocene to Pliocene supradetachment sedimentary succession that records the unroofing of the Panamint metamorphic core complex, west of Death Valley, California. Basin stratigraphy reflects the evolution of sedimentation processes from landslide emplacement during basin initiation to the development of alluvial fans composed of reworked, uplifted sections of the basin fill. ⁴⁰Ar/³⁹Ar geochronology of volcanic units in middle and lower parts of the sequence provide age control on the tectonic and depositional evolution of the basin and, more generally, insights regarding the rate of change of depositional environments in supradetachment basins. Our work, along with earlier research, indicate basin deposition from 11.38 Ma to 3.35 Ma. The data imply sedimentation rates, uncorrected for compaction, of ~100 m Myr⁻¹ in the lower, high-energy part to ~1000 m Myr⁻¹ in the middle part characterized by debris-flow fan deposition. The observed variation in sediment flux rate during basin evolution suggests that supradetachment basins have complex depositional histories involving rapid transitions in both the style and rate of sedimentation.     

Uplift and heat flow following the injection of magmas into the lithosphere

Summary. The thermal effect of a rapid injection of hot magmas into the lower part of the lithosphere is modelled as an increase in heat production through the invaded region. The change in surface heat flow and the uplift resulting from the thermal expansion are determined in three-dimensional axially symmetric geometry: they are expressed as the space time convolutions of a Green's function with the anomalous heat production.
The anomalies with shorter wavelength (compared to the lithospheric thickness) are attenuated. This filtering affects the surface uplift more than the heat flow anomaly; the attenuation effect is larger when only the lower part of the lithosphere is invaded.
The uplift time constant is of the same order as the heat conduction time if the lower lithosphere is invaded by magmas at a moderate rate (i.e. the rate of injection does not exceed the equivalent of 0.1 per cent of the lithospheric volume in 10⁶yr). Fifty per cent of the total uplift takes place in about 80 × 10⁶yr for a lithosphere 100 km thick. The uplift is slightly faster when the whole lithosphere is invaded. The heat flow anomaly is delayed when the lower part of the lithosphere is invaded.
The spatial extent and the timing of the uplift and heat flow anomalies are critical in determining the mechanism's feasibility. Magma injections explain rapid uplifts [> 100 m (10⁶ yr)⁻¹] only if the magma is supplied at a very high rate (i.e. at least 10 per cent of the lithosphere volume per 10⁶yr). It is a feasible mechanism for uplifts that occur over longer periods of time (≊ 30 × 10⁶yr) such as those that seem to have occurred when the African plate came to rest with respect to the mantle.