Search USGS⌕ Search

SEARCH · Search USGS

Results for “Journal of Water Process 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.

48 records · Page 3Linked to original sources
USGS2020

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↗
USGS2011

A distribution-based parameterization for improved tomographic imaging of solute plumes

Difference geophysical tomography (e.g. radar, resistivity and seismic) is used increasingly for imaging fluid flow and mass transport associated with natural and engineered hydrologic phenomena, including tracer experiments, in situ remediation and aquifer storage and recovery. Tomographic data are collected over time, inverted and differenced against a background image to produce ‘snapshots’ revealing changes to the system; these snapshots readily provide qualitative information on the location and morphology of plumes of injected tracer, remedial amendment or stored water. In principle, geometric moments (i.e. total mass, centres of mass, spread, etc.) calculated from difference tomograms can provide further quantitative insight into the rates of advection, dispersion and mass transfer; however, recent work has shown that moments calculated from tomograms are commonly biased, as they are strongly affected by the subjective choice of regularization criteria. Conventional approaches to regularization (Tikhonov) and parametrization (image pixels) result in tomograms which are subject to artefacts such as smearing or pixel estimates taking on the sign opposite to that expected for the plume under study. Here, we demonstrate a novel parametrization for imaging plumes associated with hydrologic phenomena. Capitalizing on the mathematical analogy between moment-based descriptors of plumes and the moment-based parameters of probability distributions, we design an inverse problem that (1) is overdetermined and computationally efficient because the image is described by only a few parameters, (2) produces tomograms consistent with expected plume behaviour (e.g. changes of one sign relative to the background image), (3) yields parameter estimates that are readily interpreted for plume morphology and offer direct insight into hydrologic processes and (4) requires comparatively few data to achieve reasonable model estimates. We demonstrate the approach in a series of numerical examples based on straight-ray difference-attenuation radar monitoring of the transport of an ionic tracer, and show that the methodology outlined here is particularly effective when limited data are available.

Geophysical Journal International↗
USGS2004

Hydraulic and geochemical framework of the Idaho National Engineering and Environmental Laboratory vadose zone

Questions of major importance for subsurface contaminant transport at the Idaho National Engineering and Environmental Laboratory (INEEL) include (i) travel times to the aquifer, both average or typical values and the range of values to be expected, and (ii) modes of contaminant transport, especially sorption processes. The hydraulic and geochemical framework within which these questions are addressed is dominated by extreme heterogeneity in a vadose zone and aquifer consisting of interbedded basalts and sediments. Hydraulically, major issues include diverse possible types of flow pathways, extreme anisotropy, preferential flow, combined vertical and horizontal flow, and temporary saturation or perching. Geochemically, major issues include contaminant mobility as influenced by redox conditions, the concentration of organic and inorganic complexing solutes and other local variables, the interaction with infiltrating waters and with the contaminant source environment, and the aqueous speciation of contaminants such as actinides. Another major issue is the possibility of colloid transport, which inverts some of the traditional concepts of mobility, as sorbed contaminants on mobile colloids may be transported with ease compared with contaminants that are not sorbed. With respect to the goal of minimizing aquifer concentrations of contaminants, some characteristics of the vadose zone are essentially completely favorable. Examples include the great thickness (200 m) of the vadose zone, and the presence of substantial quantities of fine sediments that can retard contaminant transport both hydraulically and chemically. Most characteristics, however, have both favorable and unfavorable aspects. For example, preferential flow, as promoted by several notable features of the vadose zone at the INEEL, can provide fast, minimally sorbing pathways for contaminants to reach the aquifer easily, but it also leads to a wide dispersal of contaminants in a large volume of subsurface material, thus increasing the opportunity for dilution and sorption.

Vadose Zone Journal↗
USGS2004

Object-based inversion of crosswell radar tomography data to monitor vegetable oil injection experiments

