John S.Armstrong-Altrin

Abstract Zircons are abundant in the beach sediments. In this study, surface microtexture, mineralogy, bulk sediment geochemistry, trace element composition and U-Pb isotopic geochronology of detrital zircons collected from the Riachuelos and Palma Sola beach areas, southwestern Gulf of Mexico were performed to infer the sediment provenance and palaeoenvironment. The zircon microtexture was categorized as mechanically- and/or chemicallyinduced features. The weathering index values for the Riachuelos (～72-77) and Palma Sola (～71-74) beach sediments indicated moderate weathering of both of the two source areas. The major and trace element data of bulk sediments suggested passive margin settings for the two areas. The trace elemental ratios and chondritenormalized rare earth element (REE) patterns of bulk sediments revealed that the sediments were likely sourced by felsic and intermediate igneous rocks. And the zircon Th/U ratios (mostly more than 0.2) and zircon REE patterns(with negative Eu and positive Ce anomalies) suggested a magmatic origin for both of the beach sediments from these two areas. Two distinct zircon age peaks respectively belonging to the Paleozoic and the Cenozoic were identified both in the Riachuelos and Palma Sola beach sediments. Zircon geochronology comparison research between the Riachuelos-Palma Sola beach sediments and potential source areas in SW Gulf of Mexico revealed that the source terrane supplied the Paleozoic zircons of this study was identified as the Mesa Central Province(MCP), and the Cenozoic zircons were transported from the nearby Eastern Alkaline Province (EAP). Moreover,although the Precambrian zircons were very few in the studied sediments, their geochronology and geochemistry results still could infer that they were contributed by the source terranes of Grenvillian igneous suites in the Oaxaca and the Chiapas Massif Complexes.

Keywords: Detrital zircon, Beach sediment,U-Pb dating, Zircon grain morphology, Microtexture, Mineralogy,Geochemistry, Geochronology, Gulf of Mexico

1 Introduction

Zircon is a common accessory mineral in clastic sediments, which retains U-Pb isotopic signatures during erosion, recycling, and transport; differences in zircon age populations in sediments are linked to the nature and age of source terranes (Ji et al. 2019; Joy et al. 2019;Liu et al. 2019; Chaudhuri et al. 2020). Hence, detrital zircon U-Pb age geochronology is reliable to infer provenance, as well as helpful to discriminate sediments derived from various source terranes and to understand their sediment-transport pathways (An et al. 2016;Madhavaraju et al. 2018; Wang et al. 2018a). Although U-Pb ages of detrital zircons have been used widely in different studies to investigate sediment provenance, few studies recommended that the combination of zircon age data and trace element data is more reliable to better understand the depositional history of a complex depositional basin (Wang et al. 2017; Li et al. 2019).

On the other hand, the importance of mineralogical and geochemical composition of clastic sediments to infer intensity of weathering, provenance, and tectonic setting of the source region is well established in the literature (Basu et al. 2016; Carranza-Edwards et al.2019; Casse et al. 2019; Karudu 2019; Men et al. 2019;Kettanah et al. 2020). Although clastic sediments are affected by factors such as weathering, diagenesis, and mineral fractionation during transport, their mineralogical-geochemical composition is mainly related to the parent rock types and the influence factors of the source(McLennan et al. 1993; Basu 2017; Ramos-Vázquez et al.2017, 2018; Bansal et al. 2018; Ndjigui et al. 2019;Rivera-Gómez et al. 2020). Therefore, the geochemical composition of clastic sediments is also reliable to infer the tectonic setting of a sedimentary basin. The immobile trace elements and their elemental ratios differ widely between felsic and mafic rocks,hence are particularly helpful to discriminate the sediments derived from felsic (high in Zr, Nb, Hf, Y, and Th) or mafic (high in Cr, V, Ni, and Sc) sources (Cullers and Podkovyrov 2002; Tawfik et al. 2017). Similarly, different source rocks can be identified by the value of Eu anomaly, for instance, sediments derived from mafic igneous rocks,especially basalt, exhibit positive Eu anomalies or no Eu anomaly (Basu 2017; Wang et al. 2018b; L wen et al.2020). In addition, the mobile elements like K, Na, and Ca among major elements are useful to infer the weathering intensity and other palaeoenvironmental conditions of the source area (Barros dos Santos et al. 2019).

Textural characteristics and geochemical compositions of beach sediments along the Gulf of Mexico have been studied by some authors (Rosales-Hoz et al. 2008;Tapia-Fernandez et al. 2017; Hernández-Hinojosa et al.2018; Ramos-Vázquez et al. 2018). And many studies have been focussed on the zircon U-Pb ages of volcanic rocks from various terranes in Mexico, like Chiapas,Oaxaca, Zacatecas, and Xolapa, with little consideration to the detrital zircons in beach sediments (Keppie et al.2003; Weber et al. 2009, 2018; Talavera-Mendoza et al.2013; Escalona-Alcázar et al. 2016; Cisneros de León et al. 2017; Pe?a-Alonso et al. 2017, 2018; Ramírez-Pe?a and Chávez-Cabello 2017). Recently, based on the geochemical composition of sandstone and the detrital zircon U-Pb geochronology, Wengler et al. (2019)interpreted the provenance of a sedimentary unit in the Mesa Central region, central Mexico. However, investigations based on the comprehensive research of mineralogy, detrital zircon U-Pb geochronology, major and trace element geochemistry of bulk sediment in the Gulf of Mexico beach area to infer the sediment provenance,are meagre.

The degree of grain roundness and microtextures on quartz grain surfaces are considered as powerful tools to infer sediment provenance and have been applied in various studies to reconstruct palaeoenvironments (Margolis and Krinsley 1974; Mahaney 2002; Madhavaraju et al. 2009; Mahaney et al. 2012; Vos et al. 2014;Kaliń?ka-Narti?a et al. 2018; Chmielowska and Woronko 2019). However, microtextures, especially on zircon grains, are not studied sufficiently, possibly because of the difficulties in separating zircon grains from rocks and sediments. To fulfil this gap, the microtextures of detrital zircons from the Riachuelos and Palma Sola beach sediments are surveyed, and the provenance information based on the types of surface features is briefly discussed.

Fig. 2 Map modified after Centeno-García (2017) and Tapia-Fernandez et al. (2017) showing the major terranes of Mexico. AH=Anegada High;CMC=Chiapas Massif Complex; CP=Chiconquiaco-Palma Sola; CUI=Cuicateco; EAP=Eastern Alkaline Province; GC=Guichicovi Complex; LT=Los Tuxtlas;MCP=Mesa Central Province;TTA=Tlanchinol-Tantima-Alamo

