Design and development of information systems for the geosciences: An application to the Middle East

Publisher's version archived with permission from publisher. http://www.gulfpetrolink.net/publication/geoarabia.htmAs our understanding grows of how the Earth functions as a complex system of myriad interrelated mechanisms, it becomes clear that a revolutionary and novel approach is needed to study and understand it. In order to take advantage of an ever-growing number of observations and large data sets and to employ them efficiently in multidisciplinary studies aimed at solving earth system science problems, we are developing a comprehensive Solid Earth Information System (SEIS). The complex nature of the solid earth sciences raises serious challenges for geoscientists in their quest to understand the nature and the dynamic mechanisms at work in the planet. SEIS forms a first step in developing a broader and more comprehensive information system for earth system sciences designed for the needs of the geoscientists of the 21st century. In a way, SEIS is a step towards the Digital Earth. Application of SEIS to the complex tectonics of the Middle East shows that information systems are crucial in multidisciplinary research studies and open new avenues in research efforts. SEIS includes an Internet module that provides open access to anyone interested. Researchers as well as educators and students can access this knowledge and information system at http://atlas.geo.cornell.edu

Vertical constraint on mantle anisotropy from shear wave splitting in the Isparta Angle, Turkey [abstract]

Abstract only availableThe Isparta Angle in southwestern Turkey is the terrestrial expression of the off shore intersection of the Hellenic and Cyprian arcs. It has been suggested that there is a tear in the down-going lithospheric slab. If a tear does exist, it will be evidenced in the mantle flow beneath the angle. When shear seismic waves travel through the mantle they can become polarized in the direction parallel to flow. On a three component seismometer, the S-wave will be recorded earlier on the horizontal component corresponding to the polarization direction. The time between the polarized and non-polarized horizontal components and the direction of polarization are both computed and plotted on a map. We are interested in the polarization directions to determine the direction of flow beneath the Isparta Angle. We studied two different types of shear waves. We analyzed local events with a maximum depth of 250 km. We also analyzed SKS events in which shear waves travel through the entire mantle and into the outer core, meaning that the flow or anisotropy causing splitting could be anywhere in the mantle. By comparing our local results that are confined to the upper 250 km of the mantle and the SKS results we have found that there is a large jump in lag time somewhere beneath 250 km. We found the average lag time for local events to be about 0.6 seconds whereas the SKS average about 1.9 seconds. This means the majority of anisotropy is in the lower portions of the mantle. Data for this research is being collected by a temporary array of seismic stations deployed around the Isparta area. The project is a NSF funded collaboration of the University of Missouri-Columbia, Kandilli Observatory-Istanbul, and Suleyman Demirel University-Isparta.National Science Foundation; Summer Arts & Science Undergraduate Research Mentorship Progra

Shear-Wave Splitting and Mantle Flow Beneath the Colorado Plateau and its Boundary with the Great Basin

Shear-wave splitting measurements from SKS and SKKS phases show fast polarization azimuths that are subparallel to North American absolute plate motion within the central Rio Grande Rift (RGR) and Colorado Plateau (CP) through to the western rim of the CP, with anisotropy beneath the CP and central RGR showing a remarkably consistent pattern with a mean fast azimuth of 4 degrees +/- degrees 6 E of N. Approaching the rim from the southeast, fast anisotropic directions become north-northeast-south-southwest (NNE-SSW), rotate counter clockwise to north-south in the CP-GB transition, and then to NNW-SSE in the western Great Basin ( GB). This change is coincident with uppermost mantle S-wave velocity perturbations that vary from +4% beneath the western CP and the eastern edge of the Marysvale volcanic field to about -8% beneath the GB. Corresponding delay times average 1.5 sec beneath the central CP, decrease to approximately 0.8 sec near the CP-GB transition, and increase to about 1.2 sec beneath the GB. For the central CP, we suggest anisotropy predominantly controlled by North American plate motion above the asthenosphere. The observed pattern of westward-rotating anisotropy from the western CP through the CP-GB transition may be influenced to asthenospheric flow around a CP lithospheric keel and/or by vertical flow arising from edge-driven small-scale convection. The anisotropic transition from the CP to the GB thus marks a first-order change from absolute plate motion dominated lithosphere-asthenosphere shear to a new regime controlled by regional flow processes. The NNW-SSE anisotropic fast directions of split SKS waves in the eastern GB area are part of a broad circular pattern of seismic anisotropic fast direction in the central GB that has recently been hypothesized to be due to toroidal flow around the sinking Juan de Fuca-Gorda slab.National Science Foundation EAR 9706094, 9707188, 9707190, 0207812Los Alamos National Laboratory Institute of Geophysics and Planetary PhysicsNational Science Foundation Cooperative EAR-000430Department of Energy National Nuclear Security AdministrationGeological Science

