Search USGS⌕ Search

SEARCH · Search USGS

Results for “Natural Hazards”

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.

486 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↗
USGS2024

Distinguishing natural sources from anthropogenic events in seismic data

As seismic data are increasingly used to investigate a diverse range of subsurface phenomena beyond regular fast-rupturing earthquakes (Peng and Gomberg, 2010; Beroza and Ide, 2011), it is important to acknowledge that human-generated ground vibrations may be mistaken for naturally generated subsurface processes (Larose et al., 2015; Li et al., 2018). Correct discrimination of natural processes from anthropogenic noise is especially pressing given the trend in seismic detection research toward automated algorithms and machine learning methods (Yoon et al., 2015; Kong et al., 2019;Mousavi and Beroza, 2022) and the growth in seismic data collection in new environments such as urban and industry settings (e.g., Díaz et al.,2017).

Seismological Research Letters↗
USGS2017

Field-trip guide to Mount Hood, Oregon, highlighting eruptive history and hazards

This guidebook describes stops of interest for a geological field trip around Mount Hood volcano. It was developed for the 2017 International Association of Volcanology and Chemistry of the Earth’s Interior (IAVCEI) Scientific Assembly in Portland, Oregon. The intent of this guidebook and accompanying contributions is to provide an overview of Mount Hood, including its chief geologic processes, magmatic system, eruptive history, local tectonics, and hazards, by visiting a variety of readily accessible localities. We also describe coeval, largely monogenetic, volcanoes in the region. Accompanying the field-trip guidebook are separately authored contributions that discuss in detail the Mount Hood magmatic system and its products and behavior (Kent and Koleszar, this volume); Mount Hood earthquakes and their relation to regional tectonics and the volcanic system (Thelen and Moran, this volume); and young surface faults cutting the broader Mount Hood area whose extent has come to light after acquisition of regional light detection and ranging coverage (Madin and others, this volume). The trip makes an approximately 175-mile (280-kilometer) clockwise loop around Mount Hood, starting and ending in Portland. The route heads east on Interstate 84 through the Columbia River Gorge National Scenic Area. The guidebook points out only a few conspicuous features of note in the gorge, but many other guides to the gorge are available. The route continues south on the Mount Hood National Scenic Byway on Oregon Route 35 following Hood River, and returns to Portland on U.S. Highway 26 following Sandy River. The route traverses rocks as old as the early Miocene Eagle Creek Formation and overlying Columbia River Basalt Group of middle Miocene age, but chiefly lava flows and clastic products of arc volcanism of late Miocene to Holocene age.

Oregon↗
USGS2024

Geologic input databases for the 2025 Puerto Rico – U.S. Virgin Islands National Seismic Hazard Model update: Crustal faults component

The last National Seismic Hazard Model (NSHM) for Puerto Rico and the U.S. Virgin Islands (PRVI) was published in 2003. In advance of the 2025 PRVI NSHM update, we created three geologic input databases to summarize new onshore and offshore fault source information in the northern Caribbean region between 62°–70° W and 16°–21° N. These databases, of fault sections, fault‐zone polygons, and geologic estimates of fault activity (fault‐slip rate and earthquake recurrence intervals) at specific sites, document updates to fault parameters used in prior seismic hazard models in PRVI. Fault sources were reviewed from published studies since 2003, which document substantial changes to the understanding of fault location, geometry, or activity. New fault section sources were added for features that meet the criteria of (1) length ≥7 km, (2) unequivocal evidence of recurrent tectonic Quaternary activity, and (3) documentation that is publicly available in a peer‐reviewed source. In addition, we revised several broad areal sources, such as the Mona and Anegada extensional zones. The 2003 model included three fault sections and two fault‐zone polygons (areal sources). These databases include 35 fault sections, 6 fault‐zone polygons, and 51 earthquake geology sites. To characterize fault activity rates, slip‐rate bins were assigned based on landscape expression and paleoseismic trench observations for faults without published slip‐rate sites. Additional fault sources were evaluated but not included in these databases due to a lack of published information about fault location, geometry, or recurrent Quaternary activity. The PRVI NSHM 2025 geologic input databases describe crustal faulting; the geometries and coupling of Puerto Rico subduction zone and Muertos Trough models are considered in a separate database. Updates to the fault sections, fault‐zone polygons, and earthquake geology databases can help inform the location and recurrence rate of damaging earthquakes in the PRVI NSHM implementation.

