Storm-Related Rhodolith Deposits from the Upper Pleistocene and Recycled Coastal Holocene on Sal Island (Cabo Verde Archipelago)

: This project examines the role of tropical storms in the northeast Atlantic Ocean related to the post-mortem deposition of rhodoliths in coastal settings during Neogene to Holocene time with primary emphasis on Sal Island in the Cabo Verde Archipelago located 600 km o ﬀ the coast of Senegal in northwest Africa. Fossil rhodoliths from 10 to 15 cm in diameter are equal in size to contemporary rhodoliths that survive for a century or more at water depths undisturbed by all but the most energetic storms. The shape of large rhodoliths makes them susceptible to rare disturbances with su ﬃ cient energy to export them beyond their preferred habitat into extreme environments that include supratidal settings. The methodology of this study gauges the relative sphericity of rhodoliths based on measurements across three axes perpendicular to one another, plots size variations on bar graphs, and considers whether or not individual nodules are nucleated around rock cores eroded from proximal rocky shores. Sal Island is impacted on a steady basis by wave swell generated from the Northeast Trade Winds, but Pleistocene and Holocene deposits with large rhodoliths on the Island’s windward coast are interpreted as the result of major storms of hurricane intensity. Comparison of Sal Island rhodoliths with Pliocene and Miocene examples from other insular localities in the Northeast Atlantic considers evidence for displacement of the Inter-Tropical Convergence Zone (ITCZ) into more northern latitudes as an inﬂuence on past hurricane tracks that are less common today.


Introduction
Rhodoliths are free-rolling nodules formed by coralline red algae. These algae secrete high-magnesium calcite in their cell walls in a concentric pattern of outward growth expressed by surface crustose growth or as radial branches resulting in nodules of spheroidal to ellipsoidal shape. Fifty-nine genera currently recognized in the subclass Corallinophycidae within the division Rhodophyta trace an evolutionary history back to the early Cretaceous [1]. Global in distribution on oceanic shelves bordering many coastlines, the latitudinal range of rhodoliths extends from places as distant as New Zealand in the southern hemisphere to Spitsbergen in the northern hemisphere [2]. Penetration of sunlight through marine water is the primary limiting factor to the sustainable growth of rhodoliths, although these have been found living at the limit or even slightly below the photic zone. Ranging in diameter between that of a large grape to that of a large grapefruit, live rhodoliths occur patterns of long-term atmospheric circulation and storm dynamics in the region through Neogene time. To this end, Sal Island Pleistocene and Holocene rhodoliths are compared with older Pliocene and Miocene rhodoliths farther afield in the Cabo Verde Archipelago and neighboring islands in the northeast Atlantic Ocean [12][13][14][15] as well as mainland Portugal [16] in order to underscore features that are peculiar to Sal Island.

Geographical and Geological Setting
The Cabo Verde Archipelago is the most southern of several volcanic archipelagos that extend from the Azores at a latitude bordering N 40 • approximately 1500 km due west of Portugal to the leeward islands of Cabo Verde at a latitude bordering 15 • N (Figure 1a). Sal is one of the windward islands in the northern sector of the archipelago, which sits 600 km west of Senegal off the northwest coast of Africa (Figure 1b). It is the smallest of the populated islands in the group [18] with an area of 216 km 2 , circumference of 83 km, and maximum elevation of 406 m (Figure 1c). The island occupies a basalt core complex that dates to 26 Ma [19] and the enclosing 50-m isobath approaches nearest to the north end, where strong surf erodes the rocky coast due to direct impact of the persistent northeast trade winds as well as strong North Atlantic swell from the northwest [20]. The Northeast Trade Winds reaching Sal from the Canary Islands commonly register 5 to 6 on the Beaufort Scale, equivalent to speeds between 10.5 and 48.3 m/s (= 38 to 48.3 km/h) [21] and capable of generating surface swells with an amplitude of 3.5 m. Unimpeded by low topography, the trade winds sweep over Sal from a northeasterly direction and deflate carbonate-rich beach sand to form a succession of small dunes that climb through narrow valleys such as the Ribeira do Lavrador (Figure 1c), enriched from 50 to 67% in grains derived from crushed rhodoliths [11]. Except for the rugged northern shores of the island, fossil and present-day rhodoliths are found to occur along all other shores around the island. The focus of this study, however, is on the larger rhodoliths from the windward east coast.
Geosciences 2020, 10, x FOR PEER REVIEW 3 of 18 highlight patterns of long-term atmospheric circulation and storm dynamics in the region through Neogene time. To this end, Sal Island Pleistocene and Holocene rhodoliths are compared with older Pliocene and Miocene rhodoliths farther afield in the Cabo Verde Archipelago and neighboring islands in the northeast Atlantic Ocean [12][13][14][15] as well as mainland Portugal [16] in order to underscore features that are peculiar to Sal Island.