High resolution regional seismic attenuation tomography in eastern Tibetan Plateau and adjacent regions

The Q of regional seismic phases Lg and Pg within the crust is assumed as a proxy for crustal Qβ and Qα, which is used as a constraint of crustal rheology. We measure regional‐phase Q of the eastern Tibetan Plateau and adjacent areas. This method eliminates contributions from source and site responses and is an improvement on the Two‐Station Method (TSM). We have generated tomographic images of crustal attenuation anomalies with resolution as high as 1°. In general we observe low Q in the northernmost portions of the Tibetan Plateau and high Q in the more tectonically stable regions such as the interior of the Qaidam basin. The calculated site responses appear to correlate with topography or sediment thickness. Furthermore the relationship between earthquake magnitudes and calculated source terms suggest that the RTM method effectively removes the source response and may be used as an alternative to source magnitude

Crustal Velocity Structure of the Northeastern Tibetan Plateau from Ambient Noise Surface-Wave Tomography and its Tectonic Implications

Broadband seismic data from the regional seismic network operated by the China Earthquake Administration and 32 temporary seismic stations are used to image the crustal velocity structure in the northeast Tibetan plateau. Empirical Rayleigh‐ and Love‐wave Green’s functions are obtained from interstation cross correlation of continuous seismic records. Group velocity dispersion curves for Rayleigh and Love waves between 10 and 50 s are obtained using the multiple‐filter analysis method with phase‐matched processing. The group velocity variations of Rayleigh and Love waves overall correlate well with the major geologic structures and tectonic units in the study region. Shear‐wave velocity structures were then inverted from Rayleigh‐ and Love‐wave dispersion maps. The results show that the Songpan–Ganzi terrane is associated with a low velocity at depth greater than 20 km. The northern Qilian orogen, with higher elevation and thicker crust compared to the southern Qilian orogen, is also dominated by low velocity at depth greater than ∼25  km. However, there is no clear evidence of the low‐velocity mid‐to‐lower crust beneath the southern Qilian orogen as the crustal flow model predicts. The low‐velocity zone (LVZ) beneath the northern Qilian orogen may suggest that the crustal thickening and surface uplift of the northern Qilian orogen are related to the LVZ, and the LVZ may be considered as an intracrustal response to bear the ongoing deformation in the northern Qilian orogen

Structure of the crust and African slab beneath the central Anatolian plateau from receiver functions: New insights on isostatic compensation and slab dynamics

The central Anatolian plateau in Turkey is a region with a long history of subduction, continental collision, accretion of continental fragments, and slab tearing and/or breakoff and tectonic escape. Central Anatolia is currently characterized as a nascent plateau with widespread Neogene volcanism and predominantly transtensional deformation. To elucidate the present-day crustal and upper mantle structure of this region, teleseismic receiver functions were calculated from 500 seismic events recorded on 92 temporary and permanent broadband seismic stations. Overall, we see a good correlation between crustal thickness and elevation throughout central Anatolia, indicating that the crust may be well compensated throughout the region. We observe the thickest crust beneath the Taurus Mountains (>40 km); it thins rapidly to the south in the Adana Basin and Arabian plate and to the northwest across the Inner Tauride suture beneath the Tuz Gölü Basin and Kırşehir block. Within the Central Anatolian Volcanic Province, we observe several low seismic velocity layers ranging from 15 to 25 km depth that spatially correlate with the Neogene volcanism in the region, and may represent crustal magma reservoirs. Beneath the central Taurus Mountains, we observe a positive amplitude, subhorizontal receiver function arrival below the Anatolian continental Moho at ∼50–80 km that we interpret as the gently dipping Moho of the subducting African lithosphere abruptly ending near the northernmost extent of the central Taurus Mountains. We suggest that the uplift of the central Taurus Mountains (∼2 km since 8 Ma), which are capped by flat-lying carbonates of late Miocene marine units, can be explained by an isostatic uplift during the late Miocene–Pliocene followed by slab breakoff and subsequent rebound coeval with the onset of faster uplift rates during the late Pliocene–early Pleistocene. The Moho signature of the subducting African lithosphere terminates near the southernmost extent of the Central Anatolian Volcanic Province, where geochemical signatures in the Quaternary volcanics indicate that asthenospheric material is rising to shallow mantle depths

