Santa Barbara and Foothill Groundwater Basins Geohydrology and Optimal Water Resources Management— Developed using Density Dependent Solute Transport and Optimization Models

U.S. Geological Survey (USGS) | July 10th, 2018

Groundwater has been a part of the city of Santa Barbara’s water-supply portfolio since the 1800s; however, since the 1960s, the majority of the city’s water has come

Groundwater has been a part of the city of Santa Barbara’s water-supply portfolio since the 1800s; however, since the 1960s, the majority of the city’s water has come from local surface water, and the remainder has come from groundwater, State Water Project, recycled water, increased water conservation, and as needed, seawater desalination. Although groundwater from the Santa Barbara and Foothill groundwater basins only accounts for a small percentage of the long-term supply, it is an important source of supplemental water during times of surface-water shortages. During the late 1980s and early 1990s, production wells extracted additional groundwater to compensate for drought related water-delivery shortfalls from other sources; in response, water levels declined substantially in the Santa Barbara and Foothill groundwater basins (below sea level in the Santa Barbara groundwater basin). In coastal basins that have groundwater extraction near shore, seawater intrusion is often a problem. Seawater intrusion in the Santa Barbara groundwater basin is thought to be more limited than in other coastal basins because of an offshore fault that acts as a partial barrier to groundwater flow. During the late 1980s and early 1990s, seawater intrusion was observed in the Santa Barbara groundwater basin, as indicated by increased chloride concentrations at several monitoring wells that ranged from 200 ft to 1,300 ft from the ocean and as close as 2,900 ft to the nearest pumping well. This demonstrated that seawater can intrude into the Santa Barbara groundwater basin when groundwater levels fall below sea level near the coast. The city of Santa Barbara is interested in developing a better understanding of the sustainability of its groundwater supplies. In 2014, California adopted historic legislation to manage its groundwater: the Sustainable Groundwater Management Act (SGMA). The SGMA requires the development and implementation of “Groundwater Sustainability Plans” in 127 priority groundwater basins; although Santa Barbara was not a designated priority basin, the city is taking steps to achieve sustainability. Sustainability was defined in the SGMA in terms of avoiding undesirable results: significant and unreasonable groundwater-level declines, reduction in groundwater storage, seawater intrusion, water-quality degradation, land subsidence, and surface-water depletion. In this project, a cooperative study between the U.S. Geological Survey (USGS) and the city of Santa Barbara, sustainable yield is defined as the volume of groundwater that can be pumped from storage without causing water-level drawdowns and the associated increases in seawater intrusion (as indicated by increases in measured chloride concentrations) at selected wells. In order to estimate the sustainability of Santa Barbara’s groundwater basins, a three-dimensional density-dependent groundwater-flow and solute-transport model (the Santa Barbara Flow and Transport Model, or SBFTM) was developed on the basis of an existing groundwater-flow model. To simulate seawater intrusion to the Santa Barbara Basin under various management strategies, the SBFTM uses the USGS code SEAWAT to simulate salinity transport and variable-density flow. The completed SBFTM was coupled with a management optimization tool, in this case a multi-objective evolutionary algorithm, to determine optimal pumping strategies that maximize the sustainable yield and at the same time satisfy user-defined drawdown and chloride-concentration constraints. As part of this study, a three-dimensional hydrogeologic framework model was developed to quantify the extent and hydrogeologic characteristics of the Santa Barbara and Foothill groundwater basins and to help define the discretization and hydraulic properties used in the SBFTM. The development of the hydrogeologic framework model required the collection and reconciliation of geologic and geophysical data from existing maps, reports, and databases, along with geologic and hydrologic data from recently drilled wells. These data were integrated into a three-dimensional hydrogeologic framework model that defines the stratigraphy and geometry of the aquifer zones and the major geologic structures in the basin. The hydrogeologic framework model also quantifies the variation in sediment grain size within each aquifer zone as the percentage of coarse-grained sediment. Previous studies indicated that there are two principal water-producing zones in the Santa Barbara groundwater basin, the upper and lower producing zones; an additional thin, productive zone was identified as part of this study. This “middle producing zone” is not as areally extensive as the upper and lower producing zones and only exists in the coastal part of Storage Unit I. These producing zones are bounded at depth by less productive shallow, middle, and deep zones. Santa Barbara and Foothill Groundwater Basins Geohydrology and Optimal Water Resources Management—Developed using Density Dependent Solute Transport and Optimization Models Two versions of the SBFTM were constructed: an initial-condition model and a modern transient model. The initial-condition model is a long-term transient model that simulates flow and solute-transport conditions during a period with limited anthropogenic influences preceeding the modern transient model. The simulation-transient model simulates flow and transport conditions from 1929 through 2013; however, because of data availability, the focus of the model calibration was 1972–2013. The SBFTM was calibrated to measured groundwater levels and drawdown, as well as measured chloride concentrations and change in concentrations, using a combination of automated and trial-and-error parameter-estimation techniques. A sensitivity analysis indicated that, in general, the SBFTM was most sensitive to recharge- and pumping-distribution parameters, specifically those controlling the amount of small-catchment recharge and the distribution of water extraction by hydrogeologic layer for production wells. The model was also sensitive to parameters controlling stream-recharge rates, horizontal and vertical hydraulic conductivity, and porosity. From 1929 to 1971, most of the water entering the area represented by the SBFTM was from creek and small-catchment recharge, and the majority of water leaving the SBFTM area was from pumping, discharge to creeks, and drains. In addition, about 37 percent of the total pumpage came from a net reduction in groundwater storage. From 1972 to 2013, the amount of water entering and leaving the SBFTM was fairly similar as that from 1929 to 1971, except the reduction in pumpage added about 17,000 acre-ft of water to storage. During this later period, there were also times of storage loss. For example, during July 1990, a month when approximately 705 acre-ft of groundwater was pumped in the study area, the pumpage was much greater than all sources of recharge combined, and about 382 acre-ft of water was removed from groundwater storage. Simulated hydraulic heads replicated the observed data to an acceptable matching of the measured water-level, flow direction, and vertical gradients. Simulated hydrographs for selected wells were in good agreement with the measured data, with an average residual of -2.7 ft and a standard deviation of 14.5 ft, indicating that the simulated heads, on average, underestimated the observed water levels. An examination of the model fit indicated that most of the discrepancies were lower simulated heads at wells proximal to production well sites. The simulated chloride concentrations reasonably matched the rising limbs of the measured breakthrough curves in terms of timing and magnitude; however, the simulation overestimated the chloride concentrations on the falling limbs. The overestimation of low chloride concentrations was attributed to the model overestimating the advance of the chloride front during periods of heavy pumping and underestimating the retreat of the chloride front during periods of low pumping. These simulation errors would result in a conservative response by local water managers to seawater intrusion. The SBFTM was used to develop a collection of predictive simulations optimized to produce pumping schedules that maximize yield, subject to a set of constraints and competing objectives. The simulations were grouped as scenarios that differed in their time horizon, initial conditions for groundwater levels and chloride concentrations, as well as precipitation, which was incorporated into the model through simulated recharge. Overall, five scenarios were developed in a multi-objective framework to obtain optimal pumping rates for all of the wells managed by the city, while minimizing excessive drawdown and seawater intrusion. For the current study, complexities in the simulation model and the optimization formulation required additional considerations. Incorporating the solute-transport equations to simulate chloride transport added a highly nonlinear process that is solved iteratively in each time step of the groundwater-flow model. These nonlinearities, coupled with the highly refined grid in the current model, creates challenges for many traditional optimization methods. Therefore, an optimization method was needed that could address nonlinear relationships as well as a very large problem size. Lastly, the optimization problem was reformulated to include multiple objectives without requiring convergence to a single solution. This approach, guided by the city’s objectives, allowed the maximum extraction of information from the complex simulation. Borg, a multi-objective evolutionary algorithm, was chosen as the optimization algorithm for this study for several reasons: (1) it is very computationally efficient; (2) it can run in parallel; (3) it requires little user input; and (4) it can solve for multiple competing objectives. The first three points allow the algorithm to proceed toward the optimal solutions at the fastest possible rate. The fourth point is advantageous for large, complex optimization problems because it is difficult to formulate the optimization problem in a way that produces only one optimal solution. The problem formulation consisted of four competing objectives and a constraint set in accordance with the main concerns of the city. The objectives were maximizing total pumpage, minimizing seawater intrusion, minimizing total drawdown in production wells, and minimizing the maximum drawdown. The constraints were pump capacity, meeting drinking-water standards for chloride, maintaining a specified minimum flowrate to a groundwater treatment plant, and maintaining minimum water levels in pumping wells. The decision variables either were quarterly pumpage by well or total pumpage by basin. Five optimization scenarios were developed that allow thedecision makers to evaluate a range of optimal solutions for a variety of water levels and chloride concentrations as well as potential future climatic conditions. Three scenarios (1, 2, and 5) were multi-objective optimization formulations that allowed for variations in management preferences and climatic conditions. The other two scenarios (3 and 4) were designed to examine the optimization results to answer specific questions. Scenario 1 described the best-case sustainable yield assuming a “full” basin (that is, high initial water levels) and typical climate conditions for 10 years. Scenario 2 also started with a “full” basin; however, this was followed by a 10-year drought. Scenario 3 determined if an “empty” basin (that is, low initial water levels) would recover to full conditions (1998 conditions) given climate assumptions and optimal pumping schedules from scenarios 1 and 2. Scenario 4 was designed to produce decision rules that can be used by water managers to help choose an optimal pumping schedule based on measured water-level or chloride data. Scenario 5 identified future pumping schedules based on short-term climate variations during a 2-year management horizon. The results from scenarios 1 and 2 described the differences in maximum pumpage in the basin under typical and dry long-term climate projections, respectively. The scenario 1 results indicated the maximum 10-year pumpage of the basin was about 31,300 acre-ft under typical conditions and controlling simulated seawater intrusion and drawdowns. For scenario 2, less recharge over the 10-year dry climate produced a maximum pumpage estimate of 30,000 acre-ft to control seawater intrusion and drawdowns. The larger pumpage for scenario 1 resulted in more seawater intrusion, but less total drawdown, compared to that of scenario 2. Results for scenarios 3 and 4 showed the basin’s response to management actions combined with climate projections. Both scenarios used the optimal pumping schedules and the 10-year climates from scenarios 1 and 2. The scenario 3 results showed that under minimal pumping, the basin did not fully recover to 1998 water levels within 10 years under either climate scenario. The relatively larger recharge from the typical climate resulted in less drawdown at coastal monitoring wells after the 10-year recovery period than that from the dry climate. The location of the seawater intrusion front was not appreciably different between the scenarios, however. Scenario 4 used the optimal results from scenarios 1 and 2 to produce decision-rule curves that illustrated the pumpage for each basin, given measured levels of chloride concentration or drawdown. This allowed the use of additional measurements at monitoring wells to assess future management decisions on the basis of the sensitivity of observations of drawdown and seawater intrusion to various pumping rates. Scenario 5 allowed managers to investigate the effects of short-term climate variations on optimal pumping schedules. Three specific 2-year simulations were optimized: typical-to-dry (scenario 5A), dry-to-typical (scenario 5B), and dry-to-dry (scenario 5C). The most noteable result from scenario 5 was the overall reduction in optimal pumpage for most schedules in scenario 5C, when the climate is simulated as dry-to-dry. There are also many optimal pumping schedules that produced an overall increase in waterlevels over the two-year simulation period, regardless of climatic condition. Similar to scenario 2, the scenario 5C results represents conservative yield estimates under a minimal-precipitation climatic condition.