Furthermore, in this study, the mineralogy and geochemistry of bulk sediments, U-Pb and trace element data of detrital zircons collected from the Riachuelos and Palma Sola beach areas are analysed, to investigate the source areas supplying sediments to the beach areas and to infer the palaeoenvironments.

2 Geological setting

The Riachuelos and Palma Sola beach areas are located at the Veracruz State, SW Gulf of Mexico (20°25′16.88″N-96°57′28.00″W and 19°46′24.73″N-96°25′19.30″W,respectively; Fig. 1).

The geology of coastal regions of the Gulf of Mexico is shown in Fig. 1. The outcrops along the Gulf of Mexico coast are mainly composed of the Precambrian-Paleozoic metamorphic rocks comprising schist and gneiss, the Mesozoic-Cenozoic clastic and calcareous sedimentary rocks, the Cenozoic mafic and intermediate volcanic rocks, and the Quaternary volcanic rocks and clastic sediments (Keppie et al. 2003; Kasper-Zubillaga et al. 2019).

The Gulf of Mexico littoral zone belongs to the wavedominated transgressive marginal sea coast with narrow coastal plain (Inman and Nordstrom 1971; Davis 1988;Boyd et al. 1992; Carranza-Edwards 2011). The Riachuelos coast is wider than the Palma Sola coast. The easternmost side of the Trans-Mexican Volcanic Belt(TMVB) ends at the Veracruz State area, which formed a cliff at northern Palma Sola beach (Verma et al. 2016),and conduced longshore surface currents of northward flow direction in summer with an average velocity of about 4.5 cm/s. Higher wind velocities were recorded during summer, varying from 3.0 m/s to 5.4 m/s (Ya?ez-Arancibia and Day Jr. 1982). However, during winter,longshore surface currents flow towards south with an average velocity of 6 cm/s. Monreal-Gómez and Salas de León (1990) documented that Gulf of Mexico water circulation and hydrodynamic conditions are controlled by loop currents and anticyclonic rings. The climate for Gulf of Mexico is considered as sub-humid to humid conditions (Tamayo 1991).

Fig.3 X-Ray diffraction analysis result patterns for the Riachuelos (a,b) and Palma Sola(c, d) beach sediments

The major rivers feeding sediment to the western and southern areas of the Gulf of Mexico are Pánuco,Tecolutla, Jamapa, Papaloapan, Coatzacoalcos, Grijalva and Usumacinta (Fig. 1). River Pánuco is 510 km long with a drainage area of ～98,227 km2, which drains into the Gulf of Mexico at Tampico State (Fig. 1). The drainage area of River Pánuco covers parts of Trans-Mexican Volcanic Belt (TMVB) and the Sierra Madre Oriental(SMOr) and its sediment lithology consists of Upper Jurassic-Lower Cretaceous carbonate rocks, Upper Cretaceous marine flysch-type clastic sediments and Pliocene-Holocene marine clastic sediments. River Tecolutla is one of the major rivers in Mexico, starts from the Xolapa uplift that is the easternmost extension of Mexican Volcanic Belt (MVB; Figs. 1 and 2), and,flows through and covers a large volcanic terrane (Self 1975) of an Neogene to Quaternary geological province extending through Central Mexico for about 1000 km(Verma 2009, 2015). River Jamapa flows for about 368 km from the Pico de Orizaba (19°02′N and 97°16′W;volcano Citlalépetl) eastward emptying into the Gulf of Mexico at Boca del Rio, Veracruz State. The sediment lithology of River Jamapa draining area is dominated by the Upper Cretaceous marine flysch, and the Pliocene-Holocene volcanics and marine clastic sediments (Fig.1). River Papaloapan joins in the Gulf of Mexico at Alvarado, Veracruz State, with an annual water discharge of 39,175,000 m3(Tamayo 1991). The drainage area of River Papaloapan is the second largest hydrological basin (17°-19°N and 95°-97°W) in Mexico,which partly covers the Cuicateco terrane (39,189 km2; Fig. 2).River Coatzacoalcos originates in the Oaxaca State area with a mean annual discharge of 32,732 hm3(Tamayo 1991), and generates a drainage basin of between 17°46′-18°10′N and 92°25′-94°31′W. The sediment lithology of the River Coatzacoalcos draining area consists mainly of intrusive igneous rocks and metasedimentary rocks of the Oaxacan Complex. River Grijalva flows for about 640 km from the Chiapas State area of southeastern Mexico. River Usumacinta, originated from the northwestern region of neighbouring Guatemala, is considered as the longest river of Mexico with ～1100 km length, and ranks second among the freshwater rivers flowing into the Gulf of Mexico (Mu?oz-Salinas et al.2017). The Grijalva and Usumacinta Rivers join upstream in about 50 km before they enter the Gulf of Mexico together (Mu?oz-Salinas et al. 2017).

Fig. 4 Energy dispersive spectrums of a Apatite in Sample RIA1; b Zircon in Sample RIA14; c Magnetite in Sample RIA15 of the Riachuelos beach sediments; and d Titanomagnetite in Sample PAL5;e Ilmenite in Sample PAL13;f Quartz in Sample PAL20 of the Palma Sola beach sediments

3 Material and methods

In total, 35 bulk sediment samples (～2 kg each) were collected from the Riachuelos (number of samples n=15) and Palma Sola (n=20) beach areas, where the waves reach the coast during high tide.

The mineralogy of 10 sediment samples in different size fractions (fine- and medium-grained) was identified using the Shimadzu XRD-6000 diffractometer at the Institute of Geology, Universidad Nacional Autónoma de México (UNAM). X-Ray diffraction (XRD) was operated with an accelerating voltage of 40 kV and a filament current of 30 mA, using CuKα radiation and a graphite monochromator.

Scanning Electron Microscopy (SEM) was used to infer the morphological features and shapes of the zircon grains. Fifteen zircon grains were randomly handpicked from each sediment sample under a stereomicroscope,and then 150 zircon grains were selected for microtexture analysis. Grains coated with gold and palladium were mounted on 3 mm-thick carbon discs and examined by JEOL-JSM-6360LV SEM with general magnifications between ×180 and ×5500 at Instituto de Ciencias del Mar y Limnología,UNAM. The quantitative compositional data of randomly-selected zircon grains (7 grains in each sample) were determined by JEOL-JXA-8900R SEM equipped with Energy Dispersive X-Ray Spectrometry (SEM-EDS) at UNAM.

Fig. 5 Cathodoluminescence images of representative detrital zircon grains collected from the Riachuelos and Palma Sola beach sediments. The yellow circles on the images represent the LA-ICP-MS analytical sites(33 μm)