Crosswell radar methods can be used to dynamically image ground-water flow and mass transport associated with tracer tests, hydraulic tests, and natural physical processes, for improved characterization of preferential flow paths and complex aquifer heterogeneity. Unfortunately, because the raypath coverage of the interwell region is limited by the borehole geometry, the tomographic inverse problem is typically underdetermined, and tomograms may contain artifacts such as spurious blurring or streaking that confuse interpretation. We implement object-based inversion (using a constrained, non-linear, least-squares algorithm) to improve results from pixel-based inversion approaches that utilize regularization criteria, such as damping or smoothness. Our approach requires pre- and post-injection travel-time data. Parameterization of the image plane comprises a small number of objects rather than a large number of pixels, resulting in an overdetermined problem that reduces the need for prior information. The nature and geometry of the objects are based on hydrologic insight into aquifer characteristics, the nature of the experiment, and the planned use of the geophysical results. The object-based inversion is demonstrated using synthetic and crosswell radar field data acquired during vegetable-oil injection experiments at a site in Fridley, Minnesota. The region where oil has displaced ground water is discretized as a stack of rectangles of variable horizontal extents. The inversion provides the geometry of the affected region and an estimate of the radar slowness change for each rectangle. Applying petrophysical models to these results and porosity from neutron logs, we estimate the vegetable-oil emulsion saturation in various layers. Using synthetic- and field-data examples, object-based inversion is shown to be an effective strategy for inverting crosswell radar tomography data acquired to monitor the emplacement of vegetable-oil emulsions. A principal advantage of object-based inversion is that it yields images that hydrologists and engineers can easily interpret and use for model calibration.

Journal of Environmental & Engineering Geophysics↗
USGS2009

Hydrochemical analysis of groundwater using multivariate statistical methods - The Volta region, Ghana

Q and R-mode multivariate statistical analyses were applied to groundwater chemical data from boreholes and wells in the northern section of the Volta region Ghana. The objective was to determine the processes that affect the hydrochemistry and the variation of these processes in space among the three main geological terrains: the Buem formation, Voltaian System and the Togo series that underlie the area. The analyses revealed three zones in the groundwater flow system: recharge, intermediate and discharge regions. All three zones are clearly different with respect to all the major chemical parameters, with concentrations increasing from the perceived recharge areas through the intermediate regions to the discharge areas. R-mode HCA and factor analysis (using varimax rotation and Kaiser Criterion) were then applied to determine the significant sources of variation in the hydrochemistry. This study finds that groundwater hydrochemistry in the area is controlled by the weathering of silicate and carbonate minerals, as well as the chemistry of infiltrating precipitation. This study finds that the ??D and ?? 18 O data from the area fall along the Global Meteoric Water Line (GMWL). An equation of regression derived for the relationship between ??D and ?? 18 O bears very close semblance to the equation which describes the GMWL. On the basis of this, groundwater in the study area is probably meteoric and fresh. The apparently low salinities and sodicities of the groundwater seem to support this interpretation. The suitability of groundwater for domestic and irrigation purposes is related to its source, which determines its constitution. A plot of the sodium adsorption ratio (SAR) and salinity (EC) data on a semilog axis, suggests that groundwater serves good irrigation quality in the area. Sixty percent (60%), 20% and 20% of the 67 data points used in this study fall within the medium salinity - low sodicity (C2-S1), low salinity -low sodicity (C1-S1) and high salinity - low sodicity (C3-S1) fields, which ascribe good irrigation quality to groundwater from this area. Salinities range from 28.1 to 1956 ??S/cm, whilst SAR values fall within the range 0-3. Extremely low sodicity waters of this kind, with salinities lower than 600 ??S/cm, have the tendency to affect the dispersive properties of irrigation soils when used for irrigation. About 50% of the groundwater in the study area fall within this category and need prior treatment before usage. ?? 2009 Korean Society of Civil Engineers and Springer-Verlag Berlin Heidelberg GmbH.

KSCE Journal of Civil Engineering↗
USGS2003

The saturated zone at Yucca Mountain: An overview of the characterization and assessment of the saturated zone as a barrier to potential radionuclide migration

