Search USGSSearch

SEARCH · Search USGS

Results for “Mathematical and Computer Modelling”

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 109 records · Page 6Linked to original sources

Optimal exploitation strategies for an animal population in a stochastic serially correlated environment

Optimal exploitation strategies were studied for an animal population in a stochastic, serially correlated environment. This is a general case and encompasses a number of important cases as simplifications. Data on the mallard (Anas platyrhynchos) were used to explore the exploitation strategies and test several hypotheses because relatively much is known concerning the life history and general ecology of this species and extensive empirical data are available for analysis. The number of small ponds on the central breeding grounds was used as an index to the state of the environment. Desirable properties of an optimal exploitation strategy were defined. A mathematical model was formulated to provide a synthesis of the existing literature, estimates of parameters developed from an analysis of data, and hypotheses regarding the specific effect of exploitation on total survival. Both the literature and the analysis of data were inconclusive concerning the effect of exploitation on survival. Therefore, alternative hypotheses were formulated: (1) exploitation mortality represents a largely additive form of mortality, or (2 ) exploitation mortality is compensatory with other forms of mortality, at least to some threshold level. Models incorporating these two hypotheses were formulated as stochastic dynamic programming models and optimal exploitation strategies were derived numerically on a digital computer. Optimal exploitation strategies were found to exist under rather general conditions. Direct feedback control was an integral component in the optimal decision-making process. Optimal exploitation was found to be substantially different depending upon the hypothesis regarding the effect of exploitation on the population. Assuming that exploitation is largely an additive force of mortality, optimal exploitation decisions are a convex function of the size of the breeding population and a linear or slightly concave function of the environmental conditions. Optimal exploitation under this hypothesis tends to reduce the variance of the size of the population. Under the hypothesis of compensatory mortality forces, optimal exploitation decisions are approximately linearly related to the size of the breeding population. Environmental variables may be somewhat more important than the size of the breeding population to the production of young mallards. In contrast, the size of the breeding population appears to be more important in the exploitation process than is the state of the environment. The form of the exploitation strategy appears to be relatively insensitive to small changes in the production rate. In general, the relative importance of the size of the breeding population may decrease as fecundity increases. The optimal level of exploitation in year t must be based on the observed size of the population and the state of the environment in year t unless the dynamics of the population, the state of the environment, and the result of the exploitation decisions are completely deterministic. Exploitation based on an average harvest, harvest rate, or designed to maintain a constant breeding population size is inefficient.

Dissertation Abstracts International

Wave propagation in layered soil deposits

Recent advances in the general theory of viscoelastic waves and rays in layered media provide a rigorous mathematical framework for site-specific, soil-response models used for earthquake resistant design. The advances provide general closed-form anelastic solutions for the classic problems of the response of a stack of soil layers to S and P waves, ray theory for reflected and refracted waves, Rayleigh- and Love-Type surface waves, and head waves. These general solutions valid for anelastic media regardless of the amount of material damping yield new insights regarding the characteristics of seismic waves and their ray paths that are not provided by conventional models. They provide corresponding numerical ground-response models and ray-tracing computation algorithms that account for changes in velocity and attenuation of anelastic waves associated with changes in inhomogeneity of the waves induced by anelastic soil and soil-rock boundaries. These anelastic effects manifest themselves as variations in amplitude response, amplitude attenuation, raypath location, and travel time as observed at the Earth’s surface. Implications of these anelastic effects for soil-response models used for earthquake resistant design are provided herein.

Conference Paper

Calibration of imperfect geophysical models by multiple satellite interferograms with measurement bias

Model calibration consists of using experimental or field data to estimate the unknown parameters of a mathematical model. The presence of model discrepancy and measurement bias in the data complicates this task. Satellite interferograms, for instance, are widely used for calibrating geophysical models in geological hazard quantification. In this work, we used satellite interferograms to relate ground deformation observations to the properties of the magma chamber at K i¯lauea Volcano in Hawai‘i. We derived closed-form marginal likelihoods and implemented posterior sampling procedures that simultaneously estimate the model discrepancy of physical models, and the measurement bias from the atmospheric error in satellite interferograms. We found that model calibration by aggregating multiple interferograms and downsampling the pixels in the interferograms can reduce the computation complexity compared to calibration approaches based on multiple data sets. The conditions that lead to no loss of information from data aggregation and downsampling are studied. Simulation illustrates that both discrepancy and measurement bias can be estimated, and real applications demonstrate that modeling both effects helps obtain a reliable estimation of a physical model’s unobserved parameters and enhance its predictive accuracy. We implement the computational tools in the RobustCalibration package available on CRAN.

