Astronomical variations in tidal magnitude can strongly modulate the severity of coastal flooding on daily, monthly, and interannual timescales. Here, we present a new quasi-nonstationary skew surge joint probability method (qn-SSJPM) that estimates interannual fluctuations in flood hazard caused by the 18.6 and quasi 4.4-year modulations of tides. We demonstrate that qn-SSJPM-derived storm tide frequency estimates are more precise and stable compared with the standard practice of fitting an extreme value distribution to measured storm tides, which is often biased by the largest few events within the observational period. Applying the qn-SSJPM in the Gulf of Maine, we find significant tidal forcing of winter storm season flood hazard by the 18.6-year nodal cycle, whereas 4.4-year modulations and a secular trend in tides are small compared to interannual variation and long-term trends in sea-level. The nodal cycle forces decadal oscillations in the 1% annual chance storm tide at an average rate of ±13.5 mm/y in Eastport, ME; ±4.0 mm/y in Portland, ME; and ±5.9 mm/y in Boston, MA. Currently (in 2020), nodal forcing is counteracting the sea-level rise-induced increase in flood hazard; however, in 2025, the nodal cycle will reach a minimum and then begin to accelerate flood hazard increase as it moves toward its maximum phase over the subsequent decade. Along the world’s meso-to-macrotidal coastlines, it is therefore critical to consider both sea-level rise and tidal non-stationarity in planning for the transition to chronic flooding that will be driven by sea-level rise in many regions over the next century.

Of the major coastal land change mechanisms responsible for relative sea-level change, tectonic subsidence is generally quoted as ranging from < mm/yr to 1 cm/yr. However, we documented coseismic and ongoing post-earthquake surface displacements from continuous GPS and tide gauge/altimetry data that indicated rapid subsidence on two of the major Samoan Islands of 12 - 20 cm during and following the 8.1 2009 Tonga-Samoa earthquake. Earlier results and our modeling of GRACE-derived gravimetric data provided a preliminary forecast of future relative sea-level rise through rapid land subsidence [Han et al., 2019]. Of course these numerical forecasts of time-dependent deformation are only as good as our input observations and our assumed rheological models. As part of our current NASA Earth Surface and Interior study, we are obtaining a wider range of data to constrain and test alternate models of ongoing postseismic deformation across American Samoa and Upolu, Samoa: (1) times series of altimetry plus tide gauge data processed to complement the cGPS data available to provide high-temporal resolution, point measurements of uplift/subsidence, (2) InSAR derived observations of surface deformation across the highly vegetated Samoan Islands, (3) evaluating and using NASA satellite lidar data (ICESat-I & ICESat-II, GEDI) for fusion with multi-source topographic data sets and for estimating topographic change on the decadal time scale. We are evaluating and using these new observations to better understand and separate out local, island-wide, and multi-island subsidence patterns and to evaluate the high impact of rising sea-level in a tectonically active region.