NASA Geodynamics Program

Activities and achievements for the period of May 1983 to May 1984 for the NASA geodynamics program are summarized. Abstracts of papers presented at the Conference are inlcuded. Current publications associated with the NASA Geodynamics Program are listed

The crustal structure of the Anatolian Plate from receiver functions and implications for the uplift of the Central and Eastern Anatolian plateaus

Understanding the crustal structure of the Anatolian Plate has important implications for its formation and evolution, including the extent to which its high elevation is maintained isostatically. However, the numerous teleseismic receiver function studies from which Anatolian Moho depths have been obtained return results that differ by ≤21 km at some seismograph stations. To address this issue, we determine Moho depth and bulk crustal VP/VS ratio (κ) at 582 broadband seismograph stations, including ∼100 for which H-κ results have not been reported previously. We use a modified H-κ stacking method in which a final solution is selected from a suite of up to 1000 repeat H-κ measurements, each calculated using randomly-selected receiver functions and H-κ input parameters. Ten quality control criteria that variously assess the final numerical result, the receiver function data set, and the extent to which the results are clustered tightly, are used to determine station quality. By refining Moho depth constraints, including identifying 182 stations, analysed previously, where H-κ stacking yields unreliable results (particularly in Eastern Anatolia and the rapidly-uplifting Taurides), our new crustal model (ANATOLIA-HK21) provides fresh insight into Anatolian crustal structure and topography. Changes in Moho depth within the Anatolian Plate occur on a shorter length-scale than has sometimes previously been assumed. For example, crustal thickness decreases abruptly from >40 km in the northern Kirsehir block to <32 km beneath the Central Anatolian Volcanic Province and Tuz Golu basin. Moho depth increases from 30-35 km on the Arabian Plate to 35-40 km across the East Anatolian Fault into Anatolia, in support of structural geological observations that Arabia-Anatolia crustal shortening was accommodated primarily on the Anatolian, not Arabian, Plate. However, there are no consistent changes in Moho depth across the North Anatolian Fault, whose development along the Intra-Pontide and İzmir-Ankara-Erzincan suture zones was more likely the result of contrasts in mantle lithospheric, not crustal, structure. While the crust thins from ∼45 km below the uplifted Eastern Anatolian Plateau to ∼25 km below lower-lying western Anatolia, Moho depth is generally correlated poorly with elevation. Residual topography calculations confirm the requirement for a mantle contribution to Anatolian Plateau uplift, with localised asthenospheric upwellings in response to slab break-off and/or lithospheric dripping/delamination example candidate driving mechanisms

Cenozoic Epeirogeny of 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

Crustal seismic structure of the eastern Mediterranean: evidence from broadband seismology

Understanding the crustal structure of the Anatolian Plate has important implications for its formation and evolution, including the extent to which its high elevation is maintained isostatically. However, previous receiver function studies of Anatolian Moho depths return results differing by <21km. H-K stacking is used routinely to infer crustal thickness and Vp/Vs ratio from teleseismic receiver functions, assuming the largest amplitude P-to-S conversions beneath the seismograph station are generated by a sharp velocity contrast at the Moho. However, synthetic seismogram analysis demonstrates that H-K results are strongly dependent on the choice of stacking input parameters. To address this issue of parameter sensitivity, an H-K approach is developed in which cluster analysis selects a final solution from 1000 results, each calculated using randomly selected input parameters via bootstrapping. When the Moho is sharp, H-K results cluster tightly and the method is reliable; in areas of more complex crustal structure, H-K analysis is often unreliable. The new crustal model for the Eastern Mediterranean (ANATOLIA-HK21) provides fresh insight into Anatolian crustal structure and topography. While the crust thins from ~45km below the uplifted Eastern Anatolian Plateau to ~25km below lower-lying western Anatolia, Moho depth is generally correlated poorly with elevation. Residual topography calculations confirm the requirement for a mantle contribution to Anatolian Plateau uplift, with localised asthenospheric upwellings in response to slab break-off and/or lithospheric dripping/delamination example candidate driving mechanisms. In the absence of good-quality H-K results on Cyprus, Rayleigh wave group velocity inversions place new constraints on the island's crustal structure. Resulting tomographic images reveal high velocities at short periods directly beneath the surface expression of the Troodos Ophiolite; anomaly amplitudes decrease at longer periods (from +30% at 8s to +5% at 14s) and shift northeastward, corroborating studies that consider the ophiolite a mid-to-upper crustal feature.Open Acces

