Curie Point Depth Estimates from Aeromagnetic Data from Death Valley and Surrounding Regions, California

Aeromagnetic data were analyzed to determine the Curie point depth (CPD) by power density spectral and three-dimensional inversion methods within and surrounding Death Valley in southern California. We calculated the CPD for 0.5° regions using 2D power density spectral methods and found that the CPDs varied between 8 and 17 km. However, the 0.5° region may average areas that include shallow and deep CPDs, and because of this limitation, we used the 3D inversion method to determine if this method may provide better resolution of the CPDs. The final 3D model indicates that the depth to the bottom of the magnetic susceptible bodies varies between 5 and 23 km. Even though both methods produced roughly similar results, the 3D inversion method produced a higher lateral resolution of the CPDs. The shallowest CPDs occur within the central and southern Death Valley, Panamint Valley, Coso geothermal field and the Tecopa hot springs region. Deeper (>15 km) CPDs occur over outcropping granitic and Precambrian lithologies in the Panamint Range, Grapevine Mountains, Black Mountains and the Argus Range. The shallowest CPD occurs within the central Death Valley over a possible seismically imaged magma body and slightly deeper values occur within the Panamint Valley, southern Death Valley and Tecopa Hot Springs. The shallow CPD values suggest that partially molten material may also be found in these latter regions. The CPD computed heat flow values for the region suggest that the entire area has high heat flow values (>100 mW m^?2), on the other hand, locally extremely high values (>200 mW m^?2) occur within the Panamint Valley, the southern and central Death Valley and Tecopa Hot Springs region. These locally high heat flow values may be related to midcrustal magma bodies; but additional geophysical experiments are needed to determine if the magma bodies exist.     

Fast wavelet-based stochastic simulation using training images

Spatial uncertainty modelling is a complex and challenging job for orebody modelling in mining, reservoir characterization in petroleum, and contamination modelling in air and water. Stochastic simulation algorithms are popular methods for such modelling. In this paper, discrete wavelet transformation (DWT)-based multiple point simulation algorithm for continuous variable is proposed that handles multi-scale spatial characteristics in datasets and training images. The DWT of a training image provides multi-scale high-frequency wavelet images and one low-frequency scaling image at the coarsest scale. The simulation of the proposed approach is performed on the frequency (wavelet) domain where the scaling image and wavelet images across the scale are simulated jointly. The inverse DWT reconstructs simulated realizations of an attribute of interest in the space domain. An automatic scale-selection algorithm using dominant mode difference is applied for the selection of the optimal scale of wavelet decomposition. The proposed algorithm reduces the computational time required for simulating large domain as compared to spatial domain multi-point simulation algorithm. The algorithm is tested with an exhaustive dataset using conditional and unconditional simulation in two- and three-dimensional fluvial reservoir and mining blasted rock data. The realizations generated by the proposed algorithm perform well and reproduce the statistics of the training image. The study conducted comparing the spatial domain filtersim multiple-point simulation algorithm suggests that the proposed algorithm generates equally good realizations at lower computational cost.     

Shear Wave Velocity Structure of the Sinai Peninsula from Rayleigh Wave Analysis

The lithospheric structure of the Sinai Peninsula is shown by means of nine shear velocity profiles for depths ranging from zero to 50 km, determined from the Rayleigh wave analysis. The traces of 30 earthquakes, which occurred from 1992 to 1999 in and around the study area, have been used to obtain Rayleigh wave dispersion. These earthquakes were registered by a broadband station located in Egypt (KEG station). The dispersion curves were obtained for periods between 3 and 40 s, by digital filtering with a combination of MFT and TVF filtering techniques. After that, all seismic events were grouped in source zones to obtain a dispersion curve for each source-station path. These dispersion curves were inverted according to generalized inversion theory, to obtain shear wave velocity models for each source-station path, which is the main goal of this study. The shear velocity structure obtained for the Sinai Peninsula is shown through the shear velocity distributions with depth. These results agree well with the geology and other geophysical results, previously obtained from seismic and gravity data. The obtained velocity models suggest the existence of lateral and vertical heterogeneity. The shear velocity increases generally with depth for all paths analyzed in the study area. Nevertheless, in some paths a small low velocity channel in the upper or lower crust occurs. Along these profiles, it is found that the crustal structure of the Sinai Peninsula consists of three principal layers: upper crust with a sedimentary layer and lower crust. The upper crust has a sedimentary cover of 2 km thick with an average S-velocity of 2.53 km/s. This upper crust has a variable thickness ranging from 12 to 18 km, with S-wave velocity ranging from 3.24 to 3.69 km/s. The Moho discontinuity is located at a depth of 30 km, which is reflected by a sharp increase in the S-velocity values that jump from 3.70–4.12 to 4.33–4.61 km/s.     

Another look at the 1993 and 1995 Gulf of Aqaba earthquakes from the analysis of teleseismic waveforms

In this work we estimated the source process and the source parameters of the 1993 and 1995 Gulf of Aqaba earthquakes. To investigate the source process of both events, we inverted the teleseismic P- and SH-waveform data using an iterative technique in which the rupture is modeled as a series of point source subevents with varying mechanisms. The main source process of the 1993 event can be explained by two subevents with essentially the same mechanism. These two subevents had different focal depths. The second subevent was about 5 s later than the first one and about 70% of the moment rate released with the second one. The total rupture duration time was about 12 s. Our solution for the 1995 Gulf of Aqaba earthquake indicated that the event consists of three subevents with various fault geometries and about 60% of the moment released by the second subevent. The first subevent occurred at a depth of 13 km, and was followed after 9 s by the second one at a depth of 11 km; the third subevent, at a depth of 9 km occurred 19 s after the second one. The total duration of the earthquake rupture process was about 18.7 s. The obtained mechanisms for the 1993 and 1995 Gulf of Aqaba earthquakes are well correlated with the structural setting of the Gulf of Aqaba.     

Effects of diagenesis on paleoclimate reconstructions from modern and young fossil corals

The integrity of coral-based reconstructions of past climate variability depends on a comprehensive knowledge of the effects of post-depositional alteration on coral skeletal geochemistry. Here we combine millimeter-scale and micro-scale coral Sr/Ca data, scanning electron microscopy (SEM) images, and X-ray diffraction with previously published δ¹⁸O records to investigate the effects of submarine and subaerial diagenesis on paleoclimate reconstructions in modern and young sub-fossil corals from the central tropical Pacific. In a 40-year-old modern coral, we find secondary aragonite is associated with relatively high coral δ¹⁸O and Sr/Ca, equivalent to sea-surface temperature (SST) artifacts as large as −3 and −5 °C, respectively. Secondary aragonite observed in a 350-year-old fossil coral is associated with relatively high δ¹⁸O and Sr/Ca, resulting in apparent paleo-SST offsets of up to −2 and −4 °C, respectively. Secondary Ion Mass Spectrometry (SIMS) analyses of secondary aragonite yield Sr/Ca ratios ranging from 10.78 to 12.39 mmol/mol, significantly higher compared to 9.15 ± 0.37 mmol/mol measured in more pristine sections of the same fossil coral. Widespread dissolution and secondary calcite observed in a 750-year-old fossil coral is associated with relatively low δ¹⁸O and Sr/Ca. SIMS Sr/Ca measurements of the secondary calcite (1.96-9.74 mmol/mol) are significantly lower and more variable than Sr/Ca values from more pristine portions of the same fossil coral (8.22 ± 0.13 mmol/mol). Our results indicate that while diagenesis has a much larger impact on Sr/Ca-based paleoclimate reconstructions than δ¹⁸O-based reconstructions at our site, SIMS analyses of relatively pristine skeletal elements in an altered coral may provide robust estimates of Sr/Ca which can be used to derive paleo-SSTs.     

Source mechanisms of the June 2004 Tabuk earthquake sequence, Eastern Red Sea margin, Kingdom of Saudi Arabia

A sequence of earthquakes took place in June 2004 approximately 60 km southeast of Tabuk, Saudi Arabia. The first felt event (M _W = 3.9) occurred on June 9 and caused minor damage in the epicentral area according to the National Earthquake Information Center and the local reports. Another moderate size event occurred on June 22 (M _W = 5.1) and was followed by a few felt aftershocks without any reported damage. This earthquake sequence caused considerable alarm at Tabuk and highlights the fact that damaging earthquakes can occur in this region away from the major plate boundary in the Red Sea. Being the largest well-recorded event in the area for which the digital and broadband records from Saudi Arabia, Egypt, Israel, Jordan, Turkey, Cyprus, and Kuwait are available, it provides an excellent opportunity to study the tectonic process and present day stress field acting on this area. The digital records from these regional networks were used to relocate the largest three events of this sequence. Focal mechanisms were obtained from full waveform inversion and indicate normal faulting mechanisms with two nodal planes oriented NW–SE in parallel to the faults bounding the Tabuk graben and the Red Sea rift axis. These events originated at shallow focal depths of 4–5 km, possibly contributing to the widely felt ground motions. These events offer a unique opportunity to study the active tectonics of the region as well as inform future studies of seismic hazard in northwestern Saudi Arabia, the Gulf of Aqaba, and eastern Egypt.     

Imaging fracture distributions of the Al-Khuff Formation outcrops using GPR and ERT geophysical techniques,Al-Qassim area,Saudi Arabia

Khuff Formation is of utmost importance in Saudi Arabia for oil and gas reservoir although it is composed mainly of limestone. This reason refers to the existence of intensive fractures which play a vital role in the increase in porosity and permeability of this formation. The fracture pattern in the study area was verified through 2D and 3D ground penetrating radar (GPR)-defined and electrical resistivity tomography (ERT)-defined surveys. In this respect, ten of 2D GPR surveys were collected along an intersected grid of profiles covering the study area while ERT data were collected along three profiles of the GPR grid. The results were interpreted in light of the field-based structural and stratigraphic assessment of the outcropping rocks. The analysis of the inverted ERT and filtered GPR sections revealed the presence of fractures. Several resistivity and electromagnetic reflection anomalies were laterally and vertically identified across the measured sections clarifying fractures that extend to a depth of 24 m in the limestone. Most fractures are oriented vertical to sub-vertical dipping both east-west and north-south.     

Seismic hazard assessment and rheological implications: a case study selected for cities of Saudi Arabia along the eastern coast of Red Sea

A probabilistic approach is used to evaluate the seismic hazard for 12 strategic cities in Saudi Arabia along the eastern coast of Red Sea. The focal depth variations controlled by rheological characteristics are taken into account for hazard calculations, and its creditability is tested through sensitivity analysis for hazard results. This study presents a neo-probabilistic seismic hazard assessment methodology in which the focal depth distribution of earthquakes within seismogenic layer is divided into three depth slices. These depth slices are based upon rheological characteristic of seismogenic layer. The hazard results are obtained using this depth-slice methodology and conventional approach in which uniform distribution of seismicity within seismogenic layer is assumed. The sensitivity analysis culminated in underestimation of hazard values in higher frequencies for uniform distribution of seismicity within seismogenic layer. Foregoing the observations recorded above, it can be concluded that the exploitation of depth-slices biased by the rheology to calculate hazard is relatively preferable in the situations demanding safety measures.     

On the relationship between climatic variables and pressure systems over Saudi Arabia in the winter season

The synoptic circulation over Saudi Arabia is complicated and frequently governed by the effect of large-scale pressure systems. In this work, we used NCEP–NCAR global data to illustrate the relationship between climatic variables and the main pressure systems that affect the weather and climate of Saudi Arabia, and also to investigate the influence of these pressure systems on surface air temperature(SAT) and rainfall over the region in the winter season. It was found that there are two primary patterns of pressure that influence the weather and climate of Saudi Arabia. The first occurs in cases of a strengthening Subtropical High(Sub H), a weakening Siberian High(Sib H), a deepening of the Icelandic Low(Ice L), or a weakening of the Sudanese Low(Sud L). During this pattern, the Sub H combines with the Sib H and an obvious increase of sea level pressure(SLP) occurs over southern European, the Mediterranean, North Africa, and the Middle East. This belt of high pressure prevents interaction between midlatitude and extratropical systems, which leads to a decrease in the SAT,relative humidity(RH) and rainfall over Saudi Arabia. The second pattern occurs in association with a weakening of the Sub H, a strengthening of the Sib H, a weakening of the Ice L, or a deepening of the Sud L. The pattern arising in this case leads to an interaction between two different air masses: the first(cold moist) air mass is associated with the Mediterranean depression travelling from west to east, while the second(warm moist) air mass is associated with the northward oscillation of the Sud L and its inverted V-shape trough. The interaction between these two air masses increases the SAT, RH and the probability of rainfall over Saudi Arabia, especially over the northwest and northeast regions.