USGS - science for a changing world

Northern Prairie Wildlife Research Center

  Home About NPWRC Our Science Staff Employment Contacts Common Questions About the Site

Influences of Fathead Minnows on Phosphorus and Nitrogen Cycling in a Prairie Wetland

Methods


We estimated prey consumption and nutrient budgets for fathead minnows by using the bioenergetic models of Hewett and Johnson (1987) and Kraft (1992). This technique involves modeling individual cohorts through time using growth and mortality rates specific for each cohort. Variables of interest (i.e. consumption and excretion rates) are then modeled for each cohort during the sampling period. Results for all cohorts can then be compiled to estimate variables of interest at the total population level. Use of these models requires collection of field data on population characteristics such as mortality, growth, diet composition, and water temperature, as well as estimating species-specific physiological variables in a laboratory setting.
GIF-Study location in Minnesota.
Figure 1. Location of the study site in west-central Minnesota.

Collection of Field Data

Our study wetland, Sagebraten, is located on a USFWS Waterfowl Production Area (WPA) in west-central Minnesota (Figure 1) and has a surface area of 3.71 ha and a maximum depth of 2.46 m. We established two sampling strata consisting of deep (water depth > 1.25 m) and shallow (depth < 1.25 m) water. We used an area-density method to estimate the population size of fathead minnows (Duffy 1998), and this technique requires estimation of the surface area of both strata. A Topcon® total station was used to generate data for a basin-morphometry map for the study site in mid-July, and the area of each depth-stratum was determined using Surfer software (Golden Software 1997).

Fathead minnows were sampled every two weeks from 19 May through 14 August 1997 using square-shaped pop nets (Dewey et al. 1989) with a sampling area of 1 m2 and fitted with 0.5 mm (bar measure) mesh. On each date 15 samples were taken at random locations at wadeable depths in each stratum and the total number of fish captured per set was determined. We then determined the average number of fish per m2 for each stratum from the 15 samples, and estimated the total population size by summing the average density multiplied by area for each stratum.

A substantial portion (39%) of our study wetland was too deep to be sampled without use of a boat, and considerable error can be introduced into the population estimates if densities vary between the sampled and unsampled area of the deep stratum. Therefore, we sampled this deeper water on 15 July to determine if densities were significantly different than those observed in shallower water of the deep stratum. Fifteen samples were taken and compared to 15 samples taken on the same day in shallower water of the deep stratum. As no difference was found between the two depths (T-test; t=2.2, 28 df, P=0.35), we concluded that the shallower samples accurately estimate densities for the entire wetland.

Use of bioenergetic models requires estimating the growth of fish during the time period being modeled, which is usually done by identifying and tracking the growth of individual cohorts. However, standard techniques cannot be used to identify cohorts in fathead minnow populations, as they are fractious spawners and produce multiple cohorts throughout the summer (Duffy 1998). Thus, we used size-frequency data to identify and track the growth of cohorts throughout the sampling season. On each sampling date the total length and wet-mass were determined for 300 randomly selected minnows. Modal lengths were identified, and fish within + 20% of a mode were associated as a cohort (Duffy 1998). The average weight of fish in each cohort was determined on each date, and growth of individuals in each cohort was determined by the increase in average weight between one sampling date and the next.

We analyzed the diet of fathead minnows by collecting 15 adult (total length > 40 mm) and 15 juvenile (total length < 40 mm) fish on each sampling date and immediately preserving them in 10% buffered formalin. In the lab, the anterior 1/3 of each stomach was removed, and it's contents enumerated and identified to the lowest feasible taxonomic group. Lengths of prey were determined and converted to biomass with length-mass equations of McCauley (1984) and Smit et al. (1993). The percent of total diet by weight was then determined for each prey type. Results for fish collected in each month were pooled, giving an average diet composition for the months of May, June, July and August. We were unable to determine diets of fish < 30 mm in length, but Held and Peterka (1974) found that fish < 20 mm in length fed mainly on cladocerans and copepods. Therefore, we assumed all YOY fathead minnows fed exclusively on copepod prey (cladocerans were rarely found in stomachs of fish or in the water column during our study) during the first two weeks of the model (when their total lengths ranged from 13-19 mm) and utilized the observed diet data after two weeks (when their total length was > 19 mm).

