AgentSkillsCN

nws-flood-thresholds

从美国国家气象局(NWS)下载洪水等级阈值。适用于为USGS站点确定洪水水位、获取行动级/轻微级/中度级/重度级洪水阶段,或为各站点匹配对应的洪水阈值时使用。

SKILL.md
--- frontmatter
name: nws-flood-thresholds
description: Download flood stage thresholds from NWS (National Weather Service). Use when determining flood levels for USGS stations, accessing action/minor/moderate/major flood stages, or matching stations to their flood thresholds.
license: MIT

NWS Flood Thresholds Guide

Overview

The National Weather Service (NWS) maintains flood stage thresholds for thousands of stream gages across the United States. These thresholds define when water levels become hazardous.

Data Sources

Option 1: Bulk CSV Download (Recommended for Multiple Stations)

code
https://water.noaa.gov/resources/downloads/reports/nwps_all_gauges_report.csv

Option 2: Individual Station Pages

code
https://water.noaa.gov/gauges/<station_id>

Example: https://water.noaa.gov/gauges/04118105

Flood Stage Categories

CategoryCSV ColumnDescription
Action Stageaction stageWater level requiring monitoring, preparation may be needed
Flood Stage (Minor)flood stageMinimal property damage, some public threat. Use this to determine if flooding occurred.
Moderate Flood Stagemoderate flood stageStructure inundation, evacuations may be needed
Major Flood Stagemajor flood stageExtensive damage, significant evacuations required

For general flood detection, use the flood stage column as the threshold.

Downloading Bulk CSV

python
import pandas as pd
import csv
import urllib.request
import io

nws_url = "https://water.noaa.gov/resources/downloads/reports/nwps_all_gauges_report.csv"

response = urllib.request.urlopen(nws_url)
content = response.read().decode('utf-8')
reader = csv.reader(io.StringIO(content))
headers = next(reader)
data = [row[:43] for row in reader]  # Truncate to 43 columns
nws_df = pd.DataFrame(data, columns=headers)

Important: CSV Column Mismatch

The NWS CSV has a known issue: header row has 43 columns but data rows have 44 columns. Always truncate data rows to match header count:

python
data = [row[:43] for row in reader]

Key Columns

Column NameDescription
usgs idUSGS station ID (8-digit string)
location nameStation name/location
stateTwo-letter state code
action stageAction threshold (feet)
flood stageMinor flood threshold (feet)
moderate flood stageModerate flood threshold (feet)
major flood stageMajor flood threshold (feet)

Converting to Numeric

Threshold columns need conversion from strings:

python
nws_df['flood stage'] = pd.to_numeric(nws_df['flood stage'], errors='coerce')

Filtering by State

python
# Get stations for a specific state
state_stations = nws_df[
    (nws_df['state'] == '<STATE_CODE>') &
    (nws_df['usgs id'].notna()) &
    (nws_df['usgs id'] != '') &
    (nws_df['flood stage'].notna()) &
    (nws_df['flood stage'] != -9999)
]

Matching Thresholds to Station IDs

python
# Build a dictionary of station thresholds
station_ids = ['<id_1>', '<id_2>', '<id_3>']
thresholds = {}

for _, row in nws_df.iterrows():
    usgs_id = str(row['usgs id']).strip()
    if usgs_id in station_ids:
        thresholds[usgs_id] = {
            'name': row['location name'],
            'flood': row['flood stage']
        }

Common Issues

IssueCauseSolution
Column mismatch errorCSV has 44 data columns but 43 headersTruncate rows to 43 columns
Missing thresholdsStation not in NWS databaseSkip station or use alternative source
Value is -9999No threshold definedFilter out these values
Empty usgs idNWS-only stationFilter by usgs id != ''

Best Practices

  • Always truncate CSV rows to match header count
  • Convert threshold columns to numeric before comparison
  • Filter out -9999 values (indicates no threshold defined)
  • Match stations by USGS ID (8-digit string with leading zeros)
  • Some stations may have flood stage but not action/moderate/major