Geology topics

Research about Kingdom of Saudi Arabia

Source-linked reports with geographic coverage including Kingdom of Saudi Arabia.

At least 19 recordsLinked to original sources

USGS2023

Active volcanism on the Arabian Shield—Geology, volcanology, and geophysics of northern Harrat Rahat and vicinity, Kingdom of Saudi Arabia

Active volcanic systems pose serious hazards to people and property including inundation and incineration by lava, blanketing by tephra (volcanic ash), exposure to noxious volcanic gases, and damage from shallow earthquakes triggered by ascending molten material (magma). To improve understanding of volcanism and associated seismicity on the western Arabia Plate, the Saudi Geological Survey and the U.S. Geological Survey conducted a multi-year investigation of the northern Harrat Rahat volcanic field adjacent to the city of Al Madīnah al Munawwarah, Kingdom of Saudi Arabia. Project components included creation of a high-resolution digital topographic base; interpretation of eruptive history supported by detailed geologic mapping, paleomagnetism, and abundant high-precision geochronology of volcanic deposits; assessments of eruptive styles and volcanic hazards by physical volcanology; investigation of the origins of magmas in the mantle and of their differentiation in the crust revealed by chemical and isotopic petrology; gravity and magnetotelluric surveys to reveal crustal structures and to search for magma reservoirs; and regional and local seismic tomography and analyses of seismic hazards. Project results are presented in this Professional Paper as chapters written for technical scientific audiences.

Professional Paper↗

USGS2023

The Saudi Geological Survey-U.S. Geological Survey northern Harrat Rahat project—Styles, rates, causes, and hazards of volcanism near Al Madīnah al Munawwarah, Kingdom of Saudi Arabia

Al Madīnah al Munawwarah↗

USGS2023

Cenozoic tectonics of the western Arabia Plate related to harrat magmatism near Al Madīnah, Kingdom of Saudi Arabia

Sprawling volcanic fields, or harrats, in western Saudi Arabia have been emplaced during the past 30 million years following effusions of flood basalts in Ethiopia and Yemen. Although broadly associated with volcanism in three rift valleys (Red Sea, Gulf of Aden, and East African Rift Zone) radiating from the Afar depression, harrat abundance on the Arabian Peninsula indicates that volcanism is distinct from processes that opened the Red Sea and drive its magmatism. Harrats primarily lie unconformably on the Arabian Shield and locally on Paleozoic platform deposits, but several are conformable upon or interbedded with Paleogene supracrustal strata. Harrats erupted from linear chains of vents; some are oriented parallel to the Red Sea (N. 30° W.), but most trend northerly or north-northeasterly. Harrats consist mainly of weakly alkalic basalts and lesser amounts of their differentiation products, but some contain strongly alkalic rocks such as basanites and phonolites. Early harrats predate initial opening of the Red Sea 24–20 million years ago (Ma) and uplift of its flank 20–11 Ma, but younger harrats postdate the uplift and are larger and more abundant. Seismic studies identify substantially thinned continental lithosphere beneath the most voluminous harrats, but little crustal thinning is apparent at the surface, an indication that the lithosphere was eroded from below. The combination of plume impingement, rift valley formation, and erosion of the lithosphere beneath the Arabian Shield since 30 Ma indicates that the Afar plume thermally weakened the region, enabling the Arabia Plate to separate from the Africa Plate along the Red Sea and Gulf of Aden. The plume also thinned interior parts of the newly formed Arabia Plate, not just its margins, promoting mantle melting and volcanism.

Professional Paper↗

USGS2023

Eruptive history within the vicinity of Al Madīnah in northern Harrat Rahat, Kingdom of Saudi Arabia