Major element concentration of total 35 bulk sediment samples collected from the Riachuelos and Palma Sola beach areas were analysed by Rigaku RIX-3000 X-Ray Fluorescence at the Institute of Geology, UNAM. Accuracy of major element analysis was monitored by the standard JGB1 (GSJ). The precision of major element data was better than 5%. Loss on ignition was obtained by weighing after 1 h combustion at 1000°C. Trace element concentration of 30 bulk sediment samples were determined by a VG Elemental PQII plus ICP-MS and the operation procedure was similar as details in Jarvis(1988). The United States Geological Survey Standard BCR-2 (Basalt, Columbia River) was used for trace element data calibration, and the analytical precision was less than 5% in general.

Two bulk sediment samples,one from Riachuelos beach area (RIA4) and another from Palma Sola beach area(PAL20) were selected for detrital zircon U-Pb dating. In total, 263 detrital zircon grains (146 grains from the Riachuelos and 117 grains from the Palma Sola) were handpicked from bulk sediments,mounted in epoxy resins and polished. Cathodoluminescence images were taken using an ELM 3R luminoscope to reveal zircon internal texture.All 263 zircon grains were analysed for obtaining information of trace element (Nb, Hf, Th, U, Pb) and REE geochemistry,and U-Pb isotope geochronology.

Zircon U-Pb dating and trace element concentration analyses were simultaneously conducted using a LAICP-MS coupled with Thermo Xii Series quadrupole mass spectrometry, followed by the methodology described by Solari et al. (2018). ANIST 610 glass standard was used to recalculate the trace element concentration,by normalizing them with29Si. U and Th concentrationsare calculated by employing an external standard zircon as mentioned in Paton et al. (2010).207Pb/206Pb ratios,ages and errors are calculated according to Petrus and Kamber (2012). Individual analyses with ＞10%uncertainty, ＞20% discordance, or ＞5% reverse discordance were not considered for interpretation. The U-Pb concordia-age plots, the probability density distribution histograms were plotted using ISOPLOT software(Ludwig 2008). Further detailed analytical methods followed to measure zircon U-Pb isotopes by LA-ICPMS were the same as descriptions in Solari et al. (2018).

Table 1 Microtextures of mechanical and chemical features identified on the zircon grain surfaces in the Riachuelos and Palma Sola beach sediments

Fig.6 Scanning Electron Microscopy(SEM) images showing surface microtextures of detrital zircons from the Riachuelos beach sediments.a Dual striated zircon (dsz) well abraded at both ends; b Euhedral zircon with one side broken edge (bez); c Enlarged view of the broken edge of grain b showing crescentic gouge (crg); d Grain showing arc-shaped steps (as) and linear steps (ls); e Euhedral zircon with one side broken edge (bez),probably due to collision (bec); f Zircon grain with bulbous edge (ble) and reworked conchoidal fracture (rcf); g Grain showing collision impact point (clp), collision fracture (cf) with straight groove (sgr); h Grain with broken edge due to collision (bec); i Enlarged view of broken edge of grain h showing arc-shaped steps (as); j Grain with cavities (gc), adhered particles (ad), solution and precipitation features (s/p); k K-feldspar (kf)adhering to a zircon grain, also showing collision fracture (cf), and solution and precipitation features (s/p); l Grain with silica pellicle (sp), adhered particle(ad), and solution and precipitation features(s/p)

Comparison research between the detrital zircon UPb dating ages tested in this study with the previously published age data (Cantagrel and Robin 1979; López-Infanzón 1991; Nelson and González-Caver 1992; Ferrari et al. 2005; Weber et al. 2009, 2018; Solari et al. 2011;Ortega-Obregón et al. 2014; Escalona-Alcázar et al.2016; Cisneros de León et al. 2017; Wengler et al. 2019)obtained from the nearby and distal major source terranes along the coastal regions of the Gulf of Mexico(such as Chiapas Massif Complex, Eastern Alkaline Province, Mesa Central Province, and Oaxacan Complex in Fig. 2) was conducted to investigate the provenance of the Riachuelos and Palma Sola beach sediments.

Fig. 7 Scanning Electron Microscopy (SEM) images showing surface microtextures of detrital zircons from the Palma Sola beach sediments. a Dual striated zircon(dsz)well abraded at both ends; b Zircon grain with an abraded edge(abe)and adhered particles(ad); c Zircon grain with an abraded edge (abe) and silica pellicle (sp); d Angular grain with circular solution pits (csp) on fracture plane and silica pellicle (sp); e Collision fracture (cf) with adhered particles (ad); f A sub-rounded zircon grain with adhered particles (ad); g A euhedral zircon grain with collision fracture(cf); h Zircon grain showing delamination (dl), meandering ridge (mr), and adhered particle appears to be silica globule (ads); i Enlarged view of grain h showing delamination (dl); j V-shaped percussion cracks (vs) and arc-shaped steps (as); k Reworked conchoidal fracture (rcf) with linear steps(ls);l Zircon grain with adhered particles(ad),silica globule(sg), merged silica globule(msg),and silica flower(sf),due to precipitation

Table 2 Major element concentrations (in wt.%) of the Riachuelos and Palma Sola beach sediments

Table 3 Trace element concentrations (in ppm)of the Riachuelos and Palma Sola beach sediments

Table 3 Trace element concentrations (in ppm)of the Riachuelos and Palma Sola beach sediments (Continued)

4 Results

4.1 Mineralogy of bulk sediments

X-Ray diffraction analysis reveals that zircon,quartz,magnetite,and ilmenite are the major and chromite,geikielite,spinel, and labradorite are the minor minerals in the Riachuelos(Fig.3a,b)and Palma Sola(Fig.3c,d)beach sediments. Based on the quantitative compositional data(Additional File 1:Supplementary Information 1)obtained by SEM-EDS analysis, heavy minerals such as apatite(Fig.4a),zircon(Fig.4b),and magnetite(Fig.4c)are identified in Riachuelos beach sediments, and titanomagnetite(Fig.4d), ilmenite (Fig. 4e), and quartz (Fig.4f)are identified in Palma Sola beach sediments.

4.2 Zircon CL images

Zircons collected from the Riachuelos and Palma Sola beach sediments are mostly prismatic and showing oscillatory zoning with luminescent overgrowths, which supports for a magmatic origin (Fig. 5). More cathodoluminescence images for the zircon sequences of the Riachuelos (Supplementary Information 8) and Palma Sola (Supplementary Information 9) beach sediments are provided in Additional File 2 to better recognize their grain morphology. In general, the morphology differences between zircons from the Riachuelos and the Palma Sola beach sediments are not significant (Fig. 5).

4.3 Zircon surface microtextures

To infer the sedimentary provenance, transportation medium and depositional history of the beach sediments,microtextures on detrital zircon grains from the Riachuelos and Palma Sola beach sediments are investigated and interpreted based on the traditional classification procedures followed by various researchers (Mahaney 2002;Madhavaraju et al.2009;Mahaney et al.2012;Armstrong-Altrin and Natalhy-Pineda 2014; Vos et al. 2014; Mohammad et al. 2020). Microtextures of mechanical and chemical features identified on the Riachuelos and Palma Sola detrital zircons are listed in Table 1,and their SEM images are shown in Figs.6 and 7,respectively.

