USGS Science⌕ Search

SEARCH · USGS Science

Results for “Journal of Structural Engineering”

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 55 records · Page 3Linked to original sources

Advances in a distributed approach for ocean model data interoperability

An infrastructure for earth science data is emerging across the globe based on common data models and web services. As we evolve from custom file formats and web sites to standards-based web services and tools, data is becoming easier to distribute, find and retrieve, leaving more time for science. We describe recent advances that make it easier for ocean model providers to share their data, and for users to search, access, analyze and visualize ocean data using MATLAB® and Python®. These include a technique for modelers to create aggregated, Climate and Forecast (CF) metadata convention datasets from collections of non-standard Network Common Data Form (NetCDF) output files, the capability to remotely access data from CF-1.6-compliant NetCDF files using the Open Geospatial Consortium (OGC) Sensor Observation Service (SOS), a metadata standard for unstructured grid model output (UGRID), and tools that utilize both CF and UGRID standards to allow interoperable data search, browse and access. We use examples from the U.S. Integrated Ocean Observing System (IOOS®) Coastal and Ocean Modeling Testbed, a project in which modelers using both structured and unstructured grid model output needed to share their results, to compare their results with other models, and to compare models with observed data. The same techniques used here for ocean modeling output can be applied to atmospheric and climate model output, remote sensing data, digital terrain and bathymetric data.

Journal of Marine Science and Engineering↗

Characterization of the structural–stratigraphic and reservoir controls on the occurrence of gas hydrates in the Eileen Gas Hydrate Trend, Alaska North Slope

One of the most studied permafrost-associated gas hydrate accumulations in Arctic Alaska is the Eileen Gas Hydrate Trend. This study provides a detailed re-examination of the Eileen Gas Hydrate Trend with a focus on the gas hydrate accumulation in the western part of the Prudhoe Bay Unit. This integrated analysis of downhole well log data and published geophysical data has provided new insight on structural, stratigraphic, and reservoir controls on the occurrence of gas hydrates in the Eileen Gas Hydrate Trend. This study revealed the relatively complex nature of the gas hydrate occurrences in the Eileen Gas Hydrate Trend, with gas hydrates present in a series of coarsening upward, laterally pervasive, mostly fine-grained sand beds exhibiting high gas hydrate saturations. Most of the gas hydrate-bearing reservoirs in the Eileen Gas Hydrate Trend are laterally segmented into distinct northwest- to southeast-trending fault blocks, occur in a combination of structural–stratigraphic traps, and are only partially hydrate filled with distinct down-dip water contacts. These findings suggest that the traditional parts of a petroleum system (i.e., reservoir, gas source, gas migration, and geologic timing of the system formation) also control the occurrence of gas hydrates in the Eileen Gas Hydrate Trend.

Alaska↗

Groundwater use: Equilibrium between social benefits and potential environmental costs

In many countries groundwater resources are under-appreciated and, therefore, underutilizied; whereas, in some areas they are inappropriately exploited and, therefore, over-utilized. “Over utilization” can lead to depletion in quantity or a degradation in quality or both. Obstacles to effective management include: (1) lack of knowledge of basic principles of groundwater science among water planners, (2) in many, if not in most countries, ownership of groundwater is in the private domain with the result that codependence is unrecognized, and (3) a misunderstanding by water planners of the concepts of “overexploitation%rdquo; and conjunctive use. The economic, social, and hydrologic constraints and procedures for management for sustainable development of groundwater are significantly different from those for surface water because these differences result from such things as (1) groundwater development is not dependent on large scale collective projects (unlike the utilization of surface water that requires engineering structures for diverting, regulating and transporting water), (2) the activities of many different groups can affect the quality of water, and (3) users of groundwater often are not aware of their co-dependence on the groundwater heritage in which each participates. Hydrogeologists should try to identify those governmental policies that have a detrimental environmental effect, promote those policies that are beneficial, and demonstrate the need for a policy in matters where a policy is lacking.

Hydrogeology Journal↗

Variability in wood-frame building damage using broad-band synthetic ground motions: a comparative numerical study with recorded motions

