Provenance of sediments from Barra del Tordo and Tesoro beaches,Tamaulipas State,northwestern Gulf of Mexico

Mayla A.Ramos-Vázquez and John S.Armstrong-Altrin

Abstract

Keywords:Geochemistry,Sedimentary provenance,Zircon U–Pb,Beach sediment,Rare earth elements,Mexico

1 Introduction

The geochemistry of beach sediments in the Gulf of Mexico have been studied by various authors(Rosales-Hoz et al.2015;Machain-Castillo et al.2019;Ayala-Pérez et al.2021;Kasper-Zubillaga et al.2021),whereas the U–Pb ages of detrital zircons from the Gulf of Mexico coastal sedimentshave not been studied extensively.Nevertheless,U–Pb ages of zircon grains in volcanic rocks from various terranes in Mexico such as Cuicateco,Oaxacan Complex,Xolapa Complex,Mesa Central Province,and Eastern Alkaline Province were widely studied(Solari et al.2004;Rubio-Cisneros and Lawton 2011;Lawton and Molina-Garza 2014;Juárez-Arriaga et al.2019;Sieck et al.2019;Barboza-Gudi?o et al.2020;Torres-Sánchez et al.2020;Vermaet al.2021).

The mineralogical and geochemical compositions of detrital sediments are commonly used to investigate the provenance,weathering,and tectonic setting of source area(e.g.Cullers et al.1988;Chaudhuri et al.2018,2020;Critelli 2018;Al-Kaaby and Albadran 2020;Banerjee et al.2020;Rivera-Gómez et al.2020;Critelli et al.2021;Karlik et al.2021).Although detrital sediments are influenced by diagenesis and recycling,their mineral and geochemistry composition is primarily depended on source rocks.The immobile trace elements like REE,Cr,Hf,Nb,Ni,Th,and V are highly reliable indicators to differentiate sediments derived from felsic and/or mafic igneous rocks(Critelli et al.1997,2003;Zeng et al.2019;Arribas et al.2000;Cullers 2000).Similarly,chondrite-normalized REE patterns and the Eu anomaly are also extensively utilized in various studies to infer source rock characteristics(e.g.Bankole et al.2020;Chaudhuri et al.2020).

Detrital zircon U–Pb age has become a common method to investigate zircon origin and to further discriminate source terranes(e.g.Lee et al.2015;McRivette et al.2019;Wang et al.2020;Armstrong-Altrin et al.2021).Zircon is a common mineral in acidic rocks,and due to its resistance capacity against weathering and thermal alteration, zircon retains U–Pb isotopic signatures related to their origin and parent rocks(Potter-McIntyre et al.2018;Dew et al.2019).Numerous studies proved that detrital zircon U–Pb geochronology is a powerful tool to investigate sediment provenance and their transport pathway(e.g.G?rtner et al.2017;Hoskin and Ireland 2000;Al-Juboury et al.2020).Besides,a few studies also documented the importance of the combination of the zircon U–Pb age and sediment geochemistry to infer provenance and to locate source terranes(Turzewski et al.2020;Zeng et al.2020).

In this study,the mineralogy,bulk sediment composition,and U–Pb ages of detrital zircons retrieved from the Barra del Tordo and Tesoro beaches,northwestern Gulf of Mexico are analyzed for the aims(1)to investigate the sediment provenance and(2)to identify the source terranes supplying sediments to the beach areas.

2 Study area and geology

Forty sediment samples(3kg each)were collected in the Barra del Tordo(23°07′30.78″N–97°45′53.04″W)and Tesoro(22°29′42.39″N–97°51′7.99″W)beaches,located at the Tamaulipas State,the northwestern part of the Gulf of Mexico(Fig.1).A 1-cm-thick top layer of the beach sand was removed before collecting the samples.Approximately a 100 m interval was maintained between sample locations.Metamorphic rocks from the Precambrian to the Paleozoic,marine sedimentary rocks from the Paleozoic to the Quaternary,igneous rocks of the Cenozoic age,and the unconsolidated Quaternary sediments are outcropped in the Tamaulipas State(Demant and Robin 1975;Rubio-Cisneros and Lawton 2011).The sediments in the coastal plain in the northwestern Gulf of Mexico mainly consist of(1)Mesozoic clastic(sandstone,siltstone,and shale)and calcareous(limestone)sedimentary rocks,and(2)Cenozoic volcanic rocks,which belong to the Sierra Madre Oriental Province(Demant and Robin 1975;Hudson 2003).The Pánuco,Tamesí,and Soto la Marina River drainage basins consist of arid to semi-arid Mesa Central Province,Oaxacan Complex,the north-south trending Sierra Madre Oriental,and the Gulf Coastal Plain.

In summer,longshore surface currents have a northward flow direction with an average velocity of 4.5 cm/s.In winter,longshore surface currents flow towards the south with an average velocity of 6 cm/s.Higher wind velocities were recorded during summer,which vary from 3.0 m/s to 5.4 m/s(Yá?ez-Arancibia et al.2009).Monreal-Gómez et al.(1992)documented that the water circulation and hydrodynamic condition of the Gulf of Mexico are controlled by loop currents and anticyclonic rings.There are“northers”associated with weather conditions of a short time scale(2–3 days)with highpressure system that originates frequently in the Rocky Mountains of the United States and travels from high latitudes to the Gulf of Mexico.

