Irrigated agriculture has expanded greatly in the water-rich U.S. northern lake states during the past half century. Source water there is usually obtained from glacial aquifers strongly connected to surface waters, so irrigation has a potential to locally decrease base flows in streams and water levels in aquifers, lakes, and wetlands. During the nascent phase of the irrigation expansion, water availability was explored in works of some fame in the Wisconsin central sands by Weeks et al. (1965) on the Little Plover River and Weeks and Stangland (1971) on "headwater area" streams and lakes. Four decades later, and after irrigation has grown to a dominant landscape presence, we revisited irrigation effects on central sands hydrology. Irrigation effects have been substantial, on average decreasing base flows by a third or more in many stream headwaters and diminishing water levels by more than a meter in places. This explains why some surface waters have become flow and stage impaired, sometimes to the point of drying, with attendant losses of aquatic ecosystems. Irrigation exerts its effects by increasing evapotranspiration by an estimated 45 to 142 mm/year compared with pre-irrigated land cover. We conclude that irrigation water availability in the northern lake states and other regions with strong groundwater-surface water connections is tied to concerns for surface water health, requiring a focus on managing the upper few meters of aquifers on which surface waters depend rather than the depletability of an aquifer.

Runoff models such as the Curve Number (CN) model are dependent upon land use and soil type within the watershed or contributing area. In regions with internal drainage (e.g. wetlands) watershed delineation methods that fill sinks can result in inaccurate contributing areas and estimations of runoff from models such as the CN model. Two methods to account for this inaccuracy have been 1) to adjust the initial abstraction value within the CN model; or 2) to improve the watershed delineation in order to better account for internal drainage. We used a combined approach of examining the watershed delineation, and refining the CN model by the incorporating of dual hydrologic soil classifications. For eighteen watersheds within Wisconsin, we compared the CN model results of three watershed delineation methods to USGS gaged values. We found that for large precipitation events (>100 mm) the CN model estimations are closer to observed values for watershed delineations that identify the directly connected watershed and use the undrained hydrologic soil classification.

Digital Elevation Models (DEMs) are spatial grids which are used to automate watershed boundary determination. Sinks are present within most DEMs. In order to easily process the watershed boundary, the sinks are reassigned to elevation equivalent to an adjacent cell. The derived DEM is called a "filled" DEM. Due to its relative simplicity, the use of the "filled" DEM is one of the most widely used methods to delineate watershed boundaries and works well in about 70 percent of the watersheds in the US. In landscapes with internal drainage, sinks may accurately represent these depressions. In this study, we compare two delineation methods that do not fill in sinks to another method that does fill in sinks. We examined ten gaged watersheds in Wisconsin and Minnesota. We found the spatial extent of the watersheds from the three methods were significantly different. To evaluate the delineation methods, we modeled ten runoff events using the Curve Number (CN) method and compared them to USGS gage discharge for each watershed. For small storms we found that there were no significant differences in the modeled runoff for three delineation methods. For large storms, we found the no-fill methods had a smaller error, but overall the difference was insignificant. This research suggests that capturing internal drainage by the delineation does not have much of an impact on the widely used CN model.

Haucke, Jessica and Katherine A. Clancy, 2011. Stationarity of Streamflow Records and Their Influence on Bankfull Regional Curves. Journal of the American Water Resources Association (JAWRA) 47(6):1338–1347. DOI: 10.1111/j.1752‐1688.2011.00590.x Abstract:  Bankfull regional curves, which are curves that establish relationships among channel morphology, discharge, drainage area, are used extensively for stream restoration. These curves are developed upon the assumption that streamflows maintain stationarity over the entire record. We examined this assumption in the Driftless Area of southwestern Wisconsin where agricultural soil retention practices have changed, and precipitation has increased since the 1970s. We developed a bankfull regional curve for this area using field surveys of bankfull channel performed during 2008‐2009 and annual series of peak streamflows for 10 rivers with streamflow records ranging from the 1930s to 2009. We found bankfull flows to correlate to a 1.1 return period. To evaluate gage data statistics, we used the sign test to compare our channel morphology to historic 1.5 return period discharge (Q1.5) for five time periods: 1959‐1972, 1973‐1992, 1993‐2008, 1999‐2008, and the 1959‐2008 period of record. Analysis of the historic gage data indicated that there has been a more than 30% decline in Q1.5 since 1959. Our research suggests that land conservation practices may have a larger impact on gaging station stationarity than annual precipitation changes do. Additionally, historic peak flow data from gages, which have records that span land conservation changes, may need to be truncated to represent current flow regimes.

Macholl, Jacob A., Katherine A. Clancy, and Paul M. McGinley, 2011. Using a GIS Model to Identify Internally Drained Areas and Runoff Contribution in a Glaciated Watershed. Journal of the American Water Resources Association (JAWRA) 47(1):114‐125. DOI: 10.1111/j.1752‐1688.2010.00495.x Abstract:  Glaciated watersheds are not easily delineated using geographic information systems’ elevation‐based algorithms, especially where stream networks are disconnected and there are large regions of internally drained areas. This paper presents the results of an analysis using the Potential Contributing Source Area (PCSA) model to identify potential contributing areas, defined as areas from which runoff is physically capable of reaching a drainage network. The investigation was conducted to define the potential contributing areas in a glaciated region of northwest Wisconsin. The curve number (CN) method was used to predict runoff volumes in the watershed. The streamflows of four tributaries were measured and the runoff portion of the hydrograph quantified to be compared with runoff estimates calculated using the potential contributing areas and the traditional catchment area. Runoff producing events occurred, but the use of area‐weighted CN values was unsuccessful in modeling runoff due to all precipitation depths during the study period falling below the initial abstraction. A distributed CN approach provided runoff estimates that were generally better using the potential contributing areas compared with using the traditional catchment area. The extent of the minimum contributing area, estimated for a range of precipitation events, was found to be substantially less than the potential contributing areas, suggesting that the PCSA model delimits the maximum boundary of potential contributing areas.