Lateral variation of crust and upper mantle structures in NW Iran derived from surface wave analysis

To obtain the shear velocity structure across North-West of Iran and surrounding areas to a depth of 160 km, we performed a namely Hedgehog nonlinear inversion on Rayleigh wave group velocity dispersion curves in the period range from 7 to 60 s. The distributed dispersion curves are the results of our surface waves dispersion tomography using the data of 280 local and regional seismic events, recorded by the medium and broad band seismic stations in the region. We outline different crust and upper mantle structures for the study area based on calculated group and shear velocities. Our results reveal relatively low velocities at the shorter periods (7 - 10 s) in the presence of sedimentary basins (e.g. South Caspian Basin) and for eastern Anatolia and relatively high velocities along the Sanandaj-Sirjan Metamorphic zone, Alborz, Talesh and the Lesser Caucasus Mountains. By depth inversion of group velocities, we observed a 14 km thick sediments in South Caspian Basin and Kura Depression. Based on our maps at 20 s, we outline different crustal models for the region and highlight the differences between South Caspian Basin and NW Iran, on one side, and the similarities between the South Caspian Basin and Kura Depression, that extend beneath Talesh, Alborz and Lesser Caucasus, on the other. Comparing the shear velocity of lower crust in South Caspian Basin and Kura Depression with that of NW Iran proves different origination of lower crust in the basin, probably oceanic source, because of its significant higher shear velocity rather than NW Iran. The extension of lower crust beneath Talesh is more than middle crust while in Alborz and Lesser Caucasus the amount of extension for middle and lower crust is the same The analysis of group velocities at longer periods ( 65 35 s) and obtained shear velocity models allows us to outline different lithospheric structures and crustal depth in the region. The high group velocities in Talesh, South Caspian Sea and Lesser Caucasus on one side and Zagros Folding and Thrust Belt on the other, beside the result of shear velocity models suggest the presence of a stable and thick mantle lid that seems to be thin or absent in the eastern Anatolia and much of NW Iran. The shallowest Moho and Lithosphere Asthenosphere boundary depth of 37 and 63 km, were observed in Easter Anatolian Accretionary Complex. The thin mantle lid in this region has affected the whole crust in such a way that we observed the lowest shear velocities inside the crust in this region. We observed a significant thickening of both crust and lithosphere in Sanandaj-Sirjan Metamorphic zone comparing to Urmieh Dokhtar Magmatic Arc and Zagros Folding and Thrust Belt on its two sides

The 2008 Methoni earthquake sequence: the relationship between the earthquake cycle on the subduction interface and coastal uplift in SW Greece

Seismological, GPS and historical data suggest that most of the 40 mm yr−1 convergence at the Hellenic Subduction Zone is accommodated through aseismic creep, with earthquakes of MW ≲ 7 rupturing isolated locked patches of the subduction interface. The size and location of these locked patches are poorly constrained despite their importance for assessment of seismic hazard. We present continuous GPS time-series covering the 2008 MW 6.9 Methoni earthquake, the largest earthquake on the subduction interface since 1960. Post-seismic displacements from this earthquake at onshore GPS sites are comparable in magnitude with the coseismic displacements; elastic-dislocation modelling shows that they are consistent with afterslip on the subduction interface, suggesting that much of this part of the interface is able to slip aseismically and is not locked and accumulating elastic strain. In the Hellenic and other subduction zones, the relationship between earthquakes on the subduction interface and observed long-term coastal uplift is poorly understood. We use cGPS-measured coseismic offsets and seismological body-waveform modelling to constrain centroid locations and depths for the 2008 Methoni MW 6.9 and 2013 Crete MW 6.5 earthquakes, showing that the subduction interface reaches the base of the seismogenic layer SW of the coast of Greece. These earthquakes caused subsidence of the coast in regions where the presence of Pliocene–Quaternary marine terraces indicates recent uplift, so we conclude that deformation associated with the earthquake cycle on the subduction interface is not the dominant control on vertical motions of the coastline. It is likely that minor uplift on a short length scale (∼15 km) occurs in the footwalls of normal faults. We suggest, however, that most of the observed Plio-Quaternary coastal uplift in SW Greece is the result of thickening of the overriding crust of the Aegean by reverse faulting or distributed shortening in the accretionary wedge, by underplating of sediment of the Mediterranean seafloor, or a combination of these mechanisms