Crustal structure of the Arabian Plate: New constraints from the analysis of teleseismic receiver functions

An edited version of this paper was published by Elsevier Science. Copyright 2005, Elsevier Science. See also: http://dx.doi.org/10.1016/j.epsl.2004.12.020; http://atlas.geo.cornell.edu/SaudiArabia/publications/Al-Damegh%202005.htmReceiver functions for numerous teleseismic earthquakes recorded at 23 broadband and mid-band stations in Saudi Arabia and Jordan were analyzed to map crustal thickness within and around the Arabian plate. We used spectral division as well as time domain deconvolution to compute the individual receiver functions and receiver function stacks. The receiver functions were then stacked using the slant stacking approach to estimate Moho depths and Vp/Vs for each station. The errors in the slant stacking were estimated using a bootstrap re-sampling technique. We also employed a grid search waveform modeling technique to estimate the crustal velocity structure for seven stations. A jackknife re-sampling approach was used to estimate errors in the grid search results for three stations. In addition to our results, we have also included published receiver function results from two temporary networks in the Arabian shield and Oman as well as three permanent GSN stations in the region. The average crustal thickness of the late Proterozoic Arabian shield is 39 km. The crust thins to about 23 km along the Red Sea coast and to about 25 km along the margin of the Gulf of Aqaba. In the northern part of the Arabian platform, the crust varies from 33 to 37 km thick. However, the crust is thicker (41?53 km) in the southeastern part of the platform. There is a dramatic change in crustal thickness between the topographic escarpment of the Arabian shield and the shorelines of the Red Sea. We compared our results in the Arabian shield to nine other Proterozoic and Archean shields that include reasonably well determined Moho depths, mostly based on receiver functions. The average crustal thickness for all shields is 39 km, while the average for Proterozoic shields is 40 km, and the average for Archean shields is 38 km. We found the crustal thickness of Proterozoic shields to vary between 33 and 44 km, while Archean shields vary between 32 and 47 km. Overall, we do not observe a significant difference between Proterozoic and Archean crustal thickness. We observed a dramatic change in crustal thickness along the Red Sea margin that occurs over a very short distance. We projected our results over a cross-section extending from the Red Sea ridge to the shield escarpment and contrasted it with a typical Atlantic margin. The transition from oceanic to continental crust of the Red Sea margin occurs over a distance of about 250 km, while the transition along a typical portion of the western Atlantic margin occurs at a distance of about 450 km. This important new observation highlights the abruptness of the breakup of Arabia. We argue that a preexisting zone of weakness coupled with anomalously hot upper mantle could have initiated and expedited the breakup

Crustal and uppermost mantle shear-wave velocity structure beneath the Middle East from surface-wave tomography

