Next Article in Journal
Fluid-Rock Interactions in a Paleo-Geothermal Reservoir (Noble Hills Granite, California, USA). Part 1: Granite Pervasive Alteration Processes away from Fracture Zones
Next Article in Special Issue
A Long-Term Record of Quaternary Facies Patterns and Palaeonvironmental Trends from the Po Plain (NE Italy) as Revealed by Bio-Sedimentary Data
Previous Article in Journal
Cloud Masks Derived from the Mars Daily Global Maps and an Application to the Tropical Cloud Belt on Mars
Previous Article in Special Issue
Distribution of Holocene Marine Mud and Its Relation to Damage from the 1923 Earthquake Disaster in the Tokyo Metropolitan Area, Japan
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Holocene Sea-Level Changes in Southern Brazil Based on High-Resolution Radar Stratigraphy

by
Eduardo Guimarães Barboza
1,2,*,
Sergio Rebello Dillenburg
1,
Matias do Nascimento Ritter
2,
Rodolfo José Angulo
3,
Anderson Biancini da Silva
4,
Maria Luiza Correaa da Camara Rosa
1,
Felipe Caron
2 and
Maria Cristina de Souza
3
1
Centro de Estudos de Geologia Costeira e Oceânica, Instituto de Geociências, Universidade Federal do Rio Grande do Sul, Porto Alegre 91509-900, RS, Brazil
2
Centro de Estudos Costeiros, Limnológicos e Marinhos, Campus do Litoral, Universidade Federal do Rio Grande do Sul, Imbé 95625-000, RS, Brazil
3
Laboratório de Estudos Costeiros, Departamento de Geologia, Universidade Federal do Paraná, Curitiba 81530-000, PR, Brazil
4
Gerência de Licenciamento Ambiental de Atividades Estratégicas, Instituto do Meio Ambiente de Santa Catarina, Florianópolis 88020-300, SC, Brazil
*
Author to whom correspondence should be addressed.
Submission received: 11 June 2021 / Revised: 24 July 2021 / Accepted: 29 July 2021 / Published: 31 July 2021

Abstract

:
This paper focuses on high-resolution coastal stratigraphy data, which were revealed by the Ground Penetrating Radar (GPR) system. Surveys performed with GPR on the surface of prograded barriers reveal patterns of reflections that allow the interpretation of the geometry and stratigraphy of coastal deposits in a continuous mode. At the Curumim prograded barrier in southern Brazil (29°30′ S–49°53′ W), a two-dimensional transverse GPR survey revealed, with high precision, a serial of contacts between aeolian deposits of relict foredunes and relict beach deposits that have a strong correlation with sea level. In a 4 km GPR profile, a total of 24 of these contacts were identified. The high accurate spatial positioning of the contacts combined with Optical Stimulated Luminescence dating resulted in the first confident sea-level curve that tells the history of sea-level changes during the last 7 ka on the southernmost sector of the Brazilian coast. The curve shows that sea-level was still rising before 6 ka BP, with a maximum level of 1.9 m reached close to 5 ka BP; after that, sea-level started to falling slowly until around 4 ka BP when fall accelerated.

1. Introduction

The pursuit of a confident sea-level curve representing oscillations of the sea level during the middle and late Holocene has been a constant goal in the last decades. In Brazil, since 1979, when the first most detailed curve was proposed [1], and later contested [2], two groups have debated the existence of high-frequency oscillations (2–3 m), operating in the scale of centuries (500–600 years), that could have occurred during the overall sea-level fall established after a maximum level of few meters was reached between 6–5 ka BP [3,4]. Since then, high-frequency oscillations have lost credibility. Moreover, one reason for that is that strong (convincing) dated coastal geomorphological records were never formally presented, confirming such oscillations. From all indicators of paleo sea-levels used by those authors, such as marine (beach), lagoon and mangrove deposits, vermetid incrustations and coral reefs, and even archaeological indicators such as “Sambaquis”, the vermetid incrustations got more credibility and confidence as a sea-level indicator [4].
In South Brazil, to the south of the Santa Marta Cape (28°35′ S-48°49′ W–Figure 1), a formal and confident curve was never proposed due to the lack of confident indicators of paleo sea level along with this almost entirely sandy coast. For instance, vermetid incrustations occur fixed on coastal rocks, which do not occur to the south of the Santa Marta Cape, except for Torres (29°20′ S–49°44′ W–Figure 1) beach, where a rock promontory occurs, which has no occurrence of fossil vermetid. A very general curve was proposed but not supported by any of the indicators mentioned above of paleo sea levels [5]. Later, a lagoonal terrace positioned at around +3 m asl was related to the maximum sea level reached between 6–5 ka BP [6], and the contact between aeolian and beach deposits in cores of the prograded Curumim barrier (29°38′ S–Figure 1) was correlated with paleo sea levels [7]. These contacts were also identified in a GPR record of a prograded barrier located 50 km to the north of Curumim barrier at the Passo de Torres (29°15′ S–49°42′ W–Figure 1) and correlated with paleo sea levels [8]. Only recently, in South Brazil, better attention was given to this paleo sea-level indicator (aeolian/beach contact), along the prograded Cassino barrier (32°10′ S-52°11′ W–Figure 1), where a maximum sea-level of 2.0 m reached at 6–5 ka BP and after that followed by an overall sea-level fall was determined [9]. This indicator of paleo sea level was detected in eight cores along a 15 km long transect on the prograded barrier of Cassino. The contact was defined with a precision of 0.5 m due to the remarkably high contrast of compaction between relict aeolian sand deposits of foredunes and backshore sands occurring underneath. The chronology of barrier progradation was determined by a combination of radiocarbon (14C) and Optical Stimulated Luminescence (OSL) datings. The limitation at the Cassino barrier was the necessity of a high-cost drilling system and the consequent discontinuity of such records along a transect.
After the pioneer study that combined the dune-beach facies boundary and the dating of the dune base as an indicator of paleo sea-level [10], some new studies have used the relict top of the berm facies [11] or the raised beach strata [12] as indicators of paleo sea levels. More recently, converging methods of paleo sea-level reconstructions were presented [13,14], based on the contact between relict foredunes and backshore deposits along with prograded barriers. The advantage of these methods is the more precise identification of the foredune/backshore deposits contact in GPR records. The continuity of the contact in GPR records makes it possible to determine a great number of spatial positions of such contacts. Combined with OSL datings (informing the chronology of the development of a prograded barrier), the method allows the reconstruction of sea-level behavior during the period of barrier progradation.
In a sandy coastal region as southern Brazil, showing two large prograded barriers (Cassino and Curumim), we choose the Curumim barrier to extract relict records of paleo sea levels. By combining high-quality geophysical records obtained in GPR surveys, OSL dating, and determinations of altitudes by a Global Navigation Satellite System (GNSS) system, the reconstruction of sea-level behavior during 4 km progradation of the Curumim barrier was made possible.
Figure 1. Location of the Curumim barrier study area (dashed red rectangle) showing the general geology of the coastal plain of southern Brazil–datum WGS84 (modified after [15]).
Figure 1. Location of the Curumim barrier study area (dashed red rectangle) showing the general geology of the coastal plain of southern Brazil–datum WGS84 (modified after [15]).
Geosciences 11 00326 g001

2. Regional Setting

