Water-quality data for 06610000 Missouri River at Omaha, Nebr.

Description

Scatterplots of water-quality data for 06610000 Missouri River at Omaha, Nebr.

Simazine concentrations in the
                                          Missouri River at Omaha, Nebr.
Prometon concentrations in the
                                          Missouri River at Omaha, Nebr.
Cyanazine concentrations in the
                                          Missouri River at Omaha, Nebr.
Metolachlor concentrations in the
                                          Missouri River at Omaha, Nebr.
Alachlor concentrations in the
                                          Missouri River at Omaha, Nebr.
Metribuzin concentrations in the
                                            Missouri River at Omaha, Nebr.
Trifluralin concentrations in the
                                            Missouri River at Omaha, Nebr.
EPTC concentrations in the
                                      Missouri River at Omaha, Nebr.
Trifluralin concentrations in the
                                            Missouri River at Omaha, Nebr.
Boxplot of data using regression on order statistics (Helsel, 2005; Lee, 2010) to estimate censored values. Boxplots of chemical concentrations 
                                          in the Missouri River at Omaha, Nebr.

Usage

1

Format

A data frame containing 115 water-quality samples for eight chemical constituents. There are 20 variables.

staid character USGS Station identification number
dates date Date water-quality sample collected
times numeric Time sample was collected
R04035 character Remark code (blank, _, <, or E)
P04035 numeric Simazine, water, filtered, recoverable, micrograms per liter
R04037 character Remark code (blank, _, <, or E)
P04037 numeric Prometon, water, filtered, recoverable, micrograms per liter
R04041 character Remark code (blank, _, <, or E)
P04041 numeric Cyanazine, water, filtered, recoverable, micrograms per liter
R39415 character Remark code (blank, _, <, or E)
P39415 numeric Metolachlor, water, filtered, recoverable, micrograms per liter
R46342 character Remark code (blank, _, <, or E)
P46342 numeric Alachlor, water, filtered, recoverable, micrograms per liter
R82630 character Remark code (blank, _, <, or E)
P82630 numeric Metribuzin, water, filtered, recoverable, micrograms per liter
R82661 character Remark code (blank, _, <, or E)
P82661 numeric Trifluralin, water, filtered (0.7 micron glass fiber filter), recoverable, micrograms per liter
R82668 character Remark code (blank, _, <, or E)
P82668 numeric EPTC, water, filtered (0.7 micron glass fiber filter), recoverable, micrograms per liter

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; <, indicating a censored value that is less than the number reported in the associated concentration field; and E, indicating that the value has been estimated. See Oblinger Childress and others (1999) for information on the remark codes used by the U.S. Geological Survey.