The US Department of Energy is pursuing Yucca Mountain, Nevada, for the development of a geologic repository for the disposal of spent nuclear fuel and high-level radioactive waste, if the repository is able to meet applicable radiation protection standards established by the US Nuclear Regulatory Commission and the US Environmental Protection Agency (EPA). Effective performance of such a repository would rely on a number of natural and engineered barriers to isolate radioactive waste from the accessible environment. Groundwater beneath Yucca Mountain is the primary medium through which most radionuclides might move away from the potential repository. The saturated zone (SZ) system is expected to act as a natural barrier to this possible movement of radionuclides both by delaying their transport and by reducing their concentration before they reach the accessible environment. Information obtained from Yucca Mountain Site Characterization Project activities is used to estimate groundwater flow rates through the site-scale SZ flow and transport model area and to constrain general conceptual models of groundwater flow in the site-scale area. The site-scale conceptual model is a synthesis of what is known about flow and transport processes at the scale required for total system performance assessment of the site. This knowledge builds on and is consistent with knowledge that has accumulated at the regional scale but is more detailed because more data are available at the site-scale level. The mathematical basis of the site-scale model and the associated numerical approaches are designed to assist in quantifying the uncertainty in the permeability of rocks in the geologic framework model and to represent accurately the flow and transport processes included in the site-scale conceptual model. Confidence in the results of the mathematical model was obtained by comparing calculated to observed hydraulic heads, estimated to measured permeabilities, and lateral flow rates calculated by the site-scale model to those calculated by the regional-scale flow model. In addition, it was confirmed that the flow paths leaving the region of the potential repository are consistent with those inferred from gradients of measured head and those independently inferred from water-chemistry data. The general approach of the site-scale SZ flow and transport model analysis is to calculate unit breakthrough curves for radionuclides at the interface between the SZ and the biosphere using the three-dimensional site-scale SZ flow and transport model. Uncertainties are explicitly incorporated into the site-scale SZ flow and transport abstractions through key parameters and conceptual models. ?? 2002 Elsevier Science B.V. All rights reserved.

Journal of Contaminant Hydrology↗
USGS2015

Ground-truthing electrical resistivity methods in support of submarine groundwater discharge studies: Examples from Hawaii, Washington, and California

Submarine groundwater discharge (SGD) is an important conduit that links terrestrial and marine environments. SGD conveys both water and water-borne constituents into coastal waters, where these inflows may impact near-shore ecosystem health and sustainability. Multichannel electrical resistivity techniques have proven to be a powerful tool to examine scales and dynamics of SGD and SGD forcings. However, there are uncertainties both in data aquisition and data processing that must be addressed to maximize the effectiveness of this tool in estuarine or marine environments. These issues most often relate to discerning subtle nuances in the flow of electricity through variably saturated media that can also be highly conductive ( i.e., seawater). Three contrasting field sites were examined for this study to assess the effectiveness of electrical resistivity techniques in varying coastal settings by comparing resistivity data to direct salinity and resistivity observations, quantifying changes in lithology and beach geomorphology, and fine-tuning inversion protocols. The three study sites all have substantial (up to 85 cm day −1 ) submarine groundwater discharge rates, but the hydrologic, oceanographic, and geologic characteristics of the sites are all very different. At a site in Pelekane Bay on the Big Island of Hawaii, seasonal flooding introduces very high concentrations of fine to coarse sediment into the bay. Near-shore circulation is limited in Pelekane Bay, so this newly introduced sediment can become deposited in the bay where it accumulates over time. At a site in Hood Canal, a fjord within Puget Sound, Washington, SGD rates can be high because of the large tidal range, abundant recharge, and steep hydrologic gradients. At Younger Lagoon in northern California, the flow of groundwater towards the coast is much more parsimonious, but here marine processes, including recirculated seawater, are important in controlling the flow of material towards the coast. Rigorous ground-truthing at each field site showed that multi-channel electrcial resistivity techniques can reproduce the scales and dynamics of a seepage field when such data are correctly collected, and when the model inversions are tuned to field site characteristics. Such information can provide a unique perspective on the scales and dynamics of exchange processes within a coastal aquifer—information essential to scientists and resource managers alike.

California, Hawaii, Washington↗
USGS2026

Response of juvenile Chinook Salmon to an experimental reservoir operation to facilitate downstream migration

Objectives High-head dams and storage reservoirs impact ecological processes, impeding recovery efforts of Pacific salmon Oncorhynchus spp. in the western United States. These projects can delay migration, increase risks of predation and disease, and induce mortality for juvenile Pacific salmon during downstream passage. To mitigate these negative effects, large-scale changes to dam operations have been proposed as interim measures or alternatives to engineered fish passage solutions. A novel deep drawdown operation of the Lookout Point Project (two dams and reservoirs) in Oregon serves as a case study for learning from a large-scale management experiment. Beginning in July 2024, Lookout Point Reservoir was slowly lowered to attain minimal levels from mid-November to mid-December with the intent of expediting downstream passage through the project for juvenile Chinook Salmon Oncorhynchus tshawytscha . We evaluated the impacts of the months-long operation on two life histories of juvenile Chinook Salmon with different rearing environments and growth rates. Methods Acoustic telemetry was used to monitor juvenile Chinook Salmon in the Lookout Point Project. Environmental data, including water temperature and hydrological conditions, were collected. We fitted a time-integrated migration survival model to these data to estimate the effects of management-influenced covariates and fish length on apparent survival and residence time in each reach of the Lookout Point Project. Results Larger fish had higher survival than smaller fish in each reach of the Lookout Point Project. Apparent survival through each dam of the Lookout Point Project was highest and forebay residence times lowest during the latter part of the drawdown (late November). However, early periods of the drawdown (September) were associated with a transition from a thermally stratified reservoir to a warm isothermal reservoir that negatively impacted reservoir survival. Conclusions The deep drawdown of the Lookout Point Project achieved the management intent to improve downstream passage for juvenile Chinook Salmon at the lowest reservoir elevations but also caused unintended changes in the reservoir rearing environment prior to full drawdown. Further research could resolve remaining uncertainties regarding the overall effectiveness of drawdown operations with respect to the population dynamics of different juvenile life history strategies. Understanding these factors would help optimize future management actions aimed at Pacific salmon recovery above the Lookout Point Project.

