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

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

Geodynamic evolution of the lithosphere and upper mantle beneath the Alboran region of the western Mediterranean: Constraints from travel time tomography

An edited version of this paper was published by the American Geophysical Union. Copyright 2000, AGU. See also: http://www.agu.org/pubs/crossref/2000/2000JB900024.shtml; http://atlas.geo.cornell.edu/morocco/publications/calvert2000.htmA number of different geodynamic models have been proposed to explain the extension that occurred during the Miocene in the Alboran Sea region of the western Mediterranean despite the continued convergence and shortening of northern Africa and southern Iberia. In an effort to provide additional geophysical constraints on these models, we performed a local, regional, and teleseismic tomographic travel time inversion for the lithospheric and upper mantle velocity structure and earthquake locations beneath the Alboran region in an area of 800 x 800 km^2. We picked P and S arrival times from digital and analog seismograms recorded by 96 seismic stations in Morocco and Spain between 1989 and 1996 and combined them with arrivals carefully selected from local and global catalogs (1964-1998) to generate a starting data set containing over 100,000 arrival times. Our results indicate that a N-S line of intermediate depth earthquakes extending from crustal depths significantly inland from the southern Iberian coat to depths of over 100 km beneath the center of the Alboran Sea coincided with a W to E transition from high to low velocities imaged in the uppermost mantle. A high-velocity body, striking approximately NE-SW, is imaged to dip southeastwards from lithospheric depths beneath the low-velocity region to depths of ~350 km. Between 350 and 500 km the imaged velocity anomalies become more diffuse. However, pronounced high-velocity anomalies are again imaged at 600 km near an isolated cluster of deep earthquakes. In addition to standard tomographic methods of error assessment, the effects of systematic and random errors were assessed using block shifting and bootstrap resampling techniques, respectively. We interpret the upper mantle high-velocity anomalies as regions of colder mantle that originate from lithospheric depths. These observations, when combined with results from other studies, suggest that delamination of a continental lithosphere played an important role in the Neogene and Quaternary evolution of the region

Propagation of regional seismic phases (Lg and Sn) and Pn velocity structure along the Africa-Iberia plate boundary zone

An edited version of this paper was published by Blackwell Publishing. Copyright 2000, Blackwell Publishing. See also: http://www.blackwell-synergy.com/doi/abs/10.1046/j.1365-246x.2000.00160.x; http://atlas.geo.cornell.edu/morocco/publications/calvert2000GJI.htmWe used over 1000 regional waveforms recorded by 60 seismic stations located in northwest Africa and Iberia to map the efficiency of L g and Sn wave propagation beneath the Gulf of Cadiz, Alboran Sea and bounding Betic, Rif and Atlas mountain belts. Crustal attenuation is inferred from the tomographic inversion of L g/Pg amplitude ratios. Upper mantle attenuation is inferred from maps of Sn propagation efficiency derived by inversion of well-defined qualitative efficiency assignments based on waveform characteristics. Regions of L g attenuation correlate well with areas of thinned continental or oceanic crust, significant sedimentary basins, and lateral crustal variations. Comparison of the Sn efficiency results with velocities obtained from an anisotropic Pn traveltime inversion shows a fairly good correlation between regions of poor Sn efficiency and low Pn velocity. A low Pn velocity (7.6?7.8 km s-1) and significant Sn attenuation in the uppermost mantle is imaged beneath the Betics in southern Spain, in sharp contrast to the relatively normal Pn velocity (8.0?8.1 km s-1) and efficient Sn imaged beneath the Alboran Sea. Slow Pn velocity anomalies are also imaged beneath the Rif and Middle Atlas in Morocco. We do not identify any conclusive evidence of lithospheric-scale upper mantle attenuation beneath the Rif, although the crust in the Gibraltar region appears highly attenuating, making observations at stations in this region ambiguous. Paths crossing the Gulf of Cadiz, eastern Atlantic and the Moroccan and Iberian mesetas show very efficient Sn propagation and are imaged with high Pn velocities (8.1?8.2 km s-1). The spatial distribution of attenuation and velocity anomalies lead us to conclude that some recovery of the mantle lid beneath the Alboran Sea must have occurred since the early Miocene episode of extension and volcanism. We interpret the low velocity and attenuating regions beneath the Betics and possibly the Rif as indicating the presence of partial melt in the uppermost mantle which may be underlain by faster less attenuating mantle. In the light of observations from other geophysical and geological studies, the presence of melt at the base of the Betic crust may be an indication that delamination of continental lithosphere has played a role in the Neogene evolution of the Alboran Sea region