Source

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

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

  1
  2
  3
  4
  5
  6
  7
  8
  9
 10
 11
 12
 13
 14
 15
 16
 17
 18
 19
 20
 21
 22
 23
 24
 25
 26
 27
 28
 29
 30
 31
 32
 33
 34
 35
 36
 37
 38
 39
 40
 41
 42
 43
 44
 45
 46
 47
 48
 49
 50
 51
 52
 53
 54
 55
 56
 57
 58
 59
 60
 61
 62
 63
 64
 65
 66
 67
 68
 69
 70
 71
 72
 73
 74
 75
 76
 77
 78
 79
 80
 81
 82
 83
 84
 85
 86
 87
 88
 89
 90
 91
 92
 93
 94
 95
 96
 97
 98
 99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
  data(swData)
  
  # summary of water-quality concentrations
  apply(qwMoRivOmaha[,grep("P[[:digit:]]",dimnames(qwMoRivOmaha)[[2]])], 
  2, summary)
  
  # simple boxplot of water-quality concentrations
  rosBoxPlot(qwMoRivOmaha, qwcols=c("R", "P"))
  
  # same boxplot function with many additional plotting arguments
  rosBoxPlot(qwMoRivOmaha, site="06610000 Missouri River at Omaha, Nebr.", 
             log="y", yaxt="n", ylim=c(0.000001, 10), 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.000001, 0.00001, 0.0001, 0.001, 0.01, 0.1, 1, 10), 
       labels=c("0.000001", "0.00001", "0.0001", "0.001", "0.01", 
                "0.1", "1", "10"), cex.axis=0.7)
  
  # scatter plot of Simazine concentrations
  cenScatPlot(qwMoRivOmaha, pname="04035")
  
  # scatter plot with many additional plotting arguments
  par(las=1, tcl=0.5)
  cenScatPlot(qwMoRivOmaha, pname="04035", 
              site="06610000 Missouri River at Omaha, Nebr.",
              ylabel="Simazine concentration, in micrograms per liter", 
              legcex=0.7, qwcols=c("R", "P"),
              ylim=c(0,0.1), yaxs="i", cex.lab=0.9, cex.axis=0.9,
              xlim=c(as.Date("1996-01-01"), as.Date("2004-01-01")), 
              xaxs="i", xaxt="n")
  axdates<-c("1996-01-01", "1998-01-01", "2000-01-01", "2002-01-01", 
            "2004-01-01")
  axis(1, as.Date(axdates), labels=c("1996", "1998", "2000", "2002", 
                                     "2004"), cex.axis=0.9)
  
  # Prometon scatter plot
  cenScatPlot(qwMoRivOmaha, pname="04037", 
              site="06610000 Missouri River at Omaha, Nebr.",
              ylabel="Prometon concentration, in micrograms per liter", 
              legcex=0.7, qwcols=c("R", "P"),
              ylim=c(0,0.1), yaxs="i", cex.lab=0.9, cex.axis=0.9,
              xlim=c(as.Date("1996-01-01"), 
                     as.Date("2004-01-01")), xaxs="i", 
              xaxt="n")
  axdates<-c("1996-01-01", "1998-01-01", "2000-01-01", "2002-01-01", 
             "2004-01-01")
  axis(1, as.Date(axdates), labels=c("1996", "1998", "2000", "2002", 
                                     "2004"), cex.axis=0.9)

  # Cyanazine scatter plot
  cenScatPlot(qwMoRivOmaha, pname="04041", 
              site="06610000 Missouri River at Omaha, Nebr.",
              ylabel="Cyanazine concentration, in micrograms per liter",
              legcex=0.7, qwcols=c("R", "P"), 
              ylim=c(0.001,5), yaxs="i", cex.lab=0.9, cex.axis=0.9,
              xlim=c(as.Date("1996-01-01"), 
                     as.Date("2004-01-01")), xaxs="i", 
              xaxt="n", log="y")
  axdates<-c("1996-01-01", "1998-01-01", "2000-01-01", "2002-01-01", 
             "2004-01-01")
  axis(1, as.Date(axdates), labels=c("1996", "1998", "2000", "2002", 
                                     "2004"), cex.axis=0.9)

  # Metolachlor scatter plot
  cenScatPlot(qwMoRivOmaha, pname="39415", 
              site="06610000 Missouri River at Omaha, Nebr.",
              ylabel="Metolachlor concentration, in micrograms per liter",
              legcex=0.7, qwcols=c("R", "P"),
              ylim=c(0.001,5), yaxs="i", cex.lab=0.9, cex.axis=0.9,
              xlim=c(as.Date("1996-01-01"), 
                     as.Date("2004-01-01")), xaxs="i", 
              xaxt="n", log="y", legpos="bottomleft")
  axdates<-c("1996-01-01", "1998-01-01", "2000-01-01", "2002-01-01", 
             "2004-01-01")
  axis(1, as.Date(axdates), labels=c("1996", "1998", "2000", "2002", 
                                     "2004"), cex.axis=0.9)

  # Alachlor scatter plot
  cenScatPlot(qwMoRivOmaha, pname="46342", 
              site="06610000 Missouri River at Omaha, Nebr.",
              ylabel="Alachlor concentration, in micrograms per liter", 
              legcex=0.7, qwcols=c("R", "P"),
              ylim=c(0, 0.1), yaxs="i", cex.lab=0.9, cex.axis=0.9,
              xlim=c(as.Date("1996-01-01"), 
                     as.Date("2004-01-01")), xaxs="i", 
              xaxt="n")
  axdates<-c("1996-01-01", "1998-01-01", "2000-01-01", "2002-01-01", 
             "2004-01-01")
  axis(1, as.Date(axdates), labels=c("1996", "1998", "2000", "2002", 
                                     "2004"), cex.axis=0.9)
  
  # Metribuzin scatter plot
  cenScatPlot(qwMoRivOmaha, pname="82630", 
              site="06610000 Missouri River at Omaha, Nebr.",
              ylabel="Metribuzin concentration, in micrograms per liter",
              legcex=0.7, qwcols=c("R", "P"),
              ylim=c(0, 0.1), yaxs="i", cex.lab=0.9, cex.axis=0.9,
              xlim=c(as.Date("1996-01-01"), 
                     as.Date("2004-01-01")), xaxs="i", 
              xaxt="n")
  axdates<-c("1996-01-01", "1998-01-01", "2000-01-01", "2002-01-01", 
             "2004-01-01")
  axis(1, as.Date(axdates), labels=c("1996", "1998", "2000", "2002", 
                                     "2004"), cex.axis=0.9)

  # Trifluralin scatter plot
  cenScatPlot(qwMoRivOmaha, pname="82661", 
              site="06610000 Missouri River at Omaha, Nebr.",
              ylabel="Trifluralin concentration, in micrograms per liter",
              legcex=0.7, qwcols=c("R", "P"),
              ylim=c(0,0.03), yaxs="i", cex.lab=0.9, cex.axis=0.9,
              xlim=c(as.Date("1996-01-01"), 
                     as.Date("2004-01-01")), xaxs="i", 
              xaxt="n")
  axdates<-c("1996-01-01", "1998-01-01", "2000-01-01", "2002-01-01", 
             "2004-01-01")
  axis(1, as.Date(axdates), labels=c("1996", "1998", "2000", "2002", 
                                    "2004"), cex.axis=0.9)

  # EPTC scatter plot
  cenScatPlot(qwMoRivOmaha, pname="82668", 
              site="06610000 Missouri River at Omaha, Nebr.",
              ylabel="EPTC concentration, in micrograms per liter", 
              legcex=0.7, qwcols=c("R", "P"),
              ylim=c(0.001,1), yaxs="i", cex.lab=0.9, cex.axis=0.9,
              xlim=c(as.Date("1996-01-01"), 
                     as.Date("2004-01-01")), xaxs="i", 
              xaxt="n", log="y")
  axdates<-c("1996-01-01", "1998-01-01", "2000-01-01", "2002-01-01", 
             "2004-01-01")
  axis(1, as.Date(axdates), labels=c("1996", "1998", "2000", "2002", 
                                     "2004"), cex.axis=0.9)

Want to suggest features or report bugs for rdrr.io? Use the GitHub issue tracker.