Geology Reports⌕ Search

SEARCH · Geology Reports

Results for “PLoS ONE”

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.

1,186 records · Page 6Linked to original sources
USGS2026

Controls on natural hydrogen generation during serpentinization of mantle rocks

Mantle rocks undergoing serpentinization can generate significant amounts of natural hydrogen, yet the rates and controlling processes remain poorly understood. Here, we constrain the possible hydrogen generation rates in two distinct mantle rock types, the fertile lherzolites of the Western Pyrenees and the depleted harzburgites of Northern California, to relatively low rates of ~0.1 to ~0.5 tonnes H₂ yr⁻¹ km⁻³ of reactive rock. When integrated over the full reactive volumes, this corresponds to total production rates of ~300 to ~600 tonnes H₂ yr⁻¹. By combining three-dimensional geophysical inversion with numerical modelling of fluid-rock processes, we show that hydrogen generation rates are mainly limited by H₂ saturation in the fluid and reaction kinetics. Under these constraints, hydrogen generation in mantle-derived serpentinization systems proceeds slowly, making rapid large-scale replenishment unlikely and suggesting that large, economically relevant accumulations, would require timescales of thousands to tens of thousands of years to develop.

northern California, western Pyrenees↗
USGS2026

On-demand global Landsat evapotranspiration product: Development, evaluation, and dissemination

Global actual evapotranspiration (ET) is one of the essential climate variables needed to understand and manage the relationships among food, energy, and water resources. The U.S. Geological Survey Earth Resources Observation and Science (EROS) Center launched a provisional ET product in 2020, offering on-demand, field-scale global coverage derived from Landsat data through the EROS Science Processing Architecture (ESPA) platform. The ESPA interface provides ET data for cloud-free Landsat overpasses starting in 1982 with Landsat 4 through the current Landsat 9. The ET data are delivered as a Provisional Level-3 Science product created using the Operational Simplified Surface Energy Balance (SSEBop) model. Landsat surface temperature and reference ET are the main model drivers along with vegetation index and net radiation for model parameterization. A large volume of Landsat-based ET orders (e.g., over 1,200,000 images from June 2020 through December 2025) around the world indicate increasing awareness and application of the ET data. The ESPA platform enables land and water resource managers and researchers to access a first-order ET product without requiring advanced knowledge of remote sensing technology or evapotranspiration modeling. We present the methodology and workflow of the on-demand Landsat ET product and its performance evaluations over diverse hydro-climatic settings. The product can help estimate field-scale consumptive water use and thus quickly and consistently assess historical water use, allocation, and budget to inform water management under changing environments. Future ET data aggregated to monthly and seasonal time scales are expected to enhance integration with decision-making tools and procedures.

Remote Sensing of Environment↗
USGS2026

Evaluating the effects of a multi-modal deterrent on upstream passage by invasive carp and native species

Invasive carp, including Silver Carp ( Hypophthalmichthys molitrix ) and Grass Carp ( Ctenopharyngodon idella ), are increasing their range within the Mississippi River Basin, threatening native fishes and ecosystem function. We conducted a 3-year field test of a multi-modal deterrent that uses sound, bubbles, and light to test its ability to impede upstream passage by invasive carp. The test was performed at a lock and dam on the Cumberland River, Kentucky, USA, using acoustic telemetry to track the movements of Silver Carp, Grass Carp, and several native fish species. When the deterrent was operating, the estimated probability of upstream lock passage was 53% lower (95% CI: 41%–63%) for Silver Carp but 74% higher (95% CI: 3%–190%) for Grass Carp compared to times when the deterrent was not operating. Paddlefish ( Polyodon spathula ) passage was reduced by 43% (95% CI: 71% reduction to 9% increase), although this was not statistically significant. We detected no significant effects on Freshwater Drum ( Aplodinotus grunniens ) or Smallmouth Buffalo ( Ictiobus bubalus ). However, sample sizes were small for these native species, limiting our ability to identify effects. Water temperature, tailwater elevation, and vessel lockages affected passage rates, although this varied among species.

