doi:10.1016/j.epsl.2005.04.016
Copyright © 2005 Elsevier B.V. All rights reserved.
Gravity anomalies of sedimentary basins and their mechanical implications: Application to the Ross Sea basins, West Antarctica
References and further reading may be available for this article. To view references and further reading you must
purchase this article.
Garry D. Karner
,
, Michael Studinger
and Robin E. Bell
Lamont–Doherty Earth Observatory of Columbia University, 61 Route 9W, Palisades, NY 10964, USA
Received 19 October 2004;
revised 8 April 2005;
accepted 11 April 2005.
Editor: V. Courtillot.
Available online 6 June 2005.
Abstract
In general, sedimentary basins are characterized by negative free-air and Bouguer gravity anomalies. However, the extensional basins of the Ross Sea are paradoxical in that positive gravity anomalies overlay the Victoria Land Basin, Northern Basin, Central Trough and Northern Central Trough while basement highs are associated with negative gravity anomalies. Measured basement densities from DSDP basement cores give values between 2600–2800 kg/m3 while bulk sediment densities range from 1210–2200 kg/m3, indicating a normal density relationship between basement and sediment infill. In contrast, the relatively young and narrow Terror Rift is associated with negative free-air and Bouguer gravity anomalies, but has a different geological history as compared to the larger Ross Sea basins. Process-oriented gravity modeling indicates that magmatic underplating and crustal intrusions are inconsistent with the observed gravity and basement geometry of the Ross Sea basins. The magma volume necessary to account for the distribution and amplitude of the positive gravity anomaly of the Central Basin and be isostatically balanced would need to be comparable to the tholeiitic flood basalt volume of the Columbia River province—it is thus unlikely that the volume of Neogene volcanics of the Ross Sea region is sufficient to explain the observed gravity relationship by modifying the bulk density of the crust.
We demonstrate that positive free-air and Bouguer gravity anomalies over extensional basins are the consequence of a relatively low flexural strength of the lithosphere during rifting being contrasted by higher flexural strengths later during sedimentation. As the difference between the rigidity of the lithosphere during sedimentation increases relative to the rigidity of the rifted lithosphere, the gravity over the basin becomes progressively more positive but only for a limited range of wavelengths. The narrow width of the Terror Rift precludes it from having a positive gravity anomaly while the opposite is true for the large Ross Sea basins. For the Ross Sea region, such a loading scenario requires a significant delay between extension and the timing of sediment infilling of the basins, consistent with the late Cretaceous extension of the Ross Sea region and the sedimentary succession being dominated by large-scale late Eocene–Neogene glaciogenic progradational sequences. Sediment source was presumably from the denudation of the Transantarctic Mountains, which commenced in the late Paleogene. The time delay between the late Cretaceous formation of the Transantarctic Mountains, late Paleogene exhumation, and the generation of significant Paleogene paleobathymetry requires either the Ross Sea region to be sub-aerial and sediment starved for most of the Paleogene and/or the Paleogene climate was ineffective in producing clastics until the onset of glaciation in the late Eocene–early Oligocene.
Keywords: gravity anomalies; isostasy; lithospheric flexure; extensional basins; Paleogene climate; Antarctica
Fig. 1. General bathymetry of the Ross Sea region [60]. Contour interval is 50 m between 0–1000 m and 500 m between 1000–5000 m. Geological features discussed in the text are located, along with the position of the Cape Roberts and DSDP drilling sites. The thin, white lines locate the main sedimentary basins and basement highs comprising the Ross Sea. Apart from the bathymetric escarpment separating the transition from the Ross Sea to the Pacific Ocean (and the approximate location of the ocean–continent boundary), the most conspicuous features of the bathymetry are northeast–southwest glacial scars.
Fig. 2. Bouguer gravity anomaly of the Ross sea region [60]. Reduction density is 2400 kg/m3. Contour interval is 10 mGal. The location of the main sedimentary basins (thin white lines) and interbasinal highs are superposed and clearly show a general anti-correlation between Bouguer gravity high and the location of sedimentary basins. Basement highs are associated with gravity lows. Location of the late Eocene–Neogene Terror Rift is shown by the bold yellow line while the surface outcrop of alkali basalts is shown by the thin red line. Seismic profiles shown in Fig. 3 and Fig. 4 are also located.
Fig. 3. Seismic reflection profile and corresponding Bouguer gravity anomaly profile across the Victoria Land Basin, Central Trough, and Eastern Basin of the Ross Sea [24]. Clearly shown is the anti-correlation between basin depocenters and positive Bouguer gravity anomalies.
Fig. 4. Depth to magnetic basement (gray circles) superimposed on basement (bold black line) mapped from seismic reflection data across the western Ross Sea. The Bouguer gravity [60] and magnetic data from the GANOVEX IV mission [42] are shown for reference. The gray circles are the Werner deconvolution solutions for depth to magnetic basement with the size of the circles being scaled to magnetic susceptibility.
 |
