Geomorphic Response of the Georgia Bight Coastal Zone to Accelerating Sea Level Rise, Southeastern USA

: Synthesis of geologic and chronologic data generated from Holocene sedimentary sequences recovered along the inner continental shelf, shoreface, and modern coastal zone of the Georgia Bight reveal a synchronous sequence of paleoenvironmental events that occurred in response to rate of sea level rise tipping points. During the early Holocene (11.7–8.2 cal kyr BP), the paleoshoreline was overstepped and submerged by rapidly rising seas that averaged ~5 mm yr − 1 . Rates of rise during the middle Holocene (8.2–4.2 cal kyr BP) averaged ~2 mm yr − 1 and this deceleration resulted in the formation of coastal environments and sedimentary sequences that were subsequently reworked as the shoreface continued its landward and upward migration. The modern coastal zone emerged commensurate with the late Holocene (4.2–0 cal kyr BP), when the rate of sea level rise averaged <1 mm yr − 1 . Analysis of water level data collected at six NOAA tide gauge stations located along the Georgia Bight coast indicates the rate of relative sea level rise has increased from a historical average of 3.6 ± 0.2 mm yr − 1 (<1972 to 2022) to 6.6 ± 0.8 (1993 to 2022) and during the 21st century it has averaged 9.8 ± 0.3 mm yr − 1 (2003 to 2022). The current rate of sea level rise is nearly double the early Holocene rate of rise. Based upon a novel application of the principle of uniformitarianism (i.e., the past is the key to the future), the likely geomorphic trajectory of the Georgia Bight coastal zone under conditions of 21st century accelerating sea level rise will be one of increasing instability (e.g., coastal erosion) and flooding (e.g., overwash, breaching). Evidence of an emerging instability within the coastal zone has been previously reported throughout the region and supports the trajectory of geomorphic change proposed herein. This will ultimately result in the submergence of existing landscapes and replacement by estuarine and marine environments, which may hasten in pace and scale given the current rate of sea level rise is expected to continue accelerating throughout this century. These findings have not been previously reported and should be considered by coastal practitioners responsible for conceptualizing risk, as well as the formulation and implementation of adaptation action plans designed to mitigate threats to the built and natural environment induced by climate change.


Introduction
This investigation was designed to assist coastal practitioners, responsible for the formulation and implementation of climate change adaptation action plans, with a scientific basis for envisioning future conditions by answering the following question-what is the probable geomorphic response of the Georgia Bight coastal zone to accelerating sea level rise during the 21st century?The projections are based on a novel application of the principle of uniformitarianism; the past is the key to the future.The work began by conducting a regional analysis of geological investigations previously performed on the modern coastal zone and inner continental shelf to generate a unified model that associated widespread and synchronous changes in sedimentology, stratigraphy, and paleoenvironment with contemporaneous rates of sea level rise throughout the Holocene Epoch.Three distinct facies transitions are recognized, each corresponding to a substantial deceleration in the rate of sea level rise.These "tipping points" were compared with historical and recent rates of sea level rise observed at six NOAA tide gage stations located within the study domain to inform projections of geomorphic change likely to unfold along the Georgia Bight to the year 2050.While the findings of this study are likely applicable beyond mid-century, projections extending only a few decades into the future present less ambiguity to coastal zone managers, decision makers, and other practitioners [1].

Background
The origin of emergent marine and coastal strata observed globally along modern shorelines was the subject of considerable debate through the mid-20th century.Some argued a novel mechanism of long-term changes in mean sea level, while others preferred to invoke the historically accepted endogenous process of tectonism [2,3].However, this debate was promptly resolved upon the advent of the radiocarbon method as a reliable means of age determination (c.f., [4]).When used in tandem with bathymetric, sedimentologic, and stratigraphic data obtained on the continental shelves of the Gulf of Mexico [5][6][7][8] and the western Atlantic Ocean [9][10][11][12][13][14][15], investigators presented unequivocal evidence that a sea level rise of more than 100 m had occurred during the late Quaternary as a result of post-glacial melting.
Over the next several decades geologists, working along the eastern seaboard of the continental United States, often coupled their observations of sedimentology, stratigraphy, and geochronology with contemporaneous Holocene sea level rise to generate paleoenvironmental models of coastal evolution.For example, the presence of a thin and discontinuous sand sheet on the inner continental shelf was attributed to the reworking of older strata in response to a relatively rapid rise in sea level during early to middle Holocene [16][17][18][19][20]. Beneath modern coastal environments (e.g., barrier island, back barrier estuary, or bay), investigations often reported the presence of a transgressive facies sequence consisting of Pleistocene fluvial deposits overlain by estuarine muds and capped by barrier island sands (e.g., overwash, flood tidal shoals, aeolian) [21][22][23][24].This sequence was interpreted to have formed by the deposition and preservation of coastal sediments during the middle to late Holocene in response to a relatively slow rate of sea level rise.Located between the inner shelf sand sheet and modern coastal zone is the shoreface and early workers (cf, [23]) typically reported the presence of a transgressive facies sequence that was partially (nearshore) to completely reworked (offshore) by the surf zone and associated ravinement surface as it migrated both landward and upward in response to Holocene sea level rise.These observations and those of the many studies that followed formed the database and context from which the findings of this study are grounded.