3 Materials and methods

3.1 Granulometry

Twenty sediment samples from the Barra del Tordo beach and 20 from the Tesoro beach were air-dried and sieved by a Ro-Tap sieve Shaker using American Standard Test Materials sieves from+7 to+270 mesh sizes with 0.50?interval for 40min(?being a grade-scale that is inverse of the grain diameter;Wentworth 1922).Cumulative curves were constructed and used to calculate the statistical grain-size parameters(mean size and sorting values)by applying the equations of Folk and Ward(1957).The Ro-Tap sieve Shaker is located at the Sedimentology Laboratory,Instituto de Ciencias del Mar y Limnología(ICML),Universidad Nacional Autónoma de México(UNAM),Mexico City.

3.2 SEM–EDS

Sand grain composition was obtained randomly by wave dispersive analyses(WDS)using a JEOL JXA-8900R electron microprobe housed in the Institute of Geophysics,UNAM,Mexico City.Peak counting times were 40s for each element,except for Na and K with 10s.

3.3 Thin-section petrography

Twenty thin-sections(10 from Tesoro beach and 10 from Tordo beach)were prepared following a common procedure.The modal component was studied based on the Gazzi-Dickinson point counting methodology(Gazzi 1966;Dickinson 1970).Combined organic and inorganic stains specific for iron-rich calcite(Katz and Friedman 1965)were adopted to identify the mineralogical variations.

3.4 Sediment composition

Fig.1 Map showing beach location and simplified geology of the Gulf of Mexico coastal areas(Source:Carta Geologica;Scale:1:1000000;Servicio Geológico Mexicano 2008).SMOr:Sierra Madre Oriental.Y:Mesoproterozoic;Pz:Paleozoic;C:Carboniferous;PzCP:Carboniferous–Permian;P:Permian;Mz:Mesozoic;TR:Triassic;J:Jurassic;K:Cretaceous;JK:Jurassic–Cretaceous;KPE:Cretaceous–Paleogene;T:Tertiary;PE:Paleogene;Tpa:Paleocene;Te:Eocene;PEeo:Eocene–Oligocene;To:Oligocene;N:Neogene;M:Miocene;Nmp:Neogene–Miocene–Pliocene;Tp:Pliocene;PQ:Pliocene–Quaternary;Q:Quaternary

Twenty bulk sediment samples(10 from Tesoro beach and 10 from Tordo beach)were powdered by an agate mortar,and the major element concentrationswere determined using a Thermo Scientific Niton FXL 950 X-Ray Fluorescence(XRF)analyzer(Balaram 2021).Accuracy of the major element analysis was monitored by an International Standard Gabbro JGB1(GSJ).The precision of major element data was better than 5%.Loss on ignition was obtained by weighing after combusting 1 h at 1000°C.Similarly,the trace and REE concentrations for 20 sediment samples were determined by a VG Elemental PQII Plus ICP–MSand the operation procedure was similar as the detailed in Jarvis(1988).For data calibration,the United States Geological Survey Standard BCR-2(Basalt,Columbia River)was used.In general,the analytical precision was less than 5%.Eu and Ce anomalies are calculated as Eu/Eu*=EuCN/[(SmCN)·(GdCN)]1/2and Ce/Ce*=CeCN/[(LaCN)·(PrCN)]1/2,respectively (CNmeans chondritenormalized valuesfrom Taylor and McLennan 1985).

3.5 Detrital zircon U–Pb dating

Four hundred zircon grains were analyzed for U–Pb ages,200 from Tordo beach(sample numbers PBT1 and PBT13)and 200 from Tesoro beach(sample numbers PT1 and PT19).The detrital zircon grains were handpicked under a binocular microscope and then mounted in an epoxy disc.The cathodoluminescence images of detrital zircon grains from the Barra del Tordo and Tesoro beach sediments analyzed in this study are provided in Additional file 1(Supplementary Information 1 and 2,respectively).Zircon U–Pb geochronology was performed by a Laser Ablation Inductively Coupled Plasma Mass Spectrometry(LA–ICPMS)at the Laboratorio de Estudios Isotópicos(LEI),Centro de Geociencias,UNAM.

The laser ablation system at LEI consists of a Resonetics M050 workstation,equipped with an LPX 220 excimer laser and an S-155 two-volume cell(Müller et al.2009;Paton et al.2010;Petrus and Kamber 2012).The laser workstation is connected to either a Thermo ICap Qcquadrupole ICPMS or a Thermo Neptune Plus multi-collector ICPMS.A“squid”signal homogenizer is used right after the ablation cell before the ablated material enters the plasma.350ml of He is used as carrier gas,mixed downstream with 4.5ml of N2.A frequency of 5Hz was employed,with a fluence of 6J/cm2.