The mechanical features identified in the Riachuelos zircons include well-abraded euhedral grains with broken edges (Fig. 6a, b, e, h, i), crescentic gouge(Fig. 6c, d), arc-shaped steps and linear steps (Fig. 6d, i),bulbous edge and reworked conchoidal fracture (Fig. 6f),and, collision fractures with straight groove (Fig. 6g, h, k).Similarly, the microtextures of chemical origin identified in the Riachuelos zircons are cavities, adhered particles(ad), solution and precipitation features (s/p), and silica pellicle(sp)(Fig.6j-l).

The microtextures identified in the Palma Sola zircons are like the Riachuelos zircons which exhibit both mechanical and chemical features. The mechanical features include abraded edges (abe) (Fig. 7a), collision fractures(cf) (Fig. 7e, g), V-shaped percussion cracks (vs) and arcshaped steps (as) (Fig. 7j), meandering ridge (mr) (Fig.7h), and reworked conchoidal fracture with linear steps(Fig. 7k). The chemical features include adhered particles (ad) (Fig. 7b, e, f), silica pellicle (sp) (Fig. 7c, d),circular solution pits (csp) (Fig. 7d), delamination (dl)(Fig. 7h, i), silica globule (sg) (Fig. 7h, l), and silica flower (sf) (Fig. 7l).

4.4 Geochemical composition of bulk sediments

The major element concentrations and their elemental ratios are provided in Table 2 (in wt.%). The Riachuelos beach sediments are enriched in Al2O3(the mean value of 14±1 wt.% with one-standard-deviation), TiO2(mean 0.98±0.18 wt.%), Fe2O3(mean 6.72±0.89 wt.%), and P2O5(mean 0.16±0.09 wt.%) contents and depleted in SiO2content (the mean value of 59.6±2.4 wt.%) relative to the Palma Sola sediments (Table 3). Compared with the average upper continental crust (UCC; Taylor and McLennan 1985), the Riachuelos beach sediments are slightly enriched in TiO2, Fe2O3, MnO, MgO, and CaO contents, whereas K2O, Na2O, and P2O5contents of both beach sediments are depleted (Fig. 8a).

Fig. 8 Upper continental crust (UCC) and chondrite-normalized a major element, b trace element, and c rare earth element distribution patterns for the Riachuelos and Palma Sola beach sediments.Normalization values are from Taylor and McLennan(1985)

The trace element and REE concentrations of the Riachuelos and Palma Sola beach sediments are listed in Tables 3 and 4, respectively. Compared with average UCC values, the Riachuelos beach sediments are moderately enriched in V, Sr, Nb, and Ba contents, whereas in Palma Sola beach sediments, other elements except Cr are moderately to strongly depleted (Fig. 8b). Among high field strength elements (HFSE), Zr and Hf are largely immobile during weathering and metamorphic processes (McLennan et al. 1993; Men et al. 2019). Theconcentrations of HFSE are similar between the Riachuelos (～110-173 ppm) and Palma Sola (～80-178 ppm) sediments. The average concentration of transition trace elements like Sc, V, Co, and Ni are higher in Riachuelos beach sediments (16±5 ppm, 139±22 ppm,18±3 ppm,and 34±4 ppm,respectively)than in the Palma Sola beach sediments (9±3 ppm, 73±16 ppm, 10±2 ppm,and 12±1 ppm, respectively) (Table 3). The ∑REE content is also higher in the Riachuelos beach sediments (～109-116 ppm) than in the Palma Sola beach sediments (～47-68 ppm)(Table 4).The chondrite-normalized(CN)REE patterns for the Riachuelos and Palma Sola beach sediments are similar, enriched in LREE (LaCN/SmCN=～2.78-3.93 and ～2.66-3.95, respectively), depleted in HREE (GdCN/YbCN=～2.0-2.9 and ～1.3-2.1, respectively), with low negative to positive Eu anomalies (Eu/Eu*=～0.90-1.17 and ～0.81-1.31,respectively)(Fig.8c).

Table 4 Rare earth element concentrations (in ppm) of the Riachuelos and Palma Sola beach sediments

Table 4 Rare earth element concentrations (in ppm) of the Riachuelos and Palma Sola beach sediments(Continued)

4.5 Geochemical composition of detrital zircons

The trace (Nb, Hf, Th, and U) and REE concentrations of detrital zircons from the Riachuelos and Palma Sola beach sediments are listed in Additional File 1 (Supplementary Information 2 and 3, respectively).

The variation in Y content between the Riachuelos zircons (～221-9300 ppm, n=146) and the Palma Sola zircons (～266-9400 ppm, n=117) is not significant.However, the variations in Th and U contents of the Palma Sola zircons (～21-2900 ppm and ～28-2290 ppm,respectively) are slightly higher than those of the Riachuelos zircons (～16-2017 ppm and ～31-1451 ppm). The Th/U ratios of detrital zircons are generally used to differentiate igneous and/or metamorphic origins: Th/U ratio for zircons of igneous origin is ＞0.3 and for zircons of metamorphic origin is ＜0.1 (Rubatto 2002). The Th/U ratios of the Riachuelos and Palma Sola detrital zircons, respectively varying between 0.14 and 1.56 (0.67±0.29, n=146) and varying between 0.13 and 2.27 (0.68±0.36, n=117), are both consistent with a magmatic source (Fig. 9).

Variations in ∑REE and HREE contents of the Riachuelos zircons (～215-2275 ppm and ～209-2132 ppm,respectively) and of the Palma Sola zircons (～256-2697 ppm and ～250-2608 ppm, respectively) are large. The REE patterns between the Riachuelos and Palma Sola zircons are similar, with depleted LREE and enriched HREE (GdCN/YbCN=2-3 and 1.3-2.1, respectively), and mostly with pronounced negative Eu and positive Ce anomalies (Figs. 10 and 11, respectively). A few zircons enriched in LREE lack Ce anomalies. Moreover, the zircon REE patterns for the Proterozoic, Paleozoic, Mesozoic, and Cenozoic sediments in the Riachuelos (Fig. 10)and Palma Sola (Fig. 11) beach areas are similar.

Fig. 9 Plot of Th/U ratios versus U-Pb ages for the Riachuelos and Palma Sola detrital zircons.Fields are after Rubatto(2002)

Fig. 10 Rare earth element patterns of the detrital zircons collected from the a Proterozoic(n=3);b Paleozoic(n=30);c Mesozoic(n=9);d Cenozoic(n=87)sediments in the Riachuelos beach area.Chondrite normalization values are from Taylor and McLennan(1985).n=Number of zircon grains