Study Area
This study focused on an area of the southeast United States known as the Georgia Bight (Figure 1) and was selected to conform with the boundaries of previous regional assessments [24,25].It extends from Cape Lookout, North Carolina, southward to Fernandina Beach, Florida, and seaward across the inner continental shelf to a depth of ~25 m.This depth roughly coincides with the elevation of sea level at the onset of the Holocene Epoch (~11.7 cal kyr BP).The study area is subdivided into four coastal sectors: southern North Carolina, northern South Carolina, southern South Carolina, and the Georgia-Florida border.Sector boundaries roughly coincide with shoreline compartments II-V of Hayes [25] and are based upon prominent geomorphic features and physical oceanographic conditions (Table 1).

Sea Level Rise
Since the last glacial maximum, the globally averaged eustatic sea level has undergone a rise of at least 120 m, with the rate of ascent decelerating during the Holocene Epoch [26][27][28].Submergence records indicate that the initial transgression was marked by several meltwater pulses [29,30], characterized by extremely rapid rises exceeding 40 mm yr −1 [31][32][33][34], until ~11 cal kyr BP or the commencement of the Holocene Epoch.Subsequently, the long-term (spanning centuries) average trajectory of relative sea level rise within the study area is conventionally represented as a smooth curve illustrating a deceleration to the present [35][36][37] (Figure 2).The deceleration is attributed to a reduction of meltwater input and diminishing response of the Earth's mantle to glacial-isostatic adjustment [27][28][29].In addition to global eustatic, steric, and glacial-isostatic contributions to relative sea level rise, the submergence history of the Georgia Bight has also been influenced by two regional neotectonic features; the Cape Fear Arch and Ocala Uplift (Figure 1) [25,30].Cape Fear Arch uplift is estimated to be 0.24 mm yr −1 [31].Vertical motion associated the Ocala Uplift appears to be <0.01 mm yr −1 [32] and thus its contribution to relative sea level rise (hereafter sea level rise) during the Holocene Epoch is considered negligible for the purposes of this study.Supporting information is provided in Supplementary Materials Table S1.Superimposed are the three submergence curves considered in this study: North Carolina (NC) [33], northeast Florida (FL) [35], and the wider Caribbean (WC) [34].Holocene Epoch subdivisions after Walker et al. [36].SC = South Carolina.Supporting information is provided in Supplementary Materials Table S1.Superimposed are the three submergence curves considered in this study: North Carolina (NC) [33], northeast Florida (FL) [35], and the wider Caribbean (WC) [34].Holocene Epoch subdivisions after Walker et al. [36].SC = South Carolina.
Three regionally relevant reconstructions were used in this study to characterize sea level rise during the Holocene marine transgression [33][34][35] (Figure 2, Table 2).While differences between them are obvious and attributed to the phenomena described previously, the synchronous deceleration is apparent.When paired with the stratigraphical subdivisions Coasts 2024, 4 5 (aka stages) of the Holocene Epoch [36], three distinct intervals of decelerating submergence are evident: (1) an interval of rapid rise during the early Holocene (11.7 to 8.2 cal kyr BP) that averaged ~5.2 mm yr −1 , (2) an intermediate rate of rise averaging ~2.4 mm yr −1 that occurred during the middle Holocene (8.2 to 4.2 cal kyr BP), and (3) a relatively slow late Holocene rise that averaged ~0.8 mm yr −1 .These trends are consistent with the results of other sea level studies conducted along the Atlantic coast of the United States [29,37,38].

Response to Holocene Sea Level Rise
The methodological framework employed in this investigation is rooted in the seminal contributions of Fischer (1961) and Curray (1964) [5,9].They recognized that the stratigraphic record of marine transgressions and regressions is influenced by two fundamental factors: (a) fluctuations in sea level and (b) the net sedimentation rate.For instance, shoreline transgression occurs when the rate of sea level rise is faster than the rate of net sedimentation (i.e., sediment deficit).In the context of the Holocene marine transgression along the Georgia Bight coast, the rate of sea level rise was not constant but rather exhibited a decelerating trend over time.This deceleration coincided with changing patterns of sedimentation and stratigraphy.Consequently, if variations in the rate of Holocene sea level rise correspond with regionally significant shifts in paleoenvironmental evolution, these rates or tipping points can be used to model the anticipated response or geomorphic trajectory of the modern coastal zone to accelerating 21st century sea level rise.
The project began by establishing a regional conceptual framework of coastal evolution.This was accomplished by synthesizing the findings of previously published studies that (a) were executed along the contemporary coastline and/or inner continental shelf, (b) specifically focused on the Holocene stratigraphic sequence, (c) provided detailed characterizations of sedimentology, stratigraphy, and geochronology, and (d) incorporated paleoenvironmental descriptions pertaining to coastal evolution within the framework of concurrent sea level fluctuations.The evolution framework was established in two steps.First, a comparison of the aggregates studies conducted within each of the coastal sectors was performed to identify shared paleo-environmental events as evidenced in the stratigraphic record (e.g., deposition of a transgressive sand sheet).The geochronologic boundaries (e.g., start, end) of those shared events were then used to determine synchroneity within each sector.Once the initial, sector-based evaluation was completed, a between-sector comparison was undertaken to generate a regional, comprehensive understanding of how the Georgia Bight coastal zone evolved under conditions of Holocene sea level rise.
The geochronology of regionally synchronous paleoenvironmental events was constrained on the basis of (a) radiocarbon dates (aka absolute) generated from samples collected from a specific sediment type (e.g., saltwater peat) or environment of deposition (e.g., estuary) and (b) inferred ages derived by superimposing the depth(s) of in situ coastal sedimentary sequences onto the suite of sea level curves considered in this study (hereafter depth-age relationships).All radiocarbon ages are expressed in cal kyr BP.Conventional dates were calibrated using the online IntCal20 model by OxCal [39][40][41][42][43].The reported dates are relative to 1950 BP.All the geochronologic data used in this study are provided in Supplementary Materials Table S1.These were then merged with three Holocene sea level histories to identify tipping points that corresponded to regionally synchronous shifts in the geomorphic trajectory.