Curumim gives name to the prograded (or regressive) barrier that was formed in a very gentle embayment between Tramandaí and Torres [15] (Figure 1). The coastal plain where the Curumim barrier occurs is characterized by coastal sandy barriers and lagoons formed during the Quaternary [16], and its location in an intra-plate position suggests a tectonically stable condition [17]. There is no indication that basin subsidence and sediment compaction were relevant during the Holocene [18]. Right at Curumim, the barrier shows its widest dimension (5 km) due to barrier progradation during the last 7.3 ka [19]. The prograded plain of the barrier is completely and dominantly covered by aeolian deposits of transgressive sand sheets (TSS), and rare transgressive dunefields. At least 12 phases of TSS were formed during barrier progradation, and their mechanism of formation seems to be related to cyclic processes of coastal erosion that have occurred during barrier progradation [20]. Each phase is connected with the destruction of at least one foredune ridge by storm waves. The Curumim barrier belongs to the barrier-lagoon system IV, as referenced in southern Brazil, formed during the Holocene [21]. This barrier-lagoon system corresponds to a high-frequency depositional sequence [22,23]. The modern beach of the Curumim barrier is formed by fine sand, showing dissipative to intermediate morphodynamic stages, and submitted to a microtidal regime with a mean range of only 0.3 m, and a maximum high spring tide of 0.53 m [24].
Consequently, sediment transport and deposition along the open coast at the Curumim barrier is primarily dominated by moderate wave energy with a mean significant wave height (Hs) of 1.25 m [25]. During autumn and winter storms (April–July) the wave height may frequently exceed 2.0 m, and the sea-level can surge up to 1.3 m [26]. Expected extreme events can increase sea-level up to 3.2 m [24]. The climate is subtropical, constantly humid (Cfa), without a dry season, with hot summer [27]. Rainfall ranges from 1600 to 1900 mm and is evenly distributed throughout the year. NE winds predominate along the entire coast but combined with the secondary winds from the S, a resultant drift direction to WNW is produced [28].

Sea-Level History

The postglacial sea-level history of the study area extends back approximately 17.5 ka BP when the sea level was approximately 120 to 130 m lower [29]. After that time, the sea level rose at an average rate of 1.2 cm/year and reached the present level for the first time at around 6 ka BP (this work). As mentioned before, there are no reliable data on sea-level behavior during the middle to late Holocene time along this coast. However, sea-level curves for areas to the north indicate that at the culmination of the Postglacial Marine Transgression (PMT) at around 5 ka BP, the sea level was approximately 1–3 m above its present level, after which it slowly fell [4] (Figure 2A). Relict vermetid incrustations records of a coastal area located 200 km to the north of the Curumim barrier indicate a maximum sea level of +2.1 m at around 5 ka BP [30]. To the south (~700 km), in Uruguay, a sea-level envelope indicates a maximum level varying between +3.2 and +2.1 m occurred at around 6 ka BP, followed by an overall fall [31] (Figure 2B). Also, from Uruguay comes the closest record of modern sea-level rise, which indicates a rate of 4.7 mm/years, from 1961 to 2014 [32]. Records of middle-late Holocene punctual sea-level positions obtained in studies along of the coastal plain of the Rio Grande do Sul [6,9], fit the envelope curve of Figure 2A [4].

3. Materials and Methods

3.1. Ground Penetrating Radar (GPR) Profile

The two-dimensional high-resolution stratigraphic GPR survey was performed on the northern coastal plain of the Rio Grande do Sul State (Figure 1) at the Arroio do Sal (29°30′ S-49°53′ W). The survey was carried out using a four 4WD vehicle, in which the acquisition and positioning units were installed. The utilized GPR system was composed of a GSSI™ (Geophysical Survey Systems, Inc., Nashua, NH, USA) SIR-3000 data collector for monostatic ground-coupled antenna 200 MHz. The GPR profile was collected using the Common Offset method with a two-way-travel time (TWTT) range of 450 ns, which allows for penetration up to 20 m of depth. A vertical IIR low pass (500 MHz), vertical IIR high pass (60 MHz), stacking (3), and gain filters were applied during the time of data acquisition.
The data were post-processed with the Radan™ and Prism2® software packages after background removal and the application of bandpass frequency filters, Ormsby bandpass filter, gain equalization, zero-point, topographic corrections, and time-to-depth conversion. The trace analysis was done, and a dielectric constant of 10 for wet sand was used to convert travel time to depth, which represents a velocity of 0.09 m/ns [33,34,35]. This constant was validated using lithological data obtained from SPT drill holes [9,36,37].
The GPR profile interpretation was based on the seismostratigraphy method [38] adapted for GPR [39]. The method is based on reflections termination (onlap, downlap, toplap, and truncations), geometry, and pattern of reflections [40,41,42,43]. These properties of the reflectors together with their continuity characteristics and signal amplitude were the criteria to define the point of contact between the base of the foredune and the beach (backshore) deposits with better precision. Records of short-term extreme events of sea-level rise are easily identified as erosional surfaces truncating foredune internal reflections.

3.2. Global Navigation Satellite System (GNSS) Profile

The GPR profile was topographically corrected using GNSS post-processed elevation data points collected along the profile lines at interval times of 1 s. These data were acquired using a GNSS Trimble® (Sunnyvale, CA, USA) ProXRT (datum: WGS84) and analyzed in a Geographic Information System (GIS). The data were post-processed with the Trimble® Pathfinder® Office software (Sunnyvale, CA, USA), resulting in a horizontal and a vertical precision of 0.1 and 0.3 m, respectively.

3.3. OSL Ages

Four samples from aeolian deposits of TSS, with each sample representing one phase of TSS formation, were collected by horizontal insertion of a pipe into the sidewalls cut for road construction, following the classical protocol to avoid any exposure to sunlight. They were dated at the Luminescence Dating Laboratory, Victoria University of Wellington, NZ. OSL ages estimated were obtained on quartz grains, whose concentration followed standard procedures and included wet sieving to isolate the 90–200 µm grain size fraction followed by chemical treatments with hydrogen peroxide (10%) and hydrochloric acid (10%) to eliminate organic compounds and carbonates, respectively [44]. Quartz grains (2.58–2.65 g/cm3) were subsequently isolated from heavy minerals and feldspars using lithium polytungstate solutions. Quartz concentrates were then treated with hydrofluoric acid (40% HF) to eliminate the outer rind of grains affected by alpha radiation and remnant feldspars.
Equivalent doses were determined using the single-aliquot regenerative dose (SAR) protocol applied to multigrain aliquots of quartz [45]. Luminescence measurements were done using a standard Riso TL-DA15 measurement system, equipped with 3 Hoya U340 optical filters (stack thickness ~9 mm) to select the ultraviolet luminescence band centered around 340 nm. Stimulation was done by blue diodes, α-irradiations were done on the Riso TL-DA15 90Sr,90Y β-irradiator, calibrated against the Risø National Laboratory, Denmark, to about 3% accuracy. The palaeodoses were estimated by use of the Single Aliquot Regenerative method (SAR). In the SAR method, a number of aliquots are subjected to a repetitive cycle of irradiation, preheat, and measurement. In the first cycle, the natural luminescence output is measured, in all following cycles an artificial dose is applied. Usually, four or five of these dose points are used to build the luminescence growth curve (β-induced luminescence intensity versus added dose) and bracket the natural luminescence output. This allows interpolation of the equivalent dose (the β-dose equivalent to the palaeodose). In order to correct for potential sensitivity changes from cycle to cycle, a test dose is applied between the cycles, preheated (‘cut heat’), and measured. For the samples reported here 16 aliquots were measured, preheat temperature was 260 °C for 10 s, cut heat was 220 °C for 10 s, and measurement time 40 s (which resets the luminescence signal to a negligible residual).
Dose rates were calculated through concentrations of uranium (U), thorium (Th), and potassium (K). The dry, ground, and homogenised soil samples were encapsulated in airtight perspex containers and stored for at least four weeks. This procedure minimizes the loss of the short-lived noble gas 222Rn and allows 226Ra to reach equilibrium with its daughters 214Pb and 214Bi. The samples were counted using high-resolution gamma spectrometry with a broad energy Ge detector for a minimum time of 24 h. The spectra were analysed using GENIE2000 software, and cosmic dose rates were calculated [46]. The dose rate calculation is based on the activity concentration of the nuclides 40K, 208Tl, 212Pb, 228Ac, 210Pb, 214Bi, 214Pb, 226Ra (Table 1).

