USGS Streamgages Linked to the Medium Resolution NHD

David W. Stewart, Alan Rea and David M. Wolock April 2006 USGS Streamgages Linked to the Medium Resolution NHD vector digital data Data Series DS-195 Reston, Va U.S. Geological Survey http://water.usgs.gov/GIS/metadata/usgswrd/XML/Streamgage.xml USGS_Streamgages-NHD_Locations The locations of approximately 23,000 current and historical U.S. Geological Survey (USGS) streamgages in the United States and Puerto Rico (with the exception of Alaska) have been snapped to the medium resolution National Hydrography Dataset (NHD). The NHD contains geospatial information about mapped surface-water features, such as streams, lakes, and reservoirs, etc., creating a hydrologic network that can be used to determine what is upstream or downstream from a point of interest on the NHD network. An automated snapping process made the initial determination of the NHD location of each streamgage. These initial NHD locations were comprehensively reviewed by local USGS personnel to ensure that streamgages were snapped to the correct NHD reaches. About 75 percent of the streamgages snapped to the appropriate NHD reach location initially and 25 percent required adjustment and relocation. This process resulted in approximately 23,000 gages being successfully snapped to the NHD. This dataset contains the latitude and longitude coordinates of the point on the NHD to which the streamgage is snapped and the location of the gage house for each streamgage. A process known as indexing may be used to create reference points (event tables) to the NHD reaches, expressed as a reach code and measure (distance along the reach). Indexing is dependent on the version of NHD to which the indexing is referenced. These data are well suited for use in indexing because nearly all the streamgage NHD locations have been reviewed and adjusted if necessary, to ensure they will index to the appropriate NHD reach. Flow characteristics were computed from the daily streamflow data recorded at each streamgage for the period of record. The flow characteristics associated with each streamgage include: First date (year, month, day) of streamflow data Last date (year, month, day) of streamflow data Number of days of streamflow data Number of days of non-zero streamflow data Minimum and maximum daily flow for the period of record (cubic feet per second) Percentiles (1, 5, 10, 20, 25, 50, 75, 80, 90, 95, 99) of daily flow for the period of record (cubic feet per second) Average and standard deviation of daily flow for the period of record (cubic feet per second) Mean annual base-flow index (BFI) computed for the period of record (fraction, ranging from 0 to 1) Year-to-year standard deviation of the annual base-flow index computed for the period of record (fraction) Number of years of data used to compute the base-flow index (years) The streamflow data used to compute flow characteristics were copied from the NWIS-Web historical daily discharge archive (nadww01.er.usgs.gov:/www/htdocs/nwisweb/data/discharge) on June 15, 2005. The purpose of this dataset is to document the location of USGS active and historical streamgages on the NHD version 2002 medium resolution (1:100,000-scale). This dataset can also be used to index these USGS streamgages to other versions of the NHD as well. The flow characteristics tied to each of the streamgages attached to the NHD will greatly facilitate hydrologic analysis and can be used for such things as the development of regional statistical models of streamflow at ungaged locations. en 2005 Source information retrieved from NWIS April 2004 Complete As needed -159.750305 -65.302467 48.999947 17.956547 -159.750305 -65.302467 17.956547 48.999947 Gages ISO 19115 Topic Category Streamgages Gaging Stations USGS Streamgage Network Locations for NHD Version June 2002 Medium Resolution Streamflow statistics Streamflow data Base flow index Hydrography National Hydrography Dataset NHD Streamgages inlandWaters United States of America United States USA Alabama Arizona Arkansas California Colorado Connecticut Deleware District of Columbia Florida Georgia Hawaii Idaho Illinois Indiana Iowa Kansas Kentucky Louisiana Maine Maryland Massachusetts Michigan Minnesota Mississippi Missouri Montana Nebraska Nevada New Hampshire New Jersey New Mexico New York North Carolina North Dakota Ohio Oklahoma Oregon Pennsylvania Rhode Island South Carolina South Dakota Tennessee Texas Utah Vermont Virginia Washington West Virginia Wisconsin Wyoming Puerto Rico None None US Gelogical Survey mailing address

437 USGS National Center

Reston Va. 22092 USA http://answers.usgs.gov/cgi-bin/gsanswers David W. Stewart, Alan Rea and David M. Wolock None Unclassified None ESRI ArcCatalog 9.1.0.722 Shapefile none none The positional accuracy of the NHD locations is very good. Most locations were reviewed by USGS Water Science Center personnel who knew the precise locations for the vast majority of all the streamgages. These reviewed locations have the REVIEWED field set to 'Y', which stands for "yes, these locations have been reviewed." In general, the positional accuracy of the active streamgage location and streamgages that were active in the past 30 years could be considered slightly more accurate then the older historical streamgages. Most NHD locations were reviewed by USGS Water Science Center personnel. The NHD locations of 17,482 are within a few meters of the delineated stream reach, since they were automatically snapped to the nearest latitude/longitude vertices in the NHD arc. The remaining 4,271 NHD locations were manually adjusted during the review and moved to the correct locations. Although the new locations were not snapped exactly to NHD vertices, reviewers were instructed to place the point so that it would be close enough to snap to the proper position on the NHD. The NHD locations not reviewed by USGS Water Science Center personnel are also included. These locations have the REVIEWED field set to 'N' which stands for "not reviewed. These locations were snapped to the NHD in an automated process, but no review has been done to ensure they snapped to the proper NHD reach. A total of 2,042 gages in the dataset were not reviewed." U.S. Geological Survey NWIS-Web streamflow data Digital Data 20050615 ground condition USGS-NWIS Point locations and hydrologic data Creating the streamgage dataset: The NWISWeb database was used exclusively to create the initial streamgage dataset. NWISWeb consolidates all the needed data from all USGS Water Science Center (WSC) NWIS databases in one place, making it easy to use without having to query each WSC's NWIS database individually. NWISWeb only contains information from NWIS db_01. Stations in other NWIS databases were not considered and are not part of the streamgage dataset. There are a total of 25,768 stations in the streamgage dataset for areas that are covered by NHD, which are, the conterminous US, plus Hawaii and Puerto Rico. Stations in Alaska, US Virgin Islands, Canada, and other US territories are not included because these areas do not have NHD coverage, at the time of the indexing. The following outlines how the streamgage dataset was constructed: 1. Retrieved all daily-values (DV) for the entire period of record from the NWISWeb database on April 20, 2004. The results are in 48 Water Science Center-specific files. 2. Each district DV file was examined. Stations containing at least 1 day of data associated with one or more of the following parameters were flagged: - parm = 00060 - Streamflow, ft/s - parm = 00065 - Gage height, ft - parm = 00054 - Reservoir storage, acre-ft - parm = 00053 - Surface Area, acres - parm = 00055 - Stream velocity, ft/s - parm = 00059 - Streamflow rate, gal/min - parm = 00062 - Reservoir elevation, ft - parm = 50042 - Discharge, gallons per minute - parm = 50051 - Flow rate, instantaneous, million gallons per day - parm = 72022 - Reservoir contents, M/gal - parm = 99020 - Elevation above NGVD 1929, ft - parm = 99060 - Streamflow, m/s - parm = 99065 - Gage height, m 3. The flagged stations were appended and duplicate stations removed. 4. Stations with data in water-years 2003 or 2004 were identified as "active". (Not all daily-values for the current water-year (2004) are in the database, and that is why water-year 2003 data was also considered. This means a few stations that were recently discontinued will be incorrectly flagged as "active".) 5. All 25,768 stations in the streamgage dataset were batch indexed using the US Environmental Protection Agency's (EPA) medium-resolution NHD database. 6. EPA provided the results of the indexing, containing a latitude/longitude coordinates of point location on the NHD stream that is associated with each station; known as an "event point". This "event point" and the latitude/longitude coordinate of the station were converted to the Albers Equal Area projection so a distance, in meters, between the station and the "event point" could be computed and added as an attribute. 7. The final streamgage dataset is in geographic coordinates using the NAD83 datum. Water Science Center (WSC) personnel then reviewed the batch-indexing results using a provided streamgage dataset and GIS/ArcMap custom application. This process not only tied USGS streamgages to the NHD, but also represents an important opportunity to use an efficient process to correct station location errors in NWIS. After the review was completed, the results were then consolidated creating this product. The reviewers examined 23,945 streamgages and found - 17,482 to be batch indexed correctly and the event point was not moved. ~ 73%. - 962 gages were flagged as "Do Not Use" and removed from this dataset, gages that do not contain any meaningful information. ~ 4% - 1380 gages were flagged as "No 100k NHD Reach" and removed from this dataset ~ 5% - 4271 gages have modified/moved event points - moved, on average, 795 meters. ~ 18% - 4111 gages have modified/moved site locations (gage house) - moved, on average, 1123 meters ~ 17% 8. The non-reviewed streamgages were added to the final dataset and these streamgages were flagged as "not reviewed" and the reviewed streamgages flagged as "reviewed". 9. The streamgages that were not included in the initial streamgages database but had streamflow characteristics were also added to the final streamgage database. The nearest NHD reach was automatically determined and the nearest point on the nearest NHD reach was calculated. This nearest NHD point became the NHD location. These streamgages have not been reviewed. The final streamgage database contains the location of 23,426 streamgages. Started with 25,768 streamgages, 962 were removed because they were identified during the review as "do not use". Another 1380 streamgages were removed because they were located on streams not delineated on the 100K NHD. Of the total 23,426 streamgages, 21,384 were reviewed and 2,042 were not. Estimating the flow characteristics: 1. Daily streamflow data for all streamgages were copied from the files located in the subdirectories of nadww01.er.usgs.gov:/www/htdocs/nwisweb/data/discharge on June 15, 2005. 2. The Perl program qstats.pl (Xiaodong Jian, USGS, written commun., 2005) was used to compute the time period and flow statistics for each streamgage. 3. Base-flow index estimation: A base-flow index (BFI) was computed for the streamgages using a Fortran program written by Tony Wahl (Bureau of Reclamation. U.S. Department of the Interior) and Ken Wahl (U.S. Geological Survey) (http://www.usbr.gov/pmts/hydraulics_lab/twahl/bfi/index.html). The computer program name is bfi41.f and was the most current version available at the time. The computer program was compiled on a Sun Unix system and run with the default parameter values. The BFI Web page (http://www.usbr.gov/pmts/hydraulics_lab/twahl/bfi/index.html) states: "The BFI program was developed to make the base-flow separation process less tedious and more objective. The program implements a deterministic procedure proposed in 1980 by the British Institute of Hydrology. The method combines a local minimums approach with a recession slope test. The program estimates the annual base-flow volume of unregulated rivers and streams and computes an annual base-flow index (BFI, the ratio of base flow to total flow volume for a given year) for multiple years of data at one or more gage sites. Although the method may not yield the true base flow as might be determined by a more sophisticated analysis, the index has been found to be consistent and indicative of base flow, and thus may be useful for analysis of long term base-flow trends. Users should be very cautious about using methods such as this for short-term storm events or for locations where streamflow is affected by upstream regulation, such as reservoir releases. In general, the method interprets most regulated releases as base flow. If the program is used for regulated streams, the effects of regulation must be carefully accounted for through manual adjustment of the program output." BFI references: Wahl, K.L., and Wahl, T.L., 1995, Determining the flow of Comal Springs at New Braunfels, Texas, In Proceedings of Texas Water '95, American Society of Civil Engineers, August 16-17, 1995 San Antonio, Texas: American Society of Civil Engineers, p. 77-86, information available on the World Wide Web, accessed March 25, 2003 , at URL http://www.usbr.gov/pmts/hydraulics_lab/twahl/bfi/texaswater95/comalsprings.html Wahl, K.L., and Wahl, T.L., 1988, Effects of regional ground-water declines on streamflows in the Oklahoma Panhandle, In Proceedings of Symposium on Water-Use Data for Water Resources Management: Tucson, Arizona, American Water Resources Association, p. 239-249, information available on the World Wide Web, accessed March 25, 2003, at URL http://www.usbr.gov/pmts/hydraulics_lab/twahl/bfi/bfi_beaver_river.pdf The use of firm, trade, and brand names is for identification purposes only and does not constitute endorsement by the U.S. government. Unknown Metadata imported. withheld Metadata imported. withheld Metadata imported. withheld Metadata imported. withheld Metadata imported. withheld Metadata imported. withheld 20060613 10222300 Dataset copied. withheld 20060613 10340700 Vector Simple FALSE 8 TRUE FALSE North American Datum of 1983 Geodetic Reference System 80 6378137.000000 298.257222 GCS_North_American_1983 0.000010 Explicit elevation coordinate included with horizontal coordinates Streamgages_USGS Feature Class 8 GOTQ Flag indicating flow data (1) or no flow data (0) OID 4 OBJECTID 0 0 Esri Sequential unique whole numbers that are automatically generated. NWISWEB URL to the station's NWISWeb home page. String 75 NWISWEB 0 0 Esri Sequential unique whole numbers that are automatically generated. HUC Hydrologic Cataloging Unit (8 digit HUC) of the gage String 8 HUC 0 0 Esri Sequential unique whole numbers that are automatically generated. HUC_REG Hydrologic Region (2 digit HUC) of the gage String 2 HUC_REG 0 0 Esri Sequential unique whole numbers that are automatically generated. SUBREGIONS Hydrologic sub-region (4 digit HUC) of the gage String 4 SUBREGIONS 0 0 Esri Sequential unique whole numbers that are automatically generated. FID_SITES NHD2GAGE_D Double 8 Internal Identifer 0 0 Esri Sequential unique whole numbers that are automatically generated. NHD2GAGE_D ACCOUNTING String 6 0 0 Distance between the NHD location and the gage house location, in meters Esri Sequential unique whole numbers that are automatically generated. ACCOUNTING Hydrologic accounting unit (6 digit HUC) of the gage Shape Geometry 0 0 0 Esri Coordinates defining the features. Shape Feature geometry. ESRI Coordinates defining the features. FID_SITES Integer 4 0 0 FID_sites Internal identifier String 5 AGENCY_CD 0 0 AGENCY_CD Gov. Agency responsible for the streamgage String 15 SITE_NO 0 0 SITE_NO USGS Station/Site identification number String 60 STATION_NM 0 0 STATION_NM Station name String 2 STATE_CD 0 0 STATE_CD 2 digit state FIPS code of the WSC maintaining the gage. Puerto Rico is listed as a state. String 2 STATE 0 0 STATE 2 char. state postal abrev. of the WSC maintaining the gage. Puerto Rico is listed as a state. String 1 SITESTATUS 0 0 SITESTATUS Active (A) or Inactive (I) where active streamgages contain streamflow data in 2003 and/or 2004 Double 8 DA_SQ_MILE 0 0 DA_SQ_MILE Reported drainage area in square miles. Stations with DA of -999999 means there is no reported DA in NWIS. Double 8 LON_SITE 0 0 LON_SITE Longitude of the streamgage (site) location - gage house in decimal degrees, NAD83 Double 8 LAT_SITE 0 0 LAT_SITE Latitude of the streamgage (site) location - gage house in decimal degrees, NAD83 Double 8 LON_NHD 0 0 LON_NHD Longitude of the NHD location in decimal degrees, NAD83 Double 8 LAT_NHD 0 0 LAT_NHD Latitude of the NHD location in decimal degrees, NAD83 String 1 REVIEWED 0 0 REVIEWED Flag to indicate review status - Y or N Integer 4 BFIYRS 0 0 BFIYRS Number of years used in the base-flow index computation Double 8 BFI_AVE 0 0 BFI_AVE Average annual base-flow index value Double 8 BFI_STDEV 0 0 BFI_STDEV Standard deviation of annual base-flow index Integer 4 GOTBFI 0 0 GOTBFI Data flag for BFI Double 8 DAY1 0 0 DAY1 First date of flow data (yyyymmdd) Double 8 DAYN 0 0 DAYN Last date of flow data (yyyymmdd) Integer 4 NDAYS 0 0 NDAYS Number of days of flow Integer 4 NDAYSGT0 0 0 NDAYSGT0 Number of days of non-zero flow Double 8 MIN_ 0 0 MIN_ Minimum daily flow for the period of record Double 8 P1 0 0 P1 1st percentile of daily flow for the period of record. Negative values indicate reverse flow; tidal or backwater Double 8 P5 0 0 P5 5th percentile of daily flow for the period of record. Negative values indicate reverse flow; tidal or backwater Double 8 P10 0 0 P10 10th percentile of daily flow for the period of record. Negative values indicate reverse flow; tidal or backwater Double 8 P20 0 0 P20 20th percentile of daily flow for the period of record. Negative values indicate reverse flow; tidal or backwater Double 8 P25 0 0 P25 25th percentile of daily flow for the period of record. Negative values indicate reverse flow; tidal or backwater Double 8 P30 0 0 P30 30th percentile of daily flow for the period of record. Negative values indicate reverse flow; tidal or backwater Double 8 P40 0 0 P40 40th percentile of daily flow for the period of record. Negative values indicate reverse flow; tidal or backwater Double 8 P50 0 0 P50 50th percentile of daily flow for the period of record. Negative values indicate reverse flow; tidal or backwater Double 8 P60 0 0 P60 60th percentile of daily flow for the period of record. Negative values indicate reverse flow; tidal or backwater Double 8 P70 0 0 P70 70th percentile of daily flow for the period of record. Negative values indicate reverse flow; tidal or backwater Double 8 P75 0 0 P75 75th percentile of daily flow for the period of record. Negative values indicate reverse flow; tidal or backwater Double 8 P80 0 0 P80 80th percentile of daily flow for the period of record. Negative values indicate reverse flow; tidal or backwater Double 8 P90 0 0 P90 90th percentile of daily flow for the period of record. Negative values indicate reverse flow; tidal or backwater Double 8 P95 0 0 P95 95th percentile of daily flow for the period of record. Negative values indicate reverse flow; tidal or backwater Double 8 P99 0 0 P99 99th percentile of daily flow for the period of record. Negative values indicate reverse flow; tidal or backwater Double 8 MAX_ 0 0 MAX_ Maximum daily flow for the period of record Double 8 AVE 0 0 AVE Average daily flow for the period of record Double 8 STDEV 0 0 STDEV Standard deviation of flow for daily the period of record Integer 4 GOTQ 0 0 Downloadable Data 3.44MB 18.6MB http://water.usgs.gov/GIS/dsdl/USGS_Streamgages-NHD_Locations_Shape.zip 20060620 REQUIRED: The organization responsible for the metadata information. REQUIRED: The person responsible for the metadata information. REQUIRED: The mailing and/or physical address for the organization or individual. REQUIRED: The city of the address. REQUIRED: The state or province of the address. REQUIRED: The ZIP or other postal code of the address. REQUIRED: The telephone number by which individuals can speak to the organization or individual. FGDC Content Standards for Digital Geospatial Metadata FGDC-STD-001-1998 local time http://www.esri.com/metadata/esriprof80.html ESRI Metadata Profile en 20060613 10340700 FALSE 20200110 13052600 20200110 13052600 {F9E503BD-32CA-48B1-B51F-BA46D18F9DEA} Streamgages_USGS 002 -12470236.164200 -11771486.103000 4267724.739400 5331162.213000 1 Projected GCS_WGS_1984 Linear Unit: Meter (1.000000) WGS_1984_Web_Mercator_Auxiliary_Sphere <ProjectedCoordinateSystem xsi:type='typens:ProjectedCoordinateSystem' xmlns:xsi='http://www.w3.org/2001/XMLSchema-instance' xmlns:xs='http://www.w3.org/2001/XMLSchema' xmlns:typens='http://www.esri.com/schemas/ArcGIS/2.4.0'><WKT>PROJCS["WGS_1984_Web_Mercator_Auxiliary_Sphere",GEOGCS["GCS_WGS_1984",DATUM["D_WGS_1984",SPHEROID["WGS_1984",6378137.0,298.257223563]],PRIMEM["Greenwich",0.0],UNIT["Degree",0.0174532925199433]],PROJECTION["Mercator_Auxiliary_Sphere"],PARAMETER["False_Easting",0.0],PARAMETER["False_Northing",0.0],PARAMETER["Central_Meridian",0.0],PARAMETER["Standard_Parallel_1",0.0],PARAMETER["Auxiliary_Sphere_Type",0.0],UNIT["Meter",1.0],AUTHORITY["EPSG",3857]]</WKT><XOrigin>-20037700</XOrigin><YOrigin>-30241100</YOrigin><XYScale>10000</XYScale><ZOrigin>-100000</ZOrigin><ZScale>10000</ZScale><MOrigin>-100000</MOrigin><MScale>10000</MScale><XYTolerance>0.001</XYTolerance><ZTolerance>0.001</ZTolerance><MTolerance>0.001</MTolerance><HighPrecision>true</HighPrecision><WKID>102100</WKID><LatestWKID>3857</LatestWKID></ProjectedCoordinateSystem> Esri ArcGIS 12.4.2.19948 Streamgages_USGS -159.750305 -65.302467 48.999947 17.956547 1 -159.750305 -65.302467 48.999947 17.956547 1 1 -112.022037 -105.745059 43.125949 35.760561 ISO 19115 Geographic Information - Metadata DIS_ESRI1.0 dataset 002 withheld Local Area Network 0.626 Shapefile File Geodatabase Feature Class GCS_North_American_1983 EPSG 8.8(9.3.1.2) 8 20200110