Sea Level Rise in the 20th and 21st Centuries
The investigation into historical and contemporary sea level behavior along the coast of the Georgia Bight involved the utilization of six tide gauges (Figure 1), as detailed in the supplementary data accompanying the 2022 technical report by the National Oceanic and Atmospheric Administration (NOAA) [44].NOAA's mean sea level monthly data products (i.e., monthly mean, minimum, and maximum values) [45] were employed for this analysis.The data for each station were initially transformed from the mean sea level datum to mean sea level at the year 2000 by computing the mean sea level value for the year 2000 and subtracting this value from all observations.This normalization was undertaken to maintain consistency with other studies [44].Trends in sea level rise were enumerated through a multi-parameter linear least-square estimator to determine optimal-fit parameters (i.e., intercept, slope) and their associated uncertainties using the Python programming language's scipy.optimize.curve_fitfunction.Additionally, best-fit calculations were performed using a linear model with a periodic (sinusoidal) component, necessitating the inclusion of two supplementary parameters (amplitude and phase) associated with the 18.61-year lunar tide periodicity, which can exert a substantial influence on mean and maximum sea level values [46].However, the outcomes derived from both linear and periodic analyses indicated only marginal deviations and, consequently, only the linear results were considered in this study.The best-fit analysis was conducted using three distinct data subsets, each characterized by a different start date and a common end date in 2022.The first subset encompassed the entire dataset (i.e., t start -2022).The second and third subsets focused on station data collected during the periods 1993-2022 and 2003-2022, respectively.The three observational-based means were then extrapolated to the year 2050 and compared to the scenario-based trajectories calculated by the National Oceanic and Atmospheric Administration (NOAA) [44].

21st Century Geomorphic Trajectories
Conceptualizations of the likely geomorphic response or trajectory of the Georgia Bight coastal zone to 21st century accelerating sea level rise were based upon the premise that the past is the key to the future.Specifically, the trajectory can be elucidated by comparison between (a) sea level rise tipping points documented during the Holocene marine transgression and (b) rates of sea level rise derived from the six NOAA tide gages.The projections of this study are considered conservative since the rate of sea level rise is expected to continue accelerating throughout the century [44,[47][48][49].The authors have utilized this methodology in other studies undertaken to identify the Holocene rate of sea level rise tipping points [50,51] and construct predictive conceptual models of environmental change [52-54].

Sedimentology, Stratigraphy, and Geomorphology
Observations and generalizations regarding the evolution of the Georgia Bight coastal zone in response to the Holocene marine transgression are based upon a review of over 150 publications describing geological investigations conducting within the study domain.The sedimentologic, stratigraphic, and geomorphic data contained in 32 of those (Tables 3 and 4) were used to reconstruct the paleoenvironmental evolution of the region, as described below, and generate four regional geologic cross-sections considered representative of each coastal sector.Each cross-section is an illustrative summary of the sediments and sequences that were preserved along the pathway of the transgressing sea.The location of each cross-section was selected after reviewing the density and distribution of all data within each sector to ensure it was constructed using the largest data set possible.The spatial distribution of data used to construct each cross-section is illustrated in Figure 3 and the generalized geologic cross-sections are illustrated in Figure 4.
Table 3. Published information (x) utilized to construct geologic cross-sections considered representative of the four coastal sectors (see Figures 1 and 3).Definition of terms and abbreviations provided in Table 4. Matching cell colors highlight link between area of study and observations of original authors.Of the 32 publications, 21 were conducted on the inner continental shelf and all reported the presence of a thin (~1 m) and discontinuous transgressive sand sheet that unconformably overlies unconsolidated Pleistocene sands and older strata (Table 3).Local bathymetric relief is generally minimal (e.g., meters), although in some areas (e.g., location f in Figure 3) bedrock scarps exceed 5 m.None of the investigations conducted on the continental shelf report the presence of in situ coastal sediments (e.g., basal peat, estuarine mud).

Map Inset Key
Key to location of data shown in Figure 3.In some cases, the original studies were conducted in more than one coastal sector as defined in this investigation.Thus, some of the original studies are listed more than once.In those cases, the total aerial extent of that study is indicated by a numeric value (km 2 ).

Authors
Citation identifying author(s) whose published data and observations were used in this investigation to construct the generalized geologic cross-sections shown in Figure 4. Seven of the thirty-two studies were conducted over the shoreface (Table 3) and only one included a detailed description of subsurface geology [18].Based upon geophysical surveys conducted in the southern North Carolina sector, these authors report the presence of an inner shelf sand sheet that grades landward into a transgressive facies sequence consisting of basal Pleistocene sediments, overlain by Holocene lagoon and tidal inlet strata, in turn truncated by the modern ravinement surface of the Bogue Banks shoreface (Figure 4).
Eighteen of the thirty-two investigations were conducted in the modern coastal zone (e.g., barrier island, back barrier bay, mainland).This region exhibits the highest degree of geomorphic and stratigraphic variability within the study domain because (1) it is subject to the effects of physical oceanographic processes (e.g., waves and tides; Table 1) and (2) may host Pleistocene high stand deposits along the mainland shoreline or embedded in the Holocene back barrier bay sediment succession (Figure 4) [24,67,82,83].The alongshore variation in the presence and elevation of Pleistocene high stand deposits is attributed to regional differences in neotectonics [32].In all coastal sectors of the Georgia Bight except Georgia-Florida, investigators reported the presence of a transgressive facies sequence beneath the barrier island and back barrier bay environments.No investigations appear to Coasts 2024, 4 9 have been conducted on the subsurface geology of barrier islands located in the Georgia-Florida sector (Table 3), so the presence of a similar stratigraphy has yet to be determined.3.  3.

Geochronology
The literature review yielded 12 papers and 115 radiocarbon dates considered relevant to the geochronologic element of this study (Table 5, Supplementary Materials Table S1).The radiocarbon dates obtained from the inner shelf transgressive sand sheet vary widely (Figure 2) and often yield geochronological inversions because the stratigraphic sequence was generated by the reworking of Holocene, Pleistocene, and older materials along the leading edge of the transgression (c.f., [77]).Depth-age relationships suggest the sand sheet was deposited between ~10 (25 m) and ~8 (10 m) cal kyr BP; Figures 2 and 4) or early Holocene, as has been suggested by others [56].Long et al. [71] recovered a core in water depths of ~13 m at a site located just seaward of the modern shoreface in the northern South Carolina sector.The core penetrated ~1.5 m of sand that unconformably overlay Pleistocene muds.Eight radiocarbon dates were generated from the sand layer and no ages were inverted.The oldest sample, collected at the base of the sand layer, dated 8.5 cal kyr BP or terminal early Holocene.3. Circled question mark (?) indicates zone of uncertainty due to limited or unavailable data.Vertical exaggeration is ~700×.Transect locations shown in Figure 1.
Table 3. Published information (x) utilized to construct geologic cross-sections considered representative of the four coastal sectors (see Figures 1 and 3).Definition of terms and abbreviations provided in Table 4. Matching cell colors highlight link between area of study and observations of original authors.3. Circled question mark (?) indicates zone of uncertainty due to limited or unavailable data.Vertical exaggeration is ~700×.Transect locations shown in Figure 1.

