USGS Science⌕ Search

SEARCH · USGS Science

Results for “Hydrology Journal”

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.

At least 865 records · Page 48Linked to original sources
USGS2013

Hysteresis of unsaturated hydromechanical properties of a silty soil

Laboratory tests to examine hysteresis in the hydrologic and mechanical properties of partially saturated soils were conducted on six intact specimens collected from a landslide-prone area of Alameda County, California. The results reveal that the pore-size distribution parameter remains statistically unchanged between the wetting and drying paths; however, the wetting or drying state has a pronounced influence on the water-entry pressure, the water-filled porosity at zero suction, and the saturated hydraulic conductivity. The suction stress values obtained from the shear-strength tests under both natural moisture and resaturated conditions were mostly bounded by the suction stress characteristic curves (SSCCs) obtained from the hydrologic tests. This finding experimentally confirms that the soil-water retention curve, hydraulic conductivity function, and SSCC are intrinsically related.

Journal of Geotechnical and Geoenvironmental Engin↗
USGS1999

Tributary stream infiltration as a source of herbicides in an alluvial aquifer

Where Walnut Creek flows across the South Skunk River alluvial aquifer, it provides a potential source of herbicides and herbicide metabolites. This straightened reach of the creek loses water and dissolved contaminants to the alluvial aquifer through a layer of fine-grained flood plain deposits. Estimates of potential flux of chemicals were based on measurements taken during baseflow in April 1994 before herbicides were applied to the watershed and in June 1994 after chemical application and when stream discharge included runoff and tile-drainage water. Hydraulic head measurements between the creek and flood plain deposits and between the creek and aquifer confirmed the potential for downward groundwater flow during both sampling periods. Hydraulic conductivity estimates from slug tests were used to calculate an average linear groundwater velocity of 0.5 m d −1 in the fine-grained flood plain deposits. At this velocity, contaminants could be advectively transported to the aquifer within 6 d. The potential for atrazine (2-chloro-4-ethylamino-6-isopropylamino- s -triazine) flux to the aquifer from the creek was estimated to be between 60 and 3000 µg d −1 m −2 . This rate is one to three orders of magnitude greater than the estimated flux via leaching beneath a typical field. If the process of vertical stream leakage occurs in many hydrologic settings, it may constitute a substantial source of herbicides to shallow alluvial aquifers in many areas of the Midwest.

Iowa↗
USGS2014

Nitrate fate and transport through current and former depressional wetlands in an agricultural landscape, Choptank Watershed, Maryland, United States

Understanding local groundwater hydrology and geochemistry is critical for evaluating the effectiveness of wetlands at mitigating agricultural impacts on surface waters. The effectiveness of depressional wetlands at mitigating nitrate (NO 3 ) transport from fertilized row crops, through groundwater, to local streams was examined in the watershed of the upper Choptank River, a tributary of Chesapeake Bay on the Atlantic Coastal Plain. Hydrologic, geochemical, and water quality data were collected from January of 2008 through December of 2009 from surface waters and networks of piezometers installed in and around current or former depressional wetlands of three major types along a gradient of anthropogenic alteration: (1) natural wetlands with native vegetation (i.e., forested); (2) prior-converted croplands, which are former wetlands located in cultivated fields; and (3) hydrologically restored wetlands, including one wetland restoration and one shallow water management area. These data were collected to estimate the orientation of groundwater flow paths and likely interactions of groundwater containing NO 3 from agricultural sources with reducing conditions associated with wetlands of different types. Natural wetlands were found to have longer periods of soil saturation and reducing conditions conducive to denitrification compared to the other wetland types studied. Because natural wetlands are typically located in groundwater recharge areas along watershed divides, nitrogen (N) from nearby agriculture was not intercepted. However, these wetlands likely improve water quality in adjacent streams via dilution. Soil and geochemical conditions conducive to denitrification were also present in restored wetlands and prior-converted croplands, and substantial losses of agricultural NO 3 were observed in groundwater flowing through these wetland sediments. However, delivery of NO 3 from agricultural areas through groundwater to these wetlands resulting in opportunities for denitrification were limited, particularly where reducing conditions did not extend throughout the entire thickness of the surficial aquifer allowing NO 3 to pass conservatively beneath a wetland along deeper groundwater flow paths. The complexity of N fate and transport associated with depressional wetlands complicates the understanding of their importance to water quality in adjacent streams. Although depressional wetlands often contribute low NO 3 water to local streams, their effectiveness as landscape sinks, for N from adjacent agriculture varies with natural conditions, such as the thickness of the aquifer and the extent of reducing conditions. Measurement of such natural geologic, hydrologic, and geochemical conditions are therefore fundamental to understanding N mitigation in individual wetlands.

