ABSTRACT

Gravity data are sensitive to local vertical offsets across high-angle faults, where rocks with different densities are juxtaposed. Yet high densities in some Mesozoic sedimentary rocks just above the basement may smear out the subtle gravity signatures of basement faults. At this study the gravity data processing tends to avoid ill-described “black-box” techniques. The study area is situated in the Palomas site, Cuatrociénegas region, Coahuila, NE Mexico. The San Marcos Fault is at least 300 km long and has WNW-ESE trend from the central part of Nuevo León State through Coahuila, and finally to the eastern part of Chihuahua State. Gravimetric data shows that the lowest values of free air and Bouguer anomalies are in the southern part of the area, and the highest values are in the western and central part of the area. Between these parts exists a zone of high horizontal gravity gradient. Configuration of linear elements of gravity field (gradient zones) delimited the San Marcos Fault in the San Marcos valley below thickness of recent sedimentary cover. Two density models were carried out, which showed that the Cretaceous rocks are in discordant contact with the Paleozoic rocks that can be related to the San Marcos Fault. The density was determinate using to Nettleton’s method, which results highlight the presence of the San Marcos Fault. Density models showed that the Quaternary sediments are in direct contact with the San Marcos Fault.

1. Introduction

Steep, straight faults are commonly expressed as subtle potential-field lineaments, which can be gradient zones, alignments of separate local anomalies of various types and shapes, aligned breaks or discontinuities in the anomaly pattern, and so on [1-4]. Many large magnetic and gravity anomalies represent the ductile, ancient, healed basement structures, obscuring the desirable subtle features [1-8].

Subtlety of the desirable lineaments necessitates detailed data processing, using a wide range of anomalyenhancement techniques and display parameters. The final choice of processing steps depends on which aspects of the anomalies one aims to enhance, as well as on experimentation with various techniques [1,3,5,9,10].

Not a panacea, data processing is a necessary evil. Because the signal and noise anomaly characteristics commonly overlap, complete separation between them may be impossible. Noise artifacts may actually be introduced [3].

Unexpected consequences and processing side-effects are normal. As well, it may be hard to know in advance which of the many anomalies are desirable.

The processed and enhanced anomalies should ideally be easy to relate back to the original anomaly shapes. We kept the data processing to a minimum, avoided ill-described “black-box” techniques, and relied on mathematically simple and intuitive procedures [3,5-11].

The San Marcos Fault (SMF), was defined by Charleston [12,13], and is outcropping in Central Coahuila in northeast Mexico correspond to one of the three major zones of west-northwest trending lineaments interpreted as basement faults that controlled the tectono-sedimentary evolution since the Triassic-Jurassic in the region [14,15] (Figures 1 and 2), the others structural lineaments are: The La Babia fault (LBF) [12,13], located in Northern Coahuila; and the Mojave-Sonora Megashear (MSM) [14-16]. From these basement faults, only the SMF has a

trace well defined from stratigraphic and structural evidence, which can be followed since Potrero La Gavia, through San Marcos Valley (SMV), to Potrero Colorado and to the Sierra Mojada [17-21]. These authors mentioned that SMF is distinguished by: contrasting structural styles between the Coahuila block and the Coahuila folded belt; the existence of a pre-Tithonian and Neocomian major sedimentary clastic wedge related to its activity and deposited in the hanging wall of the fault; and the existence of fault contacts between Permian, Jurassic, and Cretaceous rocks in the SMV and Potrero Colorado areas (Figure 2). In areas where the SMF is buried, such as the Camargo volcanic field has been inferred from structural features in the surface [22-24].