Thermal Perturbations beneath the Incipient Okavango Rift Zone, Northwest Botswana

We used aeromagnetic and gravity data to investigate the thermal structure beneath the incipient Okavango Rift Zone (ORZ) in northwestern Botswana in order to understand its role in strain localization during rift initiation. We used three-dimensional (3-D) inversion of aeromagnetic data to estimate the Curie Point Depth (CPD) and heat flow under the rift and surrounding basement. We also used two-dimensional (2-D) power-density spectrum analysis of gravity data to estimate the Moho depth. Our results reveal shallow CPD values (8-15 km) and high heat flow (60-90 mW m-2) beneath a ∼60 km wide NE-trending zone coincident with major rift-related border faults and the boundary between Proterozoic orogenic belts. This is accompanied by thin crust ( \u3c 30 km) in the northeastern and southwestern parts of the ORZ. Within the Precambrian basement areas, the CPD values are deeper (16-30 km) and the heat flow estimates are lower (30-50 mW m-2), corresponding to thicker crust (∼40-50 km). We interpret the thermal structure under the ORZ as due to upward migration of hot mantle fluids through the lithospheric column that utilized the presence of Precambrian lithospheric shear zones as conduits. These fluids weaken the crust, enhancing rift nucleation. Our interpretation is supported by 2-D forward modeling of gravity data suggesting the presence of a wedge of altered lithospheric mantle centered beneath the ORZ. If our interpretation is correct, it may result in a potential paradigm shift in which strain localization at continental rift initiation could be achieved through fluid-assisted lithospheric weakening without asthenospheric involvement

Thermal Perturbations beneath the Incipient Okavango Rift Zone, Northwest Botswana

We used aeromagnetic and gravity data to investigate the thermal structure beneath the incipient Okavango Rift Zone (ORZ) in northwestern Botswana in order to understand its role in strain localization during rift initiation. We used three-dimensional (3-D) inversion of aeromagnetic data to estimate the Curie Point Depth (CPD) and heat flow under the rift and surrounding basement. We also used two-dimensional (2-D) power-density spectrum analysis of gravity data to estimate the Moho depth. Our results reveal shallow CPD values (8-15 km) and high heat flow (60-90 mW m-2) beneath a ∼60 km wide NE-trending zone coincident with major rift-related border faults and the boundary between Proterozoic orogenic belts. This is accompanied by thin crust ( \u3c 30 km) in the northeastern and southwestern parts of the ORZ. Within the Precambrian basement areas, the CPD values are deeper (16-30 km) and the heat flow estimates are lower (30-50 mW m-2), corresponding to thicker crust (∼40-50 km). We interpret the thermal structure under the ORZ as due to upward migration of hot mantle fluids through the lithospheric column that utilized the presence of Precambrian lithospheric shear zones as conduits. These fluids weaken the crust, enhancing rift nucleation. Our interpretation is supported by 2-D forward modeling of gravity data suggesting the presence of a wedge of altered lithospheric mantle centered beneath the ORZ. If our interpretation is correct, it may result in a potential paradigm shift in which strain localization at continental rift initiation could be achieved through fluid-assisted lithospheric weakening without asthenospheric involvement