Fig. 5. Kinematic and isostatic modeling procedure for extensional basins and passive margins. (a) Crustal stretching factor, β(x), for single and multiple faulting. After faulting, the resultant crustal thickness is tc(1 − 1 / β(x)). The form of β(x) mirrors the distribution and amplitude of accommodation generated across the region. The basin-bounding fault of the system is at 800 km and after 1200 km, the crustal extension is a constant β(x) of 2. (b) Kinematic description of simple slip along a single border fault that produces a basin by the collapse of the hanging wall block. The geometry of the fault governs the size and shape of the resulting basin. (c) Kinematic description of multiple faulting of the crust. (d) Local isostatic adjustment of the crustal extension shown in (c) assuming water as the infill material and generation of the compensating Moho topography. (e) Local isostatic adjustment of the crustal extension shown in (c) assuming sediment as the infill material and generation of the compensating Moho topography.
 |
Fig. 6. Modeled gravity across a series of rift depocenters simulating the general geometry and distribution of basins comprising the Ross Sea region (scheme shown in Fig. 5d). Water, sediment, crustal and mantle densities are assumed to be 1030, 2400, 2800 and 3300 kg/m3, respectively. (a) Basin geometry and modeled gravity assuming local isostasy and no sedimentation. Major depocenters are characterized by a negative gravity effect. (b) Basin geometry and modeled gravity assuming local isostasy and sediment infilling. The resulting deflection locally depresses and subdues the Moho topography. Again, major depocenters are characterized by a negative gravity effect. (c) Basin geometry and modeled gravity assuming that the infilling sediments load a crust with finite flexural strength (Te of 30 km). In this case, the deflection regionally depresses the Moho topography. In particular, the Moho relief is preserved rather than attenuated as in the low-strength case. Now, the major depocenters are associated with a gravity high while basement highs are characterized by gravity lows.
 |
Fig. 7. Modeled free-air and Bouguer gravity across a series of rift depocenters simulating the general geometry and distribution of basins comprising the Ross Sea region (scheme shown in Fig. 5d). Water, sediment, crustal and mantle densities are assumed to be 1030, 2400, 2800 and 3300 kg/m3, respectively. (a) Basin geometry and modeled free-air (bold black line) and Bouguer (thin red line) gravity assuming a rift effective elastic thickness, Te, of 5 km and no sedimentation. Major depocenters are characterized by a negative free-air gravity. In contrast, the Bouguer gravity is dominated by the local and regional variations of the Moho. (b) Basin geometry and modeled free-air and Bouguer gravity assuming an effective elastic thickness of 5 km during both rifting and sedimentation. Major depocenters are characterized by negative free-air and Bouguer gravity. (c) Basin geometry and modeled free-air and Bouguer gravity assuming that the infilling sediments load a crust with a Te of 30 km while the rifting event is associated with a Te of 5 km. The major depocenters are characterized by gravity highs (both free-air and Bouguer) while the basement highs are characterized by gravity lows.
 |
Fig. 8. Process-oriented gravity modeling investigating the contribution of magmatic underplating and crustal intrusions to the gravity anomaly. Local isostasy is assumed during crustal thinning. During sedimentation, magmatic underplating and crustal intrusions the effective elastic thickness is assumed to be either zero (local isostasy; left-hand column) or a Te of 30 km (right-hand column). (a) Gravity effect across a modeled sedimentary basin assuming local isostasy during extension and either local isostasy (left-hand column) or a Te of 30 km (right-hand column) during sedimentation. For locally compensated systems, the gravity negative correlates with the localization of the basin depocenter. In contrast, for sedimentation with Te of 30 km, the gravity across the basin depocenter is positive. (b) Gravity effect across a modeled sedimentary basin assuming local isostasy during extension and either local isostasy (left-hand column) or a Te of 30 km (right-hand column) during magmatic underplating. Underplating induces uplift and for local isostasy, inverts the basin center while adding to the crustal thickness. The resultant gravity is positive but correlates with a zone of major basement uplift within the basin. In contrast, for magmatic underplating when Te is relatively large, the basin is regionally uplifted and the gravity over the basin depocenter is negative. (c) Gravity effect across a modeled sedimentary basin assuming local isostasy during extension and either local isostasy (left-hand column) or a Te of 30 km (right-hand column). Crustal intrusions effectively increase the density of the lower crust and induce subsidence. Two intrusion densities were assumed, 2900 and 3000 kg/m3 (red and black lines, respectively). For local isostasy, basin subsidence is amplified and the gravity effect is negative. In contrast, if Te is significantly large during the intrusion process, the basin is forced to regionally subside. The gravity effect is positive over the basin depocenter.
 |
Fig. 9. Admittance function Z(k) between the Moho and basin gravity as a function of varying rift and sedimentation flexural rigidities. Two families of admittance functions are shown depending on the effective elastic thickness of the extended lithosphere; Ter of 5 km (black curves), Ter of 30 km (red curves). For Ter of 5 km, sediment loading effective elastic thickness, Tes, ranges from 5–30 km while for a Ter of 30 km, sediment loading effective elastic thickness ranges from 30–80 km. For Z(k) < 1, the gravity effect over the basin is dominated by the density variation of the basin. For Z(k) > 1, the gravity is dominated by the Moho and the basin is associated with a positive gravity effect. When Ter and Tes are equal, the gravity is always negative across the basin (Z(k) < 1). As the difference between rift and sediment loading effective elastic thickness increases, the gravity over the basin becomes progressively more positive but only for a limited range of wavelengths. The wide late Cretaceous basins of the Ross Sea should be characterized by a positive gravity anomaly while the narrow late Paleogene–Neogene Terror Rift should be associated with a negative gravity anomaly.
Table 1.
Ross Sea basins modeling parameters


Corresponding author. Tel.: +1 845 365 8355.