Puerto Rico, U.S. Virgin Islands↗
USGS2026

Beaver dam analogs as nature-based solutions to mitigate snowpack loss in northern New Mexico

Reductions in snow have left many streams in northern New Mexico dry or with very low flows during the summer. Base flow, defined as the contribution of groundwater to streamflow, can sustain streamflow during the summer and during drought conditions. Beaver dam analogs can be used to “slow the flow,” or increase the infiltration of rain-based runoff to replenish shallow groundwater reservoirs and increase base flow and potentially increase summertime streamflow.

New Mexico↗
USGS2025

Comparative assessment of a restored and natural wetland using 13C-DNA SIP reveals a higher potential for methane production in the restored wetland

Wetlands are the largest natural source of methane (CH 4 ), a potent greenhouse gas produced by methanogens. Methanogenesis rates are controlled by environmental factors such as redox potential, temperature, and carbon and electron acceptor availability and are presumably dependent on the composition of the active methanogen community. We collected intact soil cores from a restored and natural freshwater depressional wetland on Maryland’s Delmarva Peninsula (USA) to assess the effects of wetland restoration and redox shifts on microbial processes. Intact soil cores were incubated under either saturated (anoxic) or unsaturated (oxic) conditions and amended with 13 C-acetate for quantitative stable isotope probing (qSIP) of the 16S rRNA gene. Restored wetland cores supported a distinct community of methanogens compared to natural cores, and acetoclastic methanogens putatively identified in the genus Methanosarcina were among the most abundant taxa in restored anoxic and oxic cores. The active microbial communities in the restored wetland cores were also distinguished by the unique presence of facultatively anaerobic bacteria belonging to the orders Firmicutes and Bacteroidetes . In natural wetland incubations, methanogen populations were not among the most abundant taxa, and these communities were instead distinguished by the unique presence of aerobic bacteria in the phyla Acidobacteria , Actinobacteria , and class Alphaproteobacteria . Iron-reducing bacteria, in the genus Geobacter , were active across all redox conditions in both the restored and the natural cores, except the natural oxic–anoxic condition. These findings suggest an overall higher potential for methanogenesis in the restored wetland site compared to the natural wetland site, even when there is evidence of Fe reduction.

Maryland↗
USGS2025

A crosswalk of the 2015 World Terrestrial Ecosystems to the International Union for the Conservation of Nature Global Ecosystem Typology Framework

To support ecosystem mapping and accounting applications, we aligned the 2015 U.S. Geological Survey/Esri/The Nature Conservancy-World Terrestrial Ecosystems (WTEs) with the International Union for Conservation of Nature Global Ecosystem Typology (GET) framework. This process, known as “crosswalking,” enabled the development of a global map of GET level 3 Ecosystem Functional Groups (EFGs) at a 250-meter spatial resolution. Crosswalking involved manually assigning 1,781 biogeographically stratified WTEs to their most probable EFG based on similarities in climate, terrain, vegetation, and geographic distribution. We compared attributes of the WTE dataset with summary characteristics of the EFGs. The resulting crosswalked global map of International Union for Conservation of Nature GET ecosystems is intended to be useful for standardizing ecosystem classification and reporting under frameworks such as the Kunming-Montreal Global Biodiversity Framework and the United Nations System of Environmental-Economic Accounting. We discuss key challenges in reconciling non-identical classifications, such as many-to-one relationships and variation in data quality.

Open-File Report↗
USGS2026

The Great Lakes Geologic Mapping Coalition—Working collaboratively to understand the geology of the Great Lakes Region

