M. J. Hoggard

Spatial and temporal patterns of Cenozoic dynamic topography around Australia

[1]xa0Despite its importance, the spatial and temporal pattern of dynamic topography generated by mantle convective circulation is poorly known. We present accurate estimates of dynamic topography from oceanic basins and continental margins surrounding Australia. Our starting point is measurement of residual depth anomalies on the oldest oceanic floor adjacent to the continental shelf. These anomalies were determined from a combined dataset of ~200 seismic reflection and wide-angle images of well-sedimented oceanic crust. They have amplitudes of between −1u2009km and +0.5u2009km, and their spatial variation is broadly consistent with long-wavelength free-air gravity and shallow seismic tomographic anomalies. Along the Northwest Shelf, a regional depth anomaly of −300 to −700u2009m intersects the adjacent continental shelf. The temporal evolution of this anomaly was determined by analyzing the stratigraphic architecture of an extensive carbonate platform, which fringes the shelf and records a dramatic switch from progradation to aggradation during Neogene times. Three-dimensional seismic mapping calibrated by boreholes was used to calculate water-loaded subsidence histories at rollover points of clinoforms along the shelf. At 9u2009±u20093u2009Ma, the rate of subsidence increases from 5 to up 75u2009mu2009Myr−1, generating a subsidence anomaly of −300 to −700u2009m. The amplitude of this anomaly varies along the shelf and cannot be generated by glacio-eustatic sea-level variation. Instead, we propose that a regional subsidence episode, which affects both the proximal shelf and the distal oceanic basin, was generated by convective drawdown. By combining our results with other published estimates of uplift and subsidence, a map of Australia, which shows the spatial and temporal pattern of dynamic topography is presented. Most, but not all, of Australias epeirogeny can be attributed to rapid northward motion of the Australian plate over a pre-existing pattern of convective circulation.

Cenozoic epeirogeny of the Arabian Peninsula from drainage modeling

It is generally accepted that the Arabian Peninsula has been uplifted by subcrustal processes. Positive residual depth anomalies from oceanic crust in the Red Sea and in the Gulf of Aden suggest that a region surrounding this peninsula is dynamically supported. Admittance calculations, surface wave tomography studies, and receiver function analyses all imply that regional topography is generated and maintained by some combination of mantle convective circulation and lithospheric thickness changes. Despite these significant advances, the spatial and temporal uplift rate history of the Arabian Peninsula is not well known. Here we show that a regional uplift rate history can be obtained by jointly inverting 225 longitudinal river profiles that drain this peninsula. Our strategy assumes that shapes of individual river profiles are controlled by uplift rate history and moderated by erosional processes. We used local measurements of incision rate to calibrate the relevant erosional parameters. In our inverse algorithm, uplift rate is permitted to vary smoothly as a function of space and time but upstream drainage area remains invariant. We also assume that knickzone migration is not lithologically controlled. Implications of these important assumptions have been investigated. Our results suggest that the Arabian Peninsula underwent two phases of asymmetric uplift during the last 20–30 Ma at rates of 0.05–0.1 mm a−1. The southwestern flank of the peninsula has been uplifted by 1.5–2.5 km. Regional stratigraphic constraints, the age and composition of volcanism, paleosol formation, incised peneplains, emergent marine terraces, and thermochronometric measurements corroborate our calculated patterns of uplift. Progressive development of three domal swells along the western margin of the peninsula is consistent with localized upwelling of hot asthenospheric mantle.

Oceanic Residual Depth Measurements, the Plate Cooling Model and Global Dynamic Topography

