USGS Science⌕ Search

SEARCH · USGS Science

Results for “Interfaces”

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 937 records · Page 52Linked to original sources
USGS2004

Crustal deformation measurements in Guerrero, Mexico

GPS measurements of crustal deformation in Guerrero, southern Mexico, include surveys collected between 1992 and 2001 as well as continuous GPS measurements at a few sites. These geodetic observations are used to calculate interseismic deformation rates and assess the presence and possible location of transient deformation during the period encompassing 1992.25 to 2001.75. The data are used to examine transient deformation in 1998 previously described from data at a single site by Lowry et al. [2001]. Survey measurements and continuous data from a site near Popocate??petl volcano confirm the 1998 transient, and survey data also suggest another transient occurred following the 14 September 1995 (Mw = 7.3) Copala earthquake. All of the available GPS position estimates have been inverted for a combined model of slip during each event plus the steady state slip on the plate interface. Modeling of the steady state deformation rates confirms that the Guerrero seismic gap is partially frictionally locked at depths shallower than about 25 km and accumulating strain that may eventually be released in a great earthquake. The data also suggest that there is frictional coupling to much greater (>40 km) depths, which releases more frequently in aseismic slip events. The locations and sizes of the transient events are only partially constrained by the available data. However, the transient models which best fit the GPS coordinate time series suggest that aseismic slip was centered downdip of the seismogenic portion of the plate-bounding thrust in both events, and the moment release had equivalent magnitudes Mw = 7.1 + 1.3/-1.0 in 1995-1996 and 7.1 + 0.4/-0.1 in 1998. Copyright 2004 by the American Geophysical Union.

Journal of Geophysical Research B: Solid Earth↗
USGS2011

Proposed initiative would study Earth's weathering engine

At the Earth's surface, a complex suite of chemical, biological, and physical processes combines to create the engine that transforms bedrock into soil (Figure 1). Earth's weathering engine provides nutrients to nourish ecosystems and human society mediates the transport of toxic components within the biosphere, creates water flow paths that carve and weaken bedrock, and contributes to the evolution of landscapes at all temporal and spatial scales. At the longest time scales, the weathering engine sequesters CO 2 , thereby influencing long-term climate change. Despite the importance of soil, our knowledge of the rate of soil formation is limited because the weathering zone forms a complex, ever-changing interface, and because scientific approaches and funding paradigms have not promoted integrated research agendas to investigate such complex interactions. No national initiative has promoted a systems approach to investigation of weathering science across the broad array of geology, soil science, ecology and hydrology. Such a program is certainly needed, and this article describes a platform on which to build the initiative to answer the following question: How does the Earth weathering engine break down rock to nourish ecosystems, carve terrestrial landscapes, and control carbon dioxide in the global atmosphere?

Eos Science News↗
USGS2005

The wrinkle-like slip pulse is not important in earthquake dynamics

A particular solution for slip on an interface between different elastic materials, the wrinkle-like slip pulse, propagates in only one direction with reduced normal compressive stress. More general solutions, and natural earthquakes, need not share those properties. In a 3D dynamic model with a drop in friction and heterogeneous initial stress, the wrinkle-like slip pulse is only a small part of the solution. Rupture propagation is determined primarily by the potential stress drop, not by the wrinkle-like slip pulse. A 2D calculation with much finer resolution shows that energy loss to friction might not be significantly reduced in the wrinkle-like slip pulse. Copyright 2005 by the American Geophysical Union.

Geophysical Research Letters↗
USGS2005

Coseismic slip distribution of the 1923 Kanto earthquake, Japan

The slip distribution associated with the 1923 M = 7.9 Kanto, Japan, earthquake is reexamined in light of new data and modeling. We utilize a combination of first-order triangulation, second-order triangulation, and leveling data in order to constrain the coseismic deformation. The second-order triangulation data, which have not been utilized in previous studies of 1923 coseismic deformation, are associated with only slightly smaller errors than the first-order triangulation data and expand the available triangulation data set by about a factor of 10. Interpretation of these data in terms of uniform-slip models in a companion study by Nyst et al. shows that a model involving uniform coseismic slip on two distinct rupture planes explains the data very well and matches or exceeds the fit obtained by previous studies, even one which involved distributed slip. Using the geometry of the Nyst et al. two-plane slip model, we perform inversions of the same geodetic data set for distributed slip. Our preferred model of distributed slip on the Philippine Sea plate interface has a moment magnitude of 7.86. We find slip maxima of ???8-9 m beneath Odawara and ???7-8 m beneath the Miura peninsula, with a roughly 2:1 ratio of strike-slip to dip-slip motion, in agreement with a previous study. However, the Miura slip maximum is imaged as a more broadly extended feature in our study, with the high-slip region continuing from the Miura peninsula to the southern Boso peninsula region. The second-order triangulation data provide good evidence for ???3 m right-lateral strike slip on a 35-km-long splay structure occupying the volume between the upper surface of the descending Philippine Sea plate and the southern Boso peninsula. Copyright 2005 by the American Geophysical Union.