Introduction The Great Lakes Geologic Mapping Coalition (GLGMC), commonly referred to as the “Coalition,” is a partnership between the U.S. Geological Survey (USGS), the U.S. States of Illinois, Indiana, Michigan, Minnesota, New York, Ohio, Pennsylvania, and Wisconsin and the Canadian province of Ontario. The member States receive funding for geologic mapping work from the USGS National Cooperative Geologic Mapping Program (NCGMP), whereas Ontario participates as a nonfunded partner. The mission of the GLGMC is to produce three-dimensional (3D) geologic maps that depict unconsolidated sediments and near-surface bedrock in the Great Lakes region of North America. Geologic maps are the basis of most earth science investigations and help support resource exploration (energy, minerals, groundwater), natural hazard mitigation, infrastructure development, and land-use planning, all of which can be used to advance economic development and strengthen national security in the Great Lakes region. During the last few million years, the Great Lakes region has experienced repeated glacial advances and retreats, leaving behind extensive sediments, abundant natural resources, and widespread effects on the underlying bedrock geology (Swezey and others, 2022). Linked by shared histories of past glaciations, industrial agriculture, and legacy automotive, coal, steel, and manufacturing industries, the GLGMC member States collaborate to improve the understanding of the 3D distribution of the sediments overlying the region’s bedrock (fig. 1). Developing a comprehensive subsurface 3D framework of this glaciated terrain can provide earth science data to policymakers at all levels. These insights facilitate informed decisions on the exploration, use, and protection of vital resources, such as critical minerals, industrial materials, and aquifers, thereby supporting economic prosperity and the well-being of the citizens of this region. Since its inception in 1998, the Coalition has completed more than 100 geologic mapping projects across the Great Lakes region. Each project aims to deliver geologic maps, 3D datasets, and other information that improves understanding of the geology of the Great Lakes region, with an emphasis on economic and water resources. Key deliverables include 3D geologic maps and models typically portraying sediment thickness, often derived from top-of-bedrock and borehole data. These products are developed through a combination of fieldwork, subsurface modeling, and the collection and analysis of rock and sediment cores. To support Coalition goals, member States collaborate with scientists working on related STATEMAP, EDMAP, and FEDMAP projects. Coalition scientists also engage with Tribal Nations in the Great Lakes region to ensure that Tribal interests pertaining to Coalition work are addressed. Through this collaboration, the Coalition unites the efforts of State, Federal, and Tribal Nation stakeholders to advance geologic data production and enhance understanding of the geologic resources of the Great Lakes region.

Illinois, Indiana, Michigan, Minnesota, New York, ↗
USGS2025

Assessment of western Oregon debris-flow hazards in burned and unburned environments

In the steep and mountainous environment of western Oregon, debris flows pose a considerable threat to property, infrastructure and life. Wildfire is commonly known to increase the susceptibility of steep slopes to debris flows, but the extent of this process in the western Cascades is not well understood. The US Geological Survey (USGS) currently estimates postfire debris-flow likelihood and triggering rainfall thresholds using a model calibrated to a southern California inventory of debris flows generated by excess runoff within the first year after fire. Because of a lack of available data, this model has not been tested in western Oregon, or in locations where postfire debris flows initiate via other mechanisms (e.g., shallow landslides or in-channel failures). Using repeat field observations and aerial imagery, we developed two new debris-flow inventories within and adjacent to the perimeters of five 2020 wildfires in western Oregon: Archie Creek, Holiday Farm, Beachie Creek, Lionshead and Riverside. The first inventory focuses on postfire debris flows (2020–2022); the second focuses on debris flows prior to fires (1995–2020). Our inventories of prefire and postfire debris flows were used to document initiation mechanisms in Oregon's western Cascades and to evaluate the effects of wildfire. We found that wildfire changed the distribution of debris-flow initiation mechanisms in the western Cascades. After the wildfires, annual rates of runoff-generated debris flows increased by 22% and the number of shallow landslide-initiated debris flows decreased by 17% relative to before the wildfires. Despite this shift, shallow landsliding was the dominant debris-flow initiation mechanism in both unburned and burned environments. We found the performance of the current USGS debris-flow likelihood model was degraded relative to other previously tested locations across the intermountain western United States. Our results highlight the need for improved postfire hazard assessment in western Oregon based on regional model calibration that is tuned to the dominant debris-flow initiation mechanisms.

Oregon↗
USGS2026

Development of liquefaction-targeted design basis in U.S. seismic provisions

