The conjunctive use model consists of three loosely-coupled sub-models: 1) a surface water supply (SWS) model, 2) an unsaturated zone water budget (UZWB) model, and 3) a groundwater flow model. The base period of the study covers the fiscal water years of 1970-99. The purpose of the SWS model is to calculate the surface water balance for the source and diversion channels in the inter-district channel network. For each modeled surface water channel, the SWS model computes surface water deliveries from it to each district and conveyance losses from it due to evaporation and channel seepage. The primary model outputs are monthly surface water deliveries to each district and monthly seepage rates from modeled channels. The surface water deliveries became input for the UZWB model. The channel seepage became input for the groundwater flow model as localized aquifer recharge. The allocation of surface water within each district, via the implicitly modeled intra-district surface water distribution system, is estimated by the UZWB model.
The total imported surface water for 1970-99 from the CVP and the Success Reservoir are 13,329,262 and 4,653,501 acre-feet (af), respectively. The SWP and the Kings River imported the lesser amounts of 88,625 and 7,332 af, respec- tively. Annual CVP diversions varied from 125,970 af in 1977 to 679,298 af in 1993 with a 30-year annual average of 444,309 af. The Tule River and Pioneer Ditch both receive regulated releases from Success Reservoir. Tule River annual imports varied from 11,034 af in 1977 to 607,154 af in 1983 while the Pioneer Ditch varied from 3,445 af in 1973 to 5,874 af in 1990. The total natural runoff from the Deer Creek and White River from 1970-99 were 703,444 and 219,098 af, respectively. Deer Creek runoff varied from 4,082 af in 1992 to 103,716 af in 1983 while the White River runoff varied from 422 af in 1977 to 37,985 af in 1998.
From 1970-99, a total of 15 million af of surface water was applied by the service districts in the study area. The applied surface water varied from a low of 135,482 af in 1977 to a high of 708,293 af in 1996. The Lower Tule River Ir- rigation District and the Delano-Earlimart Irrigation District together account for 59% of the total applied surface water while occupying approximately 40% of the incorporated area in the study area. Over the 30-year base period, an estimated total of 3.5 million af of seepage conveyance loss occurred in all sur- face water channels. Seepage in the Tule River, Deer Creek, and White River accounted for 85% of the total seepage. Total annual seepage varied from a low of 8,128 af in 1977 to 467,084 af in 1983.
The UZWB model then calculates the monthly water storage changes in the soil root zone and deep vadose zone of each land unit, where the land unit is the UZWB model scale of resolution. It also models the intra-district surface water distribution system by estimating the monthly allocation of surface water to individual land units within each district. The main model outputs were the recharge to the unconfined aquifer from surface applied water and precipitation, and the groundwater pumping demand from the unconfined and confined aquifers. The recharge and groundwater pumping rates became input for the groundwater flow model.
The total annual agricultural and urban consumptive use ranged from 865,800 af in 1970 to 1,246,700 af in 1999. The estimated total pumping ranged from 148,100 af in 1978 to 570,000 af in 1990. As expected, pumping was heaviest during the droughts of 1975-77 and 1987-92, and lightest during the wet years of 1973, 1978, 1982-83, 1995, and 1998. Precipitation totals varied from 177,800 af in 1990 to 974,400 af in 1998. Diffuse recharge from surface applied water ranged from 64,800 af in 1992 to 350,100 af in 1983.
The net aquifer recharge for the entire study area was computed by ag- gregating the aquifer recharge and groundwater pumping of each land unit to this scale and adding the contribution to aquifer recharge from channel seepage. The monthly net recharge was then summed to produce a cumulative annual net recharge from 1970 to each fiscal water year from 1971-99. The water balance computed for the entire study area neglects horizontal groundwater inflows and outflows through its vertical boundaries. Groundwater fluxes undoubtedly exist along these boundaries. However, net fluxes are likely small in comparison to the total changes in storage due to vertical stresses applied to the entire study area (e.g. groundwater pumping, evapotranspiration, applied surface water, channel seepage). Horizontal groundwater flow on the inter-land unit and inter-district scales is expected to be more significant. For computing a total water balance, however, we made the simplifying assumption that the study area behaves as a relatively closed system where the net horizontal groundwater inflows through its vertical boundaries are small. Invoking this assumption, we then use the cumulative net recharge as an estimate of the cumulative groundwater storage change in the aquifer system. Ideally, verification of these estimates is performed by comparing them with an objective measure of the study area aquifer storage changes. However, changes in groundwater storage are not directly observable and must always be estimated using non-direct measures. As such, an objec- tive measure for verification does not exist. As an alternative, we compare the water balance model results with those produced by the water-table fluctuation (WTF) method.
The trends in cumulative annual groundwater storage changes computed from the water balance and the WTF method from 1970-99 were quite similar. The minimum and maximum differences between them were 2,450 af (1980) and 752,387 af (1991), respectively. From 1970, the maximum amount of ground- water accumulation occurred in the spring of 1987 with the WTF method and the water balance estimating positive storage changes of 1,146,286 and 898,128 af, respectively. The maximum groundwater overdraft occurred in 1993 with the WTF method and the water balance estimating negative storage changes of 1,610,210 and 1,218,566 af, respectively. The 1987 and 1993 fiscal water years marked the beginning and ending of a major 6-year drought in California, respectively.
Finally, the groundwater flow model calculates the changes in water levels in the aquifer system subject to transient groundwater recharge and pump- ing stresses. A post-processing routine calculates the cumulative groundwater storage changes over each district and the entire study area for each stress pe- riod. An automated calibration of the groundwater flow model was performed to refine the conceptual model of the hydrogeology and to estimate the spa- tial distributions of the aquifer system horizontal hydraulic conductivity. The calibration period of the groundwater flow model is 1970-85 and the validation period is 1986-99.
Three different conceptual models of the aquifer system horizontal hydraulic conductivity, Kh, structure were evaluated in the calibration process: 1) Khas an exponential function of the specific yield, Sy, distribution, 2) Kh as a linear function of the saturated hydraulic conductivity of the soil survey map- ping units, and 3) division of the model domain into square zones of uniform size. The models were calibrated against both spatially distributed hydraulic head targets and cumulative groundwater storage change targets for seven of the largest districts. The discretization of the model domain into uniform square zones provided the most robust Kh structure and produced the most reason- able estimates of hydraulic head and district groundwater storage changes from the three conceptual models over the 1971-85 calibration period. The calibrated model was then used to compute the annual net inter-district groundwater fluxes between adjacent districts. In general, groundwater flux directions were con- sistent with the large-scale hydraulic gradients. Annual inter-district net fluxes between adjacent districts ranged from negligibly small ( < 100 af) to as much as 80,000 af (e.g. net flux from Lower Tule River ID to Pixley ID). Net inter-district fluxes were generally a function of the local transmissivity, the length of the shared border between adjacent districts, and the differences in their surface water supplies.