3.4. Geostatistic

The first curve of Holocene paleo sea levels for southern Brazil (29°30′ S-49°53′ W) is based on 24 positions, in which the distance from shore (m) and altimetric values relative to the present mean sea-level (MSL) are indicated (Table 2). A modern point at (0, 0) was also included to indicate the present mean sea level. The positions were then modeled linearly using the function “lm” with a 5th-order polynomial employing the argument “poly”. The predicted intervals are based on 0.99 level confidence, which represents here the envelope curve. It means a 99% confidence that the following new observation (distance of shore) will fall within this range. The adjusted coefficient of determination R2 was determined with a significance level gamma threshold value of 0.05. All analyses were carried out in the open-source R language [47].

4. Results

The analysis and interpretation of the GPR record obtained along a 4.3 km long profile made possible the identification of different patterns of reflectors, which were separated as distinct radarfacies (Rdf) (Figure 3).
Rdf I presents irregular, non-continuous, and sub-parallel reflectors in a general segmented undulating pattern, commonly found at the upper shoreface. Rdf II shows a downlap termination over radarfacies I and a toplap termination with radarfacies III. The downlap and toplap terminations characterize an oblique-parallel progradational configuration (clinoforms) formed under a high-energy depositional condition, typical of the backshore/foreshore zone. Rdf III shows reflectors in downlap termination over the clinoforms of Rdf II for both directions (landwards and seawards), a marked characteristic of foredunes. The upper limit of this radarfacies is marked by a continuous and undulating high amplitude reflector. Rdf V has continuous, sub-parallel reflectors with a low-angle downlap termination in the landward direction, corresponding to the aeolian TSS that dominate at the barrier surface [48,49,50].
As demonstrated by many authors, the contact between foredune and backshore deposits is strongly correlated with sea-level position, e.g., [9,10,12,13,14,51]. On coasts such as the one of Curumim barrier, the low micro-tidal range (0.3 m) makes wave energy (maximum runup during regular storms) the main factor controlling the vertical distance of this contact in relation to the present sea-level position. On this southern coast of Brazil, this vertical distance ranges from 1.0 to 1.5 m [9]. At the Curumim barrier, the modern contact between foredune and backshore deposits can be identified by the downlap and toplap terminations of reflectors representing foredune and backshore deposits, respectively (see Figure 3 and Figure 4). Eleven (11) measurements made at the modern (active) dune-beach system of the Curumim barrier in 2003 and 2004 (unpublished data) indicated an average vertical distance of 1.5 ± 0.4 m between this contact and the mean sea level. This vertical distance is also found 50 km to the north of the Curumim barrier [8]. In contrast, on coasts such as Troia Peninsula (Portugal) (macro-tidal range) and at The Granites (South Australia) (moderate wave energy but macro-tidal range), this vertical distance is around 4.6 and 3.0 m, respectively [13,52]. Thus, it is clear the great influence of the tidal range on the altitude of the contact between foredune and backshore deposits. With the assumption that this vertical distance has remained almost the same during the Holocene evolution of prograded barriers (meaning relatively constant wave energy and tidal range), it is possible to access past sea-level positions by subtracting the vertical distance from the heights of the contact of relict foredune and backshore deposits identified in GPR records [53].
A total of 24 points precisely marking the contact of downlap and toplap terminations between foredune and backshore radarfacies were determined (Figure 3 and Figure 5). Figure 3 presents 250 m only of the 4.3 km long GPR profile performed at Arroio do Sal on the Curumim barrier (Figure 1 and Figure 5). The portions along of the GPR profile in Figure 5, where the contacts between the foredune and the backshore deposits were not identified, are related to places eroded by channelized waters.
Relative ages of these contacts were obtained in the form of four OSL ages (Table 1), whose projections are presented together with the GPR profile in Figure 5. Each analyzed sample provided the age of a starting point of a TSS, which in turn approximately corresponds to the location of a former foredune ridge. The four ages inform the chronology of barrier progradation and consequently inform the relative chronology of the formation of the 24 contacts between foredune and backshore radarfacies presented in Figure 5.
After subtracting the 1.5 m vertical distance (the vertical distance between the modern contact and the present position of sea-level–Figure 4B), 24 paleo sea levels were determined (Table 2).
The fitted sea-level variation curve, presented in Figure 6, shows a promising and statistically significant outcome (R2 = 0.966, p < 0.001). The predicated confidence intervals at 99% level (Figure 6; black dashed lines) agree with the envelopes [4,28] presented in Figure 2. It is noteworthy that the curve, notwithstanding be bidimensonally reached, using the distance from shore versus sea-level, encompasses an age-related tridimensional understanding, including age-estimates prior to the maximum sea-level reached at about 5 ka BP.

5. Discussion

