Mark E. Gettings

Saudi Arabian seismic-refraction profile: A traveltime interpretation of crustal and upper mantle structure

Abstract The crustal and upper mantle compressional-wave velocity structure across the southwestern Arabian Shield has been investigated by a 1000-km-long seismic refraction profile. The profile begins in Mesozoic cover rocks near Riyadh on the Arabian Platform, trends southwesterly across three major Precambrian tectonic provinces, traverses Cenozoic rocks of the coastal plain near Jizan, and terminates at the outer edge of the Farasan Bank in the southern Red Sea. More than 500 surveyed recording sites were occupied, and six shot points were used, including one in the Red Sea. Two-dimensional ray-tracing techniques, used to analyze amplitude-normalized record sections indicate that the Arabian Shield is composed, to first order, of two layers, each about 20 km thick, with average velocities of about 6.3 km/s and 7.0 km/s, respectively. West of the Shield-Red Sea margin, the crust thins to a total thickness of less than 20 km, beyond which the Red Sea shelf and coastal plain are interpreted to be underlain by oceanic crust. A major crustal inhomogeneity at the northeast end of the profile probably represents the suture zone between two crustal blocks of different composition. Elsewhere along the profile, several high-velocity anomalies in the upper crust correlate with mapped gneiss domes, the most prominent of which is the Khamis Mushayt gneiss. Based on their velocities, these domes may constitute areas where lower crustal rocks have been raised some 20 km. Two intracrustal reflectors in the center of the Shield at 13 km depth probably represent the tops of mafic intrusives. The Mohorovicic discontinuity beneath the Shield varies from a depth of 43 km and mantle velocity of 8.2 km/s in the northeast to a depth of 38 km and mantle velocity of 8.0 km/s depth in the southwest near the Shield-Red Sea transition. Two velocity discontinuities occur in the upper mantle, at 59 and 70 km depth. The crustal and upper mantle velocity structure of the Arabian Shield is interpreted as revealing a complex crust derived from the suturing of island arcs in the Precarnbrian. The Shield is currently flanked by the active spreading boundary in the Red Sea.

Texture Analysis for Automated Classification of Geologic Structures

Texture present in aeromagnetic anomaly images offers an abundance of useful geological information for discriminating between rock types, but current analysis of such images still relies on tedious, human interpretation. This study is believed to be the first effort to quantitatively assess the performance of texture-based digital image analysis for this geophysical exploration application. We computed several texture measures and determined the best subset using automated feature selection techniques. Pattern classification experiments measured the ability of various texture measures to automatically predict rock types. The classification accuracy was significantly better than a priori probability and prior weights-of-evidence results. The accuracy rates and choice of texture measures that minimize the error rate are reported

A predictive penetrative fracture mapping method from regional potential field and geologic datasets, southwest Colorado Plateau, U.S.A.

Some aquifers of the southwest Colorado Plateau, U.S.A., are deeply buried and overlain by several impermeable units, and thus recharge to the aquifer is probably mainly by seepage down penetrative fracture systems. This purpose of this study was to develop a method to map the location of candidate deep penetrative fractures over a 120,000 km2 area using gravity and aeromagnetic anomaly data together with surficial fracture data. The resulting database constitutes a spatially registered estimate of recharge location. Candidate deep fractures were obtained by spatial correlation of horizontal gradient and analytic signal maxima of gravity and magnetic anomalies vertically with major surficial lineaments obtained from geologic, topographic, side-looking airborne radar, and satellite imagery. The maps define a sub-set of possible penetrative fractures because of limitations of data coverage and the analysis technique. The data and techniques employed do not yield any indication as to whether fractures are open or closed. Correlations were carried out using image processing software in such a way that every pixel on the resulting grids was coded to uniquely identify which datasets correlated. The technique correctly identified known deep fracture systems and many new ones. Maps of the correlations also define in detail the tectonic fabrics of the southwestern Colorado Plateau.

Monte Carlo Simulations of Product Distributions and Contained Metal Estimates

Estimation of product distributions of two factors was simulated by conventional Monte Carlo techniques using factor distributions that were independent (uncorrelated). Several simulations using uniform distributions of factors show that the product distribution has a central peak approximately centered at the product of the medians of the factor distributions. Factor distributions that are peaked, such as Gaussian (normal) produce an even more peaked product distribution. Piecewise analytic solutions can be obtained for independent factor distributions and yield insight into the properties of the product distribution. As an example, porphyry copper grades and tonnages are now available in at least one public database and their distributions were analyzed. Although both grade and tonnage can be approximated with lognormal distributions, they are not exactly fit by them. The grade shows some nonlinear correlation with tonnage for the published database. Sampling by deposit from available databases of grade, tonnage, and geological details of each deposit specifies both grade and tonnage for that deposit. Any correlation between grade and tonnage is then preserved and the observed distribution of grades and tonnages can be used with no assumption of distribution form.

New method to integrate remotely sensed hydrothermal alteration mapping into quantitative mineral resource assessments

Hydrothermal alteration data mapped using the Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER) were compiled into hydrothermal alteration polygons for use in an assessment of porphyry copper mineral resource potential in the southwestern United States. Hydrothermal alteration polygons along with geochemistry, gravity and magnetic, lithologic, and deposit and prospects data were compiled in a GIS to produce a quantitative set of physical properties for each polygon that were effectively used in making estimates of undiscovered deposits for each permissive tract. Results show a higher estimate of potential undiscovered deposits (17 vs 14) for permissive tracts when ASTER alteration data were used in the assessment.

Magnetic susceptibility data for some exposed bedrock in the western conterminous United States

..........................................................................................................................................................