Kentucky↗
USGS2024

The impact of source time function complexity on stress drop estimates

Earthquake stress drop—a key parameter for describing the energetics of earthquake rupture—can be estimated in several different, but theoretically equivalent, ways. However, independent estimates for the same earthquakes sometimes differ significantly. We find that earthquake source complexity plays a significant role in why theoretically (for simple rupture models) equivalent methods produce different estimates. We apply time‐ and frequency‐domain methods to estimate stress drops for real earthquakes in the SCARDEC (Seismic source ChAracteristics Retrieved from DEConvolving teleseismic body waves, Vallée and Douet, 2016 ) source time function (STF) database and analyze how rupture complexity drives stress‐drop estimate discrepancies. Specifically, we identify two complexity metrics—Brune relative energy (BRE) and spectral decay—that parameterize an earthquake’s complexity relative to the standard Brune model and strongly correlate with the estimate discrepancies. We find that the observed systematic magnitude–stress‐drop trends may reflect underlying changes in STF complexity, not necessarily trends in actual stress drop. Both the decay and BRE parameters vary systematically with magnitude, but whether this magnitude–complexity relationship is real remains unresolved.

Bulletin of the Seismological Society of America↗
USGS2026

Linking distribution and return-on-investment models to optimize woody management for prairie grouse in Nebraska

Grasslands in Nebraska, USA, face threats from agricultural conversion, urban development, and woody encroachment, all of which negatively affect prairie grouse ( Tympanuchus spp.) populations. To optimize conservation planning, Nebraska wildlife agencies developed probabilistic area-based surveys for greater prairie-chicken ( T. cupido ) and sharp-tailed grouse ( T. phasianellus ) to sample landscapes across a range of environmental conditions. This design improves historical surveys and enables the development of distribution models that quantitatively define habitat associations and support scenario-based conservation planning. Using survey data collected during 2020–2022, we modeled prairie grouse occurrence and abundance as functions of land cover, topography, and climate using Bayesian logistic and zero-inflated negative binomial models with regularized horseshoe priors. We then conducted a maximum potential return-on-investment analysis of woody cover treatments, assuming sustained treatment success, relative to projected impacts of woody encroachment on prairie grouse populations by 2050. Among modeled associations were a positive association with grasslands having low woody cover and a negative association with grasslands having high woody cover. Across the 3-year period, median estimated annual populations were 142,380 for greater prairie-chicken (range of 95% CIs across years = 68,821–277,615) and 64,154 for sharp-tailed grouse (range of 95% CIs across years = 27,550–146,559). Under projected woody encroachment, mean predicted population declines were 10% for greater prairie-chicken (range of 95% CIs = 7–14%) and 6% for sharp-tailed grouse (range of 95% CIs = 5–7%). Areas with high prairie grouse density and low treatment costs in 2021, and high projected woody encroachment and population loss by 2050, offered the greatest return on investment for woody management. Return on investment was greatest in the northwestern Shortgrass Prairie ecoregion (northwestern Nebraska) for sharp-tailed grouse and the eastern Sandhills ecoregion (central Nebraska) for both species. These models underscore the value of evidence-based, quantitative approaches for prioritizing conservation actions on working lands. Scenario-based modeling could be extended to guide other treatments, such as optimizing restoration (e.g., Conservation Reserve Program) or incentivizing grassland persistence in areas with predicted climate resilience.

Nebraska↗
USGS1999

Effects of snowmobile use on snowpack chemistry in Yellowstone National Park, 1998

Snowmobile use in Yellowstone National Park has increased substantially in the past three decades. In areas of greatest snowmobile use, elevated levels of by-products of gasoline combustion such as ammonium and benzene have been detected in snowpack samples. Annual snowpacks and snow-covered roadways trap deposition from local and regional atmospheric emissions. Snowpack samples representing most of the winter precipitation were collected at about the time of maximum annual snow accumulation at a variety of locations in the park to observe effects of a range of snowmobile traffic levels. Concentrations of organic and inorganic compounds in snow samples from pairs of sites located directly in and off snow-packed roadways used by snowmobiles were compared to concentrations in samples collected at nearby off-road sites. Concentrations of ammonium were 2 to 5 times higher for the in-road snow compared to off-road snow for each pair of sites. Thus, concentrations decreased rapidly with distance from roadways. In addition, concentrations of ammonium, nitrate, sulfate, benzene, and toluene in snow were positively correlated with snowmobile use.