Oregon↗
USGS2014

High-resolution delineation of chlorinated volatile organic compounds in a dipping, fractured mudstone: depth- and strata-dependent spatial variability from rock-core sampling

Synthesis of rock-core sampling and chlorinated volatile organic compound (CVOC) analysis at five coreholes, with hydraulic and water-quality monitoring and a detailed hydrogeologic framework, was used to characterize the fine-scale distribution of CVOCs in dipping, fractured mudstones of the Lockatong Formation of Triassic age, of the Newark Basin in West Trenton, New Jersey. From these results, a refined conceptual model for more than 55 years of migration of CVOCs and depth- and strata-dependent rock-matrix contamination was developed. Industrial use of trichloroethene (TCE) at the former Naval Air Warfare Center (NAWC) from 1953 to 1995 resulted in dense non-aqueous phase liquid (DNAPL) TCE and dissolved TCE and related breakdown products, including other CVOCs, in underlying mudstones. Shallow highly weathered and fractured strata overlie unweathered, gently dipping, fractured strata that become progressively less fractured with depth. The unweathered lithology includes black highly fractured (fissile) carbon-rich strata, gray mildly fractured thinly layered (laminated) strata, and light-gray weakly fractured massive strata. CVOC concentrations in water samples pumped from the shallow weathered and highly fractured strata remain elevated near residual DNAPL TCE, but dilution by uncontaminated recharge, and other natural and engineered attenuation processes, have substantially reduced concentrations along flow paths removed from sources and residual DNAPL. CVOCs also were detected in most rock-core samples in source areas in shallow wells. In many locations, lower aqueous concentrations, compared to rock core concentrations, suggest that CVOCs are presently back-diffusing from the rock matrix. Below the weathered and highly fractured strata, and to depths of at least 50 meters (m), groundwater flow and contaminant transport is primarily in bedding-plane-oriented fractures in thin fissile high-carbon strata, and in fractured, laminated strata of the gently dipping mudstones. Despite more than 18 years of pump and treat (P&T) remediation, and natural attenuation processes, CVOC concentrations in aqueous samples pumped from these deeper strata remain elevated in isolated intervals. DNAPL was detected in one borehole during coring at a depth of 27 m. In contrast to core samples from the weathered zone, concentrations in core samples from deeper unweathered and unfractured strata are typically below detection. However, high CVOC concentrations were found in isolated samples from fissile black carbon-rich strata and fractured gray laminated strata. Aqueous-phase concentrations were correspondingly high in samples pumped from these strata via short-interval wells or packer-isolated zones in long boreholes. A refined conceptual site model considers that prior to P&T remediation groundwater flow was primarily subhorizontal in the higher-permeability near surface strata, and the bulk of contaminant mass was shallow. CVOCs diffused into these fractured and weathered mudstones. DNAPL and high concentrations of CVOCs migrated slowly down in deeper unweathered strata, primarily along isolated dipping bedding-plane fractures. After P&T began in 1995, using wells open to both shallow and deep strata, downward transport of dissolved CVOCs accelerated. Diffusion of TCE and other CVOCs from deeper fractures penetrated only a few centimeters into the unweathered rock matrix, likely due to sorption of CVOCs on rock organic carbon. Remediation in the deep, unweathered strata may benefit from the relatively limited migration of CVOCs into the rock matrix. Synthesis of rock core sampling from closely spaced boreholes with geophysical logging and hydraulic testing improves understanding of the controls on CVOC delineation and informs remediation design and monitoring.