The northernmost part of the Harrat Rahat volcanic field contains early Pleistocene to Holocene mafic eruptive products within the vicinity of the city of Al Madīnah, Kingdom of Saudi Arabia. A detailed geologic investigation into the eruptive history of a 570 square kilometer (km 2 ) area covering Al Madīnah and the surrounding area has yielded 33 mapped Quaternary volcanic units consisting of lava flows, scoria cones, and shield volcanoes. These eruptive products consist of continental, intraplate alkalic and minor transitional basalts, hawaiites, and a single mugearite that were emplaced from at least 1,014±14 thousand years ago (ka) to a single Holocene eruption in 1256 C.E. Lava flows are generally 10 to 15 kilometers (km) long (but can reach 23 km long), 1 to 3 km wide, and at least 10 meters thick. Most of the mapped units erupted episodically between 400 and 340 ka and 180 and 100 ka. Despite small individual volumes (less than 1 cubic kilometers dense rock equivalent), each unit represents eruption of a distinct magma batch that was strongly influenced by clinopyroxene, olivine, and plagioclase fractionation. Some of these units are interpreted to have undergone magma mixing pre- and (or) syneruptively. Integrating eruption ages, geochemistry, and paleomagnetic data yields evidence that some eruptions were temporally and (or) spatially clustered. Aligned scoria cones and elongate vent edifices were constructed atop fissure vent systems that reflect the local stress field, which controls dike ascent through the middle and upper crust.

Professional Paper↗

USGS2023

Eruptive history of northern Harrat Rahat—Volume, timing, and composition of volcanism over the past 1.2 million years

Harrat Rahat, one of several large, basalt-dominated volcanic fields in the western part of the Kingdom of Saudi Arabia, is a prime example of continental, intraplate volcanism. Excellent exposure makes this an outstanding site to investigate changing volcanic flux and composition through time. We present 93 40 Ar/ 39 Ar ages and 6 36 Cl surface-exposure ages for volcanic deposits throughout northern Harrat Rahat that, integrated with a new geologic map, define 12 eruptive stages. Exposed volcanic deposits in the study area erupted less than 1.2 million years ago (Ma), and 214 of 234 identified eruptions occurred less than 570 thousand years ago (ka). Two eruptions were in the Holocene, including a historically described basaltic eruption in 1256 C.E. and a trachyte eruption newly recognized as Holocene (4.2±5.2 ka). An estimated approximately 82 cubic kilometers (km 3 ; dense rock equivalent) of volcanic products can be documented as having erupted since 1.2 Ma, though this is a lower limit because of concealment of deposits older than 570 ka. Over the last 570 thousand years (k.y.), the average eruption rate was 0.14 cubic kilometers per thousand years (km 3 /k.y.), but volcanism was episodic with periods alternating between low (0.04–0.06 km 3 /k.y.) and high (0.1–0.3 km 3 /k.y.) effusion rates. Before 180 ka, eruptions vented from the volcanic field’s dominant eastern vent axis and from a subsidiary, diffuse, western vent axis. After 180 ka, volcanism focused along the eastern vent axis, and the composition of volcanism varied systematically along its length from basalt dominated in the north to trachyte dominated in the south. We hypothesize that these compositional variations younger than 180 k.y. reflect the growth of a mafic intrusive complex beneath the southern part of the vent axis, which led to the development of evolved magmas. Lastly, these new age data allow for a reassessment of the volcanic recurrence interval at northern Harrat Rahat. Based on available data, volcanism in northern Harrat Rahat over the last 180 k.y. is poorly described using a Poisson distribution with a single recurrence interval. Instead, data for northern Harrat Rahat are better described using a mixed exponential distribution that is applicable for volcanic systems characterized by two different eruptive states, where one state with a longer recurrence interval corresponding to periods of low eruption frequency and one state with a shorter recurrence interval corresponding to periods of high eruption frequency. The preferred model for northern Harrat Rahat over the last 180 k.y. uses a long recurrence interval of 4.0 k.y. and a short recurrence interval of 0.22 k.y.

Professional Paper↗

USGS2023

Lava flow emplacement in Harrat Rahat with implications for eruptions in mafic volcanic fields

