Hostname: page-component-848d4c4894-v5vhk Total loading time: 0 Render date: 2024-06-21T18:48:37.714Z Has data issue: false hasContentIssue false

Theoretical Gravity Anomalies of Glaciers Having Parabolic Cross-Sections

Published online by Cambridge University Press:  30 January 2017

Charles E. Corbató*
Affiliation:
Department of Geology, University of California, Los Angeles, California, U.S.A.
Rights & Permissions [Opens in a new window]

Abstract

Equations and a graph are presented for calculating gravity anomalies on a two-dimensional glacier model having a horizontal upper boundary and a lower boundary which is a parabola with a vertical axis of symmetry.

Résumé

Résumé

On donne des équations et un graphique pour le calcul des anomalies de gravité sur la base d’un modèle bi-dimentionnel de glacier ayant une surface horizontale et une base de forme parabolique à axe de symétrie vertical.

Zusammenfassung

Zusammenfassung

Für die Berechnung der Schwereanomalien an einem zweidimensionalen Gletschermodell werden Gleichungen und eine graphische Darstellung angegeben. Das Modell ist oben durch eine Horizontale, unten durch eine Parabel mit vertikaler Symmetrieachse begrenzt.

Type
Research Article
Copyright
Copyright © International Glaciological Society 1964

Gravity measurements on valley glaciers have become more common during the last decade as investigators have attempted to determine subglacial bedrock configurations. The early work utilized very simple methods of analysis and only recently have more sophisticated procedures of two- and three-dimensional analysis been attempted. Reference KanasewichKanasewich (1963) found that a two-dimensional parabolic shape for the bedrock was a good first approximation to use for comparison with the observed gravity values; however, he used an approximation method to determine the gravity variations across this model.

The gravity effects of a two-dimensional glacier model can be calculated exactly for the case where the upper boundary is horizontal and the lower one a parabola with a vertical axis of symmetry. Let a be the half width and b the maximum depth of the glacier (Fig. 1 ). The vertical component of gravity, g z, at a point on the upper boundary of this body (Reference HeilandHeiland 1940, p. 146) can be expressed as

where γ is the gravitational constant, ρ is the density difference between the glacier and bedrock, ξ = x/a, ζ = z/a, β = b/a, and χ = x′/a, x′ being the distance of the gravity station from the center line of the glacier. Integrating with respect to ζ yields

Fig. 1. Notation used for the calculation of the gravity anomaly at point x′. The upper boundary of the glacier is horizontal and the lower boundary is a parabola with a vertical axis of symmetry

The final integration can easily be performed if the left-hand integrand is factored into

where

. By evaluating the integrals and re-arranging terms the following result is obtained when χ≠1

and when χ = 1,

The equations are shown graphically in Figure 2 for 0.1≤β≤2.0 and 0≤χ≤1.5 (The calculation of the 320 points used in determining these curves was done on an IBM 7090 computer in less than one minute.) Comparison of the curves with observed gravity anomalies allows one to make a rapid first approximation for the depth of a valley glacier.

Fig. 2. The junction gz/γρa plotted against x′/a for values of b/a from 0.1 to 2.0

Acknowledgement

The author wishes to thank the Computing Facility of the University of California, Los Angeles, for the generous use of its IBM 7090 high-speed digital computer.

References

Heiland, C. A. 1940. Geophysical exploration. New York, Prentice-Hall, Inc.Google Scholar
Kanasewich, E. R. 1963. Gravity measurements on the Athabaska Glacier, Alberta, Canada. Journal of Glaciology, Vol. 4, No. 35, p. 61731.CrossRefGoogle Scholar
Figure 0

Fig. 1. Notation used for the calculation of the gravity anomaly at point x′. The upper boundary of the glacier is horizontal and the lower boundary is a parabola with a vertical axis of symmetry

Figure 1

Fig. 2. The junction gz/γρa plotted against x′/a for values of b/a from 0.1 to 2.0