Search USGSSearch

Geology topics

M. E. Gettings

Publications and source records attributed to M. E. Gettings.

34 records · Page 2Linked to original sources

Saudi Arabian seismic deep-refraction profiles; final project report

In February 1978 a seismic deep-refraction profile was recorded by the U.S. Geological Survey along a 1000-km line across the Arabian Shield in western Saudi Arabia. The line begins in Mesozoic cover rocks near Riyadh on the Arabian Platform, leads southwesterly across three major Precambrian tectonic provinces, traverses Cenozoic rocks of the coastal plain near Jizan (Tihamat-Asir), and terminates at the outer edge of the Farasan Bank in the southern Red Sea. More than 500 surveyed recording sites were occupied, including 19 in the Farasan Islands. Six shot points were used: five on land, with most charges placed below the water table in drill holes, and one at sea, with charges placed on the sea floor and detonated from a ship. Slightly more than 61 metric tons of explosives were used in 19 discrete firings. Seismic energy was recorded by 100 newly-developed portable seismic stations deployed in approximately 200 km-long arrays for each firing. Each station consisted of a standard 2-Hz vertical component geophone coupled to a self-contained analog recording instrument equipped with a magnetic-tape cassette. In this final report, we fully document the field and data-processing procedures and present the final seismogram data set as both a digital magnetic tape and as record sections for each shot point. Record sections include a normalized set of seismograms, reduced at 6 km/s, and a true-amplitude set, reduced at 8 km/s, which have been adjusted for amplifier gain, individual shot size, and distance from the shot point. Appendices give recorder station and shot information, digital data set descriptions, computer program listings, arrival times used in the interpretation, and a bibliography of reports published as a result of this project. We used two-dimensional ray-tracing techniques in the data analysis, and our interpretation is based primarily on horizontally layered models. The Arabian Shield is composed, to first-order, of two layers, each about 20 km thick, with average velocities of 6.3 km/s and 7.0 km/s, respectively. At the western shield margin the crust thins to less than 20 km total thickness, beyond which the Red Sea shelf and coastal plain are interpreted to be underlain by oceanic crust. A major crustal lateral velocity inhomogeneity northeast of Sabhah in the Shammar Tectonic Province is interpreted as the suture zone of two crustal blocks of different composition. Several high-velocity anomalies in the upper crust correlate with mapped gneissic dome structures. Two intra-crustal reflectors at13 km depth are interpreted as the tops of mafic intrusives. The Mohorovicic discontinuity beneath the shield varies from 43 km depth in the northeast with 8.2 km/s mantle velocity to 38 km depth in the southwest with 8.0 km/s mantle velocity. Two velocity discontinuities are identified in the upper mantle, at 59 and 70 km depth. We suggest further work, including refined analyses of the data employing filtering and synthetic seismogram techniques, as well as consideration of attenuation properties. Extension of the seismic refraction profile to the Arabian Gulf and some short profiles perpendicular to the existing profile would be fruitful areas for future field work.

Open-File Report

Selected physical properties of rocks from the Baid al Jimalah West tungsten deposit, Kingdom of Saudi Arabia, and recommendations for geophysical surveys

Bulk density and magnetic susceptibility of 11 outcrop samples representing the Proterozoic lithologic units at the Baid al Jimalah West tungsten deposit, Kingdom of Saudi Arabia were measured. Induced polarization response, apparent resistivity, and electromagnetic conductivity were determined for four specimens of the sample suite. Measurements show that there is a negative density contrast of about -0.17 g-cm^-3 between metasedimentary rocks of the Murdama group and the Bald al Jimalah graaite and that this contrast decreases with increasing mineralization of the granite. Similarly, the bulk magnetic susceptibility of the granite is about one-third that of the Murdama rocks for this sample suite; however, magnetic susceptibility increases with increasing mineralization in the granite specimens. Electromagnetic conductivities are uniformly low, in part because the specimens are weathered, but probably also because intense silicification accompanies the mineralization. Induced polarization chargeability increases in the granitic specimens with increasing mineralization and reflects higher percentages of sulfide minerals. Chargeability for the mineralized rocks is about four times higher than for the Murdama host rocks, and apparent resistivity values are about one-fifth the values of host rocks. Based on these results, it is recommended that during reconnaissance exploration of the area 15 detailed high-precision gravity profiles at 10 m to 50 m station spacing and eight induced polarization dipole-dipole profiles at 25 m dipole spacing and maximum 'n' of 6 be measured. To help define subsurface structure, a high-precision, ground-magnetic survey (map at 2-gamma contour interval) and a four-channel gamma ray spectrometric survey on a 25x50 m grid covering the area of the profiles are recommended.