Santa Cruz Mid-County Groundwater Basin Groundwater Sustainability Plan

Santa Cruz Mid-County Groundwater Basin Groundwater Sustainability Agency | November 21st, 2019

Due to size, this file has been divided Part 1, Part 2, Part 3, Part 4, Part 5

Due to size, this file has been divided Part 1, Part 2, Part 3, Part 4, Part 5

Satellite-based monitoring of groundwater depletion in California’s Central Valley

Nature Portfolio (Springer Nature) | November 5th, 2019

Range change data, obtained from Synthetic Aperture Radar satellites, form the basis for estimates of aquifer volume change in california’s central Valley. th

Range change data, obtained from Synthetic Aperture Radar satellites, form the basis for estimates of aquifer volume change in california’s central Valley. the estimation algorithm incorporates a function penalizing changes far from known well locations, linking the aquifer volume changes to agricultural, industrial, and municipal pumping within the tulare basin. We show that the range changes are compatible with the hypothesis that the source of aquifer volume changes are variations in effective pressure around documented wells. Specifically, inclusion of the well distance penalty does not degrade the fit to the observations, inversions with and without it both give variance reductions of 99.6%. The patterns of aquifer volume change vary significantly from the drought year, between October 2015 and October 2016, to a wet year in 2017, and into 2018, a year with near average rainfall. The 2.3 million acre-feet of estimated volume reduction, a lower bound on the amount of water extracted from the basin between October 2015 and 2016, agrees with independent estimates of 1.8 and 2.3 million acre- feet. The aquifer volume reduction is also compatible with a loss of 3.1 km3 (2.5 million acre-feet) in groundwater volume derived from Gravity Recovery and Climate Experiment (GRACE) satellite data.