Cross-Section Location
As noted previously, only one publication described the stratigraphy of the modern shoreface.However, the authors [18] did not provide any geochronological data.Depth-age relationships (Figures 2 and 4) suggest the lagoon strata, located at the base of transgressive facies sequence they encountered, was deposited ~8 cal kyr BP and therefore during the onset of the middle Holocene.
All samples collected from beneath the modern back barrier coastline and adjacent bay areas dated <7 cal kyr BP (Figure 2).These data plot along a concave upward curve that emulates contemporaneous sea level, and their ages suggest these depositional systems formed during the intermediate and terminal stages of middle Holocene and have persisted throughout the late Holocene marine transgression.A scatter of the data within each sector is apparent and attributed to (a) errors associated with estimating sample age and elevation (both in situ and relative to contemporaneous sea level), (b) the thickness of the sample interval (e.g., bulk peat samples vs. an individual shell), and (c) post-depositional autocompaction and/or sediment consolidation associated with the coring process (cf, [39,69]).Between sector differences are also apparent, with the data collected from the northern two sectors tacking more closely to the North Carolina sea level reconstruction of Kopp et al. [35] and the southern two more tightly with the Florida curve of Hawkes et al. [37].These trends are attributed to glacial-isostatic processes.When the geology and geochronology of the study area are considered in tandem with contemporaneous rates of Holocene sea level rise, three synchronous depositional events are revealed.The inner shelf transgressive sand sheet was generated as the coastline was overstepped in response to a relatively rapid rise in sea level (~5.2 mm yr −1 ) during the early Holocene.During the middle Holocene, sea level rise slowed to an average rate of ~2.4 mm yr −1 and although the data are limited, this corresponds to a transition from coastal overstep to the formation of an emergent, but landward migrating, erosional coastline.Sediment cores recovered from the modern back barrier coastline and bays yielded dates <7 cal kyr BP and record a period of continuous sedimentation and geomorphic stability that has persisted throughout the late Holocene sea, during which the rate of level rise averaged ~0.8 mm yr −1 .