Maryl↗
USGS2026

CrustCheck: A conceptual framework for interpreting biocrust trait–function relationships in drylands

Introduction Biocrusts are increasingly recognized for their contribution to soil stability, hydrology, nutrient cycling, and energy balance in drylands, yet practical field methods for assessing functional roles are limited. Existing approaches typically rely on cover estimates, taxonomic identification, or laboratory analyses, providing little insight into functional variation among biocrust types. To address this gap, we developed CrustCheck, a literature-derived conceptual framework that synthesizes potential relationships between easily observable biocrust traits and key ecological functions. The framework provides a foundation for future empirical validation of these relationships. Objectives We aimed to: (1) identify visual characteristics of biocrusts potentially associated with key ecological functions, and (2) synthesize these relationships into a literature-derived conceptual framework for comparative interpretation of potential biocrust contributions to landscape function. Methods We conducted a structured review of global literature on biocrusts' contribution to landscape function, drawing on multiple databases, targeted journals, and identifying relevant scientific publications from reference lists. We screened sources for information on relationships between visual traits and metrics of landscape function, synthesized relationships thematically, and used them to develop a multi-trait, function-oriented conceptual framework to inform further testing of these hypotheses in the field. Results Our review revealed that visual traits such as thickness/height, patchiness, cover, and surface roughness emerged as potential indicators of surface stability, hydrological function, and nutrient accumulation. Community composition was also frequently associated with these functions, while color reflected energy balance and carbon fixation potential. Conclusions Literature reveals that visual traits of biocrusts may serve as proxies for landscape function. CrustCheck is presented here as a conceptual synthesis framework intended to support interpretation of potential trait–function relationships rather than as a quantitatively validated predictive tool.

Restoration Ecology↗
USGS2018

Interseismic ground deformation and fault slip rates in the greater San Francisco Bay Area from two decades of space geodetic data

The detailed spatial variations of strain accumulation and creep on major faults in the northern San Francisco Bay Area (North Bay), which are important for seismic potential and evaluation of natural hazards, remain poorly understood. Here we combine interferometric synthetic aperture radar data from the ERS‐1/2 and Envisat satellites between 1992 and 2010 with continuous and campaign GPS data to obtain a high spatial and temporal coverage of ground deformation of the North Bay. The SAR data from both ascending and descending orbits are combined to separate horizontal and vertical components of the deformation. We jointly invert the horizontal component of the mean velocities derived from these data to infer the deep strike‐slip rates on major locked faults. We use the estimated deep rates to simulate the long‐wavelength deformation due to interseismic elastic strain accumulation along these locked faults. After removing the long‐wavelength signal from the InSAR horizontal mean velocity field, we estimate fault‐parallel surface creep rates of up to 2 mm/year along the central section of the Rodgers Creek fault and surface creep rates ranging between 2 and 4 mm/year along the Concord fault. No surface creep is geodetically resolved along the West Napa and Green Valley fault zones. We identified characteristically repeating earthquakes on the Rodgers Creek fault, the West Napa fault, the Green Valley fault, and the Concord fault. Nontectonic deformation in the Geysers geothermal field and in Late Cenozoic basins (Rohnert Park and Sonoma basins) are also observed, likely due to hydrological and sediment‐compaction processes, respectively.

Journal of Geophysical Research B: Solid Earth↗
USGS2007

Altered stream-flow regimes and invasive plant species: The Tamarix case