Scenarios of Future California Water Demand Through 2050

California Department of Water Resources (DWR) | August 14th, 2014

This paper describes the application of future scenarios for Water Plan Update 2013 (DWR 2013) to quantify future water demand for California and our 10 hydrologic region

This paper describes the application of future scenarios for Water Plan Update 2013 (DWR 2013) to quantify future water demand for California and our 10 hydrologic regions out to the year 2050. Starting with the publication of California Water Plan Update 2005 (DWR 2005), DWR has applied the concept of scenarios to describe future water management conditions. Scenarios are used to capture a broad range of uncertain factors that affect water management that are considered to be beyond the control of water managers. For Update 2013 scenarios were organized around narrative themes of growth and climate change, which were then quantified using factors about future population growth, housing densities, land use patterns, and climate.  For Updates 2005 and 2009, DWR attempted only to quantify future water demand management strategies at the hydrologic region scale. The hydrologic regions are large areas of tens of thousands square miles, too large to compute meaningful information on regional water management actions. But in parallel to the hydrologic region scale analysis as in current Update 2013, DWR has been pursuing a longer term project, to perform a more complete hydrologic analysis at the finer planning area scale. Results and analysis of planning area water management strategies is given in a separate report.  In addition to this paper, background information on the application of scenarios can be found in Chapter 5 and 6 of Update 2013. Chapter 5 provides information how scenarios are used to quantify water management uncertainty. Chapter 6 describes DWR’s efforts to effect long term improvements in data and analytical tools to support the Water Plan. 