Open-File Report

Crustal structure of southwestern Saudi Arabia

The southwestern Arabian Shield is composed of uplifted Proterozoic metamorphic and plutonic rocks. The Shield is bordered on the southwest by Cenozoic sedimentary and igneous rocks of the Red Sea paar and on the east by the Arabian Platform, an area of basin sedimentation throughout Phanerozoic time. The Shield appears to have been formed by successive episodes of island arc volcanism and sea-floor spreading, followed by several cycles of compressive tectonism and metamorphism. An interpretation and synthesis of a deep-refraction seismic profile from the Riyadh area to the Farasan Islands, and regional gravity, aeromagnetic, heat flow, and surface geologic data have yielded a self-consistent regional-scale model of the crust and upper mantle for this area. The model consists of two 20 km-thick layers of crust with an average compressional wave velocity in the upper crust of about 6.3 km/s and an average velocity in the lower. crust of about 7.0 km/s. This crust thins abruptly to less than 20 km near the southwestern end of the profile where Precambrian outcrops abut the Cenozoic rocks and to 8 km beneath the Farasan Islands. The data over the coastal plain and Red Sea shelf areas are fit satisfactorily by an oceanic crustal model. A major lateral velocity inhomogeneity in the crust is inferred about 25 km northeast of Sabhah and is supported by surface geologic evidence. The major velocity discontinuities occur at about the same depth across the entire Shield and are interpreted to indicate horizontal metamorphic stratification of the Precambrian crust. Several lateral inhomogenities in both the upper and lower .crust of the . Shield are interpreted, to indicate bulk compositional variations. The subcrustal portion of the model is composed of a hot, low-density lithosphere beneath the Red Sea which is systematically cooler and denser to the northeast. This model provides a mechanism which explains the observed topographic uplift, regional gravity pattern, heat flow, and mantle compressional wave velocities. Such a lithosphere could be produced by upwelling of hot asthenosphere beneath the Red Sea which then flows laterally beneath the lithosphere of the Arabian Plate.

Open-File Report

A simple Bouguer gravity anomaly map of southwestern Saudi Arabia and an initial interpretation