The coastal region between the Santa Marta Cape (28°35′ S-48°49′ W–Figure 1) and La Coronilla (33°55′ S-53°30′ W–Figure 1) is 750 km long, representing the emergent portion of the Pelotas oceanic marginal basin, showing essentially Quaternary age deposits. In the last 7 ka, a coastal barrier was formed as a shore parallel geomorphological structure, continuous for 750 km [15]. Immediately to the north of the Santa Marta Cape, at Imbituba-Laguna (28°14′ S-48°39′ W), spatial positions of relict vermetid incrustations and corresponding 14C datings indicated that +2.1 m was the highest position reached by sea level during middle Holocene (6–5 ka BP), followed after that by an overall sea-level fall [29]. Later, a sea-level envelope was presented for the Imbituba-Laguna coastal region, indicating a maximum sea level varying from +3 to +1 m, reached at around 5 ka BP [4].
The precision of the spatial position of the 24 indicator points of past sea-level, determined in GPR records, is controlled by the operational basis (fundamentals) of the GNSS system and of the GPR system itself. In regard to the GPR system, precision is basically related to the velocity and frequency of the electrical signal in the unconsolidated sediments (sand-dominated in the case of this study). However, as the depths of the target contacts are less than 5 m and considering that contacts were also determined directly in 14 cores obtained from drill-holes [36], a much more real dielectric constant was applied (10, for wet sand), and thus precision seems to be remarkably high. The frequency of the GPR system (200 MHz) corresponds to a reflector resolution of 0.25 to 0.30 m [54,55]. A precision of 0.1 m (horizontal) and 0.3 m (vertical) was determined regarding the GNSS system. Thus, the altitudes of the target points are subjected to a maximum error of ±0.6 m. Indeed, the less precision value used to build the proposed sea-level curve is the vertical distance of 1.5 m, subtracted from the altitudes of the GPR reflectors representing the 24 relict contacts between foredune and backshore deposits identified along the prograded plain of the Curumim barrier. As mentioned above, this 1.5 m value was obtained from average measurements of 11 surveys performed at the modern dune-beach system of the Curumim barrier in 2003 and 2004, and from a survey made 50 km to the north of the Curumim barrier [8]. Moreover, as this vertical altitude is mainly controlled by wave energy acting on the Curumim barrier, due to the small average tidal range of 0.3 m, it was assumed that wave energy and consequently this average value of 1.5 m were kept constant during 7 ka of barrier progradation.
An exercise of validation of the sea-level curve obtained from the Curumim barrier was done by comparison with nearby data of sea-level changes during the Holocene. At Imbituba-Laguna (160 km to the north of Arroio do Sal), the maximum sea level of +2.1 m, around 5 ka BP [30] was obtained from fossil vermetid records whose precision is indicated by a standard deviation of ±1 m. This standard deviation results from variations of the altitude of the vermetid bioconstructions produced by variations of wave energy and tidal range [56,57]. The variations also result from the morphology of the rocky coast, the homologous reference level considered, and the preserved part of the bioconstruction [58].
Table 2 shows that at Arroio do Sal, between 5.99 and 5.03 ka BP, altitudes of sea-level positions varied between 1.9 and 0.4 m. Thus, independently of the margin of errors of the two methods of paleo sea-level determination, the maximum sea levels that occurred at Arroio do Sal (1.9 m) and at Imbituba-Laguna (2.1 m), at around 5 ka BP, are remarkably similar. To the south of Arroio do Sal, the nearest sea-level data are distant by approximately 700 km, embracing a large coastal area of the Atlantic Ocean (Argentina and Uruguay) and the La Plata River [31,59,60]. Those data are mainly based on sea-level indicators extracted from beach ridges and beach storm deposits, and they have been intensely debated in regard of their meaning as sea-level indicators for the last 7 ka of this coastal region [31,59,60]. One of the above studies interpreted a maximum level of +4 m at around 6 ka BP and two high-frequency oscillations of sea level that may have occurred at around 4.3 and 2.0 ka BP [59], while the other two studies interpreted a continuous sea-level curve, progressively falling to the present position, however, they disagree between each other by around 1 m in regard of the maximum level reached at around 6 ka BP [31,60]. Thus, despite the disagreement between researchers studying sea-level behavior at Argentina and Uruguay, some general similarities with the sea-level record obtained at Arroio do Sal could be highlighted, such as the maximum level of few meters reached around 6–5 ka BP and followed by an overall sea-level fall to the present position (not considering high-frequency oscillations). However, when comparing the Arroio do Sal sea-level data with results of Uruguay [31] only, similarities are worthy of note. For the Uruguayan coast a maximum sea level varying between +3.2 to +2.1 m (average +2.65 m), at around 6 ka BP, followed by a continuous and progressively fall [31].
This paper does not discuss the drivers of sea-level behavior in the last 7 ka. However, considering that the Curumim barrier is located in the central part of a marginal sedimentary basin (Pelotas Basin) and that, for instance, the Atlantic coast of Uruguay occurs at a structural high that corresponds to the southern limit of Pelotas Basin, some geological factors such as the differential compaction of sedimentary deposits, and tectonics could come to explain the difference between the maximum sea-levels reached at around 6–5 ka BP, recorded at Arroio do Sal and at Argentina and Uruguayan coasts. For instance, regarding tectonics is worth of mention that in a south direction, the eastern continental margin of South America gets closer to its counterpart tectonic-active western continental margin. These factors have not already deserved appropriate investigation by coastal researchers from Argentina, Uruguay, and Brazil. As a second hypothesis, the difference between the maximum sea levels reached at both regions (Uruguay and South Brazil) around 6–5 ka BP, could be related to the adjustment (correction) of the sea level indicators used to infer the paleo-sea levels at the Uruguayan coast. The altitude of the paleo-sea levels was determined by the subtraction of 1.5 m (considered the storm mean altitude), but the authors informed that the extreme storms recorded in the last 100 years reached 4.7 m [31]. As extreme events commonly produce records with a high potential of preservation, it is possible that the dated storm deposits correspond to storms that reached altitudes higher than 1.5 m, and therefore the inferred altitudes could be over-estimated.
Important to mention here that these extreme storm events would have reached southern Brazil as well. It was mentioned above an expected storm surge of up to 3.2 m [24]. But differently from the paleo sea-level indicators of the Argentine and Uruguayan coast, which are mostly beach storm deposits, the indicator analysed in the progradational record of the Curumim barrier was the depositional contact between the base of the foredune and the backshore deposits, which in turn was formed under a combination of fair-weather conditions and regular (not extreme) storms. Extreme storm events are recorded in GPR profiles, but they are easily identified as erosional surfaces truncating the foredune internal reflections.
The curve extracted from the Curumim barrier shows that sea-level was still rising before 6 ka BP, with a maximum level of 1.9 m reached close to 5 ka BP; after that, sea-level started to fall slowly until around 4 ka BP when fall was accelerated. It is not the goal of this paper to discuss the geological consequences of the sea-level behavior indicated by the curve of the Curumim barrier; however, it is important to raise the attention to the possibility that the fall acceleration after 4 ka BP may have resulted in changes over coastal systems of southern Brazil. For instance, palynology studies of Holocene lagoonal deposits of southern Brazil (28°01′ S-48°38′ W and 33°43′ S-53°21′ W) [61,62], respectively, found that the lagoonal systems have changed to a freshwater dominance after 4–3 ka BP. Other potential consequences deserve future investigation.
The global sea-level rise driven by global warming was not detected in the sedimentary records of the prograded barriers of the Pelotas Basin. Historical aerial pho-tographies of the last 40 years indicated that the Curumim barrier is still prograding at a rate of 3.0 to 3.4 m/year [63]. This means that a positive sediment balance may be overcoming the ongoing contemporaneous sea-level rise. As an exercise, by taking into account the rate of sea-level rise of 4.7 mm/year acting for 100 years, an elevation of 47 cm was expected of the contact between foredune and backshore radarfacies, at the GPR records close to the shoreline. The GPR record is blind along the last 150 m due to the saltwater incursion, which made it impossible to verify this probable elevation of the contact between the radarfacies.
At the end of the last century, high-resolution seismic reflection equipment was still used on Quaternary terrains aiming the obtention of continuous stratigraphic records, but at an extremely high cost and too much hard work. Fortunately, after the establishment of the GPR as a tool to investigate geological problems, following the early background fundamentals [64,65], the GPR devices quickly turned to portable and precise operating systems, furnishing conditions to run relatively low-cost quick surveys on Quaternary coastal terrains. Thus, it is gratefully and worthy of mention that after 30 years of pursuing the production of a sea-level curve telling the sea-level history for the last 7 ka of the southernmost coast of Brazil, this goal was finally reached here, although perfectible.

6. Conclusions

The two-dimensional GPR profile performed on the prograded barrier of Curumim revealed a detailed stratigraphy, showing patterns of reflectors that allowed the identification of the geometries of foredunes and beach (backshore) deposits. The downlap and toplap reflectors termination of these two deposits allowed to define with high precision their contacts, which in turn have a strong correlation with sea level.
With great confidence, a total of 24 points representing these contacts were identified in the GPR profile, allowing the construction of the first sea-level curve that represents the changes on sea-level during the last 7 ka, which occurred at the southernmost sector of the Brazilian coast.
The curve shows that sea level was still rising before 6 ka BP, with a maximum level of 1.9 m reached close to 5 ka BP; after that, sea level started to fall slowly until around 4 ka BP when fall was accelerated.
By providing the precise and continuous visualization of the well-defined contact between foredune and backshore deposits on prograded barriers, the GPR system, when combined with high-precision spatial positioning systems and with geochronological data such as OSL, made possible the construction of a curve remarkably similar to curves constructed, for instance, using biological indicators such as fossil vermetid, which have allowed the production of sea-level curves for most of the Brazilian coast in the last 30 years.

Supplementary Materials

The following are available online at https://0-www-mdpi-com.brum.beds.ac.uk/article/10.3390/geosciences11080326/s1, Figures S1–S4: GPR Section 1–4.