Wyoming↗
USGS2026

Simulation of the impacts of spring diversions on streamflow in the Strawberry Creek watershed, San Bernardino County, California, using an integrated hydrological model

The Strawberry Creek watershed, in the San Bernardino Mountains of southern California, features a group of natural springs known as Arrowhead Springs that have been augmented with diversions in the form of sub-horizontal borings and tunnels. Understanding the effects of these structures on streamflow through groundwater capture is crucial for managing surface-water resources in this watershed. In this study, we constructed the Strawberry Creek integrated hydrological model to increase this understanding. The Strawberry Creek integrated hydrological model is an integrated surface runoff and groundwater model, GSFLOW to simulate surface runoff and infiltration and groundwater flow. The model has 263 rows and 176 columns in each of its three layers. The model area includes the Strawberry Creek and four adjacent watersheds. The precipitation-runoff modeling system was calibrated using two streamflow gaging stations, and the groundwater and surface-water flow model was calibrated to the reported spring diversion discharge and a sparse number of groundwater-level measurements. The Strawberry Creek integrated hydrological model was run with and without diversions active, and simulated streamflow was compared, indicating that in the headwaters of Strawberry Creek, about 35 percent of the diversion flow was captured from base flow.

California↗
USGS2025

Rebuilding a volcano one lava flow at a time—Visualizing the lava dome-building eruption in the crater of Mount St. Helens, 1982–1986

Between 1980 and 1986, the U.S. Geological Survey made a series of 1:2,000-scale topographic contour maps from aerial photographic surveys to monitor the eruption. These maps were made for operational purposes and were not intended for publication. Since then, advances in technology made it possible to digitize the original, highly detailed hardcopy maps and derive new digital data elevation models of the surface of the lava dome. These digital elevation models allow for the visualization of the progression of the eruption and reveal the rubbly, chaotic surface of the lava flows and dome. Additionally, these new data help fill gaps in the long-term record of topographic changes that have occurred at the volcano since the May 18, 1980, eruption.

Washington↗
USGS2024

Simulated effects of projected 2014–40 withdrawals on groundwater flow and water levels in the New Jersey Coastal Plain

Abstract Groundwater flow between 2014 through 2040 was simulated in the New Jersey Coastal Plain based on three withdrawal scenarios. Two of the scenarios were based on projected population trends and the assumption of water conservation; the nominal water-loss scenario projected a status quo in the efficiency of water loss in the delivery systems whereas the optimal water-loss scenario projected a better water-loss efficiency resulting in less withdrawals. The third scenario assumes that all wells will withdraw water at their full allocation level which is generally much more than reported withdrawals in 2013 or projected under the other two scenarios. Maps and summaries of heads and drawdowns are presented for nine confined aquifers. All the aquifers have areas with heads below sea level by 2040. Of the three scenarios, the drawdowns are most extreme in the full allocation scenarios; there are large areas of head decline greater than 20 feet in 5 of the 9 confined aquifers. The exceptions are the Vincentown aquifer, despite some areas of large drawdown in the vicinity of wells, and the three Potomac-Raritan-Magothy (PRM) aquifers where withdrawals are regulated by Critical Area restrictions. The nominal and optimal water-loss scenarios have some areas of head declines; most are less than 15 feet. The simulation of these scenarios shows some extensive areas of head recovery as well—especially in the aquifers that are regulated by the Critical Area restrictions. Budgets of inflow and outflow components were calculated for 44 hydrologic budget areas (HBAs). The budget analysis shows that the water movement is complex and varies based on the aquifer geometry and location of pumping wells. Flow components between the unconfined and confined parts of the system were summarized by HUC11 (hydrologic unit code 11) basins.

