Geology Reports⌕ Search

SEARCH · Geology Reports

Results for “Ground Water Modeling Investigation Report”

Search indexed USGS publications on groundwater, aquifers, geologic maps, mineral resources and earthquakes. Explore source records by subject and place.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 811 records · Page 45Linked to original sources
USGS1997

Review of selected features of the natural system model, and suggestions for applications in south Florida

A study was conducted to review selected features of the Natural System Model, version 4.3. The Natural System Model is a regional-scale model that uses recent climatic data and estimates of historic vegetation and topography to simulate pre-canal-drainage hydrologic response in south Florida. Equations used to represent the hydrologic system and the numerical solution of these equations in the model were documented and reviewed. Convergence testing was performed using 1965 input data, and selected other aspects of the model were evaluated. Some conclusions from the evaluation of the Natural System Model include the following observations. Simulations were generally insensitive to the temporal resolution used in the model. However, reduction of the computational cell size from 2-mile by 2-mile to 2/3-mile by 2/3-mile resulted in a decrease in spatial mean ponding depths for October of 0.35 foot for a 3-hour time step. Review of the computer code indicated that there is no limit on the amount of water that can be transferred from the river system to the overland flow system, on the amount of seepage from the river to the ground-water system, on evaporation from the river system, or on evapotranspiration from the overland-flow system. Oscillations of 0.2 foot or less in simulated river stage were identified and attributed to a volume limiting function which is applied in solution of the overland-flow equations. The computation of the resistance coefficient is not consistent with the computation of overland-flow velocity. Ground-water boundary conditions do not always ensure a no-flow condition at the boundary. These inconsistencies had varying degrees of effects on model simulations, and it is likely that simulations longer than 1 year are needed to fully identify effects. However, inconsistencies in model formulations should not be ignored, even if the effects of such errors on model results appear to be small or have not been clearly defined. The Natural System Model can be a very useful tool for estimating pre-drainage hydrologic response in south Florida. The model includes all of the important physical processes needed to simulate a water balance. With a few exceptions, these hydrologic processes are represented in a reasonable manner using empirical, semiempirical, and mechanistic relations. The data sets that have been assembled to represent physical features, and hydrologic and meteorological conditions are quite extensive in their scope. Some suggestions for model application were made. Simulation results from the Natural System Model need to be interpreted on a regional basis, rather than cell by cell. The available evidence suggests that simulated water levels should be interpreted with about a plus or minus 1 foot uncertainty. It is probably not appropriate to use the Natural System Model to estimate pre-drainage discharges (as opposed to hydroperiods and water levels) at a particular location or across a set of adjacent computational cells. All simulated results for computational cells within about 10 miles of the model boundaries have a higher degree of uncertainty than results for the interior of the model domain. It is most appropriate to interpret the Natural System Model simulation results in connection with other available information. Stronger linkages between hydrologic inputs to the Everglades and the ecological response of the system would enhance restoration efforts.

Florida↗
USGS1994

Hydrogeologic framework and simulation of shallow ground-water flow in the vicinity of a hazardous-waste landfill near Pinewood, South Carolina

The geologic units in the vicinity of a hazardous- waste landfill near Pinewood, S.C., were divided into hydrogeologic units on the basis of lithologic and hydrologic characteristics. A quasi-3- dimensional finite-difference ground-water-flow model was constructed to represent the hydrogeologic flamework. The simulation results indicate that if non-reactive constituents were released to the Lang Syne water-bearing zone underlying the central and western pans of the disposal area, the constituents would move in a southwesterly direction at a rate of about one-half to 7 feet per year. Contaminants could move from the Lang Syne water-bearing zone upward to the surficial aquifer, to streams, or to Lake Marion. Although these flow rates indicate that it would require at least 50 years for contaminants to travel between the disposal area and a nearby (400 ft) potential discharge area, the heterogeneity of the site hydrogeology imparts an uncertainty to the conclusion. Faster travel times cannot be ruled out if contamination enters an area having a higher hydraulic conductivity than those determined in this investigation. Faster arrival times at Lake Marion also could result if there are pathways shorter than about 400 feet between contaminated water and an area where it can discharge to the surficial aquifer or to streams. If contaminant releases were to occur on the eastern side of the ground-water mounds, near landfill section II and the southeastern part of land fill section I, initial flow directions would be toward the water-level depression in the eastern part of the facility. Ground water within water- level depression would flow downward, probably to the underlying lower Sawdust Landing water-beating zone. Movement of non-reactive constituents in the tower Sawdust Landing water-bearing zone would be southwestward toward Lake Marion at a rate of about 8 to 20 feet per year. Transport to the lake by this route could require more than 200 years.

South Carolina↗
USGS1981

Saline ground-water discharge to the Smoky Hill River between Salina and Abilene, central Kansas

Saline water discharges from the alluvial aquifer into the Smoky Hill and Solomon Rivers between New Cambria and Sand Springs, Kansas. During relatively stable base flow in 1976-77, the discharge was about 32 cubic feet per second. Chloride concentrations at base flow increased about 800 milligrams per liter in the Smoky Hill River and 550 milligrams per liter in the Solomon River. The source of the saline water is the underlying Wellington aquifer, a zone of dissolution, subsidence, and collapse along the eastern margin of the Wellington Formation. Locally, brine from the aquifer moves upward through collapse structures in the confining layer at the base of the alluvium. The brine discharge ranges from 0.3 to 0.8 cubic foot per second, and the chloride load ranges from 150 to 370 tons per day. Results from a mathematical model of the fl ow system indicated that recharge from periodic flooding, as in 1973, was sufficient to reverse the normal (197677) hydraulic gradient between aquifers. Although brine discharge was temporarily reduced, saline-water discharge to the rivers increased. The discharge of brine could be intercepted by pumping wells completed in the Wellington aquifer to lower the potentiometric head and reduce upward flow. The intercepted brine could be discharged into formations underlying the area at depth or piped to a storage reservoir from which the brine would be evaporated or released to the rivers during peak flows. Also, the freshwater in upstream base flow could be diverted by canal, and the saline-water discharge could be retained by low-head dams in the river channel for release during peak flows.

Kansas↗
USGS2004

Simulation of ground-water flow, contributing recharge areas, and ground-water travel time in the Missouri River alluvial aquifer near Ft. Leavenworth, Kansas

