Future coastal hazards along the U.S. North and South Carolina coasts

By Patrick L. Barnard, Kevin M. Befus, Jeffrey J. Danielson, Anita C. Engelstad, Li H. Erikson, Amy C. Foxgrover, Matthew W. Hardy, Daniel J. Hoover, Tim Leijnse, Chris Massey, Robert McCall, Norberto Nadal-Caraballo, Kees M. Nederhoff, Leonard Ohenhen, Andrea O'Neill, Kai A. Parker, Manoocher Shirzaei, Xin Su, Jennifer A. Thomas, Maarten van Ormondt, Sean F. Vitousek, Kilian Vos, and Madison C. Yawn

https://doi.org/10.5066/P9W91314


Dates

Published: Jan. 28, 2023

Summary

This product consists of several datasets that map future coastal flooding and erosion hazards due to sea level rise (SLR) and storms along the North and South Carolina coast. The SLR scenarios encompass a plausible range of projections by 2100 based on the best available, science and with enough resolution to support a suite of different planning horizons. The storm scenarios are derived with the use of atmospheric drivers from the Coupled Model Intercomparison Project Phase 6 (CMIP6) High Resolution Model Intercomparison Project (HighResMIP), a widely endorsed study that seeks to improve GCM hurricane modeling and applies a systematic approach to investigate the impact of horizontal resolution on simulating climate variables and hurricanes. Downscaled data products provided as part of this data release include: Projected Flood hazards (28 scenarios of 7 SLRs in combination with 3 storms and daily conditions), storm surge and astronomic tide time-series within the nearshore region for projected future storms (2020-2050), ​shoreline change time-series (2020-2100), depth to water table (7 SLR scenarios)​, and vertical land motion (14 years of historical observations).