Earthquake damage to light-frame wood buildings is a major concern for North America because of the volume of this construction type. In order to estimate wood building damage using synthetic ground motions, we need to verify the ability of synthetically generated ground motions to simulate realistic damage for this structure type. Through a calibrated damage potential indicator, four different synthetic ground motion models are compared with the historically recorded ground motions at corresponding sites. We conclude that damage for sites farther from the fault (>20 km) is under-predicted on average and damage at closer sites is sometimes over-predicted.

Journal of Earthquake Engineering↗

Shear-wave velocity-based probabilistic and deterministic assessment of seismic soil liquefaction potential

Shear-wave velocity (V s ) offers a means to determine the seismic resistance of soil to liquefaction by a fundamental soil property. This paper presents the results of an 11-year international project to gather new V s site data and develop probabilistic correlations for seismic soil liquefaction occurrence. Toward that objective, shear-wave velocity test sites were identified, and measurements made for 301 new liquefaction field case histories in China, Japan, Taiwan, Greece, and the United States over a decade. The majority of these new case histories reoccupy those previously investigated by penetration testing. These new data are combined with previously published case histories to build a global catalog of 422 case histories of V s liquefaction performance. Bayesian regression and structural reliability methods facilitate a probabilistic treatment of the V s catalog for performance-based engineering applications. Where possible, uncertainties of the variables comprising both the seismic demand and the soil capacity were estimated and included in the analysis, resulting in greatly reduced overall model uncertainty relative to previous studies. The presented data set and probabilistic analysis also help resolve the ancillary issues of adjustment for soil fines content and magnitude scaling factors.

Journal of Geotechnical and Geoenvironmental Engin↗

Mercury adsorption properties of sulfur-impregnated adsorbents

Carbonaceous and noncarbonaceous adsorbents were impregnated with elemental sulfur to evaluate the chemical and physical properties of the adsorbents and their equilibrium mercury adsorption capacities. Simulated coal combustion flue gas conditions were used to determine the equilibrium adsorption capacities for Hg0 and HgCl2 gases to better understand how to remove mercury from gas streams generated by coal-fired utility power plants. Sulfur was deposited onto the adsorbents by monolayer surface deposition or volume pore filling. Sulfur impregnation increased the total sulfur content and decreased the total and micropore surface areas and pore volumes for all of the adsorbents tested. Adsorbents with sufficient amounts of active adsorption sites and sufficient microporous structure had mercury adsorption capacities up to 4,509 ??g Hg/g adsorbent. Elemental sulfur, organic sulfur, and sulfate were formed on the adsorbents during sulfur impregnation. Correlations were established with R2>0.92 between the equilibrium Hg0/HgCl2 adsorption capacities and the mass concentrations of elemental and organic sulfur. This result indicates that elemental and organic sulfur are important active adsorption sites for Hg0 and HgCl2.

Journal of Environmental Engineering↗

Management of fluid mud in estuaries, bays, and lakes. II: Measurement, modeling, and management

Techniques for measurement, modeling, and management of fluid mud are available, but research is needed to improve them. Fluid mud can be difficult to detect, measure, or sample, which has led to new instruments and new ways of using existing instruments. Multifrequency acoustic fathometers sense neither density nor viscosity and are, therefore, unreliable in measuring fluid mud. Nuclear density probes, towed sleds, seismic, and drop probes equipped with density meters offer the potential for accurate measurements. Numerical modeling of fluid mud requires solving governing equations for flow velocity, density, pressure, salinity, water surface, plus sediment submodels. A number of such models exist in one-, two-, and three-dimensional form, but they rely on empirical relationships that require substantial site-specific validation to observations. Management of fluid mud techniques can be classified as those that accomplish: Source control, formation control, and removal. Nautical depth, a fourth category, defines the channel bottom as a specific fluid mud density or alternative parameter as safe for navigation. Source control includes watershed management measures to keep fine sediment out of waterways and in-water measures such as structures and traps. Formation control methods include streamlined channels and structures plus other measures to reduce flocculation and structures that train currents. Removal methods include the traditional dredging and transport of dredged material plus agitation that contributes to formation control and/or nautical depth. Conditioning of fluid mud by dredging and aerating offers the possibility of improved navigability. Two examples—the Atchafalaya Bar Channel and Savannah Harbor—illustrate the use of measurements and management of fluid mud.

Journal of Hydraulic Engineering↗

Mississippi River streamflow measurement techniques at St. Louis, Missouri