Metabolic rates of fish are temperature dependant, making water temperature data necessary for bioenergetic models. We measured water temperature during the sampling period with temperature-data loggers (Hobo Inc., Procasse, MA). The data logger was deployed 50 cm off the bottom at 1 m depth, and recorded water temperature every 4 hr. The average daily water temperature was then incorporated into the bioenergetics model.

Physiological Variables and Bioenergetic Modeling

Use of the bioenergetic model to estimate nutrient dynamics in fish requires constructing energy and nutrient budgets. The energy budget of fish is based on the energy mass-balance equation

G=C-(R+S+F+U)
where G is observed growth, C is energy ingested, R is energy used in metabolism, S is specific dynamic heat increment, F is egestion, and U is excretion (Hewett and Johnson 1992). Values of variables in the equation are based on a number of physiological parameters that are species specific. We used values developed specifically for fathead minnows by Duffy (1998), and these are shown in Table 1. Use of this equation also requires data on predator and prey energy concentrations. Energy concentrations for fathead minnow were taken from Duffy (1998), and invertebrate prey from Cummins and Wuychek (1971), Hanson et al. (1997), and Schindler and Eby (1997). When values were not available for specific prey, we used values for the most closely related taxonomic group.

In a similar fashion, consumption, allocation to growth, and excretion of both phosphorus and nitrogen can be estimated with the mass-balance equation

C=G+F+U
where C is mass of nutrients consumed, G is nutrients allocated to growth, F is nutrients lost in feces, and U is nutrients lost in urine (Kraft 1992). We used phosphorus and nitrogen assimilation efficiencies of 0.72 and 0.80, respectively (Nakashima and Leggett 1980), and values of F and U reported in Duffy (1998). We assumed that 50% of phosphorus in fish tissue would be quickly released upon fish death and decomposition, with the remaining tied up for longer periods in bone and scales (Kitchell et al. 1975, Smith et al. 1977). Nutrient modeling also requires data on nitrogen and phosphorus concentrations of prey, and we used data in Hanson et al. (1997) for nutrient concentrations in fathead minnows and invertebrate prey. We again used values for the most closely related taxonomic group if data were not available for a specific prey type.

We used sensitivity analysis to determine how errors in the physiological parameters would effect model estimates. Sensitivity was assessed with the equation

GIF-Equation used to assess sensitivity.
where Sx (p) is the sensitivity of the model estimate x to perturbations in the physiological parameter p, p is the original value of parameter p, p is the deviation of parameter p, X is the original value of the model estimate, and X is the deviation of the model estimate X due to p (Kitchell et al. 1977, Bartell et al. 1986). The analysis was performed by deviating individual physiological parameters +10% the value used in our model and determining the percent change in model estimates of phosphorus consumption and excretion relative to our original estimates.

All modeling was performed with Fish Bioenergetics 3.0 © (Hanson et al. 1997). All calculations were performed on a daily interval, and models for each cohort were combined for total population estimates. The model was used to estimate consumption rates of prey, nutrients contained in the fish population, nutrients released via fish mortality, and nutrient consumption, allocation to growth, and excretion.


Previous Section -- Introduction
Return to Contents
Next Section -- Results

Accessibility FOIA Privacy Policies and Notices

Take Pride in America logo USA.gov logo U.S. Department of the Interior | U.S. Geological Survey
URL: http://www.npwrc.usgs.gov/projects/wetland/studies/nitro/methods.htm
Page Contact Information: Webmaster
Page Last Modified: Thursday, 31-Aug-2006 16:49:49 EDT
Reston, VA [vaww55]