The SMF, LBF and the MSM play an important role to understand the tectonic and stratigraphic evolution of northeast Mexico. Several authors [14-16,23,25,26] propose that these fault systems (SMF and LBF) were reactivated at different times. McKee et al. (1984) found evidence of movement of the San Marcos fault during the Jurassic, Early Cretaceous, and Tertiary. Additionally, Aranda-Gómez et al. (2005) suggest that SMF have been reactivated during the Late Pliocene-Quaternary by normal faulting in the border of Chihuahua-Coahuila states. In central Coahuila not has been documented Quaternary activity, due to the works published have been focused mainly in the study of pre and syn-laramide reactivation of the fault. As mentioned before, the study of the SMF is related only with geologic information, geophysics data not are available concerning to characterization of this basement fault. One goal of interest about the possible reactivation of the SMF at Plio-Quaternary time is the seismic risk associated with potential earthquakes related to this fault. Northeast Mexico is characterized by low seismicity and a lack of strong ground motion records [27-30].

Large intraplate earthquakes in the relatively stable interior regions of continents are rare in comparison with those that originate in plate boundary regions. However, these occasional events can be extremely devastating, because cities located in continental interior regions are often built without seismic design criteria. Crone et al. (2003) mention that major intraplate earthquakes can cause widespread damage because the attenuation of seismic energy from large earthquakes is relatively low in plate interiors. Galván-Ramírez and Montalvo-Arrieta (2008) made a compilation of the historical seismicity in northeast Mexico; they found that some of the epicenters locations are overlaped or near to the traces of the major regional faults trending west-northwest postulated in northeast Mexico: LBF, SMF, and MMS. Figure 3 shows the historic seismicity for northeast Mexico. These authors hypothesize that seismic activity reported for this region, which lies near or on the trace of these faults,

could be related to reactivation of these major faults. However, to probe this hypothesis is necessary to install a seismic network to confirm or refute the presence of seismicity associated with these fault systems. Additionally, these authors mentioned that a possible critical scenario would represent a rupture (M_W = 6.5) in the south segment of the San Marcos fault in Central Coahuila. The importance of this scenario is the settlement of three of the most populated centers in northeast Mexico (Monterrey, Saltillo and Monclova with a total population of more than six million) located in a radius less than 150 km from the fault source.

Some evidence of recent seismicity could be obtained from visual observation of deformation of latest sediments covering fault zone. One way to identify buried faults systems in sedimentary valleys is by means of gravity data. The goal of this study is identifies the SMF zone in the SMV in central Coahuila from gravity data and have the first stage to identify if there is evidence of displacements in Quaternary sediments. To develop thiswe obtain for first time precision gravity data for five north-south trending lines across the Palomas area (Figure 4) in SMV.

2. Geologic Framework

The Palomas area is located to the NW of San Marcos y Pinos valley and to the SE of Sierra El Granizo, between 26˚24'N and 26˚26'N, and 101˚50'W and 101˚57'W (Figures 2 and 4). In this area there are some isolate hills composed by limestone of Upper Tamaulipas Formation that can be associated to displacement of the SMF, we recollected gravity data distributed in five profiles, and two of them cross the La Pasta hill (Figure 4).

SMV is located in central Coahuila and represents an area where there is excellent evidence of the SMF [14,19, 31].

The SMV are compose by pre-Jurassic outcrops composed by small lenses of granite rocks at NW of the valley and slates and metasandstones in fault contact with both Cretaceous and Jurassic rocks at northwest of the Las Palomas area. The Jurassic sediments composed by conglomeratic sandstone of Las Palomas beds and sandstone of the Sierra El Granizo beds are exposed along the foothills of the Sierra El Granizo and in the valley where a gradational contact with overlying Neocomian San Marcos Formation is exposed. Tanque Cuatro Palmasbeds consist of nearly 100 m of fine-grained marine sedimentary rocks without conglomeratic units [14]. The Cretaceous sediments are defined by the San Marcos Formation composed by continental alluvial deposits represented by conglomeratic units of pebbles or cobbles of volcanic rocks with subordinate quartz, and plutonic fragments, overlying in a discordant contact is limestone of the Cupido Formation, the sedimentary sequence are complete by shale and marble of the La Peña Formation, follow by limestone of Tamaulipas Superior Formation, overlying by limestone, shale and banded chert of Cuesta del Cura Formation and interbeded shales with marble and gysum veins of Indidura Formation, the marine sequence is overlie discordantly by Quaternary deposits [14, 19,31].