The Missouri River alluvial aquifer near Ft. Leavenworth, Kansas, supplies all or part of the drinking water for Ft. Leavenworth; Leavenworth, Kansas; Weston, Missouri; and cooling water for the Kansas City Power and Light, Iatan Power Plant. Ground water at three sites within the alluvial aquifer near the Ft. Leavenworth well field is contaminated with trace metals and organic compounds and concerns have been raised about the potential contamination of drinking-water supplies. In 2001, the U.S. Geological Survey, U.S. Army Corps of Engineers, and the U.S. Army began a study of ground-water flow in the Missouri River alluvial aquifer near Ft. Leavenworth. Hydrogeologic data from 173 locations in the study area was used to construct a ground-water flow model (MODFLOW-2000) and particle-tracking program (MODPATH) to determine the direction and travel time of ground-water flow and contributing recharge areas for water-supply well fields within the alluvial aquifer. The modeled area is 28.6 kilometers by 32.6 kilometers and contains the entire study area. The model uses a uniform grid size of 100 meters by 100 meters and contains 372,944 cells in 4 layers, 286 columns, and 326 rows. The model represents the alluvial aquifer using four layers of variable thickness with no intervening confining layers. The model was calibrated to both quasi-steady-state and transient hydraulic head data collected during the study and ground-water flow was simulated for five well-pumping/river-stage scenarios. The model accuracy was calculated using the root mean square error between actual measurements of hydraulic head and model generated hydraulic head at the end of each model run. The accepted error for the model calibrations were below the maximum measurement errors. The error for the quasi-steady-state calibration was 0.82 meter; for the transient calibration it was 0.33 meter. The shape, size, and ground-water travel time within the contributing recharge area for each well or well field is affected by changes in river stage and pumping rates and by the location of the well or well field with respect to the major rivers, alluvial valley walls, and other pumping wells. The shapes of the simulated contributing recharge areas for the well fields in the study area are elongated in the upstream direction for all well-pumping/river-stage scenarios. The capture of ground water by the pumping wells as it moved downgradient toward the Missouri River caused the long up-valley extent of the contributing recharge areas. Recharge to the Iatan and Weston well fields primarily is from precipitation and surface runoff from the surrounding uplands because the contributing recharge area does not intersect the Missouri River for any well-pumping/river-stage scenarios. Recharge to the Leavenworth and Ft. Leavenworth well fields is from precipitation, surface runoff from the surrounding uplands, and the Missouri River because the contributing recharge area intersects these boundaries for all well-pumping/river-stage scenarios. Particle tracking analysis indicated ground water from the three contaminated sites was captured by the Ft. Leavenworth well field for all well-pumping/river-stage scenarios. Ground-water travel times to the Ft. Leavenworth well field for average well-pumping/river-stage scenario ranged from about 33 years for the closest contamination site to about 71 years for the farthest contamination site. Ground-water flow was induced below the Missouri River by the Ft. Leavenworth and Leavenworth well fields for all well-pumping/river-stage scenarios.

Scientific Investigations Report↗
USGS2015

The Everglades Depth Estimation Network (EDEN) surface-water model, version 2

The Everglades Depth Estimation Network (EDEN) is an integrated network of water-level gages, interpolation models that generate daily water-level and water-depth data, and applications that compute derived hydrologic data across the freshwater part of the greater Everglades landscape. The U.S. Geological Survey Greater Everglades Priority Ecosystems Science provides support for EDEN in order for EDEN to provide quality-assured monitoring data for the U.S. Army Corps of Engineers Comprehensive Everglades Restoration Plan. The EDEN surface-water model, version 2 (V2), interpolates water-level data from a network of 240 gages to generate gridded daily water-level surfaces for the freshwater domain of the Everglades. When these spatiotemporal continuous surfaces are combined with EDEN’s digital elevation model of ground surface, derived hydrologic data provide scientists and water managers working in the Everglades with data necessary to analyze ecological and biotic responses to hydrologic changes in the Everglades. Derived datasets include water depth, recession rates, days since last dry, water-surface slopes, and hydroperiod. The V2 model includes enhancements from the previous model (version 1; V1) to accommodate changes in the water-level gage network, adjustments to water-level data, improved understanding of the flow dynamics (particularly near canals), and installation of an elevation benchmark network. Enhancements to the V2 model included Expansion of the EDEN domain: The model domain was expanded to include a part of southern Big Cypress National Preserve and northwestern Everglades National Park upstream of the marsh mangrove wetlands, thus completing the coastal connection along the southwestern boundary of the model; and Development of subdomain models: To account for insufficient water-control structure gage data at some subbasin boundaries, subdomain models were developed for five subdomains, and the resulting water-level surfaces were merged to generate the final water-level surface. Model performance statistics show a general improvement in the V2 model as compared to the V1 model. Overall, the root mean squared error (RMSE) was reduced by 2.42 centimeters (cm) to 4.68 cm. In Water Conservation Area 3A North and Water Conservation Area 3B, the RMSE was reduced by 10.88 and 9.15 cm, respectively. In addition to evaluating model performance statistics, 2-cm water-level maps were generated and evaluated for irregular contours that would indicate a potential problem either with data input or water-level estimates. Three applications of the EDEN-modeled water surfaces and other EDEN datasets are presented in the report to show how scientists and resource managers are using EDEN datasets to analyze biological and ecological responses to hydrologic changes in the Everglades. The biological responses of two important Everglades species, alligators and wading birds, to changes in hydrology are described. The effects of hydrology on fire dynamics in the Everglades are also discussed.

Florida↗
USGS1982

Two-dimensional, steady-state model of ground-water flow, Nevada Test Site and vicinity, Nevada-California

A two-dimensional, steady-state model of ground-water flow beneath the Nevada Test Site and vicinity has been developed using inverse techniques. The area is underlain by clastic and carbonate rocks of Precambrian and Paleozoic age and by volcanic rocks and alluvium of Tertiary and Quaternary age that have been juxtaposed by normal and strike-slip faulting. Aquifers are composed of carbonate and volcanic rocks and alluvium. Characteristics of the flow system are determined by distribution of low-conductivity rocks (barriers); by recharge originating in the Spring Mountains, Pahranagat, Timpahute, and Sheep Ranges, and in Pahute Mesa; and by underflow beneath Pahute Mesa from Gold Flat and Kawich Valley. Discharge areas (Ash Meadows, Oasis Valley, Alkali Flat, and Furnace Creek Ranch) are upgradient from barriers. Sensitivities of simulated hydraulic heads and fluxes to variations in model parameters were calculated to guide field studies and to help estimate errors in predictions from transport modeling. Hydraulic heads and fluxes are very sensitive to variations in the greater magnitude recharge/discharge terms. Transmissivity at a location may not be the most important transmissivity for determining flux there. Transmissivities and geometries of large barriers that impede flow from Pahute Mesa have major effects on fluxes elsewhere; as their transmissivities are decreased, flux beneath western Jackass Flats and Yucca Mountains is increased as water is diverted around the barriers. Fortymile Canyon is underlain by highly transmissive rocks that cause potentiometric contours to vee upgradient; increasing their transmissivity increases flow through them, and decreases it beneath Yucca Mountain.