New Jersey↗
USGS2025

Revisiting an enigma on California's north coast: The Mw6.5 Fickle Hill earthquake of 21 December 1954

Many earthquakes occur along the North Coast of California in the vicinity of the Mendocino Triple Junction (MTJ), where the Pacific, Gorda, and North American (NA) plates meet, and on the adjacent plate boundaries. The MTJ marks the nexus of the Mendocino and San Andreas faults with the Cascadia subduction zone (CSZ). Historically, most large earthquakes around the MTJ have been within the offshore Gorda plate and its subducted portion beneath the NA plate. North of the MTJ, active faults mapped in the NA plate are part of the CSZ fold‐and‐thrust belt. Although some events have been detected in the NA plate, no large historic events have been associated with mapped surface faults. The 21 December 1954 M w 6.5 earthquake in Humboldt County is one possible exception. Using published data from catalogs and articles, unpublished data from Berkeley’s archives, and S ‐ P times interpreted from two U.S. Coast and Geodetic Survey (USCGS) accelerometers, we determine a probability cloud for the earthquake’s hypocenter using NonLinLoc. The highest probability location lies beneath Fickle Hill just east of the city of Arcata, California, at 40.87° N, 124.03° W, and ∼11 km depth. Using P ‐wave polarities from Berkeley stations and the digitized waveforms from the accelerometers, we find that the focal mechanism most consistent with the data indicates thrust movement with strike, dip, and rake of 350°, 10°, and 90°, respectively, at a depth of 14 km. Given the depth uncertainties of both this event and the megathrust, this implies that the earthquake most likely took place on the subduction interface rather than on the mapped faults in the Mad River fault zone that trend 322° and dip to the northeast. The revisited intensity in the epicentral region also supports a location beneath Fickle Hill to the east of the city of Arcata, California.

California, Idaho, Nevada, Oregon, Washington↗
USGS2025

To heal or not to heal?: 1. The effect of pore fluid pressure on the frictional healing behavior of lithologies in Oklahoma

The competition between fault healing (i.e., re-strengthening) and fault loading determines the timing and magnitude of fault failure within the seismic cycle. Repeating earthquakes can give observational estimates of fault healing rates, however, it is difficult to link laboratory studies of frictional healing and observed healing rates from repeating earthquakes in part because of uncertainty in lithology at depth. Due to well-constrained and relatively simple geology, earthquakes in Oklahoma can be linked to the granitic basement rock and to the Arbuckle Group, which is primarily composed of dolomite at earthquake depths. Here, we conduct friction experiments to measure healing rates of the two earthquake-bearing lithologies at confining pressures representative of earthquake depths and pore pressures ranging from 0% to 80% of the confining pressure. We measure frictional healing by executing slide-hold-slide tests with hold times ranging from 3 s to 3000 s. The friction experiments on the Troy Granite indicate that pore fluid pressure does not greatly affect healing rate. On the other hand, the dolomite of the Arbuckle Group exhibits decreased healing with increased pore fluid pressure, with weakening at the highest pore pressure. We hypothesize that this is due to an increase in dissolution of dolomite at high pore pressures/low effective normal stress. These healing rates are used in the companion paper to understand the moment-recurrence time behavior of repeating earthquakes in Prague, Oklahoma. This work has implications for possible enhanced dissolution and weakening behavior of the Arbuckle Group during wastewater injection activities.

Oklahoma↗
USGS2026

Stress states on the eve of past earthquakes inform earthquake rupture through fault complexity along the San Andreas and San Jacinto faults

