Long-term groundwater availability in the Waihe‘e, ‘Īao, and Waikapū aquifer systems, Maui, Hawai‘i

Kolja Rotzoll; Delwyn S. Oki; Adam G. Johnson; William R. Souza

doi:10.3133/sir20215113

Home/Archive/USGS

USGS · sir20215113

Long-term groundwater availability in the Waihe‘e, ‘Īao, and Waikapū aquifer systems, Maui, Hawai‘i

Kolja Rotzoll·Delwyn S. Oki·Adam G. Johnson·William R. Souza

View original source ↗Read source document ↗Export metadata ↓

Abstract

Groundwater levels have declined since the 1940s in the Wailuku area of central Maui, Hawai‘i, on the eastern flank of West Maui volcano, mainly in response to increased groundwater withdrawals. Available data since the 1980s also indicate a thinning of the freshwater lens and an increase in chloride concentrations of pumped water from production wells. These trends, combined with projected increases in demand for groundwater in central Maui, have led to concerns over groundwater availability and have highlighted a need to improve general understanding of the hydrologic effects of proposed groundwater withdrawals in the Waihe‘e, ‘Īao, and Waikapū areas of central Maui. A numerical groundwater model was constructed to simulate the flow and salinity of groundwater in central Maui. The model simulates the effects of changes in groundwater withdrawals and recharge on water levels, freshwater-lens thicknesses, and chloride concentrations of pumped water from production wells. The model incorporates updated water-budget estimates of groundwater recharge from infiltration and direct recharge, seepage in stream channels, and inflow from inland areas. Mean annual groundwater recharge from infiltration and direct recharge was estimated using a daily water-budget model and the most current data, including the distributions of monthly rainfall and potential evapotranspiration, for the study area for nine historical periods from 1926 through 2012: 1926–69, 1970–79, 1980–84, 1985–89, 1990–94, 1995–99, 2000–04, 2005–09, and 2010–12. The water-budget model also estimated groundwater recharge based on one hypothetical scenario that used 1980–2010 rainfall and 2017 land cover. For the nine historical periods, estimated recharge from infiltration and direct recharge within the area of the groundwater model ranged from 30.4 million gallons per day (Mgal/d) during 2010–12 to 98.7 Mgal/d during 1926–69. Variability in recharge during these periods mainly reflects changes in rainfall and irrigation over time. Between 2010 and 2014, streamflow restoration in previously diverted streams resulted in an estimated increase in recharge from seepage in stream channels of about 12.5 Mgal/d. Average groundwater inflow of about 39.6 Mgal/d from inland, dike-intruded areas to the main area of interest was estimated from an existing island-wide numerical groundwater-flow model, which is at a larger scale and incorporates a greater number of simplifying assumptions. The numerical groundwater model developed for this study was calibrated to 1926–2012 transient water levels, vertical salinity profiles, and chloride concentrations of water pumped by production wells in the study area. The model was then used to evaluate one future recharge and six selected withdrawal scenarios, developed in consultation with the Maui Department of Water Supply, in terms of long-term changes in water level and 50-percent ocean-water salinity surface. The groundwater model was also used to simulate the future salinity of water withdrawn by existing and proposed production wells. The simulations were run to steady-state conditions, providing an estimate of the long-term effects of changes in withdrawal and recharge on the groundwater resource. Results of the simulated future withdrawal scenarios indicate that, relative to 2017–18 rates, the scenarios’ long-term effect of increased withdrawals ultimately leads to lower water levels and a higher 50-percent ocean-water salinity surface indicating a thinning of the freshwater lens. Results also indicate that the increased withdrawals produce some groundwater with chloride concentration below 250 milligrams per liter and some groundwater with higher chloride concentration. The amount of drawdown near production wells and the quality of water withdrawn from production wells is dependent on the rate and spatial distribution of the withdrawals. The model was also used to evaluate how groundwater availability may be affected for a drier recharge scenario based on a published study of future climate. Model results of the future recharge scenario indicate that the rate of groundwater recharge is a controlling factor for (1) water levels, (2) the 50-percent ocean-water salinity surface, and (3) the quality of water withdrawn from production wells in the Wailuku area. Coupled with reduced groundwater recharge (with all other factors remaining equal), the modeled future withdrawals in the scenario would tend to cause lower water levels, a higher 50-percent ocean-water salinity surface, and increased salinity of water withdrawn from production wells. The three-dimensional numerical groundwater model developed for this study utilizes the latest available hydrologic and geologic information and is a useful tool for understanding the long-term hydrologic effects of additional groundwater withdrawals in central Maui. The model has several limitations, including its non-uniqueness and inability to account for local-scale heterogeneities. Short-term effects of changes in recharge and withdrawals—and optimization of pumping rates to meet increased demand for water with acceptable salinity—are possible conditions for future simulation analyses.

Explore related subjects

Scientific Investigations Report Report

Source-reported bounding extent: 20.56593890346526° to 21.056307701901847° latitude; -156.7254638671875° to -155.9564208984375° longitude. This indicates report coverage, not an exact sampling location. View area on OpenStreetMap.

Keep this discovery

Explore connections, maps & timelines

BibTeX RIS

Kolja Rotzoll, Delwyn S. Oki, Adam G. Johnson, William R. Souza. 2022-06-03. Long-term groundwater availability in the Waihe‘e, ‘Īao, and Waikapū aquifer systems, Maui, Hawai‘i. https://doi.org/10.3133/sir20215113