Imaging the lithosphere beneath NE Tibet: Teleseismic P and S body wave tomography incorporating surface wave starting models

The northeastern margin of the Tibetan Plateau, which includes the Qiangtang and Songpan-Ganzi terranes as well as the Kunlun Shan and the Qaidam Basin, continues to deform in response to the ongoing India–Eurasia collision. To test competing hypotheses concerning the mechanisms for this deformation, we assembled a high-quality data set of approximately 14 000 P- and 4000 S-wave arrival times from earthquakes at teleseismic distances from the International Deep Profiling of Tibet and the Himalaya, Phase IV broad-band seismometer deployments. We analyse these arrival times to determine tomographic images of P- and S-wave velocities in the upper mantle beneath this part of the plateau. To account for the effects of major heterogeneity in crustal and uppermost mantle wave velocities in Tibet, we use recent surface wave models to construct a starting model for our teleseismic body wave inversion. We compare the results from our model with those from simpler starting models, and find that while the reduction in residuals and results for deep structure are similar between models, the results for shallow structure are different. Checkerboard tests indicate that features of ~125km length scale are reliably imaged throughout the study region. Using synthetic tests, we show that the best recovery is below ~300km, and that broad variations in shallow structure can also be recovered. We also find that significant smearing can occur, especially at the edges of the model. We observe a shallow dipping seismically fast structure at depths of ~140–240km, which dies out gradually between 33°N and 35°N. Based on the lateral continuity of this structure (from the surface waves) we interpret it as Indian lithosphere. Alternatively, the entire area could be thickened by pure shear, or the northern part could be an underthrust Lhasa Terrane lithospheric slab with only the southern part from India. We see a deep fast wave velocity anomaly (below 300?km), that is consistent with receiver function observations of a thickened transition zone and could be a fragment of oceanic lithosphere. In NE Tibet, it appears to be disconnected from faster wave velocities above (i.e. it is not downwelling or subducting here). Our models corroborate results of previous work which imaged a relatively slow wave velocity region below the Kunlun Shan and northern Songpan-Ganzi Terrane, which is difficult to reconcile with the hypothesis of southward-directed continental subduction at the northern margin. Wave velocities in the shallow mantle beneath the Qaidam Basin are faster than normal, and more so in the east than the west.This work was supported by a Natural Environment Research Council studentship (grant NE/H52449X/1)This version of record of this article can be found in Geophysical Journal International (March, 2014) 196 (3): 1724-1741. doi: 10.1093/gji/ggt47

Tomographic Pn velocity and anisotropy structure beneath the Anatolian plateau (eastern Turkey) and the surrounding regions

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.../2003GL017391.shtml; http://atlas.geo.cornell.edu/turkey/publications/Al-Lazki-et-al_2003.htmWe use Pn phase travel time residuals to invert for mantle lid velocity and anisotropy beneath northern Arabia eastern Anatolia continent-continent collision zone. The primary phase data were obtained from the temporary 29- station broadband PASSCAL array of the Eastern Turkey Seismic Experiment. These data were supplemented by phase data from available stations of the Turkish National Seismic Network, the Syrian National Seismic Network, the Iranian Long Period Array, and other stations around the southern Caspian Sea. In addition, we used carefully selected catalog data from the International Seismological Centre and the National Earthquake Information Center bulletins. Our results show that low (<8 km/s) to very low (<7.8 km/s) Pn velocity zones underlie the Anatolian plateau, the Caucasus, and northwestern Iran. Such low velocities are used to infer the presence of partially molten to absent mantle lid beneath these regions. In contrast, we observed a high Pn velocity zone beneath northern Arabia directly south of the Bitlis-Zagros suture indicating the presence of a stable Arabian mantle lid. This sharp velocity contrast across the suture zone suggests that Arabia is not underthrusting beneath the Anatolian plateau and that the surface suture extends down to the uppermost mantle. Pn anisotropy orientations within a single plate (e.g. Anatolia plate) show a higher degree of lateral variation compared to Pn velocity. Areas of coherent Pn anisotropy orientations are observed to continue across major fault zones such as the EAF zone

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