4.6 Detrital zircon U-Pb dating ages

Concordia plots and age probability density histograms of the Riachuelos and Palma Sola detrital zircons are shown in Fig. 12. In total, 146 U-Pb ages for the detrital zircons from the Riachuelos beach sediments are available and 129 involved provide the concordant ages(Fig.12a,b; Additional File 1: Supplementary Information 4). Two significant age populations are identified as: (1) the Paleozoic (～356.5-252.5 Ma; n=30), including four Carboniferous ages and 26 Permian ages;and,(2)the Cenozoic(～33.1-0.4 Ma; n=87), including four Oligocene ages, 34 Miocene ages, 12 Pliocene ages and 37 Pleistocene ages. In addition, nine zircon U-Pb ages of ～194-66 Ma represent the Mesozoic,among which,three are of Early Jurassic ages(～194-179 Ma) and six are of Cretaceous (～141.9-66.1 Ma); and three U-Pb ages of ～1370-903.8 Ma represent the Proterozoic Period.

In total, 117 detrital zircons from the Palma Sola beach sediments are analysed and 102 spots provide concordant ages (Fig. 12c, d; Additional File 1: Supplementary Information 5). The Palma Sola zircons exhibit two prominent age peaks of the Paleozoic (～480-255 Ma; n=15) and the Cenozoic (～64.9-0.4 Ma; n=69),like the Riachuelos zircon ages. Besides, minor age peaks are concentrated in the Proterozoic (～1652-630 Ma;n=9) and the Mesozoic (～185.5-72.4 Ma; n=9)(Fig. 12c, d).

5 Discussion

5.1 Statistical analysis of geochemistry data

5.1.1 Pearson's correlation coefficient

Fig. 11 Rare earth element patterns of the detrital zircons collected from the a Proterozoic (n=9); b Paleozoic (n=15); c Mesozoic (n=9); d Cenozoic (n=69) sediments in the Palma Sola beach area. Chondrite normalization values are from Taylor and McLennan (1985). n=Number of zircon grains

The correlation technique is applied in this study to measure the magnitude and direction of the association of elements between the Riachuelos and Palma Sola beach sediments. The correlation of SiO2versus Al2O3,Fe2O3, MnO, CaO, MgO, and TiO2is statistically not significant for the Riachuelos (r=0.39, -0.83, -0.80, -0.72, -0.73, and-0.62, respectively; n=15) and Palma Sola (r=-0.53, -0.54, -0.38, -0.50, -0.29, and-0.40,respectively; n=20) beach sediments (critical t value for 99% confident level is 0.537; Verma 2005), reflecting the quartz abundance effect. For the Riachuelos beach sediments, SiO2exhibits a significant positive correlation with K2O and Na2O (r=0.66 and 0.63, respectively),whereas this correlation is not significant for the Palma Sola sediments (r=-0.08 and-0.47, respectively),which is probably due to variations in the relative proportions of K-feldspar and plagioclase, respectively(Bhattacharjee et al. 2018; Madhavaraju et al. 2020). This is also revealed in differences in K2O and Na2O contents between the Riachuelos (～1.33-2.26 wt.% and ～2.08-3.17 wt.%, respectively) and Palma Sola (～1.19-1.74 wt.% and ～1.37-2.22 wt.%, respectively) beach sediments(Table 2). A significant positive correlation between Al2O3and K2O for the Riachuelos (r=0.81) and Palma Sola (r=0.87) beach sediments indicates that a considerable proportion of K is associated with fine-grained sediments or phyllosilicates (Etemad-Saeed et al. 2015;Greggio et al. 2018; Udayanapillai et al. 2020). Al2O3against TiO2reveals a negative correlation for the Riachuelos (r=-0.36, n=15) and Palma Sola (r=-0.01,n=20) beach sediments, indicating that TiO2is not associated with aluminous clays or its association with titano-magnetite and ilmenite (Nagarajan et al. 2015,2017; Papadopoulos et al. 2019). Similarly, a weak correlation between Al2O3and P2O5for the Riachuelos (r=-0.70, n=15) and Palma Sola (r=0.13, n=20) beach sediments suggests P2O5content is associated with detrital minerals, probably apatite (Fig. 4a).

The correlation between Al2O3and HFSE like Zr, Y,and Hf is statistically not significant for the Riachuelos(r=0.29, -0.28, and 0.05, respectively; n=15) and Palma Sola (r=0.17, -0.29, and 0.12, respectively; n=15) beach sediments, which indicates that HFSEs are not associated with phyllosilicates but with accessory minerals. Similarly,Sc, V, Co, Cr, and Ni are all negatively correlated with Al2O3for the Riachuelos (r=-0.31, -0.40, -0.60, -0.73,and-0.67, respectively; n=15) and Palma Sola beach sediments (r=-0.51, -0.09, -0.42, -0.45, and-0.19, respectively;n=15),suggesting the association of these trace elements with accessory components.

Fig. 12 U-Pb Concordia age diagrams and age probability density histograms of detrital zircons from the Riachuelos (a, b) and Palma Sola (c, d)beach sediments. n=Number of zircon grains

The correlation of ∑REE with SiO2, Al2O3, K2O, and Na2O is statistically not significant for the Riachuelos(r=-0.73, -0.44, -0.41, and-0.37, respectively) and Palma Sola (r=-0.60, 0.20, -0.07, and 0.15, respectively) beach sediments. In addition, the positive correlation of ∑REE versus Fe2O3, TiO2, P2O5, Cr, Ni, V, Zn,and Y is significant for the Riachuelos sediments (r=0.86, 0.70, 0.53, 0.83, 0.77, 0.80, 0.81, and 0.67, respectively), which indicates the association of REE with accessory heavy minerals. In contrast, the variations in correlation coefficient observed in the Palma Sola beach sediments (r=0.02, -0.14, 0.76, -0.43, 0.07, -0.08,0.54, and 0.72, respectively) suggests that the REE are partially associated with accessory minerals.

5.1.2 Principal component analysis (PCA)