Author Contributions

Conceptualization, E.G.B. and S.R.D.; investigation, E.G.B., S.R.D., M.L.C.d.C.R. and A.B.d.S.; data processing, E.G.B., M.L.C.d.C.R., M.d.N.R., F.C.; writing—original draft preparation, E.G.B., S.R.D., M.d.N.R. and R.J.A.; writing—review and editing, E.G.B., S.R.D., M.C.d.S., R.J.A., M.L.C.d.C.R., A.B.d.S., F.C. and M.d.N.R. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by MCTI/CNPQ/Universal 14/2014, grant number 446060/2014-3.

Institutional Review Board Statement

The study was approved by the Institutional Review Board of Geoscience Institute of UFRGS (Project No.: 29099).

Informed Consent Statement

Not applicable.

Data Availability Statement

The study did not report any organized data sets.

Acknowledgments

E.G.B., S.R.D., R.J.A. and M.L.C.d.C.R. thank the Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq) for the provision of their research fellowships.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Martin, L.; Flexor, J.M.; Vilas-Boas, G.S.; Bittencourt, A.C.S.P.; Guimarães, M.M.M. Courbe de variation du niveau relative de la mer au cours des 7000 dernie’res anne’es sur un secteur homoge’ne du littoral bre’ silien (nord de Salvador—Bahia). In Proceedings of the International Symposium on Coastal Evolution in the Quaternary, São Paulo, Brazil, 11–18 September 1978; Suguio, K., Fairchild, T.R., Martin, L., Flexor, J.M., Eds.; USP: São Paulo, Brasil, 1979; pp. 264–295. [Google Scholar]
  2. Angulo, R.J.; Lessa, G. The Brazilian sea level curves: A critical review with emphasis on the curves from Paranaguá’ and Cananéia regions. Mar. Geol. 1997, 140, 141–166. [Google Scholar] [CrossRef] [Green Version]
  3. Martin, L.; Dominguez, J.M.L.; Bittencourt, A.C.S.P. Fluctuating Holocene sea levels is eastern and southeastern Brazil: Evidence from a multiple fossil and geometric indicators. J. Coast. Res. 2003, 19, 101–124. [Google Scholar]
  4. Angulo, R.J.; Lessa, G.C.; Souza, M.C. A critical review of Mid- to Late-Holocene sea-level fluctuations on the eastern Brazilian coastline. Quat. Sci. Rev. 2006, 25, 486–506. [Google Scholar] [CrossRef]
  5. Tomazelli, L.J.; Villwock, J.A. Processos Erosivos na Costa do Rio Grande do Sul, Brasil: Evidências de uma Provável Tendência Contemporânea de Elevação do Nível Relativo do Mar. In Proceedings of the 2nd Congresso da Associação Brasileira de Estudos do Quaternário, Rio de Janeiro, Brazil, 10–16 July 1989; UFRJ: Rio de Janeiro, Brazil, 1989; p. 16. [Google Scholar]
  6. Barboza, E.G.; Tomazelli, L.J. Erosional features of the eastern margin of the Patos Lagoon, southern Brazil: Significance for Holocene history. J. Coast. Res. 2003, 35, 260–264. [Google Scholar]
  7. Dillenburg, S.R.; Tomazelli, L.J.; Hesp, P.A.; Barboza, E.G.; Clerot, L.C.P.; Silva, D.B. Stratigraphy and evolution of a prograded, transgressive dunefield barrier in southern Brazil. J. Coast. Res. 2006, 39, 132–135. [Google Scholar]
  8. Biancini da Silva, A. Mapeamento e Caracterização dos Depósitos em Subsuperfície do Setor Meridional da Planície Costeira sul de Santa Catarina. Master’s Thesis, PPGGEO/UFRGS, Porto Alegre, Brazil, 2009; 66p. Available online: http://hdl.handle.net/10183/220100 (accessed on 20 May 2021).
  9. Dillenburg, S.R.; Barboza, E.G.; Rosa, M.L.C.C.; Caron, F.; Sawakuchi, A. The complex prograded Cassino barrier in Southern Brazil: Geological and morphological evolution and records of climatic, oceanographic and sea-level changes in the last 7-6 ka. Mar. Geol. 2017, 390, 106–119. [Google Scholar] [CrossRef]
  10. van Heteren, S.; Huntley, D.J.; van de Plassche, O.; Lubberts, R.K. Optical dating of dune sand for the study of sea-level change. Geology 2000, 28, 411–414. [Google Scholar] [CrossRef]
  11. Brooke, B.P.; Huanga, Z.; Nicholas, W.A.; Oliver, T.S.N.; Tamura, T.; Woodroffe, C.D.; Nichol, S.L. Relative sea-level records preserved in Holocene beach-ridge strandplains–An example from tropical northeastern Australia. Mar. Geol. 2019, 411, 107–118. [Google Scholar] [CrossRef]
  12. Oliver, T.S.N.; Murray-Wallace, C.V.; Woodroffe, C.D. Holocene shoreline progradation and coastal evolution at Guichen and Rivoli Bays, southern Australia. Holocene 2020, 30, 106–124. [Google Scholar] [CrossRef]
  13. Costas, S.; Ferreira, O.; Plomaritis, T.A.; Leorri, E. Coastal barrier stratigraphy for Holocene high-resolution sea-level reconstruction. Sci. Rep. 2016, 6, 38726. [Google Scholar] [CrossRef] [Green Version]
  14. Dougherty, A.J. Extracting a record of Holocene storm erosion and deposition preserved in the morphostratigraphy of a prograded coastal barrier. Cont. Shelf Res. 2014, 86, 116–131. [Google Scholar] [CrossRef]
  15. Dillenburg, S.R.; Barboza, E.G. The Strike-Fed Sandy Coast of Southern Brazil. In Sedimentary Coastal Zones from High to Low Latitudes: Similarities and Differences; Martini, I.P., Wanless, H.R., Eds.; Geological Society: London, UK, 2014; Special Publications 388; pp. 333–352. [Google Scholar] [CrossRef]
  16. Fernandez, G.B.; Rocha, T.B.; Barboza, E.G.; Dillenburg, S.R.; Rosa, M.L.C.C.; Angulo, R.J.; Souza, M.C.; Oliveira, L.H.S.; Dominguez, J.M.L. Natural Landscapes Along Brazilian Coastline. In The Physical Geography of Brazil-Environment, Vegetation and Landscape; Salgado, A.A.R., Santos, L.J.C., Paisani, J.C., Eds.; Springer: Cham, Switzerland, 2019; pp. 199–218. [Google Scholar] [CrossRef]
  17. Rosa, M.L.C.C.; Tomazelli, L.J.; Costa, A.F.U.; Barboza, E.G. Integração de métodos potenciais (gravimetria e magnetometria) na caracterização do embasamento da região sudoeste da Bacia de Pelotas, sul do Brasil. Rev. Bras. Geofís. 2009, 27, 641–657. [Google Scholar] [CrossRef] [Green Version]
  18. Contreras, J.; Zühlke, R.; Bowman, S.; Bechstädt, T. Seismic stratigraphy and subsidence analysis of the southern Brazilian margin (Campos, Santos and Pelotas basins). Mar. Pet. Geol. 2010, 27, 1952–1980. [Google Scholar] [CrossRef]
  19. Dillenburg, S.R.; Barboza, E.G. Long and short term progradation of a regressive barrier in southern Brazil. J. Coast. Res. 2009, 56, 599–601. [Google Scholar]
  20. Hesp, P.A. Conceptual models of the evolution of transgressive dunefield systems. Geomorphology 2013, 199, 138–149. [Google Scholar] [CrossRef]
  21. Villwock, J.A.; Tomazelli, L.J.; Loss, E.L.; Dehnhardt, E.A.; Horn Filho, N.O.; Bachi, F.A.; Denhardt, B.A. Geology of the Rio Grande do Sul Coastal Province. In Quaternary of South America and Antarctic Peninsula; Rabassa, J., Balkema, A.A., Eds.; CRC Press: Rotterdam, The Netherlands, 1986; Volume 4, pp. 79–97. [Google Scholar] [CrossRef]
  22. Rosa, M.L.C.C.; Barboza, E.G.; Dillenburg, S.R.; Tomazelli, L.J.; Ayup-Zouain, R.N. The Rio Grande do Sul (southern Brazil) shoreline behavior during the Quaternary: A cyclostratigraphic analysis. J. Coast. Res. 2011, 64, 686–690. [Google Scholar]
  23. Rosa, M.L.C.C.; Barboza, E.G.; Abreu, V.S.; Tomazelli, L.J.; Dillenburg, S.R. High-frequency sequences in the Quaternary of Pelotas Basin (coastal plain): A record of degradational stacking as a function of longer-term base-level fall. Braz. J. Geol. 2017, 47, 183–207. [Google Scholar] [CrossRef] [Green Version]
  24. Andrade, M.M.; Toldo, E.E.; Nunes, J.C.R. Tidal and subtidal oscillations in a shallow water system in southern Brazil. Brazil. J. Oceanog. 2018, 66, 245–254. [Google Scholar] [CrossRef]
  25. Cecilio, R.O.; Dillenburg, S.R. An ocean wind-wave climatology for the Southern Brazilian Shelf. Part I: Problem presentation and model validation. Dyn. Atmos. Ocean 2020, 89, 101101. [Google Scholar] [CrossRef]
  26. Barletta, R.C.; Calliari, L.J. Determinação da intensidade das tempestades que atuam no litoral do Rio Grande do Sul, Brasil. Pesqui. Geociênc. 2001, 28, 117–124. [Google Scholar] [CrossRef] [Green Version]
  27. Alvares, C.A.; Stape, J.L.; Sentelhas, P.C.; Gonçalves, J.L.M.; Sparovek, G. Köppen’s climate classification map for Brazil. Meteorol. Z. 2014, 22, 711–728. [Google Scholar] [CrossRef]
  28. Martinho, C.T.; Hesp, P.A.; Dillenburg, S.R. Morphological and temporal variations of transgressive dunefields of the northern and mid-littoral Rio Grande do Sul coast, Southern Brazil. Geomorphology 2010, 117, 14–32. [Google Scholar] [CrossRef]
  29. Corrêa, I.C.S. Les variations du niveau de la merdurant les derniers 17.500 ans BP. L’exemple de la plateforme continentale du Rio Grande do Sul-Bresil. Mar. Geol. 1996, 130, 163–178. [Google Scholar] [CrossRef]
  30. Angulo, R.J.; Giannini, P.C.F.; Suguio, K.; Pessenda, L.C.R. The relative sea-level changes in the last 5500 years southern Brazil (Laguna-Imbituba region, Santa Catarina State) based on vermetid 14C ages. Mar. Geol. 1999, 159, 327–339. [Google Scholar] [CrossRef]
  31. Martínez, S.; Rojas, A. Relative sea level during the Holocene in Uruguay. Palaeogeogr. Palaeoclimatol. Palaeoecol. 2013, 374, 123–131. [Google Scholar] [CrossRef]
  32. Verocai, J.E.; Hagy, G.J.; Bidegain, M. Seal-level trends along freshwater and seawater mixing in the Uruguayan Rio de la Plata estuary and Atlantic Ocean coast. Int. J. Mar. Sci. 2016, 6, 1–18. [Google Scholar] [CrossRef]
  33. Leandro, C.G.; Barboza, E.G.; Caron, F.; Jesus, F.A.N. GPR trace analysis for coastal depositional environments of southern Brazil. J. Appl. Geophys. 2019, 162, 1–12. [Google Scholar] [CrossRef]
  34. Daniels, J.; Roberts, R.; Vendl, M. Ground penetrating radar for the detection of liquid contaminants. J. Appl. Geophys. 1995, 33, 195–207. [Google Scholar] [CrossRef]
  35. Martinez, A.; Byrnes, A.P. Modeling Dielectric-Constant Values of Geologic Materials: An Aid to Ground-Penetrating Radar Data Collection and Inter-pretation. Current Research in Earth Sciences, Kansas Geological Survey 2001, Bulletin 247, Part 1. Available online: http://www.kgs.ku.edu/Current/2001/martinez/martinez1.html (accessed on 22 May 2021).
  36. Lima, L.G.; Parise, C.K. Holocene coastal evolution of the transition from transgressive to regressive barrier in southern Brazil. Catena 2020, 185, 104263. [Google Scholar] [CrossRef]
  37. Barboza, E.G.; Dillenburg, S.R.; Lopes, R.P.; Rosa, M.L.C.C.; Caron, F.; Abreu, V.S.; Manzolli, R.P.; Nunes, J.C.R.; Weschenfelder, J.; Tomazelli, L.J. Geomorphological and Stratigraphic Evolution of a Fluvial Incision in the Coastal Plain and Inner Continental Shelf in Southern Brazil. Mar. Geol. 2021, 437, 106514. [Google Scholar] [CrossRef]
  38. Payton, C.E. Seismic Stratigraphy—Applications to Hydrocarbon Exploration; EUA: Tulsa, OK, USA, 1977; Volume 26, 516p. [Google Scholar] [CrossRef]
  39. Neal, A. Ground-penetrating radar and its use in sedimentology: Principles, problems and progress. Earth-Sci. Rev. 2004, 66, 261–330. [Google Scholar] [CrossRef]
  40. Mitchum, R.M.; Vail, P.R.; Sangree, J.B. Seismic Stratigraphy and Global Changes of Sea Level, Part 6: Stratigraphy interpretation of seismic reflection patterns in depositional sequences. In Seismic Stratigraphy—Applications to Hydrocarbon Exploration; Payton, C.E., Ed.; EUA: Tulsa, OK, USA, 1977; Volume 26, pp. 117–133. [Google Scholar] [CrossRef]
  41. Abreu, V.S.; Neal, J.; Vail, P.R. Integration of Sequence Stratigraphy concepts. In Sequence Stratigraphy of Siliciclastic Systems–The ExxonMobil Methodology: Atlas of Exercises; Abreu, V.S., Neal, J., Bohacs, K.M., Kalbas, J.L., Eds.; EUA: Houston, TX, USA, 2010; pp. 209–224. [Google Scholar]
  42. Barboza, E.G.; Rosa, M.L.C.C.; Dillenburg, S.R.; Biancini da Silva, A.; Tomazelli, L.J. Stratigraphic analysis applied on the recognition of the interface between marine and fluvial depositional systems. J. Coast. Res. 2014, 70, 205–210. [Google Scholar] [CrossRef]
  43. Neal, J.E.; Abreu, V.; Bohacs, K.M.; Feldman, H.R.; Pederson, K.H. Accommodation succession (δA/δS) sequence stratigraphy: Observational method, utility and insights into sequence boundary formation. J. Geol. Soc. 2016, 173, 803–816. [Google Scholar] [CrossRef]
  44. Aitken, M.J. An Introduction to Optical Dating; Oxford University Press: New York, NY, USA, 1998; 280p. [Google Scholar]
  45. Murray, A.S.; Wintle, A.G. The single aliquot regenerative dose protocol: Potential for improvements in reliability. Radiat. Meas. 2003, 37, 377–381. [Google Scholar] [CrossRef]
  46. Prescott, J.R.; Hutton, J.T. Cosmic ray contributions to dose rates for luminescence and ESR dating: Large depths and long-term time variations. Radiat. Meas. 1994, 23, 497–500. [Google Scholar] [CrossRef]
  47. R Core Team. R: A Language and Environment for Statistical Computing; R Foundation for Statistical Computing: Vienna, Austria, 2021; Available online: https://www.R.-project.org/ (accessed on 20 April 2021).
  48. Hesp, P.A.; Dillenburg, S.R.; Barboza, E.G.; Tomazelli, L.J.; Ayup-Zouain, R.N.; Esteves, L.S.; Gruber, N.L.S.; Toldo, E.E., Jr.; Tabajara, L.L.C.A.; Clerot, L.C.P. Beach Ridges, Foredunes or Transgressive Dunefields? Definitions and an Examination of the Torres to Tramandaí Barrier System, Southern Brazil. An. Acad. Bras. Ciênc. 2005, 77, 493–508. [Google Scholar] [CrossRef]
  49. Hesp, P.A.; Dillenburg, S.R.; Barboza, E.G.; Clerot, L.C.; Tomazelli, L.J.; Ayup-Zouain, R.N. Morphology of the Itapeva to Tramandaí transgressive dunefield barrier system and mid-to late Holocene sea level change. Earth Surf. Process. Landf. 2007, 32, 407–414. [Google Scholar] [CrossRef]
  50. Barboza, E.G.; Rosa, M.L.C.C.; Dillenburg, S.R.; Tomazelli, L.J. Preservation Potential of Foredunes in the Stratigraphic Record. J. Coast. Res. 2013, 65, 1265–1270. [Google Scholar] [CrossRef]
  51. Hesp, P.A. A 34 year record of foredune evolution, Dark Point, NSW, Australia. J. Coast. Res. 2013, 65, 1295–1300. [Google Scholar] [CrossRef]
  52. Dillenburg, S.R.; Hesp, P.A.; Keane, R.; Miot da Silva, G.; Sawakuchi, A.O.; Moffat, I.; Barboza, E.G.; Bitencourt, V.J.B. Geochronology and evolution of a complex barrier, Younghusband Peninsula, South Australia. Geomorphology 2020, 354, 107044. [Google Scholar] [CrossRef]
  53. Dougherty, A.J.; Choi, J.H.; Turney, C.S.M.; Dosseto, A. Optimizing the utility of combined GPR, OSL, and Lidar (GOaL) to extract paleoenvironmental records and decipher shoreline evolution. Clim. Past. 2019, 15, 389–404. [Google Scholar] [CrossRef] [Green Version]
  54. Jol, H. Ground Penetrating Radar: Theory and Applications; Elsevier: Amsterdam, The Netherlands, 2009; 544p. [Google Scholar]
  55. Venkateswarlu, B.; Tewari, V.C. Geotechnical Applications of Ground Penetrating Radar (GPR). J. Ind. Geol. Cong. 2014, 6, 35–46. [Google Scholar]
  56. Laborel, J. Fixed marine organisms as biological indicator for the study of recent sea level and climatic variations along the Brazilian tropical coast. In Proceedings of the International Symposium on Coastal Evolution in the Quaternary, São Paulo, Brasil, 11–18 September 1978; Suguio, K., Fairchild, T.R., Martin, L., Flexor, J.M., Eds.; USP: São Paulo, Brasil, 1979; pp. 193–211. [Google Scholar]
  57. Laborel, J. Vermetid gastropods as sea-level indicators. In Sea-Level Research: A Manual for the Collection and Evaluation of Data; Van de Plassche, O., Ed.; Geo Books: Norwich, UK, 1986; pp. 281–310. [Google Scholar]
  58. Toniolo, T.F.; Giannini, P.C.F.; Angulo, R.J.; Souza, M.C.; Pessenda, L.C.R.; Spotorno-Oliveira, P. Sea-level fall and coastal water cooling during the Late Holocene in Southeastern Brazil based on vermetid bioconstructions. Mar. Geol. 2020, 428, 106281. [Google Scholar] [CrossRef]
  59. Bracco, R.; García-Rodríguez, F.; Inda, H.; del Puerto, L.; Castiñeira, C.; Panario, D. Niveles relativos del mar durante el Pleistoceno final–Holoceno y las costas de Uruguay. In El Holoceno en la zona costera del Uruguay; García-Rodríguez, F., Ed.; CSIC-Universidad de la República, Facultad de Ciencias: Montevideo, Uruguay, 2011; pp. 65–92. [Google Scholar]
  60. Prieto, A.R.; Mourelle, D.; Peltier, W.R.; Drummond, R.; Vilanova, I.; Ricci, L. Relative sea-level changes during the Holocene in the Río de la Plata, Argentina and Uruguay: A review. Quat. Int. 2017, 442, 35–49. [Google Scholar] [CrossRef]
  61. Kuhn, L.A.; Souza, P.A.; Cancelli, R.R.; Silva, W.G.; Macedo, R.B. Paleoenvironmental evolution of the coastal plain of Southern Brazil: Palynological data from a Holocene core in Santa Catarina State. An. Acad. Bras. Ciênc. 2017, 89, 2581–2595. [Google Scholar] [CrossRef] [Green Version]
  62. Masetto, E.; Lorscheitter, M.L. Vegetation dynamics during the last 7500 years on the extreme southern Brazilian coastal plain. Quat. Int. 2019, 524, 48–56. [Google Scholar] [CrossRef]
  63. Esteves, L.S.; Dillenburg, S.R.; Toldo, E.E., Jr. Alongshore Patterns of Shoreline Movements in Southern Brazil. J. Coast. Res. 2006, 39, 215–219. [Google Scholar]
  64. Annan, A.P.; Davis, J.L. Impulse radar sounding in permafrost. Radio Sci. 1976, 11, 383–394. [Google Scholar] [CrossRef]
  65. Olhoeft, G.R. Applications and limitations of ground penetrating radar. In Proceedings of the 54th Annual International Meeting, Atlanta, GA, USA, 2–6 December 1984; Society of Exploration Geophysicists: Tulsa, OK, USA; pp. 147–148. [Google Scholar]
