Abstract
Beds of free-living coated nodules (coralline algae, bryozoans, acervulinid foraminifera) create shallow-water carbonate biogenic benthic habitats, which are sensitive to human disturbance and slow to recover. Holocene bryoliths, ranging from sub-spheroidal, sub-discoidal to sub-ellipsoidal in shape, were found scattered in the foredunes in ca. 30-km stretch along the Hermenegildo and Concheiros do Albardão beaches on the southernmost coast of Brazil (Santa Vitória do Palmar municipality, Rio Grande do Sul State). The dominating bryozoan species forming the bryolith is Biflustra holocenica Vieira, Spotorno-Oliveira and Tâmega sp. nov. The inner bryolith arrangement, generally asymmetrical, shows multilamellar and circumrotatory growth of colonies that envelop the bivalve Ostrea puelchana. Bryozoans and subordinate corals characterize the outer bryolith surfaces. The ichnogenera Gastrochaenolites (made by the boring bivalve Lithophaga patagonica) and Caulostrepsis occur throughout the bryoliths, from the inner part up to the outer surface. The studied bryoliths, originated in a shoreface setting at ca. 7910–7620 cal. yr BP and during subsequent storm waves, were resedimented onto the foreshore and foredunes (to ca. 5700 cal. yr BP) where the bryoliths were finally fossilized.
Introduction
Beds of free-living nodules are built by marine encrusting calcareous organisms growing around a hard nucleus on mobile substrates and are overturned by waves, currents or bioturbation (e.g. Moissette et al., 2010; Prager and Ginsburg, 1989). Macroids (i.e. macro-oncoids or large oncoids), organic, centimetre-sized (>10 mm) coated grains, are formed predominantly by encrusting metazoans or protozoans, sometimes associated with calcareous algae (Bassi et al., 2012; Hottinger, 1983; Peryt, 1983). The beds provide a stable habitat on which seaweeds, coralline red algae and invertebrates can thrive. Fossil and modern free-living nodules made predominantly of a single species of an encrusting sclerobiont include rhodoliths formed by coralline algae (e.g. Aguirre et al., 2017; Baarli et al., 2012; Bosence, 1983), coralliths by corals (e.g. Glynn, 1974), bryoliths by bryozoans (ectoproctaliths, bryoncoids, bryooids; for example, Nebelsick, 1996; Scholz, 2000) and ostreoliths by oysters (Wilson et al., 1998).
Bryoliths formed by cheilostome and cyclostome bryozoans are found on sand and shell-gravel substrates (James et al., 2006). Among cheilostomes, the main taxa producing bryoliths are antroporids (Nishizawa, 1985), celleporids (El-Sorogy, 2015; Nebelsick, 1996; Spjeldnaes and Moissette, 1997), membraniporids (Kidwell and Gyllenhaal, 1998; Klicpera et al., 2013), microporids (Moissette et al., 2010), schizoporellids (Zabin et al., 2010) and smittinids (James et al., 2006). Several bryolith-forming taxa have been reported from the Miocene: Cellepora sp. (El-Sorogy, 2015; Nebelsick, 1996), Celleporina sp., Cerioporidae sp., Tetrocycloecia sp. (El-Sorogy, 2015, as Tretocycloecia) and Turbicellepora sp. (El-Sorogy, 2015; Nebelsick, 1996). Conopeum commensale (Kirkpatrick and Metzelaar, 1922) has been reported both from the Pliocene (Cuffey and Johnson, 1997; Kidwell and Gyllenhaal, 1998; Moissette et al., 2010) and the present day (Klicpera et al., 2013).
The modern Brazilian continental shelf from 3° to 22° S is considered the largest area in the world supporting rhodolith beds (Foster, 2001), ranging in depth from 20 to 250 m (Henriques et al., 2014; Kempf, 1980; Kempf et al., 1969). Pleistocene-to-Recent bryozoan counterparts of the rhodolith beds have been recently studied. In Rio de Janeiro State, Tâmega et al. (2013, 2014) reported a number of bryozoan taxa associated with fossil Holocene (4 spp.) and recent rhodolith beds (20 spp.) surveyed from 94 to 105 m water depth. Barbosa (1967) recognized one single undated bryolith formed by Conopeum commensale collected in the Concheiros do Albardão beach, Santa Vitória do Palmar, Rio Grande do Sul State (RS), Brazil. The bryolith is associated with Pleistocene and Holocene vertebrate remains (e.g. fishes and mammals). Ramalho and Calliari (2015) described eight bryozoan species occurring on the RS shoreface from 12 to 21 m water depth.
In this study, we report for the first time the occurrence of Holocene bryoliths from the southern Brazilian coast at 33° S. The bryoliths are characterized in terms of growth forms, inner arrangement and taxonomic composition, including the description of a new bryozoan species. The palaeoecological analysis of these fossil bryoliths assesses their palaeoenvironmental dynamics recorded in the Hermenegildo and Concheiros do Albardão sedimentary deposits over the past 7900 years.
Study sites
This study was carried out at the Hermenegildo (33°32′49″ S, 53°06′28″ W) and Concheiros do Albardão (33°31′49″ S, 53°05′13″ W) beaches, Santa Vitória do Palmar, RS, Brazil (Figure 1). These beaches are characterized by fine to medium sand in summer (aggregation profile) and mixed composition of coarse sand and bioclastic gravel in winter (erosion profile; Asp et al., 1999; Dillenburg et al., 2004).