Mafic volcanic fields are widespread, but few have erupted in historical times, providing limited observations of the magnitudes, dynamics, and timescales of lava flow emplacement in these settings. The Harrat Rahat volcanic field in western Saudi Arabia offers a good opportunity to study eruptions in such a setting, with a historical eruption in 1256 C.E. (654 in the year of the Hijra) and numerous well-preserved late Pleistocene lava flows. We combine historical observations and rheological and morphological analyses of the youngest flows with analytical models to reconstruct eruptive histories and lava flow emplacement conditions in Harrat Rahat. Petrologic analysis of samples for emplacement temperatures and crystallinities show cooling trends from vent to toe of ~1,140 to ~1,090 degrees Celsius (°C) at rates of 2 to 7 °C per kilometer, crystallinities increasing from 0.5 to 60 volume percent, and apparent viscosities increasing from 10 2 to 10 9 pascal seconds. High-resolution topographic data facilitate quantitative analysis of morphology and interpolation of preeruptive surfaces to measure flow thicknesses, channels, and levees, and enable calculation of eruptive volumes. Analytical models relating flow morphology to emplacement conditions are applied to estimate effusion rates. Within the suite of studied flows, minimum volume estimates range from 0.07 to 0.42 cubic kilometers dense rock equivalent, with effusion rates on the order of tens to hundreds of cubic meters per second and durations from 1 to 15 weeks. These integrated analyses quantify past lava flow emplacement conditions and dynamics in Harrat Rahat, improving our understanding and observations of fundamental parameters and controls of effusive eruptions in Harrat Rahat and other mafic volcanic fields.

northeastern Harrat Rahat↗

USGS2023

The duration and characteristics of magmatic differentiation from basalt to trachyte within the Matan volcanic center, northern Harrat Rahat, Kingdom of Saudi Arabia

A fundamental goal of igneous petrology is to quantify the durations required to produce evolved magmas following influx of basalt into the crust. However, in many cases, complex field relations and (or) the presence of a long-lived magmatic system make it difficult to assess how basaltic inputs relate to more evolved magmas, therefore precluding calculation of meaningful timescales. Here we present field relations, geochemistry, 40 Ar/ 39 Ar ages, and 36 Cl surface-exposure ages for volcanic rocks from the Matan volcanic center, located in the northern part of the Harrat Rahat volcanic field, in the Kingdom of Saudi Arabia. These data document a systematic and repeated temporal progression from alkali basalt to trachyte for the youngest eruptive products. From ~155–17 thousand years ago, the following eruptive sequence occurred four times: (1) alkali basalt, (2) hawaiite, mugearite, or benmoreite, and (3) trachyte. We interpret each eruptive sequence to result from injection of basalt into the crust, and its subsequent differentiation and eruption of progressively evolved magmas. We use the interval time between successive eruptions within a given sequence to calculate the timespans required to produce trachyte from alkali basalt. Differentiation from alkali basalt to intermediate compositions (hawaiite, mugearite, and benmoreite) took ≤3 thousand years (k.y.) on average. Differentiation from intermediate compositions to trachyte took a maximum of 6.7±3.6 to 22.9±1.7 k.y. Thus, the total duration of differentiation was ~10–25 k.y. Timescales presented here are independent of the processes evoked to drive differentiation because they are based solely on the ages and compositions of eruptive products from a system characterized by a simple, repeated differentiation sequence.

northern Harrat Rahat↗

USGS2023

Explosive trachyte eruptions from the Al Efairia volcanic center in northern Harrat Rahat, Kingdom of Saudi Arabia