Seismic design criteria for new buildings in the United States have historically focused on life safety performance objectives through minimizing the potential for structural collapse. Development of design criteria to meet this performance objective has evolved over time, leading to the current, risk-targeted maximum considered earthquake (MCE R ) design basis. Corresponding MCE R seismic design loads incorporate the full range of probabilistic ground motion hazard information available from U.S. Geological Survey National Seismic Hazard Models, as well as a representation of the uncertainty in structural collapse capacity. Despite these considerable advances in structural design criteria, guidelines for geotechnical hazards such as liquefaction-induced damage, which continue to be based on uniform ground-shaking hazard (i.e., MCE ground motions), are not well-connected with risk-based, collapse prevention performance objectives. This study, undertaken as part of Building Seismic Safety Council efforts to update recommended seismic provisions through the National Earthquake Hazards Reduction Program, presents probabilistic, liquefaction-targeted design objectives and corresponding ground motion parameters, with a focus on improving consistency in first-level screening criteria for liquefaction hazard assessment. This study identifies a potential maximum acceptable annualized probability of liquefaction triggering ( P L,max ) of 2.5% in 50 years, above which liquefaction consequences and potential mitigation measures would be considered. The P L,max -based design basis is calibrated to effective probabilistic liquefaction hazard levels obtained in practice using current MCE peak ground accelerations, and could be implemented via a new liquefaction-targeted maximum considered earthquake (MCE L ) peak ground acceleration ( PGA L ). The proposed PGA L is derived from probabilistic liquefaction hazard curves and includes uncertainties inherent to liquefaction hazard modeling. The potential design effects of the proposed PGA L are demonstrated via comparison with current MCE PGA levels, effects on first-level liquefaction screening analyses, and improvements in the consistency of P L,max levels across the United States, and for a wide range of site conditions.

conterminous United States↗
USGS1976

Uranium analysis of single drops of natural waters using the fission-track technique

The fission-track technique has been applied to uranium analysis of natural waters using a single drop of water. Samples may be collected in the field by evaporating a drop of water on a detector, and samples may then be analyzed by returning the detector to the laboratory. A sample collecting kit has been assembled, which facilitates the designed sampling method. Some unique advantages of this technique are picogram sensitivity (10 -12 g) for uranium on single drops of water, maintenance of sample integrity during transport and storage, and the ability to discern whether the uranium is in solution or in the suspended material. Counting the tracks is the most time-consuming part of the analysis, but it can be reduced by using automatic counting systems. The analytical error is generally ± 10 percent over the entire range of concentrations.

Open-File Report↗
USGS2025

Metaproteomics and metagenomics reveal microbial pathways of organic matter degradation and methanogenesis in a marginally producing natural gas well

The expansion of natural gas production and utilization worldwide has led to the decline of many once-productive wells, eventually resulting in costly well-plugging and unused infrastructure. However, in areas like the Michigan basin, MI, where the majority of natural gas is biogenically produced, microbial communities could potentially be stimulated to generate additional methane, increasing gas supply and reducing the need to drill new wells. In this study, we performed metaproteomic, metagenomic, and geochemical analyses of Antrim Shale formation water from a marginally producing natural gas well to evaluate resident microbial community functions in the context of potential bioenergy production. Functional proteins involved in methanogenesis, degradation/catabolism (including organic matter degradation), biosynthesis, energy utilization, transmembrane transport, and stress response were among the most commonly identified groups. Three metagenome-assembled genomes (MAGs) were characterized, including Methanomicrobiaceae, Methanothrix , and Smithella . For each, the identified proteins involved in methanogenesis and the degradation of diverse organic compounds, strongly suggest their role in utilizing shale-derived organic matter. These findings provide an increased understanding of the microorganisms and their metabolisms generating natural gas in the Antrim Shale and establish a foundation for future stimulation efforts aimed at enhancing biogenic methane production in marginal gas wells.

Michigan↗
USGS2026

Structural evolution and slip rate variations through time of the Puente-Hills blind-thrust fault beneath Los Angeles: Implications for seismic hazard and folding kinematics

