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
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary-alt:20170925211724-11483-mediumThumb-S0022143000028860_eqn1.jpg?pub-status=live)
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
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary-alt:20170925211724-46797-mediumThumb-S0022143000028860_eqn2.jpg?pub-status=live)
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary-alt:20170925211724-17763-mediumThumb-S0022143000028860_fig3g.jpg?pub-status=live)
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
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary-alt:20170925211724-60800-mediumThumb-S0022143000028860_eqn4.jpg?pub-status=live)
where
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20170925050843300-0248:S0022143000028860:S0022143000028860_eqn5.gif?pub-status=live)
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary-alt:20170925211724-94351-mediumThumb-S0022143000028860_eqn6.jpg?pub-status=live)
and when χ = 1,
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary-alt:20170925211724-05815-mediumThumb-S0022143000028860_eqn7.jpg?pub-status=live)
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.
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary-alt:20170925211724-52851-mediumThumb-S0022143000028860_fig8g.jpg?pub-status=live)
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.