SEARCH · Geology Reports

Results for “Nature Methods”

Search indexed USGS publications on groundwater, aquifers, geologic maps, mineral resources and earthquakes. Explore source records by subject and place.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,099 records · Page 61Linked to original sources

USGS2008

An exact solution for ideal dam-break floods on steep slopes

The shallow‐water equations are used to model the flow resulting from the sudden release of a finite volume of frictionless, incompressible fluid down a uniform slope of arbitrary inclination. The hodograph transformation and Riemann's method make it possible to transform the governing equations into a linear system and then deduce an exact analytical solution expressed in terms of readily evaluated integrals. Although the solution treats an idealized case never strictly realized in nature, it is uniquely well‐suited for testing the robustness and accuracy of numerical models used to model shallow‐water flows on steep slopes.

Water Resources Research↗

USGS1977

Additional results on palaeomagnetic stratigraphy of the Koobi Fora Formation, east of Lake Turkana (Lake Rudolf), Kenya

The magnetostratigraphy of the hominid-bearing sediments exposed east of Lake Turkana has been strengthened by new palaeomagnetic results. Ages obtained from several tuffs by the 40Ar/39Ar method suggest an approxmate match between the observed magnetozones and the geomagnetic polarity time scale; however, the palaeomagnetic results are also compatible with a younger chronology suggested by conventional K-Ar dating of the KBS Tuff. ?? 1977 Nature Publishing Group.

Nature↗

USGS2020

Postfire growth of seeded and planted big sagebrush - Strategic designs for restoring Greater Sage-grouse nesting habitat

Wildfires change plant community structure and impact wildlife habitat and population dynamics. Recent wildfire‐induced losses of big sagebrush ( Artemisia tridentata ) in North American shrublands are outpacing natural recovery and leading to substantial losses in habitat for sagebrush‐obligate species such as Greater Sage‐grouse. Managers are considering restoration strategies that include planting container‐grown sagebrush to improve establishment within areas using more conventional seeding methods. Although it is thought that planting sagebrush provides initial structural advantages over seeding, empirical comparisons of sagebrush growth are lacking between individuals established post‐fire using both methods. Using a Bayesian hierarchical approach, we evaluated sagebrush height and canopy area growth rates for plants established in 26 seeded and 20 planted locations within the Great Basin. We then related recovery rates to previously published nesting habitat requirements for sage‐grouse. Under average weather conditions, planted or seeded sagebrush will require 3 or 4 years, respectively, and a relatively high density (≥ 2 plants/m 2 ) to achieve the minimum recommended canopy cover for sage‐grouse (15 %). Sagebrush grown in warmer and drier conditions met this cover goal months earlier. Although planted sagebrush reached heights to meet sage‐grouse nesting requirements (30 cm) one year earlier than seeded plants, seeded individuals were ~19 cm taller with 410 cm 2 more canopy area than planted sagebrush after 8 years. However, big sagebrush establishment from seed is unreliable. Strategically planting small, high density patches of container‐grown sagebrush in historic sage‐grouse nesting habitat combined with lower density seedings in larger surrounding areas may accelerate sage‐grouse habitat restoration.

California, Idaho, Nevada, Oregon, Utah↗

USGS1973

A review of benthic faunal surveys in San Francisco Bay