Approximately 2,200 gravity stations on a 10-km2 grid were used to construct a simple Bouguer gravity anomaly map at 1:2,000,000 scale along a 150-km-wide by 850-km-long strip of the Arabian Peninsula from Sanam, southwest of Ar Riyad, through the Farasan Islands and including offshore islands, the coastal plain, and the Hijaz-Asir escarpment from Jiddah to the Yemen border. On the Precambrian Arabian Shield, local positive gravity anomalies are associated with greenstone belts, gneiss domes, and the Najd fault zones. Local negative gravity anomalies correlate with granitic plutonic rocks. A steep gravity gradient of as much as 4 mgal-km-1 marks the continental margin on the coastal plain near the southwestern end of the strip. Bouguer gravity anomaly values range from -10 to +40 mgal southwest of this gradient and from -170 to -100 mgal in a 300-km-wide gravity minimum northeast of the gradient. Farther northeast, the minimum is terminated by a regional gradient of about 0.1 mgal-km-1 that increases toward the Arabian Gulf. The regional gravity anomaly pattern has been modeled by using seismic refraction and Raleigh wave studies, heat-flow measurements, and isostatic considerations as constraints. The model is consistent with the hypothesis of upwelling of hot mantle material beneath the Red Sea and lateral mantle flow beneath the Arabian plate. The model yields best-fitting average crustal densities of 2.80 g-cm-3 (0-20 km depth) and 3.00 g-cm-3 (20-40 km depth) southwest of the Nabitah suture zone and 2.74 g-cm-3 (0-20 km depth) and 2.94 g-cm-3 (20-40 km depth) northeast of the suture zone. The gravity model requires that the crust be about 20 km thick at the continental margin and that the lower crust between the margin and Bishah (lat 20? N., long 42.5? E.) be somewhat denser than the lower crust to the northeast. Detailed correlations between 1:250,000- and 1:500,000-scale geologic maps and the gravity anomaly map suggest that the greenstone belts associated with gravity highs contain a large proportion of gabbroic and dioritic intrusive rocks and that the bulk density of the upper crust associated with some of the batholithic complexes has been lowered by the large-scale intrusion of granitic material at depth, as well as by that exposed at the surface. A comparison of known base and precious metals occurrences with the Bouguer gravity anomaly field shows, in some cases, a correlation between such occurrences and the features of the gravity anomaly map. Several areas were identified between known mineral occurrences along gravity-defined structures that may contain mineral deposits if the lithologic environment is favorable.

Open-File Report

Near-station terrain corrections for gravity data by a surface-integral technique

A new method of computing gravity terrain corrections by use of a digitizer and digital computer can result in substantial savings in the time and manual labor required to perform such corrections by conventional manual ring-chart techniques. The method is typically applied to estimate terrain effects for topography near the station, for example within 3 km of the station, although it has been used successfully to a radius of 15 km to estimate corrections in areas where topographic mapping is poor. Points (about 20) that define topographic maxima, minima, and changes in the slope gradient are picked on the topographic map, within the desired radius of correction about the station. Particular attention must be paid to the area immediately surrounding the station to ensure a good topographic representation. The horizontal and vertical coordinates of these points are entered into the computer, usually by means of a digitizer. The computer then fits a multiquadric surface to the input points to form an analytic representation of the surface. By means of the divergence theorem, the gravity effect of an interior closed solid can be expressed as a surface integral, and the terrain correction is calculated by numerical evaluation of the integral over the surfaces of a cylinder, The vertical sides of which are at the correction radius about the station, the flat bottom surface at the topographic minimum, and the upper surface given by the multiquadric equation. The method has been tested with favorable results against models for which an exact result is available and against manually computed field-station locations in areas of rugged topography. By increasing the number of points defining the topographic surface, any desired degree of accuracy can be obtained. The method is more objective than manual ring-chart techniques because no average compartment elevations need be estimated ?

Open-File Report

Preliminary geological investigation of the Bir Hayzan diatomite deposit, Kingdom of Saudi Arabia, with a section on selected physical properties and implications for future geophysical exploration