Technometrics

Steady and intermittent slipping in a model of landslide motion regulated by pore-pressure feedback

This paper studies a parsimonious model of landslide motion, which consists of the one-dimensional diffusion equation (for pore pressure) coupled through a boundary condition to a first-order ODE (Newton's second law). Velocity weakening of sliding friction gives rise to nonlinearity in the model. Analysis shows that solutions of the model equations exhibit a subcritical Hopf bifurcation in which stable, steady sliding can transition to cyclical, stick-slip motion. Numerical computations confirm the analytical predictions of the parameter values at which bifurcation occurs. The existence of stick-slip behavior in part of the parameter space is particularly noteworthy because, unlike stick-slip behavior in classical models , here it arises in the absence of a reversible (elastic) driving force. Instead, the driving force is static (gravitational), mediated by the effects of pore-pressure diffusion on frictional resistance.

SIAM Journal on Applied Mathematics

Forecasting runout of rock and debris avalanches

Physically based mathematical models and statistically based empirical equations each may provide useful means of forecasting runout of rock and debris avalanches. This paper compares the foundations, strengths, and limitations of a physically based model and a statistically based forecasting method, both of which were developed to predict runout across three-dimensional topography. The chief advantage of the physically based model results from its ties to physical conservation laws and well-tested axioms of soil and rock mechanics, such as the Coulomb friction rule and effective-stress principle. The output of this model provides detailed information about the dynamics of avalanche runout, at the expense of high demands for accurate input data, numerical computation, and experimental testing. In comparison, the statistical method requires relatively modest computation and no input data except identification of prospective avalanche source areas and a range of postulated avalanche volumes. Like the physically based model, the statistical method yields maps of predicted runout, but it provides no information on runout dynamics. Although the two methods differ significantly in their structure and objectives, insights gained from one method can aid refinement of the other.

Conference Paper

A suggestion for computing objective function in model calibration

A parameter-optimization process (model calibration) is usually required for numerical model applications, which involves the use of an objective function to determine the model cost (model-data errors). The sum of square errors (SSR) has been widely adopted as the objective function in various optimization procedures. However, ‘square error’ calculation was found to be more sensitive to extreme or high values. Thus, we proposed that the sum of absolute errors (SAR) may be a better option than SSR for model calibration. To test this hypothesis, we used two case studies—a hydrological model calibration and a biogeochemical model calibration—to investigate the behavior of a group of potential objective functions: SSR, SAR, sum of squared relative deviation (SSRD), and sum of absolute relative deviation (SARD). Mathematical evaluation of model performance demonstrates that ‘absolute error’ (SAR and SARD) are superior to ‘square error’ (SSR and SSRD) in calculating objective function for model calibration, and SAR behaved the best (with the least error and highest efficiency). This study suggests that SSR might be overly used in real applications, and SAR may be a reasonable choice in common optimization implementations without emphasizing either high or low values (e.g., modeling for supporting resources management).

Ecological Informatics

Use of the truncated shifted Pareto distribution in assessing size distribution of oil and gas fields

The truncated shifted Pareto (TSP) distribution, a variant of the two-parameter Pareto distribution, in which one parameter is added to shift the distribution right and left and the right-hand side is truncated, is used to model size distributions of oil and gas fields for resource assessment. Assumptions about limits to the left-hand and right-hand side reduce the number of parameters to two. The TSP distribution has advantages over the more customary lognormal distribution because it has a simple analytic expression, allowing exact computation of several statistics of interest, has a "J-shape," and has more flexibility in the thickness of the right-hand tail. Oil field sizes from the Minnelusa play in the Powder River Basin, Wyoming and Montana, are used as a case study. Probability plotting procedures allow easy visualization of the fit and help the assessment. ?? 1988 International Association for Mathematical Geology.

Mathematical Geology

A transient laboratory method for determining the hydraulic properties of 'tight' rocks-I. Theory

Transient pulse testing has been employed increasingly in the laboratory to measure the hydraulic properties of rock samples with low permeability. Several investigators have proposed a mathematical model in terms of an initial-boundary value problem to describe fluid flow in a transient pulse test. However, the solution of this problem has not been available. In analyzing data from the transient pulse test, previous investigators have either employed analytical solutions that are derived with the use of additional, restrictive assumptions, or have resorted to numerical methods. In Part I of this paper, a general, analytical solution for the transient pulse test is presented. This solution is graphically illustrated by plots of dimensionless variables for several cases of interest. The solution is shown to contain, as limiting cases, the more restrictive analytical solutions that the previous investigators have derived. A method of computing both the permeability and specific storage of the test sample from experimental data will be presented in Part II.

International Journal of Rock Mechanics and Mining

Hydrograph simulation models of the Hillsborough and Alafia Rivers, Florida: a preliminary report

Mathematical (digital) models that simulate flood hydrographs from rainfall records have been developed for the following gaging stations in the Hillsborough and Alafia River basins of west-central Florida: Hillsborough River near Tampa, Alafia River at Lithia, and north Prong Alafia River near Keysville. These models, which were developed from historical streamflow and and rainfall records, are based on rainfall-runoff and unit-hydrograph procedures involving an arbitrary separation of the flood hydrograph. These models assume the flood hydrograph to be composed of only two flow components, direct (storm) runoff, and base flow. Expressions describing these two flow components are derived from streamflow and rainfall records and are combined analytically to form algorithms (models), which are programmed for processing on a digital computing system. Most Hillsborough and Alafia River flood discharges can be simulated with expected relative errors less than or equal to 30 percent and flood peaks can be simulated with average relative errors less than 15 percent. Because of the inadequate rainfall network that is used in obtaining input data for the North Prong Alafia River model, simulated peaks are frequently in error by more than 40 percent, particularly for storms having highly variable areal rainfall distribution. Simulation errors are the result of rainfall sample errors and, to a lesser extent, model inadequacy. Data errors associated with the determination of mean basin precipitation are the result of the small number and poor areal distribution of rainfall stations available for use in the study. Model inadequacy, however, is attributed to the basic underlying theory, particularly the rainfall-runoff relation. These models broaden and enhance existing water-management capabilities within these basins by allowing the establishment and implementation of programs providing for continued development in these areas. Specifically, the models serve not only as a basis for forecasting floods, but also for simulating hydrologic information needed in flood-plain mapping and delineating and evaluating alternative flood control and abatement plans.

Florida

Documentation of a deep percolation model for estimating ground-water recharge

A deep percolation model, which operates on a daily basis, was developed to estimate long-term average groundwater recharge from precipitation. It has been designed primarily to simulate recharge in large areas with variable weather, soils, and land uses, but it can also be used at any scale. The physical and mathematical concepts of the deep percolation model, its subroutines and data requirements, and input data sequence and formats are documented. The physical processes simulated are soil moisture accumulation, evaporation from bare soil, plant transpiration, surface water runoff, snow accumulation and melt, and accumulation and evaporation of intercepted precipitation. The minimum data sets for the operation of the model are daily values of precipitation and maximum and minimum air temperature, soil thickness and available water capacity, soil texture, and land use. Long-term average annual precipitation, actual daily stream discharge, monthly estimates of base flow, Soil Conservation Service surface runoff curve numbers, land surface altitude-slope-aspect, and temperature lapse rates are optional. The program is written in the FORTRAN 77 language with no enhancements and should run on most computer systems without modifications. Documentation has been prepared so that program modifications may be made for inclusions of additional physical processes or deletion of ones not considered important. (Author 's abstract)

Open-File Report

Hydraulic modeling development and application in water resources engineering

The use of modeling has become widespread in water resources engineering and science to study rivers, lakes, estuaries, and coastal regions. For example, computer models are commonly used to forecast anthropogenic effects on the environment, and to help provide advanced mitigation measures against catastrophic events such as natural and dam-break floods. Linking hydraulic models to vegetation and habitat models has expanded their use in multidisciplinary applications to the riparian corridor. Implementation of these models in software packages on personal desktop computers has made them accessible to the general engineering community, and their use has been popularized by the need of minimal training due to intuitive graphical user interface front ends. Models are, however, complex and nontrivial, to the extent that even common terminology is sometimes ambiguous and often applied incorrectly. In fact, many efforts are currently under way in order to standardize terminology and offer guidelines for good practice, but none has yet reached unanimous acceptance. This chapter provides a view of the elements involved in modeling surface flows for the application in environmental water resources engineering. It presents the concepts and steps necessary for rational model development and use by starting with the exploration of the ideas involved in defining a model. Tangible form of those ideas is provided by the development of a mathematical and corresponding numerical hydraulic model, which is given with a substantial amount of detail. The issues of model deployment in a practical and productive work environment are also addressed. The chapter ends by presenting a few model applications highlighting the need for good quality control in model validation.

Book chapter

Parallel Fortran-MPI software for numerical inversion of the Laplace transform and its application to oscillatory water levels in groundwater environments

A parallel Fortran-MPI (Message Passing Interface) software for numerical inversion of the Laplace transform based on a Fourier series method is developed to meet the need of solving intensive computational problems involving oscillatory water level's response to hydraulic tests in a groundwater environment. The software is a parallel version of ACM (The Association for Computing Machinery) Transactions on Mathematical Software (TOMS) Algorithm 796. Running 38 test examples indicated that implementation of MPI techniques with distributed memory architecture speedups the processing and improves the efficiency. Applications to oscillatory water levels in a well during aquifer tests are presented to illustrate how this package can be applied to solve complicated environmental problems involved in differential and integral equations. The package is free and is easy to use for people with little or no previous experience in using MPI but who wish to get off to a quick start in parallel computing. ?? 2004 Elsevier Ltd. All rights reserved.

Environmental Modelling and Software

Conversion and comparison of the mathematical, three-dimensional, finite-difference, ground-water flow model to the modular, three-dimensional, finite-difference, ground-water flow model for the Tesuque aquifer system in northern New Mexico

The three-dimensional finite-difference groundwater model (using a mathematical groundwater flow code) of the Tesuque aquifer system in northern New Mexico was converted to run using the U.S. Geological Survey 's modular groundwater flow code. Results from the final versions of the predevelopment and 1947 to 2080 transient simulations of the two models are compared. A correlation coefficient of 0.9905 was obtained for the match in block-by-block head-dependent fluxes for predevelopment conditions. There are, however, significant differences in at least two specific cases. In the first case, a difference is associated with the net loss from the Pojoaque River and its tributaries to the aquifer. The net loss by the river is given as 1.134 cu ft/sec using the original groundwater model, which is 38.1% less than the net loss by the river of 1.8319 cu ft/sec computed in this study. In the second case, the large difference is computed for the transient decline in the hydraulic head of a model block near Tesuque Pueblo. The hydraulic-head decline by 2080 is, using the original model, 249 ft, which is 14.7% less than the hydraulic head of 292 ft computed by this study. In general, the differences between the two sets of results are not large enough to lead to different conclusions regarding the behavior of the system at steady state or when pumped. (USGS)

Open-File Report

Ground-water and surface-water flow and estimated water budget for Lake Seminole, southwestern Georgia and northwestern Florida

Lake Seminole is a 37,600-acre impoundment formed at the confluence of the Flint and Chattahoochee Rivers along the Georgia?Florida State line. Outflow from Lake Seminole through Jim Woodruff Lock and Dam provides headwater to the Apalachicola River, which is a major supply of freshwater, nutrients, and detritus to ecosystems downstream. These rivers,together with their tributaries, are hydraulically connected to karst limestone units that constitute most of the Upper Floridan aquifer and to a chemically weathered residuum of undifferentiated overburden. The ground-water flow system near Lake Seminole consists of the Upper Floridan aquifer and undifferentiated overburden. The aquifer is confined below by low-permeability sediments of the Lisbon Formation and, generally, is semiconfined above by undifferentiated overburden. Ground-water flow within the Upper Floridan aquifer is unconfined or semiconfined and discharges at discrete points by springflow or diffuse leakage into streams and other surface-water bodies. The high degree of connectivity between the Upper Floridan aquifer and surface-water bodies is limited to the upper Eocene Ocala Limestone and younger units that are in contact with streams in the Lake Seminole area. The impoundment of Lake Seminole inundated natural stream channels and other low-lying areas near streams and raised the water-level altitude of the Upper Floridan aquifer near the lake to nearly that of the lake, about 77 feet. Surface-water inflow from the Chattahoochee and Flint Rivers and Spring Creek and outflow to the Apalachicola River through Jim Woodruff Lock and Dam dominate the water budget for Lake Seminole. About 81 percent of the total water-budget inflow consists of surface water; about 18 percent is ground water, and the remaining 1 percent is lake precipitation. Similarly, lake outflow consists of about 89 percent surface water, as flow to the Apalachicola River through Jim Woodruff Lock and Dam, about 4 percent ground water, and about 2 percent lake evaporation. Measurement error and uncertainty in flux calculations cause a flow imbalance of about 4 percent between inflow and outflow water-budget components. Most of this error can be attributed to errors in estimating ground-water discharge from the lake, which was calculated using a ground-water model calibrated to October 1986 conditions for the entire Apalachicola?Chattahoochee?Flint River Basin and not just the area around Lake Seminole. Evaporation rates were determined using the preferred, but mathematically complex, energy budget and five empirical equations: Priestley-Taylor, Penman, DeBruin-Keijman, Papadakis, and the Priestley-Taylor used by the Georgia Automated Environmental Monitoring Network. Empirical equations require a significant amount of data but are relatively easy to calculate and compare well to long-term average annual (April 2000?March 2001) pan evaporation, which is 65 inches. Calculated annual lake evaporation, for the study period, using the energy-budget method was 67.2 inches, which overestimated long-term average annual pan evaporation by 2.2 inches. The empirical equations did not compare well with the energy-budget method during the 18-month study period, with average differences in computed evaporation using each equation ranging from 8 to 26 percent. The empirical equations also compared poorly with long-term average annual pan evaporation, with average differences in evaporation ranging from 3 to 23 percent. Energy budget and long-term average annual pan evaporation estimates did compare well, with only a 3-percent difference between estimates. Monthly evaporation estimates using all methods ranged from 0.7 to 9.5 inches and were lowest during December 2000 and highest during May 2000. Although the energy budget is generally the preferred method, the dominance of surface water in the Lake Seminole water budget makes the method inaccurate and difficult to use, because surface water makes up m

Florida, Georgia

Are equilibrium shoreline models just convolutions?

Yes. Equilibrium shoreline models, which simulate wave-driven cross-shore erosion and accretion, are mathematically equivalent to a discrete convolution (i.e., a weighted, moving average) of a time series of wave-forcing conditions with a parameterized memory-decay kernel function. The direct equivalence between equilibrium shoreline models and convolutions reveals key theoretical aspects of equilibrium behavior. Convolutions (representing quasi-low-pass filter operations) provide an intuitive theoretical description of shoreline erosion and accretion behavior in response to waves: that is, shoreline position often mirrors the weighted moving average of wave time series. Model-convolution equivalence also provides a conceptual basis to interpret, evaluate, and construct data-driven Machine-Learning/Deep-Learning (ML/DL) models that use convolutions to extract features from data and then apply them for prediction (e.g., Convolutional Neural Networks (CNNs)). Finally, our findings provide a methodological pathway (based on Fourier transforms) for future understanding of wave-driven shoreline change, which can be used to interpret the coherence between the frequency spectrum of the processes of waves and shoreline change and construct more computationally efficient and effective shoreline-modeling approaches.

JGR Earth Surface

Multi-model approach to petroleum resource appraisal using analytic methodologies for probabilistic systems

The geologic appraisal model that is selected for a petroleum resource assessment depends upon purpose of the assessment, basic geologic assumptions of the area, type of available data, time available before deadlines, available human and financial resources, available computer facilities, and, most importantly, the available quantitative methodology with corresponding computer software and any new quantitative methodology that would have to be developed. Therefore, different resource assessment projects usually require different geologic models. Also, more than one geologic model might be needed in a single project for assessing different regions of the study or for cross-checking resource estimates of the area. Some geologic analyses used in the past for petroleum resource appraisal involved play analysis. The corresponding quantitative methodologies of these analyses usually consisted of Monte Carlo simulation techniques. A probabilistic system of petroleum resource appraisal for play analysis has been designed to meet the following requirements: (1) includes a variety of geologic models, (2) uses an analytic methodology instead of Monte Carlo simulation, (3) possesses the capacity to aggregate estimates from many areas that have been assessed by different geologic models, and (4) runs quickly on a microcomputer. Geologic models consist of four basic types: reservoir engineering, volumetric yield, field size, and direct assessment. Several case histories and present studies by the U.S. Geological Survey are discussed. ?? 1988 International Association for Mathematical Geology.

Mathematical Geology

Effective groundwater model calibration: With analysis of data, sensitivities, predictions, and uncertainty

Methods and guidelines for developing and using mathematical models Turn to Effective Groundwater Model Calibration for a set of methods and guidelines that can help produce more accurate and transparent mathematical models. The models can represent groundwater flow and transport and other natural and engineered systems. Use this book and its extensive exercises to learn methods to fully exploit the data on hand, maximize the model's potential, and troubleshoot any problems that arise. Use the methods to perform: Sensitivity analysis to evaluate the information content of data Data assessment to identify (a) existing measurements that dominate model development and predictions and (b) potential measurements likely to improve the reliability of predictions Calibration to develop models that are consistent with the data in an optimal manner Uncertainty evaluation to quantify and communicate errors in simulated results that are often used to make important societal decisions Most of the methods are based on linear and nonlinear regression theory. Fourteen guidelines show the reader how to use the methods advantageously in practical situations. Exercises focus on a groundwater flow system and management problem, enabling readers to apply all the methods presented in the text. The exercises can be completed using the material provided in the book, or as hands-on computer exercises using instructions and files available on the text's accompanying Web site. Throughout the book, the authors stress the need for valid statistical concepts and easily understood presentation methods required to achieve well-tested, transparent models. Most of the examples and all of the exercises focus on simulating groundwater systems; other examples come from surface-water hydrology and geophysics. The methods and guidelines in the text are broadly applicable and can be used by students, researchers, and engineers to simulate many kinds systems.

Book

Dynamics of seismogenic volcanic extrusion resisted by a solid surface plug, Mount St. Helens, 2004-2005

The 2004-5 eruption of Mount St. Helens exhibited sustained, near-equilibrium behavior characterized by nearly steady extrusion of a solid dacite plug and nearly periodic occurrence of shallow earthquakes. Diverse data support the hypothesis that these earthquakes resulted from stick-slip motion along the margins of the plug as it was forced incrementally upward by ascending, solidifying, gas-poor magma. I formalize this hypothesis with a mathematical model derived by assuming that magma enters the base of the eruption conduit at a steady rate, invoking conservation of mass and momentum of the magma and plug, and postulating simple constitutive equations that describe magma and conduit compressibilities and friction along the plug margins. Reduction of the model equations reveals a strong mathematical analogy between the dynamics of the magma-plug system and those of a variably damped oscillator. Oscillations in extrusion velocity result from the interaction of plug inertia, a variable upward force due to magma pressure, and a downward force due to the plug weight. Damping of oscillations depends mostly on plug-boundary friction, and oscillations grow unstably if friction exhibits rate weakening similar to that observed in experiments. When growth of oscillations causes the extrusion rate to reach zero, however, gravity causes friction to reverse direction, and this reversal instigates a transition from unstable oscillations to self-regulating stick-slip cycles. The transition occurs irrespective of the details of rate-weakening behavior, and repetitive stick-slip cycles are, therefore, robust features of the system’s dynamics. The presence of a highly compressible elastic driving element (that is, magma containing bubbles) appears crucial for enabling seismogenic slip events to occur repeatedly at the shallow earthquake focal depths (<1 km) observed during the 2004-5 eruption. Computations show that fluctuations in magma pressure accompanying such slip events are <3 kPa, indicating that deviations from mechanical equilibrium are slight and that coseismic force drops are <10 8 N. These results imply that the system’s self-regulating behavior is not susceptible to dramatic change--provided that the rate of magma ascent remains similar to the rate of magma accretion at the base of the plug, that plug surface erosion more or less compensates for mass gain due to basal accretion, and that magma and rock properties do not change significantly. Even if disequilibrium initial conditions are imposed, the dynamics of the magma-plug system are strongly attracted to self-regulating stick-slip cycles, although this self-regulating behavior can be bypassed on the way to runaway behavior if the initial state is too far from equilibrium.

Washington