During the past 60 years, considerable effort has been expended in studies of the relations of the biotic community and physicochemical characteristics of San Francisco Bay water. In very recent years these studies have emphasized the relations between the 'state of health' of bottom-living invertebrates (the benthos) and the levels of pollutants in the bay. Benthic organisms, generally sessile, are unable to escape deleterious environmental changes, and they reflect these changes in alterations of normal species composition of assemblages and species abundance. Data that expands understanding of these relations in urbanized areas such as San Francisco Bay are critical. Because of the implications of such data in control of water quality, the U.S. Geological Survey undertook a review of the results and major conclusions of San Francisco Bay benthic surveys. The size and species composition of faunal assemblages are largely controlled by the salinity of the water, the texture of the bottom sediments, and locally by wastes discharged into the bay. Efforts to describe the structure and function of benthic communities of the bay and to quantify the effects of waste discharge on them have been hampered by inconsistent and often faulty sampling methodology and species identification. Studies made show a lack of information on the normal life processes of the organisms concerned. The diversity index (a mathematical expression of the number of kinds of organisms present at a location), commonly used to describe the 'health' of the benthic community, has been employed without regard for the need for standardizing methodology and species identifications or for understanding natural biological processes that affect such mathematical indices. There are few reliable quantitative data on the distribution of benthic organisms in San Francisco Bay with which future assessments of the 'health' of the benthic community might be compared. Methods for study of the benthos must be standardized, identifications of species verified by trained taxonomists, and new field and laboratory studies undertaken before we can expect to obtain an accurate description of the benthic fauna and its relations with the environment.

Circular↗

USGS2010

Noninvasive methods for monitoring bear population trends

The U.S. Geological Survey began a grizzly bear research project in 2009 in the Northern Continental Divide Ecosystem (NCDE) of northwestern Montana. This work uses hair collection and DNA analysis methods similar to those used in the 2004 Northern Divide Grizzly Bear Project. However, instead of producing a snapshot of population size, the objectives of this new work are to estimate population growth rates by collecting hair at natural bear rubs along trails, roads, and fence and power lines. This approach holds promise of providing reliable estimates of population trends in an efficient, cost-effective, and unobtrusive way.

Fact Sheet↗

USGS2009

Abundance Trends and Status of the Little Colorado River Population of Humpback Chub: An Update Considering Data From 1989-2008

Mark-recapture methods have been used for the past two decades to assess trends in adult abundance and recruitment of the Little Colorado River (LCR) population of humpback chub. These methods indicate that the adult population declined through the 1980s and early 1990s but has been increasing for the past decade. Recruitment appears also to have increased, particularly in the 2003-4 period. Considering a range of assumed natural mortality-rates and magnitude of ageing error, it is unlikely that there are currently less than 6,000 adults or more than 10,000 adults. Our best estimate of the current adult (age 4 years or more) population is approximately 7,650 fish. Recent humpback chub assessments using the Age-Structured Mark-Recapture model (ASMR) and reported in 2006 (Melis and others, 2006) and 2008 (Coggins, 2008a,b) have provided abundance and recruitment trend estimates that have changed progressively over time as more data are considered by the model. The general pattern of change implies a less severe decline in adult abundance during the late 1980s through early 1990s, with attendant changes in recruitment supporting this demographic pattern. We have been concerned that these changes are not indicative of the true population and may be associated with a 'retrospective' bias as additional data are included in the ASMR model. To investigate this possibility, we developed a realistic individual-based simulation model (IBM) to generate replicate artificial data sets with similar characteristics to the true humpback chub data. The artificial data have known abundance trends and we analyzed these data with ASMR. On the basis of these simulations, we believe that errors in assigning age (and therefore brood-year) to fish based on their length are likely to have caused the retrospective bias pattern seen in the assessments and to have caused both less severe trends in the adult abundance estimates and progressively more severe downward bias in estimates of adult mortality-rates. This 'smearing', or assignment of fish from a single brood-year into multiple incorrect brood-years, is a result of variation in growth rates. The IBM simulations indicate that as a result of this error source, the best estimates of abundance and recruitment for any calendar year are those obtained from data collected previous to and within a year or two after each calendar year.

Open-File Report↗

USGS2014

Multi-elemental analysis of aqueous geological samples by inductively coupled plasma-optical emission spectrometry