The background average was subtracted from the net intensity measured for each isotope.The acquisition involved an alternation of 2 analyses of Ple?ovice reference zircon(～337Ma,Sláma et al.2008),NIST 610 standard glasses and 5 unknown zircons,using standardunknown bracketing method to allow downhole fractionation corrections to be performed with an inhouse developed software(Solari et al.2010).Precision on measured ratios of207Pb/206Pb,206Pb/238U, and208Pb/232Th was ～0.7% with 1σ relative standard deviation.Replicate analyses of the Ple?ovice zircon indicate an external reproducibility of 0.8%,0.7%,and 1.6%on the measured207Pb/206Pb,206Pb/238U and208Pb/232Th ratios,respectively.These errors are quadratically included in the quoted uncertainties for individual analyses of the analyzed zircons.Concordia,probability density plots and cumulative proportion curves were generated from the corrected207Pb/206Pb ratios and the integrated concordant and discordant ages through ISOPLOT 3.70 Software(Ludwig 2003).The Tuff-Zirc algorithm combined in the same software was used to calculate the mean206Pb/238U ages and their errors,as well as to filter outliers,which are preferred for grains younger than 1000Ma(Ludwig 1998;Ludwig 2003).For grains with ages of<1000Ma,the analysis was considered concordant if the206Pb/238U and207Pb/235U ages differed by less than 10%.For the grains with ages>1000Ma,the same test wascarried out considering206Pb/238U and207Pb/206Pb ages.

4 Results

4.1 Bulk sediment

4.1.1 Textural parameters

The textural parameters for the Barra del Tordo and Tesoro beach sediments are listed in Additional file 2(Supplementary Information 3).The Tordo beach sediments are predominately medium-grained,moderately sorted,near symmetrical,and very leptokurtic.On the other hand,Tesoro beach sediments are predominantly fine-grained,moderately sorted,near symmetrical,and very leptokurtic.

4.1.2 Mineralogy

The minerals detected by SEM–EDS are mainly quartz,ilmenite,magnetite,titanite,zircon,and anorthite in both the Tordo and Tesoro beach sediments(Fig.2).In addition,the thin section study reveals some accessory minerals,such as amphibole,biotite,chlorite,epidote,feldspar,plagioclase,pyroxene,monocrystalline quartz,lithic volcanic,lithic sediments,and bioclast in the Tordo beach sediments(Fig.3a–f),and biotite,chlorite,plagioclase,and lithic volcanic in the Tesoro beach sediments(Fig.3g–l).

4.1.3 Geochemistry

Fig.2 SEM–EDSspectrum for the a–c Barra del Tordo and d–f Tesoro beach sediments,northwestern Gulf of Mexico.a Quartz spectrum;b Ilmenite spectrum;c Magnetite spectrum;d Titanite spectrum;e Zircon spectrum;and f Anorthite spectrum

The major element concentrations of Tordo and Tesoro beach sediments are reported in the Supplementary Information 4 and their UCC-normalized patterns are showed in Fig.4a.SiO2content is>73wt.%in Tordo and Tesoro beach sediments,which varies between 73.2–91.4 wt.% and 82.7–84.1wt.% respectively.The SiO2content shows negative correlation with TiO2,Al2O3,Fe2O3,MgO,MnO,and CaO for Tordo beach sediments(r=?0.99,?0.19,?0.99,?0.98,?0.99,and?0.63,respectively;n=10) and for Tesoro beach sediments(r=?0.19,?0.17,?0.06,?0.37,?0.43,and?0.52,respectively;n=10).Significantly positive correlations between TiO2vs.P2O5(r=0.74,n=10)and Al2O3vs.K2O(r=0.72,n=10)are observed for the Tesoro beach sediments.The differences in major element concentrations between the Tordo and Tesoro beach sediments are probably related with grain size variations(Mean grain size=2.26±0.29 and 3.11±0.11,respectively;Supplementary Information 3).In Tordo beach sediments,compared to the characteristics of UCC,an enrichment of SiO2content is observed,whereas Al2O3,MgO,CaO,Na2O,K2O and P2O5are depleted(Fig.4a).However,in Tesoro beach sediments,UCCnormalized major element patterns are basically different,except for the enrichment in SiO2content relative to UCC(Fig.4a).

The trace element concentrations of the Tordo and Tesoro beach sediments are reported in Additional file 2(Supplementary Information 5). On the UCCnormalized trace element plot(Fig.4b),the trace element concentrations of both beach sediments are depleted relative to the UCC except Cu content.Correlations between SiO2and V,Sr,Rb,Th,Sc,Zr,and REE contents are not significant in both Tordo(r=?0.99,?0.16,0.25,?0.96,?0.99,?0.99,and?0.97,respectively;n=10)and Tesoro beach sediments(r=0.13,?0.33,?0.09,?0.23,?0.38,0.098,and?0.03,respectively;n=10).While Al2O3is negatively correlated with Co,Ni,and Zr in the samples from Tordo(r=?0.218,?0.22,and?0.21,respectively)and Tesoro(r=?0.49,?0.31,and?0.58,respectively)beaches.Positive correlations can be observed between Cu and Ca,Sr,Rb,Ba and Zn(r=0.63,0.68,0.70,0.71,and 0.91,n=20)in the samples from the Tordo and Tesoro beaches.