Journal of Geophysical Research B: Solid Earth↗
USGS2006

Seismic characteristics of central Brazil crust and upper mantle: A deep seismic refraction study

[1] A two‐dimensional model of the Brazilian central crust and upper mantle was obtained from the traveltime interpretation of deep seismic refraction data from the Porangatu and Cavalcante lines, each approximately 300 km long. When the lines were deployed, they overlapped by 50 km, forming an E‐W transect approximately 530 km long across the Tocantins Province and western São Francisco Craton. The Tocantins Province formed during the Neoproterozoic when the São Francisco, the Paranapanema, and the Amazon cratons collided, following the subduction of the former Goiás ocean basin. Average crustal V P and V P /V S ratios, Moho topography, and lateral discontinuities within crustal layers suggest that the crust beneath central Brazil can be associated with major geological domains recognized at the surface. The Moho is an irregular interface, between 36 and 44 km deep, that shows evidences of first‐order tectonic structures. The 8.05 and 8.23 km s −1 P wave velocities identify the upper mantle beneath the Porangatu and Cavalcante lines, respectively. The observed seismic features allow for the identification of (1) the crust has largely felsic composition in the studied region, (2) the absence of the mafic‐ultramafic root beneath the Goiás magmatic arc, and (3) block tectonics in the foreland fold‐and‐thrust belt of the northern Brasília Belt during the Neoproterozoic. Seismic data also suggested that the Bouguer gravimetric discontinuities are mainly compensated by differences in mass distribution within the lithospheric mantle. Finally, the Goiás‐Tocantins seismic belt can be interpreted as a natural seismic alignment related to the Neoproterozoic mantle domain.

Journal of Geophysical Research B: Solid Earth↗
USGS2006

Impacts of the 2004 tsunami on groundwater resources in Sri Lanka

The 26 December 2004 tsunami caused widespread destruction and contamination of coastal aquifers across southern Asia. Seawater filled domestic open dug wells and also entered the aquifers via direct infiltration during the first flooding waves and later as ponded seawater infiltrated through the permeable sands that are typical of coastal aquifers. In Sri Lanka alone, it is estimated that over 40,000 drinking water wells were either destroyed or contaminated. From February through September 2005, a team of United States, Sri Lankan, and Danish water resource scientists and engineers surveyed the coastal groundwater resources of Sri Lanka to develop an understanding of the impacts of the tsunami and to provide recommendations for the future of coastal water resources in south Asia. In the tsunami‐affected areas, seawater was found to have infiltrated and mixed with fresh groundwater lenses as indicated by the elevated groundwater salinity levels. Seawater infiltrated through the shallow vadose zone as well as entered aquifers directly through flooded open wells. Our preliminary transport analysis demonstrates that the intruded seawater has vertically mixed in the aquifers because of both forced and free convection. Widespread pumping of wells to remove seawater was effective in some areas, but overpumping has led to upconing of the saltwater interface and rising salinity. We estimate that groundwater recharge from several monsoon seasons will reduce salinity of many sandy Sri Lankan coastal aquifers. However, the continued sustainability of these small and fragile aquifers for potable water will be difficult because of the rapid growth of human activities that results in more intensive groundwater pumping and increased pollution. Long‐term sustainability of coastal aquifers is also impacted by the decrease in sand replenishment of the beaches due to sand mining and erosion.

Water Resources Research↗
USGS2007

Anomalous deepening of a belt of intraslab earthquakes in the Pacific slab crust under Kanto, central Japan: Possible anomalous thermal shielding, dehydration reactions, and seismicity caused by shallower cold slab material