Estimating the evolving state of stress along active fault systems can provide insight into the conditions that generated past ground‐rupturing earthquakes and influenced their ability to propagate through areas of geometric complexity, such as fault branches and stepovers. We use quasi‐static forward numerical models that incorporate the 3D complex configuration of active faults in southern California to estimate shear tractions on the geometrically complex southern San Andreas and San Jacinto faults from 1000 to 1900 C.E. These tractions include interseismic accumulation of traction due to tectonic loading, viscoelastic relaxation of shear stress within the upper crust between earthquakes, and effects of other earthquakes on the fault network. We simulate ground‐rupturing earthquakes based on the along‐strike earthquake extents modeled by Scharer and Yule (2020) , assuming that stress drop is complete in each earthquake. We use Monte Carlo simulations to estimate uncertainty in evolving shear tractions due to uncertainties in earthquake timing and in upper‐crustal viscosity. Pre‐earthquake shear tractions typically do not exceed ∼2 MPa. Although ruptures with length <200 km have pre‐earthquake shear tractions that range from near zero to ∼1.75 MPa, these tractions are not less than ∼0.4 MPa for earthquakes with rupture length >200 km. Earthquakes with long (>200 km) ruptures occur only in the single‐stranded part of the system, whereas those with short (<125 km) rupture length and high pre‐earthquake shear traction occur near fault stepovers and branches. This suggests that high accumulated shear traction encourages longer rupture propagation, but may not be sufficient to overcome geometric complexities. This modeling approach informs our understanding of rupture propagation and provides estimates of fault shear tractions that are unavailable from direct measurements.

Callifornia↗
USGS2026

Northern bobwhites select for native grasses on working grazing land

Northern bobwhite ( Colinus virginianus ; bobwhite) populations have experienced an 85% decline across most of its range since the 1960s. The most drastic decreases have been in the southeastern United States where biologists attribute the decline to widespread habitat loss, including conversion of native grasslands to tall fescue ( Lolium arundinaceum ) and other exotic cool-season grass (CSG) pasture and hayfields. Including agricultural lands in conservation programs could improve habitat conditions on a regional scale. One working-lands conservation strategy involves the use of native warm-season grasses (NWSGs) rather than cool-season exotic grasses as economically viable cattle forage. To evaluate this management practice for creating bobwhite habitat, we conducted a field study on a bobwhite population in eastern Kentucky, USA, where grazed NWSG and burned NWSG fields were adjacent to exotic grazed and exotic hayed CSG fields. Between April 2019 and October 2022, we used radio telemetry to track bobwhites to evaluate resource selection at second- and third-order scales. Bobwhites used traditionally managed CSGs less than would be expected by chance at the second-order scale in both breeding and non-breeding seasons. Bobwhite use of grazed NWSG and burned NWSG was greater than would be expected by chance at the second-order scale in both breeding and non-breeding seasons. Similarly, at the third-order scale, bobwhites used CSGs less than would be expected by chance, whereas NWSGs and woody stems were used more than by chance alone. Our study suggests that bobwhites select NWSGs at multiple spatial scales despite broad-scale availability of managed exotic CSGs. Under a working-lands framework, the integration of NWSGs into working lands could create bobwhite habitat in the southeastern United States, especially if woody stems are present.

Kentucky↗
USGS2024

Analyzing spatial distributions and alignments of pitted cone features in Utopia Planitia on Mars

Martian geomorphology and surface features provide links to understanding past geologic processes such as fluid movement, local and regional tectonics, and feature formation mechanisms. Pitted cones are common features in the northern plains basins of Mars. They have been proposed to have formed from upwelling volatile-rich fluids, such as magma or water-sediment slurries. In this study, we map the spatial distributions of pitted cone features across the Utopia Planitia (UP) basin. Using the Average Nearest Neighbor technique, we find that pitted cone features appear to be clustered across the basin, occurring in a narrow band around the circumferential rim of UP. Additionally, we find that pitted cone features also appear to be aligned in chains that are sub-parallel with the basin margin. Parallel bands of cones generally follow elevation contours of the UP basin, which suggests elevation or a correlated factor plays a major role in pitted cone and cone chain formation. We propose that pitted cones and cone chains may be related to vertical fractures formed around the UP basin rim from subsidence of infilling basin material.

Icarus↗
USGS2025

