Variations in the geomagnetic field are shown in aeromagnetic survey maps, in which magnetic anomalies are interpreted as being the result of fossilization of magnetic rocks. Thus, rocks exhibit the characteristics of magnetic fields. This observation is often explained using magnetic profiles and maps of magnetic lineaments  . Magnetic minerals can be mapped from the surface to greater depths in crustal rocks depending on the dimension, shape, and magnetic properties of the rock  . Cameroon is underlain by Precambrian rocks, Cretaceous sediments and Cenozoic sediments and volcanic formations  . Meso and Neoprotozoic rocks are found in the Southeastern part of the country. The Cameroon rock basement is divided into two units: the Congo Craton in the South and the Central African Mobile Zone (CAMZ) in the North. CAMZ is a domain of the Pan-African  . It consists of micaschist, plagioclase, bearing and micaceous gneisses, and migmatites intruded by quartz diorite and granodiorite. The study area is in the southern region of Cameroon (Central Africa). It is situated in the northern hemisphere, between the meridians 13˚50' and 14˚20' of longitudes East and parallels 2˚10' and 2˚35' of latitudes North, with an average altitude of 850 m. In a geoelectrical study carried out to the north of this region, it is showed that the intense activities of gold washers encountered in the studied area do attest to the presence of clay mineral concentrations   . To the east of this study area, recent geophysical study  has underlined the major network of lineaments identified in the area under study, which extend from SW to NE with an ENE-WSW major trend and can also be connected to those identified NE-SW by  . The aim of this work was to interpret aeromagnetic data for geological and structural mappings in the southeastern Cameroon. This study also focused on the determination of certain mineral resources potential in the area.
Cameroon is a part of CAMZ, with a geological basement of sedimentary rock. These sedimentary rocks are restricted to the southwestern and northern part of the country. In southwest Cameroon, the oldest sedimentary rocks are massive cross-bedded sandstones and conglomerates  . The sandstones are overlained by fossiliferous shales of early upper Cretaceous series ages  . The geological mapping of the study area (Figure 1) was first done by  . The investigated area lied in the southern part of Cameroon, is covered by Neoproterozoic formations and includes rocks resulting from the intermediate series, made up of schist and quartzite formed following an epizonal metamorphism   .  considered the Ayos, Mbalmayo and Bengbis formations as being of a Precambrian mean age. The area is also included in the Dja inferior series and is comprised of the tillitic, doleritic and Bek complexes. It belongs indeed to the meta-sediments series of Dja-Ayos-Mbalmayo-Bengbis  . The geology of the area is dominated by the extension of the Archean Congo craton (Ntem Complex ~3 Ga). The composition is made up of chlorite-greenschist, mica-schist with muscovite, and the lentilles of quartzite interstratifies. The doleritic complex (sills and dykes) is composed of greenschist, the mylonitics of quartzites, the dolerites and the gabbros. The Bek complex is made up of sandstones of quartzite ochre, clay-schist and blue schist. The base complex series of garnet is mainly composed of ectinites (gneiss and amphibolite), migmatites (granite of magmatic) and intrusive rocks (granite and granodiorite). The study area is a part of the Congo Craton (CC). Its formation dates to the rejuvenation during the Pan-African orogeny about 550 Ma    . That major tectonic features lie under geological covers. Geophysical studies show that the boundary of the CC and the Pan-African is around 4˚N    . An aeromagnetic study done by  , estimates the northern boundary of the Congo Craton as starting from 3˚7'N of West to 3˚75'N of East. Its depth is estimated around 2.6 km for deep and 0.1 km for shallow while the direction is ENE-WSW and the NW slope varies from 30˚ to 60˚. Also, an audio-magnetotelluric study carried
Figure 1. Geology map (modified).
out by  shows many discontinuities in the topography of the subsurface. This topography presents a major deep-seated fault with E-W direction. Also,  showed the morphological difference and the tectonically subdivision into two tectonic sectors corresponding to the Congo Craton in the south, Pan-African in the north and helped identify the tectonic boundary separating them at depth. The tectonic features (faults and folds) are fitted into the CAFB’s deformation history and could be due to the Trans-Saharan east-west collision system  . Some buried faults have been confirmed in geophysical studies   . The Archaean metamorphic rock assembly has an E-W trending foliation (S1) and intrusions that are usually aligned along a NE-SW trending foliation (S2). The S1 foliation has near vertical dips to the north and is locally deformed into mesoscopic isoclinal D2 folds. The S2 foliation is a regional, steeply dipping planar fabric with variably oriented stretching lineation and large-scale open folds that are associated with N-S trending sinistral and dextral strike-slip faults and mylonitic (S3) foliation. The S2 foliation is well observed in the greenstone units of the Ntem complex and its development is linked to dome-and-basin tectonics related to diapiric movements  .
2. Materials and Methods
2.1. Data Source and Acquisition
The aeromagnetic data set used in this study was from the recent acquisition by the Compagnie Minière du Cameroun SA (CMC SA) on 2012. Airborne geophysics with 200 m line spacing was conducted on Upper Nyong Division of the East Region of Cameroon in February 2012. The Survey was carried out by an AS350B2 helicopter, with a total of 11˚305 lines km of magnetic data and the flight was NE-SW trending. The survey was carried out at a nominal terrain clearance of 100 m which was monitored by a radar altimeter with an accuracy of ±20 m.
2.2.1. Total Magnetic Intensity Reduction to the Equator and Its Residual Anomaly
After extraction, the aeromagnetic data survey was corrected and controlled on Excel. The anomaly of the total magnetic field (TMI), is the difference between the magnetic field extracted International Geomagnetic Reference Field (IGRF) BRef to each station at the date February 15, 2012. This magnetic anomaly is given by (1):
The TMI map was reduced to the equator (RTE) (2) to avoid the problems associated with low-latitude magnetic data  . RTE is a complementary filter to the Reduction-to-Pole (RTP) procedure. Any magnetic anomaly point depends on the inclination and the declination of the main magnetic field of the earth.
It transforms an anomaly of non-zero inclination into an anomaly that would be observed on the same body with zero inclination. The RTE can be expressed as:
where I is the geomagnetic inclination, D is the geomagnetic declination, sin I is the amplitude component, and is the phase component.
The software Geosoft of Oasis Montaj v.8.4 using the convolution of Fourier transformation made it possible to obtain the map of TMI reduced to the equator and other relevant maps. To realize this map, this computed the values of the Inclination (I) and the Declination (D). These values are: I = −22.77˚ and D = −1.08˚ respectively.
The residual anomaly ( ) is calculated in each point of the regular grid by taking the difference between the anomaly of TMI RTE ( ) and Upward Continue to TMI RTE at 2 km ( ). It is given by (3):
2.2.2. Tilt Angle Approach
Also known as Variation of inclination, the tilt angle is a transformation which includes the first vertical derivative and the module of the first horizontal derivative of the anomaly of residual TMI reduced to the equator. The advantage of the tilt angle is that, compared with the other methods, it does not require the knowledge of parameters such as (density, magnetic susceptibility, structural index etc.). The other advantages of this transformation of the tilt angle are that, by computing an angle, all the forms are represented in a similar way; such that the anomaly either has minimum or maximum amplitude and that it also allows mapping the features with high resolution. Indeed, the arc-tangent function has as effect, to distribute the signal calculated between −90˚ and +90˚. The tilt angle applied to the anomaly of the magnetic field reduced to the equator permits to estimate the depth of the upper end of the sources. According to some authors    , it is given by (4):
where . Putting and , then . Thus for and for
In other words, the estimated depth of the upper end of the source is obtained by measuring the perpendicular distance between contours θ = 0rad and
to the feature. The map of tilt angle is obtained by using the software Geosoft of Oasis Montaj v.8.4. This transformed map presents a correlation (or analogy) with the geological structure of the basement in this study area with a description of the zones of contact of the geological formations.
2.2.3. Euler Deconvolution Approach
Euler deconvolution   is a method of estimating the depth and slopes of subsurface magnetic anomalies and can be applied to any homogeneous field of magnetic data  . In Euler’s method, the structural index must be assumed as prior information because the quality of the depth estimation depends mainly on the choice of the proper structural index, which is a function of the geometry of the causative bodies and characterizes the rate of the variation of the anomaly intensity with a distance. Some authors   showed that the optimum structural index usually yields the tightest clustering of the solutions. The depth estimates from magnetic data are more accurate with the equator-reduced magnetic field. The estimated Euler sources for the single-points are based on computing Euler’s homogeneity Equation (4) and result in clusters used to constrain the overall geometry of the model. The 3D form of Euler’s Equation (5) can be defined  as:
where B is the regional value of the total magnetic field and (x0, y0, z0) is the position of the magnetic source, which produces the total magnetic field F measured at (x, y, z). N is the Structural Index on how to characterize the source.
Therefore, it has assigned a value of 1.0 as a structural index to locate the possible magnetic contacts because it is particularly good at delineating the subsurface contacts. An overlapping moving window is 10 km by 10 km, a tolerance of 10% for Euler residual field and 15% for Euler Tilt angle and a proportioned symbol base of 100.
3. Results & Discussion
3.1. Anomaly of Total Magnetic Intensity in the Studied Area
The total magnetic field (TMI) is the response produced by rocks containing magnetic minerals. In Figure 2(a) & Figure 2(b), the TMI anomalies are shown as 2D and 3D maps, varying from −495.7 nT to 93.7 nT and are distributed in a bipolar way on the study area. The correlation between the TMI anomalies map and geological contacts is weak. The magnetic bandings which leads to magnetic bifurcation and a quasi-horizontal gradient at YANBOT are notices. On either side of the DJADOM and ETA axes, present the quasi-vertical gradient. At certain place, appear peak in the strongly magnetics formations. These observations underline the intense activities between the CC and the Pan-African. The DJADOM axis shows a lateral extension and longitudinal heterogeneous anomalies dominated by the positive anomalies, with a maximal value of 93.7 nT and a long wavelength of 33176.9 m.
Figure 2. (a) Anomaly map of total magnetic Intensity of the study area. (b) Three-dimension Anomaly of total magnetic intensity map of the study area with even color bar.
Geologically, the preceding observations are due to sandstone ochre quartz and the schist of the Bek complex, the dolerite of the doleritic complex, and the silver micaschist and ore quartzite in the base complex. On the TMI anomalies map, several places show high susceptibility contrasts, which is an indication of strong magnetization. Geological indicators point to inferred magnetite, dolerite and ochre schist quartzite which have a strong magnetization in this zone. The presence of weakly magnetized anomalies would be due to the migmatites of the base complex series.
3.2. Anomaly Total Magnetic Intensity Map Reduced to Equator
The anomaly total magnetic intensity map reduced to equator (TMI-RTE) in Figure 3 shows the anomalies varying from -465.8 to 134.5 nT. Compared to the TMI, the anomalies preserve their forms. Positive anomalies show the response of a strongly magnetized base. The longer-wavelength anomalies reflect deeper magnetic responses, whereas shortest-wavelength anomalies cause shallow magnetic responses  . The greater the wavelength, the deeper the response. These features characterize the earlier mentioned geological formations. The positive anomalies in the ESE-WNW and ENE-WSW directions tend to increase in amplitude and wavelength. The earlier mentioned gradients and the peaks are noticed here. These observations highlight the event of subduction of the CC onto the Pan-African.
The strong circular anomalies in this zone show the presence of accumulated minerals of strong susceptibilities in the major faults. This is due to the presence on the one hand of diamagnetic minerals such as ochre quartz and on the other hand, ferromagnetic minerals of strong magnetization remnants such as the ectinites, dolerites in the series of the basic complex and even magnetite.
3.3. Total Magnetic Intensity Residual of the Study Area
The values of the anomalies vary from −204 nT to + 121.7 nT, with a 12.8 nT reduction compared to the total field RTE. The map of residual anomaly in Figure 4 shows the inferred magnetic bodies. The effects of the surface structures are masked by those of the underlying structures. The magnetic anomalies of the DJADOM axis have disappeared. The distributions of these anomalies are more refined and indicate the local maxima of the ESE-WNW and ENE-WSW directions as previously noted on the map of the TMI and its transforms. This
Figure 3. Anomaly map Reduced to Equator of the total magnetic intensity of the study area.
Figure 4. Anomaly of the residual Total magnetic intensity map of the study area.
orientation of positive anomalies makes it possible to identify the directions of major structures and to locate them. The magnetic peaks and the bandings are observed to burst and break in several areas, and these reflect the distribution of magnetization from the top to the bottom soil. The observed magnetization could be attributed to the intrusion of banded iron, magnetite, dolerite, mylonitic quartzite with dolerite and amphibole in the northern part and to the area situated southwest of DJADOM, then BIF, schist and sandstone quartzite ochre, and diamond in the southern part. The opening formation accompanied by virgations to the southwest of DJADOM is remark, giving rise to two great geological undulations. Since the amplitudes and wavelengths are maximal, enable to say that the geological undulations go from the covers to the base.
The structural map leads to put in evidence a network of faults in the study area and thus to show that the network of faults in the adjacent Eastern zone of this area is prolonged into the current zone  . The continuity of this prolongation Northwards of this zone is highlighted in the audio-magnetotelluric  and in aeromagnetic  studies. The geometrical description of this structure suggests an open synclinal transposed on vertical foliations: the major fault at the DJADOM axis is quasi-parallel to the Northern limit of the CC and parallel to the Sanaga Fault (SF) and the Central Cameroon Shear Zone (CCSZ). The directions of the lineaments are show in Figure 9.
3.4. Variation of Inclination Angles (Tilt Angle)
The Tilt angle of the residual anomalies of TMI shown in Figure 5(a) and the amplitudes of the angles vary between −1.4 rad and + 1.3 rad while the lengths also vary. The tilt angle reflects the surface and subsurface distribution of structural contacts. This map shows several features on the subsurface which are not visible on the geological map. The interpretation of the aeromagnetic data shows
Figure 5. (a) Variation of magnetic inclination; (b) Lineaments of variation of magnetic inclination.
several magnetic features in this zone of study. The previous maps permitted us to observe Breakings. The positive magnetic anomalies are burst and laid out according to the features. The structures in the E-W quasi direction meet a shock the southwest of DJADOM, and then start on a curve to take the ESE-WNW, E-W and ENE-WSW directions. This curve of the structural directions indicates a fold of drive, revealing a sinister movement following the WNW-ESE direction, and dextral in the SWS-NEN direction. These can be interpreted as fractures. In this zone where these fractures are more significant (DJADOM and ETA axes), notice the compartments with fractures. These zones of significant fractures are considered as delimiting significant features which refer to as the accident of the ETA for the ESE-WNW direction.
According to the amplitude and the length (Figure 5(b)), 39 major features and 397 minor features, making a total of 436 are noted. These principal features have lengths varying from 1731.9 m to 33 176.9 m, and depths going from 5.0 m to 238.8 m. The characteristics (lengths, directions and depths) of the principal features are grouped in Table 1. The minor lineaments show that the base was affected by the tectonic events which characterize the transition between the zone from the Congo Craton and the belt from the folds of the Pan-African. Generally, foliations are of quasi-parallel texture. These foliations are controlled by Eburnean oogenesis and confirm those mentions by  .
3.5. Euler Deconvolution
3.5.1. Euler Deconvolution of Residual Field
The Euler’s solutions enable the characterization of magnetic responses and determination of the depth and geometry of intruding bodies  . The Euler’s solution (Figure 6) has as structural index N = 1, tolerance T = 10% and Nyquist Window W = 10 km × 10 km. The depth of Euler’s solution varies from 13 m to 314 m, including those of the bodies with their geometries. The interpretation of the structural map highlights the various features affecting this study area. The deep features and the limit of the intrusive bodies are distinguished. The great directions of lineaments of this area are: NE-SW, ENE-WSW, E-W and
Figure 6. Euler’s solution of TMI residual N = 1, T = 10%, W = 10 km × 10 km.
Table 1. Direction of major features, their lengths and their depths resulting from Tilt angle (depth 1), resulting from the Euler’s solution of VMI (depth 2).
ESE-WNW. The important lineament, earlier referred to as the ETA geological accident is observed on this map. Their depths vary from 13.3 m to 313.1 m.
3.5.2. Euler Deconvolution Variation of the Inclination (Tilt Angle)
Euler’s solution shown in Figure 7 presents a depth varying from 0.7 m to 238 m. More details on the depth of the intruding geological structures in the base are shown. The solution is used for the final interpretation of the geological contacts associated to the features noted previously. The extent of the number of structures on the map highlights the intense tectonic activity that this area undergoes. This is the reason of the metamorphism of the cover, renovated by the internal stress at the time of the collision.
Effectively, the series of solutions are observed which indicate the major fault at ETA. The depths of the structures in this map are proof of the fact that this study area belongs to the mobile zone of the central African mobile belt (CAMB) as stated in gravimetric investigations  . The quasi-discordance of contacts is confirmed by the disposition of Euler’s solutions (Figure 6).
Figure 7. Euler’s solution of VIM N = 1, T = 15%, W = 10 km × 10 km.
3.6. Correlation of Tilt Angle Structures and Its Euler’s Deconvolution Map
The interpretation done by superimposing maps of tilt angle and Euler’s is based on criteria elaborated by many authors     .
If the features of the Variation of the Magnetic Inclination (VMI) and those of the Euler solution are almost parallel and do not merge, the VMI features represent the contacts while those of Euler’s solutions indicate the sense and direction of their slopes for these contacts and the overlapping of the ores located in these features.
If the features of the VMI and those of the Euler’s solution are almost parallel and are merged, the VMI features represent the contacts while those of the Euler solution indicate the direction and sense for the vertical slopes of these contacts and the vertical gradients affecting the base.
The correlation map (Figure 8) shows the superposition of the features resulting from the VMI and the Euler’s solutions which are almost parallel, and their slopes are like those obtained within the framework of the geological study. There are several principal lineaments with vertical contact. These slopes are dextral in the Center and South, and sinistral in the North.
The main results obtained in this study highlight new elements which enable the improvement of the knowledge on the structural features of the studied area. The tilt angle method is used to delineate geological contacts and structures and to estimate the depth and the length.
The bifurcation, accompanied with virgations, given the fact that there are many faults with WNW-ESE, W-E and WSW-ENE major directions. There is
Figure 8. Euler’s solution of VIM N = 1, T = 15%, W = 10 km × 10 km.
Figure 9. (a) Rosacea diagram of the major lineament orientations; (b) Rosacea diagram of the minor lineaments orientations.
a faulting system in this study area. Principal lineaments are determined with the major direction being WNW-ESE, while that for the minor lineaments is W-E (Figure 9). Two major faults ESE-WNW and ENE-WSW are highlighted with the dominant direction being ESE-WNW, which could be named the geological accident of ETA. This shows the series of deformations with a principal ESE-WNW direction under conditions which are realized in extreme cases of the metamorphic state, which according to have a rather recent deformation. This study enables a better knowledge of the subsurface structure of this area. The intrusion of dolerite, mylonitic quartzite of dolerite, magnetite, quartzite schist and sandstone quartzite of ochre are remarked. A geophysical exploration of this zone would enable us to better highlight the preceding observations.
The authors are grateful to the reviewers and the readership for their kind help in making the manuscript clearer, more correct and mature for publication.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
 Toteu, S.F., Penaye, J. and Poudjom Djomani, Y. (2004) Geodynamic Evolution of the Pan-African Belt in Central Africa with Special Reference to Cameroon. Canadian Journal of Earth Sciences, 41, 73-85.
 Gouet, D.H., Meying, A., Assembe, S.P. and Ndougsa-Mbarga, T. (2013) Gold Mineralization Channels Identification in the Tindikala-Boutou Area (Eastern-Cameroon) Using Geoelectrical (DC & IP) Methods: A Case Study. International Journal of Geosciences, 4, 643-655.
 Gouet, D.H., Meying, A., Assemble, S.P. and Ndougsa-Mbarga, T. (2015) Clay Minerals Channels Identification in the Tindikala-Boutou Area (Eastern-Cameroon) along the Kadey River Using Direct Curent (DC) Method. Journal of Geoscience and Environment Protection, 3, 123-133.
 Feumoe, A.N.S., Ndougsa-Mbarga, T., Manguelle-Dicoum, E. and Derek Fairhead, J. (2012) Delineation of Tectonic Lineaments Using Aeromagnetic Data for the South-East Cameroon Area. Geofizika, 29, 175-192.
 Van Den Hende (1969) Notice Explicative sur la Région Batouri-Ouest Ouest de la Carte Géologique de Reconnaissance à l’échelle du 1/500.000, Direction des Mines et Géologie, Yaoundé, Cameroun, 23.
 Djoukan F.M.L., Shüssler, U. and Tchouankoue, J.P. (2004) Geochemistry and Petrogenesis of Neoproterozoicgranitoids within the Central African Fold Belt in the Bafoussam Area, Western Cameroon (Central Africa). University of Wuerzburg, Germany, 4.
 Ndougsa-Mbarga, T., Yufenyiu Layu, D., Quentin Yene-Atangana, J. and Tabod C.T. (2014) Delineation of the Northern Limit of the Congo Craton Based on Spectral Analysis and 2.5D Modeling of Aeromagnetic Data in the Akonolinga-Mbama Area, Cameroon. Geofísica Internacional, 53, 5-16.
 Tadjou, J.M., Nouayou, R., Kamguia, J., Kande, H.L. and Manguelle-Dicoum, E. (2009) Gravity Analysis of the Boundary between the Congo Craton and the Pan-African Belt of Cameroon. Austrian Journal of Earth Sciences, 102, 71-79.
 Meying, A., Ndougsa-Mbarga, T. and Manguelle-Dicoum, E. (2009) Evidence of Fractures from the Image of the Subsurface in the Akonolinga-Ayos Area (Cameroon) by Combining the Classical and the Bostick Approaches in the Interpretation of Audio-Magnetotelluric Data. Journal of Geology and Mining Research, 1, 159-171.
 Yandjimain, J., Ndougsa-Mbarga, T., Meying, A., Bi-Alou, M.B., Ngoumou, P.C., Assembe, S.P., Ngoh, J.D. and Owono-Amougou, O.U.I. (2017) Combination of Tilt-Angle and Euler Deconvolution Approaches to Determine Structural Features from Aeromagnetic Data Modeling over Akonolinga-Loum Area (Centre-East, Cameroon). International Journal of Geosciences, 8, 925-947.
 Mvondo, H., Owona, S., Mvondo-Ondoa, J. and Essono, J. (2007) Tectonic Evolution of the Yaoundé Segment of the Neoproterozoic Central African Orogenic Belt in Southern Cameroon. Canadian Journal of Earth Sciences, 44, 433-444.
 Ndougsa-Mbarga, T., Manguelle-Dicoum, E., Tabod, C.T. and Mbom-Abane, S. (2003) Modélisation d’anomalies gravimétriques dans la région de Mengueme-Akonolinga (Cameroun). Science, Technology and Development, 10, 64-74.
 Grauch, V.J.S., Bauer, P.W. and Kelson, K.I. (2004) Preliminary Interpretation of High-Resolution Aeromagnetic Data Collected near Taos, New Mexico. New Mexico Geological Society, Guidebook, 244-256.
 Salem, A., William, S., Fairhead, D., Ravat, D. and Smith, R. (2007) Tilt-Depth Method: A Simple Depth Estimation Method Using First-Order Magnetic Derivatives. The Leading-Edge December, Meter Reader, Vol. 150, 2-5.
 Verduzco, B., Fairhead, J.D., Green, C.M. and Mac-Kenzie, C. (2004) New Insights into Magnetic Derivatives for Structural Mapping. The Leading Edge, 23, 116-119.
 Miller, H.G. and Singh, V. (1994) Potential Field Tilt a New Concept for Location of Potential Field Sources. Journal of Applied Geophysics, 32, 213-217.
 Reid, A.B., Allsop, J.M., Granser, H., Millett, A.J. and Somerton, I.W. (1990) Magnetic Interpretation in Three Dimensions Using Euler Deconvolution. Geophysics, 55, 80-90.
 Durrheim, R.J. and Cooper, R.J. (1998) EULDEP: A Program for the Euler Deconvolution of Magnetic and Gravity Data. Computer and Geosciences, 24, 545-550.
 Langenheim, V.E. and Jachens, R.C. (2014) Aeromagnetic Data, Processing, and Maps of Fort Irwin and Vicinity, California. In: Buesch, D.C., Ed., Geology and Geophysics Applied to Groundwater Hydrology at Fort Irwin, U.S. Geological Survey Open-Fort Report 2013-1024, Chap. 1.
 Feumoé and Ndougsa-Mbarga (2017) Curie Point Depth Variations Derived from Aeromagnetic Data and the Thermal Structure of the Crust at the Zone of Continental Collision (South-East Cameroon). Geophysica, 52, 31-45.
 Phillips, J.D. (1998) Processing and Interpretation of Aeromagnetic Data for the Santa Cruz Basin—Patahonia Mountains Area, South-Central Arizona. U.S. Geological Survey Open-File Report, Arizona, 02-98.