[1] A belt of intraslab seismicity in the Pacific slab crust parallel to iso-depth contours of the plate interface has been found beneath Hokkaido and Tohoku. Hypocenter relocations have shown that this seismic belt does not run parallel to but obliquely to the iso-depth contours beneath Kanto, deepening toward the north from ∼100 km to ∼140 km depth. The depth limit of the contact zone with the overlying Philippine Sea slab is located close to and parallel to this obliquely oriented seismic belt, suggesting that the deepening of the seismic belt there is caused by the contact with the overlying slab. The contact with this cold slab hinders the heating of the Pacific slab crust by hot mantle wedge, which would cause delay of eclogite-forming phase transformations and hence deepening of the seismic belt there. The depth limit of the subducting low-velocity crust also deepens toward the north, supporting this idea.

Geophysical Research Letters↗
USGS2008

Not all water becomes wine: Sulfur inputs as an opportune tracer of hydrochemical losses from vineyards

California's widespread and economically important vineyards offer substantial opportunities to understand the interface between hydrology and biogeochemistry in agricultural soils. The common use of native sulfur (S) as a fumigant or soil additive provides a novel way to isotopically differentiate among sulfate (SO 4 2− ) pools, allowing the estimation of water and SO 4 2− budgets. The objectives of this study were (1) to characterize the near‐surface hydrological flow paths in a vineyard during irrigation and storm events and (2) to determine how those flow paths affect the fate and transport of SO 4 2− across seasons. Integrating hydrological theory with measurements of SO 4 2− concentration and sulfate‐S isotopic ratios (expressed as [SO 4 2− ] and δ 34 S, respectively) in inputs, soil water, and leachate provided a means of determining flow paths. Low [SO 4 2− ] and δ 34 S in leachate during 4‐h irrigation events reflect minimal engagement of the soil matrix, indicating that preferential flow was the dominant path for water in the near surface. In contrast, high [SO 4 2− ] and δ 34 S values during 8‐h irrigation and storm events reflect near‐complete engagement of the soil matrix, indicating that lateral flow was the dominant pathway. Because hydrologic response and SO 4 2− mobility are tightly coupled in these soils, the magnitude of water fluxes through the near surface controls S cycling both on and off site. These results indicate that preferential flow is an important loss pathway to consider in managing both water resources and water quality (reactive elements) in vineyard land use systems.

Water Resources Research↗
USGS2008

An autocorrelation method to detect low frequency earthquakes within tremor

Recent studies have shown that deep tremor in the Nankai Trough under western Shikoku consists of a swarm of low frequency earthquakes (LFEs) that occur as slow shear slip on the down-dip extension of the primary seismogenic zone of the plate interface. The similarity of tremor in other locations suggests a similar mechanism, but the absence of cataloged low frequency earthquakes prevents a similar analysis. In this study, we develop a method for identifying LFEs within tremor. The method employs a matched-filter algorithm, similar to the technique used to infer that tremor in parts of Shikoku is comprised of LFEs; however, in this case we do not assume the origin times or locations of any LFEs a priori. We search for LFEs using the running autocorrelation of tremor waveforms for 6 Hi-Net stations in the vicinity of the tremor source. Time lags showing strong similarity in the autocorrelation represent either repeats, or near repeats, of LFEs within the tremor. We test the method on an hour of Hi-Net recordings of tremor and demonstrates that it extracts both known and previously unidentified LFEs. Once identified, we cross correlate waveforms to measure relative arrival times and locate the LFEs. The results are able to explain most of the tremor as a swarm of LFEs and the locations of newly identified events appear to fill a gap in the spatial distribution of known LFEs. This method should allow us to extend the analysis of Shelly et al. (2007a) to parts of the Nankai Trough in Shikoku that have sparse LFE coverage, and may also allow us to extend our analysis to other regions that experience deep tremor, but where LFEs have not yet been identified. Copyright 2008 by the American Geophysical Union.

Geophysical Research Letters↗
USGS2009

Three‐dimensional model of Hellenic Arc deformation and origin of the Cretan uplift