Science and Decision Making: Water Management and Tree‐Ring Data in the Western United States

American Water Resources Association (AWRA) | October 5th, 2009

Growing populations, limited resources, and sustained drought are placing increased pressure on already over‐allocated water supplies in the western United States, prom

Growing populations, limited resources, and sustained drought are placing increased pressure on already over‐allocated water supplies in the western United States, prompting some water managers to seek out and utilize new forms of climate data in their planning efforts. One source of information that is now being considered by water resource management is extended hydrologic records from tree‐ring data. Scientists with the Western Water Assessment (WWA) have been providing reconstructions of streamflow (i.e., paleoclimate data) to water managers in Colorado and other western states (Arizona, New Mexico, and Wyoming), and presenting technical workshops explaining the applications of tree‐ring data for water management for the past eight years. Little is known, however, about what has resulted from these engagements between scientists and water managers. Using in‐depth interviews and a survey questionnaire, we attempt to address this lack of information by examining the outcomes of the interactions between WWA scientists and western water managers to better understand how paleoclimate data has been translated to water resource management. This assessment includes an analysis of what prompts water managers to seek out tree‐ring data, how paleoclimate data are utilized by water managers in both quantitative and qualitative ways, and how tree‐ring data are interpreted in the context of organization mandates and histories. We situate this study within a framework that examines the coproduction of science and policy, where scientists and resource managers collectively define and examine research and planning needs, the activities of which are embedded within wider social and political contexts. These findings have broader applications for understanding science‐policy interactions related to climate and climate change in resource management, and point to the potential benefits of reflexive interactions of scientists and decision makers.

Science be Dammed: How Ignoring Inconvenient Science Drained the Colorado River

University of Arizona Press | November 26th, 2019

Science Be Dammed is an alarming reminder of the high stakes in the management—and perils in the mismanagement—of water in the western United States.  It seems decep