Aim: To test the hypothesis that anthropogenic alteration of stream-flow regimes is a key driver of compositional shifts from native to introduced riparian plant species. Location: The arid south-western United States; 24 river reaches in the Gila and Lower Colorado drainage basins of Arizona. Methods: We compared the abundance of three dominant woody riparian taxa (native Populus fremontii and Salix gooddingii, and introduced Tamarix) between river reaches that varied in stream-flow permanence (perennial vs. intermittent), presence or absence of an upstream flow-regulating dam, and presence or absence of municipal effluent as a stream water source. Results: Populus and Salix were the dominant pioneer trees along the reaches with perennial flow and a natural flood regime. In contrast, Tamarix had high abundance (patch area and basal area) along reaches with intermittent stream flows (caused by natural and cultural factors), as well as those with dam-regulated flows. Main conclusions: Stream-flow regimes are strong determinants of riparian vegetation structure, and hydrological alterations can drive dominance shifts to introduced species that have an adaptive suite of traits. Deep alluvial groundwater on intermittent rivers favours the deep-rooted, stress-adapted Tamarix over the shallower-rooted and more competitive Populus and Salix. On flow-regulated rivers, shifts in flood timing favour the reproductively opportunistic Tamarix over Populus and Salix, both of which have narrow germination windows. The prevailing hydrological conditions thus favour a new dominant pioneer species in the riparian corridors of the American Southwest. These results reaffirm the importance of reinstating stream-flow regimes (inclusive of groundwater flows) for re-establishing the native pioneer trees as the dominant forest type. ?? 2007 The Authors Journal compilation ?? 2007 Blackwell Publishing Ltd.

Global Ecology and Biogeography↗
USGS2025

Wet antecedent soil moisture increases atmospheric river streamflow magnitudes non-linearly

Atmospheric rivers (ARs) drive most riverine floods on the United States (U.S.) West Coast. However, estimating flood risk based solely on AR intensity and duration is challenging because precipitation phase, antecedent conditions, and physical watershed characteristics (e.g., slope and soil depth) can influence the magnitude of floods. Here, we analyze how antecedent soil moisture (ASM) conditions contribute to variability in streamflow during AR events and how that changes across climatic regimes and physiography in 122 U.S. West Coast watersheds. We identify a robust non-linear relationship between streamflow and ASM during ARs in 89% of watersheds. The inflection point in this relationship represents a watershed-specific critical ASM threshold, above which event maximum streamflow is, on average, two to four and a half times larger. Wet ASM conditions amplify the hydrologic impacts of more frequent but weaker, lower moisture transport AR events, while dry ASM conditions attenuate the hydrologic impacts that stronger, higher moisture transport AR events could otherwise cause. Our research shows that watersheds prone to ASM-amplified streamflows have higher evaporation ratios, lower cold-season precipitation, lower snow-to-rain ratios, and shallower, clay-rich soils. Higher evaporation and lower precipitation lead to greater ASM variability during the cold season, increasing streamflow during wet periods and buffering streamflow during dry periods. Lower snow fraction and shallower soils limit the antecedent water storage capacity of a watershed, contributing to greater sensitivity of streamflow peaks to ASM variability. Incorporating ASM thresholds into hydrologic models in these regions prone to AR-amplified streamflow could improve forecasts and decrease uncertainty.

California, Nevada, Oregon, Washington↗
USGS2025

Alaskan hydrology in transition: Changing precipitation and evapotranspiration patterns are projected to reshape seasonal streamflow and water temperature by midcentury (2035-2064)

High spatial and temporal resolution models are essential for understanding future climate impacts and developing effective climate resilience plans. However, existing regional and global river models often lack the resolution needed to accurately capture local conditions. This study uses a series of high-resolution models, including the Regional Arctic System Model, mizuRoute, and the river basin model, to analyze Arctic and sub-Arctic Alaskan hydrology. We compare a historical baseline (1991–2020) with six midcentury (2035–64) futures: two pseudo–global warming scenarios based on historical meteorology and four direct dynamically downscaled global climate models. The six futures reveal significant uncertainty in future annual discharge and peak flows, although a widespread increase in discharge during April (+63%) and October (+31%) is consistently shown across models. Projected increases in rain and shifting weather patterns lead to a transition from snow to rain in spring and autumn, reducing the fraction of snowmelt contributing to river discharge. Rising evapotranspiration moderates discharge changes, particularly in autumn, by offsetting precipitation increases. Average summer river temperatures are projected to increase by approximately 1.5°C, doubling the number of river segments that experience 18°C days, a critical threshold for salmon survival, and intensifying the heat flux to the ocean adding an average of 3.3 × 10 12 MJ yr −1 . These changes in the hydrologic cycle could profoundly impact riverine and oceanic ecosystems, posing substantial challenges to communities reliant on these environments.

Alaska, Yukon↗
USGS2017

Projecting impacts of climate change on water availability using artificial neural network techniques

