:: Back to the Main Experimental STOFS Portal Page ::

Table of Contents


STOFS Disclaimer

This is an experimental web portal for graphic visualization of the operational and experimental results from the storm surge and tide forecast systems being developed and tested by the National Ocean Service.

STOFS model output is NOT total water level guidance. Actual water levels can be significantly higher than forecast due to waves and other components not presently included in STOFS.

Please check with your regional National Weather Service forecast office for the official water level forecast.




Page Description: STOFS-2D-Global

This page provides experimental web dissemination for the output and validation of the Global Surge and Tide Operational Forecast System 2-D (STOFS-2D-Global).

The Coastal & Marine Modeling Branch (CMMB) of the Coast Survey Development Laboratory (CSDL) has upgraded STOFS-2D-Global to v3.1. This STOFS-2D-Global upgrade includes:

*Inclusion of a gridded bias correction. This bias correction uses the current station bias correction approach used in operational STOFS2D-Global (v2.1), which corrects the station water level forecast guidance biases by removing the bias between the observed water levels at NOS/CO-OPS stations and most recent 5 days of nowcast water levels. This bias correction at stations is then propagated across the entire global grid using a graph diffusion approach. The result is that both station and gridded output is bias corrected.

*Addition of a few stations to the station output, as requested by National Weather Service forecast offices: Bogue Sound at Emerald Isle, NC; White Oak River at Swansboro, NC; Kings Bay MSF Pier, GA; and TEMCO Kalama Terminal, WA. Also, some redundant stations are removed from the output: Duck, NC; Fort Point, NH; Money Point, VA; and Springmaid Pier, SC. South Amboy, NJ was removed as it is no longer active and S Amboy Raritan River, NJ is the active station in that area. Finally, Naples Bay, North, FL has an updated NOS ID from 8725110 to 8725114.

The user evaluation and feedback for STOFS-2D-Global was officially collected up until January 16 2026, but we will always welcome feedback. The STOFS-2D-Global Upgrade Briefing to the NCEP Director occurred on January 22 2026. Operational implementation of the STOFS-2D-Global upgrade to v3.1 occurred at the 12 UTC cycle on August 17 2026.

Evaluation is always being sought to verify that the upgrade of STOFS-2D-Global to v3.1 improves or maintains skill in the forecast guidance it provides for your use.

Please send your feedback to:
Greg Seroka [email]
Edward Myers [email]  
Saeed Moghimi [email]
 

Operational STOFS-2D-Global Configuration

Model: ADCIRC v56, 2D depth-integrated, with internal tides, self-attraction & loading. Spherical coordinates on an unstructured global mesh

Forcing
: GFS-FV3 13km 10-m winds, mslp, seaice, tidal potential with short-period 13 harmonics

Mesh
: Global. Node count: 12,785,004. Elements count: 24,875,336. Coastal resolution: at least 1.5 km globally. Nominal coastal resolution (US): 120m. Maximal coastal resolution: 80m (30m Chetco River, OR). Open ocean resolution: 40km. Upland boundary: +10m MSL along the US East Coast, up to 20m for some of the US Pacific Islands. Flood plain spatial resolution: nominal 500m for the US East Coast, up to 80m for some of the US Islands.

Computation
: Computational time step: 6 seconds (implicit). In total, STOFS-2D-Global currently runs on 4,064 computational cores on Weather and Climate Operational Supercomputing System 2 (WCOSS2).



Model Output

Maximal forecasted water level elevation (maxele) on a native netCDF, 6-min forecast water level timeseries at 1688 point locations worldwide, hourly snapshot fields in GRIB2 format on NDFD subgrids (CONUS, North Pacific, Alaska, Hawaii, Guam, Puerto Rico).

Digital Output:
netCDF and GRIB2 output is available via
https://nomads.ncep.noaa.gov/pub/data/nccf/com/stofs/prod/.

Graphics Output:
this page provides links to reports including map plots of the maximal forecasted water level elevation (maxele), map plots of the maximal forecasted GFS winds, timeseries plots of the forecasted combined water levels along with the observed water levels where available, average and individual skill metrics for the point outputs where the observations are available. The reports are generated for operational STOFS-2D-Global.

The typical auto-generated model report provides the following information:
Name (tag) of the experiment. The tag includes datestamp of the forecast cycle. It is printed on each plot. Please verify that you are looking at the latest cycle by checking the tag.
Datestamp when the report was generated (UTC)
Date span for the statistical analysis (UTC)
Bounding box of the geographical coordinates (East longitudes/North latitudes)
The section Maximal Forecasted Water Levels containing the maps of maxele fields for the whole domain and key sub-domains.
The section Time-Series Metrics containing the output and the analysis of the water level point output.


Skill Reporting

Skill metrics are reported for the time series comparison of the point output against NOAA coastal tide gauges water level observations, where available. The below describes the different metrics. For more details, please see https://github.com/noaa-ocs-modeling/autoval#package-description.

Metrics include:
Skill (units)
Description
Best Skill Value
RMSE (meters) Root Mean Square Error between the model and the observations. RMSE=0
PEAK (meters) Under/overestimation of the maximal water level.  PEAK=0
PLAG (minutes) Time lag between the modeled and the observed peak in water level. PLAG=0
BIAS (meters) Linear bias in the modeled water level. BIAS=0
VEXP (%) Variance explained, a measure of a coherence between the model and the observations. VEXP=100
SKIL (unitless) Statistical Skill of the model against the observations. SKIL=1
RVAL (unitless) R-Value of the model against the observations. RVAL=1
NPTS (unitless) The number of 6-minute model/data pairs at the location that went into computing the above metrics.

Metrics are computed for each individual point output station. All point output is provided in the Individual Time-Series Statistics sub-section of the report.

Below is the example of the time series analysis for a given station, showing the modeled (blue) and observed (green) hydrographs, and the corresponding metrics:

ts

Individual station analysis provides time series skill for the model on the date span where the results overlap the observations. On the example above, the statistics were computed over the 17.8 hours (178 6-minute points where both model and data are valid). Skill metrics is represented as both a text table and a dashboard, with acceptable ranges for each skill shown in green bars.

The maps for each skill metrics of all valid individual stations are provided in the experiment report in the
Average Time-Series Statistics sub-section of the report.

The example below shows the map of RMSE values at each valid station. Values below a 0.2m threshold are plotted as circles. Values above 0.2m are plotted as upward-looking triangles, to help identify areas where the error is larger than a given threshold.   

rmse


Stations without the observations only report modeled nowcast/forecast time series (blue curve).

Average skill metrics for the experiment is shown in the table in the Average Time-Series Statistics sub-section of the report.



STOFS Pre- and Post-Processing Software

The pre- and post-processing software including skill computations has been developed by CSDL, in collaboration with other partners. The codes are written in Python3, and are publicly available from GitHub.

The baseline library csdllib.
The automated skill analysis application autoval.


Most Recent Model Output



 STOFS-2D-Global
[Report]
operational
* STOFS-2D-Global operational run output is broadcast to nomads server.


Maximal forecasted GFS winds for the latest operational STOFS-2D-Global cycle (in knots):

maxwvel