New Jersey, New York, Pennsylvania↗
USGS2022

Upper Rio Grande Basin water-resource status and trends: Focus area study review and synthesis

The Upper Rio Grande Basin (URGB) is a critical international water resource under pressure from a myriad of climatic, ecological, infrastructural, water-use, and legal constraints. The objective of this study is to provide a comprehensive assessment of the spatial distribution and temporal trends of selected water-budget components (snow processes, evapotranspiration (ET), streamflow processes, and groundwater storage) using integrated analyses, such as watershed modeling and water availability and use data in the URGB over the past three decades. A spatially distributed snow evolution modeling system simulated snowpack processes over 34 years (1984–2017). It highlighted snow water equivalent declines from -35 to -77 mm/decade with widespread variability across elevation zones and land cover types. Gridded actual ET data from the SSEBop model were developed and tested for the URGB and demonstrated that all land-cover types had significant decreasing trends (1986-2015) ranging from -14 to -80 mm/decade. Conductivity-mass-balance (CMB) hydrograph separation results found that baseflow forms a large component of total streamflow, ranging from 29 to 69% (49% average) of total streamflow at 17 URGB sites upstream of Albuquerque, NM. Three of 4 graphical hydrograph separation methods in the U.S. Geological Survey Groundwater Toolbox were found to be inappropriate for estimating baseflow in the URGB; the most promising method, baseflow index (BFI) Standard, was optimized using CMB data and tested at three URGB sites, with resulting overestimation of 0 to 47%. Simulated changes in groundwater storage were extracted from historical and recent groundwater-flow models of select alluvial basins (San Luis, Española, Middle Rio Grande, and Tularosa-Hueco). In general, decreases in groundwater storage were observed from 1903 to 2013 except for the San Luis alluvial basin (Colorado), where periods of recovery are observed. The PRMS hydrologic model was successfully calibrated for 9 near-native subbasins (Nash-Sutcliffe efficiency 0.47 to 0.85) and parameters translated to the remaining subbasins; compared to simulated near-native flows (with minimal influence of reservoirs or diversions), observed Rio Grande streamgage flows demonstrated reductions of 40% or more for New Mexico and Texas areas of the basin. Significant decreasing trends (1980-2015) in precipitation, snowmelt rate, streamflow, and baseflow were observed at many of the 12 streamgage basins studied, which suggests that the decreasing trends for actual ET may be related to overall decreasing water availability in the basin, with negative implications for agricultural production and groundwater abstraction. Water security concerns arise from our findings of higher fraction precipitation as rain, slower snowmelt rates leading to decreasing streamflow production, and an increasing fraction of baseflow, all of which will affect the timing and magnitude of water available for human needs in the basin.

Upper Rio Grande Basin↗
USGS2018

A 30-m landsat-derived cropland extent product of Australia and China using random forest machine learning algorithm on Google Earth Engine cloud computing platform

Mapping high resolution (30-m or better) cropland extent over very large areas such as continents or large countries or regions accurately, precisely, repeatedly, and rapidly is of great importance for addressing the global food and water security challenges. Such cropland extent products capture individual farm fields, small or large, and are crucial for developing accurate higher-level cropland products such as cropping intensities, crop types, crop watering methods (irrigated or rainfed), crop productivity, and crop water productivity. It also brings many challenges that include handling massively large data volumes, computing power, and collecting resource intensive reference training and validation data over complex geographic and political boundaries. Thereby, this study developed a precise and accurate Landsat 30-m derived cropland extent product for two very important, distinct, diverse, and large countries: Australia and China. The study used of eight bands (blue, green, red, NIR, SWIR1, SWIR2, TIR1, and NDVI) of Landsat-8 every 16-day Operational Land Imager (OLI) data for the years 2013–2015. The classification was performed by using a pixel-based supervised random forest (RF) machine learning algorithm (MLA) executed on the Google Earth Engine (GEE) cloud computing platform. Each band was time-composited over 4–6 time-periods over a year using median value for various agro-ecological zones (AEZs) of Australia and China. This resulted in a 32–48-layer mega-file data-cube (MFDC) for each of the AEZs. Reference training and validation data were gathered from: (a) field visits, (b) sub-meter to 5-m very high spatial resolution imagery (VHRI) data, and (c) ancillary sources such as from the National agriculture bureaus. Croplands versus non-croplands knowledge base for training the RF algorithm were derived from MFDC using 958 reference-training samples for Australia and 2130 reference-training samples for China. The resulting 30-m cropland extent product was assessed for accuracies using independent validation samples: 900 for Australia and 1972 for China. The 30-m cropland extent product of Australia showed an overall accuracy of 97.6% with a producer’s accuracy of 98.8% (errors of omissions = 1.2%), and user’s accuracy of 79% (errors of commissions = 21%) for the cropland class. For China, overall accuracies were 94% with a producer’s accuracy of 80% (errors of omissions = 20%), and user’s accuracy of 84.2% (errors of commissions = 15.8%) for cropland class. Total cropland areas of Australia were estimated as 35.1 million hectares and 165.2 million hectares for China. These estimates were higher by 8.6% for Australia and 3.9% for China when compared with the traditionally derived national statistics. The cropland extent product further demonstrated the ability to estimate sub-national cropland areas accurately by providing an R 2 value of 0.85 when compared with province-wise cropland areas of China. The study provides a paradigm-shift on how cropland maps are produced using multi-date remote sensing. These products can be browsed at www.croplands.org and made available for download at NASA’s Land Processes Distributed Active Archive Center (LP DAAC) https://www.lpdaac.usgs.gov/node/1282 .