Lago Loíza reservoir in east-central Puerto Rico is one of the primary sources of public water supply for the San Juan metropolitan area. To evaluate and predict the Lago Loíza water budget, an artificial neural network (ANN) technique is trained to predict river inflows. A method is developed to combine ANN-predicted daily flows with ANN-predicted 30-day cumulative flows to improve flow estimates. The ANN application trains well for representing 2007–2012 and the drier 1994–1997 periods. Rainfall data downscaled from global circulation model (GCM) simulations are used to predict 2050–2055 conditions. Evapotranspiration is estimated with the Hargreaves equation using minimum and maximum air temperatures from the downscaled GCM data. These simulated 2050–2055 river flows are input to a water budget formulation for the Lago Loíza reservoir for comparison with 2007–2012. The ANN scenarios require far less computational effort than a numerical model application, yet produce results with sufficient accuracy to evaluate and compare hydrologic scenarios. This hydrologic tool will be useful for future evaluations of the Lago Loíza reservoir and water supply to the San Juan metropolitan area.

Puerto Rico↗
USGS2007

The concept of hydrologic landscapes

Hydrologic landscapes are multiples or variations of fundamental hydrologic landscape units. A fundamental hydrologic landscape unit is defined on the basis of land-surface form, geology, and climate. The basic land-surface form of a fundamental hydrologic landscape unit is an upland separated from a lowland by an intervening steeper slope. Fundamental hydrologic landscape units have a complete hydrologic system consisting of surface runoff, ground-water flow, and interaction with atmospheric water. By describing actual landscapes in terms of land-surface slope, hydraulic properties of soils and geologic framework, and the difference between precipitation and evapotranspiration, the hydrologic system of actual landscapes can be conceptualized in a uniform way. This conceptual framework can then be the foundation for design of studies and data networks, syntheses of information on local to national scales, and comparison of process research across small study units in a variety of settings. The Crow Wing River watershed in central Minnesota is used as an example of evaluating stream discharge in the context of hydrologic landscapes. Lake-research watersheds in Wisconsin, Minnesota, North Dakota, and Nebraska are used as an example of using the hydrologic-landscapes concept to evaluate the effect of ground water on the degree of mineralization and major-ion chemistry of lakes that lie within ground-water flow systems.

Minnesota, Nebraska, North Dakota, Wisconsin↗
USGS2026

Assessing environmental drivers of denitrification in restored riverine floodplains

Restoration of impaired floodplains is an increasingly prevalent strategy for alleviating water quality concerns and reducing downstream flooding at watershed scales. Floodplains temporarily store water and slow flow velocity to promote sedimentation during overbank flooding and remove inorganic nitrogen from floodwater and groundwater via denitrification. Evaluating the impacts of different restoration strategies on denitrification can inform more strategic investments into floodplain modifications that improve water quality outcomes. Our research investigates how denitrification rates in floodplains respond to environmental factors that are actionable from an engineering perspective through design and water resources management. We seasonally measured soil denitrification enzyme activity and various environmental characteristics in 4 floodplains with different restoration design and management approaches at the confluence of the Wabash and Tippecanoe Rivers in Indiana, United States. Our results showed that denitrification rates in an agricultural floodplain were significantly lower than in restored floodplains with native vegetation. Certain soil conditions characteristic of floodplain wetlands were associated with higher denitrification, particularly elevated total nitrogen, moisture, silt, and organic matter contents. Vegetation species composition was correlated with denitrification rates. This link may reflect the direct effects of vegetation on soil conditions, such as supplying labile organic carbon, or indirect effects, such as vegetation acting as an indicator of hydrologic regime and land use. Denitrification seasonally varied, peaking in winter when nitrate supply from rivers draining agricultural watersheds in the region is also high. Substrate limitation of soil denitrification enzyme activity was most significant during the summer when overbank flooding, which replenishes soil nitrogen stocks, rarely occurs. Our findings indicate that denitrification capacity will likely be maximized in riverine floodplains that are restored as wetlands with diverse native vegetation and enhanced hydrologic connectivity. Such restoration activities promote higher denitrification rates via elevated moisture, fine sediment deposition, and soil organic matter.

Tippecanoe River, Wabash River↗
USGS2018

Tropical wetlands in the Anthropocene: The critical role of wet-dry cycles