California, Nevada↗
USGS2005

Pond-aquifer flow and water availability in the vicinity of two coastal area seepage ponds, Glynn and Bulloch Counties, Georgia

Pond-aquifer flow and water availability at excavated seepage pond sites in Glynn County and in southern Bulloch County, Georgia, were evaluated to determine their potential as sources of water supply for irrigation. Excavated seepage ponds derive water primarily from ground water seeping into the pond, in a manner similar to a dug well completed in a surficial aquifer. The availability of water from seepage ponds is controlled by the permeability of surficial deposits, the amount of precipitation recharging the ground-water system, and the volume of water stored in the pond. The viability of seepage ponds as supplies for irrigation is limited by low seepage rates and high dependence on climatic conditions. Ponds will not refill unless there is adequate precipitation to recharge the surficial aquifer, which subsequently drains (seeps) into the pond. Ground-water seepage was estimated using a water-budget approach that utilized on-site climatic and hydrologic measurements, computing pond-volume changes during pond pumping tests, and by digital simulation using steady-state and transient ground-water flow models. From August 1999 to May 2000, the Glynn County pond was mostly losing water (as indicated by negative net seepage); whereas from October 2000 to June 2001, the Bulloch County pond was mostly gaining water. At both sites, most ground-water seepage entered the pond following major rainfall events that provided recharge to the surficial aquifer. Net ground-water seepage, estimated using water-budget analysis and simulation, ranged from -11.5 to 15 gallons per minute (gal/min) at the Glynn County pond site and from -55 to 31 gal/min at the Bulloch County pond site. Simulated values during pumping tests indicate that groundwater seepage to both ponds increases with decreased pond stage. At the Glynn County pond, simulated net ground-water seepage varied between 7.8 gal/min at the beginning of the test (high pond stage and low hydraulic gradient) and 103 gal/min at the end of the test (low pond stage and high hydraulic gradient). At the Bulloch County pond site, values ranged from -17.7 gal/min at the beginning of the test to 15 gal/min at the end of the test. Results at the two pond sites indicate that the use of excavated seepage ponds as sources for irrigation supply is limited by pond-storage volume and low net ground-water seepage rates during periods of low precipitation. Pumps withdrawing 1,000 gal/min for 10 hours per dayunder climatic and hydrologic conditions similar to those observed during pond pumping tests at each sitewould drain the Glynn County pond within 30 days and the Bulloch County pond within 3.5 days. Because the two pond sites are considered to represent the extremes of likely conditions to be encountered in the coastal Georgia area, it is likely that other seepage ponds would have similar storage-depletion rates.

Georgia↗
USGS2011

Summary of November 2010 meeting to evaluate turbidite data for constraining the recurrence parameters of great Cascadia earthquakes for the update of national seismic hazard maps