Typically, 27 major, minor, and trace elements are determined in natural waters, acid mine drainage, extraction fluids, and leachates of geological and environmental samples by inductively coupled plasma-optical emission spectrometry (ICP-OES). At the discretion of the analyst, additional elements may be determined after suitable method modifications and performance data are established. Samples are preserved in 1–2 percent nitric acid (HNO3) at sample collection or as soon as possible after collection. The aqueous samples are aspirated into the ICP-OES discharge, where the elemental emission signals are measured simultaneously for 27 elements. Calibration is performed with a series of matrix-matched, multi-element solution standards.

Open-File Report↗

USGS1954

The investigation of chemical quality of water in tidal rivers

This report has been prepared for the guidance of personnel of the Water Resources Division who are engaged in water-quality investigations of tidal rivers. The study of tidal rivers is beset with many complexities not present in the investigation of non-tidal rivers. The periodic rise and fall of the tide may result in a corresponding periodic change in salinity at a sampling location on the tidal river. When the fresh water discharge is low, saline water may intrude up-river, and any factor changing the relative elevations of the ocean and the mean river level has an effect on the extent of salt-water intrusion. Variations in water composition between samples taken at several locations up or down river, at different depths, or at several locations across the stream are likely to be more pronounced than for similar sets of samples from a non-tidal stream. The nature of these variations and factors responsible for them are discussed, and the need for consideration of them in planning a sampling routine is stressed. The nature and mechanism of ocean-water intrusion in tidal rivers is discussed and sampling procedures for its detection are described. lllustrative examples - mostly from the work of the United States Geological Survey or State agencies - show various methods for correlating and presenting data from quality-of-water surveys of tidal rivers. Each tidal river presents an individual problem which can best be understood from a study of the factors involved. To that end the report is supplemented by an annotated bibliography of selected publications in the field.

Open-File Report↗

USGS1981

Comparison of the propane-area tracer method and predictive equations for determination of stream-reaeration coefficients on two small streams in Wisconsin

This study was made to identify the best predictive equations for a stream's reaeration-rate coefficient. Reaeration-rate information is needed in dissolvedoxygen modeling work, but an actual tracer measurement is not always possible. The propane-area gas-tracer method and predictive equations were compared for determination of stream-reaeration coefficients (K.2) for reaches of two small streams in Wisconsin. The study was made by the U.S. Geological Survey in cooperation with the Wisconsin Department of Natural Resources. The reaeration-rate coefficients actually measured by the propane-area tracer technique were 14.0 per day and 10.5 per day for two reaches of Honey Creek near Monroe, Wisconsin, with 6.98 per day and 0.98 per day measured at separate reaches on Mill Creek near Marshfield, Wisconsin. Of 20 predictive equations evaluated, the top five ranking equations were as follows: Tsivoglou-Neal with 34 percent mean error, Foree with 34.8 percent, Cadwallader with 45.5 percent, Isaacs-Gaudy with 45.8 percent, and Langbein- Durum with 49 percent.

Wisconsin↗

USGS1999

Estimating Concentrations of Road-Salt Constituents in Highway-Runoff from Measurements of Specific Conductance