In the face of climate change and increasing human water demands for agriculture, industry, and cities, the fate of wetland ecosystems in tropical wet-dry climates is threatened. To maximize biodiversity and ecological resilience, the value of the ecosystem services provided by tropical wetlands can be incorporated into regional land use and water management decisions. Environmental planners and resource managers can work to protect both the “dry” and “wet” phases of the wet-dry hydrologic cycles. These cycles have shaped and maintained these ecosystems in the past and they can be used to maximize biodiversity and resilience in the future.

Solutions Journal↗
USGS1996

A test of vegetation-related indicators of wetland quality in the prairie pothole region

This study was part of an effort by the U.S. Environmental Protection Agency to quantitatively assess the environmental quality or 'health' of wetland resources on regional and national scales. During a two-year pilot study, we tested selected indicators of wetland quality in the U.S. portion of the prairie pothole region (PPR). We assumed that the amount of cropland versus non-cropland (mostly grassland) in the plots containing these basins was a proxy for their quality. We then tested indicators by their ability to discriminate between wetlands at the extremes of that proxy. Amounts of standing dead vegetation were greater in zones of greater water permanence. Depth of litter was greater in zones of greater water permanence and in zones of basins in poor-quality watersheds. Amounts of unvegetated bottom were greater in basins in poor-quality watersheds; lesser amounts occurred in all wetlands during a wetter year. Greater amounts of open water occurred during a wetter year and in zones of greater water permanence. When unadjusted for areas (ha) of communities, plant taxon richness was higher in wet-meadow and shallow-marsh zones in good-quality watersheds than in similar zones in poor-quality watersheds. Wet-meadow zones in good-quality watersheds had greater numbers of native perennials than those in poor-quality watersheds. This relation held when we eliminated all communities in good-quality watersheds larger than the largest communities in poor-quality watersheds from the data set. We conclude that although amounts of unvegetated bottom and plant taxon richness in wet-meadow zones were useful indicators of wetland quality during our study, the search for additional such indicators should continue. The value of these indicators may change with the notoriously unstable hydrological conditions in the PPR. Most valuable would be indicators that could be photographed or otherwise remotely sensed and would remain relatively stable under various hydrological conditions. An ideal set of indicators could detect the absence of stressors, as well as the presence of structures or functions, of known value to major groups of organisms.

Journal of Aquatic Ecosystem Health↗
USGS2011

Climatic controls on the snowmelt hydrology of the northern Rocky Mountains

The northern Rocky Mountains (NRMs) are a critical headwaters region with the majority of water resources originating from mountain snowpack. Observations showing declines in western U.S. snowpack have implications for water resources and biophysical processes in high-mountain environments. This study investigates oceanic and atmospheric controls underlying changes in timing, variability, and trends documented across the entire hydroclimatic-monitoring system within critical NRM watersheds. Analyses were conducted using records from 25 snow telemetry (SNOTEL) stations, 148 1 April snow course records, stream gauge records from 14 relatively unimpaired rivers, and 37 valley meteorological stations. Over the past four decades, midelevation SNOTEL records show a tendency toward decreased snowpack with peak snow water equivalent (SWE) arriving and melting out earlier. Temperature records show significant seasonal and annual decreases in the number of frost days (days ≤0°C) and changes in spring minimum temperatures that correspond with atmospheric circulation changes and surface–albedo feedbacks in March and April. Warmer spring temperatures coupled with increases in mean and variance of spring precipitation correspond strongly to earlier snowmeltout, an increased number of snow-free days, and observed changes in streamflow timing and discharge. The majority of the variability in peak and total annual snowpack and streamflow, however, is explained by season-dependent interannual-to-interdecadal changes in atmospheric circulation associated with Pacific Ocean sea surface temperatures. Over recent decades, increased spring precipitation appears to be buffering NRM total annual streamflow from what would otherwise be greater snow-related declines in hydrologic yield. Results have important implications for ecosystems, water resources, and long-lead-forecasting capabilities.

Idaho, Montana↗
USGS2009

Hydrologic and biogeochemical controls of river subsurface solutes under agriculturally enhanced ground water flow

The relative influences of hydrologic processes and biogeochemistry on the transport and retention of minor solutes were compared in the riverbed of the lower Merced River (California, USA). The subsurface of this reach receives ground water discharge and surface water infiltration due to an altered hydraulic setting resulting from agricultural irrigation. Filtered ground water samples were collected from 30 drive point locations in March, June, and October 2004. Hydrologic processes, described previously, were verified by observations of bromine concentrations; manganese was used to indicate redox conditions. The separate responses of the minor solutes strontium, barium, uranium, and phosphorus to these influences were examined. Correlation and principal component analyses indicate that hydrologic processes dominate the distribution of trace elements in the ground water. Redox conditions appear to be independent of hydrologic processes and account for most of the remaining data variability. With some variability, major processes are consistent in two sampling transects separated by 100 m.