Convective circulation of the mantle causes deflections of the Earths surface that vary as a function of space and time. Accurate measurements of this dynamic topography are complicated by the need to isolate and remove other sources of elevation, arising from flexure and lithospheric isostasy. The complex architecture of continental lithosphere means that measurement of present-day dynamic topography is more straightforward in the oceanic realm. Here, we present an updated methodology for calculating oceanic residual bathymetry, which is a proxy for dynamic topography. Corrections are applied that account for the effects of sedimentary loading and compaction, for anomalous crustal thickness variations, for subsidence of oceanic lithosphere as a function of age, and for non-hydrostatic geoid height variations. Errors are formally propagated to estimate measurement uncertainties. We apply this methodology to a global database of 1,936 seismic surveys located on oceanic crust and generate 2,297 spot measurements of residual topography, including 1,161 with crustal corrections. The resultant anomalies have amplitudes of ±1 km and wavelengths of ∼1,000 km. Spectral analysis of our database using cross-validation demonstrates that spherical harmonics up to and including degree 30 (i.e. wavelengths down to 1,300 km) are required to accurately represent these observations. Truncation of the expansion at a lower maximum degree erroneously increases the amplitude of inferred long-wavelength dynamic topography. There is a strong correlation between our observations and free-air gravity anomalies, magmatism, ridge seismicity, vertical motions of adjacent rifted margins, and global tomographic models. We infer that shorter wavelength components of the observed pattern of dynamic topography may be attributable to the presence of thermal anomalies within the shallow asthenospheric mantle.

A Cenozoic Uplift History of Mexico and its Surroundings From Longitudinal River Profiles

Geodynamic models of mantle convection predict that Mexico and western North America share a history of dynamic support. We calculate admittance between gravity and topography, which indicates that the elastic thickness of the plate in Mexico is 11 km and in western North America it is 12 km. Admittance at wavelengthsu2009>u2009500 km in these regions suggests that topography is partly supported by subcrustal processes. These results corroborate estimates of residual topography from isostatic calculations and suggest that the amount of North American topography supported by the mantle may exceed 1 km. The Cenozoic history of magmatism, sedimentary flux, thermochronometric denudation estimates, and uplifted marine terraces imply that North American lithosphere was uplifted and eroded during the last 30 Ma. We jointly invert 533 Mexican and North American longitudinal river profiles to reconstruct a continent-scale rock uplift rate history. Uplift rate is permitted to vary in space and time. Erosional parameters are calibrated using incision rate data in southwest Mexico and the Colorado Plateau. Calculated rock uplift rates were 0.15–0.2 mm/yr between 25 and10 Ma. Central Mexico experienced the highest uplift rates. Central and southern Mexico continued to uplift at 0.1 mm/yr until recent times. This uplift history is corroborated by independent constraints. We predict clastic flux to the Gulf of Mexico and compare it to independent estimates. We tentatively suggest that the loop between uplift, erosion, and deposition can be closed here. Mexicos staged uplift history suggests that its dynamic support has changed during the last 30 Ma.

Cenozoic Epeirogeny of the Indian Peninsula

Peninsular India is a cratonic region with asymmetric relief manifest by eastward tilting from the 1.5 km high Western Ghats escarpment toward the floodplains of eastward-draining rivers. Oceanic residual depth measurements on either side of India show that this west-east asymmetry is broader scale, occurring over distances of >2,000 km. Admittance analysis of free-air gravity and topography shows that the elastic thickness is 10 ±3 km, suggesting that regional uplift is not solely caused by flexural loading. To investigate how Indian physiography is generated, we have jointly inverted 530 river profiles to determine rock uplift rate as a function of space and time. Key erosional parameters are calibrated using independent geologic constraints (e.g. emergent marine deposits, elevated paleosurfaces, uplifted lignite deposits). Our results suggest that regional tilt grew at rates of up to 0.1 mm a– 1 between 25 Ma and the present day. Neogene uplift initiated in the south and propagated northward along the western margin. This calculated history is corroborated by low-temperature thermochronologic observations, by sedimentary flux of clastic deposits into the Krishna-Godavari delta, and by sequence stratigraphic architecture along adjacent rifted margins. Onset of regional uplift predates intensification of the Indian monsoon at 8 Ma, suggesting that rock uplift rather than climatic change is responsible for modern-day relief. A positive correlation between residual depth measurements and shear wave velocities beneath the lithosphere suggests that regional uplift is generated and maintained by temperature anomalies of ±100°C within a 200 ±25 km thick asthenospheric channel. This article is protected by copyright. All rights reserved.

Spatial and temporal uplift history of South America from calibrated drainage analysis