Harrat Rahat is an alkali basalt, continental, intraplate volcanic field located within the central-western part of the Kingdom of Saudi Arabia. The northern quarter of Harrat Rahat contains evolved volcanic products that achieve trachyte compositions (>60 weight percent SiO 2 ). Within the Al Efairia volcanic center, pyroclastic-flow and -surge deposits that reflect explosive trachyte volcanism (and minor exposed lava domes that reflect effusive trachyte eruptions) sit at the surface as the youngest expression of volcanic activity within this part of the Harrat Rahat volcanic field. Five trachyte deposits emplaced explosively have been identified within the Al Efairia volcanic center based on geologic mapping, petrography, geochemistry, and paleomagnetism. These units are the trachytes of Um Rgaibah, Gura 5, Gura 4, Al Efairia, and Al Qayf, in descending stratigraphic order. Here, we present 14 40 Ar/ 39 Ar analyses from four of these units, which yield eruption ages of 4.2±5.2 thousand years (ka) for the trachyte of Um Rgaibah, 79.7±1.6 ka for the trachyte of Gura 5, 84.3±1.6 ka for the trachyte of Gura 4, and 88.0±1.8 ka for the trachyte of Al Efairia. The eruption age of the trachyte of Al Qayf has been constrained to between 410.3±3.4 and 418.8±1.9 ka using paleomagnetic correlations and 40 Ar/ 39 Ar ages from overlying and underlying intermediate composition lava flows. Most of these trachytes have distinct geochemical compositions, petrographic characteristics, and directions of remanent magnetization. The exceptions are for the trachytes of Gura 4 and Gura 5, which overlap in their geochemical, petrographic, paleomagnetic, and geochronologic affinities. Based on these similarities, we interpret the trachytes of Gura 4 and Gura 5 to have erupted during a closely spaced (a few decades) time interval from the same magma batch but from craters that are >2 kilometers (km) apart. The eruption of the trachyte of Al Efairia at 88.0±1.8 ka is the result of a different magma batch that erupted a few thousand years prior to the trachytes of Gura 4 and Gura 5. The Al Efairia volcanic center is remarkably different from the Matan volcanic center located ~10 km to the north, which has also erupted young (<150 ka) trachytes. The Matan volcanic center has been shown to produce trachyte compositions only after eruption of basalt followed by intermediate lava flows, whereas only trachyte compositions have erupted within the Al Efairia volcanic center over this same time interval.

Al Efairia volcanic center, Harrat Rahat↗

USGS2023

Paleomagnetism of the Harrat Rahat volcanic field, Kingdom of Saudi Arabia—Geologic unit correlations and geomagnetic cryptochron identifications

Paleomagnetic rock samples were collected from 173 drill sites in the Quaternary alkali basaltic volcanic field of northern Harrat Rahat, Kingdom of Saudi Arabia. Laboratory measurements on these samples established that lava flows and vent complexes—identified and mapped from field characteristics, rock types, and compositions as products of single or temporally close eruptions—typically record single, or very similar, directions of remanent magnetization. Correlations defined through geologic mapping, spatial association, geochemistry, geochronology, and identical mean remanent directions indicate at least 16 brief episodes of temporally clustered eruptions. These episodes had durations of a few centuries or less. Anomalous remanent magnetic directions were found for at least 13 mapped lavas of northern Harrat Rahat, which demonstrate that they were acquired during brief geomagnetic cryptochrons during the Brunhes Normal Polarity Chron. These uncharacteristic directions enhance the opportunity to identify common eruptive episodes, and to better understand and evaluate assessments of eruption ages based on 40 Ar/ 39 Ar geochronology. Combining paleomagnetic and regional archaeomagnetic results for the youngest eruptions allows us to evaluate their historical age assignments and, in one case, refute a previously assigned provisional age.

Harrat Rahat volcanic field↗

USGS2023

Mantle origin and crustal differentiation of basalts and hawaiites of northern Harrat Rahat, Kingdom of Saudi Arabia

