Palaeomagnetic distortion modelling and possible recovery by inversion

A robust finite-element technique is presented for computation of both the internal demagnetization effects and magnetic terrain effects in bodies with arbitrary shape and arbitrary susceptibility distribution. This method facilitates a flexible analysis of the palaeomagnetic deflection problem. Tests on geologically realistic settings of highly magnetic rocks demonstrate that deflections of several degrees may occur even for relatively simple two-dimensional models. Similarly, the magnetic intensity may well be biased by 5-15% by demagnetization effects. The present paper focuses on deflections and intensity variations inside the magnetized body, where we find a systematic shallowing of inclination for bodies with a horizontal elongation. Because the bodies sampled at a typical palaeomagnetic site will have a dominant direction of elongation, the magnetic deflection effect will tend to impose a systematic bias which doesn’t average out. An inversion-based procedure for elimination of the deflection effect is presented. It requires that the magnetic body is quite homogeneous and that its surface geometry is known, as may be the case for historical lava flows. Tests demonstrate that in order to recover both ambient palaeofield direction and the effective susceptibility at blocking temperature it is necessary to sample near strong topographic elements in the magnetic body. Since the surface geometry rarely is known it is proposed as an alternative to inversion that an effective susceptibility is assessed and a horizontal slab correction is applied for samples taken far from topographical features. When shape geometry is unknown and no correction applied, palaeomagnetic conclusions must take into account the possible bias from internal demagnetization and magnetic terrain effects.     

New paleomagnetic pole and magnetostratigraphy of Faroe Islands flood volcanics, North Atlantic igneous province

A paleomagnetic sampling was carried out along four sections (altogether 86 lava flows, 548 samples) in the North Atlantic Igneous Province outcropping in Faroe Islands, Denmark. The four polarity zones in the 700-m-thick exposed part of the Faroes lower formation can be correlated with the geomagnetic polarity time scale as C26n-C25r-C25n-C24r. The seven lava flows erupted during C25n indicate a very low eruption rate in the upper part of the Faroes lower formation of ∼1/70 kyr. The Faroes middle and upper formations (composite thickness ∼2300 m) are all reversely magnetized corresponding to C24r. The eruption rate at the onset of middle formation volcanism was very high as evidenced by several thick lava sequences recording essentially spot readings of the paleomagnetic field. The shift in eruption rate between the Faroes lower and middle formations and evidence that onset of the Faroes middle formation volcanism took place in C24r are of particular importance, placing onset of middle formation volcanism in close temporal relation to North Atlantic continental break-up and the late Paleocene thermal maximum. After grouping flows recording the same field directions, we obtained 43 independent readings of the paleomagnetic field, yielding a paleomagnetic pole with coordinates 71.4°N, 154.7°E (A₉₅=6.0°, K=14, N=43); age 55-58 Ma. The pole is supported by a positive reversal test. Paleosecular variation, estimated as the angular standard deviation of the virtual geomagnetic pole distribution 21.7°+3.9°/−2.8°, is close to expected for the given age and paleolatitude. Our new Faroes paleomagnetic pole is statistically different from the majority of previously published poles from the British and Faroes igneous provinces, and we suggest that these older data should be used with care.     

Electrical properties of Titan's surface from Cassini RADAR scatterometer measurements

We report regional-scale low-resolution backscatter images of Titan's surface acquired by the Cassini RADAR scatterometer at a wavelength of 2.18-cm. We find that the average angular dependence of the backscatter from large regions and from specific surface features is consistent with a model composed of a quasi-specular Hagfors term plus a diffuse cosine component. A Gaussian quasi-specular term also fits the data, but less well than the Hagfors term. We derive values for the mean dielectric constant and root-mean-square (rms) slope of the surface from the quasi-specular term, which we ascribe to scattering from the surface interface only. The diffuse term accommodates contributions from volume scattering, multiple scattering, or wavelength-scale near-surface structure. The Hagfors model results imply a surface with regional mean dielectric constants between 1.9 and 3.6 and regional surface roughness that varies between 5.3° and 13.4° in rms-slope. Dielectric constants between 2 and 3 are expected for a surface composed of solid simple hydrocarbons, water ice, or a mixture of both. Smaller dielectric constants, between 1.6 and 1.9, are consistent with liquid hydrocarbons, while larger dielectric constants, near 4.5, may indicate the presence of water-ammonia ice [Lorenz, R.D., 1998. Icarus 136, 344-348] or organic heteropolymers [Thompson, W.R., Squyres, S.W., 1990. Icarus 86, 336-354]. We present backscatter images corrected for angular effects using the model residuals, which show strong features that correspond roughly to those in 0.94-μm ISS images. We model the localized backscatter from specific features to estimate dielectric constant and rms slope when the angular coverage is within the quasi-specular part of the backscatter curve. Only two apparent surface features are scanned with angular coverage sufficient for accurate modeling. Data from the bright albedo feature Quivira suggests a dielectric constant near 2.8 and rms slope near 10.1°. The dark albedo feature Shangri-La is best fit by a Hagfors model with a dielectric constant close to 2.4 and an rms slope near 9.5°. From the modeled backscatter curves, we find the average radar albedo in the same linear (SL) polarization to be near 0.34. We constrain the total-power albedo in order to compare the measurements with available groundbased radar results, which are typically obtained in both senses of circular polarization. We estimate an upper limit of 0.4 on the total-power albedo, a value that is significantly higher than the 0.21 total albedo value measured at 13 cm [Campbell, D., Black, G., Carter, L., Ostro, S., 2003. Science 302, 431-434]. This is consistent with a surface that has more small-scale structure and is thus more reflective at 2-cm than 13-cm. We compare results across overlapping observations and observe that the reduction and analysis are repeatable and consistent. We also confirm the strong correlations between radar and near-infrared images.     