ISPRS Journal of Photogrammetry and Remote Sensing↗
USGS2010

Mineral Commodity Summaries 2010

Each chapter of the 2010 edition of the U.S. Geological Survey (USGS) Mineral Commodity Summaries (MCS) includes information on events, trends, and issues for each mineral commodity as well as discussions and tabular presentations on domestic industry structure, Government programs, tariffs, 5-year salient statistics, and world production and resources. The MCS is the earliest comprehensive source of 2009 mineral production data for the world. More than 90 individual minerals and materials are covered by two-page synopses. For mineral commodities for which there is a Government stockpile, detailed information concerning the stockpile status is included in the two-page synopsis. National reserves information for most mineral commodities found in this report, including those for the United States, are derived from a variety of sources. The ideal source of such information would be comprehensive evaluations that apply the same criteria to deposits in different geographic areas and report the results by country. In the absence of such evaluations, national reserves estimates compiled by countries for selected mineral commodities are a primary source of national reserves information. Lacking national assessment information by governments, sources such as academic articles, company reports, presentations by company representatives, and trade journal articles, or a combination of these, serve as the basis for national reserves information reported in the mineral commodity sections of this publication. A national estimate may be assembled from the following: historically reported reserves information carried for years without alteration because no new information is available; historically reported reserves reduced by the amount of historical production; and company reported reserves. International minerals availability studies conducted by the U.S. Bureau of Mines (USBM), before 1996, and estimates of identified resources by an international collaborative effort (the International Strategic Minerals Inventory) are the basis for some reserves estimates. The USGS collects information about the quantity and quality of mineral resources but does not directly measure reserves, and companies or governments do not directly report reserves to the USGS. Reassessment of reserves is a continuing process, and the intensity of this process differs for mineral commodities, countries, and time period. Throughout the history of Mineral Commodity Summaries and its predecessor prior to 1978, Commodity Data Summaries, the presentation of resource data has evolved. Although world resources have been discussed each year, presentation of reserves and reserve base data varied. From 1957 through 1979, only reserves information was published in the reports, but from 1980 through 1987, only estimates of reserve base, a concept introduced by the U.S. Bureau of Mines (USBM) and the USGS in 1980, were published. Beginning in 1988, both reserves and reserve base information were listed for each mineral commodity where applicable and available. Prior to 1996, the minerals availability studies conducted by the USBM and work with international collaborators were the basis for reserve base data reported in Mineral Commodity Summaries. When the USBM was closed in 1996, this function was discontinued. Since that time, reserve base estimates have been updated to be consistent with changes in reserves, but the nonreserves component of the information upon which the reserve base data were estimated is not current enough to support defensible reserve base estimates. For that reason, publication of reserve base estimates was discontinued for Mineral Commodity Summaries 2010. Abbreviations and units of measure, and definitions of selected terms used in the report, are in Appendix A and Appendix B, respectively. A resource/reserve classification for minerals, based on USGS Circular 831 (published with the U.S. Bureau of Mines) is Appendix C, and a directory of USGS minerals information country specialists and their responsibilities is Appendix D. The USGS continually strives to improve the value of its publications to users. Constructive comments and suggestions by readers of the MCS 2010 are welcomed.

Mineral Commodity Summaries↗