Quaternary volcanic rocks of northern Harrat Rahat, Kingdom of Saudi Arabia, are chiefly alkali basalts with subordinate transitional basalts, hawaiites, mugearites, benmoreites, and trachytes. Geochemical and isotopic results indicate that crystallization-differentiation, mixing, and cumulate reassimilation within the magmatic system produced most of its compositional diversity, with only minor involvement of Neoproterozoic crust. With increasing evolution, crystal assemblages that separated from and drove basalt-to-hawaiite differentiation passed from (1) dunitic or troctolitic to (2) olivine gabbroic to (3) titanomagnetite-bearing olivine gabbroic, with typical hawaiites representing about 20 weight percent residual liquids from an estimated primary alkali basaltic parent. Crystallization-differentiation simulations for midcrustal pressures yield the closest compositional match to the basalt-hawaiite suite, and stagnation in the midcrustal area (near 20 kilometers [km] depth) may result from density trapping beneath the lower density upper continental crust. Least differentiated alkali basalts have fractionation-adjusted major-oxide compositions that are consistent with separation from the deeper parts (60–70 km) of the spinel-lherzolite stability field at pressures that are close to the local lithosphere-asthenosphere boundary (60–80 km). Mantle potential temperature estimates are strongly sensitive to modeling approach, but potential temperatures (1,345–1,390 degrees Celsius [°C]) are not discernably greater than for midocean ridge basalts (MORB; 1,350–1,410 °C) if adherence to spinel-lherzolite melting relations is required. Inversion of the trace-element concentrations of the lesser fractionated basalts indicates a depleted mantle source, similar to MORB-source estimates, but one that is enriched in Sr and includes a greater relative proportion of melting in the garnet-peridotite stability field. These geochemical and thermal relations, as well as radiogenic isotopes, point to a dominantly depleted mantle asthenospheric source for Harrat Rahat basalts, admixed with subordinate materials either from the Afar mantle plume or an enriched MORB component in the ambient asthenosphere. The lithosphere-asthenosphere boundary is shallower beneath the belt of major volcanic fields on the Arabia Plate, and restoration of rifting across the Red Sea and Gulf of Aden places the south end of this belt adjacent to the northern part of the Afar region, suggesting a once-continuous structure, possibly an arch, a weakness, or a discontinuity along the base of the lithosphere. Magma generation can be ascribed to focused upwelling and decompression melting, perhaps driven by a magmatic-feedback mechanism whereby basaltic intrusions into the deep lithosphere solidify as eclogites, causing lithospheric foundering and further asthenospheric upwelling and decompression melting in a restricted region.

northern Harrat Rahat↗

USGS2023

Isotopic and geochemical evidence for the source of volcanism at Harrat Rahat, Kingdom of Saudi Arabia