This report summarizes a meeting of geologists, marine sedimentologists, geophysicists, and seismologists that was held on November 18–19, 2010 at Oregon State University in Corvallis, Oregon. The overall goal of the meeting was to evaluate observations of turbidite deposits to provide constraints on the recurrence time and rupture extent of great Cascadia subduction zone (CSZ) earthquakes for the next update of the U.S. national seismic hazard maps (NSHM). The meeting was convened at Oregon State University because this is the major center for collecting and evaluating turbidite evidence of great Cascadia earthquakes by Chris Goldfinger and his colleagues. We especially wanted the participants to see some of the numerous deep sea cores this group has collected that contain the turbidite deposits. Great earthquakes on the CSZ pose a major tsunami, ground-shaking, and ground-failure hazard to the Pacific Northwest. Figure 1 shows a map of the Pacific Northwest with a model for the rupture zone of a moment magnitude M w 9.0 earthquake on the CSZ and the ground shaking intensity (in ShakeMap format) expected from such an earthquake, based on empirical ground-motion prediction equations. The damaging effects of such an earthquake would occur over a wide swath of the Pacific Northwest and an accompanying tsunami would likely cause devastation along the Pacifc Northwest coast and possibly cause damage and loss of life in other areas of the Pacific. A magnitude 8 earthquake on the CSZ would cause damaging ground shaking and ground failure over a substantial area and could also generate a destructive tsunami. The recent tragic occurrence of the 2011 M w 9.0 Tohoku-Oki, Japan, earthquake highlights the importance of having accurate estimates of the recurrence times and magnitudes of great earthquakes on subduction zones. For the U.S. national seismic hazard maps, estimating the hazard from the Cascadia subduction zone has been based on coastal paleoseismic evidence of great earthquakes over the past 5,000 years. The instrumental catalog of earthquakes is of little use for constraining the hazard of the CSZ, because there are virtually no recorded earthquakes on most of the plate interface of the CSZ. There are no historical accounts in the past 150 years of large earthquakes on most of the CSZ. Until about 20 years ago, some interpreted this lack of recent and historical earthquakes as an indicator that the subduction zone was slipping aseismically and could not produce a great earthquake. The work of Brian Atwater and others, in the late 1980s and the 1990s (Atwater, 1987, 1992; Atwater and others, 1995; Nelson and others, 1996; Clague, 1997; Atwater and Hemphill-Haley, 1997; Atwater and others, 2004) demonstrated that submerged forests, buried soils, tsunami deposits, and liquefaction along and near the coast were compelling evidence of repeated great Cascadia earthquakes over at least the past 5,000 years. Atwater and Hemphill-Haley (1997) concluded from paleoseismic evidence at Willapa Bay, Washington, that great earthquakes ruptured the CSZ with an average recurrence time of about 500 years. The date of the last great CSZ earthquake, January 26, 1700, was established from historical records of the so-called orphan tsunami in Japan that is inferred to have been produced by this earthquake (Satake and others, 1996, 2003; Atwater and others, 2005) and is consistent with tree-ring data from drowned forests in Washington and Oregon. From modeling the observations of the tsunami, Satake and others (2003) estimated a moment magnitude of about 9.0 for this earthquake. Many other paleoseismic sites have been investigated along the Pacific Northwest coast from Vancouver Island to northern California and show evidence of great CSZ earthquakes. Nelson and others (2006) summarized the dates found from these studies and proposed correlations between sites indicating the extent of rupture for individual events. Dating of inferred tsunami deposits in Bradley Lake, Oregon by Kelsey and others (2005), as well as tsunami and subsidence evidence from Six Rivers, Oregon (Kelsey and others, 2002) and Coquille River (Witter and others, 2003), indicates that there were probably M w 8 ruptures in the southern portion of the CSZ in addition to the M w 9 events that rupture the whole length of the CSZ (Nelson and others, 2006). A parallel development over the past 20 years or more is the use of deep-sea turbidite deposits for identifying and dating great Cascadia earthquakes over the past 10,000 years (Adams, 1990; Goldfinger and others, 2003, 2008, in press; Goldfinger, 2011). Turbidites are sediment deposits in the deep ocean from turbidity currents, which are energetic flows of sediment and water along the continental shelf and slope. Adams (1990), using the counts of turbidites in deep-sea cores off the coast of Oregon and Washington collected and analyzed by Griggs (1969) and Griggs and others (1969), proposed that these turbidites were caused by the shaking of great Cascadia earthquakes. Part of his reasoning was that the number (13) of turbidite deposits that occurred since deposition of the Mazama Ash 7,000 years ago gave a recurrence time of about 500 years, consistent with that derived from the coastal submergence data. Adams (1990) also proposed the “confluence test” which evaluates the number of turbidites for submarine channels that form a confluence. He reported that the number of turbidites in the single downstream channel equaled the number in each of the tributary channels. He reasoned that this indicated that the turbidites in each tributary were simultaneously triggered and were, therefore, caused by a common forcing agent. He concluded that shaking from extended ruptures of great Cascadia earthquakes was the most likely cause of these turbidites. Based on the paleoseismic evidence of past great earthquakes, the hazard from the Cascadia subduction zone was included in the 1996 U.S. NSHM (Frankel and others, 1996), which were the basis for seismic provisions in the 2000 International Building Code. These hazard maps used the paleoseismic studies to constrain the recurrence rate of great CSZ earthquakes. Goldfinger and his colleagues have since collected many more deep ocean cores and done extensive analysis on the turbidite deposits that they identified in the cores (Goldfinger and others, 2003, 2008, in press; Goldfinger, 2011). Using their dating of the sediments and correlation of features in the logs of density and magnetic susceptibility between cores, they developed a detailed chronology of great earthquakes along the CSZ for the past 10,000 years (Goldfinger and others, in press). These correlations consist of attempting to match the peaks and valleys in logs of density and magnetic susceptibility between cores separated, in some cases, by hundreds of kilometers. Based on this work, Goldfinger and others (2003, 2008, in press) proposed that the turbidite evidence indicated the occurrence of great earthquakes (M w 8) that only ruptured the southern portion of the CSZ, as well as earthquakes with about M w 9 that ruptured the entire length of the CSZ. For the southernmost portion of the CSZ, Goldfinger and others (in press) proposed a recurrence time of M w 8 or larger earthquakes of about 230 years. This proposed recurrence time was shorter than the 500 year time that was incorporated in one scenario in the NSHM’s. It is important to note that the hazard maps of 1996 and later also included a scenario or set of scenarios with a shorter recurrence time for M w 8 earthquakes, using rupture zones that are distributed along the length of the CSZ (Frankel and others, 1996; Petersen and others, 2008). Originally, this scenario was meant to correspond to the idea that some of the 500-year averaged ruptures seen in the paleoseismic evidence could have been a series of M w 8 earthquakes that occurred over a short period of time (a few decades), rather than M w 9 earthquakes. Figure 2 shows the logic tree for the CSZ used in the 2008 NSHM’s (Petersen and others, 2008). This logic tree includes whole CSZ rupture earthquakes (M w 8.8–9.2) and partial CSZ rupture earthquakes (M w 8.0–8.7). In this latest version of the NSHM’s, the effective recurrence time of earthquakes on the CSZ with moment magnitudes greater than or equal to 8.0 over the various models is about 270 years (Petersen and others, 2008). This recurrence time applies to the entire CSZ, so that the hazard from great earthquakes was approximately equal along the whole zone, although the hazard estimates taper on the northern and southern ends of the CSZ, because of the way rupture zones of M w 8 earthquakes were distributed along the strike of the CSZ. The NSHM will be updated in 2013, as part of the standard update cycle that corresponds to the update cycle of the national model building codes that are based on the seismic hazard maps. A meeting was necessary to assemble a wide group of experts to hear Dr. Goldfinger explain his methodology for dating and correlating the turbidites and for developing the earthquake chronology. The overall goal of the workshop was to evaluate observations of turbidite deposits to provide constraints on the recurrence times and rupture extents of great Cascadia subduction zone earthquakes for the next update of the NSHM. Before the meeting, participants were supplied with the U.S. Geological Survey (USGS) Professional Paper of Goldfinger and others (in press), as well as material from Brian Atwater and Alan Nelson. The agenda of the meeting was developed by Art Frankel, with assistance from Chris Goldfinger, Brian Atwater, Alan Nelson, Mark Petersen, and Craig Weaver. The meeting was hosted by Chris Goldfinger of Oregon State University. We stress that it is difficult to evaluate in a two-day meeting the large amount of work that Goldfinger and his colleagues have done over the past 15 years or more. This meeting is the first step in a process that develops the inputs to the update of the national maps. The conclusions of this workshop will be discussed and possibly modified at the regional Pacific Northwest workshop for the hazard maps to be held in early 2012. Vetting new research results using informed expert opinion is an integral part of updating the national maps and does not reflect on the veracity of these results.

Cascadia↗
USGS1998

Estimate of aquifer properties by numerically simulating ground-water/surface-water interactions, Fort Wainwright, Alaska