Sea Level Rise in the 20th and 21st Centuries
The sea level data acquired from the six NOAA tide gauge stations (Figure 5, Table 6) reveal an acceleration in the average rate of sea level rise.The mean rates of historical sea level rise, computed utilizing the entire dataset, ranged from 2.22 ± 0.01 mm yr −1 (Fernandina Beach, Florida) to 4.27 ± 0.05 mm yr −1 (Beaufort, North Carolina), with an overall average of 3.26 ± 0.21 mm yr −1 .The mean rate of sea level rise since 1993 ranged from 4.80 ± 0.57 mm yr −1 (Fernandina Beach, Florida) to 6.32 ± 0.43 mm yr −1 (Charleston, South Carolina) and averaged 5.63 ± 0.81 mm yr −1 .The mean rate of sea level rise during the 21st century (2003-2022) ranged from 8.24 ± 1.83 mm yr −1 (Fernandina Beach, Florida) to 10.96 ± 1.40 mm yr −1 (Charleston, South Carolina) and averaged 9.82 ± 0.31 mm yr −1 .Differences in the rates of sea level rise among stations is attributed to variations in glacialisostatic adjustment [39].neotectonics (i.e., Cape Fear Arch), surface temperature [84], prevailing winds, semi-annual (e.g., sunny-day) and long-term (e.g., nodal) tidal cycles [46], the regional coastal sea level fingerprint [85], climate variability [86], and oceanic circulation (e.g., the Gulf Stream) [87] unique to each location.However, in the context of this study, the numerous processes that contribute to local variations in the rate of relative sea level rise are not relevant.Rather, it is their aggregate that determines the rate of rise to which the geomorphology has, is, and will respond.Observation-based extrapolations of sea level elevation in 2050 were calculated for each of the three time periods considered in this study (i.e., tstart-2022, 1993-2022, and 2003-2022).These are compared to the scenario-based trajectories of Sweet et al. [44] in Table 7.

21st Century Geomorphic Trajectories
The paleo-environmental observations culled from geological investigations conducted within the study domain reveal a shared regional and synchronous response to changing rates of Holocene sea level rise.An initial interval of coastal overstep, demonstrated by the presence of a transgressive sand sheet, occurred during the Early Holocene (11.7-8.2 cal kyr BP) when the rate of sea level rise exceeded 5 mm yr −1 .Towards the modern shoreface, the sand sheet thickens and is interpreted as a transgressive facies sequence.This interval of shoreline retreat occurred during the Middle Holocene (8.2-4.2 cal kyr BP), when the rate of SLR averaged ~2 mm yr −1 .The modern coastal zone emerged and stabilized during the Late Holocene (4.2 cal kyr BP-present), when sea level was rising at a rate of <1 mm yr −1 .Local differences in sedimentology, stratigraphy, and geochronology within and between each coastal sector are evident, but these differences did not mask the Observation-based extrapolations of sea level elevation in 2050 were calculated for each of the three time periods considered in this study (i.e., t start -2022, 1993-2022, and 2003-2022).These are compared to the scenario-based trajectories of Sweet et al. [44]