Temporal and spatial comparison of coal mine ventilation methane emissions and mitigation quantified using PRISMA satellite data and on-site measurements

Emission monitoring at the facility level (bottom-up, BU) is key for accurate reporting of coal mine methane (CMM) emissions. Recent advances in aerial and satellite observations (top-down, TD) indicate that these methods have the potential to support CMM emissions monitoring and reporting of, as well as track the effectiveness of, mitigation actions. That said, studies have shown discrepancies between BU and TD estimations. Performing TD monitoring with concurrent BU measurements at the same mines may help address the observed discrepancies and improve quantification methods to narrow the gap between BU and TD data. This paper presents a comparison of concurrent BU-TD methane emission fluxes monitored from two ventilation shafts at a southwestern Virginia, USA, longwall mine to complement the existing body of studies on satellite-based monitoring of coal mines by incorporating continuous ground-based monitoring with concurrent TD monitoring of methane emission fluxes to address the gap and provide valuable insights into temporal emission patterns. The shafts were monitored on multiple dates between 2020 and 2023. BU monitoring was performed at the exhaust fans, while TD quantifications used PRISMA hyperspectral satellite data and two different wind reanalysis datasets (i.e., GEOS-FP and ERA5). This mine also offered a unique opportunity for BU-TD comparisons before and after ventilation air methane (VAM) oxidizer operation, which was installed at one of the shafts in 2022. The results showed that TD-estimated mean fluxes were generally lower than BU data, which were attributed to quantification difficulties associated with the low albedo caused by heavy vegetation and the terrain of the area. However, despite the discrepancies in mean emissions, both the interquartile range and the data range of the distributions generally overlapped, and the estimates correctly showed the emission trends.

Virginia↗
USGS2026

Effects of human activity on Khumbu Glacier: Towards a sustainable Everest Base Camp

The Everest Base Camp (EBC), a critical staging area for mountaineers attempting to summit Mount Everest, has seen a significant increase in human activity over the past few decades. The increasing number of climbers, trekkers, and support staff, particularly at the campsite area, has intensified environmental impacts on the Khumbu Glacier. This study investigates the effects of human activity at EBC, focusing on three contributors to glacier melt: global climate change, local fossil fuel consumption, and human urinary discharge. Through field data collected during the spring 2023 climbing season, this study examines the energy released through the use of liquefied petroleum gas (LPG), kerosene, and petrol for cooking, heating, and electricity, as well as the heat generated by human urinary discharge. The total heat energy released from these activities at EBC in the spring of 2023 was 849,174 ± 179,774 MJ, which would be sufficient to melt 2492 ± 528 tons of glacier ice and snow. Trends in land surface temperature at EBC and adjacent environments during the 1991–2023 period were calculated from Landsat 5, 7, and 9 satellite data. The 32-year Landsat record reveals that EBC surface temperature increased by 0.28 °C year −1 , which was roughly twice the rate of warming at the surface of Khumbu Glacier adjacent to EBC and about 15% higher than the temperature increase of the debris-covered terrain immediately to the north of EBC. The findings suggest that addressing these anthropogenic influences could help to preserve the Khumbu Glacier and support the sustainability of mountaineering in the region. This paper also takes a transformative approach and explores stable and safer locations in case there is a necessity to relocate the current EBC. Two sites have been identified southwest of the existing EBC, which, unlike the current supraglacial site, are situated on stable ground.

Everest Base Camp, Khumbu Glacier↗
USGS2025

The influence of postfire root strength decay on shallow landslide susceptibility in western Oregon