Geographic (left) and geological (right) location on the studied area (Hermenegildo and Concheiros do Albardão beaches), Rio Grande do Sul, Brazil. 1, alluvial lake system; 2, Holocene/Pleistocene lagoonal system; 3, Pleistocene barrier; 4, Holocene barrier (modified from Rosa et al., 2017; Tomazelli and Villwock, 1996).
These studied sites belong to the Rio Grande do Sul State coastal plain (RSCP), which is Neogene–Quaternary in age and broadens along the 615 km length (Tomazelli et al., 2000).
In the northern and southern ends of the RSCP, shorefaces are completely different. The southern part of the Lagoa dos Patos inlet (Figure 1) has a gradually decreasing slope towards the sea (Fachin, 1998). The extensive and shallow upper shoreface has a boundary at 10 m water depth with predominantly gravel and sand deposits, whereas the lower shoreface has a boundary at 40 m depth with predominantly fine sand deposits (Toldo et al., 2006). The average significant wave height and period are 1 m and 10–11 s, respectively (Tozzi and Calliari, 2000). During autumn and winter storms, wave height frequently exceeds 2 m and storm surges can reach up to 1.3 m above the modern mean sea level (MSL; Calliari et al., 1998; Parise et al., 2009).
The RS shoreface and inner shelf are characterized by gently steepened bedding, the most prominent features of which are linear sand shoals, marine terraces and exposed hard bottoms of beach rocks. There is no current sediment input, only biogenic contribution by biofouling and platform mud (Buchmann et al., 2001; Buchmann and Tomazelli, 2003). The submerged Pleistocene–Holocene beach rocks are constituted essentially by carbonate rocks. The RSCP is characterized by several marine and continental fossil deposits located both on the continental shelf (Buchmann, 1994; Figueiredo, 1975) and in continental areas landwards of the present-day coastline (Lopes, 2013). On the inner shelf, predominantly Holocene mollusc shell deposits (Ritter et al., 2017) are exposed during autumn and winter erosive events by strong storm waves produced by the passage of extratropical cyclones and transported to the beach (Cruz et al., 2016).
Methods
Sampling methods, bryolith shapes and composition
The bryoliths were found scattered along the foreshore and foredunes, between Hermenegildo and Concheiros do Albardão beaches (Figure 1). A total of 36 bryoliths were collected. To determine the biotic differences from the inner to the outer bryolith parts, the relative coverage of bryoliths and faunal components were estimated by point-counting (50 random points) on photographed 36 bryolith specimens using the image analysis program Coral Point Count with Excel extensions (CPCe; Kohler and Gill, 2006). The measurements of the long, intermediate and short orthogonal axes were taken from each bryolith to determine the sphericity, illustrated by the Triangular Diagram Plotting Spreadsheet (TRI-PLOT; Graham and Midgley, 2000). The amount of boring with respect to the original bryolith arrangement was defined by the bioerosion index (BI; Bassi et al., 2012).
The taxonomic composition and relative cover of the bryolith builders were assessed by one-way analysis of variance (ANOVA; Sokal and Rohlf, 1995). Significant differences were investigated via the Tukey test (Zar, 1999). Bryozoan taxonomic identification was carried out by scanning electron microscopy (SEM) analysis at the Centro de Microscopia Eletrônica do Sul of the Universidade Federal do Rio Grande (CEME-Sul, FURG), Brazil. The studied material was compared with reference bryozoans and types occurring in different bryozoan collections at the Universidade Federal de Pernambuco, Brazil (UFPE); Museu de Zoologia da Universidade Federal da Bahia, Brazil (MZUFBA); Museu de Zoologia da Universidade de São Paulo, Brazil (MZUSP); Smithsonian Institution National Museum of Natural History, USA (NMNH); and the Natural History Museum, UK (NHMUK).
Voucher specimens are deposited at the Coleção Paleontológica do Laboratório de Geologia e Paleontologia, Universidade Federal do Rio Grande (LGP-FURG), Brazil; Museu Oceanográfico Prof. Eliézer de Carvalho Rios, Universidade Federal do Rio Grande (MORG-FURG), Brazil; Coleção de Invertebrados Fósseis do Museu Nacional, Universidade Federal do Rio de Janeiro (MNRJ), Brazil; Coleção Científica de Zoologia–Museu Jóias da Natureza (MJNC), Brazil; Coleção Paleontológica do Laboratório de Estratigrafia e Paleontologia, Universidade Estadual Paulista (LEP-UNESP), Brazil; UFPE, Brazil. Comparative specimens are deposited at NHMUK, London.
Radiocarbon dating
The radiocarbon measurements were carried out by Accelerator Mass Spectrometry (AMS). Two samples of bryolith colonies (LGP-FURG L-016) were dated: the innermost layer and the external layer. Samples were physically and chemically pretreated at the Radiocarbon Laboratory of the Universidade Federal Fluminense (LAC-UFF; Brazil). The thorough bryolith sampling rules out the exogenous material, since the AMS technique requires ca. 20 mg of sample material. The subsamples were treated with 0.1 M hydrochloric acid in order to remove 50% of the sample and any possible contamination. Cleaned solid samples were placed in vials with septum lids in order to be pumped out through a needle. The samples were then hydrolysed with 1.0 mL phosphoric acid and left overnight at room temperature. The tubes were connected to the vacuum line and the gas was purified by means of dry ice/ethanol and liquid nitrogen traps. CO2 samples were then transferred to 9 mm Pyrex tubes previously prepared with titanium hydride and zinc on the bottom and iron powder inside a 6-mm inner tube. Tubes were sealed with an oxy-acetylene torch and heated to 550°C for 7 h so that CO2 could be reduced to graphite (Macario et al., 2015, 2017; Xu et al., 2007). Graphite samples were pressed in aluminium cathodes and their radiocarbon concentration was measured at the Center for Applied Isotope Studies at the University of Georgia (Cherkinsky et al., 2010).
Radiocarbon ages were calibrated using the Marine13 curve (Reimer et al., 2013) in order to account for the marine reservoir effect (MRE; see, for example, Alves et al., 2018). A value of 17 ± 56 14C years was used for the local offset from the average global marine reservoir (Alves et al., 2015). Apart from the independent individual calibration, the two dated samples can be related by growth order, as the inner sample can be assumed to be older than outer one. This prior knowledge allows for the use of sequential Bayesian modelling (Bronk Ramsey, 2008). Therefore, a simple sequence model was used for calibration of the set of two samples within the OxCal software (Bronk Ramsey, 2008).
X-ray diffraction analyses
Four samples were collected for mineralogical analyses: two from the outer surface and two from the inner part of two bryoliths (LGP-FURG L-012 and L-016). Approximately 3 g of powdered samples was submitted for X-ray diffraction (XRD) analysis at CEME-Sul of FURG, using a Brucker AXS D8 Advance model (Cu Kα radiation, 40 kV, 40 mA). This was operated in a Bragg–Brentano θ/θ configuration, with the diffraction patterns being collected in a flat geometry with steps of 0.05 degrees and accumulation time of 2.0 s per step using a PSD detector (Bruker AXS LynexEye model). The XRD data were refined following the Rietveld method with the GSAS-II software (Toby and Von Dreele, 2013).
Results
Bryolith shapes and composition
The studied bryoliths show encrusting growth forms (Figure 2(a)–(g)).

Outer bryolith surfaces showing the encrusting growth forms and the ichnotaxa Caulostrepsis (Ca) and Gastrochaenolites (Ga) along with the identified taphonomic signatures. The bryoliths (a–e) are mainly sub-spheroidal in shape with subordinate sub-ellipsoidal forms (f and g). Note the highly abraded (abr) bryolith outer surfaces. Bry: bryozoans; cor: corals (Astrangia rathbuni); Disc: discoidal; Ellip: ellipsoidal; Sph: spheroidal. Scale bar represents 4 cm.
The maximum bryolith diameter ranges from 7.5 to 31.0 cm (mean = 19.1 cm; n = 36) and the minimum diameter ranges from 3 to 18 cm (mean = 10 cm; n = 36) with sub-spheroidal, sub-discoidal and sub-ellipsoidal shapes (Figure 2). Four samples consist of sub-spheroidal coral cobbles with a mean maximum diameter of 8.25 cm (n = 4). The coral cobbles are composed of a single zooxanthellate coral species, Astrangia rathbuni Vaughan, 1906.
All collected bryolith specimens are composed mainly by a single encrusting bryozoan species herein identified as a new species: Biflustra holocenica sp. nov. (Figures 5–7). Two subordinate bryozoan taxa were also recognized: Conopeum sp. and Schizoporella sp.
The inner bryolith arrangement is generally asymmetrical, with multilamellar and concentric growths enveloping oyster bivalve shells, sometimes articulated. Shells of Ostrea puelchana d’Orbigny, 1842, were found preserved as the bryolith nucleus (Figure 3(a)–(c), (e)–(g)). In the inner bryolith part, subordinate components (less than 2%) are represented by the bivalves Sheldonella bisulcata (Lamarck, 1819), Glycymeris longior (Sowerby, 1833), Lithophaga patagonica (d’Orbigny, 1847), Mactraisabelleana (d’Orbigny, 1846), Nucula sp., Ostrea equestris Say, 1834, Phlyctiderma semiaspera (Philippi, 1836), Sphenia fragilis (Adams and Adams, 1854), Plicatula gibbosa (Lamarck, 1801) and by Polychaeta (Serpulidae). No constructional voids (sensu Aguirre et al., 2017) are present.

