United States
                    Environmental Protection
                    Agency
                                            Office Of Water
                                            (WH-553)
EPA841-F-93-004
November 1992
Number 4
                    TMDL  Case  Study
                    Nomini  Creek  Watershed
 Key Feature:
 Project Name:
 Location:
.- .- r • ',*, 'jf\ -. . •'.•:"„" ,:.-.
 Scope/Size:    '*';*.•••
 Land Type: •-'•-

 Type of Activity:   "
 Pollutants: ..-...;•--,.
 TMDL Development:
 Data Sources:
 Data Mechanisms:

 Monitoring Plan:

 Control Measures:
                 Use of a geographic information system
                 and watershed models to identify areas of
                 critical nonpoint pollution.

                 Nomini Creek Watershed GIS Study
                 USEPA Region Ill/Westmoreland
                -County, Virginia/Potomac River
                - Small watershed, 1505 hectares
                 Ecoregion 65, Southeastern   „
                 plains;        -..•;.'
                 Agriculture   ->.  •     •• .. .-":..
                 Nutrients, sediment      ,
                                 '
                              .             .
                 State, locairfederal    ~J :,-: ,-'
                 Modeling (SLOSS, PHOSPH); GIS
                 (VirGIS)
                 Yes, long-term BMP effectiveness
                 monitoring:         ;:
                 BMPs   :. , .       V ."-...;;.
Summary: Using the Nomini Creek watershed (Figure 1),
Virginia's Department of Conservation and Recreation and its
Polytechnic and State University demonstrated how geographic
information system (GIS) technology can be used to (1) prioritize
and target waterbodies with multiple water quality concerns, and
(2) target BMPs to critical nonpoint source loading areas to meet
load allocations more effectively.
                                                    Nomjni Greek
                                                      Watershed
  Chesapeake Bay
              Inset
              Area
                                                  FIGURE 1. Location of the Nomini Creek watershed
                                                  in Virginia  ,                    \   '
The Department of Conservation and Recreation selected the Nomini Creek watershed as an area in which to evaluate and
monitor the effectiveness of best management practices (BMPs) for the Chesapeake Bay Program. To identify the critical
phosphorus and sediment loading areas within the watershed so that BMPs could be sited effectively, the Division tested
the feasibility of integrating VkGIS, a state-run GIS, with two simple pollutant yield models (SLOSS and PHOSPH).
Because Virginia's data base was sufficiently large, VkGIS  was able to provide the data requked to run the models.  The
output from these models successfully identified critical areas of nonpoint source loading.  BMPs were sited on these areas
and an intensive water quality monitoring program is currently in place to evaluate BMP effectiveness and to verify the
estimated pollutant loads.                                              .
Contact:
           Jf, Michael Flagg, Vkgiaia DCR-DSWC, 203 Governor Street, Suite 206, Richmond, VA 23219, phone
           {804)786-3959                      ;        .-   -"

-------
 BACKGROUND

 Programmatic Issues

 EDITOR'S NOTE:  Water quality management in the
 Chesapeake Bay employs the essential elements of a
 TMDL assessment.  Stakeholders recognized water
 quality problems, identified the causes and sources  of
 these problems, set an achievable target or endpoint, and
 then targeted controls for specific point and nonpoint
 pollution sources to decrease their contributions.
 Technical efforts to address difficulties or obstacles to
 effective water quality management in the Chesapeake
 basin therefore have a unique place within the TMDL
 Case Study series.

: la the early 1970s, water quality managers recognized
 that the water quality of the Chesapeake had become
 substantially degraded and prevented the Bay from  fully
 meeting its designated uses.  Subsequent studies of the
 Chesapeake Bay, sponsored by USEPA, identified
, nutrients as a primary cause of the depleted oxygen and
. other water quality problems. USEPA, the District of
 Columbia, and the states of Pennsylvania, Maryland and
 Virginia agreed that reducing total nutrient loads to the
 Bay by 40 percent was an achievable goal and would
 improve the health of the ecosystem.  As  the primary
 source of nutrients in the basin, agricultural activities
 were targeted for implementation of BMPs.

 To begin to address these nonpoint source problems,
 Virginia, Maryland, Pennsylvania, and the District  of
 Columbia have initiated watershed programs to
 demonstrate and monitor the benefits of various BMPs.
 The new information on BMP effectiveness will
 eventually be used to encourage more widespread use of
 BMPs throughout the Chesapeake drainage to achieve the
 40 percent reduction in nutrient loads.

  In 1985, the Virginia Department of Conservation  and
  Recreation, Division of Soil and Water Conservation
  (DCR-DSWC or the Division) selected the Nomini Creek
  watershed as an area in which to evaluate and monitor
  BMP effectiveness.  The watershed was chosen because
  of its proximity to the Bay and because agriculture is its
  predominant land use activity.  The watershed is also
  typical of row-cropped agricultural areas  hi the Virginia
  coastal plain. Agricultural loadings of nitrogen,
  phosphorus,  and sediment from Nomini Creek enter the
  lower Potomac River, which is a main tributary to the
  Chesapeake  Bay.

  To site BMPs effectively, the Division sought to identify
  critical nutrient loading areas within the watershed.  The
  Virginia Geographic Information System (VirGIS), in
  conjunction  with simple pollutant loading models,  was
  proposed to accomplish this task.  Although the GIS was