Figure 2. (A) Sea-level envelope for the Brazilian coast south of 28° S (modified after [4]); (B) Sea-level envelope for the Uruguayan coast (modified after [31]).
Figure 2. (A) Sea-level envelope for the Brazilian coast south of 28° S (modified after [4]); (B) Sea-level envelope for the Uruguayan coast (modified after [31]).
Geosciences 11 00326 g002
Figure 3. The short interval of the GPR profile above was performed at the Arroio do Sal (see location in Figure 1). (A) Processed section showing, in the dashed red rectangle, the contact between the foredune and the backshore (zoom in for a better observation); (B) Interpreted section showing the four radarfacies. The arrows represent the reflection terminations (toplap: orange, downlap: green, truncations: black). See insertion of this GPR interval in the large and stacked GPR record of Figure 5B.
Figure 3. The short interval of the GPR profile above was performed at the Arroio do Sal (see location in Figure 1). (A) Processed section showing, in the dashed red rectangle, the contact between the foredune and the backshore (zoom in for a better observation); (B) Interpreted section showing the four radarfacies. The arrows represent the reflection terminations (toplap: orange, downlap: green, truncations: black). See insertion of this GPR interval in the large and stacked GPR record of Figure 5B.
Geosciences 11 00326 g003
Figure 4. (A) The morphological transition between the foredune and the backshore at the Curumim barrier and (B) detail of the sharp contact (black arrow) between backshore sands (light grey color) and foredune sands (white to pale yellow). At this place the contact is located 1.5 m above mean sea level (picture orientation is NW-SE).
Figure 4. (A) The morphological transition between the foredune and the backshore at the Curumim barrier and (B) detail of the sharp contact (black arrow) between backshore sands (light grey color) and foredune sands (white to pale yellow). At this place the contact is located 1.5 m above mean sea level (picture orientation is NW-SE).
Geosciences 11 00326 g004
Figure 5. (A) Location of the GPR profile along the Arroio do Sal (UTM–Datum: WGS84). (B) Correspondent GPR section. The 24 red dots indicate the positions of the contacts between the foot of the foredune and the backshore. The large intervals along the profile showing no dots correspond to places eroded by channelized waters (black arrows). The orange points represent a projection of the location of OSL ages. The black dashed rectangle corresponds to the location of the GPR profile of Figure 3. Four enlarged (detailed) GPR sections from where the 24 contacts between the foot of the foredune and the backshore were extracted are presented as Supplementary Material.
Figure 5. (A) Location of the GPR profile along the Arroio do Sal (UTM–Datum: WGS84). (B) Correspondent GPR section. The 24 red dots indicate the positions of the contacts between the foot of the foredune and the backshore. The large intervals along the profile showing no dots correspond to places eroded by channelized waters (black arrows). The orange points represent a projection of the location of OSL ages. The black dashed rectangle corresponds to the location of the GPR profile of Figure 3. Four enlarged (detailed) GPR sections from where the 24 contacts between the foot of the foredune and the backshore were extracted are presented as Supplementary Material.
Geosciences 11 00326 g005
Figure 6. The sea-level curve of the southernmost Brazilian coast, this work, based on 24 foredune-backshore contact points identified by GPR at Arroio do Sal (29°30′ S-49°53′ W–Figure 5). These paleo sea levels were adjusted to the current mean sea level (MSL). The 5th-order polynomial curve (R2 = 0.966, p < 0.001) and the envelope (black dashed lines), this work, represent the lower and upper limit of the predicted intervals (CI, 99%). OSL ages are displayed at the superior x-axis. We have also included a modern point at (0, 0) to indicate the present mean sea level. The red dotted and dashed lines represent the envelopes of Imbituba-Laguna [4] and Uruguay [31], respectively. These two envelopes were plotted according to axes Y (sea level altitude) and X (age).
Figure 6. The sea-level curve of the southernmost Brazilian coast, this work, based on 24 foredune-backshore contact points identified by GPR at Arroio do Sal (29°30′ S-49°53′ W–Figure 5). These paleo sea levels were adjusted to the current mean sea level (MSL). The 5th-order polynomial curve (R2 = 0.966, p < 0.001) and the envelope (black dashed lines), this work, represent the lower and upper limit of the predicted intervals (CI, 99%). OSL ages are displayed at the superior x-axis. We have also included a modern point at (0, 0) to indicate the present mean sea level. The red dotted and dashed lines represent the envelopes of Imbituba-Laguna [4] and Uruguay [31], respectively. These two envelopes were plotted according to axes Y (sea level altitude) and X (age).
Geosciences 11 00326 g006
Table 1. Radionuclide, water contents, measured dose rate, equivalent dose, and luminescence ages.
Table 1. Radionuclide, water contents, measured dose rate, equivalent dose, and luminescence ages.
SampleWater Content δ 1U (μg/g) from 234ThU (μg/g) 2 from 226Ra, 214Pb, 214BiU (μg/g) from 210PbTh (μg/g)2 from 208Tl, 212Pb, 228AcK (%)dD/dt (Gy/ka) #De (Gy)OSL-Age (ka)
WLL3801.0830.30 ± 0.090.45 ± 0.010.29 ± 0.091.57 ± 0.030.62 ± 0.010.918 ± 0.0296.14 ± 0.926.69 ± 1.02
WLL3811.1610.35 ± 0.080.47 ± 0.010.51 ± 0.091.38 ± 0.030.35 ± 0.010.641 ± 0.0313.84 ± 0.665.99 ± 1.07
WLL3831.1940.26 ± 0.080.33 ± 0.010.37 ± 0.090.95 ± 0.030.38 ± 0.010.593 ± 0.0342.98 ± 0.515.03 ± 0.91
WLL4041.1920.27 ± 0.050.31 ± 0.010.42 ± 0.050.97 ± 0.020.32 ± 0.010.543 ± 0.0302.38 ± 0.554.38 ± 1.04
1 Ratio wet sample to dry sample weight. Errors assumed 50% of (δ-1). 2 U and Th-content is calculated from the error weighted mean of the isotope equivalent contents. # The doserate dD/dt was calculated using the conversion factors of Adamiec and Aitken (1998).
Table 2. Table showing the distances of the 24 points identified, with their respective mean and corrected altitudes, associated with the correlated OSL ages.
Table 2. Table showing the distances of the 24 points identified, with their respective mean and corrected altitudes, associated with the correlated OSL ages.
Distance (m)Altitude (m asl) GNSS68% Accuracy HZ (m)68% Accuracy VT (m)Altitude of Paleo Sea Levels (m asl)Correlated OSL Age (ka BP)
165.302.00.10.30.5
265.802.30.10.30.8
364.702.50.10.31.0
420.802.40.10.30.9
523.202.30.10.30.8
657.303.00.10.31.5
743.503.00.10.31.5
1000.003.00.10.31.54.38 ± 1.04
1737.263.20.10.31.7
2320.003.30.10.31.85.03 ± 0.91
2411.563.30.10.31.8
2457.563.30.10.31.8
2557.163.40.10.31.9
2596.163.30.10.31.8
2733.933.30.10.31.8
2988.933.20.10.31.7
2997.833.10.10.31.6
3070.632.40.10.30.9
3301.432.20.10.30.7
3458.831.90.10.30.4
3580.001.00.10.3−0.55.99 ± 1.07
3705.83−0.20.10.3−1.7
3917.53−0.40.10.3−1.9
3967.73−0.70.10.3−2.2
3975.00−1.70.10.3−3.26.69 ± 1.02
3986.63−1.90.10.3−3.4
4019.63−2.20.10.3−3.7
4041.03−2.60.10.3−4.1
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Barboza, E.G.; Dillenburg, S.R.; do Nascimento Ritter, M.; Angulo, R.J.; da Silva, A.B.; da Camara Rosa, M.L.C.; Caron, F.; de Souza, M.C. Holocene Sea-Level Changes in Southern Brazil Based on High-Resolution Radar Stratigraphy. Geosciences 2021, 11, 326. https://0-doi-org.brum.beds.ac.uk/10.3390/geosciences11080326

AMA Style

Barboza EG, Dillenburg SR, do Nascimento Ritter M, Angulo RJ, da Silva AB, da Camara Rosa MLC, Caron F, de Souza MC. Holocene Sea-Level Changes in Southern Brazil Based on High-Resolution Radar Stratigraphy. Geosciences. 2021; 11(8):326. https://0-doi-org.brum.beds.ac.uk/10.3390/geosciences11080326

Chicago/Turabian Style

Barboza, Eduardo Guimarães, Sergio Rebello Dillenburg, Matias do Nascimento Ritter, Rodolfo José Angulo, Anderson Biancini da Silva, Maria Luiza Correaa da Camara Rosa, Felipe Caron, and Maria Cristina de Souza. 2021. "Holocene Sea-Level Changes in Southern Brazil Based on High-Resolution Radar Stratigraphy" Geosciences 11, no. 8: 326. https://0-doi-org.brum.beds.ac.uk/10.3390/geosciences11080326

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop