<?xml version='1.0' encoding='UTF-8'?>
<metadata xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance">
  <idinfo>
    <citation>
      <citeinfo>
        <origin>Noah M Schmadel</origin>
        <origin>John Gala</origin>
        <origin>Cristiana Figueroa-Kaminsky</origin>
        <origin>Zachary C. Johnson</origin>
        <origin>Daniel Wise</origin>
        <pubdate>20250131</pubdate>
        <title>Nutrient Calibration Load Estimation Data for the Puget Sound Region, 2005 through 2020</title>
        <geoform>tabular digital data</geoform>
        <pubinfo>
          <pubplace>https://www.sciencebase.gov</pubplace>
          <publish>U.S. Geological Survey</publish>
        </pubinfo>
        <othercit>Any use of trade, firm, or product names is for descriptive purposes only and does not imply endorsement by the U.S. Government.</othercit>
        <onlink>https://doi.org/10.5066/P9LY1PQF</onlink>
        <lworkcit>
          <citeinfo>
            <origin>Noah Schmadel</origin>
            <origin>Cristiana Figueroa-Kaminsky</origin>
            <origin>Daniel Wise</origin>
            <origin>Jamie K. Wasielewski</origin>
            <origin>Zachary Johnson</origin>
            <origin>Robert W. Black</origin>
            <pubdate>2026</pubdate>
            <title>Simulated seasonal loads of total nitrogen and total phosphorus by major source from watersheds draining to Washington waters of the Salish Sea, 2005 through 2020</title>
            <geoform>publication</geoform>
            <pubinfo>
              <pubplace>n/a</pubplace>
              <publish>US Geological Survey</publish>
            </pubinfo>
            <onlink>https://doi.org/10.3133/sir20265001</onlink>
          </citeinfo>
        </lworkcit>
      </citeinfo>
    </citation>
    <descript>
      <abstract>This data release contains estimated seasonal total nitrogen (TN) and total phosphorus (TP) loads during 2005 through 2020 in selected streams of the watersheds draining to Washington waters of the Salish Sea (Puget Sound region). Two load estimation methods were used: Fluxmaster (Schwarz and others, 2006, https://doi.org/10.3133/tm6B3) and Weighted Regression on Time, Discharge, and Season (WRTDS) with Kalman filtering (Hirsch and De Cicco, 2015, https://doi.org/10.1016/j.envsoft.2015.07.017). Together with the water-quality data, streamflow data, and estimated nitrogen and phosphorus loads, this release also provides copies of computer programming script that can be used to recreate the load estimates or modified for similar applications (see Process Steps). Seasonal TN and TP loads paired with streamflow were estimated at 50 locations. For convenience, seasonal summaries are also provided of observed streamflow at 128 locations and of observed water temperature at 29 locations.</abstract>
      <purpose>Nutrient load estimates were generated for use as the dependent variable in calibrating dynamic SPARROW (Spatially Referenced Regressions on Watershed attributes) models (Schwarz and others, 2006, https://doi.org/10.3133/tm6B3) of the watersheds draining to Washington waters of the Salish Sea (Puget Sound region). For convenience, as they may provide additional explanatory or validation data in SPARROW modeling efforts, seasonal summaries of observed streamflow and water temperature are included. Funding was provided by Washington State Department of Ecology and U.S. Geological Survey cooperative matching.

This data release is intended to allow users to recreate the load estimates using included input data and computer script and to support other potential uses of the estimates and data inputs. The hundreds of individual Fluxmaster and WRTDS model files included have been thoroughly reviewed for each individual stream location. However, there are still possibilities that estimated loads for any particular individual stream or period may be incorrect or unreasonable due to a combination of factors including poor model fit, poor data quality, or some unique local condition such as point source or backwater or tidal influences. Users interested in loads in a specific stream, for example, may be best served by using included data to re-calibrate WRTDS or Fluxmaster models which can then be additionally reviewed or evaluated for their intended purposes.</purpose>
      <supplinf>Copies of the Fluxmaster ("fluxmaster_puget.zip") and WRTDS ("wrtds_puget.zip") models are provided, which contain all individual discrete water-quality values, daily streamflow collected at gages, and R (https://www.r-project.org/ version 4.3.0 used) (for WRTDS) and SAS (www.sas.com) (for Fluxmaster) programming scripts. Hundreds of individual files are used and created, with files generally in comma-delimited (CSV) or tab-delimited format, but both models start at the same highest level by using the "wq_q_list.csv" as a key that matches water-quality and streamflow data. Water-quality data were combined at several locations that were originally entered into the public data repository as similar yet different site identification numbers; an additional list "wq_dup_list.csv" is provided to document which sites contained duplicated (but any duplicates were removed) or combined water-quality data. For convenience, a seasonal summary of the total nitrogen and total phosphorus loads estimated via Fluxmaster and WRTDS are provided ("TN_TP_loads_puget.csv"). There were more locations of streamflow gaging stations compared to water-quality stations. Therefore, a list of streamflow gages with mostly complete records from 2005 through 2020 ("flow_list.csv"), their upstream stream reaches identified ("key_upstream_gage.csv"), and a seasonal summary of those streamflows ("streamflow_seasonal_puget.csv") are also provided. For convenience, a seasonal summary of observed water temperature co-located with water-quality data at 29 locations is provided ("wtemp_seasonal_puget.csv").        

References Cited:
Hirsch, R. M. and L. A. De Cicco, 2015. User Guide to Exploration and Graphics for River Trends (Egret) and Dataretrieval: R Packages for Hydrologic Data.  Report 4-A10, 104 p. http://pubs.er.usgs.gov/publication/tm4A10.
Schwarz, G. E., A. B. Hoos, R. B. Alexander and R. A. Smith, 2006. Section 3. The Sparrow Surface Water-Quality Model: Theory, Application and User Documentation. USGS Techniques and Methods 6-B3, 248 p. https://doi.org/10.3133/tm6B3.</supplinf>
    </descript>
    <timeperd>
      <timeinfo>
        <rngdates>
          <begdate>20050101</begdate>
          <enddate>20201231</enddate>
        </rngdates>
      </timeinfo>
      <current>ground condition</current>
    </timeperd>
    <status>
      <progress>Complete</progress>
      <update>None planned</update>
    </status>
    <spdom>
      <bounding>
        <westbc>-124.7069</westbc>
        <eastbc>-120.6476</eastbc>
        <northbc>49.0050</northbc>
        <southbc>46.5895</southbc>
      </bounding>
    </spdom>
    <keywords>
      <theme>
        <themekt>ISO 19115 Topic Category</themekt>
        <themekey>inlandWaters</themekey>
      </theme>
      <theme>
        <themekt>None</themekt>
        <themekey>nitrogen</themekey>
        <themekey>phosphorus</themekey>
        <themekey>stream load</themekey>
        <themekey>Washington State</themekey>
      </theme>
      <theme>
        <themekt>USGS Metadata Identifier</themekt>
        <themekey>USGS:65860a21d34eff134d437203</themekey>
      </theme>
      <place>
        <placekt>U.S. Geological Survey</placekt>
        <placekey>Washington State</placekey>
        <placekey>water quality</placekey>
        <placekey>nitrogen</placekey>
        <placekey>phosphorus</placekey>
      </place>
    </keywords>
    <accconst>None.  Please see 'Distribution Info' for details.</accconst>
    <useconst>These data are in the public domain in accordance with Creative Commons Zero v1.0 Universal Public Domain Dedication (CC0-1.0) and have no use constraints.  

Users are advised to read the dataset's metadata thoroughly to understand appropriate use and data limitations. Use of this data release is considered acceptance of the limitations of this data and that the user has read and understood this metadata prior to its use in any form.</useconst>
    <ptcontac>
      <cntinfo>
        <cntperp>
          <cntper>Noah M Schmadel</cntper>
          <cntorg>U.S. Geological Survey | Earth System Processes Division</cntorg>
        </cntperp>
        <cntpos>Hydrologist</cntpos>
        <cntaddr>
          <addrtype>mailing address</addrtype>
          <address>601 SW 2nd Ave. Suite 1950</address>
          <city>Portland</city>
          <state>OR</state>
          <postal>97204</postal>
          <country>US</country>
        </cntaddr>
        <cntvoice>9714421691</cntvoice>
        <cntemail>nschmadel@usgs.gov</cntemail>
      </cntinfo>
    </ptcontac>
    <datacred>Washington State Department of Ecology</datacred>
    <native>The estimates were produced using SAS (www.sas.com) and the R software (https://www.r-project.org/ version 4.3.0 used) programming languages on a Microsoft Windows platform. R scripts were implemented using the RStudio (https://posit.co/download/rstudio-desktop/) editor. Any use of trade, firm, or product names is for descriptive purposes only and does not imply endorsement by the U.S. Government.</native>
    <crossref>
      <citeinfo>
        <origin>Noah M. Schmadel</origin>
        <origin>Cristiana Figueroa-Kaminsky</origin>
        <origin>Daniel R. Wise</origin>
        <origin>Jamie K. Wasielewski</origin>
        <origin>Zachary C. Johnson</origin>
        <origin>Robert Black</origin>
        <pubdate>2025</pubdate>
        <title>Simulated seasonal loads of total nitrogen and total phosphorus by major source from watersheds draining to Washington waters of the Salish Sea, 2005 through 2020</title>
        <geoform>publication</geoform>
        <pubinfo>
          <pubplace>Publication Warehouse</pubplace>
          <publish>U.S. Geological Survey</publish>
        </pubinfo>
      </citeinfo>
    </crossref>
    <crossref>
      <citeinfo>
        <origin>Noah M. Schmadel</origin>
        <origin>Cristiana Figueroa-Kaminsky</origin>
        <origin>Daniel R. Wise</origin>
        <origin>Jamie K. Wasielewski</origin>
        <origin>John Gala</origin>
        <origin>Zachary C. Johnson</origin>
        <pubdate>2025</pubdate>
        <title>Model application and calibration load data for seasonally dynamic total nitrogen and total phosphorus SPARROW models developed for watersheds draining to Washington waters of the Salish Sea, 2005 through 2020</title>
        <geoform>publication</geoform>
        <pubinfo>
          <pubplace>https://www.sciencebase.gov</pubplace>
          <publish>U.S. Geological Survey</publish>
        </pubinfo>
        <onlink>https://www.sciencebase.gov/catalog/item/6515efc7d34e469cabfce175</onlink>
      </citeinfo>
    </crossref>
  </idinfo>
  <dataqual>
    <attracc>
      <attraccr>Data were extensively checked and corrected for obvious errors (see Process steps); however, errors in datasets may still exist. For example, there are still possibilities that estimated loads for any particular individual stream or period may be incorrect or unreasonable due to a combination of factors including poor model fit, poor data quality, or some unique local condition such as point source or backwater or tidal influences.</attraccr>
    </attracc>
    <logic>No formal logical accuracy tests were conducted. All values are within expected range.</logic>
    <complete>Data are considered complete for the information presented, as described. Users are advised to read the rest of the metadata record for additional details.</complete>
    <posacc>
      <horizpa>
        <horizpar>No formal positional accuracy tests were conducted</horizpar>
      </horizpa>
      <vertacc>
        <vertaccr>No formal positional accuracy tests were conducted</vertaccr>
      </vertacc>
    </posacc>
    <lineage>
      <srcinfo>
        <srccite>
          <citeinfo>
            <origin>Water Quality Portal</origin>
            <pubdate>2022</pubdate>
            <title>Water-quality data downloaded from the Water-Quality Portal</title>
            <geoform>tabular digital data</geoform>
            <onlink>www.waterqualitydata.us</onlink>
          </citeinfo>
        </srccite>
        <typesrc>Digital</typesrc>
        <srctime>
          <timeinfo>
            <rngdates>
              <begdate>20050101</begdate>
              <enddate>20201231</enddate>
            </rngdates>
          </timeinfo>
          <srccurr>Data are current as of 7/1/2022</srccurr>
        </srctime>
        <srccitea>WQP</srccitea>
        <srccontr>Contributed to water-quality data</srccontr>
      </srcinfo>
      <srcinfo>
        <srccite>
          <citeinfo>
            <origin>David A. Saad</origin>
            <origin>Gregory E. Schwarz</origin>
            <origin>Denise M. Argue</origin>
            <origin>David W. Anning</origin>
            <origin>Scott W. Ator</origin>
            <origin>Anne B. Hoos</origin>
            <origin>Stephen D. Preston</origin>
            <origin>Dale M. Robertson</origin>
            <origin>Daniel Wise</origin>
            <pubdate>20190819</pubdate>
            <title>Estimates of long-term mean daily streamflow and annual nutrient and suspended-sediment loads considered for use in regional SPARROW models of the conterminous United States, 2012 base year</title>
            <geoform>publication</geoform>
            <onlink>doi.org/10.5066/F7DN436B</onlink>
          </citeinfo>
        </srccite>
        <typesrc>Digital</typesrc>
        <srctime>
          <timeinfo>
            <rngdates>
              <begdate>19991001</begdate>
              <enddate>20140930</enddate>
            </rngdates>
          </timeinfo>
          <srccurr>Data are current as of the original publication date.</srccurr>
        </srctime>
        <srccitea>Saad and others (2019)</srccitea>
        <srccontr>Contributed to water-quality data included in the folder, "saad2019".</srccontr>
      </srcinfo>
      <srcinfo>
        <srccite>
          <citeinfo>
            <origin>Washington State Department of Ecology</origin>
            <pubdate>2019</pubdate>
            <title>EIM - environmental monitoring data. Environmental Assessment Program, Washington State Department of Ecology, Olympia, WA</title>
            <geoform>application/service</geoform>
            <onlink>https://ecology.wa.gov/Research-Data/Data-resources/Environmental-Information-Management-database</onlink>
          </citeinfo>
        </srccite>
        <typesrc>Digital and/or Hardcopy</typesrc>
        <srctime>
          <timeinfo>
            <rngdates>
              <begdate>20050101</begdate>
              <enddate>20201231</enddate>
            </rngdates>
          </timeinfo>
          <srccurr>observed</srccurr>
        </srctime>
        <srccitea>EC</srccitea>
        <srccontr>Contributed water-quality and streamflow data</srccontr>
      </srcinfo>
      <srcinfo>
        <srccite>
          <citeinfo>
            <origin>King County</origin>
            <pubdate>2024</pubdate>
            <title>Hydrologic Information Center Data Download, King County Water and Land Resource Division, WA</title>
            <geoform>application/service</geoform>
            <onlink>https://green2.kingcounty.gov/hydrology/Data.aspx</onlink>
          </citeinfo>
        </srccite>
        <typesrc>Digital and/or Hardcopy</typesrc>
        <srctime>
          <timeinfo>
            <rngdates>
              <begdate>20050101</begdate>
              <enddate>20201231</enddate>
            </rngdates>
          </timeinfo>
          <srccurr>observed</srccurr>
        </srctime>
        <srccitea>KC</srccitea>
        <srccontr>Contributed water-quality and streamflow data</srccontr>
      </srcinfo>
      <srcinfo>
        <srccite>
          <citeinfo>
            <origin>U.S. Geological Survey</origin>
            <pubdate>1994</pubdate>
            <title>USGS Water Data for the Nation</title>
            <geoform>application/service</geoform>
            <pubinfo>
              <pubplace>http://waterdata.usgs.gov</pubplace>
              <publish>U.S. Geological Survey</publish>
            </pubinfo>
            <onlink>https://doi.org/10.5066/f7p55kjn</onlink>
          </citeinfo>
        </srccite>
        <typesrc>Digital and/or Hardcopy</typesrc>
        <srctime>
          <timeinfo>
            <rngdates>
              <begdate>20050101</begdate>
              <enddate>20201231</enddate>
            </rngdates>
          </timeinfo>
          <srccurr>observed</srccurr>
        </srctime>
        <srccitea>NWIS</srccitea>
        <srccontr>Provided streamflow data</srccontr>
      </srcinfo>
      <srcinfo>
        <srccite>
          <citeinfo>
            <origin>Schwarz, G.E.</origin>
            <pubdate>2019</pubdate>
            <title>E2NHDPlusV2_us: Database of Ancillary Hydrologic Attributes and Modified Routing for NHDPlus Version 2.1 Flowlines</title>
            <geoform>dataset</geoform>
            <pubinfo>
              <pubplace>https://www.sciencebase.gov</pubplace>
              <publish>U.S. Geological Survey</publish>
            </pubinfo>
            <onlink>https://doi.org/10.5066/p986kzem</onlink>
          </citeinfo>
        </srccite>
        <typesrc>Digital and/or Hardcopy</typesrc>
        <srctime>
          <timeinfo>
            <sngdate>
              <caldate>20190101</caldate>
            </sngdate>
          </timeinfo>
          <srccurr>publication date</srccurr>
        </srctime>
        <srccitea>e2NHD</srccitea>
        <srccontr>Provided digital stream network navigation attributes</srccontr>
      </srcinfo>
      <srcinfo>
        <srccite>
          <citeinfo>
            <origin>U.S. Environmental Protection Agency</origin>
            <pubdate>20190313</pubdate>
            <title>NHDPlus National Data</title>
            <geoform>tabular digital data</geoform>
            <pubinfo>
              <pubplace>Online</pubplace>
              <publish>Water Data and Tools</publish>
            </pubinfo>
            <onlink>https://www.epa.gov/waterdata/nhdplus-national-data</onlink>
          </citeinfo>
        </srccite>
        <typesrc>Digital and/or Hardcopy</typesrc>
        <srctime>
          <timeinfo>
            <sngdate>
              <caldate>20190313</caldate>
            </sngdate>
          </timeinfo>
          <srccurr>publication date</srccurr>
        </srctime>
        <srccitea>NHDPlusV2</srccitea>
        <srccontr>Provided flowline, waterbody, and catchment shapefiles</srccontr>
      </srcinfo>
      <srcinfo>
        <srccite>
          <citeinfo>
            <origin>Robert M. Hirsch</origin>
            <origin>Stacey A. Archfield</origin>
            <origin>Laura A. De Cicco</origin>
            <pubdate>201511</pubdate>
            <title>A bootstrap method for estimating uncertainty of water quality trends</title>
            <geoform>publication</geoform>
            <serinfo>
              <sername>Environmental Modelling &amp;amp; Software</sername>
              <issue>vol. 73</issue>
            </serinfo>
            <pubinfo>
              <pubplace>n/a</pubplace>
              <publish>Elsevier BV</publish>
            </pubinfo>
            <othercit>ppg. 148-166</othercit>
            <onlink>https://doi.org/10.1016/j.envsoft.2015.07.017</onlink>
          </citeinfo>
        </srccite>
        <typesrc>Digital and/or Hardcopy</typesrc>
        <srctime>
          <timeinfo>
            <sngdate>
              <caldate>201511</caldate>
            </sngdate>
          </timeinfo>
          <srccurr>publication date</srccurr>
        </srctime>
        <srccitea>WRTDS</srccitea>
        <srccontr>Citation of WRTDS model used to estimate loads</srccontr>
      </srcinfo>
      <srcinfo>
        <srccite>
          <citeinfo>
            <origin>Gregory Schwarz</origin>
            <origin>Anne B. Hoos</origin>
            <origin>R.B. Alexander</origin>
            <origin>R.A. Smith</origin>
            <pubdate>2006</pubdate>
            <title>Section 3. The SPARROW Surface Water-Quality Model—Theory, application and user documentation</title>
            <geoform>publication</geoform>
            <pubinfo>
              <pubplace>n/a</pubplace>
              <publish>US Geological Survey</publish>
            </pubinfo>
            <onlink>https://doi.org/10.3133/tm6B3</onlink>
          </citeinfo>
        </srccite>
        <typesrc>Digital and/or Hardcopy</typesrc>
        <srctime>
          <timeinfo>
            <sngdate>
              <caldate>2006</caldate>
            </sngdate>
          </timeinfo>
          <srccurr>publication date</srccurr>
        </srctime>
        <srccitea>Fluxmaster</srccitea>
        <srccontr>Citation for Fluxmaster model used to estimate loads</srccontr>
      </srcinfo>
      <srcinfo>
        <srccite>
          <citeinfo>
            <origin>National Oceanic and Atmospheric Administration</origin>
            <pubdate>20220101</pubdate>
            <title>National Water Model 2022 CONUS Retrospective Dataset</title>
            <geoform>tabular digital data</geoform>
            <pubinfo>
              <pubplace>Online</pubplace>
              <publish>National Oceanic and Atmospheric Administration</publish>
            </pubinfo>
            <othercit>accessed February 23, 2022</othercit>
            <onlink>https://registry.opendata.aws/nwm-archive/</onlink>
          </citeinfo>
        </srccite>
        <typesrc>Digital and/or Hardcopy</typesrc>
        <srctime>
          <timeinfo>
            <rngdates>
              <begdate>20050101</begdate>
              <enddate>20201231</enddate>
            </rngdates>
          </timeinfo>
          <srccurr>ground condition</srccurr>
        </srctime>
        <srccitea>WRF</srccitea>
        <srccontr>Provided simulated daily streamflow at a few select stream reaches</srccontr>
      </srcinfo>
      <srcinfo>
        <srccite>
          <citeinfo>
            <origin>U.S. Geological Survey</origin>
            <pubdate>20230901</pubdate>
            <title>Discrete water-quality data compiled for use in load estimation</title>
            <geoform>tabular digital data</geoform>
          </citeinfo>
        </srccite>
        <typesrc>Digital and/or Hardcopy</typesrc>
        <srctime>
          <timeinfo>
            <rngdates>
              <begdate>20050101</begdate>
              <enddate>20201231</enddate>
            </rngdates>
          </timeinfo>
          <srccurr>observed</srccurr>
        </srctime>
        <srccitea>WQ Data</srccitea>
        <srccontr>Water-quality data in this data release, product or Process Step 1</srccontr>
      </srcinfo>
      <srcinfo>
        <srccite>
          <citeinfo>
            <origin>U.S. Geological Survey</origin>
            <pubdate>20231010</pubdate>
            <title>Streamflow data compiled for use in load estimation and streamflow adjustment</title>
            <geoform>tabular digital data</geoform>
          </citeinfo>
        </srccite>
        <typesrc>Digital and/or Hardcopy</typesrc>
        <srctime>
          <timeinfo>
            <rngdates>
              <begdate>20050101</begdate>
              <enddate>20201231</enddate>
            </rngdates>
          </timeinfo>
          <srccurr>observed</srccurr>
        </srctime>
        <srccitea>Q Data</srccitea>
        <srccontr>Streamflow data in this data release, product of Process Step 2</srccontr>
      </srcinfo>
      <srcinfo>
        <srccite>
          <citeinfo>
            <origin>U.S. Geological Survey</origin>
            <pubdate>20240110</pubdate>
            <title>Matched pairs of water-quality and streamflow locations</title>
            <geoform>tabular digital data</geoform>
          </citeinfo>
        </srccite>
        <typesrc>Digital and/or Hardcopy</typesrc>
        <srctime>
          <timeinfo>
            <sngdate>
              <caldate>20241231</caldate>
            </sngdate>
          </timeinfo>
          <srccurr>publication date</srccurr>
        </srctime>
        <srccitea>wq_q_list.csv</srccitea>
        <srccontr>Spatial linking of water-quality and streamflow locations in this data release, product of Process Step 3</srccontr>
      </srcinfo>
      <srcinfo>
        <srccite>
          <citeinfo>
            <origin>U.S. Geological Survey</origin>
            <pubdate>20240110</pubdate>
            <title>Fluxmaster and WRTDS seasonal TN and TP load summary</title>
            <geoform>tabular digital data</geoform>
          </citeinfo>
        </srccite>
        <typesrc>Digital and/or Hardcopy</typesrc>
        <srctime>
          <timeinfo>
            <rngdates>
              <begdate>20050101</begdate>
              <enddate>20201231</enddate>
            </rngdates>
          </timeinfo>
          <srccurr>ground condition</srccurr>
        </srctime>
        <srccitea>TN_TP_loads_puget.csv</srccitea>
        <srccontr>TN and TP seasonal loads in this data release, product of Process Step 4</srccontr>
      </srcinfo>
      <srcinfo>
        <srccite>
          <citeinfo>
            <origin>U.S. Department of Agriculture</origin>
            <pubdate>2024</pubdate>
            <title>NorWEST Stream Temp,</title>
            <geoform>tabular digital data</geoform>
            <othercit>accessed December 01, 2023</othercit>
            <onlink>https://www.fs.usda.gov/rm/boise/AWAE/projects/NorWeST.html</onlink>
          </citeinfo>
        </srccite>
        <typesrc>Digital and/or Hardcopy</typesrc>
        <srctime>
          <timeinfo>
            <sngdate>
              <caldate>2024</caldate>
            </sngdate>
          </timeinfo>
          <srccurr>publication date</srccurr>
        </srctime>
        <srccitea>norwest</srccitea>
        <srccontr>Observed water temperature</srccontr>
      </srcinfo>
      <srcinfo>
        <srccite>
          <citeinfo>
            <origin>Scott W Ator</origin>
            <origin>David A Saad</origin>
            <origin>Lauren E Koenig Snyder</origin>
            <origin>Zachary C Johnson</origin>
            <origin>Gretchen P Oelsner</origin>
            <pubdate>2024</pubdate>
            <title>Estimated seasonal nitrogen and phosphorus loads in selected streams of the conterminous United States, 1999 - 2020</title>
            <geoform>dataset</geoform>
            <pubinfo>
              <pubplace>https://www.sciencebase.gov</pubplace>
              <publish>U.S. Geological Survey</publish>
            </pubinfo>
            <onlink>https://doi.org/10.5066/p13ysthm</onlink>
          </citeinfo>
        </srccite>
        <typesrc>Digital and/or Hardcopy</typesrc>
        <srctime>
          <timeinfo>
            <sngdate>
              <caldate>2024</caldate>
            </sngdate>
          </timeinfo>
          <srccurr>publication date</srccurr>
        </srctime>
        <srccitea>Ator</srccitea>
        <srccontr>Fluxmaster options</srccontr>
      </srcinfo>
      <procstep>
        <procdesc>1. Retrieval, Filtering, and Harmonization of Water-Quality Data 

Available water-quality data (sampled instream nitrogen and phosphorus concentrations) within the Puget Sound region were retrieved by Washington State Department of Ecology staff in June 2022. Data were retrieved from numerous sources including EC, KC, WQP, and NWIS. The EC datasouce is a comprehensive data repository that includes Tribal and county data, and also uploads to the WQP. The WQP was checked for any additional data not already included in EC, KC, or NWIS.     

Upon initial data retrieval, several of the water-quality data station names were different yet similar (such as names KCM-0438 and 0438) for data collected at the same location (or very near same location); therefore, any water-quality data collected at the same location were grouped and renamed to the most frequently used identification number. Regular water surface samples were retrieved to exclude quality control data; however, grouping by similar name caused some duplicate records and a few cases where one date had more than one value. For those cases, the other value was near zero, implying a quality control record. Therefore, keeping maximum by date was performed to retain only field samples and also remove any duplicates. 

All data were screened and standardized following methods adapted from Saad and others (2019). For example, data were screened to limit observations to those of total (unfiltered) forms of nitrogen or phosphorus under ambient conditions in non-tidal streams and to remove observations with insufficient metadata to reasonably verify such conditions. Concentrations of total nitrogen were computed where possible from reported concentrations of individual nitrogen species; however, nearly all nitrogen data used were originally reported as unfiltered total form.</procdesc>
        <srcused>WQP</srcused>
        <srcused>Saad and others (2019)</srcused>
        <srcused>EC</srcused>
        <srcused>KC</srcused>
        <srcused>NWIS</srcused>
        <procdate>20230901</procdate>
        <srcprod>WQ Data</srcprod>
      </procstep>
      <procstep>
        <procdesc>2. Retrieval of Streamflow Data

Daily streamflow data used to estimate loads were retrieved June 2022 from EC and KC and September 2023 from NWIS for 2005 to 2020. At a few select stream reaches (n=5) that contained water-quality data but not a streamflow gage, daily streamflow estimates from WRF was used as a streamflow proxy.  

For convenience, seasonal and mean summaries of observed streamflow are provided in "streamflow_seasonal_puget.csv" where at least 85 streamflow values had to exist per season to keep a seasonal estimate.</procdesc>
        <srcused>NWIS</srcused>
        <srcused>EC</srcused>
        <srcused>KC</srcused>
        <srcused>WRF</srcused>
        <procdate>20231010</procdate>
        <srcprod>Q Data</srcprod>
      </procstep>
      <procstep>
        <procdesc>3. Continuous streamflow estimated at each water-quality station

The nearest NHDPlusV2 (same for e2NHD) common identifier (COMID) was assigned to each water-quality and streamflow station based on the latitudes and longitudes or each station. Every COMID-station assignment was visually checked for correctness by viewing in maps and inspecting the closest NHDPlusV2 flowline reach to the station point.

Many water-quality and streamflow stations are co-located, but some are not. To determine if a nearby streamflow station could be used as a representative streamflow at a nearby water-quality station, ratios of drainage areas (provided by e2NHD) were checked. Streamflow was considered representative at a water-quality station if the ratio of drainage area between the water-quality and streamflow station was greater than 0.75 (water-quality station upstream) or less than 1.33 (water-quality station downstream), or no more than a 25% difference in areas. 

Locations relative to other features such as diversions, reservoirs, point sources, and stormwater outfalls were also checked by viewing maps. Ideally, the collection location of water-quality and streamflow would be the same or with nothing substantial in between, but pairing to locations near by was done to include more calibration stations. For example, instead of combining closely located water-quality stations 09A080 and KCM-3106 on the Green River that were originally assigned to the same COMID reach, station 09A080 was re-assigned to the next upstream COMID because several stormwater outfalls are located in-between those stations.      

The final station list is provided in "wq_q_list.csv" where 50 water-quality station locations were identified that paired with 47 streamflow gages. There were 127 streamflow gage locations identified where 80 of those gages were not paired with water-quality stations. For convenience, seasonal streamflow was summarized at those additional gages and provided here, together with the COMID of every upstream stream reach that flows into each streamflow gage in "key_upstream_gage.csv" and a list of all gages in "flow_list.csv," to allow for possible estimation of streamflow upstream and downstream of gages at more reach locations.</procdesc>
        <srcused>NHDPlusV2</srcused>
        <srcused>e2NHD</srcused>
        <procdate>20240110</procdate>
        <srcprod>wq_q_list.csv</srcprod>
      </procstep>
      <procstep>
        <procdesc>4. Load Estimation 

Total nitrogen and phosphorus loads were estimated using Fluxmaster and WRTDS similar to Ator, and compiled seasonally into "TN_TP_loads_puget.csv" for convenience.  

Fluxmaster runs in the SAS computing environment (https://www.sas.com/en_us/home.html, version SAS 9.4 TS Level 1M7, X64_10PRO platform). A 5-parameter (intercept, streamflow, season, trend, serial correlation) regression method was selected and requires three input files: 1) discretely measured constituent concentrations at a station, 2) daily observed (or simulated) streamflow at a station, and 3) a file describing the the spatial linking of water-quality and streamflow locations including drainage area ratios. WRTDS runs within the R software (https://www.r-project.org/ version 4.3.0 used) programming environment using the RStudio (https://posit.co/download/rstudio-desktop/) editor and was completed using Kalman smoothing. WRTDS also requires 1) constiuent data, 2) daily streamflow, and 3) station information. Streamflow was scaled by the drainage area ratio to better match water-quality location (note that many ratios were 1). Both modeling tools (Fluxmaster and WRTDS) were executed using computer code included in this data release (see below). 

Fluxmaster results can be reproduced as follows:

1. Read "control" files into the SAS Program Editor:
a.fluxmaster_puget/control/fluxmaster_control_Puget_TN.sas (for TN load estimates)
b.fluxmaster_puget/control/fluxmaster_control_Puget_TP.sas (for TP load estimates)

2. In the control file, update the directory "home_dir" to the to the user's local computer pathname where fluxmaster_puget was saved. The control file will automatically locate the provided Fluxmaster scripts (fluxmaster_puget/master_100418).  

3. Import to SAS datasets:
a: fluxmaster_puget/data/sitegage.txt 
b: fluxmaster_puget/data/wq_data_TN.txt (for TN load estimates)
c. fluxmaster_puget/data/wq_data_TP.txt (for TP load estimates)
(note: a SAS script "make_data.sas" is provided to convert text files to SAS-readable, but copies in SAS-ready format are provided for convenience)

4. Run the .sas code (control files) in the Program Editor by clicking on the "running person" icon. An warning message in the Console may occur, which is expected and innocuous. The control file will use code within the zipped directory: master_100418.zip. The following output files are produced:
a. fluxmaster_puget/results/tn/predict_flux_all_sub_periods_TN.txt (for TN load estimates)
b. fluxmaster_puget/results/tn/wq_summary_betaest_all_TN.txt (for TN load estimates)
c. fluxmaster_puget/results/tp/predict_flux_all_sub_periods_TP.txt (for TP load estimates)
d. fluxmaster_puget/results/tp/wq_summary_betaest_all_TP.txt (for TP load estimates)

WRTDS results can be reproduced as follows:

1. To run WRTDS and recreate results in wrtds_puget/eLists: 
a. read wrtds_puget/targets1.R (same script for TN and TP) into RStudio. The "working directory" (wd) should automatically set to the user's local computer pathname where wrtds_puget was saved.
b. Run targets1.R by clicking "Source" in RStudio. To generate TN load estimates, on line 23 to 27, set "varmn = "TN" " and "path_wq = "/wq_files/tn" ". To generate TP load estimates, set "varmn = "TP" " and "path_wq = "/wq_files/tp" " and re-run. Any notes or warnings in the Console are expected and innocuous. Model output is stored as daily "eLists" and file names are automatically saved in \wrtds_puget\eLists\TN or TP by station name. 

2. To aggregate daily eLists to seasonal loads
a. read wrtds_puget/targets2.R (same script for TN and TP) into RStudio. The "working directory" (wd) should automatically set to the user's local computer pathname where wrtds_puget was saved.
b. To generate seasonal TN load estimates, on line 22 and 23, set "varmn = "TN" ". To generate seasonal TP load estimates, set "varmn = "TP" " and re-run. Seasonal load results are automatically stored in \wrtds_puget\seasonal_fluxes\TN or TP.</procdesc>
        <srcused>Fluxmaster</srcused>
        <srcused>WRTDS</srcused>
        <srcused>WQ Data</srcused>
        <srcused>Q Data</srcused>
        <srcused>wq_q_list.csv</srcused>
        <srcused>Ator</srcused>
        <procdate>20240110</procdate>
        <srcprod>TN_TP_loads_puget.csv</srcprod>
      </procstep>
    </lineage>
  </dataqual>
  <eainfo>
    <detailed>
      <enttyp>
        <enttypl>wq_q_list.csv</enttypl>
        <enttypd>Comma Separated Value (CSV) file containing a list key of water-quality stations paired with streamflow gages.</enttypd>
        <enttypds>Producer Defined</enttypds>
      </enttyp>
      <attr>
        <attrlabl>fluxmaster</attrlabl>
        <attrdef>Flag Y = Fluxmaster TN load used instead of WRTDS TN load, N otherwise</attrdef>
        <attrdefs>Producer Defined</attrdefs>
        <attrdomv>
          <edom>
            <edomv>N</edomv>
            <edomvd>WRTDS TN load used instead of Fluxmaster TN load</edomvd>
            <edomvds>Producer defined</edomvds>
          </edom>
        </attrdomv>
        <attrdomv>
          <edom>
            <edomv>Y</edomv>
            <edomvd>Fluxmaster TN load used</edomvd>
            <edomvds>Producer defined</edomvds>
          </edom>
        </attrdomv>
      </attr>
      <attr>
        <attrlabl>target_site_previous</attrlabl>
        <attrdef>Previous identification number temporarily assigned to water-quality stations with different name but similar location, see wq_dup_list.csv to lookup which data were combined with "wq_dup_" followed by this number.</attrdef>
        <attrdefs>Producer Defined</attrdefs>
        <attrdomv>
          <udom>Unique character identifiers for each site</udom>
        </attrdomv>
      </attr>
      <attr>
        <attrlabl>target_site</attrlabl>
        <attrdef>Water-quality station identification number, which sometimes have lead or trailing zeros, from various data sources (see Process steps)</attrdef>
        <attrdefs>Producer Defined</attrdefs>
        <attrdomv>
          <udom>Unique character identifiers for each site</udom>
        </attrdomv>
      </attr>
      <attr>
        <attrlabl>target_nm</attrlabl>
        <attrdef>Water-quality station name from various data sources (see Process steps)</attrdef>
        <attrdefs>Producer Defined</attrdefs>
        <attrdomv>
          <udom>Unique character name for each site</udom>
        </attrdomv>
      </attr>
      <attr>
        <attrlabl>target_comid</attrlabl>
        <attrdef>Stream reach common identifier (COMID) from e2NHD associated with water-quality station</attrdef>
        <attrdefs>Producer Defined</attrdefs>
        <attrdomv>
          <rdom>
            <rdommin>23955772</rdommin>
            <rdommax>947110181</rdommax>
          </rdom>
        </attrdomv>
      </attr>
      <attr>
        <attrlabl>target_lat</attrlabl>
        <attrdef>Water-quality station latitude in decimal degrees from various data sources (see Process steps)</attrdef>
        <attrdefs>Producer Defined</attrdefs>
        <attrdomv>
          <rdom>
            <rdommin>47.01176</rdommin>
            <rdommax>48.9241</rdommax>
            <attrunit>Decimal degrees</attrunit>
          </rdom>
        </attrdomv>
      </attr>
      <attr>
        <attrlabl>target_lon</attrlabl>
        <attrdef>Water-quality station longitude in decimal degrees from from various data sources (see Process steps)</attrdef>
        <attrdefs>Producer Defined</attrdefs>
        <attrdomv>
          <rdom>
            <rdommin>-123.9540923</rdommin>
            <rdommax>-121.42902</rdommax>
            <attrunit>Decimal degrees</attrunit>
          </rdom>
        </attrdomv>
      </attr>
      <attr>
        <attrlabl>wtemp_source_target</attrlabl>
        <attrdef>Data source of water temperature observed at water-quality station, EC or norwest</attrdef>
        <attrdefs>Producer Defined</attrdefs>
        <attrdomv>
          <edom>
            <edomv>blank</edomv>
            <edomvd>No data</edomvd>
            <edomvds>Producer defined</edomvds>
          </edom>
        </attrdomv>
        <attrdomv>
          <edom>
            <edomv>EC</edomv>
            <edomvd>Washington State Department of Ecology observed water temperature</edomvd>
            <edomvds>Producer defined</edomvds>
          </edom>
        </attrdomv>
        <attrdomv>
          <edom>
            <edomv>norwest</edomv>
            <edomvd>U.S. Department of Agriculture NorWEST Stream Temperature</edomvd>
            <edomvds>Producer defined</edomvds>
          </edom>
        </attrdomv>
      </attr>
      <attr>
        <attrlabl>q_lat</attrlabl>
        <attrdef>Streamflow gage latitude in decimal degrees from NWIS, EC, or KC data sources</attrdef>
        <attrdefs>Producer Defined</attrdefs>
        <attrdomv>
          <rdom>
            <rdommin>46.93343245</rdommin>
            <rdommax>48.93194444</rdommax>
            <attrunit>Decimal degrees</attrunit>
          </rdom>
        </attrdomv>
      </attr>
      <attr>
        <attrlabl>q_lon</attrlabl>
        <attrdef>Streamflow gage longitude in decimal degrees from NWIS, EC, or KC data sources</attrdef>
        <attrdefs>Producer Defined</attrdefs>
        <attrdomv>
          <rdom>
            <rdommin>-123.9658333</rdommin>
            <rdommax>-121.4298499</rdommax>
            <attrunit>Decimal degrees</attrunit>
          </rdom>
        </attrdomv>
      </attr>
      <attr>
        <attrlabl>q_site</attrlabl>
        <attrdef>Streamflow gage identification number, which sometimes have leading or trailing zeros, from NWIS, EC, or KC data sources</attrdef>
        <attrdefs>Producer Defined</attrdefs>
        <attrdomv>
          <udom>Unique character identifiers for each site</udom>
        </attrdomv>
      </attr>
      <attr>
        <attrlabl>q_site_nm</attrlabl>
        <attrdef>Streamflow gage name from NWIS, EC, or KC data sources</attrdef>
        <attrdefs>Producer Defined</attrdefs>
        <attrdomv>
          <udom>Unique character name for each site</udom>
        </attrdomv>
      </attr>
      <attr>
        <attrlabl>q_comid</attrlabl>
        <attrdef>Stream reach common identifier (COMID) from e2NHD associated with streamflow gage</attrdef>
        <attrdefs>Producer Defined</attrdefs>
        <attrdomv>
          <rdom>
            <rdommin>23955772</rdommin>
            <rdommax>24538380</rdommax>
          </rdom>
        </attrdomv>
      </attr>
      <attr>
        <attrlabl>q_CumAreaKm2</attrlabl>
        <attrdef>Accumulated drainage area from e2NHD source data</attrdef>
        <attrdefs>Producer Defined</attrdefs>
        <attrdomv>
          <rdom>
            <rdommin>3.996</rdommin>
            <rdommax>8026.9668</rdommax>
            <attrunit>Square kilometers</attrunit>
          </rdom>
        </attrdomv>
      </attr>
      <attr>
        <attrlabl>drain_ratio</attrlabl>
        <attrdef>Drainage area ratio of the water-quality station to the streamflow gage where greater than 1 indicates station is downstream of gage</attrdef>
        <attrdefs>Producer Defined</attrdefs>
        <attrdomv>
          <rdom>
            <rdommin>0.780680745</rdommin>
            <rdommax>1.317580248</rdommax>
            <attrunit>Dimensionless</attrunit>
          </rdom>
        </attrdomv>
      </attr>
      <attr>
        <attrlabl>dist_km</attrlabl>
        <attrdef>Distance along e2NHD stream reaches between water-quality station and streamflow gage</attrdef>
        <attrdefs>Producer Defined</attrdefs>
        <attrdomv>
          <rdom>
            <rdommin>0.0</rdommin>
            <rdommax>34.162</rdommax>
            <attrunit>Kilometers</attrunit>
          </rdom>
        </attrdomv>
      </attr>
      <attr>
        <attrlabl>q_source</attrlabl>
        <attrdef>Data source of streamflow, NWIS, EC, KC, or WRF</attrdef>
        <attrdefs>Producer Defined</attrdefs>
        <attrdomv>
          <edom>
            <edomv>NWIS</edomv>
            <edomvd>U.S. Geological survey observed streamflow</edomvd>
            <edomvds>Producer defined</edomvds>
          </edom>
        </attrdomv>
        <attrdomv>
          <edom>
            <edomv>EC</edomv>
            <edomvd>Washington State Department of Ecology observed streamflow</edomvd>
            <edomvds>Producer defined</edomvds>
          </edom>
        </attrdomv>
        <attrdomv>
          <edom>
            <edomv>WRF</edomv>
            <edomvd>Simulated streamflow by the National Water Model (WRF-Hydro)</edomvd>
            <edomvds>Producer defined</edomvds>
          </edom>
        </attrdomv>
        <attrdomv>
          <edom>
            <edomv>KC</edomv>
            <edomvd>King County observed streamflow</edomvd>
            <edomvds>Producer defined</edomvds>
          </edom>
        </attrdomv>
      </attr>
      <attr>
        <attrlabl>wtemp_source_q</attrlabl>
        <attrdef>Data source of water temperature observed at streamflow gage, NWIS, EC, norwest, or WQP</attrdef>
        <attrdefs>Producer Defined</attrdefs>
        <attrdomv>
          <edom>
            <edomv>blank</edomv>
            <edomvd>No data</edomvd>
            <edomvds>Producer defined</edomvds>
          </edom>
        </attrdomv>
        <attrdomv>
          <edom>
            <edomv>NWIS</edomv>
            <edomvd>U.S. Geological survey observed water temperature</edomvd>
            <edomvds>Producer defined</edomvds>
          </edom>
        </attrdomv>
        <attrdomv>
          <edom>
            <edomv>EC</edomv>
            <edomvd>Washington State Department of Ecology observed water temperature</edomvd>
            <edomvds>Producer defined</edomvds>
          </edom>
        </attrdomv>
        <attrdomv>
          <edom>
            <edomv>norwest</edomv>
            <edomvd>U.S. Department of Agriculture NorWEST Stream Temperature</edomvd>
            <edomvds>Producer defined</edomvds>
          </edom>
        </attrdomv>
        <attrdomv>
          <edom>
            <edomv>WQP</edomv>
            <edomvd>Water Quality Portal</edomvd>
            <edomvds>Producer defined</edomvds>
          </edom>
        </attrdomv>
      </attr>
      <attr>
        <attrlabl>LevelPathI</attrlabl>
        <attrdef>Level path (Hydroseq number identifier for most downstream stream reach COMID in path) from e2NHD</attrdef>
        <attrdefs>USGS Defined</attrdefs>
        <attrdomv>
          <rdom>
            <rdommin>53879</rdommin>
            <rdommax>796070</rdommax>
          </rdom>
        </attrdomv>
      </attr>
    </detailed>
    <detailed>
      <enttyp>
        <enttypl>wq_dup_list.csv</enttypl>
        <enttypd>Comma Separated Value (CSV) file containing list of water-quality stations with combined data.</enttypd>
        <enttypds>Producer Defined</enttypds>
      </enttyp>
      <attr>
        <attrlabl>site</attrlabl>
        <attrdef>Raw water-quality station site identification character as entered into the public data repository</attrdef>
        <attrdefs>Producer Defined</attrdefs>
        <attrdomv>
          <udom>Unique character identifiers for each site</udom>
        </attrdomv>
      </attr>
      <attr>
        <attrlabl>dup_no</attrlabl>
        <attrdef>Assigned number that identifies stations with different identification but same location, where duplicated numbers indicate same location</attrdef>
        <attrdefs>Producer Defined</attrdefs>
        <attrdomv>
          <rdom>
            <rdommin>6</rdommin>
            <rdommax>140</rdommax>
          </rdom>
        </attrdomv>
      </attr>
      <attr>
        <attrlabl>comid</attrlabl>
        <attrdef>Stream reach common identifier (COMID) from e2NHD associated with water-quality site</attrdef>
        <attrdefs>USGS Defined</attrdefs>
        <attrdomv>
          <rdom>
            <rdommin>23956550</rdommin>
            <rdommax>947110180</rdommax>
          </rdom>
        </attrdomv>
      </attr>
    </detailed>
    <detailed>
      <enttyp>
        <enttypl>TN_TP_loads_puget.csv</enttypl>
        <enttypd>Comma Separated Value (CSV) file containing a seasonal summary of TN and TP loads.</enttypd>
        <enttypds>Producer Defined</enttypds>
      </enttyp>
      <attr>
        <attrlabl>year</attrlabl>
        <attrdef>Calendar year</attrdef>
        <attrdefs>Producer Defined</attrdefs>
        <attrdomv>
          <rdom>
            <rdommin>2005</rdommin>
            <rdommax>2020</rdommax>
          </rdom>
        </attrdomv>
      </attr>
      <attr>
        <attrlabl>quarter</attrlabl>
        <attrdef>Quarter year or season, 1 = winter Jan Feb Mar, 2 = spring Apr May Jun, 3 = summer Jul Aug Sep, 4 = fall Oct Nov Dec</attrdef>
        <attrdefs>Producer Defined</attrdefs>
        <attrdomv>
          <rdom>
            <rdommin>1</rdommin>
            <rdommax>4</rdommax>
            <attrunit>Quarter year</attrunit>
          </rdom>
        </attrdomv>
      </attr>
      <attr>
        <attrlabl>comid</attrlabl>
        <attrdef>Stream reach common identifier (COMID) from e2NHD associated with water-quality station</attrdef>
        <attrdefs>Producer Defined</attrdefs>
        <attrdomv>
          <rdom>
            <rdommin>23955772</rdommin>
            <rdommax>947110181</rdommax>
          </rdom>
        </attrdomv>
      </attr>
      <attr>
        <attrlabl>TN_kfluxkg</attrlabl>
        <attrdef>Seasonal TN load estimated using WRTDS with Kalman filtering</attrdef>
        <attrdefs>Producer Defined</attrdefs>
        <attrdomv>
          <edom>
            <edomv>&lt;&lt; empty cell &gt;&gt;</edomv>
            <edomvd>No Data</edomvd>
            <edomvds>Producer defined</edomvds>
          </edom>
        </attrdomv>
        <attrdomv>
          <rdom>
            <rdommin>61.554287227</rdommin>
            <rdommax>1381406.7803</rdommax>
            <attrunit>Kilograms per season</attrunit>
          </rdom>
        </attrdomv>
      </attr>
      <attr>
        <attrlabl>TP_kfluxkg</attrlabl>
        <attrdef>Seasonal TP load estimated using WRTDS with Kalman filtering</attrdef>
        <attrdefs>Producer Defined</attrdefs>
        <attrdomv>
          <edom>
            <edomv>&lt;&lt; empty cell &gt;&gt;</edomv>
            <edomvd>No Data</edomvd>
            <edomvds>Producer defined</edomvds>
          </edom>
        </attrdomv>
        <attrdomv>
          <rdom>
            <rdommin>2.5365186299</rdommin>
            <rdommax>1302477.4058</rdommax>
            <attrunit>Kilograms per season</attrunit>
          </rdom>
        </attrdomv>
      </attr>
      <attr>
        <attrlabl>station_id</attrlabl>
        <attrdef>Water-quality station identification number, which sometimes have lead or trailing zeros</attrdef>
        <attrdefs>Producer Defined</attrdefs>
        <attrdomv>
          <edom>
            <edomv>&lt;&lt; empty cell &gt;&gt;</edomv>
            <edomvd>No Data</edomvd>
            <edomvds>Producer defined</edomvds>
          </edom>
        </attrdomv>
        <attrdomv>
          <udom>Unique character identified</udom>
        </attrdomv>
      </attr>
      <attr>
        <attrlabl>load_00600</attrlabl>
        <attrdef>Seasonal TN load estimated using Fluxmaster</attrdef>
        <attrdefs>Producer Defined</attrdefs>
        <attrdomv>
          <edom>
            <edomv>&lt;&lt; empty cell &gt;&gt;</edomv>
            <edomvd>No Data</edomvd>
            <edomvds>Producer defined</edomvds>
          </edom>
        </attrdomv>
        <attrdomv>
          <rdom>
            <rdommin>41.804791782</rdommin>
            <rdommax>1326856.7672</rdommax>
            <attrunit>Kilograms per season</attrunit>
          </rdom>
        </attrdomv>
      </attr>
      <attr>
        <attrlabl>sload_00600</attrlabl>
        <attrdef>Standard error of load_00600</attrdef>
        <attrdefs>Producer Defined</attrdefs>
        <attrdomv>
          <edom>
            <edomv>&lt;&lt; empty cell &gt;&gt;</edomv>
            <edomvd>No Data</edomvd>
            <edomvds>Producer defined</edomvds>
          </edom>
        </attrdomv>
        <attrdomv>
          <rdom>
            <rdommin>2.6617866435</rdommin>
            <rdommax>101765.29663</rdommax>
            <attrunit>Kilograms per season</attrunit>
          </rdom>
        </attrdomv>
      </attr>
      <attr>
        <attrlabl>load_00665</attrlabl>
        <attrdef>Seasonal TP load estimated using Fluxmaster</attrdef>
        <attrdefs>Producer Defined</attrdefs>
        <attrdomv>
          <edom>
            <edomv>&lt;&lt; empty cell &gt;&gt;</edomv>
            <edomvd>No Data</edomvd>
            <edomvds>Producer defined</edomvds>
          </edom>
        </attrdomv>
        <attrdomv>
          <rdom>
            <rdommin>1.7761504457</rdommin>
            <rdommax>1277678.4249</rdommax>
            <attrunit>Kilograms per season</attrunit>
          </rdom>
        </attrdomv>
      </attr>
      <attr>
        <attrlabl>sload_00665</attrlabl>
        <attrdef>Standard error of load_00665</attrdef>
        <attrdefs>Producer Defined</attrdefs>
        <attrdomv>
          <edom>
            <edomv>&lt;&lt; empty cell &gt;&gt;</edomv>
            <edomvd>No Data</edomvd>
            <edomvds>Producer defined</edomvds>
          </edom>
        </attrdomv>
        <attrdomv>
          <rdom>
            <rdommin>0.2828189635</rdommin>
            <rdommax>422621.67978</rdommax>
            <attrunit>Kilograms per season</attrunit>
          </rdom>
        </attrdomv>
      </attr>
    </detailed>
    <detailed>
      <enttyp>
        <enttypl>streamflow_seasonal_puget.csv</enttypl>
        <enttypd>Comma Separated Value (CSV) file containing a summary of seasonal streamflow at each gage.</enttypd>
        <enttypds>Producer Defined</enttypds>
      </enttyp>
      <attr>
        <attrlabl>q_site</attrlabl>
        <attrdef>Streamflow gage identification number, which sometimes have leading or trailing zeros, from NWIS, EC, or KC data sources</attrdef>
        <attrdefs>Producer Defined</attrdefs>
        <attrdomv>
          <udom>Unique character identifier of each site</udom>
        </attrdomv>
      </attr>
      <attr>
        <attrlabl>comid</attrlabl>
        <attrdef>Stream reach common identifier (COMID) from e2NHD associated with streamflow gage</attrdef>
        <attrdefs>Producer Defined</attrdefs>
        <attrdomv>
          <rdom>
            <rdommin>23955772</rdommin>
            <rdommax>971091491</rdommax>
          </rdom>
        </attrdomv>
      </attr>
      <attr>
        <attrlabl>year</attrlabl>
        <attrdef>Calendar year</attrdef>
        <attrdefs>Producer Defined</attrdefs>
        <attrdomv>
          <rdom>
            <rdommin>2005</rdommin>
            <rdommax>2020</rdommax>
          </rdom>
        </attrdomv>
      </attr>
      <attr>
        <attrlabl>quarter</attrlabl>
        <attrdef>Quarter year or season, 1 = winter Jan Feb Mar, 2 = spring Apr May Jun, 3 = summer Jul Aug Sep, 4 = fall Oct Nov Dec</attrdef>
        <attrdefs>Producer Defined</attrdefs>
        <attrdomv>
          <rdom>
            <rdommin>1</rdommin>
            <rdommax>4</rdommax>
            <attrunit>Quarter year</attrunit>
          </rdom>
        </attrdomv>
      </attr>
      <attr>
        <attrlabl>Q_obsv_cfs</attrlabl>
        <attrdef>Mean seasonal streamflow estimated from observed mean daily streamflow from NWIS, EC, and KC data sources</attrdef>
        <attrdefs>Producer Defined</attrdefs>
        <attrdomv>
          <rdom>
            <rdommin>0.073186312125304</rdommin>
            <rdommax>26785.8358045024</rdommax>
            <attrunit>Cubic feet per second</attrunit>
          </rdom>
        </attrdomv>
      </attr>
    </detailed>
    <detailed>
      <enttyp>
        <enttypl>key_upstream_gage.csv</enttypl>
        <enttypd>Comma Separated Value (CSV) file containing a list key of every stream reach upstream of each gage.</enttypd>
        <enttypds>Producer Defined</enttypds>
      </enttyp>
      <attr>
        <attrlabl>q_site</attrlabl>
        <attrdef>Streamflow gage identification number, which sometimes have leading or trailing zeros, from NWIS, EC, or KC data sources</attrdef>
        <attrdefs>Producer Defined</attrdefs>
        <attrdomv>
          <udom>Unique character identifier</udom>
        </attrdomv>
      </attr>
      <attr>
        <attrlabl>comidg</attrlabl>
        <attrdef>Stream reach common identifier (COMID) from e2NHD associated with streamflow gage, many values are repeated</attrdef>
        <attrdefs>Producer Defined</attrdefs>
        <attrdomv>
          <rdom>
            <rdommin>23955772</rdommin>
            <rdommax>947110179</rdommax>
          </rdom>
        </attrdomv>
      </attr>
      <attr>
        <attrlabl>comidu</attrlabl>
        <attrdef>Stream reach common identifier (COMID) from e2NHD of all features upstream of each streamflow gage</attrdef>
        <attrdefs>Producer Defined</attrdefs>
        <attrdomv>
          <rdom>
            <rdommin>23955772</rdommin>
            <rdommax>971091491</rdommax>
          </rdom>
        </attrdomv>
      </attr>
    </detailed>
    <detailed>
      <enttyp>
        <enttypl>flow_list.csv</enttypl>
        <enttypd>Comma Separated Value (CSV) file containing a list of streamflow gages.</enttypd>
        <enttypds>Producer Defined</enttypds>
      </enttyp>
      <attr>
        <attrlabl>q_site</attrlabl>
        <attrdef>Streamflow gage identification number, which sometimes have leading or trailing zeros, from NWIS, EC, or KC data sources</attrdef>
        <attrdefs>Producer Defined</attrdefs>
        <attrdomv>
          <udom>Unique character identifier for each site</udom>
        </attrdomv>
      </attr>
      <attr>
        <attrlabl>q_site_nm</attrlabl>
        <attrdef>Streamflow gage name from NWIS, EC, or KC data sources</attrdef>
        <attrdefs>Producer Defined</attrdefs>
        <attrdomv>
          <udom>Unique character name for each site</udom>
        </attrdomv>
      </attr>
      <attr>
        <attrlabl>q_lat</attrlabl>
        <attrdef>Streamflow gage latitude in decimal degrees from NWIS, EC, or KC data sources</attrdef>
        <attrdefs>Producer Defined</attrdefs>
        <attrdomv>
          <rdom>
            <rdommin>46.74427397</rdommin>
            <rdommax>48.9387257</rdommax>
            <attrunit>Decimal degrees</attrunit>
          </rdom>
        </attrdomv>
      </attr>
      <attr>
        <attrlabl>q_lon</attrlabl>
        <attrdef>Streamflow gage longitude in decimal degrees from NWIS, EC, or KC data sources</attrdef>
        <attrdefs>Producer Defined</attrdefs>
        <attrdomv>
          <rdom>
            <rdommin>-124.3838328</rdommin>
            <rdommax>-121.0728964</rdommax>
            <attrunit>Decimal degrees</attrunit>
          </rdom>
        </attrdomv>
      </attr>
      <attr>
        <attrlabl>q_source</attrlabl>
        <attrdef>Data source of streamflow, NWIS, EC, or KC</attrdef>
        <attrdefs>Producer Defined</attrdefs>
        <attrdomv>
          <edom>
            <edomv>NWIS</edomv>
            <edomvd>U.S. Geological Survey observed streamflow</edomvd>
            <edomvds>Producer defined</edomvds>
          </edom>
        </attrdomv>
        <attrdomv>
          <edom>
            <edomv>EC</edomv>
            <edomvd>Washington State Department of Ecology observed streamflow</edomvd>
            <edomvds>Producer defined</edomvds>
          </edom>
        </attrdomv>
        <attrdomv>
          <edom>
            <edomv>KC</edomv>
            <edomvd>King County observed streamflow</edomvd>
            <edomvds>Producer defined</edomvds>
          </edom>
        </attrdomv>
      </attr>
      <attr>
        <attrlabl>comid</attrlabl>
        <attrdef>Streamflow gage paired common identifier (COMID) of closest associated stream reach from e2NHD</attrdef>
        <attrdefs>USGS Defined</attrdefs>
        <attrdomv>
          <rdom>
            <rdommin>23955772</rdommin>
            <rdommax>971091491</rdommax>
          </rdom>
        </attrdomv>
      </attr>
    </detailed>
    <detailed>
      <enttyp>
        <enttypl>fluxmaster_puget.zip</enttypl>
        <enttypd>Fluxmaster model copy, see Process Steps</enttypd>
        <enttypds>Producer Defined</enttypds>
      </enttyp>
    </detailed>
    <detailed>
      <enttyp>
        <enttypl>wrtds_puget.zip</enttypl>
        <enttypd>WRTDS model copy, see Process Steps</enttypd>
        <enttypds>Producer Defined</enttypds>
      </enttyp>
    </detailed>
    <detailed>
      <enttyp>
        <enttypl>wtemp_seasonal_puget.csv</enttypl>
        <enttypd>Comma Separated Value (CSV) file containing a summary of seasonal water temperature at select streamflow gages and water-quality stations.</enttypd>
        <enttypds>Producer Defined</enttypds>
      </enttyp>
      <attr>
        <attrlabl>wtemp_site</attrlabl>
        <attrdef>Site identification number of observed water temperature, which sometimes have leading or trailing zeros, from NWIS, EC, norwest, or WQP data sources</attrdef>
        <attrdefs>Producer Defined</attrdefs>
        <attrdomv>
          <udom>Unique character identifiers for each site</udom>
        </attrdomv>
      </attr>
      <attr>
        <attrlabl>comid</attrlabl>
        <attrdef>Stream reach common identifier (COMID) from e2NHD associated with observed water temperature</attrdef>
        <attrdefs>USGS Defined</attrdefs>
        <attrdomv>
          <rdom>
            <rdommin>23955772</rdommin>
            <rdommax>947110181</rdommax>
          </rdom>
        </attrdomv>
      </attr>
      <attr>
        <attrlabl>year</attrlabl>
        <attrdef>Calendar year</attrdef>
        <attrdefs>Producer Defined</attrdefs>
        <attrdomv>
          <rdom>
            <rdommin>2005</rdommin>
            <rdommax>2020</rdommax>
            <attrunit>year</attrunit>
          </rdom>
        </attrdomv>
      </attr>
      <attr>
        <attrlabl>quarter</attrlabl>
        <attrdef>Quarter year or season, 1 = winter Jan Feb Mar, 2 = spring Apr May Jun, 3 = summer Jul Aug Sep, 4 = fall Oct Nov Dec</attrdef>
        <attrdefs>Producer Defined</attrdefs>
        <attrdomv>
          <rdom>
            <rdommin>1</rdommin>
            <rdommax>4</rdommax>
            <attrunit>quarter year</attrunit>
          </rdom>
        </attrdomv>
      </attr>
      <attr>
        <attrlabl>wtemp_mean</attrlabl>
        <attrdef>Mean seasonal water temperature in degrees Celsius computed from mean daily observations</attrdef>
        <attrdefs>Producer Defined</attrdefs>
        <attrdomv>
          <rdom>
            <rdommin>3.734529398</rdommin>
            <rdommax>19.02380952</rdommax>
            <attrunit>degrees Celsius</attrunit>
          </rdom>
        </attrdomv>
      </attr>
      <attr>
        <attrlabl>wtemp_sd</attrlabl>
        <attrdef>Standard deviation of daily observations used to compute mean seasonal water temperature in degrees Celsius</attrdef>
        <attrdefs>Producer Defined</attrdefs>
        <attrdomv>
          <rdom>
            <rdommin>0.25663817</rdommin>
            <rdommax>4.245584838</rdommax>
            <attrunit>degrees Celsius</attrunit>
          </rdom>
        </attrdomv>
      </attr>
      <attr>
        <attrlabl>wtemp_cnt</attrlabl>
        <attrdef>Number of daily observations available and used to compute a mean seasonal water temperature</attrdef>
        <attrdefs>Producer Defined</attrdefs>
        <attrdomv>
          <rdom>
            <rdommin>52</rdommin>
            <rdommax>92</rdommax>
            <attrunit>Integer count</attrunit>
          </rdom>
        </attrdomv>
      </attr>
    </detailed>
    <detailed>
      <enttyp>
        <enttypl>dataDictionary_Fluxmaster_WRTDS.csv</enttypl>
        <enttypd>Comma Separated Value (CSV) file containing data describing files in copies of WRTDS and Fluxmaster.</enttypd>
        <enttypds>Producer Defined</enttypds>
      </enttyp>
      <attr>
        <attrlabl>folderName</attrlabl>
        <attrdef>The name of folder containing file</attrdef>
        <attrdefs>U.S. Geological Survey</attrdefs>
        <attrdomv>
          <udom/>
        </attrdomv>
      </attr>
      <attr>
        <attrlabl>fileName</attrlabl>
        <attrdef>The name of the file</attrdef>
        <attrdefs>U.S. Geological Survey</attrdefs>
        <attrdomv>
          <udom/>
        </attrdomv>
      </attr>
      <attr>
        <attrlabl>fileDesc</attrlabl>
        <attrdef>A description of the file</attrdef>
        <attrdefs>U.S. Geological Survey</attrdefs>
        <attrdomv>
          <udom/>
        </attrdomv>
      </attr>
      <attr>
        <attrlabl>varName</attrlabl>
        <attrdef>The name of the variable within the file</attrdef>
        <attrdefs>U.S. Geological Survey</attrdefs>
        <attrdomv>
          <udom/>
        </attrdomv>
      </attr>
      <attr>
        <attrlabl>varDef</attrlabl>
        <attrdef>A description of the variable within the file</attrdef>
        <attrdefs>U.S. Geological Survey</attrdefs>
        <attrdomv>
          <udom/>
        </attrdomv>
      </attr>
      <attr>
        <attrlabl>varMin</attrlabl>
        <attrdef>The minimum value of the variable if numeric or date</attrdef>
        <attrdefs>U.S. Geological Survey</attrdefs>
        <attrdomv>
          <udom/>
        </attrdomv>
      </attr>
      <attr>
        <attrlabl>varMax</attrlabl>
        <attrdef>The maximum value of the variable if numeric or date</attrdef>
        <attrdefs>U.S. Geological Survey</attrdefs>
        <attrdomv>
          <udom/>
        </attrdomv>
      </attr>
      <attr>
        <attrlabl>varUnits</attrlabl>
        <attrdef>The units of the variable, if applicable</attrdef>
        <attrdefs>U.S. Geological Survey</attrdefs>
        <attrdomv>
          <udom/>
        </attrdomv>
      </attr>
      <attr>
        <attrlabl>varCharVal</attrlabl>
        <attrdef>The character value of the variable, if a character</attrdef>
        <attrdefs>U.S. Geological Survey</attrdefs>
        <attrdomv>
          <udom/>
        </attrdomv>
      </attr>
      <attr>
        <attrlabl>varCharValDef</attrlabl>
        <attrdef>The definition of the character value</attrdef>
        <attrdefs>U.S. Geological Survey</attrdefs>
        <attrdomv>
          <udom/>
        </attrdomv>
      </attr>
    </detailed>
  </eainfo>
  <distinfo>
    <distrib>
      <cntinfo>
        <cntperp>
          <cntper>GS ScienceBase</cntper>
          <cntorg>U.S. Geological Survey</cntorg>
        </cntperp>
        <cntaddr>
          <addrtype>mailing address</addrtype>
          <address>Denver Federal Center, Building 810, Mail Stop 302</address>
          <city>Denver</city>
          <state>CO</state>
          <postal>80225</postal>
          <country>United States</country>
        </cntaddr>
        <cntvoice>1-888-275-8747</cntvoice>
        <cntemail>sciencebase@usgs.gov</cntemail>
      </cntinfo>
    </distrib>
    <distliab>Unless otherwise stated, all data, metadata and related materials are considered to satisfy the quality standards relative to the purpose for which the data were collected. Although these data and associated metadata have been reviewed for accuracy and completeness and approved for release by the U.S. Geological Survey (USGS), no warranty expressed or implied is made regarding the display or utility of the data on any other system or for general or scientific purposes, nor shall the act of distribution constitute any such warranty. Furthermore, the database is released on condition that neither the USGS nor the U.S. Government shall be held liable for any damages resulting from its authorized or unauthorized use. Although the software has been subjected to rigorous review, the USGS reserves the right to update the software as needed pursuant to further analysis and review. No warranty, expressed or implied, is made by the USGS or the U.S. Government as to the functionality of the software and related material nor shall the fact of release constitute any such warranty. Furthermore, the software is released on condition that neither the USGS nor the U.S. Government shall be held liable for any damages resulting from its authorized or unauthorized use. Any use of trade, product, or firm names is for descriptive purposes only and does not imply endorsement by the U.S. Government.</distliab>
    <stdorder>
      <digform>
        <digtinfo>
          <formname>Digital Data</formname>
        </digtinfo>
        <digtopt>
          <onlinopt>
            <computer>
              <networka>
                <networkr>https://doi.org/10.5066/P9LY1PQF</networkr>
              </networka>
            </computer>
          </onlinopt>
        </digtopt>
      </digform>
      <fees>None</fees>
    </stdorder>
  </distinfo>
  <metainfo>
    <metd>20260409</metd>
    <metc>
      <cntinfo>
        <cntorgp>
          <cntorg>U.S. Geological Survey</cntorg>
          <cntper>Noah M. Schmadel</cntper>
        </cntorgp>
        <cntpos>Hydrologist</cntpos>
        <cntaddr>
          <addrtype>mailing and physical</addrtype>
          <address>601 SW 2nd Ave. Suite 1950</address>
          <city>Portland</city>
          <state>OR</state>
          <postal>97204</postal>
          <country>United States</country>
        </cntaddr>
        <cntvoice>9714421691</cntvoice>
        <cntemail>nschmadel@usgs.gov</cntemail>
      </cntinfo>
    </metc>
    <metstdn>FGDC Content Standard for Digital Geospatial Metadata</metstdn>
    <metstdv>FGDC-STD-001-1998</metstdv>
  </metainfo>
</metadata>
