Groundwater Storage Estimates in the Central Valley Aquifer Using GRACE Data

EarthzineArticles, Earth Observation, Original, Sections, Themed Articles, Water Availability Theme

A map outlining the Central Valley aquifer.

A map outlining the Central Valley aquifer.

Fig. 1. The Central Valley aquifer is outlined in black and is the study area for which C2VSIM calculates groundwater storage changes.

Amber Kuss, San Francisco State University;

William T. Brandt, California State University, Monterey Bay;

Joshua Randall, Arizona State University;

Bridget Floyd, University of California, Berkeley;

Abdelwahab Bourai, Cupertino High School;

Michelle Newcomer, San Francisco State University;

Cindy Schmidt, NASA Ames Research Center;

J.W. Skiles Ph.D., NASA Ames Research Center

Abstract — The role of the NASA Applied Sciences DEVELOP student internship program is to use NASA satellite missions to explore Earth-based research questions in collaboration with state and federal agencies. One recent project focused on the use of the Gravity Recovery and Climate Experiment (GRACE), to estimate changes in groundwater storage in the Central Valley aquifer in collaboration with the California Department of Water Resources (DWR). GRACE can be used to measure gravitational anomalies on Earth which are further processed to remove other gravitational effects, and are obtained as measurements of equivalent water thickness, which relate to changes in total water storage (TWS) throughout the world. For this study, GRACE TWS anomalies were obtained from October 2002 to September 2009, for two hydrologic regions, the Sacramento River Basin and the San Joaquin River Basin, including the Tulare Lake Basin, encompassing the Central Valley aquifer. To calculate monthly groundwater storage estimates, additional variables such as soil moisture, snowpack, and surface water storage were combined with GRACE TWS values to estimate groundwater storage anomalies. The GRACE-derived changes in groundwater storage at the basin and regional level or the two basins combined, were then compared to modeled values calculated using the California Department of Water Resources’ Central Valley Groundwater-Surface Water Simulation Model (C2VSIM) and a Geographical Information Systems Change in Storage Tool (GIS CST). Groundwater storage estimates from GRACE and C2VSIM were comparable for the entire Central Valley, showing similar seasonal, annual, and long-term trends. This work has the potential to improve California’s groundwater resource management and the use of existing hydrologic models for the Central Valley.

I. INTRODUCTION TO NASA DEVELOP

The NASA Applied Sciences’ DEVELOP Program sponsors paid internships located at Ames Research Center for students to extend science research to local communities. DEVELOP is a NASA Science Mission Directorate Applied Sciences Program training and development internship. Students work on Earth science research projects, mentored by science advisors from NASA and partner agencies to extend research results to local communities. The NASA Ames DEVELOP program hosts graduate, undergraduate and high school students in a ten week summer internship. Students are involved in a wide range of projects such as vegetation mapping, air quality assessments, wetland restoration, and water availability. The DEVELOP program is student run and student lead, and is supported by science advisors and mentors.

Drawing of the GRACE satellite orbiting Earth.

Fig. 2. The GRACE satellite orbiting Earth {5}.

II. CALIFORNIA’S GROUNDWATER

One pressing issue within California is the amount of available groundwater to support a growing population, agriculture and industry [1]. Groundwater, one of the most important resources on Earth, is difficult to understand and manage. The majority of groundwater in California is found in the Central Valley aquifer system (Figure 1). This aquifer supports agriculture that supplies nearly 7 percent of the United States’ food supply, with an estimated annual value of $21 billion [2]. California does not regulate groundwater pumping at the state level. Groundwater management is implemented at the local level and includes groundwater monitoring, basin management and water-use restrictions. Understanding groundwater availability is complicated and may benefit from using different methods. The California Department of Water Resources (DWR) has developed tools to assess groundwater storage changes within the Central Valley with the use of the Central Valley Groundwater-Surface Water Simulation Model (C2VSIM) and a Water Data Library (WDL) Geographic Information System (GIS) change in groundwater storage tool (CGST) [3].

The Central Valley aquifer is a basin, which is about 80 kilometers wide and 650 kilometers long, and is bounded by the coast ranges on the east and the Sierra Nevada on the west (Figure 1). The Central Valley aquifer is contained within three hydrologic regions ‰ÛÒThe Sacramento River Basin, the San Joaquin River Basin, and the Tulare Lake Basin (Figure 1). For this study, the San Joaquin River Basin and the Tulare Lake Basin were combined and will be referred to as the San Joaquin River Basin. The C2VSIM model also incorporates inflow from watersheds within the Sierra Nevada, and although the study areas from C2VSIM and GRACE data are significantly different, the Sierra Nevada is not conducive to groundwater storage, and it is assumed that the majority of changes within the GRACE defined study area are occurring within the Central Valley aquifer [3].