Geographical and Geological Setting
The Cabo Verde Archipelago is the most southern of several volcanic archipelagos that extend from the Azores at a latitude bordering N 40° approximately 1500 km due west of Portugal to the leeward islands of Cabo Verde at a latitude bordering 15° N (Figure 1a). Sal is one of the windward islands in the northern sector of the archipelago, which sits 600 km west of Senegal off the northwest coast of Africa (Figure 1b). It is the smallest of the populated islands in the group [18] with an area of 216 km 2 , circumference of 83 km, and maximum elevation of 406 m (Figure 1c). The island occupies a basalt core complex that dates to 26 Ma [19] and the enclosing 50-m isobath approaches nearest to the north end, where strong surf erodes the rocky coast due to direct impact of the persistent northeast trade winds as well as strong North Atlantic swell from the northwest [20]. The Northeast Trade Winds reaching Sal from the Canary Islands commonly register 5 to 6 on the Beaufort Scale, equivalent to speeds between 10.5 and 48.3 m/s (= 38 to 48.3 km/h) [21] and capable of generating surface swells with an amplitude of 3.5 m. Unimpeded by low topography, the trade winds sweep over Sal from a northeasterly direction and deflate carbonate-rich beach sand to form a succession of small dunes that climb through narrow valleys such as the Ribeira do Lavrador (Figure 1c), enriched from 50 to 67% in grains derived from crushed rhodoliths [11]. Except for the rugged northern shores of the island, fossil and present-day rhodoliths are found to occur along all other shores around the island. The focus of this study, however, is on the larger rhodoliths from the windward east coast.