Severe wildfires can increase the potential for postfire landslides, partly due to the loss of vegetation and root reinforcement. This study investigates the multi-year dynamics of root strength following wildfire in coniferous forests of western Oregon and evaluates implications for postfire slope stability. Root tensile strength was measured through laboratory testing of 969 roots collected from 26 test pits over four postfire timeframes: unburned conditions, 4–12 months postfire (Cedar Creek, 2022), 24–35 months postfire (Holiday Farm, 2020), and 59–67 months postfire (Eagle Creek, 2017). Variables analyzed include root diameter, tensile thread strength, stiffness, progressive tensile strength with displacement, time-dependent strength and ductility, root area ratio, and cohesion. Results indicate a 50% reduction in root strength after wildfire, with the lowest values occurring around 4 years postfire. Projections indicate it could take 10–22 years following fire for root strength to recover to 70% and 90% of the prefire original root strength, respectively. This timeline indicates a multi-year window of vulnerability for shallow landslides that reaches a maximum approximately 4 years following fire. Brittle failure behavior was observed in burned root systems compared to unburned root systems, signaling a more abrupt loss of strength at yield. Simple slope stability analyses show that reduced root strength can lead to instability on more gentle slopes and under less saturation in comparison to unburned conditions. These findings highlight the critical role of root reinforcement in postfire slope stability and the long-term implications of wildfire disturbance on landslide susceptibility.

Oregon↗
USGS2024

Effects of drought and cloud-water interception on groundwater recharge and wildfire hazard for recent and future climate conditions, Kauaʻi, Oʻahu, Molokaʻi, Maui, and the Island of Hawaiʻi

The Water-budget Accounting for Tropical Regions Model (WATRMod) code was used for Kauaʻi, Oʻahu, Molokaʻi, Maui, and the Island of Hawaiʻi to estimate the spatial distribution of groundwater recharge, soil moisture, evapotranspiration, and climatic water deficit for a set of water-budget scenarios. The scenarios included historical and future drought conditions, and a land-cover condition where shrubland and forest within the cloud zone were converted to grassland. For the historical drought condition, island-wide mean annual recharge estimates range from a decrease of 30 percent (239 million gallons per day [Mgal/d]) for Kauaʻi to a decrease of 39 percent (2,706 Mgal/d) for the Island of Hawaiʻi, relative to the reference condition consisting of 1978–2007 rainfall and 2020 land cover. For the future drought condition, estimates of island-wide mean annual recharge range from a decrease of 40 percent (477 Mgal/d) on Maui to a decrease of 51 percent (116 Mgal/day) on Molokaʻi. Complete conversion of all shrubland and forest within the cloud zone to grassland for each drought condition produces estimated land-cover-related decreases in island-wide mean annual recharge (in addition to the drought-related decreases) of 11–12 Mgal/d on Oʻahu, 119–135 Mgal/day on Maui, and 689–849 Mgal/d on the Island of Hawaiʻi. The spatial distributions of increases in conditions indicative of moisture stress and potential wildfire hazard were quantified using the relative frequency of soil moisture less than a selected threshold value (monthly mean soil moisture less than 0.074, expressed as a fraction of available water capacity), evapotranspiration less than a selected threshold value (monthly evapotranspiration less than 0.96 inches), and climatic water deficit greater than a selected threshold value (monthly climatic water deficit greater than 0.77, expressed as fraction of potential evapotranspiration). For the historical drought condition, the greatest increases in the relative frequency for the moisture-stress indicators occur across parts of east and southwest Kauaʻi; central, east, and west Oʻahu; central Molokaʻi; central Maui and low- to mid-altitude parts of West Maui volcano; and the northwestern and southern parts of the Island of Hawaiʻi. For the future drought condition, the greatest increases in the relative frequency of the moisture-stress indicators occur across parts of west Kauaʻi; central and west Oʻahu and Molokaʻi; a band of mid-altitude area on the southern slope of West Maui volcano and across the southwestern slope of Haleakalā; and mid-altitude areas of the northwestern and southern parts of the Island of Hawaiʻi. Complete conversion of all shrubland and forest within the cloud zone to grassland for each drought condition results in land-cover-related increases in the relative frequency of moisture-stress indicators around Kaʻala in the Waiʻanae Range and the southeastern part of the Koʻolau Range on Oʻahu, the southern part of West Maui volcano and the southwestern slope of Haleakalā on Maui, and the upland regions of the western and southern parts of the Island of Hawaiʻi.

Hawaii↗