Simulation Data

  • Projections of coastal flood hazards and flood potential for North Carolina and South Carolina

    Projected impacts by compound coastal flood hazards for future sea-level rise (SLR) and storm scenarios are shown for North Carolina and South Carolina. Accompanying uncertainty for each SLR and storm scenario, indicating total uncertainty from model processes and contributing datasets, are illustrated in maximum and minimum flood potential. As described by Nederhoff and others (2023), projections were made using a system of numerical models driven by output from Global Climate Models (GCMs) from the Coupled Model Intercomparison Project Phase 6 (CMIP6) and a tropical cyclone database from US Army Corp of Engineers. The resulting data products include detailed flood-hazard maps along the North Carolina and South Carolina coast due to sea level rise and plausible future storm conditions that consider the changing climate, hurricanes, and natural variability. In addition to sea-level rise, flood simulations run by these numerical models included dynamic contributions from tide, storm surge, wind, waves, river discharge, precipitation, and seasonal sea-level fluctuations. Outputs include impacts from combinations of SLR scenarios (0, 0.25, 0.5, 1.0, 1.5, 2.0, and 3.0 m) storm conditions including 1-year, 20-year and 100-year return interval storms and a background condition (no storm - astronomic tide and average atmospheric conditions).

  • Projections of coastal water depths for North Carolina and South Carolina

    Projected water depths from compound coastal flood hazards for future sea-level rise (SLR) and storm scenarios are shown for North Carolina and South Carolina. As described by Nederhoff and others (2023), projections were made using a system of numerical models driven by output from Global Climate Models (GCMs) from the Coupled Model Intercomparison Project Phase 6 (CMIP6) and a tropical cyclone database from U.S. Army Corp of Engineers. The resulting data are depths of projected flood hazards along the North Carolina and South Carolina coast due to sea level rise and plausible future storm conditions that consider the changing climate, hurricanes, and natural variability. The resulting data products include water depths that are consistent with coastal flood projections, also available in this dataset (Barnard, and others, 2023). In addition to sea-level rise, flood simulations run by these numerical models included dynamic contributions from tide, storm surge, wind, waves, river discharge, precipitation, and seasonal sea-level fluctuations. Outputs include impacts from combinations of SLR scenarios (0, 0.25, 0.5, 1.0, 1.5, 2.0, and 3.0 m), storm conditions including 1-year, 20-year, and 100-year return interval storms, and a background condition (no storm - astronomic tide and average atmospheric conditions).

  • Projected groundwater head along the North and South Carolina coasts

    Seamless unconfined groundwater heads for U.S. coastal North and South Carolina groundwater systems were modeled with homogeneous, steady-state MODFLOW simulations. The geographic extent examined was limited primarily to low-elevation (land surface less than approximately 10 m above mean sea level) areas. Steady-state MODFLOW groundwater flow models were used to obtain detailed (50-meter-scale) predictions over large geographic scales (100s of kilometers) of groundwater heads for both current and future sea-level rise (SLR) scenarios (0, 0.25, 0.5, 1, 1.5, 2, 2.5, and 3 m) using 3 spatially varying hydraulic conductivities (K); one based on published K’s, one with published K’s reduced by a factor of 10 (K*0.1), and one with published K’s increased by a factor of 10 (K*10) to assess the sensitivity of model results to K. All models had variable thicknesses corresponding to published transmissivities. The models were run with a local mean higher-high water (MHHW) marine boundary condition and with groundwater reaching the land surface removed from the model, simulating loss via natural drainage.

  • Projected water table depths along the North and South Carolina coasts

    To predict water table depths, seamless groundwater heads for unconfined coastal North and South Carolina groundwater systems were modeled with homogeneous, steady-state MODFLOW simulations. The geographic extent examined was limited primarily to low-elevation (land surface less than approximately 10 m above mean sea level) areas. Steady-state MODFLOW groundwater flow models were used to obtain detailed (50-meter-scale) predictions over large geographic scales (100s of kilometers) of groundwater heads for both current and future sea-level rise (SLR) scenarios (0, 0.25, 0.5, 1, 1.5, 2, 2.5, and 3 m) using 3 spatially varying hydraulic conductivities (K); one based on published K’s, one with published K’s reduced by a factor of 10 (K*0.1), and one with published K’s increased by a factor of 10 (K*10) to assess the sensitivity of model results to K. All models had variable thicknesses corresponding to published transmissivities. The models were run with a local mean higher-high water (MHHW) marine boundary condition, and with groundwater reaching the land surface removed from the model, simulating loss via natural drainage. Modeled groundwater heads were then subtracted from high-resolution topographic digital elevation model (DEM) data to obtain the water table depths.

  • Projected groundwater emergence and shoaling along the North and South Carolina coasts

    Groundwater emergence and shoaling extents are derived from water table depth GeoTIFFs, which are calculated as steady-state groundwater model heads subtracted from high-resolution topographic digital elevation model (DEM) land surface elevations. Results are provided as shapefiles of water table depth in specific depth ranges.

  • Nearshore water level, tide, and non-tidal residual hindcasts (1979-2016) for the North and South Carolina coasts

    A dataset of modeled nearshore water levels (WLs) was developed for the North and South Carolina coastlines. Water levels, defined for this dataset as the linear sum of tides and non-tidal residuals (NTR), were produced by Muis and others (2016) using a global tide and surge model (GTSM) forced by global atmospheric fields -. Water level outputs were extracted from the global grid at approximately 20 km resolution along the coastlines. These data were then statistically downscaled using a signal-specific set of corrections to improve skill in comparison to tide gauge observations (Parker and others, 2023). Hindcast water levels were forced by ERA5 atmospheric forcing provided by the dataset of Hersbach and others (2020). ERA5 is a reanalysis product, incorporating observations and data assimilation to best represent the experienced climate. Therefore, data from this version of the dataset are comparable to observed WLs along the study region.

  • Nearshore water level, tide, and non-tidal residual future projections (2016-2050) for the North and South Carolina coasts

    A dataset of modeled nearshore water levels (WLs) was developed for the North and South Carolina coastlines. Water levels, defined for this dataset as the linear sum of tides and non-tidal residuals (NTR), were produced by Muis and others (2016) using a global tide and surge model (GTSM) forced by global atmospheric fields. Water level outputs were extracted from the global grid at approximately 20 km resolution along the southeast Atlantic coastline. These data were then statistically downscaled using a signal-specific set of corrections to improve skill in comparison to tide gauge observations (Parker and others, 2023). Projected water levels were forced by CMIP6 future period data. This dataset provides information on how water levels are expected to change moving towards the future. Four CMIP6 climate models were selected from the High-Resolution Model Intercomparison project (highresMIP; Haarsma and others, 2016) to sample variability in climate predictions.

  • Nearshore parametric wave setup hindcast data (1979-2019) for the North and South Carolina coasts

    This dataset presents alongshore wave setup timeseries for the North and South Carolina coastlines. Wave setup was modelled using parameterization for open coast sandy beaches as presented in Stockdon and others (2006). The parameterization relates onshore wave setup to offshore wave conditions and beach characteristics. Wave conditions were extracted at approximately the 10 m depth contour and reverse shoaled to the deep-water condition. These data were then matched to cross-shore transects spaced at approximately 1 km along the coastline. Beach slope at these transects were extracted from lidar (Doran and others, 2017) and temporally averaged across all available datasets. This dataset is forced by waves taken directly from the ERA5 reanalysis, which incorporates observations and data assimilation (Hersbach and others, 2020). Therefore, data presented in this data release represent a best prediction of the observed historical wave setup along the study region. This dataset and the methodology used for its production is further detailed in Parker and others (2023). Modeled wave setup time series are presented for the hindcast period 1979 to 2019. Output includes 1-hour wave setup provided at approximately 1,600 alongshore transects at approximately 1-5 km resolution. Data are available as csv files for each transect location and are bundled by state.

  • Nearshore parametric wave setup future projections (2020-2050) for the North and South Carolina coasts

    This dataset presents alongshore wave setup timeseries for the North and South Carolina coastlines. Wave setup was modelled using parameterization for open coast sandy beaches as presented in Stockdon and others (2006). The parameterization relates onshore wave setup to offshore wave conditions and beach characteristics. Wave conditions were extracted at approximately the 10 m depth contour and reverse shoaled to the deep-water condition. These data were then matched to cross-shore transects spaced at approximately 1 km along the coastline. Beach slope at these transects were extracted from lidar (Doran and others, 2017) and temporally averaged across all available datasets. Waves were modelled using a global WAVEWATCHIII model forced by atmospheric forcing from the Coupled Model Intercomparison Project (CMIP6) future period data. This dataset provides information on how water setup is expected to change moving towards the future. The methodology used to produce this dataset is further detailed in Parker and others (2023). Data are provided for 6 CMIP6 models from the HighResMIP project (Haarsma and others, 2016). Output includes 1-hour wave setup provided at approximately 1,600 alongshore transects at approximately 1-5 km resolution. Data are available as csv files for each transect location and are bundled by state.

  • Projections of shoreline change of current and future (2005-2100) sea-level rise scenarios for North Carolina and South Carolina

    This dataset contains projections of shoreline change and uncertainty bands for future scenarios of sea-level rise (SLR). Scenarios include 25, 50, 75, 100, 150, 200, and 300 centimeters (cm) of SLR by the year 2100. Output for SLR of 0 cm is also included, reflective of conditions in 2005, in accordance with recent SLR projections and guidance from the National Oceanic and Atmospheric Administration (NOAA; see process steps). Projections were made using the Coastal Storm Modeling System - Coastal One-line Assimilated Simulation Tool (CoSMoS-COAST), a numerical model (described in Vitousek and others, 2017; 2021; 2023) run in an ensemble forced with global-to-local nested wave models and assimilated with satellite-derived shoreline (SDS) observations. Shoreline positions from models are generated at pre-determined cross-shore transects and output includes different cases covering important model behaviors (cases are described in process steps of metadata; see citations listed in the Cross References section for more details on the methodology and supporting information). This model shows change in shoreline positions along transects, considering sea level, wave conditions, along-shore/cross-shore sediment transport, long-term trends due to sediment supply, and estimated variability due to unresolved processes (as described in Vitousek and others, 2021). Variability associated with complex coastal processes (for example, beach cusps/undulations and shore-attached sandbars) are included via a noise parameter in a model, which is tuned using observations of shoreline change at each transect and run in an ensemble of 200 simulations; this approach allows for a representation of statistical variability in a model that is assimilated with sequences of noisy observations. The model synthesizes and improves upon numerous, well-established shoreline models in the scientific literature; processes and methods are described in this metadata (see lineage and process steps), but also described in more detail in Vitousek and others 2017, 2021 and 2023. KMZ data are readily viewable in Google Earth. For best display of results, it is recommended to turn off any 3D features or terrain. For technical users and researchers, shapefile and KMZ data can be ingested into geographic information system (GIS) software such as Global Mapper or QGIS.