The REEcontents in the Tordo and Tesoro beach sediments are reported in Additional file 2(Supplementary Information 5).The chondrite-normalized REE patterns of bulk sediments are shown in Fig.4c,and are more fractionated in the Tordo relative to the Tesoro beach sediments.The REE patterns are with distinct negative europium(Eu/Eu*)anomaly.The Eu/Eu*ratios for Tordo and Tesoro beach sediments range between～0.50–0.80 and～0.56–0.68 respectively.A significantly positive correlation is observed forΣREE against TiO2,Fe2O3,MnO,MgO,Th,and Zr in Tordo(r=0.96,0.98,0.98,0.98,0.99,and 0.99,respectively)and Tesoro(r=0.52,0.52,0.60,0.59,0.51,and 0.72,respectively)beach sediments.

Fig.3 Thin-section photomicrographs for the a–f Bara del Tordo and g–l Tesoro beach sediments,northwestern Gulf of Mexico.Am:Amphibole;Bio:Biotite;Cl:Chlorite;Ep:Epidote;Fsp:Feldspar;Ilm:Ilmenite;Ls:Lithic sediments;Lv:Lithic volcanic;Mag:Magnetite;Pl:Plagioclase;Px:Pyroxene;Qm:Monocrystalline quartz;Qz:Quartz;Zrn:Zircon

4.2 Detrital zircon U–Pb geochronology

4.2.1 Barra del Tordo beach

The zircon U–Pb ages analyzed for Tordo beach sediments(PBT1 and PBT13)are reported in Additional file 2(Supplementary Information 6)and their average age groups are listed in Table 1.Among 200 analyses,193 zircon grains yielded concordant ages(Fig.5a and c).The youngest concordant zircon age is 1.2±0.1Ma,and the oldest is 2094±45 Ma.The spectrums of U–Pb ages of samples PBT1 and PBT13 are similar including age populations of the Proterozoic(2094–595.8 Ma and 1732–632.2 Ma,number of zircons=34 and n=21,respectively),the Paleozoic(455.7–262 Ma and 481.2–252.5Ma,n=5 and n=8,respectively),the Mesozoic(177.7–70.3 Ma and 243.8–68.9 Ma,n=24 and n=28,respectively),and,the Cenozoic(59.5–1.2 Ma and 58.2–1.2 Ma,n=33 and n=40,respectively).The probability density plots show three major peaks(Fig.5b and d),which correspond to the Proterozoic,Mesozoic,and Cenozoic.

Fig.4 a Upper continental crust(UCC)-normalized diagram for major element concentrations of the Barra del Tordo and Tesoro beach sediments;b UCC-normalized diagram for trace element concentrations of the Barra del Tordo and Tesoro beach sediments;c Chondrite-normalized REE patterns for the Barra del Tordo and Tesoro beach sediments.The UCC and chondrite values are from Taylor and McLennan(1985)

4.2.2 Tesoro beach

The zircon U–Pb ages obtained from two samples(PT1 and PT19)are listed in Additional file 2(Supplementary Information 7).Among 200 analyses,188 zircons yielded concordant ages(Fig.6a and c).The probabilitydensity plots for samples PT1 and PT19 are shown in Fig.6b and d.The youngest concordant zircon age is 23.3±0.8 Ma,and the oldest is 2593±21Ma.Four zircon U–Pb age populations obtained from samples PT1 and PT19 are the Proterozoic(2593–560.4 Ma and 2039–669 Ma;n=25 and n=30,respectively),the Paleozoic(540–254.9Ma and 477.5–268.5 Ma,n=10 and n=10,respectively),the Mesozoic(249.3–78.4 Ma and 244–71Ma;n=30 and n=28,respectively),and the Cenozoic(64.3–25.7 Ma and 65.9–23.3Ma;n=27 and n=28,respectively).

Table 1 Detrital zircon U–Pb age groups of the Barra del Tordo and Tesoro beach sediments,northwestern Gulf of Mexico(refer Additional file 2 for complete data)

5 Discussion

5.1 Elemental variations

A negative correlation of SiO2against major(TiO2,Al2O3,Fe2O3,MgO,MnO,CaO,and P2O5)and trace(V,Sr,Rb,Th,Sc,Zr,and?REE)elements in the Barra del Tordo beach sediments indicates that silica is hosted in quartz rather than accessory minerals.Similarly,a negative correlation between Al2O3and Co,Ni,Zr for Tordo and Tesoro beach sediments indicates an association of Al2O3with accessory minerals rather than aluminosilicates.A positive correlation between TiO2and P2O5in Tordo and Tesoro beach sediments indicates a similar source and their association with accessory minerals.An enrichment in TiO2,Fe2O3,and MnO contents is noted in samples PBT1 and PBT7,indicating their association with magnetite.Similarly,a positive correlation ofΣREE versus TiO2,Fe2O3,MnO,MgO,Th,and Zr in the Tordo and Tesoro beach sediments is suggestive of the association of REE with accessory minerals.Cu reveals a positive correlation against Ca,Sr,Rb,Ba,and Zn for Tordo and Tesoro beach sediments and indicates the association of these elements with a similar source.The REE fractionation of Tordo beach sediments is probably due to the effect of weathering.