Discrete or composite samples of highway runoff may not adequately represent in-storm water-quality fluctuations because continuous records of water stage, specific conductance, pH, and temperature of the runoff indicate that these properties fluctuate substantially during a storm. Continuous records of water-quality properties can be used to maximize the information obtained about the stormwater runoff system being studied and can provide the context needed to interpret analyses of water samples. Concentrations of the road-salt constituents calcium, sodium, and chloride in highway runoff were estimated from theoretical and empirical relations between specific conductance and the concentrations of these ions. These relations were examined using the analysis of 233 highwayrunoff samples collected from August 1988 through March 1995 at four highway-drainage monitoring stations along State Route 25 in southeastern Massachusetts. Theoretically, the specific conductance of a water sample is the sum of the individual conductances attributed to each ionic species in solution-the product of the concentrations of each ion in milliequivalents per liter (meq/L) multiplied by the equivalent ionic conductance at infinite dilution-thereby establishing the principle of superposition. Superposition provides an estimate of actual specific conductance that is within measurement error throughout the conductance range of many natural waters, with errors of less than ?5 percent below 1,000 microsiemens per centimeter (?S/cm) and ?10 percent between 1,000 and 4,000 ?S/cm if all major ionic constituents are accounted for. A semi-empirical method (adjusted superposition) was used to adjust for concentration effects-superposition-method prediction errors at high and low concentrations-and to relate measured specific conductance to that calculated using superposition. The adjusted superposition method, which was developed to interpret the State Route 25 highway-runoff records, accounts for contributions of constituents other than calcium, sodium, and chloride in dilute waters. The adjusted superposition method also accounts for the attenuation of each constituent's contribution to conductance as ionic strength increases. Use of the adjusted superposition method generally reduced predictive error to within measurement error throughout the range of specific conductance (from 37 to 51,500 ?S/cm) in the highway runoff samples. The effects of pH, temperature, and organic constituents on the relation between concentrations of dissolved constituents and measured specific conductance were examined but these properties did not substantially affect interpretation of the Route 25 data set. Predictive abilities of the adjusted superposition method were similar to results obtained by standard regression techniques, but the adjusted superposition method has several advantages. Adjusted superposition can be applied using available published data about the constituents in precipitation, highway runoff, and the deicing chemicals applied to a highway. This semi-empirical method can be used as a predictive and diagnostic tool before a substantial number of samples are collected, but the power of the regression method is based upon a large number of water-quality analyses that may be affected by a bias in the data.

Water-Resources Investigations Report↗

USGS2005

Rainfall-runoff in the Albuquerque, New Mexico, area: Measurements, analyses and comparisons

Albuquerque, New Mexico, has experienced significant growth over the last 20 years like many other cities in the Southwestern United States. While the US population grew by 37% between the 1970 and 2000 censuses, the growth for Albuquerque was 83%. More people mean more development and increased problems of managing runoff from urbanizing watersheds. The U.S. Geological Survey (USGS) in cooperation with the Albuquerque Arroyo Metropolitan Flood Control Authority (AMAFCA) and the City of Albuquerque has maintained a rainfall-runoff data collection program since 1976. The data from measured precipitation events can be used to verify hydrologic modeling. In this presentation, data from a representative gaged watershed is analyzed and discussed to set the overall framework for the rainfall-runoff process in the Albuquerque area. Of particular interest are the basic relationships between rainfall and watershed runoff response and an analysis of curve numbers as an indicator of runoff function. In urbanized areas, four land treatment types (natural, irrigated lawns, compacted soil, and impervious) are used to define surface infiltration conditions. Rainfall and runoff gage data are used to compare curve number (CN) and initial abstraction/uniform infiltration (IA/INF) techniques in an Albuquerque watershed. The IA/INF method appears to produce superior results over the CN method for the measured rainfall events.

Conference Paper↗

USGS2015

North Cascades National Park Service Complex

Washington↗

USGS2014

Mount Rainier National Park

Natural Resource Condition Assessments (NRCAs) evaluate current conditions for a subset of natural resources and resource indicators in national parks. NRCAs also report on trends in resource condition (when possible), identify critical data gaps, and characterize a general level of confidence for study findings. The resources and indicators emphasized in a given project depend on the park’s resource setting, status of resource stewardship planning and science in identifying high-priority indicators, and availability of data and expertise to assess current conditions for a variety of potential study resources and indicators. Although the primary objective of NRCAs is to report on current conditions relative to logical forms of reference conditions and values, NRCAs also report on trends, when appropriate (i.e., when the underlying data and methods support such reporting), as well as influences on resource conditions. These influences may include past activities or conditions that provide a helpful context for understanding current conditions and present-day threats and stressors that are best interpreted at park, watershed, or landscape scales (though NRCAs do not report on condition status for land areas and natural resources beyond park boundaries). Intensive cause-andeffect analyses of threats and stressors, and development of detailed treatment options, are outside the scope of NRCAs. It is also important to note that NRCAs do not address resources that lack sufficient data for assessment. For Mount Rainier National Park, this includes most invertebrate species and many other animal species that are subject to significant stressors from climate change and other anthropogenic sources such as air pollutants and recreational use. In addition, we did not include an analysis of the physical hydrology associated with streams (such as riverine landforms, erosion and aggradation which is significant in MORA streams), due to a loss of staff expertise from the USGS-BRD staff conducting the work, and human disturbance landcover issues such as the effects of roads, trails, and other anthropogenic developments due to lack of funds.