III. METHODS

Image of the area of the modeled GRACE dataset compared to the outline of the Sacramento and San Joaquin hydrologic regions.

Fig. 3. Area of the modeled GRACE dataset compared to the outline of the Sacramento and San Joaquin hydrologic regions.

The GRACE satellite can detect anomalies and deviations from normal values in the Earth’s gravity field [4]. Changes in gravity over short time periods on Earth are often attributed to changes in water, ice and snow. Data were obtained from GRACE and processed to represent changes in total water storage (TWS). This included the removal of other factors such as atmospheric effects, postglacial rebound and large-scale topographic processes, with units expressed as equivalent water thickness [5] [6] [7]. The GRACE sensor operates with two identical satellites that fly in the same orbit (Figure 2). These satellites move further or closer together based on changes in the Earth’s gravitational field. These measurements can then be used with other satellites and ground measurements to calculate how much of that change may be attributed to groundwater. Previous studies from all over the world demonstrate the use of this satellite and how to calculate changes in groundwater storage [8] [9] [10] [11] [12] [13]. For this study, GRACE data were collected from the University of Colorado’s GRACE Data Analysis website [5]. The Sacramento and San Joaquin River Basins were defined using a model [6] (Figure 3).

Other sensors and field measurements can be used to detect surface water storage changes in California, and can aid in calculating changes in groundwater storage that may be due to pumping, drought or natural climate variability. The most important variables to use when converting from changes in TWS to changes in groundwater storage are changes in lakes or reservoirs, soil moisture and snowpack (Figure 4). For this study, changes in reservoir storage were obtained from the California Data Exchange Center (CDEC) [15], soil moisture can be estimated from AMSR-E) [15], and snowpack was obtained from data collected from various sources including satellite data by the National Oceanic and Atmospheric Administration (NOAA) [17]. The monthly anomalies for soil moisture (SMë±), surface water (SWë±) and snow pack (SPë±) were then subtracted from monthly TWSë±,GRACE values to calculate monthly GWë± for the Central Valley, the Sacramento River Basin and the San Joaquin River Basin.

Equation (1)

GWë± = TWSë±,GRACE – ( SWë± + SMë± + SPë± )

Where:

GWë± = groundwater storage anomaly

TWSë± = total water storage anomaly

SWë± = surface water storage anomaly

SMë± = soil moisture storage anomaly

SPë± = snowpack storage anomaly

Image showing examples of each variable of the hydrologic cycle consider in the groundwater calculation

Fig. 4. Visual example of each of the variables of the hydrologic cycle that were considered in the groundwater calculation.

Changes in storage were then calculated for each of the anomalies TWS, SW, SM, SP, and GW over the length of the study. Anomalies were plotted and a trend line was fitted. The slopes of the graphs were then converted into a total volume of water lost or gained over the course of the study.

Changes in groundwater storage in the Central Valley were also calculated using two DWR tools, the C2VSIM model and groundwater level measurements that are incorporated into a GIS CGST. The C2VSIM also calculates changes in groundwater storage using a variety of inputs including precipitation, river discharge and aquifer characteristics. The C2VSIM then calculates the change in groundwater storage for each groundwater basin. Also, the GIS CGST uses groundwater elevation measurements taken each spring to estimate changes in groundwater storage. This tool was used in the Sacramento River Basin only. The groundwater storage estimates from each method were compared.

Graph showing groundwater storage anomalies from GRACE and C2VSIIM

Fig. 5. Groundwater storage anomalies from GRACE and C2VSIM with linear trend lines {19}.

The processing of GRACE data can also have significant impacts on the groundwater storage results [17] [18]. For this study, a Gaussian smoother of 300 kilometers was used to reduce noise to an acceptable level without significantly altering the signal. GRACE also has two types of error, a measurement error and a leakage error, and for this study, the combined measurement and leakage error was found to be 1.01 kilometers3, 1.00 kilometer3, and 1.00 kilometer3 for the Sacramento River Basin, and the San Joaquin River Basin, and the Central Valley, respectively. It should be noted that this was a best-estimate approximation of the errors, and recent studies suggest the errors for smaller regions are correlated and may be larger than currently reported errors [19]. We also used an error of 15 percent to estimate the error involved in the hydrologic variables of SM, SP, and SW, and for the C2VSIM and the GIS CST.

VI. RESULTS