The PCA for the Riachuelos beach sediments demonstrates two major factors (Fig. 13; Additional File 1: Supplementary Information 6). Factor F1 reveals significant negative loading, which represents the association of a few major (TiO2, Fe2O3, MnO, and MgO), trace (V, Y,and Zn), and rare earth elements (La, Ce, Pr, Nd, Sm,Eu,Gd,Tb,Dy,Ho,Er,Tm,Yb,and Lu)related to heavy minerals like magnetite and ilmenite. Factor F2 shows positive loading with elements, which is largely associated with detrital materials. On the other hand, PCA for the Palma Sola beach sediments reveals three factors(Additional File 1: Supplementary Information 7). Factor F1 represents high positive loading with the association of alkali (K and Rb), alkaline (Ba, Mg, Nd,and Ca), and rare earth elements (Pr, Nd, Sm, Gd,Tb, Dy, Ho, Er, Tm, Pb, and Lu), related to the detrital origin. Factor F2 shows a negative loading for elements of Be, Cu, Ga, Li, Sn, Ta, Th, and Zn. Factor F3 represents a negative loading with variables of TiO2, Fe2O3, Co, Cr, Hf, V, and Zr, which are associated with heavy minerals like ilmenite, magnetite,chromite, and zircon. Therefore, the combined results obtained from the Pearson's Correlation Coefficient and the Principal Component Analysis indicate the association of major and trace elements with detrital materials rather than biogenic components.

Fig. 13 Bivariate plots of factor F1-factor F2-factor F3 for the bulk geochemistry data of the a Riachuelos and b Palma Sola beach sediments

5.2 Chemical weathering

Numerous weathering indices like Chemical Index of Alteration (CIA=100×[Al2O3/(Al2O3+CaO+Na2O+K2O)];Nesbitt and Young 1982)and Chemical Index of Weathering(CIW=100×[Al2O3/(Al2O3+CaO+Na2O)];Harnois 1988)have been widely used to infer intensity of weathering (e.g.,Ota et al. 2017; Anaya-Gregorio et al. 2018; Raza and Mondal 2018; Tawfik et al. 2018; Hossain et al. 2018; Kettanah et al. 2020). However, a few authors (Cullers 2000; Buggle et al.2011;Garzanti et al.2015)documented that sediments are affected by various factors like sorting,recycling,and having significant proportions of carbonate or phosphate minerals. Hence, care should be taken when applying CIA and CIW indices to infer the intensity of weathering, because CaO content in sediments decrease the values and may mislead interpretation. These authors further recommended to eliminate the CaO content from the CIA and CIW indices and suggested a simply modified index called CIX (CIX=100×[Al2O3/(Al2O3+Na2O+K2O)]), which eliminated CaO content(Cullers 2000;Buggle et al.2011;Garzanti et al.2014).In general,similar to other indices,the CIX values increase if weathering increases and in highly weathered sediments the values are always more than 80. The calculated CIX values for the Riachuelos and Palma Sola beach sediments vary from 74 to 77 and from 71 to 74, respectively(Table 2),suggesting a moderate intensity of chemical weathering for the source area.

In clastic sediments, Th/U, Rb/Sr, and Ba/Sr ratios can be used as indicators to infer intensity of weathering, because Th, U, Rb, Ba, and Sr are sensitive to chemical weathering and their values are always high in extremely weathered sediments (McLennan et al.1993; Xu et al. 2010). Similarly, Th/U ratios in sediments increase with increasing weathering due to oxidation and loss of uranium. Th/U ratios above four are expected to be indicative of weathering (McLennan et al. 1993). The average Rb/Sr and Ba/Sr ratios in the Riachuelos (0.09±0.01 and 1.17±0.08, respectively) and Palma Sola (0.18±0.02 and 1.79±0.18, respectively) beach sediments are lower than in UCC(0.32 and 1.57, respectively; Taylor and McLennan 1985). Similarly, average Th/U ratio in the Riachuelos and Palma Sola beach sediments are both less than 4(3.0±0.7 and 3.3±0.3, respectively; Table 3). Hence,the Rb/Sr, Ba/Sr, and Th/U ratios in the Riachuelos and Palma Sola beach sediments suggest that the intensity of weathering is moderate in the source area.

5.3 Source rock composition and tectonic setting

Trace element concentrations (e.g., Th, U, Y, Sc, Cr, Co,Ni, V, and REE) in clastic sediments are more reliable indicators of provenance than major elements and are extensively applied in various studies to infer sediment provenance (e.g., Basu 2017; Zhang et al. 2018; Casse et al. 2019). This is because incompatible trace elements are concentrated in sediments derived from felsic igneous rocks, whereas compatible elements are enriched in sediments derived from mafic igneous rocks and their weathered products. Then, elemental ratios such as La/Sc, Th/Co, and Th/Sc are significantly different between felsic and mafic rocks and can be used to infer sediment provenance (Cullers and Podkovyrov 2002). Bivariate(Fig. 14a, b) and ternary diagrams (Fig. 14c, d) show that the Riachuelos and Palma Sola beach sediments are plotted between dacite and andesite average compositions,suggesting the provenance of felsic and intermediate igneous rocks. Less contributions from mafic rocks are also revealed by the lower concentrations of Co, Ni, and V in the Riachuelos and Palma Sola beach sediments relative to in the UCC (Fig. 8b).

Fig. 14 Bivariate and ternary plots for the Riachuelos and Palma Sola beach sediments: a La/Sc versus Co/Th; b Sc versus Th/Sc; c Ni-Th*10-V; d Th-Sc-La. 1 This study; 2 The average compositions of dacite, andesite, and basalt plotted in these diagrams are compiled from Verma (2015),Verma et al. (2016),and Torres-Sánchez et al. (2019)

The Riachuelos and Palma Sola beach sediments are little fractionated (LREE/HREE=～7.16-9.54 and ～6.23-8.56,respectively) and are characterized by both negative and positive Eu anomalies (Eu/Eu*=～0.90-1.17 and ～0.81-1.31, respectively; Table 4), suggesting that the sediments were derived from both felsic and intermediate igneous rocks.Furthermore,the chondrite-normalized REE patterns of the Riachuelos and Palma Sola beach sediments are compared with potential source rocks such as rhyolite,dacite, and basaltic andesite located in the Gulf of Mexico coastal areas (Fig. 15). The similarities in REE patterns between the two beach sediments and source rocks are consistent with the interpretation that the detritus was supplied by felsic and intermediate rocks.

Fig. 15 Chondrite-normalized rare earth element patterns of the Riachuelos and Palma Sola beach sediments compared with the adjacent source rocks in the Gulf of Mexico. 1 This study; 2 The average data of basaltic andesite, dacite, and andesite are compiled from Verma (2015),Verma et al. (2016), and Torres-Sánchez et al. (2019); 3 The Upper Continental Crust(UCC)data are from Taylor and McLennan(1985)

The major and trace element concentrations of detrital sediments have been widely used in various studies to infer the tectonic setting of a sedimentary basin, because detrital sediment composition varies significantly between active and passive margin settings (Bhatia 1983;Roser and Korsch 1986; Saha et al. 2010, 2018; Verma and Armstrong-Altrin 2013; Armstrong-Altrin 2015).Although many studies continuously use the old tectonic discrimination diagrams of Bhatia (1983) and Roser and Korsch (1986), in this study, the recently proposed discriminant function diagrams of Verma and Armstrong-Altrin (2016)are considered(Fig.16).These new statistically discriminant diagrams are efficient to discriminate siliciclastic sediments of island or continental arc, continental rift, and/or collision settings.On these diagrams,the Riachuelos (Fig. 16a, c) and Palma Sola (Fig. 16b, d)beach sediments are divided well within the rift and passive margin fields, suggesting a passive margin setting for the beach areas, which is consistent with the tectonic setting of the Gulf of Mexico (Verma et al. 2016;Centeno-García 2017).