3. Geophysical Data Collection and Reduction

More than 150 gravity reading points grouped in 5 lines were carried out nearly perpendicular to the main trace of the SMF. Gravity stations spacing was from 15 - 20 up to 40 - 50 meters. Gravity meter Autograv CG-5 with 0.001 milligal (mGal) of reading resolution was used (1 Gal = 1 cm/s²). Gravity reference station was established in the Central Place of Cuatrocienegas city.

The error of measurements (RMS error) did not exceed 5 microgal (μGal). Corrections were applied in real

time for tilt errors, for long term drift, for the temperature of the sensor, and for earth tides. The seismic filter was applied too. The location of gravity stations was determined by a GPS Magellan Platinum with a horizontal accuracy of 5 - 7 m.

Elevation of the gravity reading points was measured with a Total Station Nikon DTM-551 with precision of ±2 mm. The standard corrections, e.g. instrumental drift, latitude, elevation using new standards to reducing gravity data were applied to obtained measurements [32]. As a result of field data correction and processing free air and Bouguer gravity anomalies were obtained. Data interpretation procedures include Fourier transformation, wavelength filters, upward and downward continuation, vertical and horizontal derivates, etc. [1,5,9,11,33,34]. A series of maps (Figures 5(a)-(d)) and graphs (Figures 6 and 7) show gravity anomalies interpretation technique. 2D geological-geophysical models were elaborated (Figures 8 and 9).

A sea level datum NAD 83 and the standard density 2.3 g/cm³ were used for the elevation correction. The terrain correction was calculated using the method by Hammer and was applied to each gravity station to obtain the complete Bouguer anomaly.

Gamma ray measurements were carried out additionally to gravity study. Gamma ray reading points were the same as gravity stations. GRS-500 Differential Gamma Ray Spectrometer/Scintillometer was used to measure the gamma radiation emitted by various daughter isotopes in the uranium decay series. Sample rate was selected as 10.0 seconds, and energy window was choosing detecting total count above 0.08 MeV. Observed count rates were corrected due to dead time (4 microseconds).

4. Gravity and Gamma Ray Data Analysis

The observed gravity map (Figure 5(a)) shows an isolated gravity high in the NE part of study area. Southeast of this high, an extensive gravity low is observed. The amplitude of gravity field between these extreme parts reaches 24 milligals (mGal). A high horizontal gradient zone (up to 0.01 mGal/m) is shown here (Figure 5(d)).

The free air gravity map (Figure 5(b)) shows an anomaly that divide two big blocks. The southern block present the lowest values that varying between –15.5 to –14.5 mGal. The northern block shows high gravity values. Anomalies of –9.0 - 12.5 mGal are observed here. The Bouguer anomaly (Figure 5(c)) shows the same behavior that the free air gravity, where can see two main areas. The complete Bouguer gravity (Figure 5(c)) values in the southern part of the study area are in the range of –109 up to –112 mGal and show several weak local anomalies with an amplitude of ±0.2 mGal. Gravity high with a range of –105 up to –108 mGal is shown in the northern part of area. There is a relatively gravity low (up to –110 mGal) between two extreme parts of this positive zone. A strong gravity gradient zone marks a structural boundary between these blocks. It seems that this strong anomaly is related to the main trend of the fault zone in