GRACE and C2VSIM groundwater storage anomalies exhibited similar trends for the Central Valley region from October 2002 through September 2009 (Figure 5) [20]. This finding is important as it validates the usefulness of GRACE at scales that are less than or equal to 150,000 kilometers2. It should be noted that GRACE and C2VSIM groundwater anomalies display marked differences during the seasonal peaks and troughs. The GRACE data appears to be more variable than C2VSIM, although the trends for the entire time period are similar. Also, climate variability is observed, and the drought period beginning in 2007 has a distinct negative trend in available groundwater.

Water storage anomalies produced from GRACE were also examined to understand climatic trends and long term differences in each basin. The San Joaquin River Basin exhibited the largest loss in TWS (-21.92 å± 1.92 kilometers3), snowpack (-3.90 å± 0.59 kilometers3) and groundwater storage anomalies (-16.43 å± 2.04 kilometers3) using the GRACE method. In contrast, the Sacramento River Basin displayed the largest loss in surface water storage (-4.22 å± 0.63 kilometers3) compared to the San Joaquin (-2.40 å± 0.36 kilometers3) [20]. Soil moisture storage values, however, remained relatively unchanged throughout the study period in both regions. The results presented here are consistent with the results presented by [1] and [13].

Graph showing changes in groundwater for the Sacramento River Basin, san Jouquin river basin, and the central valley aquifer

Fig. 6. Changes in groundwater for the Sacramento River Basin, the San Joaquin River Basin, and the Central Valley aquifer for GRACE and C2VSIM {19}.

The change in groundwater storage estimates from GRACE for the Sacramento River Basin and the San Joaquin River Basin were not comparable with those from C2VSIM, thus illustrating the usefulness of GRACE on very large scales rather than smaller basins [20]. The GIS CGST estimated a loss of groundwater storage for the Sacramento River Basin of -0.67 å± 0.1 kilometers3 while C2VSIM and GRACE showed losses of -7.70 å± 1.49 kilometers3 and -2.55 å± 0.38 kilometers3, respectively [20]. The apparent differences in the results by the two methods will need to be investigated in the future.

Groundwater storage loss was calculated for the entire Central Valley, the Sacramento River Basin and the San Joaquin River Basin which includes the Tulare Lake Basin using GRACE and C2VSIM. The calculations of changes in groundwater storage for each river basin using both the GRACE-derived and the C2VSIM estimates did not produce comparable results. Although there are large differences on the hydrologic region level, both tools were similar for the entire Central Valley. For the study period, GRACE calculated a total loss of groundwater storage of -14.47 å± 1.49 kilometers3 or -11.7 å± 1.2 million acre-ft, and C2VSIM calculated a total loss of -15.01 å± 2.25 kilometers3 or -12.2 å± 1.8 million acre-ft, for the Central Valley (Figure 6) [20]. For a detailed analysis of results refer to Kuss et al., 2011 [20]. The large losses in groundwater storage in the Central Valley during the study period most likely resulted from a period of extended drought from 2007 to 2009.

V. CONCLUSIONS AND IMPLICATIONS FOR GROUNDWATER MANAGEMENT

The usefulness of the GRACE satellite to monitor changes in groundwater storage for the Central Valley aquifer was examined through a comparison approach. Satellites have been shown to be useful to detect changes in groundwater storage over large areas, and could be an important tool for groundwater management on large scales. The GRACE satellite and other sensors can supplement current groundwater management techniques for the Central Valley aquifer. This study also supports DWR’s modeling efforts for the Central Valley, and illustrates the usefulness of multiple techniques. Current spatial downscaling limitations of GRACE may result in decreased usefulness for smaller scale basin management. Estimated changes in groundwater storage by the C2VSIM and GRACE are comparable, while the estimated changes in groundwater storage for the Sacramento River Basin by the GIS CGST are significantly less than the estimates by the other two methods. Understanding how groundwater storage is changing in the Central Valley is an important aspect of sustainable groundwater resource management in California. This understanding will enhance the ability to make informed decisions on proper management techniques such as curbing water use, implementing artificial recharge and maximizing conservation and water-use efficiency efforts that will aid in sustainable groundwater practices.

REFERENCES

[1] Faunt, C. C., R.T. Hanson, K. Belitz, W. Schmid, S.P. Predmore, D.L. Rewis, K. McPhearson, (2009): Groundwater availability of the Central Valley aquifer, California, US Geological Survey, pp. 246.

[2] California Department of Food and Agriculture. (2010): California Agricultural Production Statistics. Sacramento, CA: California Department of Food and Agriculture.