Fig.5 U–Pb concordia diagrams and histograms showing the relative age probability distribution of detrital zircon grains from samples PBT1(a and b)and PBT13(c and d)of the Barra del Tordo beach sediments,northwestern Gulf of Mexico.n=Number of zircon grains

5.2 Sediment weathering and sorting

The chemical composition of detrital sediments is highly useful to interpret sediment recycling and intensity of weathering(Devi et al.2017;Basu 2020;Patra and Shukla 2020;Verlekar and Kotha 2020).In this study,weathering indices like Chemical Index of Alteration(CIA;Nesbitt and Young 1982),Chemical Index of Weathering(CIW;Harnois 1988),and Plagioclase Index of Alteration(PIA;Fedo et al.1995)are applied to infer the intensity of weathering(Supplementary Information 4).The CIA,CIW,and PIA values in the Tordo and Tesoro beach sediments vary between～44.6–57.5,～47.5–74.5, and～41.9–63.8 respectively, indicating moderate to high intensity of weathering in the source area(Supplementary Information 4).

Hydraulic sorting of sediments can be evaluated by the Index of Compositional Variability(ICV),which decreases when weathering increases(Cox et al.1995).The ICV values of Tordo beach sediments vary from 0.95 to 2.41,except for two samples PBT1 and PBT7,which show higher values relative to other samples(9.4 and 5.2,respectively).ICV values of Tesoro beach sediments range between 1.61 and 1.78(Supplementary Information 4).According to Cox et al.(1995),ICV>1 indicates littleweathered detrital mineralslike feldspar,and ICV<1 represents fine-grained sediments like clay.In general,ICV values in Tesoro and Tordo beach sediments are larger than 1,indicating the dominance of detrital minerals.

5.3 Provenance

5.3.1 Geochemistry

Fig.6 U–Pb concordia diagrams and histograms showing the relative age probability distribution of detrital zircon grains from samples PT1(a and b)and PT19(c and d)of the Tesoro beach sediments,northwestern Gulf of Mexico.n=Number of zircon grains

To infer the provenance of the Barra del Tordo and Tesoro beach sediments,we employed a provenance discrimination diagram based on the major oxides(Roser and Korsch 1988),which is frequently used in various studies(e.g.Tawfik et al.2018;Bineli et al.2020;Madhavaraju et al.2020;Mustafa and Tobia 2020).On this plot,the samples are clustered in the quartzose sedimentary provenance field(Fig.7),indicating a recycled quartzose provenance.Besides,the TiO2/Zr ratio in detrital sediments is one of the best indicators on source rock types(Girty et al.1996).TiO2/Zr ratio values are higher in mafic igneous(>195)rocks than in intermediate(～55–195)and felsic(<55)igneous rocks.The TiO2/Zr ratios are very low in the Tordo and Tesoro beach sediments(～0.05–0.22 and～0.04–0.09,respectively),indicating their derivation from felsic source rocks.

Trace element concentrations of sediments,particularly immobile elements like La,Sc,Co,Cr,Th,and REE are helpful to predict the nature of source rocks,due to a compositional difference between silicic and mafic rocks(e.g.Anaya-Gregorio et al.2018;Cullers 2000).Hence,the immobile elements and their elemental ratios are utilized in various studies to differentiate source rocks(e.g.Damian et al.2019;Hossain 2019;Tobia and Shangola 2019;Patra and Shukla 2020;Ekoa Bessa et al.2021).La/Sc,La/Co,Th/Sc,and Cr/Th ratios in the Tordo and Tesoro beach sediments are compared with sediments derived from felsic and mafic sources,as well as with average UCC values(Supplementary Information 8).This comparative study reveals that these sediments were likely derived from felsic source rocks.Furthermore,the relative chondrite-normalized REE pattern and the Eu anomaly have also been used to infer the source of clastic sediments(e.g.Cullers 2000;Hernández-Hinojosa et al.2018;Kettanah et al.2021;Madhavaraju et al.2021).The REE patterns of the Tordo and Tesoro beach sediments are both with a significantly negative Eu anomaly,indicating felsic igneous rocks as their sources(Fig.4c).

5.3.2 Geochronological evidence

Fig.7 Provenance discriminant function(DF)diagram for the Barra del Tordo and Tesoro beach sediments(modified after Roser and Korsch 1988). DF 1=(?1.773·TiO2)+(0.607·Al2O3)+(0.760·Fe2O3)+(?1.500·MgO)+(0.616·CaO)+(0.509·Na2O)+(?1.224·K2O)+(?9.090);DF 2=(0.445·TiO2)+(0.070·Al2O3)+(?0.250·Fe2O3)+(?1.142·MgO)+(0.438·CaO)+(1.475·Na2O)+(1.426·K2O)+(?6.861)