21st Century Geomorphic Trajectories
The paleo-environmental observations culled from geological investigations conducted within the study domain reveal a shared regional and synchronous response to changing rates of Holocene sea level rise.An initial interval of coastal overstep, demonstrated by the presence of a transgressive sand sheet, occurred during the Early Holocene (11.7-8.2 cal kyr BP) when the rate of sea level rise exceeded 5 mm yr −1 .Towards the mod-ern shoreface, the sand sheet thickens and is interpreted as a transgressive facies sequence.This interval of shoreline retreat occurred during the Middle Holocene (8.2-4.2 cal kyr BP), when the rate of SLR averaged ~2 mm yr −1 .The modern coastal zone emerged and stabilized during the Late Holocene (4.2 cal kyr BP-present), when sea level was rising at a rate of <1 mm yr −1 .Local differences in sedimentology, stratigraphy, and geochronology within and between each coastal sector are evident, but these differences did not mask the overall uniform and synchronous trajectory of coastal evolution during the Holocene marine transgression.
The observed 21st century rate of sea level rise along the Georgia Bight coastline is nearly double the fastest rate of rise recorded during the Holocene marine transgression and tracking with the intermediate to high sea level rise scenario-based trajectories of [44].When this region was last subjected to a rate of this magnitude, the geomorphology was overstepped as demonstrated by the presence of the thin and discontinuous inner shelf transgressive sand sheet.If the past is indeed the key to the future, it is highly probable the geomorphic trajectory of the modern coastal zone to a rate of rise of this magnitude will be a shift towards increasing instability (c.f., [74][75][76][77]).This will ultimately result in erosional shoreface retreat, coastal backstepping, and landscape submergence.Early indicators of destabilization will include increasing rate and extent of (a) saltwater encroachment, landward migration of salt tolerant wetland plant communities, (b) barrier overwash and/or breaching, and (c) shoreline erosion [52, [88][89][90][91][92][93][94].

Evidence of Recent Changes and Likely Outcomes of Faster Rates of Sea Level Rise
There have been numerous studies conducted in the study area that either (1) describe historical and recent environmental changes attributed to sea level rise or (2) model future conditions using various sea level rise scenarios.For example, on the back barrier coastlines of North Carolina, Voss et al. [80] reported steep declines in coastal marsh macrophyte production and related vertical sediment accumulation which they attributed to rising sea levels.The authors concluded a growing sediment deficit will ultimately lead to marsh drowning.In South Carolina, rising seas and salinity since 2014 were shown to have negatively impacted coastal forest productivity and composition [95], leading the investigators to conclude a transition to marsh habitat was inevitable.Working in the tidal freshwater coastal floodplains of South Carolina and Georgia, Jones et al. [96] and Wang et al. [97] observed a decline in coastal forest productivity (e.g., die off, tree regeneration, growth) and marsh encroachment, which the authors linked to rising sea level and increasing salinity.Some of their observations were reported from riverine study sites located as much as 30 km upstream of the Atlantic coastline.These areas represent the leading edge of the ongoing, anthropogenically induced marine transgression.On the Georgia coastal plain, Craft [98] determined the maximum rates of vertical elevation change or soil accretion in tidal freshwater forests averaged 2.2 mm yr −1 .He noted this was lower than the recent rate of sea level rise, which he reported at the time as 3.3 mm yr −1 , and concluded submergence was likely.Submergence is even more likely now, given the 21st century rates of rise recorded at NOAA stations proximal to his study area are nearly three-times faster (i.e., Fort Pulaski, 10.62 mm yr −1 ; Fernandina Beach, 8.24 mm yr −1 ).In north Florida, Morris et al. [99] measured vertical elevation change in both marsh and mangrove habitat.The authors noted that mangroves outperformed marsh, but even so concluded rates of sea level rise >10 mm yr −1 will result in forest submergence.The rate of sea level rise in this sector of the Georgia Bight coast is currently averaging ~8 mm yr −1 (Table 6).Finally, in a study conducted in National Wildlife Refuges located throughout the study domain [100], net vertical elevation measurements were collected from a variety of coastal wetland plant communities (e.g., oligohaline and salt marsh, forested) and none exceeded 3.1 mm yr −1 .This is much lower than the 21st century average rate of sea level rise (9.82 mm yr −1 ) reported in this study using all the NOAA station data (Table 6).In summary, these ecological studies provide regional evidence that the effects of accelerating sea level rise are already occurring in fresh-and tide-water wetlands located in the study area.The pace of wetland plant community shifts and localized submergence are likely to hasten as the rate of sea level rise continues to accelerate.
No reports or publications linking accelerating sea level rise to historic or recent changes in the geomorphology of sandy beaches and barrier islands were identified during the data mining phase of this investigation.However, this was not unexpected because the natural response of these features to rising seas (e.g., shoreline retreat, more frequent overwash) is often obscured or obstructed by human activities (cf, [101,102]).These include the installation of shore protection features (e.g., sea walls, bulkheads, revetments), beach nourishment, and inlet management (e.g., jetties, sediment dredging, sand back-or bypassing).However, these anthropogenic alterations ultimately exacerbate the vulnerability of the coastline to accelerating sea level rise by eliminating relevant ecosystem services (e.g., the storm buffering capacity of wetlands), lowering elevation (e.g., removal or scalping of coastal dunes), altering the sediment budget (e.g., managed inlets), and enhancing substrate erosion (e.g., sea wall wave reflection; [103]).
The geomorphic trajectory of the Georgia Bight will also be influenced by landfall of tropical storms and cyclones.These high energy events generally result in extensive shoreline, dune and bluff erosion, island overwash and breaching, and the destruction of shore protection structures (e.g., bulkheads, revetments), buildings, and infrastructure [104][105][106][107]. Storm impacts are especially pronounced across low-lying areas, where overwash or breaching are more likely occur [108].This overtopping of the shoreline will increase proportional to sea level rise.The intensity and destructiveness of tropical cyclones are expected to increase as the climate warms [109][110][111][112] and this will further accelerate geomorphic and ecologic instability.The combination of faster rates of sea level rise, increasing storminess, and ongoing anthropogenic alterations will surely increase the pace and scale of geomorphic, ecologic, and urban change towards 2050 [113,114]).This will create substantial new challenges to coastal zone practitioners, as the management and planning instruments of governance often lack the capacity to contextualize risk at the pace at which it is changing and the financial resources necessary to implement the plans that have already been accepted [115].Table 6.Average rates of sea level rise (mm yr −1 ) recorded at NOAA six tide gage stations [45].Error expressed as uncertainty.Station locations illustrated in Figure 1.

Location
NOAA ID Latitude Longitude Observation of Period Trend t start -2022 Table 7. Observation-based linear extrapolations (this study, left column) and scenario-based median estimates of sea level elevation with likely ranges bracketed (i.e., 17th-83rd percentile [44]) in 2050 relative to a baseline of 2000.The observation-based extrapolations are derived from best-fit analysis of all data (black), three decades: 1993-2022 (red), and 21st century: 2003-2022 (purple).The two scenario-based estimates (columns) that bound the observation-based extrapolations (rows) are shown for each of the six tide gauge stations as vertical dividing lines using the observation-based color scheme (i.e., red = 1993-2022, purple = 2003-2022).Values expressed in mm.Scenario-based estimates obtained from Supplementary Data Files to NOAA's 2022 sea level rise technical report [44,116]).