Structure and tectonic evolution of the Anatolian plateau in eastern Turkey

This paper was published by the Geological Society of America (GSA). Copyright 2006, GSA. See also: http://granite.geosociety.org/bookstore/default.asp?oID=0&catID=9&pID=SPE409; http://atlas.geo.cornell.edu/turkey/publications/Barazangi-et-al_2006.htmThe Cenozoic geology and the present lithospheric and upper mantle structure of the Anatolian plateau in eastern Turkey and nearby regions are the result of the final collision and suturing of the continental Arabian plate to the Turkish terranes (i.e., micro-continents). This process of collision and suturing was strongly influenced by three active structures in the region: the Caucasus mountains, the Aegean subduction zone, and the Dead Sea fault system. Understanding these three major tectonic elements are important for the development of a robust model for the formation of the Anatolian plateau. We show that the Anatolian plateau lithosphere in eastern Turkey has no lithospheric mantle, i.e., the crust floats on a partially molten asthenosphere. The average thickness of the crust in the region is approximately 45 km. The uppermost mantle beneath this crustal block strongly attenuates Sn waves and has one of the lowest Pn velocities on earth (about 7.6 km/s). The Anatolian plateau, with an average of 2 km elevation is dissected by numerous active, seismogenic faults (mostly strike-slip and some thrust type). Neogene and Quaternary volcanism with varying composition is widespread and covers more than half of the region. We argue that the northward subduction of the northern and the southern branches of the Neo-Tethyan oceanic lithosphere since the Mesozoic has resulted in the development of arc and back-arc volcanism (i.e., the Pontide and Bitlis systems) and the development of the eastern Anatolian accretionary complex that covers a large area of eastern Turkey. The northward subduction of the southern Neo-Tethys considerably thinned and weakened the overriding Eurasian plate above the descending oceanic lithosphere of the Arabian plate. The final suturing of the continental Arabian plate with the Turkish terranes in the Miocene and the continued convergence of Arabia relative to Eurasia has resulted in the shortening of the accretionary complex both in the forearc and the back-arc regions and the development of a broad zone with numerous strike-slip faults. The mobilization of the Caucasus is also partially a consequence of this convergence. The documented major episode of widespread volcanism at about 11 Ma is probably related to the breakoff of the shallowly descending oceanic segment of the Arabian lithosphere beneath eastern Turkey. The continued convergence of Arabia relative to Eurasia has resulted in the development of the North Anatolian fault (NAF) and subsequently the East Anatolian fault (EAF) in the Pliocene. At about this time, the northern segment of the Dead Sea fault (DSF) also developed in Lebanon and northwest Syria and joined the EAF to form the Anatolian - Arabian - African triple junction in the Maras region of southern Turkey. The development of these fault systems (i.e., NAF, EAF, and DSF) provided the mechanism for the tectonic escape of the Anatolian crustal block towards the Aegean arc system

Non-monotonicity of the frictional bimaterial effect

Sliding along frictional interfaces separating dissimilar elastic materials is qualitatively different from sliding along interfaces separating identical materials due to the existence of an elastodynamic coupling between interfacial slip and normal stress perturbations in the former case. This bimaterial coupling has important implications for the dynamics of frictional interfaces, including their stability and rupture propagation along them. We show that while this bimaterial coupling is a monotonically increasing function of the bimaterial contrast, when it is coupled to interfacial shear stress perturbations through a friction law, various physical quantities exhibit a non-monotonic dependence on the bimaterial contrast. In particular, we show that for a regularized Coulomb friction, the maximal growth rate of unstable interfacial perturbations of homogeneous sliding is a non-monotonic function of the bimaterial contrast, and provide analytic insight into the origin of this non-monotonicity. We further show that for velocity-strengthening rate-and-state friction, the maximal growth rate of unstable interfacial perturbations of homogeneous sliding is also a non-monotonic function of the bimaterial contrast. Results from simulations of dynamic rupture along a bimaterial interface with slip-weakening friction provide evidence that the theoretically predicted non-monotonicity persists in non-steady, transient frictional dynamics.Comment: 14 pages, 5 figure