Fig. 16 Tectonic discrimination diagrams for the Riachuelos (blue circles) and Palma Sola (red squares) beach sediments, based on a and b major element and c and d trace element concentrations(after Verma and Armstrong-Altrin 2016). The function DF(AP)M is derived from an equation: DF(AP)M=(3.0005×ilr1TiM)+(-2.8243×ilr2AlM)+(-1.0596×ilr3FeM)+(-0.7056×ilr4MnM)+(-0.3044×ilr5MgM)+(0.6277×ilr6CaM)+(-1.1838×ilr7NaM)+(1.5915×ilr8KM)+(0.1526×ilr9PM)-5.9948.The function DF(AP)MT is derived from an equation: DF(AP)MT=(3.2683×ilr1TiMT)+(5.3873×ilr2AlMT)+(1.5546×ilr3FeMT)+(3.2166×ilr4MnMT)+(4.7542×ilr5MgMT)+(2.0390×ilr6CaMT)+(4.0490×ilr7NaMT)+(3.1505×ilr8KMT)+(2.3688×ilr9PMT)+(2.8354×ilr10CrMT)+(0.9011×ilr11NbMT)+(1.9128×ilr12NiMT)+(2.9094×ilr13VMT)+(4.1507×ilr14YMT)+(3.4871×ilr15ZrMT)-3.2088

5.4 Probable provenance based on zircon surface microtextures

The microtextures identified on zircon grain surfaces(Figs. 6, 7) and their implications on provenance are briefly discussed.

The collision fractures, V-shaped percussion cracks,meandering ridges, arc and linear steps indicate highenergy grain to grain impact and mechanical grinding by littoral transport, as well as high-energy subaqueous deposition such as in river and coastal environments (Finzel 2017; Hossain et al. 2020). However, the mechanical features on the grain surface may likely be associated with the hydrodynamic conditions during littoral transport and wave action in the nearshore zone (Kaliń?ka-Narti?a et al. 2018). The abraded grains with bulbous edges identified in a few zircons show the evidence of high-energy collision through saltation in an aeolian environment (Mahaney 2002; Chmielowska and Woronko 2019). Angular grains with broken edges indicate a short transport in littoral environment as well as a highenergy subaqueous environment, where the breaking of waves is probably the transport agent (Costa et al. 2017).Conchoidal fracture is one of the most dominant microtextures on zircon surfaces identified in both beach areas, which indicate liberation of grains from crystalline rocks (Madhavaraju et al. 2009). In addition, collisions between zircon grains and pebbles may also cause fracturing, as well as delamination (Vos et al. 2014).

The etching process is revealed by chemical features like solution pits and precipitation features associated with the percolation of seawater. The silica globules, silica pellicles, and silica flowers indicate that zircon grains were subjected to precipitation with silica saturated solutions for a considerable time. Vos et al. (2014) documented that when precipitation on grain surface continues, silica globules start to merge to form silica flowers and pellicles.

Overall, the mechanical and chemical features identified on the zircon surfaces of the Riachuelos and Palma Sola beach sediments suggest their transport processes by both littoral and aeolian currents, and a high-energy subaqueous coastal depositional environment.

5.5 Probable provenance based on detrital zircon U-Pb ages

U-Pb dating ages of the detrital zircons of the Riachuelos and Palma Sola beach sediments are compared with other zircon U-Pb ages reported by various authors from different terranes of Gulf of Mexico (Fig. 17), and are briefly discussed, to infer the provenance.

5.5.1 Relative zircon age probability distribution—Proterozoic

The Riachuelos beach sediments own three Proterozoic zircons with U-Pb dating ages varying from 1370 Ma to 903.8 Ma and the Palma Sola beach sediments own nine Proterozoic zircons with U-Pb ages varying from 1652 Ma to 629.9 Ma.The possible source terranes contributed to Proterozoic zircons of the Riachuelos and Palma Sola beach areas are the Chiapas Massif Complex and the Oaxaca Complex, where the documented zircon U-Pb ages vary from ～1500 Ma to ～450 Ma (Keppie et al. 2003;Weber et al. 2012). In addition, Escalona-Alcázar et al.(2016)and Wengler et al.(2019)also reported Proterozoic zircons in the Mesa Central Province, and they inferred that Oaxaca Complex and Chiapas Massif Complex are potential source terranes for the Mesa Central Province.

5.5.2 Relative zircon age probability distribution—Paleozoic

Fig. 17 Relative probability distribution of detrital zircon U-Pb ages of the Riachuelos and Palma Sola beach sediments in this study and from four terranes in the SW Gulf of Mexico (referring Figs. 1 and 2 for locations). a Chiapas Massif Complex(Weber et al. 2009, 2018; Cisneros de León et al. 2017); b Eastern Alkaline Province (Cantagrel and Robin 1979; López-Infanzón 1991; Nelson and González-Caver 1992; Ferrari et al. 2005;Solari et al. 2011); c Riachuelos beach area (this study); d Oaxacan Complex (Ortega-Obregón et al. 2014); e Mesa Central Province (Escalona-Alcázar et al.2016;Wengler et al.2019);f Palma Sola beach area(this study).n=Numbers of zircons

Detrital zircons of the Riachuelos and Palma Sola beach areas also show the Paleozoic ages (～480-252.5 Ma).However, Riachuelos beach sediments are more dominated by the Paleozoic zircons (n=30) than Palma Sola beach sediments (n=15). The most likely potential source terrane for the Paleozoic zircons of beach areas,is the Mesa Central Province, located near to Zacatecas City, Mexico. The Mesa Central Province consists of Zacatecas, Nazas, La Boca, and La Joya sedimentary formations (Barboza-Gudi?o et al. 2010). Wengler et al.(2019) reported that the zircon populations from the Zacatecas and Naza Formations vary between ～350-215 Ma and ～440-200 Ma, respectively. These ages are similar to the zircon ages inferred from the Riachuelos and Palma Sola beach areas. Then it can be assumed that the Mesa Central Province is also one of the source terranes, which supplied sediments to the beach areas.The Zacatecas and Nazas Formations in the Mesa Central Province are composed of marine-continental siliciclastic sediments with thin-bedded turbidites, which consist of sandstone, shale and conglomerates (Aguillón-Robles et al. 2014; Centeno-García 2017; Ramos-Vázquez and Armstrong-Altrin 2019; Torres-Sánchez et al. 2019; Verma et al. 2019).