A multi-disciplinary approach is used to analyze the Cenozoic uplift history of South America. Residual depth anomalies of oceanic crust abutting this continent help to determine the pattern of present-day dynamic topography. Admittance analysis and crustal thickness measurements indicate that the elastic thickness of the Borborema and Altiplano regions isu2009≤u200910 km with evidence for sub-plate support at longer wavelengths. A drainage inventory of 1827 river profiles is assembled and used to investigate landscape development. Linear inverse modeling enables river profiles to be fitted as a function of the spatial and temporal history of regional uplift. Erosional parameters are calibrated using observations from the Borborema Plateau and tested against continent-wide stratigraphic and thermochronologic constraints. Our results predict that two phases of regional uplift of the Altiplano plateau occurred in Neogene times. Regional uplift of the southern Patagonian Andes also appears to have occurred in Early Miocene times. The consistency between observed and predicted histories for the Borborema, Altiplano and Patagonian plateaux implies that drainage networks record coherent signals that are amenable to simple modeling strategies. Finally, the predicted pattern of incision across the Amazon catchment constrains solid sedimentary flux at the Foz do Amazonas. Observed and calculated flux estimates match, suggesting that erosion and deposition were triggered by regional Andean uplift during Miocene times.

Neogene Uplift and Magmatism of Anatolia: Insights From Drainage Analysis and Basaltic Geochemistry

It is generally agreed that mantle dynamics have played a significant role in generating and maintaining the elevated topography of Anatolia during Neogene times. However, there is much debate about the relative importance of subduction zone and asthenospheric processes. Key issues concern onset and cause of regional uplift, thickness of the lithospheric plate, and the presence or absence of temperature and / or compositional anomalies within the convecting mantle. Here, we tackle these interlinked issues by analyzing and modeling two disparate suites of observations. First, a drainage inventory of 1,844 longitudinal river profiles is assembled. This geomorphic database is inverted to calculate the 18 nvariation of Neogene regional uplift through time and space by minimizing the misfit between observed and calculated river profiles subject to independent calibration. Our results suggest that regional uplift commenced in the east and propagated westward. Secondly, we have assembled a database of geochemical analyses of basaltic rocks. Two di ff erent approaches have been used to quantitatively model this database with a view to determining the depth and degree of asthenospheric melting across Anatolia. Our results suggest nthat melting occurs at depths as shallow as 60 km in the presence of mantle potential temperatures as high as 1390◦ nC. There is evidence that potential temperatures are higher in the east, consistent with the pattern of sub-plate shear wave velocity anomalies. Our combined results are consistent with isostatic and admittance analyses and suggest that elevated asthenospheric temperatures beneath thinned Anatolian lithosphere have played a first order role in generating and maintaining regional dynamic topography and basaltic magmatism.

Research data supporting "Neogene Uplift and Magmatism of Anatolia: Insights from Drainage Analysis and Basaltic Geochemistry"

Crustal thickness, drainage inventory and igneous geochemical data for Anatolia. Dataset S1 is compilation of published receiver function estimates of Anatolian crustal thickness. Station names and locations, number of receiver functions modeled, estimated Moho depth, and technique employed are all given where available. Note that elevations are extracted from topographic grid that is low-pass filtered for wavelengths > 30 km. Dataset S2 is drainage inventory as a multi-segment .csv file. Format is: latitude (◦),longitude (◦), elevation (m), longitudinal distance (m) and upstream drainage area (m^2). River profiles were extracted from SRTM digital elevation model of using Arc-GIS flow-routing algorithms. Dataset S3 is compilation of all available geochemical analyses for Neogene magmatic rocks from Anatolia. Ages are assigned either by radiometric dating or by stratigraphic relationships. Although only samples with MgO > 5 wt% were used in this study, samples with MgO > 4 wt% are included for completeness. Major elements are normalized so that the sum of weights = 100 wt%. We also provide estimates of melt equilibration pressure and temperature.

Research data supporting "Oceanic Residual Depth Measurements, the Plate Cooling Model and Global Dynamic Topography"

CSV files containing the measurements of oceanic residual depth anomalies from the paper. These data are a proxy for dynamic topography. Circles represent the most accurate data. Downwards pointing triangles are maximum constraints and upwards are minimums. Data have been averaged within 1 and 4 degree bins on the Earths surface. Half of the points include a correction for the shape of the Earths non-hydrostatic geoid. Gridfiles of the oceanic age grid and ship-track derived residual topography compatible with GMT (Generic Mapping Tools) are also included.