[1] The Hellenic Arc of Greece is the most seismically active part of Europe, but little is know about its mechanics. We modeled deformation along the arc using a finite element model. The model was intended to capture large‐scale 3‐D structure of Nubian plate subduction beneath the Aegean block and its deformational consequences. The shape of the interface was developed using mapped traces at the surface and earthquake hypocenters at depth. Model block motions were constrained by recent compilations of GPS velocity vectors. We simulated a 10 ka period of convergence between Nubia and the Aegean and calculated the strain field in the overriding plate as well as the spatial distribution and orientation of differential stress (∣ σ 1 − σ 3 ∣). From these calculations we derived testable quantities such as the expected seismic moment rate on the interplate contact, uplift pattern, and distribution of strain modes. Our relatively simple model broadly reproduced observed uplift patterns, earthquake activity, and loci of extension and contraction. The model showed a localization of uplift near the island of Crete, where the fastest Aegean uplift rates are well documented. Comparison of calculated expected seismic moment and observed earthquake catalogs implies a nearly fully coupled interplate contact. On the basis of our modeling results, we suggest that south Aegean deformation is driven primarily by the fast moving (∼33 mm a −1 ) Aegean upper plate overriding a nearly stalled (∼5 mm a −1 ) Nubian lower plate. This tectonic setting thus more closely resembles a continental thrust than it does a typical oceanic subduction zone.

Aegean Sea, Hellenic Arc↗
USGS2008

Scaling hyporheic exchange and its influence on biogeochemical reactions in aquatic ecosystems

Hyporheic exchange and biogeochemical reactions are difficult to quantify because of the range in fluid‐flow and sediment conditions inherent to streams, wetlands, and nearshore marine ecosystems. Field measurements of biogeochemical reactions in aquatic systems are impeded by the difficulty of measuring hyporheic flow simultaneously with chemical gradients in sediments. Simplified models of hyporheic exchange have been developed using Darcy's law generated by flow and bed topography at the sediment‐water interface. However, many modes of transport are potentially involved (molecular diffusion, bioturbation, advection, shear, bed mobility, and turbulence) with even simple models being difficult to apply in complex natural systems characterized by variable sediment sizes and irregular bed geometries. In this study, we synthesize information from published hyporheic exchange investigations to develop a scaling relationship for estimating mass transfer in near‐surface sediments across a range in fluid‐flow and sediment conditions. Net hyporheic exchange was quantified using an effective diffusion coefficient ( D e ) that integrates all of the various transport processes that occur simultaneously in sediments, and dimensional analysis was used to scale D e to shear stress velocity, roughness height, and permeability that describe fluid‐flow and sediment characteristics. We demonstrated the value of the derived scaling relationship by using it to quantify dissolved oxygen (DO) uptake rates on the basis of DO profiles in sediments and compared them to independent flux measurements. The results support a broad application of the D e scaling relationship for quantifying coupled hyporheic exchange and biogeochemical reaction rates in streams and other aquatic ecosystems characterized by complex fluid‐flow and sediment conditions.

Water Resources Research↗
USGS2010

Cyclic ground tilt associated with the 2004–2008 eruption of Mount St. Helens

The long‐term behavior of the 2004–2008 effusive eruption of Mount St. Helens was characterized by a gradual decline in the rates of seismicity, dome growth, and broad‐scale ground deformation, but shallow near‐periodic “drumbeat” earthquakes over timescales of minutes indicated episodic short‐term behavior. In part to better characterize this behavior and any associated ground deformation, a network of tiltmeters was installed and recorded thousands of cyclic tilt events within the crater. The duration of these events, from minutes to hours, was too long for them to be directly caused by the drumbeat seismicity. Tilt events were characterized by a recoverable, asymmetric pattern of rapid tilt away from the vent followed by a more gradual reversal, were highly correlated between different tiltmeters in the crater, sometimes occurred in association with volcanic seismicity, and ceased at the end of the eruption. Tilt vectors converged on a point just south of the center of the preexisting 1980s lava dome, and the absence of detectable tilt outside the crater suggests a shallow source (<1 km). We examine several models, including cycles of conduit pressurization and plug slip or gas loss and stick‐slip behavior on the interface between the ascending plug and the 1980s lava dome. The small number of stations within the crater prevents a unique determination of source type or geometry, but results are consistent with a mechanism involving extrusion of the semisolid dacite plug and/or cycles of conduit pressurization.

Washington↗
USGS2011

Flow structures and sandbar dynamics in a canyon river during a controlled flood, Colorado River, Arizona