Journal of Environmental Quality↗
USGS1987

Diagenesis and fluid flow in the San Juan Basin, New Mexico - regional zonation in the mineralogy and stable isotope composition of clay minerals in sandstone.

The Westwater Canyon Member of the Upper Jurassic Morrison Formation is a relatively homogeneous, hydrologically continuous 100-m-thick sequence of massive fluvial sandstone, bounded above and below by relatively heterogeneous, hydrologically discontinuous units and has served as a primary conduit for fluids within this stratigraphic interval. Patterns of mineral-fluid reactions suggest a basinwide hydrologic regime in which warm, evolved fluids migrated up-dip from the center of the basin under the influence of a regional hydraulic head. -from Authors

American Journal of Science↗
USGS2006

Variability and regulation of denitrification in an Upper Mississippi River backwater

Sediments in the backwaters of the Upper Mississippi River (UMR) are highly organic and provide an optimal environment for N removal. We monitored an 8.6-ha UMR backwater site near La Crosse, Wisconsin, for nearly 3 y to assess temporal variability, seasonal trends, and the factors regulating denitrification. We measured rates of unamended denitrification (DEN) and denitrification enzyme activity (DEA) rates at ambient temperature and DEA at 30 degrees C (DEA30). Seasonal mean (+/- 1 SE) DEN rates ranged from 0.041 +/- 0.015 to 0.47 +/- 0.23 mu g N cm(-2) h(-1)and were highest in winter and lowest in autumn. Seasonal rates of DEA exhibited a different pattern with the highest rates in summer (25.6 +/- 3.4 mu g N cm(-2) h(-1)) and the lowest rates in winter (10.6 +/- 2.1 mu g N cm(-2) h(-1)). The overall mean DEA30 rate was 31.0 +/- 1.9 mu g N cm(-2) h(-1) but showed no significant seasonal pattern. Short-term (weekly) and seasonal variability exhibited by rates of DEN and DEA were best explained by water-column NO3- concentration and temperature, respectively. No environmental variables explained a significant amount of variability in DEA30. Our results suggest that nutrient (i.e., NO3-) availability and temperature are both regulators of denitrification, with NO3- concentration being the most important limiting factor in this system. The high DEN rates during winter were in response to elevated NO3- concentrations resulting from a chain reaction beginning with algal blooms creating oxic conditions that stimulated nitrification. Increasing hydrological connectivity in large rivers as a river management tool to reduce N flux to downstream areas may be beneficial.

Journal of the North American Benthological Societ↗
USGS2007

Investigation of the groundwater system at Masaya Caldera, Nicaragua, using transient electromagnetics and numerical simulation

The distribution of groundwater beneath Masaya Volcano, in Nicaragua, and its surrounding caldera was characterized using the transient electromagnetic method (TEM). Multiple soundings were conducted at 30 sites. Models of the TEM data consistently indicate a resistive layer that is underlain by one or more conductive layers. These two layers represent the unsaturated and saturated zones, respectively, with the boundary between them indicating the water-table elevation. A map of the TEM data shows that the water table in the caldera is a subdued replica of the topography, with higher elevations beneath the edifice in the south-central caldera and lower elevations in the eastern caldera, coinciding with the elevation of Laguna de Masaya. These TEM data, combined with regional hydrologic data, indicate that the caldera in hydrologically isolated from the surrounding region, with as much as 60??m of difference in elevation of the groundwater table across caldera-bounding faults. The water-table information and estimates of fluxes of water through the system were used to constrain a numerical simulation of groundwater flow. The simulation results indicate that basalt flows in the outer parts of the caldera have a relatively high transmissivity, whereas the central edifice has a substantially lower transmissivity. A layer of relatively high transmissivity must be present at depth within the edifice in order to deliver the observed flux of water and steam to the active vent. This hydrologic information about the caldera provides a baseline for assessing the response of this isolated groundwater system to future changes in magmatic activity. ?? 2007.

Journal of Volcanology and Geothermal Research↗