Slab surfaces of studied bryoliths showing the inner arrangement, the building components, the morphologies of common bioerosion traces and taphonomic signatures. Bal: balanids; bry: bryozoans; ost: ostreid; p: polychaetes; Ga: Gastrochaenolites; abr: abrasion; dis: dissolution. Scale bars represent (a and b, d–g) 4 cm, (c) 2 cm.
Bryozoans, corals and bioerosion traces are the dominant components on the outer bryolith surfaces with a relative cover of 54%, 10% and 33%, respectively. Significant differences in relative cover and associated fauna were detected in the sampled bryoliths (ANOVA F = 72.85, p ⩽ 0.0001; Figure 4(a)). In the outer bryolith surfaces, subordinate components are represented by the bivalves Ostrea equestris and the barnacles Amphibalanus sp. Azooxanthellate coral colonies of A. rathbuni were found associated to 25 bryoliths.

(a) Relative cover (%, mean ± SE) of benthic organisms and bioerosion traces on bryoliths surfaces sampled. Different letters above bars indicate significant differences among means detected by Tukey’s test (p < 0.05). (b) The modelled datings of the sample LGP-FURG L-016.
Two bioerosion traces were identified within the bryoliths: Gastrochaenolites Leymerie, 1842 and Caulostrepsis Clarke, 1908. The BI is 3. Gastrochaenolites is the most common trace within the samples and occurs randomly spread in all bryoliths. Gastrochaenolites (Figures 2(a), (c)–(g) and 3(a)–(c), (f) and (g)) is a clavate (club-shaped) boring trace. The aperture of the boring is slightly narrower than the main chamber (forming a neck of variable length) and has a circular to ovoid shape with a diameter up to ca. 14.44 mm. The producer of Gastrochaenolites is the bivalve Lithophaga patagonica, preserved within some borings.
Caulostrepsis occurs as cylindrical galleries usually with rectilinear trajectories. The cross-sections are generally sub-elliptical in shape, with rare eight-shaped perimeters (Figure 2(a)). The traces, up to 3 mm in length, are randomly arranged, roughly equidistant each other, locally perpendicular or oblique to the bryolith surface.
Radiocarbon dating and XRD analyses
Radiocarbon age estimates for the bryolith’s nucleus and outer surface are statistically similar (Table 1) preventing the estimation of longevity or growth rates. Considering the growth sequence, the calibrated and modelled dates (Figure 4(b)) provide similar results with up to 300 years range at 7910–7600 cal. yr BP (2σ)
Sizes and analytical data for conventional radiocarbon calibrated and modelled age estimates (95.4%) for the nucleus and outer bryolith sample.
Interm: intermediate; cal: calibrated.
The stable isotope for the inner bryolith part reflects the isotopic composition of dissolved inorganic carbon (DIC) in the marine water.
Table 2 shows the results of XRD of the inner part and outer surface of the two analysed bryoliths (LGP-FURG L-012 and L-016).
Mineral composition of the outer surface (OS) and the nucleus (N) of analysed bryoliths by weight (%).
Mg-calcite predominates in all layers. From the XRD, it was possible to determine the Mg content as 10 mol% MgCO3, classified as high magnesium calcite (HMC). Aragonite and quartz are also present.
Systematics
Conopeum commensale: Barbosa, 1967; Non Conopeum commensale Kirkpatrick and Metzelaar, 1922

Bryoliths mainly formed by Biflustra holocenica sp. nov. (a) The outer surfaces of bryolith LGP-FURGL-007 with holotype on surface. (b) Bryolith LGP-FURGL-007 showing the inner asymmetrical arrangement, with multilamellar and circumrotatory growth colonies that envelop the shell fragment. (c) Bryolith (MNRJ-4887-I) formed by bryozoan previously assigned to Conopeum commensale by Barbosa (1967) and here reassigned to the new species. (d) The inner arrangement growth colonies (MNRJ-4887-I). Scale bar represents 2 cm.