developed to facilitate these types of projects, no one had
used it yet for this purpose. The Division and the
Department of Agricultural Engineering at Virginia
Polytechnic Institute and State University therefore
conducted a sub-study to determine whether critical
nutrient loading areas could be properly identified when
VirGIS was linked with the models.

The Resource

The Nomini Creek watershed is located in Westmoreland
County in eastern Virginia. Its 3,719 acre watershed is
approximately 54 percent woodlands (there is a
significant amount of commercial forestry), 43  percent
croplands, and 3 percent homesteads and roads. There
are'no towns in the watershed. The primary agricultural
activity in the watershed is row cropping of com,  barley,
and wheat (USDA SCS, 1979).  Twenty-five farmers
cultivate the 1,646 acres of cropland in the watershed.
Thirty-three percent (546 acres) of the farmland is
farmed by the landowners/operators themselves, while
the remaining 67 percent (1100 acres) is leased for others
to farm (DCR-DSWC, 1986). One small beef cattle
operation also exists in the watershed.

Nomini Creek currently meets all of its designated uses.
Biological parameters indicate that water quality within
the watershed is good; however, preliminary water
quality sampling results indicate a high level of nutrients.
In 1986, gross soil erosion was calculated at 7,920
tons/year for the study portion of the watershed.
Approximately 1,584 tons of that sediment is delivered to
Nomini Creek, along with 10,700 pounds of phosphorus.
Eroding cropland is responsible for about 95 percent of
 the total sediment and phosphorus that is delivered to the
 creek (DCR-DSWC, 1986).  The state, in conjunction
 with the Virginia Forestry Association, currently
 monitoring to determine how much forestry practices
 contribute to the pollution load.

 Hydrologically, the watershed contains first- and second-
 order streams flowing through nearly level to gently
 sloping topography.  The soils are mostly fine sandy
 loam to loam, with slopes that range from nearly level to
 15 percent, except along stream banks where slopes
 range from 15  percent to 50 percent (USDA SCS, 1979).
 Approximately 40 percent of the cropland in the
 watershed is classified as highly erodible.
 ASSESSING AND CHARACTERIZING
 PROBLEM AREAS

 GIS and Modeling Tools

 The data necessary for characterizing nutrient loading
 patterns within the Nomini Creek watershed were

-------
 available on VirGIS, a state-run GIS. A GIS is a
 computerized system for storing and manipulating data
 that have a spatial component (i.e., data for which
 geographical location is important).  VirGIS was initiated
 in 1985 by the Division as a tool for developing
 modeling and mapping procedures that could readily
 identify land areas with nonpoint source pollution
 potential.
 can be modified to estimate potential soil erosion for
 each hydrologically homogeneous cell in a watershed.
 For Nomini Creek, a cell size of 1/9 hectare was used.

 Three main equations are used in SLOSS.The first
 equation computes soil loss per unit area of watershed
 (A,).  It is expressed as:
 VirGIS is a" modular and highly interactive program that
 consists of a large database coupled with nearly 500
 special-purpose programs that manipulate and display
 data. Basic data types, or layers, taken from 7.5-minute         where
 quadrangle maps, county soU, surveys, National High
 Altitude Program color-infraredphotos, and U.S.
 Geological Survey (USGS) digital elevation models '
 include elevation, soils, land use, surface water,  _;
 watersheds, and county boundaries.  The data are stored
 primarily in raster form (i.e., as small area! units) with
 each unit, or cell, representing from 1/9  to 1 hectare.
 These units are joined to form data layers that cover
 areas ranging from 7.3 million hectares up to 1.0.1
 million  hectares (the entire State of Virginia).  VirGIS
 can use the basic data layers to calculate "derived" data
 layers (Figure 2) for input parameters required by
 nonpoint source pollutant yield models.

 Virginia Polytechnic and State University and the state
 developed two simple nonpoint source pollutant yield
 models, SLOSS and PHOSPH, to characterize the               where
Nomini Creek watershed.

 SLOSS  is a simplified pollutant yield model designed to
estimate soil loss and sediment delivery to a stream. It is             '
based on the Universal Soil Loss Equation (USLE) and
                                                (i)