MODFLOW, a finite-difference model of ground-water flow, was used to simulate the flow of water between the aquifer and the Chena River at Fort Wainwright, Alaska. The model was calibrated by comparing simulated ground-water hydrographs to those recorded in wells during periods of fluctuating river levels. The best fit between simulated and observed hydrographs occurred for the following: 20 feet per day for vertical hydraulic conductivity, 400 feet per day for horizontal hydraulic conductivity, 1:20 for anisotropy (vertical to horizontal hydraulic conductivity), and 350 per feet for riverbed conductance. These values include a 30 percent adjustment for geometry effects. The estimated values for hydraulic conductivities of the alluvium are based on assumed values of 0.25 for specific yield and 0.000001 per foot for specific storage of the alluvium; the values assumed for bedrock are 0.1 foot per day horizontal hydraulic conductivity, 0.005 foot per day vertical hydraulic conductivity, and 0.0000001 per foot for specific storage. The resulting diffusivity for the alluvial aquifer is 1,600 feet per day. The estimated values of these hydraulic properties are nearly proportional to the assumed value of specific yield. These values were not found to be sensitive to the assumed values for bedrock. The hydrologic parameters estimated using the cross-sectional model are only valid when taken in context with the other values (both estimated and assumed) used in this study. The model simulates horizontal and vertical flow directions near the river during periods of varying river stage. This information is useful for interpreting bank-storage effects, including the flow of contaminants in the aquifer near the river.

Alaska↗
USGS2021

Development of a groundwater-simulation model in the Los Angeles Coastal Plain, Los Angeles County, California

Executive Summary The Los Angeles Coastal Plain (LACP) covers about 580 square miles and is the largest coastal plain of semiarid southern California. The LACP is heavily developed with mostly residential, commercial, and industrial land uses that rely heavily on groundwater for water supply. In 2010, the LACP was home to about 14 percent of California’s population, or about 5.4 million residents. The LACP is also a major commercial and industrial hub with industries including manufacturing, aerospace, entertainment, and tourism. There has been a heavy reliance on groundwater from the LACP for many years. An average of 305,000 acre-feet per year (acre-ft/yr) of groundwater was used annually from the LACP from 1971 to 2015. The need to replenish the groundwater basins within the LACP was recognized as far back as the 1930s, when spreading grounds were first used to replenish groundwater basins and store water underground during times of water surplus to meet demands in times of shortage. Seawater intrusion resulting from freshwater pumping was first observed in the 1940s. As a result, injection of imported water through wells at what is now the West Coast Basin Barrier Project began on an experimental basis in 1951. Managed aquifer recharge from the spreading grounds and barrier wells is now a substantial component of the LACP’s groundwater supply. The average annual recharge from water spreading from 1971 to 2015 was about 120,000 acre-ft/yr, and the average annual injection into the barrier wells was about 33,000 acre-ft/yr. Other inflows include areal recharge, underflow from San Gabriel and San Fernando Valleys, and onshore flow from the ocean. The average annual recharge from these sources was 100,000 acre-feet (acre-ft) from 1971 to 2015. Additionally, cross-boundary flow from Orange County into the western Orange County subareas of the LACP was simulated as 48,000 acre-ft from 1971 to 2015. This study, conducted in cooperation with the Water Replenishment District of Southern California (WRD), involved an assessment of the historical and present status of groundwater resources in the LACP and the development of tools to better understand the groundwater system. These efforts were built upon results from previous studies and incorporate new information and developments in modeling capabilities to provide a more detailed analysis of the aquifer systems. This study includes a comprehensive compilation of geologic and hydrologic data (Chapter A), development of a chronostratigraphic model that provides a detailed description of the LACP aquifer systems (Chapter B), characterization of the groundwater hydrology of the LACP, including a down-hole analysis of grain size using lithologic and geophysical logs (Chapter C), and development and application of the Los Angeles Coastal Plain Groundwater-flow Model (LACPGM) to simulate past groundwater conditions, estimate groundwater-budget components and flow paths, and approximate future groundwater conditions under different scenarios (Chapter D).

California↗
USGS1984

Hydrogeology and effects of tailings basins on the hydrology of Sands Plain, Marquette County, Michigan

Sands Plain, a 225-square mile area, is near the Marquette iron-mining district in Michigan's Upper Peninsula. Gribben Basin, a settling basin for disposal of waste rock particles from iron-ore concentration, is in the western part. Because Sands Plain is near iron-ore deposits, but not underlain by them, parts of the area are being considered as sites for additional tailings basins. Glacial deposits, as much as 500 feet thick, comprise the principal aquifer. Most ground water flows through the glacial deposits and discharges in a series of nearly parallel tributaries to the Chocolay River which flows into Lake Superior. Ninety-five percent of the discharge of these streams is ground-water runoff. The aquifer is recharged by precipitation at an average rate of 15 inches per year and by streamflow losses from the upper reaches of Goose Lake Outlet at an average rate of 2 inches per year. Precipitation collected at two sites had mean pH values of 4.0; rates of deposition of sulfate and total dissolved nitrogen were estimated to be 17.4 and 5.8 pounds per acre per year, respectively. Dissolved-solids concentrations in water from streams ranged from 82 to 143 milligrams per liter; sulfate ranged from 4.2 to 10 milligrams per liter. Calcium and bicarbonate were the principal dissolved substances. Highest dissolved-solids concentrations in water from wells in glacial deposits were found in a major buried valley east of Goose Lake Outlet. These concentrations ranged from 14 to 246 milligrams per liter; sulfate concentrations ranged from 0.9 to 53 milligrams per liter. Because of the high ground-water component of streamflow, mean concentrations of total nitrogen and trace metals in surface water do not differ significantly from mean concentrations in ground water. A two-dimensional digital model of ground-water flow was used to simulate water levels and ground-water runoff under steady-state and transient conditions Predictive simulations with the steady-state model were made to determine the effects of continued operation of Gribben tailings basin and construction and operation of four hypothetical tailings basins. Operation of Gribben Basin has decreased the average rate of ground-water flow to Goose Lake Outlet by 0.9 to 1.6 cubic feet per second but has increased the average rate of groundwater flow to Warner Creek by about 0.2 cubic foot per second. Continued filling of the tailings basin to its design capacity is expected to cause a slight increase in leakage from the basin to Goose Lake Outlet. Four hypothetical tailings basins, comprising a total of 11 square miles, were simulated by successively adding one more basin to the previous basin configuration. Net ground-water flow to streams was reduced by the simulated basins. The magnitude of these reductions depends on engineering decisions about the method of basin construction and a better understanding of the hydraulic properties of the materials used to seal the basin perimeters. The maximum total reduction in ground-water runoff due to construction and operation of 11 square miles of tailings basins is about 18 cubic feet per second compared to flow simulated by a steady-state simulation without tailings basins. If bottom sealing, rather than slurry wall construction, is used for one of the hypothetical basins, the total maximum reduction is 7.5 cubic feet per second. Under some assumed conditions, leakage from the tailings basins may slightly increase ground-water flow to Goose Lake Outlet and Warner Creek. The maximum probable leakage from all tailings basins is about 7 cubic feet per second; the minimum probable leakage is about 0.7 cubic foot per second.