Science Be Dammed is an alarming reminder of the high stakes in the management—and perils in the mismanagement—of water in the western United States.  It seems deceptively simple: even when clear evidence was available that the Colorado River could not sustain ambitious dreaming and planning by decision-makers throughout the twentieth century, river planners and political operatives irresponsibly made the least sustainable and most dangerous long-term decisions. Arguing that the science of the early twentieth century can shed new light on the mistakes at the heart of the over-allocation of the Colorado River, authors Eric Kuhn and John Fleck delve into rarely reported early studies, showing that scientists warned as early as the 1920s that there was not enough water for the farms and cities boosters wanted to build. Contrary to a common myth that the authors of the Colorado River Compact did the best they could with limited information, Kuhn and Fleck show that development boosters selectively chose the information needed to support their dreams, ignoring inconvenient science that suggested a more cautious approach. Today water managers are struggling to come to terms with the mistakes of the past. Focused on both science and policy, Kuhn and Fleck unravel the tangled web that has constructed the current crisis. With key decisions being made now, including negotiations for rules governing how the Colorado River water will be used after 2026, Science Be Dammed offers a clear-eyed path forward by looking back. Understanding how mistakes were made is crucial to understanding our contemporary problems. Science Be Dammed offers important lessons in the age of climate change about the necessity of seeking out the best science to support the decisions we make.

Scoping Report for Post-2026 Colorado River Reservoir Operations Upper and Lower Colorado Basin Regions

U.S. Bureau of Reclamation (USBR) | October 19th, 2023

This scoping summary was prepared to document the scoping process for the Post-2026 Colorado River Reservoir Operations (Post-2026 process). As several important reservoi

This scoping summary was prepared to document the scoping process for the Post-2026 Colorado River Reservoir Operations (Post-2026 process). As several important reservoir and water management decisional documents and agreements that govern operation of Colorado River facilities and management of Colorado River water are scheduled to expire at the end of 2026, the Post-2026 process has been initiated to develop successor domestic Colorado River operational guidelines and strategies for Lake Powell and Lake Mead.

Seasonal Peak Snow Predictability Derived From Early-Season Snow in North America

American Geophysical Union (AGU) | September 12th, 2023

Early season snow observations show promise in predicting peak snow water equivalent in locations where large-scale ocean-atmosphere patterns similarly influence fall and

Early season snow observations show promise in predicting peak snow water equivalent in locations where large-scale ocean-atmosphere patterns similarly influence fall and winter conditions, and in northern maritime locations where winter temperatures are warm enough that rain on snow and midwinter melt occur.

Second Amended Delano-Earlimart Irrigation District Groundwater Sustainability Plan 2024

Delano-Earlimart Irrigation District Groundwater Sustainability Agency (DEID GSA) | July 1st, 2024

In January of 2020, the DEID GSA submitted its Groundwater Sustainability Plan (GSP) to the California Department of Water Resources (DWR). The GSP provided technical doc

In January of 2020, the DEID GSA submitted its Groundwater Sustainability Plan (GSP) to the California Department of Water Resources (DWR). The GSP provided technical documentation of the DEID Management Area’s existing sustainability, along with addressing the needs of the Earlimart Public Utility District and Richgrove Community Services District. The GSP outlined specific management actions that will be taken in accordance with the Sustainable Groundwater Management Act. DWR’s review of the Tule Subbasin, as a whole, resulted in an incomplete determination, to which the GSAs in the Tule Subbasin had 180 days to respond. The GSA’s in the Tule Subbasin provided revised GSPs on July 27, 2022, which DWR also determined were inadequate, causing the Tule Subbasin to transition to the oversight of the State Water Resources Control Board (SWRCB). The draft GSP found below includes updates to specifically address comments received from the SWRCB and their staff in anticipation of the Tule Subbasin probationary hearing on September 17, 2024. Comments on the 2nd amended GSP can be submitted to info@deid.org and will be accepted through August 18, 2024. Due to size, this file has been divided. Part 1, Part 2, Part 3, Part 4, Part 5, Part 6

Second Supplement to State Water Resources Board Bulletin No. 11, San Joaquin County Investigation, Basic Data, Calaveras Unit, 1955-56

California Department of Water Resources (DWR) | April 20th, 1957