Streamflow measurement techniques of the Mississippi River at St. Louis have changed through time (1866–present). In addition to different methods used for discrete streamflow measurements, the density and range of discrete measurements used to define the rating curve (stage versus streamflow) have also changed. Several authors have utilized published water surface elevation (stage) and streamflow data to assess changes in the rating curve, which may be attributed to be caused by flood control and/or navigation structures. The purpose of this paper is to provide a thorough review of the available flow measurement data and techniques and to assess how a strict awareness of the limitations of the data may affect previous analyses. It is concluded that the pre-1930s discrete streamflow measurement data are not of sufficient accuracy to be compared with modern streamflow values in establishing long-term trends of river behavior.

Missouri↗

Maximum entropy derived statistics of sound speed structure in a fine-grained sediment inferred from sparse broadband acoustic measurements on the New England continental shelf

Marginal probability distributions for parameters representing an effective sound-speed structure of a fine-grained sediment are inferred from a data ensemble maximum entropy method that utilizes a sparse spatially distributed set of received pressure time series resulting from multiple explosive sources in a shallow-water ocean environment possessing significant spatial variability of the seabed. A remote sensing seabed acoustics experiment undertaken in March 2017 off the New England Shelf was designed so that multiple independent analyses could infer the statistical properties of the seabed. The current analysis incorporates the measured horizontal variability from interpretations of a subbottom profiling survey of the experimental area. An idealized range- and azimuth-dependent parameterization of the seabed is derived from identification of horizons within the seabed that define multiple sediment layers. A sparse set of explosive charges were deployed on circular tracks with radii of about 2, 4, and 6.5 km with an acoustic array at the center to correlate a set of random measurements to physical acoustic processes that characterize the seabed. The mean values of a surface sound speed ratio and a linear sound speed gradient for the fine-grained sediment layer derived from 12 data samples processed in the 25–275-Hz band provide an estimate of the effective sound-speed structure in a 130-km $^2$ area. The inferred sediment sound speed values are evaluated by predicting measured time series data not used in the statistical inference, and are also compared to historical measurements. Finally, the low-frequency maximum entropy estimate of the sediment sound speed along with physical measurements derived from piston core measurements are utilized to estimate the sediment grain bulk modulus. predictions made by the viscous grain shearing model.

Maine, Vermont, New Hampshire, Massachusetts, Rhod↗

Predicting boundary shear stress and sediment transport over bed forms

To estimate bed-load sediment transport rates in flows over bed forms such as ripples and dunes, spatially averaged velocity profiles are frequently used to predict mean boundary shear stress. However, such averaging obscures the complex, nonlinear interaction of wake decay, boundary-layer development, and topographically induced acceleration downstream of flow separation and often leads to inaccurate estimates of boundary stress, particularly skin friction, which is critically important in predicting bed-load transport rates. This paper presents an alternative methodology for predicting skin friction over 2D bed forms. The approach is based on combining the equations describing the mechanics of the internal boundary layer with semiempirical structure functions to predict the velocity at the crest of a bedform, where the flow is most similar to a uniform boundary layer. Significantly, the methodology is directed toward making specific predictions only at the bed-form crest, and as a result it avoids the difficulty and questionable validity of spatial averaging. The model provides an accurate estimate of the skin friction at the crest where transport rates are highest. Simple geometric constraints can be used to derive the mean transport rates as long as bed load is dominant.To estimate bed-load sediment transport rates in flows over bed forms such as ripples and dunes, spatially averaged velocity profiles are frequently used to predict mean boundary shear stress. However, such averaging obscures the complex, nonlinear interaction of wake decay, boundary-layer development, and topographically induced acceleration downstream of flow separation and often leads to inaccurate estimates of boundary stress, particularly skin friction, which is critically important in predicting bed-load transport rates. This paper presents an alternative methodology for predicting skin friction over 2D bed forms. The approach is based on combining the equations describing the mechanics of the internal boundary layer with semiempirical structure functions to predict the velocity at the crest of a bedform, where the flow is most similar to a uniform boundary layer. Significantly, the methodology is directed toward making specific predictions only at the bed-form crest, and as a result it avoids the difficulty and questionable validity of spatial averaging. The model provides an accurate estimate of the skin friction at the crest where transport rates are highest. Simple geometric constraints can be used to derive the mean transport rates as long as bed load is dominant.