Cite the original work for its findings. Save a collection to share your selection of sources.

KEEP EXPLORING

Related USGS reports

USGS2026

Karst hydrologic processes in carbonate fensters of the western Great Smoky Mountains, Tennessee, 2017–18

The geology of the Great Smoky Mountains in Tennessee and North Carolina is dominated by metamorphic strata. Along the western extent of the range, however, several carbonate fensters expose the Ordovician Jonesboro Limestone along the western edge of the Great Smoky Fault. The fensters (Cades Cove, Tuckaleechee Cove, and Wear Cove) create opportunities for allogenic recharge at the contact between the surrounding, relatively insoluble, strata and the underlying soluble formations, resulting in the development of karst systems. Karst environments have a high degree of connectivity between surface streams and groundwater systems; thus, karst aquifers are inherently vulnerable to contamination from surface activities. Karst environments also provide a specific environment for sensitive biota. Protection of the unique resources of karst environments requires an understanding of the behavior of karst aquifers, knowledge of karst recharge areas, and the mapping of groundwater flow paths. A study began in 2017 to characterize the hydrologic behavior of the karst systems in Great Smoky Mountains National Park through hydrologic monitoring, dye tracing, and seepage investigations. Hydrologic monitoring equipment was installed at several karst hydrogeologic features: two in-cave sites, two major springs, Abrams Creek in Cades Cove, and a karst depression wetland. Dye tracing was conducted to determine resurgences for the cave systems, to delineate spring recharge areas, and to estimate karst groundwater velocities. Seepage investigations along Abrams Creek helped determine if streams were sinking into the Jonesboro Limestone, either at the contact with insoluble strata or as streams flowed across the carbonate fenster. Study results indicate considerable differences in groundwater behavior between the conduit-dominated flow systems in Tuckaleechee Cove compared to the lower-gradient, colluvium-covered karst systems in Cades Cove. Study results also indicate that contaminants in the area of Bull sink and Whiteoak Sink could be rapidly transported through karst groundwater towards springs in Tuckaleechee Cove. Karst resources in Cades Cove are most likely to be affected by local events or activities taking place within the cove itself or in upper Abrams Creek upstream of the Great Smoky Fault. This study provides information on the vulnerability and uniqueness of the karst systems in the western Great Smoky Mountains, which are characterized by complex geology, high-gradient subterranean streams, deep caves, and carbonate fensters.

Tennessee↗

USGS2026

Analysis of groundwater geochemistry in the Lansing, Michigan, area

Groundwater resources in the Lansing, Michigan, area are protected by a wellhead protection plan, but previous studies have indicated that groundwater in this area may be susceptible to surface contamination. The U.S. Geological Survey, in cooperation with the Lansing Board of Water & Light, analyzed geochemical and aquifer characteristics to assess the potential for contamination of groundwater resources. Recharge age categories were determined from tritium concentrations for some wells, and by using ternary diagrams and principal component analysis for other wells. These age classifications, combined with hydrogeologic and geochemical data, were analyzed by using statistical methods to identify trends that may indicate areas of increased groundwater susceptibility. Wells with modern recharge tended to have higher sodium, chloride, and chloride:bromide ratios than did wells with mixed recharge, suggesting greater surface connectivity and potential influence from anthropogenic sources such as road salt. In contrast, wells in the mixed age category exhibited more stable geochemistry and lower constituent concentrations. Classification and regression tree models identified aquifer geometry—particularly the thickness of the uppermost bedrock unit and presence of confining layers—as key predictors of recharge age and geochemical indicators. The trends identified from these analyses could help to reveal locations in the Lansing area that are among the most susceptible to contamination.

Michigan↗

USGS2026

Evaluating the number of realizations (train-test splits) necessary to ensure robust estimates of conventional hydrothermal resource favorability using Monte Carlo cross-validation, the Great Basin, Western USA

Supervised machine learning approaches used to produce favorability maps in resource assessments commonly partition an input dataset into training and testing subsets (that is, through train-test splits). However, if a single split is not representative of underlying resource conditions, the resulting favorability map may be biased, particularly with sparse or heterogeneous datasets. To mitigate the potential for unrepresentative favorability maps, recent works by U.S. Geological Survey researchers have used Monte Carlo cross-validation to implement a large number ( n) of randomized train-test splits to produce n realizations. The realizations are then aggregated using cell-wise median predictions to produce the resulting final favorability map, recognizing that the median will remove extreme values that might have been estimated from a comparatively few unrepresentative realizations. Theoretically, the median prediction at every cell converges to a single estimate as n approaches infinity, and the question addressed herein is: How many realizations are necessary to get a robust estimate (that is, an estimate unaffected by additional realizations)? Therefore, convergence metrics are described to identify the number of train-test splits and realizations needed to make a reliable map using the data and supervised machine learning workflow for conventional hydrothermal resource assessment of the Great Basin of the Western United States. For this study, robust prediction occurs with n approximately equal to 600 realizations, and smaller ensembles (for example, n is approximately equal to 200 or even 100) may be suitable during preliminary or exploratory analyses. Although the resulting threshold numbers will vary depending on the data and modeling approach used, the convergence metrics proposed and applied herein are applicable to other supervised machine learning studies that use Monte Carlo cross-validation.

California, Idaho, Nevada, Oregon, Utah↗