In canyon rivers, debris fan constrictions create rapids and downstream pools characterized by secondary flow structures that are closely linked to channel morphology. In this paper we describe detailed measurements of the three‐dimensional flow structure and sandbar dynamics of two pools along the Colorado River in the Grand Canyon during a controlled flood release from Glen Canyon Dam. Results indicate that the pools are characterized by large lateral recirculation zones (eddies) resulting from flow separation downstream from the channel constrictions, as well as helical flow structures in the main channel and eddy. The lateral recirculation zones are low‐velocity areas conducive to fine sediment deposition, particularly in the vicinity of the separation and reattachment points and are thus the dominant flow structures controlling sandbar dynamics. The helical flow structures also affect morphology but appear secondary in importance to the lateral eddies. During the controlled flood, sandbars in the separation and reattachment zones at both sites tended to build gradually during the rising limb and peak flow. Deposition in shallow water on the sandbars was accompanied by erosion in deeper water along the sandbar slope at the interface with the main channel. Erosion occurred via rapid mass failures as well as by gradual boundary shear stress driven processes. The flow structures and morphologic links at our study sites are similar to those identified in other river environments, in particular sharply curved meanders and channel confluences where the coexistence of lateral recirculation and helical flows has been documented.

Journal of Geophysical Research F: Earth Surface↗
USGS2011

Mg-spinel lithology: A new rock type on the lunar farside

High-resolution compositional data from Moon Mineralogy Mapper (M 3) for the Moscoviense region on the lunar farside reveal three unusual, but distinctive, rock types along the inner basin ring. These are designated "OOS" since they are dominated by high concentrations of orthopyroxene, olivine, and Mg-rich spinel, respectively. The OOS occur as small areas, each a few kilometers in size, that are widely separated within the highly feldspathic setting of the basin rim. Although the abundance of plagioclase is not well constrained within the OOS, the mafic mineral content is exceptionally high, and two of the rock types could approach pyroxenite and harzburgite in composition. The third is a new rock type identified on the Moon that is dominated by Mg-rich spinel with no other mafic minerals detectable (<5% pyroxene, olivine). All OOS surfaces are old and undisturbed since basin formation. They are effectively invisible in image data and are only recognized by their distinctive composition identified spectroscopically. The origin of these unusual lithologies appears to be linked to one or more magmatic intrusions into the lower crust, perhaps near the crust-mantle interface. Processes such as fractional crystallization and gravity settling within such intrusions may provide a mechanism for concentrating the mafic components within zones several kilometers in dimension. The OOS are embedded within highly anorthositic material from the lunar crust; they may thus be near contemporaneous with crustal products from the cooling magma ocean. Copyright ?? 2011 by the American Geophysical Union.

Journal of Geophysical Research E: Planets↗
USGS2010

Use of electrical imaging and distributed temperature sensing methods to characterize surface water–groundwater exchange regulating uranium transport at the Hanford 300 Area, Washington

We explored the use of continuous waterborne electrical imaging (CWEI), in conjunction with fiber‐optic distributed temperature sensor (FO‐DTS) monitoring, to improve the conceptual model for uranium transport within the Columbia River corridor at the Hanford 300 Area, Washington. We first inverted resistivity and induced polarization CWEI data sets for distributions of electrical resistivity and polarizability, from which the spatial complexity of the primary hydrogeologic units was reconstructed. Variations in the depth to the interface between the overlying coarse‐grained, high‐permeability Hanford Formation and the underlying finer‐grained, less permeable Ringold Formation, an important contact that limits vertical migration of contaminants, were resolved along ∼3 km of the river corridor centered on the 300 Area. Polarizability images were translated into lithologic images using established relationships between polarizability and surface area normalized to pore volume ( S por ). The FO‐DTS data recorded along 1.5 km of cable with a 1 m spatial resolution and 5 min sampling interval revealed subreaches showing (1) temperature anomalies (relatively warm in winter and cool in summer) and (2) a strong correlation between temperature and river stage (negative in winter and positive in summer), both indicative of reaches of enhanced surface water–groundwater exchange. The FO‐DTS data sets confirm the hydrologic significance of the variability identified in the CWEI and reveal a pattern of highly focused exchange, concentrated at springs where the Hanford Formation is thickest. Our findings illustrate how the combination of CWEI and FO‐DTS technologies can characterize surface water–groundwater exchange in a complex, coupled river‐aquifer system.

Washington↗
USGS2012

Slab1.0: A three-dimensional model of global subduction zone geometries