[3] Department of Water Resources. 2003, California’s Groundwater Bulletin 118, Update 2003. Sacramento, CA. pp. 265. [4] Brush, CF, EC Dogrul, MR Moncrief, J Galef, S Shultz, M Tonkin, D Wendell, TN Kadir, and FI Chung. (2008): Estimating hydrologic flow components of the Central Valley hydrologic flow system with the California Central Valley Groundwater-Surface Water Model. In CF Brush and NL Miller, eds. Proceedings of the California Central Valley Groundwater Modeling Workshop, July 10-11, 2008, Lawrence Berkeley National Laboratory, Berkeley, California. California Water and Environmental Modeling Forum, Sacramento, CA.

[4] NASA. (2002): Studying the earth’s gravity from space: The Gravity Recovery and Climate Experiment (GRACE), National Aeronautics and Space Administration, Greenbelt, MD.

[5] CU, (2011): Real-Time GRACE Data Analysis Site, 2002 to 2009, University of Colorado Cooperative Institute for Research in Environmental Sciences. http://geoid.colorado.edu/grace/grace.php, Boulder, CO.

[6] Oki, T., and Y. C. Sud. (1998): Design of Total Runoff Integrating Pathways (TRIP)—A Global River Channel Network, Earth Interact., 2 (1), 1-37, doi:10.1175/1087-3562(1998)002<0001_DOTRIP>2.3.CO;2.

[7] Swenson, S., and J. Wahr (2006): Post-processing removal of correlated errors in GRACE data, Geophys. Res. Lett, 33(4).

[8] Rodell, M., and J. S. Famiglietti. (2002): The potential for satellite-based monitoring of groundwater storage changes using GRACE: the High Plains aquifer, Central US, Journal of Hydrology, 263(1-4), 245‰ÛÒ256.

[9] Han, S. C., C. K. Shum, C. Jekeli, and D. Alsdorf. (2005): Improved estimation of terrestrial water storage changes from GRACE, Geophys. Res. Lett, 32(7).

[10] Yirdaw, S. Z., K. R. Snelgrove, and C. O. Agboma. (2008): GRACE satellite observations of terrestrial moisture changes for drought characterization in the Canadian Prairie, Journal of Hydrology, 356(1-2), 84‰ÛÒ92.

[11] Zaitchik, B. F., M. Rodell, and R. H. Reichle. (2008): Assimilation of GRACE Terrestrial Water Storage Data into a Land Surface Model: Results for the Mississippi River Basin, J. Hydrometeor, 9(3), 535-548, doi:10.1175/2007JHM951.1.

[12] Leblanc, M. J., P. Tregoning, G. Ramillien, S. O. Tweed, and A. Fakes. (2009): Basin-scale, integrated observations of the early 21st century multiyear drought in southeast Australia, Water resources research, 45(1).

[13] Famiglietti, J. S., M. Lo, S. L. Ho, J. Bethune, K. J. Anderson, T. H. Syed, S. C. Swenson, C. R. de Linage, and M. Rodell. (2011): Satellites measure recent rates of groundwater depletion in California’s Central Valley, 38, L03403, doi:10.1029/2010GL046442.

[14] Oki, T., and Y. C. Sud. (1998): Design of Total Runoff Integrating Pathways (TRIP)—A Global River Channel Network, Earth Interact., 2 (1), 1-37, doi:10.1175/1087-3562(1998)002<0001_DOTRIP>2.3.CO;2.

[15] CDEC, (2011): Reservoir Data Summary, 1964 to 2011, Department of Water Resources California Data Exchange Center (CDEC), http://cdec.water.ca.gov/reservoir.html, Sacramento, CA, (Updated daily).

[16] Njoku, Eni. (2009): updated daily. AMSR-E/Aqua L2B Surface Soil Moisture, Ancillary Parms, & QC EASE-Grids V002, October 2002‰ÛÒSeptember 2009, Boulder, Colorado USA: National Snow and Ice Data Center. Digital media.

[17] NOAA, (2011): National Snow Data, 2002 to 2011, National Operational Hydrologic Remote Sensing Center, http://www.nohrsc.nws.gov/, Chanhassen, MN, (Updated daily).

[18] Swenson, S., and J. Wahr (2006): Post-processing removal of correlated errors in GRACE data, Geophys. Res. Lett, 33(4).

[19] Landerer, F.W. and Swensen, S.C., In Press. Accuracy of scaled GRACE terrestrial water storage estimates, Water Resources Research.

[20] Kuss, A.M., W.T. Brandt, J. Randall, B. Floyd, A. Bourai, M.E. Newcomer, C. Schmidt, J.W. Skiles. In Press. Comparison of changes in groundwater storage using GRACE data and a hydrological model in California’s Central Valley. Proceedings of the American Society for Photogrammetry and Remote Sensing (ASPRS).