SEM photos of holotype specimens of Biflustra holocenica sp. nov. (on the surface of bryolith LGP-FURG L-007): (a) Detail of the frontal surface of the colony; (b) close-up of the autozooids, showing the lateral and distal pore plates; (c) close-up of three kenozooids; (d) internal lateral walls of the autozooids, showing the multiporous mural septula (black arrows) and internal spinules (white arrow); and (e) internal view of one autozooid, showing paired hook-shaped spinules in the lateral zooidal walls. Scale bar represents (a) 500 µm, (b–d) 200 µm and (e) 50 µm.

SEM photos of living specimens of Biflustra holocenica sp. nov (MORG 52680): (a) Detail of the skeletal autozooids; (b) close-up of the proximal cryptocyst, showing the proximal denticle (white arrow); (c) internal lateral walls of the autozooids, showing the internal hook-shaped spinule (black arrow); and (d) lateral view of the autozooids, showing frontal membrane with chitinous spinules (white arrows) and internal hook-shaped spinule in the lateral zooidal walls (black arrow). Scale bar represents (a) 200 µm, (b) 50 µm, (c) 25 µm and (d) 100 µm.
Among the Biflustra species, spinules at basal and vertical walls are known only in Biflustra denticulata (Busk, 1856) and Biflustra paulensis (Almeida et al., 2017; Marcus, 1937). Biflustra denticulata (Figure 8(a)), however, differs in having cryptocystal denticles projecting into the opesia, and in having the proximal cryptocyst narrower than that of Biflustra holocenica sp. nov. Biflustra paulensis has gymnocystal tubercles (absent in Biflustra holocenica sp. nov.) and an opesia occupying three-quarters of the zooidal length (opesia occupies no more than two-thirds of the zooidal length in Biflustra holocenica sp. nov.). The well-known Acanthodesia commensale (Figure 8(b); Klicpera et al., 2013) forming bryoliths, also referred to as Biflustra commensale (Almeida et al., 2017), has large paired gymnocystal tubercles that may be coalescent, and the autozooids are outlined by a brown line, characteristics that are not seen in Biflustra holocenica sp. nov. In addition, the internal hook-shaped spinules, characteristic of Biflustra holocenica sp. nov., are absent in Biflustra commensale.

SEM photos of three morphologically related species: (a) Holotype specimen of Biflustra denticulata (Busk, 1856) (NHMUK 1892.9.6.1) from Mazatlan, Gulf of California; (b) part of the type specimen of Biflustra commensale (Kirkpatrick and Metzelaar, 1922) (NHMUK 1922.9.9.9) from Cabo Blanco, West Africa; (c) specimen of Biflustra arborescens (Canu and Bassler) (UFBA 1258) from Bahia, Brazil. Scale bar represents 200 µm.
Some Recent colonies of Biflustra collected on living gastropods from 40 km north off Hermenegildo beach (MORG 52680) have a skeletal morphology (Figure 7(a)–(d)) similar to the studied Holocene holotype and paratypes. Despite the presence of a single minute proximal denticle projecting into the opesia (Figure 7(b)) only in the recent colonies (lost in the studied fossil material during the taphonomic processes), we consider both fossil and living specimens to be conspecific. In addition, organic parts of living colonies were examined, revealing the presence of external chitinous spinules on the frontal membrane of autozooids (Figure 7(d)), a feature not seen in the fossil specimens.
A frontal membrane with spinules is also described in Biflustra arborescens (Canu and Bassler, 1928) recorded from Brazil (Almeida et al., 2017). Biflustra holocenica sp. nov. and Biflustra arborescens share other similarities in skeletal morphology, including autozooidal size, absence of gymnocyst, cryptocyst occupying one-third of zooidal length and morphology of multiporous mural septula in vertical walls. Biflustra holocenica sp. nov. differs from Biflustra arborescens, in having basal and vertical walls with spinules, and paired hook-shaped spinules at proximal third of lateral vertical walls (absent in Biflustra arborescens). The single minute proximal denticle projecting into the opesia is rarely present in Biflustra holocenica sp. nov., while these denticles may be numerous in Biflustra arborescens (Figure 8(c)).
In total, 32 Biflustra species are registered worldwide (Bock and Gordon, 2018), with 8 species distributed along the Brazilian coast, from Alagoas to Santa Catarina State. These species may be found from shallow waters to 100 m depth, being abundant between 5 and 30 m water depth (Almeida et al., 2017; Vieira et al., 2008, 2016). Biflustra is recorded here for the first time in the RS as a fossil specimen (ca. 8000 cal. yr BP forming bryoliths). Present-day Biflustra holocenica sp. nov. are found growing as sheet-like colonies on the external surface of living gastropod shells, while one colony was found on a beached styrofoam buoy.
Discussion
The studied Holocene bryoliths were found in foredunes. Their occurrence in a subaerial setting suggests that (1) the bryoliths were redeposited according to the relative sea-level (RSL) changes; (2) they formed on a coarse, soft, substrate; and (3) they were resedimented onto foredunes during the Holocene coastal dynamics that affected the RSCP area. The following discussion addresses details of the bryolith’s radiocarbon datings and compositional and morphological features to assess the RSCP palaeoenvironmental dynamics.
The radiocarbon dating
The radiocarbon ages are affected by some factors such as radiocarbon MRE, incorporation of dead carbon from limestone, possible aggregation of boring bivalves and recrystallization of aragonite and HMC. MRE contributes to the depletion of 14C concentration in marine samples leading to apparent old ages (Alves et al., 2018). Such effect was considered through calibration with the marine curve and the use of local reservoir offset estimates. Aggregation of younger bivalve shells that would result in underestimated ages is ruled out as only the bryolith layers were sampled. Incorporation of dead carbon fraction within radiocarbon-dated bryoliths would also result in an overestimated age of the bryoliths (e.g. Wiener, 2012). However, there is no evidence of limestone deposits in this region. Finally, aragonite and HMC can dissolve and recrystallize into low magnesium calcite (LMC). This process could facilitate the exchange of carbon with the environment, increasing the chances of younger carbon being incorporated into the lattice. As the XRD results show a composition of HMC and aragonite, recrystallization is unlikely to have occurred.
Assuming that no exogenous carbon was incorporated to the sample, the calibrated radiocarbon ages of the LGP-FURG L-016 bryozoan sample are ca. 7910 to 7600 cal. yr BP (Figure 9).

