diff --git a/R/readNWISunit.r b/R/readNWISunit.r
index 595ae09d0ba8e13d1c274b886ccf920e9ee16ace..5f27695e44837fa318f641d2c3bbcf0b533f6d20 100644
--- a/R/readNWISunit.r
+++ b/R/readNWISunit.r
@@ -69,7 +69,7 @@ readNWISuv <- function (siteNumbers,parameterCd,startDate="",endDate="", tz=""){
 #' Reads peak flow from NWISweb. 
 #' Data is retrieved from \url{http://waterdata.usgs.gov/nwis}. 
 #' 
-#' @param siteNumber character USGS site number.  This is usually an 8 digit number
+#' @param siteNumbers character USGS site number(or multiple sites).  This is usually an 8 digit number.
 #' @param startDate character starting date for data retrieval in the form YYYY-MM-DD. Default is "" which indicates
 #' retrieval for the earliest possible record.
 #' @param endDate character ending date for data retrieval in the form YYYY-MM-DD. Default is "" which indicates
@@ -102,9 +102,9 @@ readNWISuv <- function (siteNumbers,parameterCd,startDate="",endDate="", tz=""){
 #' }
 #' @export
 #' @examples
-#' siteNumber <- '01594440'
-#' data <- readNWISpeak(siteNumber)
-readNWISpeak <- function (siteNumber,startDate="",endDate=""){  
+#' siteNumbers <- c('01594440','040851325')
+#' data <- readNWISpeak(siteNumbers)
+readNWISpeak <- function (siteNumbers,startDate="",endDate=""){  
   
   # Doesn't seem to be a peak xml service
   url <- constructNWISURL(siteNumber,NA,startDate,endDate,"peak")
@@ -186,7 +186,7 @@ readNWISrating <- function (siteNumber,type="base"){
 #'Reads surface-water measurement data from NWISweb. Data is retrieved from \url{http://waterdata.usgs.gov/nwis}.
 #'See \url{http://waterdata.usgs.gov/usa/nwis/sw} for details about surface water.
 #'
-#' @param siteNumber character USGS site number.  This is usually an 8 digit number
+#' @param siteNumbers character USGS site number (or multiple sites).  This is usually an 8 digit number
 #' @param startDate character starting date for data retrieval in the form YYYY-MM-DD. Default is "" which indicates
 #' retrieval for the earliest possible record.
 #' @param endDate character ending date for data retrieval in the form YYYY-MM-DD. Default is "" which indicates
@@ -217,12 +217,12 @@ readNWISrating <- function (siteNumber,type="base"){
 #' }
 #' @export
 #' @examples
-#' siteNumber <- '01594440'
-#' data <- readNWISmeas(siteNumber)
-readNWISmeas <- function (siteNumber,startDate="",endDate="", tz=""){  
+#' siteNumbers <- c('01594440','040851325')
+#' data <- readNWISmeas(siteNumbers)
+readNWISmeas <- function (siteNumbers,startDate="",endDate="", tz=""){  
   
   # Doesn't seem to be a WaterML1 format option
-  url <- constructNWISURL(siteNumber,NA,startDate,endDate,"meas")
+  url <- constructNWISURL(siteNumbers,NA,startDate,endDate,"meas")
   
   data <- importRDB1(url,asDateTime=FALSE,tz=tz)
   
@@ -230,7 +230,7 @@ readNWISmeas <- function (siteNumber,startDate="",endDate="", tz=""){
     data$diff_from_rating_pc <- as.numeric(data$diff_from_rating_pc)
   }
   
-  siteInfo <- readNWISsite(siteNumber)
+  siteInfo <- readNWISsite(siteNumbers)
   
   attr(data, "siteInfo") <- siteInfo
   attr(data, "variableInfo") <- NULL
diff --git a/man/readNWISmeas.Rd b/man/readNWISmeas.Rd
index 7cb0b8b1fd6f4c3f5bb87d657800f739cafb3e12..800c18dbd511cfb0b2d69c77fcb8efc1fb031e5c 100644
--- a/man/readNWISmeas.Rd
+++ b/man/readNWISmeas.Rd
@@ -3,10 +3,10 @@
 \alias{readNWISmeas}
 \title{Reads surface-water measurement data from NWISweb.}
 \usage{
-readNWISmeas(siteNumber, startDate = "", endDate = "", tz = "")
+readNWISmeas(siteNumbers, startDate = "", endDate = "", tz = "")
 }
 \arguments{
-\item{siteNumber}{character USGS site number.  This is usually an 8 digit number}
+\item{siteNumbers}{character USGS site number (or multiple sites).  This is usually an 8 digit number}
 
 \item{startDate}{character starting date for data retrieval in the form YYYY-MM-DD. Default is "" which indicates
 retrieval for the earliest possible record.}
@@ -46,7 +46,7 @@ Reads surface-water measurement data from NWISweb. Data is retrieved from \url{h
 See \url{http://waterdata.usgs.gov/usa/nwis/sw} for details about surface water.
 }
 \examples{
-siteNumber <- '01594440'
-data <- readNWISmeas(siteNumber)
+siteNumbers <- c('01594440','040851325')
+data <- readNWISmeas(siteNumbers)
 }
 
diff --git a/man/readNWISpeak.Rd b/man/readNWISpeak.Rd
index de3d4d82058caba21435f54e5d5aafee6eacf7e4..59a601e8c017062c62c62033ccd7e8d89867e5f3 100644
--- a/man/readNWISpeak.Rd
+++ b/man/readNWISpeak.Rd
@@ -3,10 +3,10 @@
 \alias{readNWISpeak}
 \title{Reads peak flow data from NWISweb.}
 \usage{
-readNWISpeak(siteNumber, startDate = "", endDate = "")
+readNWISpeak(siteNumbers, startDate = "", endDate = "")
 }
 \arguments{
-\item{siteNumber}{character USGS site number.  This is usually an 8 digit number}
+\item{siteNumbers}{character USGS site number(or multiple sites).  This is usually an 8 digit number.}
 
 \item{startDate}{character starting date for data retrieval in the form YYYY-MM-DD. Default is "" which indicates
 retrieval for the earliest possible record.}
@@ -47,7 +47,7 @@ Reads peak flow from NWISweb.
 Data is retrieved from \url{http://waterdata.usgs.gov/nwis}.
 }
 \examples{
-siteNumber <- '01594440'
-data <- readNWISpeak(siteNumber)
+siteNumbers <- c('01594440','040851325')
+data <- readNWISpeak(siteNumbers)
 }