Conclusions
Synthesis of regional geologic and chronologic data associated with the Holocene marine transgression of the Georgia Bight coastal zone reveal a uniform, time synchronous geomorphic response to rate of sea level rise tipping points.During the early Holocene, coastal plain sediments were partially reworked and submerged by a relatively rapid sea level rise averaging at least 5 mm yr −1 .This paleoenvironmental sequence of events is evidenced by the presence of a thin and discontinuous transgressive sand sheet that unconformably overlies older strata on the inner continental shelf.The rate of sea level rise slowed to ~2 mm yr −1 during the middle Holocene and this was accompanied by 6. Conclusions Synthesis of regional geologic and chronologic data associated with the Holocene marine transgression of the Georgia Bight coastal zone reveal a uniform, time synchronous geomorphic response to rate of sea level rise tipping points.During the early Holocene, coastal plain sediments were partially reworked and submerged by a relatively rapid sea level rise averaging at least 5 mm yr −1 .This paleoenvironmental sequence of events is evidenced by the presence of a thin and discontinuous transgressive sand sheet that unconformably overlies older strata on the inner continental shelf.The rate of sea level rise slowed to ~2 mm yr −1 during the middle Holocene and this was accompanied by erosional shoreface retreat of the coastal zone and formation of a transgressive facies sequence.The modern coastline emerged during the late Holocene under conditions of a relatively slow rate of sea level rise that averaged <1 mm yr −1 .Data obtained from six NOAA tide gage stations located within the study domain indicate the rate of sea level rise has accelerated from an historical average of ~3 mm yr −1 (<1972-2022) to ~6 mm yr −1 over the past 30 years (1993-2022) and ~10 mm yr −1 during the 21st century (2003)(2004)(2005)(2006)(2007)(2008)(2009)(2010)(2011)(2012)(2013)(2014)(2015)(2016)(2017)(2018)(2019)(2020)(2021)(2022).The 21st century rate of rise is tracking with the intermediate to high scenario-based trajectories of Sweet et al. [51].At the current rate, water level elevation gain in 2050 will range between 35 to 46 cm relative to the year 2000 datum.
The 21st century rates vary between stations; however, all exceed the rate of sea level rise that accompanied the early Holocene marine transgression, a time when the coastal zone was rapidly transgressed and submerged.Recent studies have demonstrated the destabilizing effects of this acceleration on coastal geomorphology and ecology are already evident.Destabilization will likely expand in pace and scale in response to a persistent acceleration in the rate of 21st century sea level rise and especially during the latter half [44,49].By 2100, much of the modern, emergent coastal zone will have been replaced by open water marine and estuarine environments.

Figure 1 .
Figure 1.Google Earth image of Georgia Bight showing location of four coastal sectors considered in this study, NOAA tide gauge stations from which historical water level data was collected, and geologic cross-sections.Also shown are locations of two neotectonic features.

Figure 1 .
Figure 1.Google Earth image of Georgia Bight showing location of four coastal sectors considered in this study, NOAA tide gauge stations from which historical water level data was collected, and geologic cross-sections.Also shown are locations of two neotectonic features.

Figure 2 .
Figure 2. Geochronologic data plotted as a function sector location and type of material analyzed.Supporting information is provided in Supplementary Materials TableS1.Superimposed are the three submergence curves considered in this study: North Carolina (NC)[33], northeast Florida (FL)[35], and the wider Caribbean (WC)[34].Holocene Epoch subdivisions after Walker et al.[36].SC = South Carolina.

Figure 2 .
Figure 2. Geochronologic data plotted as a function sector location and type of material analyzed.Supporting information is provided in Supplementary Materials TableS1.Superimposed are the three submergence curves considered in this study: North Carolina (NC)[33], northeast Florida (FL)[35], and the wider Caribbean (WC)[34].Holocene Epoch subdivisions after Walker et al.[36].SC = South Carolina.
ObservationSedimentologic, stratigraphic, and/or paleoenvironmental observations reported by the authors and relevant to this investigation:TFSTransgressive facies sequence.Stratigraphic onlapping sequence of coastal facies indicative of the landward migration and upward translation of the coastal zone (i.e., erosional shoreface retreat).TSS Transgressive sand sheet.Thin (~1 m) and discontinuous veneer of unconsolidated sediment observed on the lower shoreface and inner shelf.These sediments unconformably overly older strata and are separated from them by a ravinement surface.HB Hardbottom.Indurated surface on the sea floor.May reflect the presence of older, lithified strata or syn-sedimentary induration at the sediment-water interface.Common on 'sediment-starved' continental shelves.PC Paleochannel.Incised 'V-shaped' valleys carved or cut into the underlying lithified strata or unconsolidated sediments, respectfully, by fluvial or tidal processes.Vertical and horizontal dimensions vary as a function of genesis and ravinement depth of cut.

Coasts 2024, 4 , 8 Figure 3 .
Figure 3. Google Earth images of the four coastal sectors illustrating the source (circles labeled in lower case letters) and location (polygons) of data harvested from previous investigations to generate geologic cross-sections shown in Figure 4. Location of NOAA tide stations (squares labeled in upper case letters) and geologic cross-section (red line) also shown.Key to data sources summarized in Table3.

Figure 3 .
Figure 3. Google Earth images of the four coastal sectors illustrating the source (circles labeled in lower case letters) and location (polygons) of data harvested from previous investigations to generate geologic cross-sections shown in Figure 4. Location of NOAA tide stations (squares labeled in upper case letters) and geologic cross-section (red line) also shown.Key to data sources summarized in Table3.