Pleistocene and Holocene basalts, hawaiites, mugearites, benmoreites, and trachytes from the northern part of the Harrat Rahat volcanic field, Kingdom of Saudi Arabia, were analyzed for Sr, Nd, Hf, and Pb isotopic compositions. Evolved trachytes with Mg number <0.1 (Mg# = Mg/[Mg+Fe 2+ ], molar) have relatively radiogenic Sr isotopic compositions indicating that they were influenced by contamination probably in the upper crust. Volcanic rocks with Mg# >0.1, consisting chiefly of alkali basalts but encompassing hawaiites, mugearites, and benmoreites, show a limited range in Hf, Nd, Sr, and Pb isotopic compositions. Although the total Pb isotope variation is only 1 percent, the Pb isotope values correlate with Mg#, where the least radiogenic Pb is in samples with the lowest Mg#. The trend formed in Pb isotope space points toward an unradiogenic Pb composition that is similar to the Pb isotopic composition of lower crust of the Precambrian Arabian-Nubian Shields, as well as to feldspars and galena in the upper crust of the western Arabian Shield. This trend is interpreted as progressive but overall minor (no more than 5 weight percent) assimilation of shield rocks, or their partial melts, during fractional crystallization. Isotopic compositions of the least evolved northern Harrat Rahat magmas are most similar among analyzed Arabian harrats to depleted spreading-ridge basalts of the active Red Sea rift, but isotopic values are displaced toward those of spreading-ridge basalts of the Gulf of Aden that are proximal to the site of the Afar mantle plume. The Pb isotopic compositions very near the Northern Hemisphere Reference Line indicate no discernable lithospheric contribution to yield the parental basalts of northern Harrat Rahat, and their isotopic compositions are consistent with derivation predominantly from depleted Northern Hemisphere asthenosphere with a subordinate (20–30 weight percent) component from the Afar mantle plume. Trace-element variations show that appreciable portions of melting were in the garnet stability field, confirming the sub-lithospheric origin of the magmas, and that melting extents were low, accounting for the alkalic, trace-element-enriched character of the suite. The presence of possible Afar mantle beneath the western part of the Arabian Shield and its absence beneath the Red Sea rift may result from capture and channelized flow along high-relief structures incised into the base of the sub-continental lithosphere, as revealed by geophysical images.

northern Harrat Rahat↗

USGS2023

Depth to basement and crustal structure of the northern Harrat Rahat volcanic field, Kingdom of Saudi Arabia, from gravity and aeromagnetic data

New gravity data reveal a prominent negative anomaly along the main vent axis of the northern Harrat Rahat volcanic field in the Kingdom of Saudi Arabia. The gravity low continues north of the volcanic field onto exposures of Proterozoic rocks, indicating that the low is caused not only by the volcanic field (and possibly underlying Cenozoic sediments), but also the underlying Proterozoic basement. An inversion of the gravity field guided by analysis of aeromagnetic data indicates (1) a broad depression of the basement surface that is deeper along the main vent axis in the eastern part and in the southwest part of the volcanic field and (2) less dense basement beneath the vent axis. Low densities within the basement most likely arise from lithologic variations in the basement, predating Cenozoic volcanism, although our analysis does not rule out small volumes of partial melt and higher temperatures or extensive fracturing at depth.

northern Harrat Rahat↗

USGS2023

Magnetotelluric investigation of northern Harrat Rahat, Kingdom of Saudi Arabia

Volcanism within the harrats (Arabic for “volcanic field”) of the Kingdom of Saudi Arabia includes at least one historical eruption occurring close to the holy city of Al Madīnah in 1256 C.E. As part of a volcanic- and seismic-hazard assessment of northern Harrat Rahat, magnetotelluric (MT) data were collected to investigate the structural setting of the area, the presence or absence of melt within the crust, and the mantle-derived magmatic source. Collected MT data were modeled in both two dimensions, where anisotropy can be estimated, and three dimensions. Interpretation of the preferred resistivity model includes a shallow sediment-filled graben beneath northern Harrat Rahat lavas, a melt-free upper crust, and a region of decompression melting in the asthenosphere below 60–70 kilometers depth. Models in two dimensions image the lower crust as anisotropic, demonstrating that a series of elongate conductivity anomalies with a strike of N. 10° E. within the lower crust of the three-dimensional model are artifacts of inverting anisotropic data with an isotropic modeling algorithm. Careful examination of the resistivity models, in combination with regional geological and geophysical data, suggests an anisotropic lower crust that is free of large zones of melt. Azimuthal anisotropy in the lower crust extends well beyond the limits of Harrat Rahat volcanic rocks, with a conductive direction oriented N. 10° E. and an anisotropy factor of 2–5 between the most and least conductive directions. Enhanced conductivity is likely caused by interconnected grain-boundary graphite, where the direction of anisotropy reflects either frozen-in fabric from the Neoproterozoic stabilization of the Arabian Shield or ductile deformation driven by channelized asthenospheric flow coupled with a thin rigid mantle lid. Asthenospheric melt is interpreted to transect the crust largely through diking, with limited melt storage and short residence times within the crustal column.

northern Harrat Rahat↗

USGS2023

Thickness of the Saudi Arabian crust

As part of a joint Saudi Geological Survey (SGS) and U.S. Geological Survey (USGS) project, we analyzed P-wave receiver functions from seismic stations covering most of the Kingdom of Saudi Arabia to map the thickness of the crust across the Arabia Plate. We present an update of crustal-thickness estimates and fill in gaps for the western Arabian Shield and the rifted margin at the Red Sea (the coastal plain), as well as the eastern Arabian Platform. We applied a conventional H-k stacking algorithm and included careful attention to stacking weights, two forms of sedimentary corrections for stations located on the Arabian Platform, and additional processing for noisy stations. We obtained useful results at 154 stations from 898 teleseismic events over a 2-year period from 1995–1997 (for non-SGS stations) and a 6-year period from 2008–2014 (for SGS stations). Average crustal thickness (that is, depth to the Mohorovičić discontinuity [Moho] below the surface) beneath the Red Sea coastal plain (the rift margin) is 29 kilometers (km), beneath the volcanic fields (known in Arabic as harra [plural] or harrat [singular]) is 35 km, beneath the Arabian Shield (excluding harrats) is 37 km, and beneath the Arabian Platform is 38 km. Crustal thinning appears not to extend east of the rift escarpment, suggesting uniform extension that is no broader at depth than at the surface. In contrast to some previous claims that the Arabian Platform crust is thicker than that of the Arabian Shield, we find no statistically significant difference between their whole crustal thicknesses. However, the average subsedimentary crustal thickness (that is, the crystalline crust) for stations on the Arabian Platform is 34 km, 3 km thinner than the crust of the Arabian Shield. Individual station P-wave (pressure) velocity and S-wave (shear) velocity ratios ( V P / V S ) are highly variable for the Arabia Plate, ranging from 1.60 to 1.97 and averaging 1.75, with a standard deviation of 0.07. There are no statistically significant differences between V P / V S ratios of the different geologic regions of Saudi Arabia. Similar V P / V S ratios, coupled with similar crustal thicknesses for harrats and the Arabian Shield, indicate that Cenozoic magmatism has contributed negligibly to crustal growth.

Professional Paper↗

USGS2023

Ambient seismic noise tomography of the Kingdom of Saudi Arabia

Harrat Rahat is a Cenozoic volcanic field in the west-central part of the Kingdom of Saudi Arabia, 150 kilometers east of the Red Sea, and is the site of the most recent eruption in the country (1256 C.E.; 654 in the year of the Hijra). The city of Al Madīnah lies at the north end of Harrat Rahat, and its volcanic and seismic risks are frequently reassessed. In 2009 C.E. an earthquake swarm at Harrat Lunayyir, 200 km west-northwest of Al Madīnah, produced significant shaking and some building damage in nearby communities, which prompted a revision of seismic hazard models for the region. A network of seismic stations on this volcanically active western side of the Arabia Plate was installed, and stations were also added in the tectonically active northern part of the country. Although regional earthquakes may be used to determine the crustal structure of the western Arabia Plate, such crustal models are hindered by insufficient numbers of earthquakes in the stable plate interior. Tomography studies can be used to infer material properties of the subsurface, such as presence of partial melt, and are beneficial for volcanic hazard assessments. We use ambient seismic noise to compute Rayleigh and Love surface-wave dispersion maps between 5 and 12 second periods for a subset of seismic stations within and near northern Harrat Rahat. The surface-wave maps are inverted to produce shear-wave velocities using a neighborhood algorithm and interpolated into a pseudo three-dimensional model. The distributions of surface-wave and shear-wave velocities are heterogenous, varying from ±3–8 percent. However, low velocities are not restricted to Harrat Rahat. We observed a difference between Rayleigh- and Love-wave velocities that extends north of the site of the 1256 C.E. eruption and coincides with a low gravity anomaly. We obtain a shear-wave velocity increase of 10–15 percent between 15- and 25-kilometers depth, which is consistent with the presence of a transition between the felsic upper crust and the mafic lower crust of the Arabian Shield. The average shear-wave velocities of the upper and lower crust are estimated to be 3.64 and 3.95 kilometers per second using Rayleigh waves and 3.53 and 4.16 kilometers per second using Love waves, which are in good agreement with the results of other geophysical surveys in this area. The modest magnitude of the low-velocity anomalies within the crust and their locations extending well beyond the limits of Harrat Rahat indicate that they are not caused by a crustal magma chamber. If magma chambers exist, they are smaller than can be imaged with our seismological method (resolution corresponds to a 15-kilometer wavelength), deeper than 30 kilometers, or shallower than 5 kilometers with a small velocity contrast. The last possibility is discounted by weak-to-absent surface geothermal activity and by a lack of shallow seismicity that characterizes other areas with known shallow magmas, such as Hawaiʻi. We then expanded our work to the entire Arabian Shield with principally the same methodology to look at large-scale regional patterns. We found modest shear-wave velocity deviations on the order of only ±3 percent, which are within expected ranges for lithological variation.

Professional Paper↗

USGS2023

Ground-motion prediction equations for the western Kingdom of Saudi Arabia

Ground-motion prediction equations (GMPEs) for the western Kingdom of Saudi Arabia are developed by employing a mixed-effects regression model to modify the Boore and others (2014) Next Generation Attenuation-West2 (NGA-West2) project GMPEs. NGA-West2 addressed several key issues concerning GMPEs for shallow crustal earthquakes in active tectonic regions. However, the NGA-West2 input data do not include many earthquakes in extensional regimes, such as those occurring in Saudi Arabia. This deficiency is redressed by calculating a magnitude scaling of the new Saudi Arabia-specific GMPEs compared to those of Boore and others (2014). Furthermore, there is a clear difference in distance scaling for the Saudi Arabian GMPEs in comparison with the NGA-West2 GMPEs. This difference is especially significant at large distances and is mainly due to lower anelastic attenuation in the crystalline crust of western Saudi Arabia. Our empirical data demonstrate that the GMPEs presented here are in good agreement with observed earthquake ground motions in western Saudi Arabia.

Professional Paper↗

USGS2023

Seismic hazard assessment for areas of volcanic activity in western Kingdom of Saudi Arabia

Earthquake swarms caused by volcanic activity, tectonic stresses, or industrial operations (oil and gas production) can pose considerable risk for nearby settlements. As a rule, a probabilistic seismic hazard assessment (PSHA) that is based on time-independent earthquakes does not take into account earthquake swarms because of their statistically time-dependent nature. We describe the technique and application of a scenario-based method for the estimation of seismic hazard from earthquake swarms that may occur in areas of volcanic activity in western Saudi Arabia. The method consists of the generation of a large number of scenario seismic swarms followed by the calculation of ground motion in a site of interest for all earthquakes in all of the swarms. The set of calculated ground-motion values permits the construction and analysis of probability distribution functions for possible ground-motion levels during the hypothetical earthquake swarms. The level of ground motion representing the hazard is selected considering an appropriate value of the probability of exceedance or non-exceedance, depending on the goal of the study. This swarm-scenario-based seismic hazard assessment may provide a valuable supplement to an ordinary PSHA and (or) deterministic seismic hazard assessment that are commonly used for emergency response and post-earthquake recovery management.

Professional Paper↗

USGS2023

Probabilistic seismic-hazard analysis for the western Kingdom of Saudi Arabia

We present a probabilistic seismic-hazard analysis (PSHA) for the west-central part of the Arabian Peninsula. Our study area includes the northern Harrat Rahat volcanic field and the nearby city of Al Madīnah, Kingdom of Saudi Arabia. This young, active volcanic field experienced one historical eruption in 1256 C.E. (654 in the year of the Hijra) that vented 20 to 22 kilometers (km) southeast of the center of Al Madīnah, which has a present population of about 1.4 million. The field also erupted numerous times in the late Pleistocene and possibly in the early Holocene. We used recently developed regional ground-motion prediction equations for Saudi Arabia to calculate the severity of ground shaking as a function of distance from the earthquake source. This information was combined with two key volcanic parameters for this area: (1) the return period of volcanic eruptions and (2) the spatial probability of the next vent opening. The calculated ground-motion levels of peak ground acceleration and peak ground velocity with a 2-percent probability of exceedance in 50 years are expected to be about 0.14 the acceleration due to gravity ( g ) and 10 centimeters per second (cm/s), respectively, at the most probable vent opening location, about 25 km southeast of Al Madīnah’s center, and about 0.07 g and 3 cm/s, respectively, in the city interior. These ground motions are higher than previous estimates that did not consider the nearby earthquakes associated with volcanic activity.

Professional Paper↗