Three major zircon U–Pb age populations are identified in the Barra del Tordo and Tesoro beach sediments(Table 1;Fig.8),i.e.,(1)Proterozoic(2094–595.8Ma and 2593–560.4 Ma respectively);(2)Mesozoic(243.8–68.9 Ma and 249.3–71 Ma respectively); and (3)Cenozoic(59.5–1.2 Ma and 65.9–23.3Ma,respectively).Detrital zircon U–Pb dating ages of this study are compared with the detrital zircon ages from potential source terranes of northwestern Mexico(Fig.8).

The Proterozoic age population in this study is similar to zircon U–Pb ages reported from the Sierra Madre Oriental,Tamaulipas State,which represent the Grenvillian(～1250–920Ma)and the Pan-African orogen period(～730–530Ma).Another likely source for the Proterozoic zircons of the Tordo and Tesoro beaches is the Oaxacan Complex,which contains extensive outcrops of the Grenvillian basement rocks(Barboza-Gudi?o et al.2010;Solari et al.2014).Weber et al.(2010)reported the zircon U–Pb ages from the Oaxacan Complex,which vary from 1500 Ma to 950Ma.Recently,Vega-Granillo et al.(2020)assigned a Proterozoic age for the zircons dated from the northern part of the Oaxacan Complex(～1300–950Ma).Furthermore,Escalona-Alcázar et al.(2016)and Barboza-Gudi?o et al.(2020)reported the Proterozoic age zircons in the Mesa Central Province and they inferred that the Oaxacan Complex is the potential source terrane which supplied Proterozoic zircons to the Mesa Central.Hence,here it is inferred that the potential source terranes,which contributed Proterozoic age zircons to the Tordo and Tesoro beaches,are Sierra Madre Oriental and Oaxacan Complex.Besides,although the numbers of Paleozoic age zircons are both rare in Tordo(n=13)and Tesoro(n=20)beach sediments,we infer that the likely source for theses Paleozoic zircons is the La Joya Formation in the Mesa Central Province(the Huizachal Group;Rubio-Cisneros and Lawton 2011)due to the resemblance of zircon U–Pb ages(～256–222 Ma;Rubio-Cisneros and Lawton 2011).

The Mesozoic age population representing the likely source terrane,which supplied Mesozoic zircons to the beach areas,is the Mesa Central Province.The Mesa Central Province consists of Zacatecas,Nazas,and La Joya sedimentary formations(Zavala-Monsiváis et al.2012).Lawton and Molina-Garza(2014)reported zircon age populations of the Zacatecas and Nazas Formations in northeastern Durango,which vary between～250–179Ma and～280–200Ma,respectively.In addition,Escalona-Alcázar et al.(2016)reported zircon U–Pb ages analyzed from the conglomerate of Zacatecas Formation in the Eastern Alkaline Province,and inferred the maximum depositional age of～81 Ma.Recently,Juárez-Arriaga et al.(2019)dated zircons from different sedimentary formations in the Mesa Central Province and reported the following ages:(1)Tezapotla(84±0.6 Ma),(2)Alamar(79.2±1Ma),(3)Jalpan-Conca(78.7±0.05 Ma),and(4)Agua Zarca(77.9±0.5Ma).Based on the comparison of zircon ages,it is identified that the Mesa Central Province is the source terrane,which supplied Mesozoic zircons to the beach areas.

Fig.8 Histograms for comparing ages of potential source terranes from northwestern Mexico with detrital zircon ages in the present study.a Sierra Madre Oriental;b Mesa Central Province;c?d Tesoro beach(thisstudy);e Oaxacan Complex;f Eastern Alkaline Province;g?h Barra del Tordo beach(this study).Age distribution of zircons from potential source terranes in the northwestern Mexico are after Barboza-Gudi?o et al.(2010),Weber et al.(2010),Rubio-Cisneros and Lawton(2011),Zavala-Monsiváis et al.(2012),Lawton and Molina-Garza(2014),Solari et al.(2014),Escalona-Alcázar et al.(2016),Juárez-Arriagaet al.(2019),Sieck et al.(2019),Barboza-Gudi?o et al.(2020),and Vega-Granillo et al.(2020)

The Cenozoic age population representsthe Eastern Alkaline Province(EAP),which includes various volcanic fields,located in the Gulf of Mexico coastal region(Verma et al.2020).The zircon ages reported from the EAP volcanic fields range from 8.0Ma to 2.82Ma.For instance,zircon agesinferred by K–Ar and Ar–Ar dating are≦9Ma in the Tlanchinol succession(7.33–2.82Ma),Sierra Tantima(～9.0Ma),and Alamo volcanic field(～6.91Ma)of the EAP(Sieck et al.2019).Another source for the Cenozoic age zircons is the Mesa Central Province(MCP),which consists of three tectonomagmatic stages(Sieck et al.2019).First,a long phase of andesitic to dacitic arc volcanism(～45–31Ma)connected with the subduction of the Farallon Plate that lasted until Eocene–Oligocene.Second is an extension-related phase(～32–25 Ma)overlapped with the Oligocene to Miocene basin.The third stage corresponds to an Early Miocene(～23–20 Ma)mafic magmatism that marked the evolution from a continental arc regime to intraplate volcanism.Hence,the similarity in zircon age populations of thisstudy with EAP and MCP reveals that these provinces are the likely sources for the Cenozoic zircons of the Tordo and Tesoro beaches.