Reconstruction of the Holocene palaeoenvironmental dynamics for the studied bryoliths in the Hermenegildo and Concheiros beaches, southern Brazil. The arrows refer to the trend of the sea level (fall or rise). The sample ages are based on radiocarbon dating (this study). The two dated bryolith nuclei (175/188) plotted with Holocene sea-level curves of Martin et al. (1979), Angulo et al. (2006), Suguio et al. (2013) and Prieto et al. (2017). (a) At ca. 7900 cal. yr BP, the studied bryoliths formed on coarse soft substrate, in the shoreface part of a barrier-lagoon depositional system (Lima et al., 2013). (b) Until the PSLM (ca. 5700 cal. yr BP), the bryoliths were resedimented onshore onto foreshore and foredunes. (c) After the end of the PSLM to the present, during a slow sea-level fall, the bryoliths were preserved in the present-day beach foredunes. MSL: mean sea level; dep syst: depositional system; PSLM: postglacial sea-level maximum (5100–5700 cal. yr BP;Angulo et al., 2006). Not to scale.
Palaeoecology of the Holocene bryoliths
At the estimated age of the studied bryoliths, the sea level was −5 m below the present sea level (e.g. Prieto et al., 2017). The RSCP has some particular characteristics concerning the shoreface bathymetry, and the inner continental shelf is a result of transgressive and regressive events caused by glacial and postglacial sea-level changes (Calliari and Klein, 2011). As the shoreface of the studied area is wide and shallow with seaward limit at 10–15 m depth (Toldo et al., 2006), at the time of the bryolith growing age, the shoreface should have been about 20 m below the present sea level.
The foredunes formed at ca. 5700 cal. yr BP (Figure 9(b)) during the postglacial sea-level maximum (PSLM) when the sea level was ca. 2 m above the present sea level (e.g. Angulo et al., 2006; Prieto et al., 2017). This time slice was ca. 2000 years after the bryolith formation. The bryoliths should not, therefore, be deeply burrowed under the substrate as they would be easily re-exhumed. Therefore, if the bryoliths formed in the shoreface setting, they should have undergone a resedimentation ranging in elevation from ca. 17 to 22 m, that took them from their factory setting onto the subaerial one.
The likelihood of a shallow burial is supported by very low sedimentation rate of the area as confirmed by moderate bioerosion (e.g. Bassi et al., 2013), three orders of magnitude difference in the shell ages of bivalve shells (Ritter et al., 2017) and the presence of the ichnotaxa Gastrochaenolites (Bromley and Asgaard, 1993), most common and produced by mytilid bivalves, and Caulostrepsis, attributed to the boring spionid polychaete Polydora (Domènech et al., 2008), eucinids (Martinell and Domènech, 2009) and cirratulids (Bromley and D’Alessandro, 1983; Taylor and Wilson, 2003). It also demonstrates that these bryoliths, as other macroids, are prone to span at least some time at the sediment–water interface before the burial due to their relatively high durability. In the studied area, the bryoliths act as spotted hard substrates because locally they are supporting a high biodiversity, being an important colonization surface for sclerobionts together with molluscs (Agostini et al., 2017).
According to Kidwell and Gyllenhaal (1998), a balance of low sedimentation rate (to avoid burial), intermittent high energy (to overturn episodically) and appreciable growth rate of the encrusting organisms (Hottinger, 1983) is required for macroids development. Fast growing organisms, like membraniporid bryozoans, are able to settle on any consolidated fragment or on bioclasts (as the oyster Ostrea puelchana) and increase rapidly in size, quickly covering the colonized substrate.
A relatively shallow setting is also argued by the morphological characteristics of the nodules and the associated identified biota.
Although no shape–water energy relationship has been so far assessed for modern and fossil bryoliths, the shape of rhodoliths is commonly used as an indicator of turbulence and how often they are overturned (e.g. Aguirre et al., 2017). Sub-spheroidal shapes, along with common sub-ellipsoidal and rare spheroidal ones, as seen in Figure 2, are indicative of easy roll on the mobile substrate. Waves and tidal currents are the most common agents for such movements, but bioturbation may also play a role (e.g. Aguirre et al., 2017; Bassi et al., 2006; Foster, 2001; Marrack, 1999). All such agents may have been active in the Holocene RSCP shallow shoreface setting. The dominating encrusting growth form suggests a constant water motion induced by wave/current bed setting as concluded for rhodoliths (e.g. Marrack, 1999).
In the studied materials, the nucleus consists of an articulated specimen of the oyster Ostrea puelchana. This oyster shows a wide depth distribution from 10 to 100 m water depth (Rios, 2009; Romero et al., 2013).
In the studied bryoliths, the identified coral Astrangia rathbuni is very common. The present-day setting for Astrangia rathbuni is related to biogenic substrates from inner to mid-continental shelf. Astrangia rathbuni is restricted to shallow water (18–81 m) being most common at 48 m depth (Kitahara et al., 2009). Considering that the studied bryoliths show almost no important breakages (e.g. Figures 2(a)–(g), 3(a)–(d), 5(a) and (b)), high abrasion (Figures 2(a)–(g), 3(c) and (d)) or size selection (Figure 2), this indicates that the bryoliths were most probably re-exhumed and transported from a shallow setting, which is in good agreement with the shallowest bathymetric distributions of Ostrea puelchana and Astrangia rathbuni. Modern bryoliths are found in subtidal marine environments at a depth of 7 m outside Punta Chivato, Gulf of California (James et al., 2006).
Biflustra holocenica sp. nov. has a bi-mineralic skeleton made up of HMC and aragonite. This composition has been so far demonstrated for a minority of bryozoan species, whereas most of them have LMC (Kuklinski and Taylor, 2009; Smith et al., 2006; Taylor et al., 2009, 2015). The colonies can therefore be attributed to the major groups of Cheilostomes (Malacostegina, Flustrina, Umbonulomorpha and Lepraliomorpha) with aragonite skeletons (Taylor et al., 2015).
Chave (1954) showed a positive correlation between the magnesium content of the Mg-calcite skeletons and the temperature of the water in which the organism lived, with higher values of Mg in warmer water. However, Smith et al. (1998) found no consistent mineral trend with temperature for Bryozoa around New Zealand, and Arctic species show no correlation with temperature (Kuklinski and Taylor, 2009). Data compilations (Kuklinski and Taylor, 2009; Taylor et al., 2015) show that a higher proportion of Bryozoa have higher Mg contents in lower latitudes; however, the absence of correlation with temperature at specific locations such as New Zealand and the Arctic suggests caution must be exercised when trying to determine palaeotemperature conditions based on just the Mg content from one sample set. The 10 mol% MgCO3 for our study bryoliths is consistent with the range of values found for Bryozoans from mid- and low-latitude locations (Taylor et al., 2015).
The climatic scenario for the marine environment along the southwestern Atlantic during the mid-Holocene indicates the sea surface temperature was slightly warmer, between ca. 32ºS and 47ºS (Aguirre et al., 2009), following a latitudinal decrease toward higher latitudes. This means that temperatures were warmer than those recorded at present at the same latitudes, as also suggested by Martínez et al. (2001, 2013) for the Uruguayan coast, allowing the macroids formation. Thus, it may be that our Bryozoan 10 mol% MgCO3 is reflective of higher temperatures at that time. However, this proposition cannot be confirmed without either experimental work or sampling along a temperature gradient, to determine the response of the Mg content to changes in temperature for this species to give our results context.
Palaeoenvironmental dynamics
The Hermenegildo sedimentary successions represent a barrier–lagoon depositional system formed about 9800 years ago during the sea-level transgression (Buchmann et al., 1998; Lima et al., 2013). The RSL had overtaken present MSL by ca. 7000 cal. yr BP (Angulo et al., 2006; Prieto et al., 2017; Figure 9).
At ca. 6700 cal. yr BP, the lagoonal floor sedimentation was active with a permanent connection to the ocean in the form of single inlets (Buchmann et al., 1998; Lima et al., 2013). The balance profile of the beach changed, making available a large amount of sedimentary material, from the reworking of the Pleistocene beaches and terraces.
The subsequent high stand occurred between 5800 and 5200 cal. yr BP (PSLM) reaching as high as +2 m above present MSL. In the studied area, this RSL rise period is characterized by transgressive foredunes moving into the lagoon. The dune transgression has been occurring at the same time as the erosion of the lagoonal barrier (Buchmann et al., 1998; Lima et al., 2013). During this RSL rise, storm surges resedimented the studied bryoliths from the shoreface setting where they formed, onto the foreshore and foredunes (Figure 9(b)). The erosive processes are attributed to storm surge and concentration of wave on the beach caused by the presence of lithified structures in the backshore and to negative sediment budget (Calliari et al., 1998; Dillenburg et al., 2004).
If bryoliths had been transported with storms from offshore, they would have been piled up in a wedge-shaped abutment against the foredune face (e.g. Johnson et al., 2012).
Resedimentation of large nodules has been interpreted for offshore and inshore rhodolith deposits (e.g. Bassi et al., 2010; Checconi et al., 2010; Johnson et al., 2012). In the studied area, resedimentation is likely to be brought about by combinations of winds and storm waves as well as the higher energy of longshore currents whose present-day southerly waves are dominant 30% of the time (Lima et al., 2001).
At the maximum level of PSLM, the transgression of the barrier during the sea-level rise was complete (Lima et al., 2013; Figure 9(b)).
At Hermenegildo and Concheiros do Albardão areas, after the end of the PSLM to the present, during a slow sea-level fall of approximately 2 m (Figure 9(c)), the progradation occurred in coastal embayments due to the onshore transfer of sands from the shoreface of RS (Dillenburg et al., 2000). It was in the foredunes formed during this time that the bryoliths were finally buried.
Concluding remarks
A new species of bryozoan has been described herein for the Southern Brazilian shelf, Biflustra holocenica Vieira, Spotorno-Oliveira and Tâmega sp. nov. This fossil bryolith occurring in the Hermenegildo and Concheiros foredunes represents a witness of a Holocene marine shoreface ecosystem that probably disappeared ca. 7600 years ago. These bryoliths are made up of the dominant component Biflustra holocenica still occurring in this area. However, although small colonies found living on shells of gastropods are still present in the area, they are not able to form bryoliths as described in this study.
Bioerosion traces of two kinds were identified within the bryoliths: Gastrochaenolites and Caulostrepsis. This demonstrates that these bioconstructions spend time at the sediment–water interface, similarly to mollusc shells, thus also acting as an important colonization surfaces for sclerobionts in a sand-dominant shelf.
During about 2000 years of sea-level rise, until the PSLM, the bryoliths were resedimented from the shoreface setting onto the transgressive foredunes, where they currently occur.
Currently, little is known about the contribution of bryoliths and rhodoliths to modern coastal dunes, or whether this type of formation is a prominent end point of bryolith/rhodolith taphonomy (Gulf of California, James et al., 2006, and Sewell et al., 2007; eastern Mediterranean, Titschack et al., 2008; Madeira Archipelago, Soares, 1973). The results of our study indicate that these coastal dunes can act as prominent end points in our study location. Further investigations on a wider range of modern and ancient coastal dune systems are required to identify the extent of and type of bryoliths/macroids/rhodoliths and extensive bioclasts derived from bryozoans, encrusting foraminifera and coralline red algae.
Footnotes
Acknowledgements
We are grateful to the Coordenação de Aperfeiçoamento de Pessoal de Nível Superior-CAPES for providing postdoc fellowships (FTST and PSO). This study was financed in part by the Coordenação de Aperfeiçoamento de Pessoal de Nível Superior - Brasil (CAPES) - Finance Code 001. We thank the Universidade Federal do Rio Grande (FURG) for field and laboratory support. We are also grateful to Luiz Rota de Oliveira for donating some bryoliths studied herein. Rudmar Krumerick and Caroline Ruas (Centro de Microscopia Eletrônica do Sul (CEME-SUL), FURG) are thanked for the scanning electron microscopy (SEM) images of Brazilian specimens; Fabrizio Scarabino for the identification of some molluscs and Paulo Matioli for the MJNCA04 bryolith information. Ana Carolina Sousa de Almeida provided SEM of Biflustra arborescens specimens. Paul Taylor and Mary Spencer Jones provided SEM images of the comparative type material of bryozoans deposited at the Natural History Museum, London. Comments and suggestions from Matias Ritter and two anonymous reviewers improved the early version of the manuscript.
Funding
KM would like to thank Brazilian financial agencies CNPq (Conselho Nacional de Desenvolvimento Científico e Tecnológico, 307771/2017-2 and INCT-FNA, 464898/2014-5) and FAPERJ (Fundação Carlos Chagas Filho de Amparo à Pesquisa do Estado do Rio de Janeiro, E-26/110.138/2014 and E26/203.019/2016) for their financial support.