5.5.3 Relative zircon age probability distribution—Mesozoic

Mesozoic zircons are identified in the Riachuelos (～194-66.1 Ma, n=9) and Palma Sola (～185.5-72.4 Ma, n=9)beach sediments. The probable source of Mesozoic zircons in Riachuelos and Palma Sola beach areas is the Chiapas Massif Complex,because similar Mesozoic zircon ages from the Chiapas Massif Complex are also reported by previous studies(e.g.,Weber et al.2009;Cisneros de León et al.2017).

5.5.4 Relative zircon age probability distribution—Cenozoic

Detrital zircon U-Pb dating results (Additional File 1:Supplementary Information 4 and 5) reveal that the Riachuelos and Palma Sola beach sediments are dominated by the Cenozoic zircons (～33.1-0.39 Ma and ～64.9-0.4 Ma, respectively; Figs. 12, 17). The potential source for these Cenozoic zircons is the Eastern Alkaline Province,located along the coastal region of the SW Gulf of Mexico and formed along the Gulf-parallel extensional fault (Robin and Tournon 1978; Robin 1982). The Eastern Alkaline Province consists of alkaline volcanics dominated by basaltic andesites. The detrital zircon ages reported for the Eastern Alkaline Province vary from 65 Ma to 0.18 Ma (Cantagrel and Robin 1979; López-Infanzón 1991; Nelson and González-Caver 1992; Ferrari et al. 2005). Thus the Eastern Alkaline Province is supposed as the provenance for the Cenozoic zircons of the Riachuelos and Palma Sola beach sediments.

6 Conclusions

Zircon, apatite, magnetite, titano-magnetite, and ilmenite were identified in the Riachuelos and Palma Sola beach sediments. The zircon grain morphology and microtextures revealed that sediments were transported by littoral and aeolian currents and deposited in a highenergy subaqueous environment. In comparison with the Palma Sola beach sediments, Riachuelos beach sediments were lower in SiO2and higher in Sc, V, Co, and Ni contents, suggesting a higher proportion of detritus derived from intermediate igneous rocks. Statistical analysis revealed that the major, trace and rare earth elements were generally associated with terrigenous rather than biogenic materials.The chondrite-normalized REE patterns for the Riachuelos and Palma Sola beach sediments were enriched in LREE, with negative and positive Eu anomalies,indicating their felsic and intermediate igneous sources. A passive margin setting inferred from the tectonic discriminant function diagrams for the beach areas was consistent with the geological setting of the Gulf of Mexico.

Th/U ratios of zircon grains were ＞0.1, indicating a magmatic origin for the Riachuelos and Palma Sola beach sediments. The Riachuelos and Palma Sola beach sediments were dominated by the Paleozoic and Cenozoic zircons, and Proterozoic zircons were limited.The Zacatecas and Nazas Formations in the Mesa Central Province were the most likely source for Paleozoic detrital zircons. The Cenozoic detrital zircons implied that the source terrane, which delivered sediments to the Riachuelos and Palma Sola beach areas, was the nearby Eastern Alkaline Province. Significant discrimination in zircon REE patterns with respect to age was not identified, which suggested that the zircon U-Pb geochronology is a better method to identify sediment provenance than zircon REE patterns and abundances. The resemblance in zircon surface microtextures and zircon age clusters between the Riachuelos and Palma Sola beach sediments demonstrated a similarity in their source rocks, transport distance from the source area,and wave energy.

7 Supplementary information

Supplementary information accompanies this paper at https://doi.org/10.1186/s42501-020-00075-9.

Additional file 1 Supplementary Information 1: Quantitative compositional data of sand grains from the Riachuelos and Palma Sola beach sediments obtained by Scanning Electron Microscopy (SEM) with Energy Dispersive X-Ray Spectroscopy (EDS). Supplementary Information 2: Trace (Nb, Hf, Th, and U) and rare earth element data (in ppm)for the detrital zircon grains in the Riachuelos (Sample No. RIA4) beach sediments, Gulf of Mexico. Supplementary Information 3: Trace (Nb,Hf, Th, and U) and rare earth element data (in ppm) for the detrital zircon grains in the Palma Sola (Sample No. PAL20) beach sediments, Gulf of Mexico. Supplementary Information 4: U-Pb ages of detrital zircons collected from the Riachuelos (Sample No. RIA4) beach sediments, Gulf of

Mexico. Supplementary Information 5: U-Pb ages of detrital zircons collected from the Palma Sola (Sample No. PAL20) beach sediments, Gulf of Mexico. Supplementary Information 6: Factor loadings for the Riachuelos beach sediments, Gulf of Mexico. Supplementary Information 7: Factor loadings for the Palma Sola beach sediments, Gulf of Mexico.Additional file 2 Supplementary Information 8: Post-ablation cathodoluminescence images of representative detrital zircons from the Riachuelos beach sediments. Supplementary Information 9: Post-ablation cathodoluminescence images of representative detrital zircons from the Palma Sola beach sediments.

Acknowledgements

Financial assistance extended by the CONACyT Ciencia Basica (A1-S-21287)and PAPIIT (IN107020) projects are highly acknowledged. I appreciate Mario Marca Castillo for sample collection along the Riachuelos and Palma Sola beach areas. I gratefully acknowledge Mayla A. Ramos-Vázquez and Vanessa del Carmen for their help in zircon separation from beach sediments.I am indebted to Susana Santiago-Perez,Edwardo Alfredo Morales,and Ricardo Martínez Domínguez for their assistance in the Sedimentology Laboratory during sample preparation. I thank the staffs Patricia Girón García, Carlos Linares-López, Laura Elena Gómez Lizárraga, and Teodoro Hernández Trevi?o for their help in XRD, SEM-EDS, SEM, and mineral identification. I also like to thank Dr. Luigi Solari, Carlos Ortega-Obregón, and Harim Arvizu Gutiérrez at the Laboratorio de Estudios Isotópicos, for their assistance in zircon U-Pb dating by LA-ICP-MS. I appreciate Edith Xiadani C. Zárate, “Ayudante de Investigador Nacional Nivel III”, for the help in the preparation of Supplementary Tables. I express my sincere gratefulness to the editor and three anonymous reviewers for their excellent suggestions, which significantly improved the presentation.

Author's contributions

AAJS wrote the manuscript and approved the final revised version.

Funding

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

Availability of data and materials

All data discussed in this study are available in the current Tables, Figures, and Supplementary material(Additional Files 1 and 2:Supplementary Information 1-9)of this manuscript.

Competing interests

The author declares that he has no competing interests; and, there are no other persons who satisfied the criteria for authorship but are not listed.

Received: 16 April 2020 Accepted: 16 September 2020