Journal of Hydraulic Engineering↗

The application of oyster reefs in shoreline protection: Are we over‐engineering for an ecosystem engineer?

Oyster reef living shorelines have been proposed as an effective alternative to traditional coastal defence structures (e.g. bulkheads, breakwaters), with the benefit that they may keep pace with sea‐level rise and provide co‐benefits, such as habitat provision. However, there remains uncertainty about the effectiveness of shoreline protection provided by oyster reefs, which limits their broader application. We draw evidence from studies along the east and gulf coasts of the United States, where much research and implementation of oyster reef restoration has occurred, to better define the existing gaps in our understanding of the use of restored oyster reefs for shoreline protection. We find potential disconnects between ecological and engineering functions of reefs. In response, we outline how engineering and ecological principles are used in the design of oyster reef living shorelines and highlight knowledge gaps where an integration of these disciplines will lead to their more effective application. Synthesis and applications . This work highlights the necessary steps to advance the application of oyster reef living shorelines. Importantly, future research should focus on appropriate designs and conditions needed for these structures to effectively protect our coasts from erosion, while supporting a sustainable oyster population, thereby providing actionable nature‐based alternatives for coastal defence to diverse end‐users.

East Coast of Continental United States↗

Extraction of near-surface properties for a lossy layered medium using the propagator matrix

Near-surface properties play an important role in advancing earthquake hazard assessment. Other areas where near-surface properties are crucial include civil engineering and detection and delineation of potable groundwater. From an exploration point of view, near-surface properties are needed for wavefield separation and correcting for the local near-receiver structure. It has been shown that these properties can be estimated for a lossless homogeneous medium using the propagator matrix. To estimate the near-surface properties, we apply deconvolution to passive borehole recordings of waves excited by an earthquake. Deconvolution of these incoherent waveforms recorded by the sensors at different depths in the borehole with the recording at the surface results in waves that propagate upwards and downwards along the array. These waves, obtained by deconvolution, can be used to estimate the P- and S-wave velocities near the surface. As opposed to waves obtained by cross-correlation that represent filtered version of the sum of causal and acausal Green's function between the two receivers, the waves obtained by deconvolution represent the elements of the propagator matrix. Finally, we show analytically the extension of the propagator matrix analysis to a lossy layered medium for a special case of normal incidence. ?? 2007 The Authors Journal compilation ?? 2007 RAS.

Geophysical Journal International↗

Hyperspectral narrowband imaging spectroscopy: A new paradigm for Earth observation