Michigan↗
USGS1991

Hydrology of Heber and Round Valleys, Wasatch County, Utah, with emphasis on simulation of ground-water flow in Heber Valley

An investigation of the hydrologic system in Heber and Round Valleys was conducted to improve understanding of the surface-water and ground-water hydrology and the effects caused by changes in recharge. Ground water is present in consolidated rocks and in unconsolidated valley-fill deposits, but the principal ground-water reservoir is in the unconsolidated valley-fill deposits. Recharge to the unconsolidated valley-fill deposits in Heber Valley from unconsumed irrigation water, stream infiltration, subsurface inflow from consolidated rocks, and precipitation is estimated to be 154 cubic feet per second. Discharge is by leakage to Deer Creek Reservoir, by springs and seeps, by seepage to the Provo River and other streams, by evapotranspiration, and by pumping from wells. Recharge to the unconsolidated valley-fill deposits in Round Valley from stream infiltration, precipitation, unconsumed irrigation water and subsurface inflow from consolidated rocks is estimated to be 11 cubic feet per second. Discharge is by springs and seeps, by evapotranspiration, and by pumping from wells. Seasonal water-level fluctuations of up to 30 feet occur primarily because of changes in recharge from unconsumed irrigation water. Water levels generally are highest during June or July when recharge from irrigation is at a maximum and lowest during the winter when irrigation is absent and recharge is at a minimum. Water levels in wells near Deer Creek Reservoir respond to changes in the reservoir level. A modular, three-dimensional, finite-difference ground-water flow model developed by McDonald and Harbaugh (1988) was used to simulate the hydrologic system in the unconsolidated valley-fill deposits of Heber Valley. Model simulations indicate that decreased recharge to the unconsolidated valley-fill deposits causes a decrease in discharge to springs and seeps, streams, and leakage to Deer Creek Reservoir. Future decreases in ground-water recharge caused by changing from flood- to sprinkler-irrigation methods will cause future decreases in ground-water discharge that will be offset to some extent by increased surface-water flows.

Utah↗
USGS1986

Simulation of ground-water flow and infiltration from the Susquehanna River to a shallow aquifer at Kirkwood and Conklin, Broome County, New York

A four-layer finite difference model was developed to simulate groundwater flow and induced infiltration to an aquifer underlying the Susquehanna River in the Towns of Kirkwood and Conklin in Broome County, NY. The aquifer consists of sand and gravel deposited in an ancestral river valley during the recession of glacial ice and is in hydraulic connection with the Susquehanna River. In 1984, he aquifer supplied 1.2 million gal/day to well fields in Kirkwood and Conklin. Horizontal hydraulic conductivity of the sand and gravel in the calibrated model ranges from 50 to 10,000 ft/day. Vertical hydraulic conductivity ranges from 1.0 to 80 ft/day. The riverbed thickness was estimated from results of piezometer tests to be 2 ft; the hydraulic conductivity of the riverbed was estimated to be 0.2 ft/day. Root-mean-square differences between computed drawdowns and drawdowns measured in observation wells and piezometers during aquifer tests at the Kirkwood well field ranged from 17% to 24%. The sizes of the well field catchment areas were estimated from a model generated flow net showing the direction and rate of groundwater flow. The Kirkwood catchment area was estimated to be 250 acres, and the Conklin catchment area was 51 acres. Groundwater budgets computed by steady-state simulations showed that 58% of the groundwater withdrawn by the Kirkwood well field is derived from the Susquehanna River during the periods of low river stage and low recharge. The factor to which induced-infiltration rate and size of well field catchment areas are most sensitive, is riverbed hydraulic conductivity.

New York↗
USGS1999

Sources and chronology of nitrate contamination in spring waters, Suwannee River basin, Florida

