There are various geological structures found in northeast Mexico, which show the features of the geological evolution in the region since the late Paleozoic. The Sabinas Basin, an area with significant mineral resources, is one of these structures. Significant accumulations of metallic and non-metallic ores and hydrocarbons are associated with the rocks that fill this basin      .
Several geological, geophysical and geochemical surveys have been carried out in this basin with the main objective of prospecting and exploring its mineral deposits  -  . Despite the number of surveys performed, the origin and evolution conditions of some levels of sedimentary rocks remain unknown.
Several previous studies have used magnetic susceptibility (κ) and natural gamma radioactivity (NGR) to assess impurities in limestone and paleoenvironmental change during the formation of sedimentary rocks  -  .
In this context, the κ of specific rocks is related to the proportion and texture (magnetic grain size and shape) in which magnetic minerals are present, with the most significant being iron oxides and hydroxides, iron sulfides and ferromagnesian silicates (e.g., magnetite, hematite, maghemite, pyrrhotite)  . These minerals may be present as accessory minerals in the sedimentary rocks  . In the case of carbonate rocks, such as limestone, this physical property is related to its impurities     .
NGR is mainly caused by the gamma radiation emitted by both the radioactive isotope K40 and radioactive isotopes of the U (U238 and U235) and Th232 series  . Therefore, the NGR is expressed in concentrations of K (%), U (ppm) and Th (ppm). Most potassium is found in rocks that form silicate minerals such as feldspars and micas, as well as some other common K-bearing minerals such as alunite and glauconite  . Th and U are constituents of accessory minerals, such as zircon, monazite, allanite, xenotime, apatite and sphene, as well as some minerals that have Th and U as a major constituent, such as uraninite (which contains U), thorite and thorianite (which contain Th), and uranothorite and uranothorianite (which contain both elements)  . The concentrations of U, Th and K can be related to the impurities in limestones    . The concentration of Th is an indicator of detrital material in limestone    , while the concentration of U may indicate the presence of in organic matter  .
As carbonate sedimentary rocks predominate in the stratigraphic sequence of the Sabinas Basin, this study presents an analysis of potential magnetic susceptibility (κ) and natural gamma radioactivity (NGR) as indirect indicators of impurities in limestones and possible relationships in the sedimentation environment of this rock sequence of the basin.
2. Geological Setting
Found in northeastern Mexico (Figure 1), the Sabinas Basin is over 5000 meters thick and comprises sedimentary rocks deposited in a marine environment. The basin initially developed on the margin of the North American Craton during the opening of the Gulf of Mexico in the early Mesozoic. The upper Tamaulipas (the Tamaulipas arc and Burro-Peyotes block) and Coahuila blocks are the main
Figure 1. Location of the study area (Sabinas Basin) and the studied sites. Structural configuration and tectonic characteristics of northeastern Mexico. Modified from  .
paleotectonic and paleogeographic limits of the basin (Figure 1), which comprises a depression formed by a series of subsided blocks and some elevations to the Northwest, which separate it from the Chihuahua basin. Small uplifted blocks of basement can be found within the Sabinas basin itself (La Mula and Monclova Island)  .
The Sabinas Basin evolved over several evolutionary phases, the first of which being a long period of Triassic to Lower Jurassic magmatic-arc activity. A rupture of the lithosphere occurred during the second phase, while the rift that occurred with opening of the Gulf of Mexico developed in the third phase  .
Several marine transgressions and regressions occurred in this region during most of the Mesozoic, with these latter processes, linked to the structural characteristics and tectonic evolution of the region, determining sedimentation patterns. Those patterns changed in the late Cretaceous, when the Laramide orogeny (late Cretaceous to Eocene) caused a tectonic uplift that obliterated the structure of the basin.
Three main depositional supersequences are associated with the rift, drift and foreland stages of the evolution of the Sabinas Basin (Figure 2). The first supersequence comprises a depositional cycle, which consists of conglomerates, salt interbedded with fine-grained red and green terrigenous rocks, anhydrite, and carbonates. The stratigraphic position of this first cycle suggests a period of deposition during the Callovian(?)-lower Oxfordian of the upper Jurassic  . The abovementioned group of rocks mostly appeared in the central part of the basin and is representative of a transgressive-regressive cycle  .
The second supersequence represents the drift stage and comprises depositional cycles II, III and IV. Cycle II is corresponding to the upper Oxfordian-Berriasian age, while Cycle III corresponds to the Valanginian-Aptian age, and Cycle IV corresponds to the Aptian-Cenomanian age.
Carbonated units in the base of Cycle II indicate a transgressive section that degrades up to coastal terrigenous facies. Black carbonaceous shales of the La
Figure 2. Generalized stratigraphic column of the Sabinas Basin and surroundings (modified from  ). Continuous solid line indicates the geological formations studied here.
Casita Formation are found above the base, with a basal member representing a maximum flooding surface that covers a broader area than the preceding cycle that occurred in the Kimmeridgian age  . This formation consists of three members in the Sabinas basin, with the lower and upper members composed of shale and the middle member composed of sandstone interbedded with carbonates. The shales, carbonates and sandstones of the La Casita Formation present continuous transgression and indicate a gradual deepening of the basin. During this period, the northwestern limit of the basin was flooded and the La Mula block was surrounded by the sea.
Terrigenous and clay-rich carbonate rocks with varying grain sizes were deposited during the Berriasian age (early Cretaceous). Represented by interbedded shales and carbonates, the Taraises Formation was deposited during this period  .
The subsidence of the basin and the accumulation of sandstones and fine-grained terrigenous rocks continued in Cycle III. During the Hauterivian-Barremian age, a gradual change in sedimentation occurred, while the deposition of dolomite began with an increase in interbedded siltstone and reddish shales towards the emerging continental areas represented by the Coahuila and Tamaulipas blocks. The Cupido Formation, constituting a barrier reef complex that restricted the entry of seawater of normal salinity into the lagoon, was deposited during these geological periods  . This geological formation was deposited as a prograding system during a sea-level highstand.
Later, a new transgressive system began (Cycle IV) in which the basal high-energy carbonates transitioned upward into low-energy carbonates, with these rocks now known as the Cupidito Limestone  .
The continuation of basin subsidence and the relative increasing rise in sea level are represented by the shales and clay-rich limestone of the La Peña Formation   , a geological formation which represents the maximum level of sea flooding that extended over the entire Sabinas Basin and the Tamaulipas block.
The deposition of basinal micritic limestones and shales began in the Albian age and continued until the Cenomanian. An extensive carbonate platform, in which carbonates and evaporites were deposited, prevailed during this period. These rock types are represented by the Aurora and Kiamichi Formations, as well as the Washita Group (Georgetown, Del Río and Buda Formations).
The regressive and prograding sequences (Supersequence 3; Cycle V), represented by the Eagle Ford, Austin, San Miguel, Olmos and Escondido Formations, originated during the late Cretaceous. The Eagle Ford Formation consists of black shales rhythmically interbedded with clayey limestone and sandstone cemented with carbonate, while the Austin Formation is composed of interbedded calcareous shales and light gray argillaceous limestone. Both geological formations present similar facie changes, which indicate a similar environmental deposition. The higher calcium carbonate content in the Austin Formation is the main difference between the two formations and is linked to climate change and/or the decrease in platform depth during the prograde cycle, which developed at a high eustatic level  .
3. Material and Methods
This study of the κ and NGR found in the carbonate sedimentary rocks of the Sabinas Basin was carried out in four stages. In the first stage, the fieldwork was planned selecting 16 outcrops for measurements of both physical properties and the sampling of the studied rocks. Such carbonate rocks (limestones) studied belong to several geological formations deposited during the Cretaceous Period  . The Taraises, Cupido, La Peña, Aurora and Kiamichi Formations belong to early Cretaceous, whereas the Washita Group and Eagle Ford and Austin Formations were deposited during the late Cretaceous (Figure 2).
In the second stage, both physical properties were measured at various outcrops of these rocks. In each of these outcrops, only one of the geological formations studied is shown, in which the measurements were taken in several points dispersed by the outcrop (Figure 3). Initially, at the selected measurement points 5 measurements of κ were performed using a KT-Plus made by Terraplus Inc., after which one NGR measurement was taken using an RS-125 super gamma ray Spectrometer made by Radiation Solutions Inc., which has 4 channels: total gamma radiation (TGR in μR/h); potassium concentration (K in %); uranium concentration (U in ppm); and, thorium concentration (Th in ppm). The measurements of NGR were performed for 120 seconds at each point, with a total of 153 samples taken at each of the sites. The KT-10 Plus and gamma ray Spectrometer were placed directly on the rock, on a flat surface. The area measured with the KT-10 Plus is approximately 7 cm2, according to the diameter of the console. The spectrometer measures an approximate space of 30 × 30 × 30 cm. In this second stage samples are also taken in the outcrops.
During the second stage, the samples were prepared in order to both measure the physical properties and perform the petrographic analysis (analysis of thin sections), while the third stage comprised the statistical processing of the data obtained for both physical properties (range, average and standard deviation) in each of the geological formations. In the final stage, all information obtained from the research (physical properties, field descriptions and petrography) was interpreted.
Although the research was carried out in outcrops in which no contacts between geological formations are observed (only a geological formation is showed) in this last stage, a virtual record of the cross-sectional measurement of the limestone sequence was also undertaken, taking the average κ and NGR values for each geological formation measures in different outcrops and considering their stratigraphic positions (see the stratigraphic column in Figure 2). The constructed stratigraphic column allows using the both physical properties to compare geological formations deposited in different geological time. Initially, an analysis of the possible impurities in limestones of these geological formations is carried out, according to the values of magnetic susceptibility and NGR. Later possible
Figure 3. Outcrops of carbonated rocks. (a) Aurora Fm; (b) Washita Group. Filled black circles indicate measurements points in the outcrops.
relationships with the sedimentation conditions of these rocks are suggested.
In the geological formations studied here, the limestones differ in terms of their κ and NGR values, with the highest κ values recorded in the Washita Group and the Taraise Fm., while the Austin and Aurora Formations recorded the lowest values (Table 1). The limestones of the Washita Group and the Taraise Fm. have the highest values of κ.
The highest NGR values (total gamma radiation) were recorded in the Washita Group, as well as in the Cupido, La Peña and Eagle Ford formations, with the lowest values recorded in the Austin, Aurora and Taraise formations (Table 2). These values are related to Th concentrations.
Table 1. Descriptive statistics of the κ × 10−3 SI in limestones of some Cretaceous formations of the Sabinas basin. N: number of outcrops and samples; SD: standard deviation; O + S: total of outcrops and samples. Av.: Average.
Generally, the limestones of the Lower Cretaceous present a decline in κ and NGR from the lowest to the highest levels of the stratigraphic sequence (Figure 4), behavior which is also observed in the Upper Cretaceous limestones. In the lower parts of both Cretaceous epochs, the limestone has high κ values.
The limestones of the Taraise Formation, located in the lower part of the Lower Cretaceous in Figure 4, present the highest κ values. These rocks have low total gamma radiation values due to their low concentrations of Th and K
Table 2. Descriptive statistics of the NGR in outcrops of some Cretaceous limestones of the Sabinas basin. M: measurement.
Figure 4. Variations of κ and NGR in the Cretaceous limestone sequences of the Sabinas Basin. Average values of κ and NGR for each geological formation are shown. Clastic sequence in the geological formation presented is not considered. Horizontal lines in the points indicate the uncertainty of the standard deviations. Chart is not scaled on the vertical. In the same are not considered the thickness of the geological formations.
and slightly higher concentrations of U.
The Cupido Formation was deposited above these limestones, as characterized by a marked decrease in κ and an increase in NGR, with the latter depending on the Th and K concentrations.
The La Peña, Aurora and Kiamichi Formations overlay the Cupido Formation, occupying the middle and upper part of the Early Cretaceous. The three former geological formations are characterized by very low κ values. The La Peña Formation presents a slight decrease in κ and total gamma radiation (TGR) in the underlying limestones. In this geological formation, the K concentration decreases, whereas the Th and U concentrations increase. The limestones of the Aurora Formation were deposited above the La Peña Formation, rocks which have lower κ and NGR values.
The analyzed levels of limestones belong to the Early Cretaceous, which ended with the deposition of the limestone of the Kiamichi Formation in which the κ and NGR values were higher.
The Late Cretaceous begins with the limestones of the Washita Group, which show an abrupt increase in the κ and NGR values, with the latter due to the K and Th concentrations. In the limestones of the Eagle Ford Formation, deposited above the Washita Group, the κ and NGR values (U and Th concentrations) decrease.
The highest level of the geological section analyzed terminates with limestones of the Austin Formation, in which the lowest values of κ are recorded, as well as a greater decrease in NGR, which is linked to Th and K concentrations.
5.1. κ and NGR as Indirect Indicators of Impurities in Cretaceous Limestones Located in the Sabinas Basin
The highest NGR values related to Th concentrations in the limestones of Washita Group and the Cupido, La Peña and Eagle Ford formations, suggesting more siliciclastic impurities   . The highest κ values recorded in the Washita Group and the Taraise Fm. indicates that the limestones of the Washita Group and the Taraise Fm. have the highest contents of paramagnetic minerals or at least medium to low ferromagnetic mineral content  . Considering the typical composition of sedimentary rocks  and the results of petrographic analyzes in this research (Figure 5 and Figure 6) and previous ones  , these minerals could be mainly iron oxides and hydroxides (e.g. magnetite, ilmenite, hematite and goethite) and iron sulphides (e.g. pyrite and chalcopyrite). The descriptions of the outcrops and the analysis of the thin section indicate some of the above mentioned minerals in the rocks studied in this research (Figure 5 and Figure 6). Other research in this region also reports these magnetic minerals on carbonate rocks  . These minerals were probably formed during the enrichment of organic matter that occurs in most of the limestones in the region   . The accumulation of organic matter could generate iron sulphides, such as pyrite  , a mineral which could subsequently undergo an oxidation process and generate iron oxides, such as magnetite  . Some limestones of the region
Figure 5. Outcrops of limestone in the Washita Group. The arrows indicate nodules of iron oxides.
Figure 6. Thin section of a sample of limestone belonging to the Cupido Formation. Microphotographs taken with reflected light and without an analyzer. Ox: iron oxides (gray color); Pr: pyrite; and, G: transformation of iron oxides to goethite.
do present iron nodules (Figure 5), thus indicating the process described.
In the Taraise Formation the low concentrations of Th and K and slightly higher concentrations of U (Figure 4) can be characteristics of limestone gray and a non-clayey. According to  , limestones with this characteristics belonging to the lower level of this geological formation. In this case, the high values of κ and the high concentrations of U probably linked to the iron sulfides and organic matter, respectively. Previous studies report both components in these rocks  . It is likely that these minerals were acquired during the accumulation of organic matter. The decrease in U concentrations and κ in the Cupido Formation deposited above the Taraise Formation indicates that, in these limestones, κ remains linked to the processes for the enrichment of organic matter, during which authigenic pyrite was formed. Both groups of physical properties suggest less pure limestones, that is, ones with higher clay content (detritus-rich limestones) and less organic matter at this stratigraphic level (Cupido Fm.).
A probable relationship between κ and the K concentration is observed in the La Peña Formation possibly linked to iron-rich and K silicates such as biotite. Increases in U and Th concentrations suggest clay with a low K content, such as kaolinite and montmorrillonite. The limestones of the Aurora Formation deposited above the La Peña Formation, have lower κ values, as well as lower U and Th concentrations indicating decreases in detrital content.
In the Kiamichi Formation, the NGR is linked to the concentrations of U and a decrease in both Th and K concentrations. This result suggests a large increase in organic matter and a decrease in detritus.
In the Washita Group, the abrupt increase in the κ and the K and Th concentrations, as well as the decrease in U concentrations, indicate a probable decrease in organic matter, and suggest that κ is related to the siliciclastic content of the rock.
In the limestones of the Eagle Ford Formation the decrease of the κ values and U and Th concentrations, as well as the small increase in K concentrations reveal a decrease in clay content, as well as increases in nonferrous silicates rich in potassium (e.g. feldspars and muscovite).
The lowest values of κ recorded in the Austin Formation, as well as the decrease in the Th and K concentrations suggests a decrease in the siliciclastic content of these limestones compared to the lower stratigraphic levels of Late Cretaceous. The slight increase in U concentrations indicates a rise in the organic matter content of these limestones, which is higher than in the limestone of the Eagle Ford Formation.
The Th concentration is an indicator of detrital material in limestone  and the U concentration indicates enrichment in organic matter  . Based on the foregoing, it is possible to use the Th/U ratio to ascertain the prevalence of increases in detrital material compared to organic matter in the lower levels of limestones in the Early Cretaceous (Aurora Formation). At the end of Early Cretaceous, the Kiamichi Formation changes the relationship pattern and reveals the accumulation of organic matter. The Late Cretaceous begins with a predominance of the accumulation of detrital material in the Washita Group. Later in the Eagle Ford Formation, the accumulation of organic matter predominates, a trend which is greatly reduced during the deposition of the Austin Formation.
Generally, in the Early Cretaceous, the κ values in the limestones are mainly related to the U concentrations and, therefore, to the content of organic matter, whereas, in the Late Cretaceous, the κ values are linked to detrital material, thus indicating the nature and characteristic of the geological formations studied here.
5.2. Possible Relationships of κ and NGR with Changes in the Sedimentation Conditions at Different Levels of Cretaceous Limestone in the Sabinas Basin
From the description and analysis of the κ and NGR behavior of Cretaceous limestones, it is possible to create a reconstruction of the formation environment of the levels of Cretaceous limestones studied here (Figure 4). Both physical properties suggest that some levels of the Taraise Formation (Early Cretaceous) were deposited in a moderately deep basin environment, dominated by precipitation of CaCO3 as well as the accumulation and transformation of organic matter. Authigenic pyrites are generated during this process, which, themselves, once oxidized, generate iron oxides and hydroxides, some of which are slightly magnetic. Subsequently, a marine transgression occurred  , causing a rise in sea levels. During this period, the limestones of the Cupido Formation were deposited in a prograding depositional environment with rising sea levels  . This environment enabled a high level of CaCO3 production, less accumulation of organic matter and a greater amount of detritus.
Another marine transgression occurred in the mid-Early Cretaceous    , generating a sea level rise and, thus, conditions for increased detritus in the limestones of the La Peña Formation. The sea level subsequently decreased in the basin, with less detritus deposited in the limestone of the Aurora Formation. Subsequently, limestones with a greater amount of organic matter and less detritus were deposited (the Kiamichi Formation).
At the beginning of the Late Cretaceous (Cenomanian-Santonian), the basin depth decreased further, leading to the deposition of a limestone sequence with little organic matter and a large amount of detritus. Both components in this rock decreased as the period progressed. The decrease in basin depth is associated with a strong uprising in either the same or a marine regression. This process caused an increase in the influx of terrigenous material into the basin, generating the limestone of the Washita Group. According to research carried out in other regions  , high increases of κ and NGR can indicate a low sea level in the basin. Subsequently, a slight subsidence began in the basin, causing the accumulation of a lower amount of terrigenous material, mainly comprising clay, and generating the limestones of the Eagle Ford Formation. This process continued, resulting in the deposition of greater amounts of pure limestone, forming the Austin Formation.
The analysis carried out in this study suggests that magnetic susceptibility and natural gamma radioactivity can be considered as indirect indicators of impurities in limestones and possibly the sedimentation environment during the deposition of the carbonated sedimentary rocks of the Sabinas Basin. The κ and NGR values suggest slight differences in the degree of impurities and the mineralogical compositions of limestones deposited in this basin during the Cretaceous. These values also suggest variability in the deposition environments of these rocks. These environments determine different compositions of paramagnetic and ferromagnetic minerals in the limestones.
The reconstruction of the formation environment of the limestone sequence of the Lower and Upper Cretaceous corroborates the potential use of κ and NGR as indirect indicators for the formation conditions of limestones, indicating the enrichment of organic matter and siliciclastic components. This finding indicates that it is possible to use these physical properties along a vertical section of a stratigraphic sequence to determine changes and conditions in the formation environment of the basin.
The authors gratefully acknowledge the University of Coahuila for the support received, and thank the Department of Crystallography Mineralogy and Mineral Deposits at the University of Barcelona for providing access to the Optical Mineralogy Laboratory.
 Eguiluz de Antuñano, S. (2001) Geologic Evolution and Gas Resources of the Sabinas Basin in Northeastern Mexico. In: Bartolini, C., Buffler, R.T. and Cantú-Chapa, A., Eds., The Western Gulf of Mexico Basin: Tectonics, sedimentary basins, and Petroleum Systems, AAPG Memoir, 75, 241-270.
 Corona-Esquivel, R., Tritlla, J., Benavides-Muñoz, M.E., Piedad-Sánchez, N. and Ferrusquía-Villafranca, I. (2006) Geología, estructura y composición de los principales yacimientos de carbón mineral en México [Geology, Structure and Composition of the Main Coal Deposits in Mexico]. Boletin de la Sociedad Geológica Mexicana, 58, 141-160.
 García-Alonso, D., Canet, C., González-Partida, E., Villanueva-Estrada, R.E., Prol-Ledesma, R.M., Alfonso, P., Caballero-Martínez, J.A. and Lozano-SantaCruz, R. (2011) The Cretaceous Sediment-Hosted Copper Deposits of San Marcos (Coahuila, Northeastern Mexico): An Approach to Ore-Forming Processes. Journal of South American Earth Sciences, 31, 432-443.
 González-Sánchez, F., González-Partida, E., Canet, C., Atudorei, V., Alfonso, P., Morales-Puente, P., Cienfuegos-Alvarado, E. and González-Ruiz, L. (2015) Geological Setting and Genesis of Stratabound Barite Deposits at Múzquiz, Coahuila in Northeastern Mexico. Ore Geology Reviews, 81, 1184-1192.
 González-Sánchez, F., Puente-Solís, R., González-Partida, E. and Camprubí, A. (2007) Estratigrafía del Noreste de México y su relación con los yacimientos estratoligados de fluorita, barita, celestina y Zn-Pb [Stratigraphy of Northeastern Mexico and Its Relationship with Stratabound Fluorite, Barite, Celestine and Zn-Pb Deposits]. Boletin de la Sociedad Geológica Mexicana, 59, 43-62.
 Adams, J.A.S. and Weaver, C.E. (1958) Thorium-to-Uranium Ratios as Indicators of Sedimentary Processes; Example of Concept of Geochemical Facies. American. Association of Petroleum Geologist Bulletin, 42, 387-430.
 Bábek, O., Kalvoda, J., Aretz, M., Cossey, P.J., Devuyst, F., Herbig, H. and Sevastopulo, G. (2010) The Correlation Potential of Magnetic Susceptibility and Outcrop Gamma-Ray Logs at Tournaisian-Viséan Boundary Sections in Western Europe. Geologica Belgica, 13, 291-308.
 Da Silva, A., Potma, K., Weissenberger, J., Whalen, M., Humblet, M., Mabille, C. and Boulvain, F. (2009) Magnetic Susceptibility Evolution and Sedimentary Environments on Carbonate Platform Sediments and Atolls, Comparison of the Frasnian from Belgium and Alberta, Canada. Sedimentary Geology, 214, 3-18.
 Ellwood, B.B., Crick, R.E., El Hassani, A., Benoist, S.L. and Young, R.H. (2000) Magnetosusceptibility Event and Cyclostratigraphy Method Applied to Marine Rocks: Detrital Input versus Carbonate Productivity. Geology, 28, 1135-1138.
 Hladil, J., Gersl, M., Strnad, L., Frana, J., Langrova, A. and Spisiak, J. (2006) Stratigraphic Variation of Complex Impurities in Platform Limestones and Possible Significance of Atmospheric Dust: A Study with Emphasis on Gamma-Ray Spectrometry and Magnetic Susceptibility Outcrop Logging (Eifelian-Frasnian, Moravia, Czech Republic). International Journal Earth Sciences, 95, 703-723.
 Doveton, J.H. and Merriam, D.F. (2004) Borehole Petrophysical Chemostratigraphy of Pennsylvanian Black Shales in the Kansas Subsurface. Chemical Geology, 3-2, 249-258.
 Lüning, S., Wendt, J., Belka, Z. and Kaufmann, B. (2004) Temporal-Spatial Reconstruction of the Early Frasnian (Late Devonian) Anoxia in NW Africa; New Field Data from the Ahnet Basin (Algeria). Sedimentary Geology, 163, 237-264.
 Chávez-Cabello, G., Aranda-Gómez, J.J., Molina-Garza, R.S., Cossío-Torres, T., Arvizu-Gutiérrez, I.R. and González-Naranjo, G.A. (2005) La falla San Marcos: una estructura jurásica de basamento multirreactivada del noreste de México. [The San Marcos Fault: A Jurassic Structure of Multireactived Basement in Northeastern Mexico]. Boletin de la Sociedad Geológica Mexicana, 57, 27-52.
 Goldhammer, R.K., Lehmann, P.J., Todd, R.G., Wilson, J.L., Ward, W.C. and Johnson, C.R. (1991) Sequence Stratigraphy and Cyclostratigraphy of the Mesozoic of Sierra Madre Oriental, Northeast Mexico. A Field Guide Book: Gulf Coast Section, Society of Economic Paleontologists and Mineralogists, Tulsa, 85 p.
 Murillo-Muñeton, G. (1999) Stratigraphic Architecture, Platform Evolution, and Mud-Mound Development in the Lower Cupido Formation (Lower Cretaceous), Northeastern Mexico. Ph.D. Dissertation, Texas A&M University, College Station.
 Wilson, J.L. and Pialli, G. (1977) A Lower Cretaceous Shelf Margin in Northern Mexico. In: Cretaceous Carbonates of Texas and Mexico: Applications to Subsurface Exploration, The University of Texas at Austin, Bureau of Economic Geology, Report of Investigations, 89, 286-298.
 Humphrey, W.E. (1949) Geology of Sierra de Los Muertos Area, Mexico (with Descriptions of Aptian Cephalopods from the La Peña Formation). Bulletin of the Geological Society of America, 60, 89-176.
 Batista-Rodríguez, J., Proenza-Fernández, J., Rodríguez-Vega, A., López-Saucedo, F. and Cázares-Carreón, K. (2017) Magnetic Susceptibility and Natural Gamma Radioactivity as Indirect Proxies for Characterization of Sandstones and Limestones of the Sabinas Basin. Geofizika, 34, 19-43.
 Merinero, R., Lunar, R. and Martínez, J. (2010) Carbonatos metanógenos y pirita framboidal autigénica: geomarcadores de la actividad de organismos quimiosintéticos en el golfo de Cádiz. [Methanogenic Carbonates and Autigenic Pyrite Framboid: Geomarkers of the Activity of Chemosynthetic Organisms in the Gulf of Cadiz]. Macla, 12, 29-37.
 Lehmann, C., Osleger, D.A., Montañez, D.A., Sliter, I.P., Arnaud-Vanneau, A. and Banner, J. (1999) Evolution of Cupido and Coahuila Carbonate Platforms, Early Cretaceous, Northeastern Mexico. Bulletin of the Geological Society of America, 111, 1010-1029.