Materials and Methods
Sampled at random, whole rhodoliths from deposits at two localities were measured (to the nearest half centimeter) across three principal axes (long, intermediate, and short). Both localities are on the east shore of Sal Island and differ in origin from consolidated (cemented) as opposed to unconsolidated (loose) deposits. Data from measurements were subject to analysis based on triangular plots among spherical, ellipsoidal, and discoidal shapes according to the standard format applied to rhodoliths by Bosence (1976Bosence ( , 1983 [22,23] as modified from Sneed and Folk (1958) [24]. In addition, these data are fitted to bar graphs to show variations in rhodolith size based on length of the long and intermediate axes correlated to the Wentworth scale for sedimentary particles in the range of pebbles and cobbles [25]. The volume of each fossil rhodolith was calculated applying the radius from the intermediate axis to the standard equation for the volume of a sphere. Few rhodoliths express the morphology of a perfect sphere, and hence the radius taken as half of the intermediate axis was used to calculate the average volume for a given rhodolith. As it was not possible to record the weight of fossil rhodolith in the field, a substitute method was applied to estimate weight. Dolomite density representative of high-magnesium carbonate was taken as 2.6 g/cm 3 [26] and this value was applied to calculate weight based on rhodolith volume. Moreover, consideration is given to the wave intensity necessary to dislodge rhodoliths from a binding limestone under storm conditions. After obtaining field measurements, each rhodolith was broken open in order to check for possible nucleation around extraneous shells or shell fragments as opposed to basalt pebbles or cobbles. Potentially, an internal basalt cobble or pebble would have the effect of increasing the overall mass of the rhodolith due to the additional density of the basaltic core at about 3.0+ g/cm 3 . Any increase in overall density related to nucleation of a given rhodolith around a basalt pebble or cobble will increase the wave energy required to move it.

Upper Pleistocene Deposit on Baía da Parda
The base map for Sal Island (Figure 1c) identifies the principal localities at Baía da Parda where the data from rhodoliths in a consolidated Pleistocene deposit and those from a loose aggregate at Praia das Gatas (in English: Beach of the Nurse Sharks), were collected, respectively. Low sea cliffs extend for approximately 2 km south from the headland bordering Baia da Parda, where rhodolithic limestone typically protrudes seaward with wave-eroded overhangs between 1 and 1.5 m above mean sea level ( Figure 2). Based on 234 U/ 238 U dating of the fossil coral Siderastrea radians at Baía Parda [27], an age of 130.4 ± 1.4 ka falls within the range of Marine Isotope Substage 5(e), or the last interglacial interval before the present. A succession of other terraces at Serra Negra ( Figure 1c) indicates that Sal Island underwent episodes of uplift beginning in earlier Pleistocene time [28]. However, upper Pleistocene deposits similarly rich in fossil rhodoliths occur at the shoreline as far south as Costa da Fragata [11].
Exposure of the same limestone strata extending a short distance south, but still within Baía da Parda offers a closer view of the dissected rhodolith deposit sitting unconformably on an unconformity above a basalt cobble base ( Figure 3). Scattered basaltic pebbles occur cemented in place by a lime-rich matrix that binds the rhodoliths together, forming a dense algal conglomerate in pockets. Where naturally broken, very few rhodoliths reveal a basaltic pebble at the core. Cementation is relatively weak and individual rhodoliths are easily extracted from the matrix. The outer surface of individual rhodoliths is worn almost smooth, but reveals lumpy growth typical of species in the genus Lithothamnium. Internally, rhodolith growth in this case is the result of algal branches tightly pressed together in an outwardly expanding pattern. Exposure of the same limestone strata extending a short distance south, but still within Baía da Parda offers a closer view of the dissected rhodolith deposit sitting unconformably on an unconformity above a basalt cobble base ( Figure 3). Scattered basaltic pebbles occur cemented in place by a lime-rich matrix that binds the rhodoliths together, forming a dense algal conglomerate in pockets. Where naturally broken, very few rhodoliths reveal a basaltic pebble at the core. Cementation is relatively weak and individual rhodoliths are easily extracted from the matrix. The outer surface of individual rhodoliths is worn almost smooth, but reveals lumpy growth typical of species in the genus Lithothamnium. Internally, rhodolith growth in this case is the result of algal branches tightly pressed together in an outwardly expanding pattern.

Upper Pleistocene Deposit at Praia das Gatas
Beginning south from the cross-island paved road between Palmeira and Salinas (Figure 1c), a dirt road along the east shore runs parallel to the coast, where a small patch of Upper Pleistocene limestone equated with the unconformity at Baía da Parda is exposed in the present-day intertidal zone (Figure 4a). On closer view (Figure 4b), the patch is shown to consist of a single layer of welded

Upper Pleistocene Deposit at Praia das Gatas
Beginning south from the cross-island paved road between Palmeira and Salinas (Figure 1c), a dirt road along the east shore runs parallel to the coast, where a small patch of Upper Pleistocene limestone equated with the unconformity at Baía da Parda is exposed in the present-day intertidal zone (Figure 4a). On closer view (Figure 4b), the patch is shown to consist of a single layer of welded rhodoliths sitting directly on a surface of basalt conglomerate. Like those at Baía da Parda, the rhodoliths are degraded by wear that smooths the lumpy surface texture typical of coralline red algae in the genus Lithothamnium. The patchy deposit also incorporates whole and broken shells belonging to the tropical conch, Thetystrombus latus (Gmelin, 1791) [29]. Some of the smaller shells are complete, but other larger shells are broken, revealing the shell's internal spire structure.

Upper Pleistocene Deposit at Praia das Gatas
Beginning south from the cross-island paved road between Palmeira and Salinas ( Figure 1c), a dirt road along the east shore runs parallel to the coast, where a small patch of Upper Pleistocene limestone equated with the unconformity at Baía da Parda is exposed in the present-day intertidal zone (Figure 4a). On closer view (Figure 4b), the patch is shown to consist of a single layer of welded rhodoliths sitting directly on a surface of basalt conglomerate. Like those at Baía da Parda, the rhodoliths are degraded by wear that smooths the lumpy surface texture typical of coralline red algae in the genus Lithothamnium. The patchy deposit also incorporates whole and broken shells belonging to the tropical conch, Thetystrombus latus (Gmelin, 1791) [29]. Some of the smaller shells are complete, but other larger shells are broken, revealing the shell's internal spire structure.

Configuration of Holocene Supratidal Deposit at Praia das Gatas
At Praia das Gatas, a massive deposit of unconsolidated rhodoliths extends along the shore for a distance of about 75 m ( Figure 5). The seaward sloping deposit occupies a supratidal zone variably 4 to 6 m in width and immediately adjacent to the intertidal zone occupied entirely by basalt boulders. Few basalt clasts are incorporated within the rhodolith deposit, which also includes broken shells belonging to the large conch, Thetystrombus latus.

Configuration of Holocene Supratidal Deposit at Praia das Gatas
At Praia das Gatas, a massive deposit of unconsolidated rhodoliths extends along the shore for a distance of about 75 m ( Figure 5). The seaward sloping deposit occupies a supratidal zone variably 4 to 6 m in width and immediately adjacent to the intertidal zone occupied entirely by basalt boulders. Few basalt clasts are incorporated within the rhodolith deposit, which also includes broken shells belonging to the large conch, Thetystrombus latus. In closer view ( Figure 6), the loosely aggregated rhodoliths from this deposit resemble the worn Pleistocene rhodoliths exposed in the intertidal zone immediately to the north and the far more substantial limestone at Baía da Parda farther to the south. The condition of individual rhodoliths exhibits exterior smoothing by abrasion against neighboring rhodoliths under surf conditions. The In closer view ( Figure 6), the loosely aggregated rhodoliths from this deposit resemble the worn Pleistocene rhodoliths exposed in the intertidal zone immediately to the north and the far more substantial limestone at Baía da Parda farther to the south. The condition of individual rhodoliths exhibits exterior smoothing by abrasion against neighboring rhodoliths under surf conditions. The surface density of rhodoliths varies from 430 to 480/m 2 based on counts that frame sample grids marked off by the corners of a folding meter stick.
Geosciences 2020, 10, x FOR PEER REVIEW 8 of 18 surface density of rhodoliths varies from 430 to 480/m 2 based on counts that frame sample grids marked off by the corners of a folding meter stick.

Comparative Variations in Rhodolith Shape
Raw data on rhodolith size collected from the sample sites at Baía da Parda and Praia das Gatas are listed in Appendix A (Tables A1 and A2). With regard to shape, individual points representing individual rhodoliths are plotted on a set of Sneed-Folk triangular diagrams (Figure 7a,b). The spread of points compared in the two plots demonstrates a high degree of similarity in the variation of shapes between those localities. Points that fall into the top triangle represent the most spherical in shape. Thirty percent of rhodoliths sampled from the Pleistocene limestone fall into this category ( Figure  7a), whereas the majority of points fall within the central part of the tier below the upper triangle. Points clustered around the core of the greater triangular represent rhodoliths for which two of the three dimensions are close in value. A few points fall into the right side of the plot, which signifies a tendency toward development of more ellipsoidal or spindle-shaped rhodoliths. The composite slope of points suggests a slightly diagonal trend, although the majority of points cluster vertically within the two central boxes. There is a total absence of points within the left side of the triangular diagram. At 60%, a higher percentage of rhodoliths sampled from the Holocene deposit at Praia das Gatas fall into the upper-most triangle of the Sneed-Folk diagram (Figure 7b), whereas a clear minority appear in the next tier directly below. As found in the Pleistocene sample (Figure 7a), very few rhodoliths from the unconsolidated Holocene deposit appear in the right side of the plot, and there are no points at all in the lower left part of the diagram (Figure 7b). The overall cluster of loose rhodoliths from the Holocene deposit is more diffuse than found in the Pleistocene sample, but also indicates a slightly diagonal tendency toward ellipsoidal shapes.

Comparative Variations in Rhodolith Shape
Raw data on rhodolith size collected from the sample sites at Baía da Parda and Praia das Gatas are listed in Appendix A (Tables A1 and A2). With regard to shape, individual points representing individual rhodoliths are plotted on a set of Sneed-Folk triangular diagrams (Figure 7a,b). The spread of points compared in the two plots demonstrates a high degree of similarity in the variation of shapes between those localities. Points that fall into the top triangle represent the most spherical in shape. Thirty percent of rhodoliths sampled from the Pleistocene limestone fall into this category (Figure 7a), whereas the majority of points fall within the central part of the tier below the upper triangle. Points clustered around the core of the greater triangular represent rhodoliths for which two of the three dimensions are close in value. A few points fall into the right side of the plot, which signifies a tendency toward development of more ellipsoidal or spindle-shaped rhodoliths. The composite slope of points suggests a slightly diagonal trend, although the majority of points cluster vertically within the two central boxes. There is a total absence of points within the left side of the triangular diagram. At 60%, a higher percentage of rhodoliths sampled from the Holocene deposit at Praia das Gatas fall into the upper-most triangle of the Sneed-Folk diagram (Figure 7b), whereas a clear minority appear in the next tier directly below. As found in the Pleistocene sample (Figure 7a), very few rhodoliths from the unconsolidated Holocene deposit appear in the right side of the plot, and there are no points at all in the lower left part of the diagram (Figure 7b). The overall cluster of loose rhodoliths from the Holocene deposit is more diffuse than found in the Pleistocene sample, but also indicates a slightly diagonal tendency toward ellipsoidal shapes.

Comparative Variations in Rhodolith Size
Rhodolith size is conveniently plotted on bar graphs as a function of maximum and intermediate diameters drawn from the original data (Tables A1 and A2). The results for field samples from the Pleistocene of Baía da Parda (Figure 8(a,b)) and Holocene of Praia das Gatas (Figure 8(c,d)) are related according to size at 1.5-cm intervals, within pebble and cobble size rhodoliths. By convention, pebbles are defined as exhibiting a range in diameter between 0.4 and 6.4 mm, whereas cobbles belong to the next larger category between 6.4 and 256 mm with boulders following thereafter [21]. The highest proportion of rhodoliths registered according to long axis (Figure 8(a)) falls into the interval with a diameter between 6.6 and 8.0 mm signifying small cobbles; however, data derived from the intermediate axes of the same rhodoliths indicates a downward shift in the largest proportion to a category in the range of large pebbles. On the whole, the size of the sampled Pleistocene rhodoliths may be characterized as bordering between large pebbles and small cobbles, but includes outliers with a few exceptionally large specimens attaining diameters as much as much as 13 and 15 cm. In contrast, the highest proportion of rhodoliths from the Holocene deposit based on the long axis (Figure 8(c)) falls below the threshold for cobbles and that proportion is even greater when measurements from the intermediate axes are taken into account (Figure 8(d)). Variations in samples from the two collections are most striking when contrasted with respect to measurements derived from the long axes of the respective rhodolith deposits (Figure 8(a,c)). More than two-thirds of the Pleistocene rhodoliths qualify as small cobbles, whereas less than half of those from the Holocene deposit meet the same standard. On the whole, the sample drawn from the Holocene deposit reflects a somewhat smaller range in size with no exceptionally large specimens greater than 11 cm in diameter.

Comparative Variations in Rhodolith Size
Rhodolith size is conveniently plotted on bar graphs as a function of maximum and intermediate diameters drawn from the original data (Tables A1 and A2). The results for field samples from the Pleistocene of Baía da Parda (Figure 8a,b) and Holocene of Praia das Gatas (Figure 8c,d) are related according to size at 1.5-cm intervals, within pebble and cobble size rhodoliths. By convention, pebbles are defined as exhibiting a range in diameter between 0.4 and 6.4 mm, whereas cobbles belong to the next larger category between 6.4 and 256 mm with boulders following thereafter [21]. The highest proportion of rhodoliths registered according to long axis (Figure 8a) falls into the interval with a diameter between 6.6 and 8.0 mm signifying small cobbles; however, data derived from the intermediate axes of the same rhodoliths indicates a downward shift in the largest proportion to a category in the range of large pebbles. On the whole, the size of the sampled Pleistocene rhodoliths may be characterized as bordering between large pebbles and small cobbles, but includes outliers with a few exceptionally large specimens attaining diameters as much as much as 13 and 15 cm. In contrast, the highest proportion of rhodoliths from the Holocene deposit based on the long axis ( Figure 8c) falls below the threshold for cobbles and that proportion is even greater when measurements from the intermediate axes are taken into account (Figure 8d). Variations in samples from the two collections are most striking when contrasted with respect to measurements derived from the long axes of the respective rhodolith deposits (Figure 8a,c). More than two-thirds of the Pleistocene rhodoliths qualify as small cobbles, whereas less than half of those from the Holocene deposit meet the same standard. On the whole, the sample drawn from the Holocene deposit reflects a somewhat smaller range in size with no exceptionally large specimens greater than 11 cm in diameter.

Rhodolith Nucleation
Among the sampled rhodoliths from the Pleistocene limestone on Baía da Parda, a single rhodolith was found to be nucleated around a basaltic pebble and another was found to have formed around a visible shell fragment (Table A1). No evidence of nucleation around basalt pebbles or recognizable shell fragments was detected among rhodoliths sampled from the Holocene deposit at Praia das Gatas (Table A2). Nucleation by coralline red algae is required before the rhodolith begins expansion in semi-spherical form, but may initiate around a sand particle not readily detected when the structure is broken open. Except for clear examples of nucleation around a shell fragment or basalt pebble, 96% of sampled rhodoliths from the two study sites initiated growth around small particles.

Storm Intensity as Function of Estimated Rhodolith Volume and Weight
The mean value of rhodolith dimensions (long, intermediate, and short axes), calculated volume and estimated weight from the two study samples collected from the exposed east coast of Sal Island are summarized in Table 1, allowing for direct comparisons across all categories.

Rhodolith Nucleation
Among the sampled rhodoliths from the Pleistocene limestone on Baía da Parda, a single rhodolith was found to be nucleated around a basaltic pebble and another was found to have formed around a visible shell fragment (Table A1). No evidence of nucleation around basalt pebbles or recognizable shell fragments was detected among rhodoliths sampled from the Holocene deposit at Praia das Gatas (Table A2). Nucleation by coralline red algae is required before the rhodolith begins expansion in semi-spherical form, but may initiate around a sand particle not readily detected when the structure is broken open. Except for clear examples of nucleation around a shell fragment or basalt pebble, 96% of sampled rhodoliths from the two study sites initiated growth around small particles.

Storm Intensity as Function of Estimated Rhodolith Volume and Weight
The mean value of rhodolith dimensions (long, intermediate, and short axes), calculated volume and estimated weight from the two study samples collected from the exposed east coast of Sal Island are summarized in Table 1, allowing for direct comparisons across all categories. In contrast, recorded values for the single, largest rhodolith in size (long, intermediate, and short axes), calculated volume and estimated weight from the same two samples are featured in Table 2. A comparison between Tables 1 and 2 brings into question the level of wave energy necessary to shift the largest rhodolith in each sample from a more offshore position to a final resting place in a more onshore position. Indeed, the largest single rhodolith drawn from the Pleistocene sample has a calculated volume 4.25 times as great as the average rhodolith from the same sample and an estimated weight 4.20 times as great. By comparison, the largest single rhodolith drawn from the Holocene sample has a calculated volume 3.25 times as great as the average rhodolith from the same sample and an estimated weight 3.3 times as great. Accumulation of rhodoliths in such evident concentrations are peculiar to each locality and dependent on a depositional model that best fits those circumstances, respectively.

Pleistocene and Holocene Depositional Models
The processes responsible for accumulation of Pleistocene and Holocene rhodoliths on the windward shores of Sal Island share certain aspects. The most striking is the relationship of organic cobble-size nodules left adjacent to or directly on top of basaltic cobbles (Figures 3 and 5). It is notable that broken rhodoliths from both deposits rarely if ever feature a basaltic core, and when that style of nucleation is evident, the basalt core is smaller than available cobbles in the underlying substrate. Another shared factor is the preponderance of highly spherical rhodoliths capable of rolling in any direction under the influence of water turbulence (Figure 7a,b). The third and most distinctive aspect is the degree to which individual rhodoliths from Pleistocene and Holocene deposits are worn as a consequence of contact abrasion. It is instructive to compare the Sal Island rhodoliths with those from the supratidal deposit on Maio Island, found to cover an area of 27,000 m 2 with roughly 25 million rhodoliths characterized by highly spherical shapes between 10 cm to 20 cm in diameter [17]. The primary difference registered by the Maio deposit is its inclusion of individual nodules with lumpy surfaces that are relatively fresh, although the overall concentration of rhodoliths per unit surface area is similar to deposits on Sal Island. Like the Sal rhodoliths, those from Maio Island lack basalt cores. A unique aspect of the Maio deposit features a sample nucleated around a worn fragment of ceramic shard, which places the deposit in the context of historical time. The hallmark of the Maio deposit appears to be a single storm event of hurricane or near-hurricane strength that swept large rhodoliths ashore with minimal degradation of surface features.
By observation of tidal flux, the loose rhodoliths at Praia das Gatas signify a massive supratidal deposit out of range of normal tides (Figures 5 and 6). Rhodoliths from the better exposed Pleistocene deposit at Baía da Parda (Figures 2 and 3) and more patchy exposure at Praia das Gatas (Figure 4a,b) differ only in the amount of carbonate matrix present. Pockets of matrix at Baía da Parda are interpreted as lime mud produced as a consequence of rhodolith abrasion. Another difference between the Pleistocene and Holocene deposits is the smaller average size range of the rhodoliths in the Holocene sample (Table 1). Because few of the Pleistocene rhodoliths and none of the Holocene rhodoliths possess a basalt core, they are sure to have accreted in an offshore setting where basalt pebbles or cobbles were mostly unavailable for nucleation. The largest rhodoliths from both deposits imply accretion over a long-time span in deeper waters normally unaffected by storms. The Pleistocene deposit cannot be interpreted as a supratidal deposit, because prior to cementation, the lime mud was more likely to be removed by the back swash of retreating waves. All factors considered, the depositional model for Pleistocene rhodoliths on the east shore points to a shallow, subtidal environment (similar to today) in a zone of rapid accumulation where surface degradation of rhodoliths was intense for a protracted period of time causing algal growth to cease. In terminology commonly applied to fossil deposits, the Pleistocene rhodoliths are attributed to a death assemblage (thanatocoenoses) marked by shoreward transport from a deeper-water habitat where they previously thrived undisturbed for as much as 100 years based on established growth rates [9,10].

Direction of Storm Impact
According to conch divers on Sal Island interviewed in March 2013 [11], large rhodoliths commonly occur as far offshore as three or more kilometers at a water depth around 40 m. Tropical cyclonic storms tracking northward between the African mainland and Sal Island would move in a direction contrary to the Northeast Trade Winds for a brief time, causing ocean swells greater than an amplitude of 3.5 m normally associated with the trade winds. The occurrence of such storms is rare and in terms of human record-keeping might be regarded as a once-in-a-century event, which corresponds to the age estimates of the longest-lived rhodoliths. Transport from deeper, offshore waters to shallow coastal waters over a distance of a few kilometers is unlikely to happen during a single super storm, but could bring larger rhodoliths to an intermediate position where lesser storms happening with greater frequency would roll those rhodoliths closer to shore. Once established as a deposit in shallower waters, these large rhodoliths were further subjected to the normal coastal surf generated by persistent Northeast Trade Winds, thereby causing post-mortem abrasion.
A different model is employed to explain how slightly smaller rhodoliths in the Holocene deposit at Praia das Gatas were recycled from the adjacent Pleistocene deposit (Figure 4a,b) as a result of selective remobilization by wave conditions only moderately more energetic than those created by the Northeast Trade Winds. Applications of mathematical equations from Nott (2003) [30] are often used to estimate the required wave height for the detachment, lifting, and removal of boulder-size clasts from joint-bound layers on rocky shores under storm attack, but this methodology is unsuited for application to fossil rhodoliths cemented in place by carbonate matrix as exposed at Baía da Parda. Instead, wave action related to squalls on a repeating, multi-decadal basis could be responsible for forming Holocene rhodolith deposits all along the east shores of Sal Island, creating supratidal accumulations. Individual rhodoliths are certain to have worked free of enclosing Pleistocene strata that protrude over the intertidal zone on Baía da Parda (Figure 2), but those particular rhodoliths were unlikely to be exported outside the bay and transported northward to Praia das Gatas against longshore currents in concert with the prevailing Northeast Trade Winds. The intertidal Pleistocene deposit of rhodoliths at Praia das Gatas would have been sufficient in scope to supply all the reworked rhodoliths in the adjacent supratidal deposit.

Inference from Historical Storms
Storm activity in the Cabo Verde Archipelago at a latitude around 15 • N is notorious for the initiation of hurricanes that typically form under conditions of unusually warm surface-water (26-27 • C) leading to cyclonal rotation of the overlying air mass. Storm tracks traced by these disturbances have been the subject of scrutiny ever since the start of weather prediction for busy shipping lanes [31] and an average September during the North Atlantic hurricane season may entail as many as a dozen disturbances that usually dissipate in the more central, extra-tropical regions of that ocean. Only a few achieve wind speeds of 119 km/h equivalent to a Category 1 hurricane or more according to the Saffir-Simpson Hurricane Wind Scale. The 2020 North Atlantic hurricane season got off to an unusual pre-season start in May with two tropical storms registering wind speeds between 63 and 118 km/h (Arthur and Bertha) and a record number of actual hurricanes (Hanna, Isaias, Laura, Marco, Nana, and Sally) crossed the Atlantic Ocean to North and Central America by September [32]. On 7 September 2020, Tropical Storm Rene brought strong winds and heavy rain to the northern islands of the Cabo Verde Archipelago, including Sal Island [32]. The disturbance subsequently intensified to a Category 1 hurricane while tracking northwest into the central North Atlantic. Subsequent systems have entered the Atlantic Ocean during this busy season to include tropical storms Teddy and Vicky. The historical record shows that very few hurricanes depart from the Cabo Verde region on a northerly track to reach the Madeira Islands and even mainland Portugal. As recently as October 2005, however, the passage of Hurricane Vince crossed Madeira as a Category 1 event and continued onward to the Iberian mainland as a weakened tropical depression. The rarity of such an event compares to a historical event of comparable intensity that struck the same region in October 1842, 163 years earlier [33]. Admittedly uncommon from a human perspective, development of Category 1 hurricanes at a frequency of one every 150 years over the last 10,000 years is potentially sufficient to disrupt the growth of large, semi-spherical rhodoliths accustomed to a habitat at water depths normally untouched by storms. Rare events in the context of a human lifetime take on a greater significance in the context of geologic time, during which even the most infrequent repetitions amount to substantial impacts on the seabed around islands in the northeast Atlantic Ocean.

Sal Island Rhodolith Deposits Compared to Others in the North Atlantic
Rhodoliths degraded as organic pebbles and cobbles are widely distributed around islands in the northeast Atlantic Ocean and can be correlated with aspects of coastal geomorphology pertinent to storms from the geologic past. The contrast between two examples of middle Miocene age from opposite shores on Ilhéu de Cima off Porto Santo in Madeira [13,14] are insightful. Both entail highly spherical rhodoliths with diameters between 10 and 20 cm. The basal rhodolith deposit at Cabeço das Laranjas has a maximum thickness of 2.6 m and occupies a rocky shelf with tall sea stacks preserved in three dimensions [13]. Some rhodoliths reveal small basalt cores, but most are without. The initial deposit is followed by three more rhodolith layers, all thin and all elevated above the unconformity surface indicative of a coastal transgression. The thick basal layer is interpreted as a major storm deposit transported shoreward from deeper water and left in a rocky-shore setting among sea stacks that acted as a baffle trapping the rhodoliths in place [13]. At Pedra de Água on the island's opposite shore, a thin rhodlolith deposit sits on an unconformity among basalt mounds with low relief occupied by corals in growth position [14]. The surrounding rhodoliths are also interpreted as part of a biocenosis preserved under life conditions. In this case, all feature a large basalt core as much as 6.5 cm in diameter and a relatively thin algal crust seldom more than 1.5 cm in thickness. The example demonstrates that large rhodoliths can form in shallow water equally conducive to the development of coral bioherms, but accreted around igneous cobbles eroded under nearshore conditions. The thin crust developed by such rhodoliths belies a shorter time span of growth compared to rhodoliths with the same diameter lacking a rock core.
Another example from the upper Miocene at Castilhano on São Nicolau in the Cabo Verde Archipelago [12] relates to smaller rhodoliths commonly 3.3 cm in diameter typically nucleated around basalt pebbles and deposited in a rocky-shore setting on an irregular basalt unconformity. A similar case may be cited from the lower Pliocene at Malbusca on Santa Maria Island in the Azores [15], where small rhodoliths from 2 to 3 cm in diameter occur with nuclei dominated by shell fragments and only rarely by basalt particles. These are found in deposits more than a meter thick entrained in surge channels eroded in pillow basalt. Based on the nucleation pattern of the Malbusca rhodoliths with few basalt cores, they are likely to have dwelled a bit father offshore than the São Nicolau rhodoliths but still in relatively shallow water. Small rhodoliths with a diameter between 4 and 5 cm are also reported in accumulations as much as a half-meter in thickness from the Lower Miocene of the Lagos-Portimão Formation in southern Portugal [16]. The exposed contact at Lagos between Miocene and Cretaceous strata is a paraconformity between limestone units. It is not surprising that the Miocene rhodoliths lack a rock core. Moreover, rhodolith beds are incorporated into the Lagos-Portimão Formation in a cyclic fashion that suggests episodic storm events. Overall, the depositional pattern shared by smaller fossil rhodoliths suggests a growth history cut short by the recurrence of lesser storms that raked sea bottoms closer to shore with greater frequency than hurricanes more capable of stirring deeper waters. Clearly, the largest deep-water rhodoliths recognized today from the Abrolhos Bank off Brazil [4] have adequate time to grow in size by accretion from smaller, incipient rhodoliths.
Hurricane incidence in the island groups of the northeast Atlantic Ocean may appear infrequent in terms of today's historical records, but a higher incidence during Miocene and Pliocene times was sure to have occurred with the Inter-Tropical Convergence Zone (ITCZ) displaced from the equator during the boreal winter to a more northern latitude beyond 20 • experienced today during the boreal summer's monsoonal climate. Analysis of Ocean-Drilling-Program (ODP) Site 659 located east of Sal Island [34], indicates a higher level of terrigenous influx off the Cabo Verde plateau between 6.2 and 4.9 million years ago than typical today. The inference is that elevated levels of storminess led to enhanced rainfall over the islands on the plateau and this factor may account for the occurrence of the extensive Miocene rhodolith deposits on São Nicolau Island [12], as well as other places like Madeira [13] and southern Portugal [16]. Such differences in the organization of the ITCZ in a more northerly position are marshalled to account for changes registered as far north as Malta in the Mediterranean Sea following the Messinian Salinity Crisis [34,35]. It remains to be seen if the pattern of ongoing climate change that contributes to increasing storminess in the Atlantic Ocean may include wayward storm tracks affecting the Canary, Madeira, and Azorean islands.

Conclusions
Rhodolith samples from Pleistocene and Holocene sites on the windward eastern shore of Sal Island in the Cabo Verde Archipelago register variations in size, volume, and weight that offer insights on the effectiveness of storm waves in the final configuration of those deposits, as well as trends from older Pliocene and Miocene deposits from other Atlantic archipelagos.

•
Rhodoliths formed by the accretion of coralline red algae in nodules taking on a semi-spherical shape are moderately large at 15 cm in diameter from the upper Pleistocene at Bahía da Parda and 10 cm in diameter from the Holocene deposit at Praia das Gatas, respectively. • Fossil rhodoliths as large as 15 cm in diameter and lacking rock cores appear to have survived for a long time, potentially as much as 100 years through continuous development in deeper, offshore waters where disruption from storm events was seldom experienced, but were conveyed by storms to a position superimposed on basalt cobbles.

•
Upper Pleistocene deposits at Praia das Gatas and Baía da Parda correspond to death assemblages dominated by rhodoliths with minor faunal elements that include fragmented corals and marine gastropods. The worn outer surface of densely branched rhodoliths typical of species in the genus Lithothamnium points to a shallow, subtidal setting in which very fine biogenic carbonate was produced as a result of contact abrasion to provide the matrix in which the rhodoliths are cemented.

•
The extensive supratidal Holocene deposit of loose rhodoliths found at Praia das Gatas was recycled in its entirety from the adjacent Pleistocene deposit at the same locality as a result of storm activity that created surf conditions more energetic than generated by the constant Northeast Trade Winds on the exposed shore of Sal Island. Rhodoliths in the Holocene deposit are somewhat smaller, overall, but with surfaces subject to additional contact abrasion.

•
Analysis of the Sal Island rhodoliths compared with older Miocene and Pliocene rhodolith deposits elsewhere in the Cabo Verde Archipelago as well as other islands in Madeira, the Azores, and mainland Portugal where deposition under conditions of increased storminess is equated with a more northern position of the ITCZ. Table A1. Quantification of rhodolith size, volume, and estimated weight from Pleistocene samples at Baía da Parda. The density magnesium-rich limestone at 2.6 gm/cm 3 is applied uniformly in order to estimate the weight of each sample based on the radius measured against the intermediate axis.