Publications

Export 4 results:
Sort by: Author Title Type [ Year  (Desc)]
2019
Avery, MS, Constable CG, Davies CJ, Gubbins D.  2019.  Spectral methods for analyzing energy balances in geodynamo simulations. Physics of the Earth and Planetary Interiors. 286:127-137.   10.1016/j.pepi.2018.10.002   AbstractWebsite

The geomagnetic field displays complicated variations over a broad range of frequencies. These variations can be decomposed by frequency and linked to physical processes using frequency domain spectral methods. These spectral methods are well developed but have not previously been applied to study the energy balance of geodynamo simulations. We illustrate their potential by analyzing output from numerical dynamo simulations that have previously been studied for their apparently Earth-like properties. We show that high coherence between variations in axial dipole energy at the outer boundary of the simulation and total magnetic energy within the fluid shell occur at frequencies below similar to 0.1 kyr(-1). This suggests that paleomagnetically-observable signals with periods exceeding 10 kyrs contain information about magnetic energy changes in the bulk core. We then use spectral analysis to investigate differences in the rate of growth and decay of the axial dipole field. This behaviour, characterised by rapid growth and slow decay, is observed when signals with frequencies higher than 0.03 kyr(-1) have been filtered out. The origin of this asymmetric growth and decay is assessed using coherence spectra between rates of change in kinetic and magnetic energy, ohmic and viscous dissipation, and work done by the buoyancy and Lorentz forces. We show that asymmetry is associated with an imbalance between ohmic dissipation and work done by the Lorentz force; when changes in magnetic energy are more coherent with ohmic dissipation the field grows rapidly and decay slowly. Variations in Ohmic dissipation reflect changes in field strength in our models, while changes in viscous dissipation are associated with amplitude fluctuations of the large-scale flow that exists on millennial timescales. Our work shows that spectral analysis coupling observable and global products of the dynamo process can elucidate the physical origin of periodic processes occurring on timescales exceeding 10 kyrs.

2006
McMillan, DG, Constable CG.  2006.  Limitations in correlation of regional relative geomagnetic paleointensity. Geochemistry Geophysics Geosystems. 7   10.1029/2006gc001350   AbstractWebsite

Time domain correlations of common features among relative paleointensity records from sedimentary cores are invaluable to paleomagnetism and paleoclimatology. Sediments with high accumulation rates might now provide millennial scale correlations of temporal variations in the geomagnetic dipole moment. Errors in the ages of paleomagnetic data samples, however, can make such correlations difficult and unreliable. We use spectral methods to assess the level of coherence expected among individual and stacked high- resolution simulated paleointensity records for the time interval 0 - 75 ka. Correlations between individual paleointensity records are systematically degraded with decreased sedimentation rate and increased magnitude of age errors. We find that with optimistic age errors and interpolation of depth sampled data to evenly spaced time series, only short period signal in high- resolution relative paleointensity is corrupted. For currently available methods of establishing chronologies, we estimate the minimum characteristic timescale of correlative features between pairs of regional stacked records at about 4.5 kyr. From an analysis of NAPIS- 75 and SAPIS data, it appears that the limit is inherent to the regional stacks and not a consequence of comparison of distant, independent data sets. A detailed comparison of the NAPIS- 75 and SAPIS stacks shows that this limit is likely larger, perhaps 6 kyr. At long periods the two regional stacks are more poorly correlated than those from our simulations, suggesting somewhat larger age errors in the individual paleointensity records.

2004
McMillan, DG, Constable CG, Parker RL.  2004.  Assessing the dipolar signal in stacked paleointensity records using a statistical error model and geodynamo simulations. Physics of the Earth and Planetary Interiors. 145:37-54.   10.1016/j.pepi.2004.02.011   AbstractWebsite