Report↗

USGS2015

An assessment of arthropod prey resources at Nakula Natural Area Reserve, a potential site of reintroduction for Kiwikiu (Pseudonestor xanthophrys) and Maui `Alauahio (Parareomyza montana)

Hawaiian forest birds have declined dramatically since humans arrived in the archipelago. Birds from all foraging guilds have been affected but insectivorous species are currently at greatest risk of extinction. On the island of Maui, populations and ranges of the insectivorous kiwikiu (Maui parrotbill; Pseudonestor xanthophrys) and Maui ‘alauahio (Maui creeper; Paroreomyza montana) have declined significantly from historic levels primarily due to habitat loss, predation,disease, and food web disruption, leading to federal listings of endangered species and species of concern, respectively. Recovery plans for these birds include reestablishment of populations in parts of their former range. Nakula Natural Area Reserve on the leeward side of HaleakalāVolcano has been targeted for release of wild-caught or captive-bred individuals. The mesic, montane koa-‘ōhi‘a (Acacia koa-Metrosideros polymorpha) forest at Nakula has been heavily impacted through grazing by feral ungulates, but recent management actions to exclude these animals are promoting forest recovery. The objective of this study was to assess the arthropod prey base at Nakula in preparation for reintroductions of kiwikiu and Maui ‘alauahio. To accomplish that goal, we compared arthropod abundances at Nakula to those at Hanawi Natural Area Reserve and Waikamoi Preserve, areas where kiwikiu and Maui ‘alauahio are currently found. We also identified diets of kiwikiu and Maui ‘alauahio from fecal samples to better understand and evaluate the prey base at Nakula. Assessment methods included clipping branch tips to sample arthropods within the foliage of koa and ‘ōhi‘a, using traps to quantify arthropods on koa and ‘ōhi‘a bark surfaces, counting exit holes to quantify abundances of beetles (Coleoptera) within dead branches of koa, and measuring the density of arthropods within the stems of ‘ākala (Rubus hawaiiensis). The diet of kiwikiu was dominated by caterpillars (Lepidoptera larvae), which comprised 90% of all prey items for 50 adult birds and 98% of all prey for two nestlings. Caterpillars were also the most important prey for Maui ‘alauahio (43% for 104 adult birds) although spiders (Araneae, 16%), beetles (12%) and true bugs, planthoppers and psyllids (Hemiptera; 12%) were also important. Caterpillars were generally the most abundant type of arthropod in the foliage of koa and ‘ōhi‘a, although spiders, beetles and hemipterans were also common. Total arthropod biomass and caterpillar biomass at Nakula was as great, or greater, than that observed at Hanawi and Waikamoi per unit of foliage of both koa and ‘ōhi‘a. Spiders generally dominated the bark fauna on both koa and ‘ōhi‘a at all sites although isopods (Isopoda), millipedes (Myriapoda: Millipeda) and lacewings (Neuroptera) were also abundant at Waikamoi and Hanawi. Total arthropod biomass on bark, as well as the biomass of several individual taxa, was significantly lower at Nakula than the other sites. Our measurement of the density of beetle exit holes in dead koa branches found no difference between Nakula and Waikamoi. Finally, no difference existed in the abundance of arthropods (primarily caterpillars and moth pupae) within ‘ākala stems among sites. With the exception of bark surfaces, our results suggest that the arthropod prey base for birds on primary foraging substrates at Nakula is similar to that found at two sites within the current range of kiwikiu and Maui ‘alauahio. However, our results should be viewed with caution because they are limited to the scale of individual branch, tree, or ‘ākala stem. To complete the assessment, our results should be scaled up to the landscape level by determining the density of each substrate within each site. Key arthropod prey of kiwikiu and Maui ‘alauahio are available at Nakula and, as habitat restoration continues, food abundance should increase to the point at which populations of these birds can be supported.

Hawaii↗

USGS2021

Food, culture and climate

The Social Sciences Coordinating Committee (SSCC) is one of multiple Interagency Groups that support the U.S. Global Change Research Program (USGCRP). USGCRP began as a Presidential initiative in 1989 and was mandated by Congress through the U.S. Global Change Research Act of 1990 “to assist the Nation and the world to understand, assess, predict, and respond to human-induced and natural processes of global change.” USGCRP is overseen by the Subcommittee on Global Change Research, composed of representatives from each of USGCRP’s 13 member agencies. The mission of the SSCC is to foster the integration of the methods, findings, and disciplinary perspectives of the social, behavioral, and economic sciences, along with interdisciplinary and transdisciplinary approaches that include these sciences, into USGCRP activities. The SSCC serves as a social science resource to other USGCRP interagency working groups, the Subcommittee on Global Change Research, and other USGCRP activities such as the National Climate Assessment. The SSCC’s Food, Culture, and Climate webinar series (held from September 14 to October 12, 2021) highlighted the ways in which social science research can elucidate the role of climate change in socio-cultural systems. It drew attention to the humanistic frameworks that underpin social scientific understanding of the ways individuals, households, and communities experience climate change. In particular, the seminars explored how the impacts of climate change are felt and understood by individuals and communities, how they interact with other stressors, and how they amplify existing inequities and vulnerabilities. This understanding is vital not only to the production of scientific knowledge, but also to the use of that knowledge in practice. This report provides a summary of the key takeaways from this webinar series. The webinar series organizers identified recurrent themes and salient points that emerged in conversations across the three events. Here, central issues of the discussions on the relationships among food, culture, and climate, as well as the role that social science plays in elucidating them, are synthesized and highlighted. Recordings of the webinar series are available online from USGCRP.

Webinar Series Report↗

USGS2013

Actual evapotranspiration modeling using the operational Simplified Surface Energy Balance (SSEBop) approach

Remote-sensing technology and surface-energy-balance methods can provide accurate and repeatable estimates of actual evapotranspiration ( ETa ) when used in combination with local weather datasets over irrigated lands. Estimates of ETa may be used to provide a consistent, accurate, and efficient approach for estimating regional water withdrawals for irrigation and associated consumptive use (CU), especially in arid cropland areas that require supplemental water due to insufficient natural supplies from rainfall, soil moisture, or groundwater. ETa in these areas is considered equivalent to CU, and represents the part of applied irrigation water that is evaporated and/or transpired, and is not available for immediate reuse. A recent U.S. Geological Survey study demonstrated the application of the remote-sensing-based Simplified Surface Energy Balance (SSEB) model to estimate 10-year average ETa at 1-kilometer resolution on national and regional scales, and compared those ETa values to the U.S. Geological Survey’s National Water-Use Information Program’s 1995 county estimates of CU. The operational version of the operational SSEB (SSEBop) method is now used to construct monthly, county-level ETa maps of the conterminous United States for the years 2000, 2005, and 2010. The performance of the SSEBop was evaluated using eddy covariance flux tower datasets compiled from 2005 datasets, and the results showed a strong linear relationship in different land cover types across diverse ecosystems in the conterminous United States (correlation coefficient [r] ranging from 0.75 to 0.95). For example, r for woody savannas (0.75), grassland (0.75), forest (0.82), cropland (0.84), shrub land (0.89), and urban (0.95). A comparison of the remote-sensing SSEBop method for estimating ETa and the Hamon temperature method for estimating potential ET ( ETp ) also was conducted, using regressions of all available county averages of ETa for 2005 and 2010, and yielded correlations of r = 0.60 and r = 0.71, respectively. Correlations generally are stronger in the Southeast where ETa is close to ETp . SSEBop ETa provides more spatial detail and accuracy in the Southwest where irrigation is practiced in a smaller proportion of the region.

Scientific Investigations Report↗

USGS2022

Incorporating streambank wells in stream mass loading studies to more effectively identify sources of solutes in stream water

Stream synoptic sampling studies that include flow estimates derived from the stream tracer dilution method are now commonly performed to identify sources and processes controlling solute transport to streams. However, a limitation of this mass-loading approach is its inability to identify the side of the stream on which a source is located in the common case where loading is largely from groundwater discharge. Such bank-specific loading information can be particularly valuable at mining-affected sites where both mining-related and natural background metal sources are often closely intermingled. Here we combine the stream mass-loading approach with data from pairs of shallow hand-installed streambank wells located on opposite sides of a gaining metal-impacted headwater stream to estimate bank-specific metal loading rates. Study results successfully identify a right-bank zone in the upper half of the primary study reach as the dominant source of loading, where groundwater discharge is elevated in metals due to natural weathering of sulfide-rich bedrock. A left-bank zone in the lower half of the primary study reach was also identified as a secondary, yet still substantial, loading source, where groundwater is likely impacted by portal discharge and/or waste piles associated with an adjacent abandoned mine. Determining the dominance of the left-bank mining-related source compared to other potential right-bank natural sulfide weathering sources would not have been possible without the streambank wells. Streambank wells also enabled collection of dissolved gas samples for groundwater dating. Computed piston-flow 3 H/ 3 He groundwater ages show little to no correlation with either metal concentrations or loading rates, suggesting that groundwater residence time and flow path variations exert little influence on groundwater discharge chemistry compared to variations in bedrock and soil composition. This study thus provides proof of concept that the proposed method of combining streambank well data with stream tracer injection and synoptic sampling can provide useful information on bank-specific mass loading rates and sources of contaminants affecting stream water quality.

Colorado↗

USGS2011

Modules based on the geochemical model PHREEQC for use in scripting and programming languages

The geochemical model PHREEQC is capable of simulating a wide range of equilibrium reactions between water and minerals, ion exchangers, surface complexes, solid solutions, and gases. It also has a general kinetic formulation that allows modeling of nonequilibrium mineral dissolution and precipitation, microbial reactions, decomposition of organic compounds, and other kinetic reactions. To facilitate use of these reaction capabilities in scripting languages and other models, PHREEQC has been implemented in modules that easily interface with other software. A Microsoft COM (component object model) has been implemented, which allows PHREEQC to be used by any software that can interface with a COM server—for example, Excel®, Visual Basic®, Python, or MATLAB". PHREEQC has been converted to a C++ class, which can be included in programs written in C++. The class also has been compiled in libraries for Linux and Windows that allow PHREEQC to be called from C++, C, and Fortran. A limited set of methods implements the full reaction capabilities of PHREEQC for each module. Input methods use strings or files to define reaction calculations in exactly the same formats used by PHREEQC. Output methods provide a table of user-selected model results, such as concentrations, activities, saturation indices, and densities. The PHREEQC module can add geochemical reaction capabilities to surface-water, groundwater, and watershed transport models. It is possible to store and manipulate solution compositions and reaction information for many cells within the module. In addition, the object-oriented nature of the PHREEQC modules simplifies implementation of parallel processing for reactive-transport models. The PHREEQC COM module may be used in scripting languages to fit parameters; to plot PHREEQC results for field, laboratory, or theoretical investigations; or to develop new models that include simple or complex geochemical calculations.

Computers & Geosciences↗