A 2.2-m-thick lake bed composed entirely of diatomite was found in the southwestern arm of the Nafud sand sea, about 67 km east of Tayma in northwestern Saudi Arabia. Its stratigraphic position above eolian sand and below a large barchanoid dune indicates deposition during an exceptional period of pluvial climate in an otherwise arid sand desert. Radiocarbon-dated carbonate lake beds found in similar stratigraphic settings in other parts of the Nafud suggest that the diatomite was deposited in the late Pleistocene epoch sometime between 32,000 and 20,000 years ago. Diatomite of the same pluvial episode has also been found in the Afar region of Ethiopia. Four samples of the diatomite were analyzed for major chemical components and 32 minor and trace elements. Two of the samples appeared to be representative of the lake bed. These samples contained about 85 percent SiO2, 2.34 percent Fe2O3, 1 4 percent Al2O3, and 8.27 percent ignition loss, with only traces of CaO, P2O5, MgO, Na2O, and TiO2. The high content of silica and low contents of alumina, iron oxide, lime, and other constituents indicate low contamination. Analyses of samples from the two atypical units indicate a small amount of contamination, probably clay or mica, in the basal unit and addition of iron oxide and carbonate by ground water in an isolated green-colored unit. Chemically the Bi'r Hayzan diatomite compares favorably with commercial deposits from other parts of the world. Low bulk densities of 0.64 g/cm3 (30-40 lbs/ft3) and an average porosity of 63 percent for five samples also suggest the deposit is of industrial grade. Further studies are recommended to determine: (1) the size and extent of the deposit and (2) the suitability of the diatomite for commercial applications. As means to these ends, electrical-resistivity surveys following winter rains are suggested to determine both the extent and thickness of the lake bed beneath the overlying sand dunes. A complete species identification of the diatom association and further laboratory testing of physical properties of the diatomite also should be made. Further exploration in the Nafud is recommended to locate other possible diatomite deposits.

Open-File Report

Heat-flow measurements at shot points along the 1978 Saudi Arabia seismic deep-refraction line; Part I, Results of the measurements

Heat-flow measurements were made at five onland shot points of the 1978 Saudi Arabian seismic deep-refraction line, which sample major tectonic elements of the Arabian Shield along a profile from Ar Riyad to the Farasan Islands. Because of the pattern drilling at each shot point, several holes (60 m deep) could be logged for temperature at each site and thus allow a better estimate of the geothermal gradient. Each site was mapped and sampled in detail, and modal and. chemical analyses of representative specimens were made in the laboratory. Thermal conductivities were computed from the modal analyses and single-mineral conductivity data. The resulting heat-flow values, combined with published values for the Red Sea and coastal plain, indicate a three-level pattern, with a heat flow of about 4.5 heat-flow unit (HFU) over the Red Sea axial trough, about 3.0 HFU over the shelf and coastal plain, and an essentially constant 1.0 HFU over the Arabian Shield at points well away from the suture zone with the oceanic crust. At three sites where the rocks are granitic, gamma-ray spectrometry techniques were employed to estimate thorium, potassium, and uranium concentrations. The resulting plot of heat generation versus heat flow suggests that in the Arabian Shield the relationship between heat flow and heat production is not linear. More heat-flow data are essential to establish or reject this conclusion.

Open-File Report

Heat-flow measurements at shot points along the 1978 Saudi Arabia seismic deep-refraction line; Part II, Discussion and interpretation

The heat-flow profile across the Arabian Shield from Ar Riyad to Ad Darb and across the Red Sea is examined for compatibility with the lithospheric structure of the area as deduced from geologic and other geophysical data. Broad continental uplift associated with Red Sea rifting is symmetric about the Red Sea axis, and geologic and geochronologic evidence indicate that uplift has occurred mainly in the interval 25-13 Ma (mega-annum) ago. Thermal-profile changes in the upper mantle resulting from an influx of hot material associated with rifting yield the correct order of magnitude of uplift, and this mechanism is suggested as the explanation for the regional doming. A lithospheric section, constructed from seismic refraction, gravity, and regional geologic data, provides the framework for construction of thermal models. Thermal gradient measurements were made in drill holes at five shot points. Geotherms for the Shield, which assume a radiogenic heat-source distribution that decreases exponentially with depth, yield temperatures of about 450?C at a depth of 40 km (base of the crust) for shot points 2 (Sabhah) and 3. The geotherm for shot point 4 (near Bishah) yields a distinctly higher temperature (about 580?C) for the same depth. Static models used to model the heat flow in the oceanic crust of the Red Sea shelf and coastal plain either yield too small a heat flow to match the observed heat flow or give lithosphere thicknesses that are so thin as to be improbable. Dynamic (solid-state accretion) models, which account for mantle flow at the base of the lithosphere, adequately match the observed heat-flow values. In the deep-water trough of the Red Sea, which is presently undergoing active sea-floor spreading, classical models of heat flow for a moving slab with accretion at the spreading center are adequate to explain the average heat-flow level. At shot point 5 (Ad Darb), the anomalous heat flow of 2 HFU (heat-flow units) can be explained in terms of a Shield component (0.8-1.0 HFU) and a component related to heating by the abutting oceanic crust a few kilometers away for periods exceeding 10 Ma. Analytical results are included for: 1) the cooling of a static sheet with an initial temperature distribution characteristic of a moving slab in a sea-floor spreading environment, and 2) the heating of a homogeneous quarter-space at its vertical boundary.