We have constructed a 3-D shear-wave velocity (Vs) model for the crust and uppermost mantle beneath the Middle East using Rayleigh wave records obtained from ambient-noise cross-correlations and regional earthquakes. We combined one decade of data collected from 852 permanent and temporary broadband stations in the region to calculate group-velocity dispersion curves. A compilation of > 54000 ray paths provides reliable group-velocity measurements for periods between 2 and 150 s. Path-averaged group velocities calculated at different periods were inverted for 2-D group-velocity maps. To overcome the problem of heterogeneous ray coverage, we used an adaptive grid parametrization for the group-velocity tomographic inversion. We then sample the period-dependent group-velocity field at each cell of a predefined grid to generate 1-D group-velocity dispersion curves, which are subsequently inverted for 1-D Vs models beneath each cell and combined to approximate the 3-D Vs structure of the area. The Vs model shows low velocities at shallow depths (5–10 km) beneath the Mesopotamian foredeep, South Caspian Basin, eastern Mediterranean and the Black Sea, in coincidence with deep sedimentary basins. Shallow high-velocity anomalies are observed in regions such as the Arabian Shield, Anatolian Plateau and Central Iran, which are dominated by widespread magmatic exposures. In the 10–20 km depth range, we find evidence for a band of high velocities (> 4.0 km/s) along the southern Red Sea and Arabian Shield, indicating the presence of upper mantle rocks. Our 3-D velocity model exhibits high velocities in the depth range of 30–50 km beneath western Arabia, eastern Mediterranean, Central Iranian Block, South Caspian Basin and the Black Sea, possibly indicating a relatively thin crust. In contrast, the Zagros mountain range, the Sanandaj-Sirjan metamorphic zone in western central Iran, the easternmost Anatolian plateau and Lesser Caucasus are characterized by low velocities at these depths. Some of these anomalies may be related to thick crustal roots that support the high topography of these regions. In the upper mantle depth range, high-velocity anomalies are obtained beneath the Arabian Platform, southern Zagros, Persian Gulf and the eastern Mediterranean, in contrast to low velocities beneath the Red Sea, Arabian Shield, Afar depression, eastern Turkey and Lut Block in eastern Iran. Our Vs model may be used as a new reference crustal model for the Middle East in a broad range of future studies

Shear wave splitting and mantle flow beneath LA RISTRA

Shear-wave splitting parameters (fast polarization direction and delay time) are determined using data from LA RISTRA (Colorado pLAteau RIo Grande Rift/Great Plains Seismic TRAnsect), a deployment of broadband seismometers extending from the Great Plains, across the Rio Grande Rift and the Jemez Lineament, to the Colorado Plateau. Results show that the fast polarization directions are sub-parallel to North American absolute plate motion. The largest deviations from the plate motion are observed within the western edge of the Great Plains and in the interior of the Colorado Plateau where lithospheric anisotropy may be significant. Delay times range from 0.8 to 1.8 seconds with an average value of 1.4 seconds; the largest values are along the Jemez Lineament and the Rio Grande Rift which are underlain by an uppermost mantle low velocity zone extending to depths of ∼200 km. The anisotropy beneath the central part of LA RISTRA shows a remarkably consistent pattern with a mean fast direction of 40° ± 6°. Seismic anisotropy can be explained by differential horizontal motion between the North American lithosphere and westerly to southwesterly flow of the asthenospheric mantle. The approximately N-S fast direction found beneath western Texas is similar to that observed beneath the southern rift and may reflect a different dynamic regime

The crustal structure of the East Anatolian plateau (Turkey) from receiver functions

An edited version of this paper was published by the American Geophysical Union (AGU). Copyright 2003, AGU. See also: http://www.agu.org/pubs/crossref/2003.../2003GL018192.shtml; http://atlas.geo.cornell.edu/turkey/publications/Zor-et-al_2003.htmThe crustal structure of the Anatolian plateau in Eastern Turkey is investigated using receiver functions obtained from the teleseismic recordings of a 29 broadband PASSCAL temporary network, i.e., the Eastern Turkey Seismic Experiment [ETSE]. The S-wave velocity structure was estimated from the stacked receiver functions by performing a 6-plane layered grid search scheme in order to model the first order features in the receiver functions with minimum trade-off. We found no significant crustal root beneath the western portion of the network, but there is some evidence of crustal thickening in the northern portion of the network. We found an average crustal thickness of 45 km and an average crustal shear velocity of 3.7 km/s for the entire eastern Anatolian plateau. Within the Anatolian plateau we found evidence of a prominent low velocity zone where the crust thickness is approximately 46 km. These results suggests that the 2 km high topography across the Anatolian plateau is dynamically supported because most of the plateau appears to be isostatically under-compensated. Also, there appears to be a region of thin crust at the easternmost edge of the Anatolian plateau that may be a relic from the accretion of island arcs to the Eurasian plate