Stacks of globally distributed relative paleointensity records from sediment cores are used to study temporal variations in the strength of the geomagnetic dipole. We assess the intrinsic accuracy and resolution of such stacks, which may be limited by errors in paleointensity, non-dipole field contributions, and the age scales assigned to each sediment core. Our approach employs two types of simulations. Numerical geodynamo models generate accurate predictions of time series of magnetic variations anywhere in the world. The predicted variations are then degraded using an appropriate statistical model to simulate expected age and paleointensity errors. A series of experiments identify the major contributors to error and loss of resolution in the resulting stacks. The statistical model simulates rock magnetic and measurement errors in paleointensity, and age errors due to finite sampling and approximations inherent in interpolation, incomplete or inaccurate tie point information, and sedimentation rate variations. Data sampling and interpolation to a designated age scale cause substantial decorrelation, and control the maximum level of agreement attainable between completely accurate records. The particular method of interpolation appears to have little effect on the coherence between accurate records, but denser tie point data improve the agreement. Age errors decorrelate geomagnetic signals, usually at shorter periods, although they can destroy coherence over a broad range of periods. The poor correlation between neighboring paleomagnetic records often observed in real data can be accounted for by age errors of moderate magnitude. In a global dataset of 20 records, modeled after the SINT800 compilation and spanning 300 kyr, our results show that dipole variations with periods longer than about 20 kyr can be recovered by the stacking process. Reasonable contributions to error in the paleointensity itself have a modest influence on the result, as do non-dipole field contributions whose effect is minor at periods longer than 10 kyr. Modest errors in the ages of tie points probably account for most of the degradation in geomagnetic signal. Stacked sedimentary paleomagnetic records can be improved by denser temporal sampling and careful selection of independent high-quality tie points. (C) 2004 Elsevier B.V. All rights reserved.

2002
McMillan, DG, Constable CG, Parker RL.  2002.  Limitations on stratigraphic analyses due to incomplete age control and their relevance to sedimentary paleomagnetism. Earth and Planetary Science Letters. 201:509-523.   10.1016/s0012-821x(02)00747-1   AbstractWebsite

A major limitation in the analysis of physical quantities measured from a stratigraphic core is incomplete knowledge of the depth to age relationship for the core. Records derived from diverse locations are often compared or combined to construct records that represent a global signal. Time series analysis of individual or combined records is commonly employed to seek quasi-periodic components or characterize the timescales of relevant physical processes. Assumptions that are frequently made in the approximation of depth to age relationships can have a dramatic and harmful effect on the spectral content of records from stratigraphic cores. A common procedure for estimating ages in a set of samples from a stratigraphic core is to assign, based on complementary data, the ages at a number of depths (tie points) and then assume a uniform accumulation rate between the tie points. Imprecisely dated or misidentified tie points and naturally varying accumulation rates give rise to discrepancies between the inferred and the actual ages of a sample. We develop a statistical model for age uncertainties in stratigraphic cores that treats the true, but in practice unknown, ages of core samples as random variables. For inaccuracies in the ages of tie points, we draw the error from a zero-mean normal distribution. For a variable accumulation rate, we require the actual ages of a sequence of samples to be monotonically increasing and the age errors to have the form of a Brownian bridge. That is, the errors are zero at the tie points. The actual ages are modeled by integrating a piecewise constant, randomly varying accumulation rate. In each case, our analysis yields closed form expressions for the expected value and variance of resulting errors in age at any depth in the core. By Monte Carlo simulation with plausible parameters, we find that age errors across a paleomagnetic record due to misdated tie points are likely of the same order as the tie point discrepancies. Those due to accumulation rate variations can be as large as 30 kyr, but are probably less than 10 kyr. We provide a method by which error estimates like these can be made for similar stratigraphic dating problems and apply our statistical model to an idealized marine sedimentary paleomagnetic record. Both types of errors severely degrade the spectral content of the inferred record. We quantify these effects using realistic tie point ages, their uncertainties and depositional parameters. (C) 2002 Elsevier Science B.V. All rights reserved.