the area. The values of horizontal gravity gradient here are about 0.0035 - 0.004 mGal/m (Figure 6) and reach 0.035 - 0.04 mGal/m near the main fault zone (Figure 5(d)). The same characteristics have the gamma ray curves, which show a close coincidence with gravity gradient zones. For example, Figure 6 shows two extreme values of horizontal gradient 0.0024 mGal/m and 0.00450 mGal/m and two gamma ray peaks (240 cps and 250 cps) correspond them. As general, observed and transformed gravity anomalies was revealed trending northwestsoutheast across the central part of Las Palomas area. An extensive gravity gradient zone trends in the same direction. This zone probably represents the San Marcos Fault, which crosses the area from northwest to southeast (Figure 2).

The gravity data suggest a complex subsurface structure in the San Marcos Valley area. It seems that two distinct systems of linear elements dominate in the basement as well as shallow structure (Figure 7). The first has a southeast direction and it coincides with the San Marcos Fault main direction, plotting along the general axe of the San Marcos Valley. The second linear system has a nearly perpendicular (north-northwest-south-southeast) direction. It seems that these faults divide the San Marcos Valley area to the several blocks of basement accordingly the main structure of the San Marcos Fault (Figure 7).

5. Modeling 2D of Gravity Data

Figure 7 shows the contour of the Bouguer map and the main geologic features in the study area. In this figure is depicted the fracture arrangement that we interpreted from gravity data too. Two sets of fractures can be interpreted from the gravity data: WNW-ESE and NNW-SSE (Figure 7). It is associated with the isolated calcareous hill (La Pasta) at center of figure, which can be correlated with the trace of the SMF (red discontinuous line).

We also interpreted secondary faults that have the same orientation that SMF and a systems of faults or fractures the cut the main trending system that corresponds to conjugate faults.

From gravity data we constructed two 2D cross sections models for profiles A-A’ and B-B’ (Figure 4). The first one (A-A’) shows the model located in the sedimentary valley (Figure 8). In this figure are identified three units with the follow density values: Red unit with a density of 1.7 g/cm³ related to alluvial deposits. The green unit has a density value of 2.3 g/cm³ correlate with calcareous materials that is in a strong discontinuity or fault contact with the grey unit that is related to material with an increase in density of 2.425 g/cm³ that correspond to Paleozoic flysch. The same figure point up the correlation between the observed and calculate gravity data by the 2D model, which shows a step geometry of the gravitational anomaly; the interpreted source of this anomaly is due to the presence of the SMF. The difference of density for each block varies 0.125 g/cm³. The denser rocks are in the handing-wall of the fault. In this profile the anomalies can be associated to: a) basement up-lift or b) density contrast between rocks. These two circumstances can be interconnected as follows: small negative anomalies can be caused by increase in the thickness of quaternary sediments in those places where the rocks are less strong and less dense; we observed local positives anomalies that can be related to denser rocks denser rocks which are less weathered and thicker than the surrounded county rocks.

Figure 9 depicts the profile B-B’, it cross-cut La Pasta hill which have a height about 40 m and where calcareous rocks are outcropped, this profile has a south-north direction. The observed and calculated gravity show a positive anomaly of 4 mGal in the La Pasta hill and a small negative anomaly nearby –0.5 mGal to the north from this hill which can be related to denser rocks (0.1 mGal). Positive anomaly in area of the La Pasta hill reflects a relative high effective density of 0.23 g/cm³ (density of adjacent rocks is about 2.53 g/cm³). Geologically, this anomaly can be correlated with stronger and denser blocks by mineralization or low grade metamorphism. Negative anomaly can be caused by a thicker package of quaternary sediments. We identified the follows density units (Figure 9): Red unit with a density of 1.7 g/cm³ related to alluvial deposits, the green unit had a density value of 2.3 g/cm³ correlate with calcareous materials, that is in a strong discontinuity or fault contact with the third unit with a density of 2.53 g/cm³, that can be related with the main trace of the SMF that include conglomeratic sandstone of the Las Palomas beds and Paleozoic flysch. To the north of main fault exists another fault or discontinuity that separates the fourth and fifth units that have density values of 2.28 g/cm³ and 2.4 g/cm³ respectively.