Figure 4 .
Figure 4. Simplified geological cross-sections of the four coastal sectors illustrating prominent stratigraphic units and associated paleoenvironment of deposition.Generalized elevation and bathymetry derived from NOAA nautical charts and USGS topographic maps (topos) as follows: Bogue Banks, NC (chart #11543 and 11541; Mansfield topo), Long Bay, SC (chart #11535 and 11531; North Island topo), Kiawah Island, SC (chart #11480 and 11521; Legareville and Kiawah topos), Fernandina Beach, FL (chart 11488; Amelia topo).Specific geological data used to construct cross-sections is summarized in Table3.Circled question mark (?) indicates zone of uncertainty due to limited or unavailable data.Vertical exaggeration is ~700×.Transect locations shown in Figure1.

Figure 4 .
Figure 4. Simplified geological cross-sections of the four coastal sectors illustrating prominent stratigraphic units and associated paleoenvironment of deposition.Generalized elevation and bathymetry derived from NOAA nautical charts and USGS topographic maps (topos) as follows: Bogue Banks, NC (chart #11543 and 11541; Mansfield topo), Long Bay, SC (chart #11535 and 11531; North Island topo), Kiawah Island, SC (chart #11480 and 11521; Legareville and Kiawah topos), Fernandina Beach, FL (chart 11488; Amelia topo).Specific geological data used to construct crosssections is summarized in Table3.Circled question mark (?) indicates zone of uncertainty due to limited or unavailable data.Vertical exaggeration is ~700×.Transect locations shown in Figure1.

Coasts 2024, 4 , 14 Figure 5 .
Figure 5. Observed sea level trends and average rate of rise (error = uncertainty) calculated for the six NOAA tide gage stations located within the study domain and illustrated as a function of duration of record.Observational-based (this study) extrapolations of sea level elevation to 2050 are shown as solid lines: black = historical average rate of rise (duration of record), red = recent average rate of rise (1993-2022), purple = 21st century average rate of rise (2003-2022).Scenario-based [44] extrapolations of sea level elevation to 2050 are shown as dashed lines: blue = intermediate-high, green = high.The shaded areas delineate the 17th-83rd percentile range of uncertainty for each of the two scenario-based projections.
Those authors considered five emissions-based sea level trajectories: low, intermediate low, intermediate, intermediate high, and high.Our sea level observation-based extrapolations indicate recent rates (1993-2022) are tracking with the low to intermediatelow scenario-based estimates.The 21st century trends (2003-2022) are tracking with the intermediate to high scenario-based estimates.At the current 21st century rate of rise, water level elevations above the 2000 datum will reach +35 cm (Fernandina Beach, Florida) to +46 cm (Wilmington, North Carolina) by 2050.

Figure 5 .
Figure 5. Observed sea level trends and average rate of rise (error = uncertainty) calculated for the six NOAA tide gage stations located within the study domain and illustrated as a function of duration of record.Observational-based (this study) extrapolations of sea level elevation to 2050 are shown as solid lines: black = historical average rate of rise (duration of record), red = recent average rate of rise (1993-2022), purple = 21st century average rate of rise (2003-2022).Scenario-based [44] extrapolations of sea level elevation to 2050 are shown as dashed lines: blue = intermediate-high, green = high.The shaded areas delineate the 17th-83rd percentile range of uncertainty for each of the two scenario-based projections.

Table 1 .
[26] coastal sectors considered in this study.Physical oceanographic data obtained from Hubbard et al.[26].

Table 1 .
[26] coastal sectors considered in this study.Physical oceanographic data obtained from Hubbard et al.[26].

Table 2 .
Average rate of sea level rise during early, middle, and late Holocene.Values derived from data shown in Figure2.AuthorRate of Sea Level Rise (mm yr−1 ) Duration (Years) Late Holocene (0-4.

Table 2 .
Average rate of sea level rise during early, middle, and late Holocene.Values derived from data shown in Figure2.

Table 4 .
Definition of Table 3 terms and abbreviations.Subtidal or intertidal environments located between the shorelines of the mainland and barrier island.Typical environments include bay, overwash, marsh, and relict or active flood tidal shoals.Bi Barrier island.Emergent sand body located between the back barrier bay and Atlantic Ocean.Sf Shoreface.Subtidal, seaward dipping sand wedge extending from the barrier island shoreline to a depth of ~10 m.Sh Shelf.Relatively flat and featureless offshore zone extending from the toe of shoreface to the shelf break (~50 m).Along the shoreface boundary and extending seaward to ~20 m, this zone may host ebb tidal deltas and linear sand shoals with significant bathymetric relief.Methods Methodological approach used by original author(s) in their investigations: H Historical. Includes aerial photographs, nautical charts, field observations, and/or grab samples.CB Core borings.Includes core borings recovered using brute force, rotary, vibratory, and gravity methods.GP Geophysical.Includes survey records obtained using seismic reflection, sidescan sonar, swath bathymetry, and ground penetrating radar.

Table 5 .
Summary of geochronologic sample data used to construct Figure2.A detailed description of these data is contained in Supplementary Materials TableS1.Sector abbreviations: sNC = southern North Carolina, nSC = northern South Carolina, sSC = southern South Carolina, G-F = Georgia-Florida.
in Table 7.Those authors considered five emissions-based sea level trajectories: low, intermediate low, intermediate, intermediate high, and high.Our sea level observation-based extrapolations indicate recent rates (1993-2022) are tracking with the low to intermediate-low scenariobased estimates.The 21st century trends (2003-2022) are tracking with the intermediate to high scenario-based estimates.At the current 21st century rate of rise, water level elevations above the 2000 datum will reach +35 cm (Fernandina Beach, Florida) to +46 cm (Wilmington, North Carolina) by 2050.

Table 7 .
[44,116]ion-based linear extrapolations (this study, left column) and scenario-based median estimates of sea level elevation with likely ranges bracketed (i.e., 17th-83rd percentile[44]) in 2050 relative to a baseline of 2000.The observation-based extrapolations are derived from best-fit analysis of all data (black), three decades: 1993-2022 (red), and 21st century: 2003-2022 (purple).The two scenario-based estimates (columns) that bound the observation-based extrapolations (rows) are shown for each of the six tide gauge stations as vertical dividing lines using the observation-based color scheme (i.e., red = 1993-2022, purple = 2003-2022).Values expressed in mm.Scenario-based estimates obtained from Supplementary Data Files to NOAA's 2022 sea level rise technical report[44,116]).