DOI: 10.5281/zenodo.5951626
title: Interagency report: Global and Regional Sea Level Rise Scenarios for the United States: Updated Mean Projections and Extreme Water Level Probabilities Along U.S. Coastlines

description: <strong>Code and data for Section 2 of the Interagency report: Global and Regional Sea Level Rise Scenarios for the United States: Updated Mean Projections and Extreme Water Level Probabilities Along U.S. Coastlines</strong> This repository contains the code and data needed to produce the trajectories, projections, and observations for the Interagency report: Global and Regional Sea Level Rise Scenarios for the United States: Updated Mean Projections and Extreme Water Level Probabilities Along U.S. Coastlines. The report can be found on https://oceanservice.noaa.gov/hazards/sealevelrise/sealevelrise-tech-report-sections.html An interactive tool to study the observations, trajectories, and scenarios can be accessed from https://sealevel.nasa.gov/task-force-scenario-tool Frequently-asked questions: https://sealevel.nasa.gov/faq/16/ <strong>Authors</strong> William V. Sweet, NOAA National Ocean Service Benjamin D. Hamlington, NASA Jet Propulsion Laboratory Robert E. Kopp, Rutgers University Christopher P. Weaver, U.S. Environmental Protection Agency Patrick L. Barnard, U.S. Geological Survey Michael Craghan, U.S. Environmental Protection Agency Gregory Dusek, NOAA National Ocean Service Thomas Frederikse, NASA Jet Propulsion Laboratory Gregory Garner, Rutgers University Ayesha S. Genz, University of Hawai‘i at Mānoa, Cooperative Institute for Marine and Atmospheric Research John P. Krasting, NOAA Geophysical Fluid Dynamics Laboratory Eric Larour, NASA Jet Propulsion Laboratory Doug Marcy, NOAA National Ocean Service John J. Marra, NOAA National Centers for Environmental Information Jayantha Obeysekera, Florida International University Mark Osler, NOAA National Ocean Service Matthew Pendleton, Lynker Daniel Roman, NOAA National Ocean Service Lauren Schmied, FEMA Risk Management Directorate William C. Veatch, U.S. Army Corps of Engineers Kathleen D. White, U.S. Department of Defense Casey Zuzak, FEMA Risk Management Directorate <strong>Contents</strong> This data and code set contains the following directories: <em>Results</em> The <code>Results</code> folder contains the resulting projections, trajectories and observations from the report. <code>TR_global_projections.nc</code>: GMSL projections, trajectory, and observations <code>TR_regional_projections.nc</code>: Regional observations, projections and trajectories <code>TR_local_projections.nc</code>: Local observations, projections and trajectories <code>TR_gridded_projections.nc</code>: Gridded projections These files are in the NetCDF forrmat. To read the NetCDF files, many free software packages are available, including ncview and Panoply. Free NetCDF packages are available to directly import the data into Julia and Python code. <em>Code</em> The <code>Code</code> folder contains all the computer code used to read and analyze the observations and the projections, and to generate the trajectories. To run this code, you need Julia. The code requires the Julia packages <code>CSV</code>, <code>Interpolations</code>, <code>JSON</code>, <code>LoopVectorization</code>, <code>MAT</code>, <code>NCDatasets</code>, <code>NetCDF</code>, <code>Plots</code>, <code>XLSX</code>, <code>LinearAlgebra</code>, and <code>Statistics</code>. They can be installed by pressing <code>]</code> at the Julia REPL and typing: <pre><code>add CSV Interpolations JSON LoopVectorization MAT NCDatasets NetCDF Plots XLSX LinearAlgebra Statistics </code></pre> This program also requires Hector. Hector needs to be installed or compiled. In the file <code>Hector.jl</code> update the path to the Hector executable on lines 30 and 104. Run <code>Run_TR.jl</code> in the REPL or run <code>julia Run_TR.jl</code> from the command line to run the projections. The projections are then written to the <code>.\Data</code> directory. The folder contains the following files: <code>Run_TR.jl</code>: This is the main routine that (eventually) calls all the functions to compute the projections. <code>ConvertNCA5ToGrid.jl</code>: Converts the original NCA5 projections to a set of netCDF files that's used throughout this code <code>ProcessObservations.jl</code>: Reads and processes the tide-gauge and altimetry observations <code>GlobalProjections.jl</code>: Reads and processes the GMSL observations and projections, and computes the trajectory <code>RegionalProjections.jl</code>: Reads and processes the regional projections and computes the trajectories <code>LocalProjections.jl</code>: Reads and processes the local projections at the tide-gauge locations and computes the trajectories <code>GriddedProjections.jl</code>: Reads the gridded NCA5 projections and add a GMSL baseline correction for the 2005 vs 2000 baseline <code>SaveFigureData.jl</code>: Reads the results and writes text files for GMT <code>Hector.jl</code>: Wrapper for Hector, used to compute trends and uncertainties. <code>Masks.jl</code>: Defines the region masks for each region. <em>Data</em> The <code>Data</code> directory contains the input data sets used during the computations. Please appropriately cite the input data if you use it. It contains the following: Directories: <code>ClimIdx</code>: Map with climate indices (NAO, PDO, MEI) used to remove internal variability. All the indices come from NOAA Physical Sciences Laboratory (PSL) and NOAA Climate Prediction Centre (CPC) <code>NCA5_projections</code> Contains the NCA5 projections for each scenario (Low, IntLow, Int, IntHigh, and High). For each scenario, the GMSL projections, projections at tide-gauge locations and on a 1-degree grid are provided. Files: <code>basin_codes.nc</code>: Map with basin codes. from Eric Leuliette/NOAA. Data provided by the NOAA Laboratory for Satellite Altimetry. <code>CDS_monthly_1993_2020.nc</code>: Monthly-mean sea level (1993-2020) from gridded altimetry. Obtained from Copernicus Climate Data Store. This dataset contains modified Copernicus Climate Change Service information [2020] <code>enso_correction.mat</code>: GMSL correction for ENSO/PDO from Hamlington, B. D., Frederikse, T., Nerem, R. S., Fasullo, J. T., & Adhikari, S. (2020). Investigating the Acceleration of Regional Sea‐level Rise During the Satellite Altimeter Era. Geophysical Research Letters. https://doi.org/10.1029/2019GL086528 <code>filelist_psmsl.txt</code>: List with PSMSL file names and PSMSL IDs. Obtained from the Permanent Service for Mean Sea Level (PSMSL), 2021, Retrieved 29 Nov 2021. Simon J. Holgate, Andrew Matthews, Philip L. Woodworth, Lesley J. Rickards, Mark E. Tamisiea, Elizabeth Bradshaw, Peter R. Foden, Kathleen M. Gordon, Svetlana Jevrejeva, and Jeff Pugh (2013) New Data Systems and Products at the Permanent Service for Mean Sea Level. Journal of Coastal Research: Volume 29, Issue 3: pp. 493 – 504. https://doi.org/:10.2112/JCOASTRES-D-12-00175.1. <code>GEBCO_bathymetry_05.nc</code>: Bathymetry map of the global oceans from the General Bathymetric Chart of the Oceans (GEBCO). Source: GEBCO Compilation Group (2021) GEBCO 2021 Grid (<code>doi:10.5285/c6612cbe-50b3-0cff-e053-6c86abc09f8f</code>) The source data have been re-gridded onto a 0.5 degree grid. <code>GIA_Caron_stats_05.nc</code>: Glacial Isostatic Adjustment estimates from Caron, L., Ivins, E. R., Larour, E., Adhikari, S., Nilsson, J., & Blewitt, G. (2018). GIA Model Statistics for GRACE Hydrology, Cryosphere, and Ocean Science. Geophysical Research Letters, 45(5), 2203–2212. https://doi.org/10.1002/2017GL076644. The source data have been re-gridded onto a 0.5 degree grid. <code>global_timeseries_measures.nc</code>: Time series of estimated 20th-century GMSL and its components, based on Frederikse, T., Landerer, F., Caron, L., Adhikari, S., Parkes, D., Humphrey, V. W., Dangendorf, S., Hogarth, P., Zanna, L., Cheng, L., & Wu, Y.-H. (2020). The causes of sea-level rise since 1900. Nature, 584(7821), 393–397. https://doi.org/10.1038/s41586-020-2591-3 <code>GMSL_ensembles.nc</code>: Ensemble GMSL reconstruction from tide-gauges based on Frederikse, T., Landerer, F., Caron, L., Adhikari, S., Parkes, D., Humphrey, V. W., Dangendorf, S., Hogarth, P., Zanna, L., Cheng, L., & Wu, Y.-H. (2020). The causes of sea-level rise since 1900. Nature, 584(7821), 393–397. https://doi.org/10.1038/s41586-020-2591-3 <code>GMSL_TPJAOS_5.0_199209_202106.txt</code>: Global Mean Sea Level Trend from Integrated Multi-Mission Ocean Altimeters TOPEX/Poseidon, Jason-1, OSTM/Jason-2, and Jason-3 Version 5.1 [Data set]. NASA Physical Oceanography DAAC. https://doi.org/10.5067/GMSLM-TJ151. This altimetry dataset uses the methods as described in Beckley, B. D., Callahan, P. S., Hancock, D. W., Mitchum, G. T., & Ray, R. D. (2017). On the “Cal-Mode” Correction to TOPEX Satellite Altimetry and Its Effect on the Global Mean Sea Level Time Series. Journal of Geophysical Research: Oceans, 122(11), 8371–8384. https://doi.org/10.1002/2017JC013090 <code>grd_1992_2020.nc</code>: Seafloor deformation due to contemporary GRD effects based on Frederikse, T., Landerer, F., Caron, L., Adhikari, S., Parkes, D., Humphrey, V. W., Dangendorf, S., Hogarth, P., Zanna, L., Cheng, L., & Wu, Y.-H. (2020). The causes of sea-level rise since 1900. Nature, 584(7821), 393–397. https://doi.org/10.1038/s41586-020-2591-3 <code>region_mask.nc</code>: Mask with the definition of all regions. <code>US_tg_monthly.xlsx</code>: Tide gauge observations from the NOAA tide gauge network <em>GMT</em> This directory contains the GMT scripts to make Figures 1.2, 2.1, 2.2, 2.6, and A.1.2 from the report. To generate the figures, make sure GMT is installed and run the Shell script in each directory.