We describe and present a new model of global subduction zone geometries, called Slab1.0. An extension of previous efforts to constrain the two-dimensional non-planar geometry of subduction zones around the focus of large earthquakes, Slab1.0 describes the detailed, non-planar, three-dimensional geometry of approximately 85% of subduction zones worldwide. While the model focuses on the detailed form of each slab from their trenches through the seismogenic zone, where it combines data sets from active source and passive seismology, it also continues to the limits of their seismic extent in the upper-mid mantle, providing a uniform approach to the definition of the entire seismically active slab geometry. Examples are shown for two well-constrained global locations; models for many other regions are available and can be freely downloaded in several formats from our new Slab1.0 website, http://on.doi.gov/d9ARbS . We describe improvements in our two-dimensional geometry constraint inversion, including the use of ‘average’ active source seismic data profiles in the shallow trench regions where data are otherwise lacking, derived from the interpolation between other active source seismic data along-strike in the same subduction zone. We include several analyses of the uncertainty and robustness of our three-dimensional interpolation methods. In addition, we use the filtered, subduction-related earthquake data sets compiled to build Slab1.0 in a reassessment of previous analyses of the deep limit of the thrust interface seismogenic zone for all subduction zones included in our global model thus far, concluding that the width of these seismogenic zones is on average 30% larger than previous studies have suggested.

Journal of Geophysical Research B: Solid Earth↗
USGS2012

Sediment entrainment by debris flows: In situ measurements from the headwaters of a steep catchment

Debris flows can dramatically increase their volume, and hence their destructive potential, by entraining sediment. Yet quantitative constraints on rates and mechanics of sediment entrainment by debris flows are limited. Using an in situ sensor network in the headwaters of a natural catchment we measured flow and bed properties during six erosive debris-flow events. Despite similar flow properties and thicknesses of bed sediment entrained across all events, time-averaged entrainment rates were significantly faster for bed sediment that was saturated prior to flow arrival compared with rates for sediment that was dry. Bed sediment was entrained from the sediment-surface downward in a progressive fashion and occurred during passage of dense granular fronts as well as water-rich, inter-surge flow. En masse failure of bed sediment along the sediment-bedrock interface was never observed. Large-magnitude, high-frequency fluctuations in total normal basal stress were dissipated within the upper 5 cm of bed sediment. Within this near surface layer, concomitant fluctuations in Coulomb frictional resistance are expected, irrespective of the influence of pore fluid pressure or fluctuations in shear stress. If the near-surface sediment was wet as it was overridden by a flow, additional large-magnitude, high-frequency pore pressure fluctuations were measured in the near-surface bed sediment. These pore pressure fluctuations propagated to depth at subsonic rates and in a diffusive manner. The depth to which large excess pore pressures propagated was typically less than 10 cm, but scaled as ( D / f i ) 0.5 , in which D is the hydraulic diffusivity and f i is the frequency of a particular pore pressure fluctuation. Shallow penetration depths of granular-normal-stress fluctuations and excess pore pressures demonstrate that only near-surface bed sediment experiences the full dynamic range of effective-stress fluctuations, and as a result, can be more easily entrained than deeper sediment. These data provide robust tests for mechanical models of entrainment and demonstrate that a debris flow over wet bed sediment will be larger than the same flow over dry bed sediment.

Colorado↗
USGS2012

Hydrologic conditions controlling runoff generation immediately after wildfire

We investigated the control of postwildfire runoff by physical and hydraulic properties of soil, hydrologic states, and an ash layer immediately following wildfire. The field site is within the area burned by the 2010 Fourmile Canyon Fire in Colorado, USA. Physical and hydraulic property characterization included ash thickness, particle size distribution, hydraulic conductivity, and soil water retention curves. Soil water content and matric potential were measured indirectly at several depths below the soil surface to document hydrologic states underneath the ash layer in the unsaturated zone, whereas precipitation and surface runoff were measured directly. Measurements of soil water content showed that almost no water infiltrated below the ash layer into the near-surface soil in the burned site at the storm time scale (i.e., minutes to hours). Runoff generation processes were controlled by and highly sensitive to ash thickness and ash hydraulic properties. The ash layer stored from 97% to 99% of rainfall, which was critical for reducing runoff amounts. The hydrologic response to two rain storms with different rainfall amounts, rainfall intensity, and durations, only ten days apart, indicated that runoff generation was predominantly by the saturation-excess mechanism perched at the ash-soil interface during the first storm and predominantly by the infiltration-excess mechanism at the ash surface during the second storm. Contributing area was not static for the two storms and was 4% (saturation excess) to 68% (infiltration excess) of the catchment area. Our results showed the importance of including hydrologic conditions and hydraulic properties of the ash layer in postwildfire runoff generation models.

Colorado↗