I am using the dataRetrieval package to download mean daily streamflow data from USGS:
data_sc_stream_flow <- readNWISdv(stateCd="AZ",service="d