Satellite Data

  • Satellite-derived shorelines for North Carolina and South Carolina (1984-2021)

    This dataset contains shoreline positions derived from available Landsat satellite imagery for North Carolina and South Carolina for the time period of 1984 to 2021. Positions were determined using CoastSat (Vos and others, 2019a and 2019b), an open-source mapping toolbox, was used to classify coastal Landsat imagery and detect shorelines at the sub-pixel scale.
    To understand shoreline evolution in complex environments and operate long-term simulations illustrating potential shoreline positions in the next century (Vitousek and others, 2017, 2021), robust historical shoreline data is necessary. Satellite-derived shorelines (SDS) offer expansive shoreline observational data over large geographic and temporal scales. Resulting shorelines for the period of 1984-2021 are presented in KMZ format. Significant uncertainty is associated with the locations of shorelines in extremely dynamic regions, including at the locations of river mouths, tidal inlets, capes, and ends of spits.
    These data are readily viewable in Google Earth. For best display of results, it is recommended to turn off any 3D viewing. For technical users and researchers, data can be ingested into Global Mapper or QGIS for more detailed analysis.

  • Vertical land motion rates for the years 2007 to 2020 for the North and South Carolina coasts

    Rates of land subsidence and uplift for the North and South Carolina coasts are derived from Sentinel-1A/B (2015-2020) and ALOS (2007-2011) synthetic aperture radar (SAR) satellites, at approximately 50-75 m resolution and mm-level precision. The data consist of vertical land motion (VLM) rates and the 1-sigma error in land motion rates and are available as csv files.

Suggested Citation

Barnard, P.L., Befus, K.M., Danielson, J.J., Engelstad, A.C., Erikson, L.H., Foxgrover, A.C., Hardy, M.W., Hoover, D.J., Leijnse, T., Massey, C., McCall, R., Nadal-Caraballo, N., Nederhoff, K.M., Ohenhen, L., O'Neill, A., Parker, K.A., Shirzaei, M., Su, X., Thomas, J.A., van Ormondt, M., Vitousek, S.F., Vox, K., and Yawn, M.C., 2023, Future coastal hazards along the U.S. North and South Carolina coasts: U.S. Geological Survey data release, https://doi.org/10.5066/P9W91314.