Fig.9 Simplified diagram showing river drainage basins and source terranes,which contributed detrital zircons to the Barra del Tordo and Tesoro beaches in the northwestern Gulf of Mexico.The longshore current illustrations are after Monreal-Gómez et al.(1992).Source rocks are inferred through the zircon U?Pb ages compiled from Weber et al.(2010),Zavala-Monsiváis et al.(2012),Lawton and Molina-Garza(2014),Escalona-Alcázar et al.(2016),Silva-Romo et al.(2018),Juárez-Arriaga et al.(2019),Sieck et al.(2019),Barboza-Gudi?o et al.(2020),and Vega-Granillo et al.(2020).SC:San Carlos;ST:Sierra de Tamaulipas;TTA:Tlanchinol?Tantima?Alamo;VH:Valle de Huizachal

Based on the above discussions,it is assumed that the San Fernando,Soto la Marina,Tamesí,and Pánuco Rivers drained along the Sierra Madre Oriental,Mesa Central Province,and Oaxacan Complex carried sediments to the Barra del Tordo and Tesoro beaches and played an important role in defining the sediment source.The probable sediment transport pathway isbriefly illustrated in Fig.9.

6 Conclusions

The mineralogy,sediment composition,and zircon U–Pb ages were analyzed to identify the sources,which supplied sediments to the Barra del Tordo and Tesoro beaches in the northwestern Gulf of Mexico.The results of this study reveal that the Barra del Tordo and Tesoro beach sedimentsaremainly composed of quartz,ilmenite,magnetite,zircon,titanite,and anorthite.The weathering indices like CIA,CIW,and PIA indicate a moderate to high intensity of weathering in the source area.The major element concentrations,trace elemental ratios of La/Sc,La/Co,Th/Sc,and Cr/Th,and REEpatterns with prominent negative Eu anomaly in the Barra del Tordo and Tesoro beach sedimentsall indicate a felsic provenance.

Zircon U–Pb dating for the Barra del Tordo and Tesoro beach sediments suggests three dominant age peaks:Proterozoic(～2593–560.4Ma),Mesozoic(～249.3–68.9Ma),Cenozoic(～69.9–1.2 Ma).Differences in U–Pb age populations between the Barra del Tordo and Tesoro beach sediments are not significant,suggesting a similar provenance.It is inferred that the Proterozoic zircons in the beach areasweretransported from the Grenvillian igneous suites in the Oaxacan Complex as well as from the Sierra Madre Oriental.The potential suppliers for the Mesozoic and Cenozoic age zircons to the beach areas are the Eastern Alkaline and Mesa Central Provinces,among which the Mesa Central Provinceisthe major contributor of sediments to the Barra del Tordo and Tesoro beach areas.Furthermore,the San Fernando,Soto la Marina,Tamesí,and Pánuco Rivers in the northwestern Gulf of Mexico are interpreted to deliver source sediments to the beach areas,which were subsequently mixed by longshore currents.

Abbreviations

∑REE:Total rare earth elements;Am:Amphibole;BCR-2:United States Geological Survey Standard Basalt (Columbia River); Bio: Biotite;C:Carboniferous;CIA:Chemical Index of Alteration;CIW:Chemical Index of Weathering;Cl:Chlorite;CMC:Chiapas Massif Complex;CN:Chondritenormalized;CONACyT:Consejo Nacional de Ciencia y Tecnología;DGAPA:Dirección General de Asuntos del Personal Académico;EAP:Eastern Alkaline Province;Ep:Epidote;Fsp:Feldspar;GSJ:Geological Society of Japan;HREE:Heavy rare earth elements;ICML:Instituto de Ciencias del Mar y Limnología;ICPMS:Inductively Coupled Plasma Mass Spectrometry;ICV: Index of Compositional Variability;Ilm:Ilmenite; J: Jurassic;JK:Cretaceous–Jurassic;JGB1:International Standard Gabbro;K:Cretaceous;KPE:Cretaceous–Paleogene;LA–ICPMS:Laser Ablation Inductively Coupled Plasma Mass Spectrometry;LEI:Laboratorio de Estudios Isotópicos;LREE:Light rare earth elements;Ls:Lithic sedimentary;Lv:Lithic volcanic;M:Miocene;Mag:Magnetite;MCP:Mesa Central Province;Mz:Mesozoic;N:Neogene;Nmp:Neogene–Miocene–Pliocene;P:Permian;PAPIIT:Programa de Apoyo a Proyectos de Investigación e Innovación Tecnología;PCML:Posgrado en Ciencias del Mar y Limnologia;PE:Paleogene;PEeo:Eocene–Oligocene;PIA:Plagioclase Index of Alteration;Pl:Plagioclase;Px:Pyroxene;Pz:Paleozoic;PzCP:Carboniferous–Permian;PQ:Pliocene–Quaternary;Q:Quaternary;Qm:Monocrystalline quartz;Qz:Quartz;R:River;REE:Rare earth element;SC:San Carlos;SEM:Scanning Electron Microscopy;SEM–EDS:Scanning Electron Microscopy–Energy Dispersive X-Ray Spectroscopy;SMOr:Sierra Madre Oriental;ST:Sierra de Tamaulipas;T:Tertiary;Te:Eocene;To:Oligocene;Tp:Pliocene;Tpa:Paleocene;TR:Triassic;TTA: Tlanchinol–Tantima–Alamo; UCC: Upper continental crust;UNAM:Universidad Nacional Autónoma de México;WDS:Wave dispersive analyses;XRF:X-ray fluorescence;Y:Mesoproterozoic;Zrn:Zircon