A multi-tracer approach, which consisted of analyzing water samples for n aturally occurring chemical and isotopic indicators, was used to better understand sources and chronology of nitrate contamination in spring wate rs discharging to the Suwannee and Santa Fe Rivers in northern Florida. Dur ing 1997 and 1998, as part of a cooperative study between the Suwannee River Water Management District and the U.S. Geological Survey, water samples were collected and analyzed from 24 springs and two wells for major ions, nutrients, dissolved organic carbon, and selected environmental isotopes [18O/16O, D/H, 13C/12C, 15N/14N]. To better understand when nitrate entered the ground-water system, water samples were analyzed for chlorofluorocarbons (CFCs; CCl3F, CCl2F2, and C2Cl3F3) and tritium (3H); in this way, the apparent ages and residence times of spring waters and water from shallow zones in the Upper Floridan aquifer were determined. In addition to information obtained from the use of isotopic and other chemical tracers, information on changes in land-use activities in the basin during 1954-97 were used to estimate nitrogen inputs from nonpoint sources for five counties in the basin. Changes in nitrate concentrations in spring waters with time were compared with estimated nitrogen inputs for Lafayette and Suwannee Counties. Agricultural activities [cropland farming, animal farming operations (beef and dairy cows, poultry, and swine)] along with atmospheric deposition have contributed large quantities of nitrogen to ground water in the Suwannee River Basin in northern Florida. Changes in agricultural land use during the past 40 years in Alachua, Columbia, Gilchrist, Lafayette, and Suwannee Counties have contributed variable amounts of nitrogen to the ground-water system. During 1955-97, total estimated nitrogen from all nonpoint sources (fertilizers, animal wastes, atmospheric deposition, and septic tanks) increased continuously in Gilchrist and Lafayette Counties. In Suwannee, Alachua, and Columbia Counties, estimated nitrogen inputs from all nonpoint sources peaked in the late 1970's corresponding to the peak use in fertilizer during this time. Fertilizer use in Columbia, Gilchrist, Lafayette, and Suwannee Counties increased substantially during 1993-97. The heavy use of fertilizers in the basin is corroborated by nitrogen isotope data. Values of d15N of nitrate (d15N-NO3) in spring waters range from 2.7 per mil (SUW718971) to 10.6 per mil (Poe Spring) with a median of 5.4 per mil. The range of values indicates that nitrate in the sampled spring waters most likely originates from a mixture of inorganic (fertilizers) and organic (animal wastes) sources; however, higher d15N values for Poe and Lafayette Blue Springs indicate that an organic source of nitrogen probably is contributing nitrate to these spring waters. Water samples from two wells sampled in Lafayette County have high d15N-NO3 values of 11.0 and 12.1 per mil, indicating the predominance of an organic source of nitrate. These two wells are located near dairy and poultry farms, where leachate from animal wastes may contribute nitrate to ground water. Dissolved-gas data (nitrogen, argon, and oxygen) indicate that denitrification has not removed large amounts of nitrate from the ground-water system. Thus, variations in d15N-NO3 values of spring waters can be attributed to variations in d15N-NO3 values of ground-water recharge, and can be used to obtain information about source(s) of nitrate. Extending the use of age-dating techniques (CFCs and 3H) to spring waters in complex karst systems required the use of several different approaches for estimating age and residence time of ground water discharging to springs. These approaches included the use of a simple reservoir model, a piston-flow model, an exponential model, and a binary-mixing model. When age data (CFC-11, CFC-113, and 3H) are combined for all springs, models that incorporate exponential mixtures seem to provide reliable estimates of average residence times of ground water discharging to springs. Whereas, data for some individual springs fit a binary-mixing model with more than 50 percent young water (recharged within the past 5 years), data from other individual springs fit a piston-flow model with a water age of about 25 years. The young ages of several spring waters (such as SUW718971, SUW725971, and Ginnie Spring) indicate the high vulnerability of the springs to contamination. For most springs, CFCs suggest that a large fraction of the water is more than 20 years old. Springs with lower flows tend to have young ages (shallow ground-water flow systems), whereas springs with higher flows tend to have older ages (deep ground-water flow systems). The chemical composition of spring waters can be used as a qualitative indicator of age and ground-water residence time. Nitrate-nitrogen concentrations and dissolved oxygen in spring waters are inversely related to the apparent ages of spring waters and ground-water residence time in the basin. Silica concentrations increase as the age of spring waters increase. Long-term trends in nitrate concentrations in selected spring waters were compared with estimated inputs of nitrogen from various sources in Suwannee and Lafayette Counties. In both counties, trends in nitrate concentrations in spring waters closely followed the estimated contribution of nitrogen from fertilizers. Decreasing nitrate concentrations in spring waters from Suwannee County followed the decrease in estimated fertilizer use from the mid-1970's to the early 1990's. Increasing nitrate concentrations in sprin g waters from Lafayette County followed the steady increase in fertilizer use from the early 1960's to the mid-1990's The relation between the concentration of nitrate in ground water and the amount of nitrogen that is added to a ground-water contributing area for a spring is controlled by complex interactions among hydrogeologic, land-use, and climatic factors, as well as several other land-management factors. Spring waters represent mixtures of converging flow paths that contain ground water with a range of ages. Even if nitrogen inputs were reduced substantially, it may take decades for nitrate concentrations in the ground-water system to return to near background levels.

Florida↗
USGS1984

A finite-element simulation model for saturated-unsaturated, fluid-density-dependent ground-water flow with energy transport or chemically- reactive single-species solute transport

SUTRA (Saturated-Unsaturated Transport) is a computer program which can be used to simulate the movement of fluid and the transport of either energy or dissolved substances in a subsurface environment. The model employs a two-dimensional hybrid finite-element and integrated-finite-difference method to approximate the governing equations that describe the two interdependent processes that are simulated by SUTRA: (1) fluid density-dependent saturated or unsaturated groundwater flow, and either (2a) transport of a solute in the groundwater, in which the solute may be subject to: equilibrium adsorption on the porous matrix, and both first-order and zero-order production or decay, or, (2b) transport of thermal energy in the groundwater and solid matrix of the aquifer. SUTRA provides, as the primary calculated results, fluid pressures and either solute concentrations or temperatures, as they vary with time, everywhere in the simulated subsurface system. SUTRA may also be used to simulate simpler subsets of the above process. SUTRA may be employed for areal and cross-sectional models of saturated groundwater flow systems, and for cross-sectional models of unsaturated zone flow. Solute transport simulation using SUTRA may be used to simulate natural or man-induced chemical transport, solute sorption, production and decay. SUTRA may be used for simulation of variable density leachate movement, and for cross-sectional simulation of salt-water intrusion in aquifers at near-well or regional scales, with either dispersed or relatively sharp transition zones between fresh water and salt water. SUTRA energy transport simulation may be employed to model thermal regimes in aquifers, subsurface heat conduction, aquifer thermal energy storage systems, geothermal reservoirs, thermal pollution of aquifers, and natural hydrogeologic convection systems. (USGS)

Water-Resources Investigations Report↗
USGS2005

Precipitation-runoff processes in the Feather River basin, northeastern California, and streamflow predictability, water years 1971-97

Precipitation-runoff processes in the Feather River Basin of northern California determine short- and long-term streamflow variations that are of considerable local, State, and Federal concern. The river is an important source of water and power for the region. The basin forms the headwaters of the California State Water Project. Lake Oroville, at the outlet of the basin, plays an important role in flood management, water quality, and the health of fisheries as far downstream as the Sacramento-San Joaquin Delta. Existing models of the river simulate streamflow in hourly, daily, weekly, and seasonal time steps, but cannot adequately describe responses to climate and land-use variations in the basin. New spatially detailed precipitation-runoff models of the basin have been developed to simulate responses to climate and land-use variations at a higher spatial resolution than was available previously. This report characterizes daily rainfall, snowpack evolution, runoff, water and energy balances, and streamflow variations from, and within, the basin above Lake Oroville. The new model's ability to predict streamflow is assessed. The Feather River Basin sits astride geologic, topographic, and climatic divides that establish a hydrologic character that is relatively unusual among the basins of the Sierra Nevada. It straddles a north-south geologic transition in the Sierra Nevada between the granitic bedrock that underlies and forms most of the central and southern Sierra Nevada and volcanic bedrock that underlies the northernmost parts of the range (and basin). Because volcanic bedrock generally is more permeable than granitic, the northern, volcanic parts of the basin contribute larger fractions of ground-water flow to streams than do the southern, granitic parts of the basin. The Sierra Nevada topographic divide forms a high altitude ridgeline running northwest to southeast through the middle of the basin. The topography east of this ridgeline is more like the rain-shadowed basins of the northeastern Sierra Nevada than the uplands of most western Sierra Nevada river basins. The climate is mediterranean, with most of the annual precipitation occurring in winter. Because the basin includes large areas that are near the average snowline, rainfall and rain-snow mixtures are common during winter storms. Consequently, the overall timing and rates of runoff from the basin are highly sensitive to winter temperature fluctuations. The models were developed to simulate runoff-generating processes in eight drainages of the Feather River Basin. Together, these models simulate streamflow from 98 percent of the basin above Lake Oroville. The models simulate daily water and heat balances, snowpack evolution and snowmelt, evaporation and transpiration, subsurface water storage and outflows, and streamflow to key streamflow gage sites. The drainages are modeled as 324 hydrologic-response units, each of which is assumed homogeneous in physical characteristics and response to precipitation and runoff. The models were calibrated with emphasis on reproducing monthly streamflow rates, and model simulations were compared to the total natural inflows into Lake Oroville as reconstructed by the California Department of Water Resources for April-July snowmelt seasons from 1971 to 1997. The models are most sensitive to input values and patterns of precipitation and soil characteristics. The input precipitation values were allowed to vary on a daily basis to reflect available observations by making daily transformations to an existing map of long-term mean monthly precipitation rates that account for altitude and rain-shadow effects. The models effectively simulate streamflow into Lake Oroville during water years (October through September) 1971-97, which is demonstrated in hydrographs and statistical results presented in this report.