This editorial introduces the Special Issue entitled “Hyperspectral Narrowband Imaging Spectroscopy: A New Paradigm for Earth Observation” in the August 2026 issue of Photogrammetric Engineering & Remote Sensing (PE&RS), the flagship journal of the American Society for Photogrammetry and Remote Sensing (ASPRS). This volume represents the fourth dedicated hyperspectral special issue published in PE&RS, following earlier contributions by Thenkabail et al. (2025, 2024a, 2024b), and continues ASPRS’s commitment to advancing cutting‑edge imaging spectroscopy research and its applications across Earth system science. Remote sensing is undergoing one of the most profound transformations in its history. The emergence of hyperspectral narrowband (HNB) imaging spectroscopy data, capable of acquiring hundreds of contiguous, narrow spectral bands, has shifted the discipline from observing Earth in a handful of broad spectral windows to capturing continuous spectral signatures of the Planet. This transition marks a decisive break from the multispectral paradigm that has dominated satellite remote sensing for nearly five decades, driven by the advent of new orbital imaging spectrometers such as EnMAP, PRISMA, and NASA’s EMIT, and by the forthcoming Surface Biology and Geology (SBG) mission ((Pires Silva et al., 2026; Bourriz et al., 2025; Thenkabail et al., 2025; Chabrillat et al., 2024; Aneece et al., 2024; Dave et al., 2024; Thenkabail et al., 2024a; Thenkabail et al., 2024b; Thenkabail, 2024a; Thenkabail, 2024b; Thompson et al., 2022; Kokaly et al., 2022; Aneece & Thenkabail, 2022; Cawse Nicholson et al., 2021; Guanter et al., 2021; Vangi et al., 2021; Thenkabail et al., 2021). These missions (e.g., Table 1) deliver unprecedented spectral fidelity, improved signal to noise ratios, and global coverage capabilities, enabling a new era of quantitative, spectroscopy based Earth observation. Where multispectral broadbands (MBBs) provide only a few discrete measurements along the electromagnetic spectrum, HNB systems deliver rich, diagnostic information that enables scientists to characterize Earth’s surface with unprecedented biochemical, biophysical, and structural detail (Figure 1a, 1b). The implications for environmental monitoring, agriculture, water resources, and mineral exploration are profound. Several overarching themes emerge: • Spectral fidelity matters. The ability to preserve subtle absorption features is essential for mineral mapping, vegetation trait retrieval, and biochemical modeling. • AI and deep learning are indispensable. From destriping to classification, modern analytics must be scalable, label‑efficient, and capable of exploiting the full spectral–spatial richness of HNB data. • Physics‑based and data‑driven approaches must converge. Radiative transfer models such as PROSAIL, enhanced with localized soil parameterizations, remain foundational for biophysical retrievals and model‑based inference. • Dimensionality reduction and feature extraction are critical. Techniques such as L1‑ISOMAP demonstrate that intelligent manifold learning can unlock the structure of fused, high‑dimensional datasets. • Next‑generation architectures must be interactive and multimodal. ICTNet exemplifies the future of hyperspectral classification: hybrid, synergistic, and capable of modeling both local textures and global spectral dependencies.

Photogrammetric Engineering and Remote Sensing (PE↗

The Eastern Gas Shales Project (EGSP) Data System: A case study in data base design, development, and application

The U.S. Geological Survey and Petroleum Information Corporation in Denver, Colorado, developed the Eastern Gas Shale Project (EGSP)Data System for the U.S. Department of Energy, Morgantown, West Virginia. Geological, geochemical, geophysical, and engineering data from Devonian shale samples from more than 5800 wells and outcrops in the Appalachian basin were edited and converted to a Petroleum Information Corporation data base. Well and sample data may be retrieved from this data system to produce (1)production-test summaries by formation and well location; (2)contoured isopach, structure, and trendsurface maps of Devonian shale units; (3)sample summary reports for samples by location, well, contractor, and sample number; (4)cross sections displaying digitized log traces, geochemical, and lithologic data by depth for wells; and (5)frequency distributions and bivariate plots. Although part of the EGSP Data System is proprietary, and distribution of complete well histories is prohibited by contract, maps and aggregated well-data listings are being made available to the public through published reports. ?? 1983 Plenum Publishing Corporation.

Journal of the International Association for Mathe↗

Assessment of modal-pushover-based scaling procedure for nonlinear response history analysis of ordinary standard bridges

The earthquake engineering profession is increasingly utilizing nonlinear response history analyses (RHA) to evaluate seismic performance of existing structures and proposed designs of new structures. One of the main ingredients of nonlinear RHA is a set of ground motion records representing the expected hazard environment for the structure. When recorded motions do not exist (as is the case in the central United States) or when high-intensity records are needed (as is the case in San Francisco and Los Angeles), ground motions from other tectonically similar regions need to be selected and scaled. The modal-pushover-based scaling (MPS) procedure was recently developed to determine scale factors for a small number of records such that the scaled records provide accurate and efficient estimates of “true” median structural responses. The adjective “accurate” refers to the discrepancy between the benchmark responses and those computed from the MPS procedure. The adjective “efficient” refers to the record-to-record variability of responses. In this paper, the accuracy and efficiency of the MPS procedure are evaluated by applying it to four types of existing Ordinary Standard bridges typical of reinforced concrete bridge construction in California. These bridges are the single-bent overpass, multi-span bridge, curved bridge, and skew bridge. As compared with benchmark analyses of unscaled records using a larger catalog of ground motions, it is demonstrated that the MPS procedure provided an accurate estimate of the engineering demand parameters (EDPs) accompanied by significantly reduced record-to-record variability of the EDPs. Thus, it is a useful tool for scaling ground motions as input to nonlinear RHAs of Ordinary Standard bridges.

Journal of Bridge Engineering↗

Assessing power of large river fish monitoring programs to detect population changes: the Missouri River sturgeon example

In 2003, the US Army Corps of Engineers initiated the Pallid Sturgeon Population Assessment Program (PSPAP) to monitor pallid sturgeon and the fish community of the Missouri River. The power analysis of PSPAP presented here was conducted to guide sampling design and effort decisions. The PSPAP sampling design has a nested structure with multiple gear subsamples within a river bend. Power analyses were based on a normal linear mixed model, using a mixed cell means approach, with variance estimates from the original data. It was found that, at current effort levels, at least 20 years for pallid and 10 years for shovelnose sturgeon is needed to detect a 5% annual decline. Modified bootstrap simulations suggest power estimates from the original data are conservative due to excessive zero fish counts. In general, the approach presented is applicable to a wide array of animal monitoring programs.

Journal of Applied Ichthyology↗

Framework for mapping liquefaction hazard–Targeted design ground motions

Liquefaction-induced ground failure poses substantial challenges to geotechnical earthquake engineering design. Current approaches for designing against liquefaction hazards, as specified in most seismic provisions, focus on estimating a liquefaction factor of safety ( 𝐹⁢𝑆𝐿 ) and typically characterize earthquake loading using design parameters based on probabilistic or deterministic ground motion levels. Because 𝐹⁢𝑆𝐿 is estimated deterministically, this basis of design neglects considerable uncertainties for estimating liquefaction triggering and its consequences and results in a lack of liquefaction-specific design criteria, particularly as structural design has advanced toward risk-targeted performance objectives. This study presents a framework for developing liquefaction-targeted design criteria based on a minimum acceptable return period of liquefaction, informed by probabilistic liquefaction hazard analysis (PLHA). PLHA quantifies annualized rates of liquefaction by considering contributions from (1) the full ground-motion probability space, and (2) uncertainties in liquefaction triggering using probabilistic models. PLHA is used in this study to characterize the current, effective return periods of 𝐹⁢𝑆𝐿 ( 𝑇𝑅,𝐹⁢𝑆 ) obtained from conventional liquefaction hazard analysis (CLHA) using uniform-hazard ground motions. 𝑇𝑅,𝐹⁢𝑆 is evaluated in a parametric study of nearly 100 sites throughout the conterminous United States. The results indicate large geographic variations in acceptable liquefaction hazard levels, with implied 𝑇𝑅,𝐹⁢𝑆 ranging between approximately 1,000 to 3,000 years. To address these inconsistencies without the computational demands of full PLHA, a framework is proposed for developing a liquefaction-targeted design peak ground acceleration, 𝑃⁢𝐺⁢𝐴𝐿 , for use in liquefaction models that result in consistent liquefaction design levels across all geographic locations. The mapped 𝑃⁢𝐺⁢𝐴𝐿 is shown to be somewhat sensitive to site-specific properties, and adjustment factors are developed and presented. The proposed 𝑃⁢𝐺⁢𝐴𝐿 mapping procedure produces 𝐹⁢𝑆𝐿 estimates that are consistent with those obtained from full PLHA at a target 𝑇𝑅,𝐹⁢𝑆 , providing a promising roadmap to incorporating PLHA concepts into current liquefaction design methods.

Journal of Geotechnical and Geoenvironmental Engin↗

Reclaimed mineland curve number response to temporal distribution of rainfall

The curve number (CN) method is a common technique to estimate runoff volume, and it is widely used in coal mining operations such as those in the Appalachian region of Kentucky. However, very little CN data are available for watersheds disturbed by surface mining and then reclaimed using traditional techniques. Furthermore, as the CN method does not readily account for variations in infiltration rates due to varying rainfall distributions, the selection of a single CN value to encompass all temporal rainfall distributions could lead engineers to substantially under- or over-size water detention structures used in mining operations or other land uses such as development. Using rainfall and runoff data from a surface coal mine located in the Cumberland Plateau of eastern Kentucky, CNs were computed for conventionally reclaimed lands. The effects of temporal rainfall distributions on CNs was also examined by classifying storms as intense, steady, multi-interval intense, or multi-interval steady. Results indicate that CNs for such reclaimed lands ranged from 62 to 94 with a mean value of 85. Temporal rainfall distributions were also shown to significantly affect CN values with intense storms having significantly higher CNs than multi-interval storms. These results indicate that a period of recovery is present between rainfall bursts of a multi-interval storm that allows depressional storage and infiltration rates to rebound. ?? 2010 American Water Resources Association.

Journal of the American Water Resources Associatio↗