seawaveQ (version 1.0.0)

IllRivValleyCty: Water-quality data for 05586100 Illinois River at Valley City, Ill.

Description

Scatterplots of water-quality data for 05586100 Illinois River at Valley City, Ill. html{

05586100Simazine.png

{Simazine concentrations in the Illinois River at Valley City, Ill.} } html{

05586100Prometon.png

{Prometon concentrations in the Illinois River at Valley City, Ill.} } html{

05586100Metribuzin.png

{Metribuzin concentrations in the Illinois River at Valley City, Ill.} } html{

05586100EPTC.png

{EPTC concentrations in the Illinois River at Valley City, Ill.} } latex{

05586100Simazine.pdf

{Simazine concentrations in the Illinois River at Valley City, Ill.} } latex{

05586100Prometon.pdf

{Prometon concentrations in the Illinois River at Valley City, Ill.} } latex{

05586100Metribuzin.pdf

{Metribuzin concentrations in the Illinois River at Valley City, Ill.} } latex{

05586100EPTC.pdf

{EPTC concentrations in the Illinois River at Valley City, Ill.} } Boxplot of data using regression on order statistics (Helsel, 2005; Lee, 2010) to estimate censored values. html{

05586100BoxPlots.png

{Boxplots of chemical concentrations in the Illinois River at Valley City, Ill.} } latex{

05586100BoxPlots.pdf

{Boxplots of chemical concentrations in the Illinois River at Valley City, Ill.} }

Usage

IllRivValleyCty

format

A data frame containing 168 water-quality samples for 4 constituents. There are 20 variables. lll{ site character Site abbreviation for study staid character USGS Station identification number dates date Date water-quality sample collected yrc numeric Year moc numeric Month dac numeric Day jdayc numeric Julian day from first day of associated streamflow data used R04035 character Remark code (blank, _,

source

Data provided by Patrick Phillips, U.S. Geological Survey, New York Water Science Center.

Details

Chemical concentration data are in the columns that start with a P and are followed by a number. Qualification codes for the concentration data are in the columns that start with an R followed by the same numbers as the associated concentration data. For example, column P04035 indicates simazine data, 04035, being the U.S. Geological Survey parameter code for simazine. The qualification codes for the simazine concentrations are found in the column R04035, indicating a U.S. Geological Survey remark code. Remark codes include _ or nothing, indicating no qualification of the value in the associated concentration field;

References

Oblinger Childress, C.J., Foreman, W.T., Connor, B.F., and Maloney, T.J., 1999, New reporting procedures based on long-term method detection levels and some considerations for interpretations of water-quality data provided by the U.S. Geological Survey Open-File Report 99--193, 19 p. (Also available at http://water.usgs.gov/owq/OFR_99-193/index.html.)

Examples

data(swData)

# summary of water-quality concentrations
apply(IllRivValleyCty[,grep("P[[:digit:]]", 
dimnames(IllRivValleyCty)[[2]])], 2, summary)

# simple boxplot of water-quality concentrations
rosBoxPlot(IllRivValleyCty)

# same boxplot function with many additional plotting arguments
rosBoxPlot(IllRivValleyCty, 
           site="05586100 Illinois River at Valley City, Ill.", log="y", 
           yaxt="n", ylim=c(0.0000001, 1), qwcols=c("R", "P"), 
           ylab=c("Concentration, micrograms per liter"), col="skyblue1", 
           cex.axis=0.7, cex.sub=0.8, par(tcl=0.5, las=1, 
                                                        yaxs="i", 
                                                        mgp=c(3,0.5,0), 
                                                        mar=c(5,5,2,2),
                                                        cex.main=0.9))
axis(2, at=c(0.0000001, 0.000001, 0.00001, 0.0001, 0.001, 0.01, 0.1, 1), 
     labels=c("0.0000001", "0.000001", "0.00001", "0.0001", "0.001", "0.01", 
              "0.1", "1"), cex.axis=0.7)

# scatter plot of simazine concentrations
cenScatPlot(IllRivValleyCty, pname="04035")

# scatter plot with many additional plotting arguments
par(las=1, tcl=0.5)
cenScatPlot(IllRivValleyCty, pname="04035", 
            site="05586100 Illinois River at Valley City, Ill.",
            ylabel="Simazine concentration, in micrograms per liter", 
            legcex=0.7, 
            ylim=c(0,0.4), yaxs="i", cex.lab=0.9, cex.axis=0.9,
            xlim=c(as.Date("1996-01-01"), as.Date("2012-01-01")), 
            xaxs="i", xaxt="n")
axdates<-c("1996-01-01", "2000-01-01", "2004-01-01", "2008-01-01",
           "2012-01-01")
axis(1, as.Date(axdates), labels=c("1996", "2000", "2004", "2008",
                                   "2012"), cex.axis=0.9)

# Prometon scatter plot
cenScatPlot(IllRivValleyCty, pname="04037", 
            site="05586100 Illinois River at Valley City, Ill.",
            ylabel="Prometon concentration, in micrograms per liter", 
            legcex=0.7, 
            ylim=c(0,0.06), yaxs="i", cex.lab=0.9, cex.axis=0.9,
            xlim=c(as.Date("1996-01-01"), 
                   as.Date("2012-01-01")), xaxs="i", 
            xaxt="n")
axdates<-c("1996-01-01", "2000-01-01", "2004-01-01", "2008-01-01", 
           "2012-01-01")
axis(1, as.Date(axdates), labels=c("1996", "2000", "2004", "2008",
                                               "2012"), cex.axis=0.9)

# Metribuzin scatter plot
cenScatPlot(IllRivValleyCty, pname="82630", 
            site="05586100 Illinois River at Valley City, Ill.",
            ylabel="Metribuzin concentration, in micrograms per liter", 
            legcex=0.7, 
            ylim=c(0,0.3), yaxs="i", cex.lab=0.9, cex.axis=0.9,
            xlim=c(as.Date("1996-01-01"), 
                   as.Date("2012-01-01")), xaxs="i", 
            xaxt="n")
axdates<-c("1996-01-01", "2000-01-01", "2004-01-01", "2008-01-01", 
           "2012-01-01")
axis(1, as.Date(axdates), labels=c("1996", "2000", "2004", "2008",
                                               "2012"), cex.axis=0.9)
# EPTC scatter plot
cenScatPlot(IllRivValleyCty, pname="82668", 
            site="05586100 Illinois River at Valley City, Ill.",
            ylabel="EPTC concentration, in micrograms per liter", 
            legcex=0.7, ylim=c(0,0.08), yaxs="i", cex.lab=0.9, 
            cex.axis=0.9, xlim=c(as.Date("1996-01-01"), 
                   as.Date("2012-01-01")), xaxs="i", xaxt="n")
axdates<-c("1996-01-01", "2000-01-01", "2004-01-01", "2008-01-01", 
           "2012-01-01")
axis(1, as.Date(axdates), labels=c("1996", "2000", "2004", "2008","2012"), 
     cex.axis=0.9)