Open-File Report

Status of Regional Gravity Project, Kingdom of Saudi Arabia, as of 1 January, 1979

Regional gravity data have now been obtained for approximately 33 percent of the area of the Arabian Shield, associated harrat areas, and the Red Sea coastal plain. Virtually no surveys have been done on the Phanerozoic platform around the edges of the shield, and although regional gravity coverage is available for the Eastern Province from the Arabian American Oil Company, it is confidential. Approximately 4,500 more gravity stations at 10 km spacing are needed to complete the coverage of the shield, and 5,250 are needed on the Phanerozoic platform to the north and east of the shield. The initial gravity maps compiled from the regional gravity project will be along the "geophysical strip", a band approximately 150 km wide across the shield from Al Qawayiyah to the Farasan Islands. As more data become available, additional maps will be compiled at 1:500,000 scale covering the same areas as the USGS Miscellaneous Investigations Maps 1-700-A to I-220-A. Computer software development and implementation for gravity data reduction, storage, retrieval, and display are 80 percent complete. Development and implementation of software for routine semi-quantitative and quantitative interpretation are 50 percent complete. If the Directorate General of Mineral Resources were to commit full support for the regional gravity project, phase 1 (data acquisition) and phase 2 (data reduction and publication of maps and principal facts) could be completed and phase 3 (interpretation) could be about 50 percent completed by the end of the Third Five Year Development Plan (1985).

Open-File Report

Magma beneath Yellowstone National Park

The Yellowstone plateau volcanic field is less than 2 million years old, lies in a region of intense tectonic and hydrothermal activity, and probably has the potential for further volcanic activity. The youngest of three volcanic cycles in the field climaxed 600,000 years ago with a voluminous ashflow eruption and the collapse of two contiguous cauldron blocks. Doming 150,000 years ago, followed by voluminous rhyolitic extrusions as recently as 70,000 years ago, and high convective heat flow at present indicate that the latest phase of volcanism may represent a new magmatic insurgence. These observations, coupled with (i) localized postglacial arcuate faulting beyond the northeast margin of the Yellowstone caldera, (ii) a major gravity low with steep bounding gradients and an amplitude regionally atypical for the elevation of the plateau, (iii) an aeromagnetic low reflecting extensive hydrothermal alteration and possibly indicating the presence of shallow material above its Curie temperature, (iv) only minor shallow seismicity within the caldera (in contrast to a high level of activity in some areas immediately outside), (v) attenuation and change of character of seismic waves crossing the caldera area, and (vi) a strong azimuthal pattern of teleseismic P-wave delays, strongly suggest that a body composed at least partly of magma underlies the region of the rhyolite plateau, including the Tertiary volcanics immediately to its northeast. The Yellowstone field represents the active end of a system of similar volcanic foci that has migrated progressively northeastward for 15 million years along the trace of the eastern Snake River Plain. Regional aeromagnetic patterns suggest that this course was guided by the structure of the Precambrian basement. If, as suggested by several investigators, the Yellowstone magma body marks a contemporary deep mantle plume, this plume, in its motion relative to the North American plate, would appear to be "navigating" along a fundamental structure in the relatively shallow and brittle lithosphere overhead. The concept that a northeastward-propagating major crustal fracture controls the migration path of the major foc

Wyoming