Magnetically moderated outbursts of WZ Sagittae

We argue that the quiescent value of the viscosity parameter of the accretion disc in WZ Sge may be  α_cold∼ 0.01  , in agreement with estimates of α_cold for other dwarf novae. Assuming the white dwarf in WZ Sge to be magnetic, we show that, in quiescence, material close to the white dwarf can be propelled to larger radii, depleting the inner accretion disc. The propeller therefore has the effect of stabilizing the inner disc and allowing the outer disc to accumulate mass. The outbursts of WZ Sge are then regulated by the (magnetically determined) evolution of the surface density of the outer disc at a radius close to the tidal limit. Numerical models confirm that the recurrence time can be significantly extended in this way. The outbursts are expected to be superoutbursts since the outer disc radius is forced to exceed the tidal (3:1 resonance) radius. The large, quiescent disc is expected to be massive, and to be able to supply the observed mass accretion rate during outburst. We predict that the long-term spin evolution of the white dwarf spin will involve a long cycle of spin-up and spin-down phases.     

Dissolved organic matter (DOM) export to a temperate estuary: seasonal variations and implications of land use

Inputs of dissolved carbon, nitrogen, and phosphorus were assessed for an estuary and its catchment (Horsens, Denmark). Seasonal patterns in the concentrations of DOM in the freshwater supply to the estuary differed depending on the soil and drainage characteristics of the area. In streams draining more natural areas the, patterns observed were largely driven by seasonal temperature fluctuations. The material exported from agricultural areas was more variable and largely controlled by precipitation events. Positive exponential relationships were found between the nitrogen and phosphorus loading, and the percentage of catchment area used for agriculture. Colored DOM (CDOM) loading measurements were found to be a good predictor of dissolved organic carbon (DOC) loading across the different subcatchments, offering a rapid and inexpensive alternative of operationally monitoring DOC export. For all the dissolved nutrient inputs to the estuary, dissolved inorganic nitrogen (DIN) and dissolved organic phosphorus dominated the loadings. Although 81% of the nitrogen annually supplied to the estuary was DIN, 83% of the nitrogen exported from the estuary was dissolved organic nitrogen (DON). Results show that increasing the area of the catchment covered by forest and natural pastures would have a positive effect on the trophic status of the estuary, leading to a considerable decrease in the phosphorus loading and a shift in the nitrogen loading from DIN to DON. Such a change in land use would also increase the export of DOC and CDOM to the estuary having the potential to increase oxygen consumption and reduce the photic depth.     

Geomagnetic variations recorded in Older (≫ 23 000 B.P.) and Younger Yoldia Clay (∼ 14 000 B.P.) at Nørre Lyngby, Denmark

Summary. Detailed palaeomagnetic results from a rapidly deposited 12.5-m Younger Yoldia Clay sequence of age around 14 000 B.P. at the coast cliff at Nørre Lyngby (northern Jutland, Denmark) and a 2-m Older Yoldia Clay sequence of age somewhere between 23 000 and 40 000 B.P. are presented and discussed. The Younger Yoldia Clay spans some 1000–1500 yr and shows swings in inclination and declination of about that period, and also more rapid oscillations which are particularly marked in inclination, showing that rapid secular variations as have occurred during historic times were indeed also present back in time. There exist easterly declinations of 80° to 90° in the upper half of the Younger Yoldia Clay which cause the virtual geomagnetic pole to migrate clockwise to around 50° away from the rotation pole. This we name the 'Nørre Lyngby declination excursion'.
In the Older Yoldia Clay, as well as secular variations in both declination and inclination, significant low inclination values are found, confirming the existence of the recently named 'Rubjerg low inclination excursion', with the virtual geomagnetic pole moving first in a clockwise then in an anticlockwise sense at 40° to 60° away from the rotation pole.
It is therefore inferred that models for the Earth's geomagnetic field should involve at least local rapid eastward as well as westward 'drift' of the non-dipole field components at various times in the past.     

Preface and Dedication

Earth, Moon, and Planets -