7 Supplementary Information

The online version contains supplementary material available at https://doi.org/10.1186/s42501-021-00101-4.

Additional file 1:Supplementary Information 1.Post-ablation cathodoluminescence images of selected detrital zircons from the Barra del Tordo beach sediments.Supplementary Information 2.Post-ablation cathodoluminescence images of selected detrital zircons from the Tesoro beach sediments.

Additional file 2:Supplementary Information 3.Graphic mean grain size(Mz)and textural parameters for the Barra del Tordo and Tesoro beach sediments,northwestern Gulf of Mexico.Grain size classifications were based on Wentworth (1922)and Folk and Ward (1957).Supplementary Information 4.Major element concentrations(in wt.%)in the Barra del Tordo and Tesoro beach sediments,northwestern Gulf of Mexico.Supplementary Information 5.Trace and rare earth element concentrations(in ppm)in the Barra del Tordo and Tesoro beach sediments,northwestern Gulf of Mexico.Supplementary Information 6.U–Pb geochronology data for zircons of samples PBT1 and PBT13 analyzed from the Barra del Tordo beach sediments,northwestern Gulf of Mexico.Supplementary Information 7.U–Pb geochronology data for zircons of samples PT1 and PT19 analyzed from the Tesoro beach sediments,northwestern Gulf of Mexico.Supplementary Information 8.Range of trace elemental ratios in the Barra del Tordo and Tesoro beach sediments,compared with ratio range of sediments derived from felsic and mafic rocks,and the Upper Continental Crust.

Acknowledgements

Armstrong-Altrin acknowledges the financial assistance provided by the CONACyTProject(No.A1-S-21287)for zircon grain U–Pb geochronology analysis.We extend our sincere thanks to Patricia Girón García,Carlos Linares-López,Laura E.Gómez Lizárraga and Teodoro Hernández Trevi?o,for XRF,SEM-EDS,SEM,and zircon separation,respectively.We appreciate Dr.Luigi Solari of Centro de Geociencias for U–Pb zircon dating.The textural parameters were analyzed by Eduardo Alfredo Morales,Ricardo Martínez Domínguez,Susana Santiago-Perez,and Arturo Ronquillo Arvizu at the Sedimentology Laboratory,ICML.Mayla A.Ramos-Vázquez is grateful to the Posgrado en Ciencias del Mar y Limnologia(PCML)postgraduate program and to Consejo Nacional de Ciencia y Tecnología(CONACyT)for a doctoral fellowship(No.595593/308610).We are also grateful to the ICMLInstitutional Project(No.616)for providing transportation facilities along the Gulf of Mexico coastal areas.Financial support extended by the Programa de Apoyo a Proyectos de Investigación e Innovación Tecnológica(PAPIIT;Project No.IN107020)is highly acknowledged.We acknowledge the Editor and three anonymous reviewers for their comments,which significantly improved the presentation.

Authors’contributions

Mayla A.Ramos-Vázquez collected sediment samples along the Barra de Tordo and Tesoro beaches,northwestern Gulf of Mexico,and prepared the samples for SEM and zircon dating studies.John S.Armstrong-Altrin participated in sample collection,examined and interpreted the data.Both authors contributed in writing,reviewing,and editing the manuscript.And both authors read and approved the final paper.

Funding

This study was supported by the 1)Consejo Nacional de Ciencia y Tecnologia(CONACyT),Project Number A1-S-21287;2)Direccíon General de Asuntos del Personal Académico(DGAPA)—Programa de Apoyo a Proyectos de Investigación e Innovación Tecnología(PAPIIT),Project Number IN107020;3)Instituto de Ciencias del Mar y Limnologia Institutional Fund Number 616.

Availability of data and materials

All data discussed in this study are available in the current Table,Figures,and Supplementary material of this manuscript(Supplementary Information Additional files 1 and 2).

Declarations

Competing interests

The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.

Author details

1Posgrado en Ciencias del Mar y Limnología,Universidad Nacional Autónoma de México,Ciudad Universitaria,04510 Ciudad de México,Mexico.

2Universidad Nacional Autónoma de México,Instituto de Ciencias del Mar y Limnología,Unidad de Procesos Oceánicos y Costeros,Ciudad Universitaria,04510 Ciudad de México,Mexico.

Received:15 January 2021 Accepted:7 July 2021