Q
P£
            = maximum number of cells;
             = soil credibility factor;
             = topographic factor;
             = land use/land cover management; and
             = support practice factor.
The SLOSS model then calculates a sediment delivery.
ratio for each cell.  This ratio relates the amount of
sediment lost, from a cell to the amount that will actually
be delivered to the stream channel.  The delivery ratio
(DR) is expressed as:
                                               (2)
        b   = land cover factor;
        Lf  = length of the flow path between cell i
                and the channel outlet; and
        Sf  = slope function.              ,.
                       Soils
                     Land Use
                    Elevation
                  Watersheds
                  Surface Waters

                   Boundaries
                      SCWD     ^
                  Boundaries
           BASE DATA LAYERS
                                                                                      Slope
                                                                                  Slope Length
                       Slope Length
                           Factor
                        Erodibility
                                                                                      Factor
                                                                                    Tolerance
                                                                                      Factor
                         Delivery
                           Factor
                      Water Quality
                      	Index
                                                                                     Erosion
                                                                                    	Index
               DERIVED DATA LAYERS
                     FIGURE 2. Basic and derived data layers used in the VirGIS database
                                                                                                            3

-------
Total sediment yield (L.) is then calculated in the
expression:
                      1-1
where the parameters are as previously defined.

PHOSPH is a simplified phosphorus loading model
developed to circumvent the intense data requirements of
more complex phosphorus models. The basic equation in
PHOSPH is:
                                                 (4)
                    l-i
 where
        TP, = total sediment-associated phosphorus
               .delivered to the stream outlet;
        Pq  * average phosphorus content of the
               surface soil layer for soil hi cell i;
        L,   s-sediment yield for each cell (eq. 3);
               and
        ER_ = phosphorus enrichment ratio.   "


 The phosphorus enrichment ratio is defined as the mass
 of phosphorus in the eroded sediment per unit mass of
•phosphorus hi the surface soil layer.  It is calculated hi
 the equation:
                                   -039
                                                 (5)
where CmiK and Cmln  are the maximum and minimum
percent clay content of the soil in each cell.

Locating NFS Problem Areas

The first step in the modeling process was to obtain the
necessary data layers from the VirGIS system.  This was
done  by creating a data window for Nomini Creek to
obtain only the relevant data from the much larger
VirGIS data base.  Once these data were extracted,
several VirGIS programs were used to convert the
.information into parameters that would be accepted by
the models (e.g., the topographic factor, LS, was
determined from VirGIS slope uiformation).  The SLOSS
model was then used with these parameters to estimate
sediment loss.  The sediment loss  estimate then was
passed to the PHOSPH model, which calculated
sediment-associated phosphorus export. This process is
summarized in Figure 3.

The spatial resolution of the VirGlS data for Nomini
Creek was 1/9 hectare. In other words, each type of data
was recorded for blocks of land 1/9 hectare in area.
Consequently, individual sediment and phosphorus export
values were calculated for each of approximately 33,470
individual cells within the watershed.  The numerical
results from both models were reprocessed by VirGIS
 into maps that allowed easy comparison of loadings
 throughout the watershed.  Figure 4 shows the GIS-
 generated jtnap for phosphorus yield in the watershed.

 It was determined that estimated sediriient yield exceeded
 the "high" threshold of 22.4 tons/hectare/year in
 approximately 15 percent  (227 hectares) of the Nomini
 Creek watershed.  Due to the lack of established values
 hi the literature, the sediment threshold was arbitrarily
        Bute Data
         Layers
                            GIS
                                         Derived Input
                                          Parameters
       Etevalkxi
        Cover
                                           Pollutant
                                          Yield Maps


Lty
i '


1
J
^


— >•




"K", "R", "LS"...




— >



    NPS Models
                                                                                                Sediment Loading
              FIGURE 3.  Conceptual framework for the integration of a GIS with water quality models

-------
      ; I?B/ACIYR
      10.00 - 200.00
         > 200.00
  Seal* 1:30000 C«ISiz«. I/SNA(057AC)
FIGURE 4.  Phosphorus yield  in the  Nomini  Creek
watershed

determined,in order to/facilitate classification into high,
medium, or low categories. For phosphorus,
approximately 21 percent of the watershed (316 hectares)
exceeded the high threshold value of 1.12
kg/hectare/year.  The phosphorus threshold was obtained
from extensive literature review.
SITING NFS CONTROLS