Using seismic reflection profiles, historical well logging data, and luminescence and radiocarbon ages, we determine a Pleistocene-Holocene slip history for the central, Santa Fe Springs segment of the Puente Hills blind-thrust fault (PHT), a major seismogenic fault situated beneath the urbanized Los Angeles metropolitan region. We analyze the geometry of correlative stratigraphic units in the forelimb and backlimb of the overlying growth-fold, the Santa Fe Springs anticline, and determine the uplift of seven age-correlative markers. Uplift measurements are converted to thrust displacements on the underlying PHT using a structural method laid out by Don et al. (2022, https://doi.org/10.1785/0120220048 ) that accounts for the geometry of the fault. These data indicate that deep thrust displacement on the PHT is partially consumed updip in the creation of a hanging-wall fault-bend fold, with forelimb growth strata recording <80% of the slip documented within the backlimb. Chronological data from growth strata yield age constraints for folding and faulting on the underlying PHT, providing a detailed incremental slip history derived from both the forelimb and backlimb folding for seven discrete growth horizons, spanning the past 1.4 million years. The resulting six incremental slip rates demonstrate that fault slip has varied through time from the middle Pleistocene to the Holocene. Moreover, these results reveal synchronous acceleration of both the central, Santa Fe Springs and western, Los Angeles segments of the fault system since late Pleistocene time (after 200 ka) and slip rates of greater than 2 mm/yr on the downdip, backlimb, fault-ramp below the anticline.

California↗
USGS2024

Methods to evaluate and improve the modeling of rupture directivity in assessment of seismic hazard

In recent years, there have been several advancements related to the modelling of near-source effects of earthquake rupture on strong ground shaking, leading to an improved characterization of ground motions and resulting seismic hazard. Some of these modifications have stemmed from physics-based numerical modelling of the earthquake rupture process, using physics-based dynamic rupture simulations. These contributions have led to a better understanding of how fault rupture characteristics, geometry, and the style of faulting can interact with the hypocenter-dependence on the path from source to site that may ultimately guide the development of seismic directivity models. Moving forward, the application of modern techniques can be used to incorporate these source characteristics and near-fault ground motion behavior that contribute to the azimuthally varying effects that result in rupture directivity. One example is the application of machine learning methods to support more automated integration of new predictor variables in model development and open more evaluation opportunities to access residuals. Here, we utilize several techniques to take advantage of the plethora of synthetic data and its ability to supplement preexisting trends observed in data. We showcase two examples of how models can be either developed, expanded upon, or constrained using artificial neural network model (ANNs). We evaluate the performance of the ANN with existing methods, comparing misfit, potential limitations, and ability to continue to improve upon these methods in the future. One approach uses a set of simulations with corresponding synthetic ground motions from the Southern California Earthquake Center (SCEC) CyberShake study to develop a ground motion model adapted to incorporate seismic directivity information using an ANN. This large database (TBs) enables us to train the model to capture magnitude, period, and distance variations and how these parameters relate to amplification from hypocenters located along finite-faults. In some cases, there is reduced misfit from better representing source features that aren’t included in base ground motion models that neglect hypocenter location (e.g. azimuthal variation, source-to-site terms). Another ANN method uses a shallow-layered neural network model to better fit a hypocenter-independent model. This method adjusts the median and aleatory variability to account for the averaged impact of various hypocenter distributions to fit the underlying directivity adjustment model. This method serves as a template to apply to other directivity models, improving computational efficiency and more readily enabling integration in hazard codes.

California↗
USGS2025

Nodal seismic deployment on Mauna Loa volcano, Hawaii: Dataset and preliminary insights

Mauna Loa is the largest active volcano on Earth, comprising ∼51% of the Island of Hawai‘i’s landmass and posing significant risks to the island’s communities, infrastructure, and natural environment. Historical eruptions have produced lava flows that have reached the ocean in as little as 3 hr. The timing and location of such lava flows in the past 200 yr underscore how critical determination of the location and geometry of magma storage and structure is for volcanic hazard assessment and eruption forecasting. Now, after nearly 38 yr of volcanic quiescence, Mauna Loa has erupted again. On 27 November 2022, fissures initiated within the summit caldera and then migrated to the northeast rift zone, where they generated a large lava flow that threatened a major highway. To improve our understanding of the geometry of this magma system, we deployed a temporary nodal array on Mauna Loa in the summer of 2024. This increased our seismic coverage sufficiently to image this magmatic system. This nodal array consists of 33 seismometers distributed on and around the volcano and was deployed for over three months to record seismic signals. The primary objective of this project is to resolve the high‐resolution seismic velocity structure and characterize seismic features associated with magma storage and ascent pathways. In this article, we present an overview of the deployment, evaluate the quality of the data, and show example recordings to evaluate the suitability of the data set for future seismic investigations, including earthquake relocation, seismic tomography, and receiver function analysis. Comparisons with nearby permanent broadband and short‐period seismic stations demonstrate that the nodal array recorded high‐quality waveforms, making it a valuable resource for constraining the magmatic system beneath Mauna Loa at multiple scales.

Hawaii↗
USGS1996

Techniques to estimate generalized skew coefficients of annual peak streamflow for natural basins in Texas

This report presents two techniques to estimate generalized skew coefficients used for log-Pearson Type III peak-streamflow frequency analysis of natural basins in Texas. A natural basin has less than 10 percent impervious cover, and less than 10 percent of its drainage area is controlled by reservoirs. The estimation of generalized skew coefficients is based on annual peak and historical peak streamflow for all U.S. Geological Survey streamflow-gaging stations having at least 20 years of annual peak-streamflow record from natural basins in Texas. Station skew coefficients calculated for each of 255 Texas stations were used to estimate generalized skew coefficients for Texas. One technique to estimate generalized skew coefficients involved the use of regression equations developed for each of eight regions in Texas, and the other involved development of a statewide map of generalized skew coefficients. The weighted mean of the weighted mean standard errors of the regression equations for the eight regions is 0.36 log 10 skew units, and the weighted mean standard error of the map is 0.35 log 10 skew units. The technique based on the map is preferred for estimating generalized skew coefficients because of its smooth transition from one region of the State to another.

Texas↗
USGS2025

Natural source zone depletion of crude oil in the subsurface: Processes controlling mass losses of individual compounds

At many petroleum hydrocarbon spill sites, residual spilled product forms a long-term source of groundwater contamination. The phrase source zone natural depletion is used to refer to the mass loss rates. Overall mass lost under environmental conditions was analyzed using conservative biomarker concentrations for a 1979 oil spill in northern Minnesota, USA. After 40–41 years, an average of 50% of the mass was lost with values ranging from 22% to 57% depending on location. It is also important to understand the composition changes in the source. To understand controls on the losses of individual compounds, concentrations of volatile hydrocarbons in oil samples were compared with aqueous solubilities, and pore-space oil saturations. The results of the comparison show that losses of the oil compounds were controlled by pore-space oil saturations, solubility, and susceptibility to degradation under methanogenic conditions. Compounds that degrade under methanogenic conditions, including toluene, o -xylene, and n -alkanes are more depleted compared to benzene, ethylbenzene, and m - and p -xylene for which losses are dominated by dissolution. These rates and compound-specific behaviors form a foundation for improved modeling approaches and risk analyses.

Minnesota↗
USGS2026

Assessing natural recharge in Indian Wells Valley, California: A Basin Characterization Model case study

The communities in Indian Wells Valley (IWV), in the northern Mojave Desert in California, rely on groundwater for domestic and agricultural use. Mountain front recharge from the surrounding Sierra Nevada is the main source of natural recharge to the valley. Increased urbanization, agricultural development, and groundwater pumping during recent decades put IWV in a state of critical overdraft. The U.S. Geological Survey Basin Characterization Model, version 8 (BCMv8) was used to evaluate historical and future climate and hydrologic conditions in IWV. The BCMv8 estimated natural recharge in IWV at 10.7 million cubic meters (Mm 3 ) per year for the period from 1981 to 2010. Future patterns of water balance variables using three future climate scenarios, hot-wet, hot-dry, and warm-moderately wet, were calculated for mid-century (2040–69) and end-of-century (2070–99) periods. Results for both wet models projected an increase in recharge in both periods, whereas the hot-dry model projected a decrease in recharge in both periods. All models reported a large increase in seasonal variability in recharge, indicating more future availability and frequent occurrences of drought years. All climate scenarios projected an increase in climatic water deficit in both periods. These increases in irrigation demand and variability of water supply highlight the importance of strategic management planning for the sustainability of water resources in IWV.

California↗