Northeastern California↗
USGS1995

Hydrogeology and simulated effects of ground-water withdrawals for citrus irrigation, Hardee and De Soto counties, Florida

The hydrogeology of Hardee and De Soto Counties in west-central Florida was evaluated, and a ground-water flow model was developed to simulate the effects of expected increases in ground-water withdrawals for citrus irrigation on the potentiometric surfaces of the intermediate aquifer system and the Upper Floridan aquifer. In 1988, total citrus acreage in Hardee and De Soto Counties was 89,041 acres. By the year 2020, citrus acreage is projected to increase to 130,000 acres. Ground water is the major source of water supply in the study area, and 94 percent of the ground-water withdrawn in the area is used for irrigation purposes. The principal sources of ground water in the study area are the surficial aquifer, the intermediate aquifer system, and upper water-yielding units of the Floridan aquifer system, commonly referred to as the Upper Floridan aquifer. The surficial aquifer is a permeable hydrogeo1ogic unit contiguous with land surface that is comprised predominately of surficial quartz sand deposits that generally are less than 100 feet thick. The intermediate aquifer system is a somewhat less permeable hydrogeologic unit that lies between and retards the exchange of water between the overlying surficial aquifer and the underlying Upper Floridan aquifer. Thickness of the intermediate aquifer system ranges from about 200 to 500 feet and transmissivity ranges from 400 to 7,000 feet squared per day. The highly productive Upper Floridan aquifer consists of 1,200 to 1,400 feet of solution-riddled and fractured limestone and dolomite. Transmissivity values for this aquifer range from 71,000 to 850,000 feet squared per day. Wells open to the Upper Floridan aquifer. the major source of water in the area, can yield as much as 2,500 gallons of water per minute. The potential effects of projected increases in water withdrawals for citrus irrigation on groundwater heads were evaluated by the use of a quasi-three-dimensional, finite-difference, ground-water flow model. The model was calibrated under steady-state conditions to simulate September 1988 heads and under transient conditions to simulate head fluctuations between September 1988 and September 1989. The calibrated model was then used to simulate hydraulic heads for the years 2000 and 2020 that might result from projected increases in pumpage for citrus irrigation. The model simulation indicated that increased pumpage might be expected to result in: A maximum decline of more than 10 feet in theintermediate aquifer system at a proposed grove in eastern De Soto County and an average decline of more than 2 feet in much of the study area. An increase in downward leakage to the intermediate aquifer system from the overlying surficial aquifer system from 178 to 183 million gallons per day. A decrease in upward leakage from the intermediate aquifer system to the surficial aquifer from 1.58 to 1.47 million gallons per day. A maximum decline of about 5 feet in the Upper Floridan aquifer at a proposed grove in eastern De Soto County and a decline of more than 2 feet in much of the model area. An increase in downward leakage to the Upper Floridan aquifer from the intermediate aquifer system from 180 to 183 million gallons per day. A decrease in upward leakage from the Upper Floridan aquifer to the intermediate aquifer system from 4.32 million gallons per day in 1989 to 3.89 million gallons per day in the year 2,000. but an increase in upward leakage to 5.10 million gallons per day by the year 2020, reflecting a change in hydraulic gradient.

Water-Resources Investigations Report↗
USGS2005

An assessment of optical properties of dissolved organic material as quantitative source indicators in the Santa Ana River basin, Southern California

The ability to rapidly, reliably, and inexpensively characterize sources of dissolved organic material (DOM) in watersheds would allow water management agencies to more quickly identify problems in water sources, and to more efficiently allocate water resources by, for example, permitting real-time identification of high-quality water suitable for ground-water recharge, or poor-quality water in need of mitigation. This study examined the feasibility of using easily measurable intrinsic optical properties' absorbance and fluorescence spectra, as quantitative indicators of DOM sources and, thus, a predictor of water quality. The study focused on the Santa Ana River Basin, in southern California, USA, which comprises an area of dense urban development and an area of intense dairy production. Base flow in the Santa Ana Basin is primarily tertiary treated wastewater discharge. Available hydrologic data indicate that urban and agricultural runoff degrades water quality during storm events by introducing pathogens, nutrients, and other contaminants, including significant amounts of DOM. These conditions provide the basis for evaluating the use of DOM optical properties as a tracer of DOM from different sources. Sample spectra representing four principal DOM sources were identified among all samples collected in 1999 on the basis of basin hydrology, and the distribution of spectral variability within all the sample data. A linear mixing model provided quantitative estimates of relative endmember contribution to sample spectra for monthly, storm, and diurnal samples. The spectral properties of the four sources (endmembers), Pristine Water, Wastewater, Urban Water, and Dairy Water, accounted for 94 percent of the variability in optical properties observed in the study, suggesting that all important DOM sources were represented. The scale and distribution of the residual spectra, that not explained by the endmembers, suggested that the endmember spectra selected did not adequately represent Urban Water base flow. However, model assignments of sources generally agreed well with those expected, based on sampling location and hydrology. The results suggest that with a fuller characterization of the endmember spectra, analysis of optical properties will provide rapid quantitative estimates of the relative contribution of DOM sources in the Santa Ana Basin.

Scientific Investigations Report↗