The critical area maps generated by VirGIS assisted in
siting BMPs where they were needed in the watershed.
The BMPs included no-till fanning, nutrient management
plans, grassed waterways, drop structures, diversions,
pasture management, and the removal of land from
production. The installation of the BMPs was
accomplished through the Virginia Agricultural BMP
Cost-Sharing Program, in cooperation with the U.S.
Department of Agriculture's Soil Conservation Service  •
and Agricultural Stabilization arid Conservation Service,
the Virginia Department of Forestry, the Cooperative
Extension Service, and the Northern Neck Soil and
Water Conservation District.
                                                            FOLLOW-UP MONITORING

                                                            Intensive water quality monitoring, which is the primary
                                                            focus of the Division's efforts in the watershed, is
                                                            continuing and will be used to verify the estimated
                                                            loadings of sediment and phosphorus.  This monitoring
                                                            includes both storm event and ambient monitoring. Storm
                                                            sampling is conducted whenever the stream level
                                                            increases by 2/10 of a fop.t Storm event parameters
                                                            include total nitrogen, total phosphorus, nitrate, ammonia,
                                                            phosphate, and total suspended solids. Ambient
                                                            monitoring occurs on a monthly basis for protozoans,
                                                            pesticides, and a suite of weather parameters.  Colifprm
                                                            bacteria are sampled every 2 weeks, and land use data
                                                            are collected twice per year. Monthly groundwater
                                                            monitoring is also conducted for pesticides, nutrients, and
                                                            water table depth.
 MANAGEMENT CONSIDERATIONS

 The maps produced by the modeling process are
 especially valuable as watershed management tools.
 First, they highlight the portions of a watershed that are
 critical hi terms of their sediment and;phosphofus loading
 potential.  (The critical loading-level can be ,set by the
 user based on literature review or professional judgment)
 This feature provides managers with a readily
 understandable means for determining areas that are in
 need of control measures.  For states, this information
 can assist in targeting areas to receive BMP cost-share
. funding, potentially increasing the cost-effectiveness'of"
 existing state cost-share programs.

 In addition, GIS output maps ate useful as an education
 tool. Landowner cooperation is sometimes a difficult
 obstacle hi the implementation of BMP cost-share
 programs. These maps could be used by agricultural
 field personnel as visual means for promoting program
 cooperation.               7

 Targeting high-priority waterbodies or watersheds for
 TMDL development often involves more than just
 technical factors. It may also involve the evaluation of
 factors related to recreational, economic, and ecological
 values such as the risk to human health and aquatic life;
 the degree of public interest and support in protecting a
 waterbody; the recreational, economic, and aesthetic
 importance of the waterbody; and its vulnerability or
 fragility as aquatic habitat  Many of these factors
 contain spatial components and can be displayed on
 maps.  Overlaying these maps on maps of potential
 pollutant yield would  illustrate which waterbodies are of
 special concern. Coupled with professional judgement,

-------
these maps could facilitate the prioritization and targeting
of watersheds with the greatest need for TMDL
development.

When sufficient data are collected, the Nornini Creek
Study will provide useful information on the
effectiveness of BMPs in reducing NFS pollution on an
entire watershed in the coastal plain. Furthermore, if the
statistical analyses show a significant instream
phosphorus reduction, this experiment will have provided
a new tool for managers to use in achieving the
Chesapeake Bay 40 percent nutrient reduction goal.  It is
important to note, however, that the success of this
technique in the coastal plain will not necessarily make it
a viable modeling tool in other parts of the Chesapeake
drainage. The ten-year monitoring period for Nomini
Creek will be complete in 1995.
REFERENCES

DCR-DSWC.  1986.  Nomini Creek Watershed Plan
Westmoreland County, Virginia. Virginia Department of
Conservation and Recreation, Division of Soil and Water
Recreation.  Revised 6/10/86.

Omernik, J.M.  1987. Ecoregions of the conterminous
United States.  Annals of the Association of American
Geographers 77(1):118-125.
Shanholz, V.O., C.J. Desai, N. Zhang, J.W. Kleene, and
C.D. Metz.  1990.  Hydrologic/water quality modeling in
a GIS environment.  Written for presentation at the 1990
International Summer Meeting sponsored by the
American Society of Agricultural Engineers, Columbus,
OH, June 24-27, 1990.  Paper No. 90-3033.

USDA SCS.  1979.  Westmoreland county soil survey.
Soil Conservation Service National Cooperative Soil
Survey, U.S. Department of Agriculture, Washington,
DC.

USEPA.  1991. Guidance for water quality-based
decisions: The TMDL process. U.S. Environmental
Protection Agency, Office of Water, Washington, DC.
EPA 440/4-91-001.
   Tb4s  case  study was prepared  by Research  Triangle
   Institute, Research; triangle Park, KG, in coajuttttioA with
   USEPA, Office  of Office  of Wetlands,  Oceans,  and
   Watersheds, Watershed Mattag0ment: Section, To obtain
   copies,  contact  your EPA  Regional  303(d)/lMDL
   Coordinatofv

-------