Geology Reports⌕ Search

SEARCH · Geology Reports

Results for “Analytical Methods”

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

The effects of flow-path modification on water-quality constituent retention in an urban stormwater detention pond and wetland system, Orlando, Florida

Changes in constituent retention in a wet stormwater-detention pond and wetland system in Orlando, Florida, were evaluated following the 1988 installation of a flow barrier which approximately doubled the flow path and increased detention time in the pond. The pond and wetland were arranged in series so that stormwater first enters the pond and overflows into the wetland before spilling over to the regional stream system. Several principal factors that contribute to constituent retention were examined, including changes in pond-water quality between storms, stormwater quality, and pond-water flushing during storms. A simple, analytical pond-water mixing model was used as the basis for interpreting changes in retention efficiencies caused by pond modification. Retention efficiencies were calculated by a modified event-mean concentration efficiency method using a minimum variance unbiased estimator approach. The results of this study generally support the hypothesis that changes in the geometry of stormwater treatment systems can significantly affect the constituent retention efficiency of the pond and wetland system. However, the results also indicate that these changes in efficiency are caused not only by changes in residence time, but also by changes in stormwater mixing and pond water flushing during storms. Additionally, the use of average efficiencies as indications of treatment effectiveness may fail to account for biases associated with sample distribution and independent physical properties of the system, such as the range and concentrations of constituents in stormwater inflows and stormwater volume. Changes in retention efficiencies varied among chemical constituents and were significantly different in the pond and wetland. Retention efficiency was related to inflow concentration for most constituents. Increased flushing of the pond after modification caused decreases in retention efficiencies for constituents that concentrate in the pond between storms (dissolved solids) and increases in retention efficiency for constituents that settle out of pond and wetland storage between storms. The greatest increase in retention efficiencies in the detention pond was observed for total lead, which increased from 19 percent before modification to 73 percent after modification. However, retention efficiencies for nutrients for nutrients and suspended constituents decreased in the wetland after modification. This was probably because of the flushing of accumulated sediments as a result of a change in flow path through the wetland. As a result, the overall effect of modification on the system (pond and wetland retention efficiencies combined) was a reduction in retention efficiency for all but two constituents (total zinc and total ammonia nitrogen).

Water-Resources Investigations Report↗

A precise method for the analysis of δ18O of dissolved inorganic phosphate in seawater

A method for preparation and analysis of the oxygen isotope composition (δ 18 O) of dissolved inorganic phosphate (DIP) has been developed and preliminary results for water samples from various locations are reported. Phosphate is extracted from seawater samples by coprecipitation with magnesium hydroxide. Phosphate is further purified through a series of precipitations and resin separation and is ultimately converted to silver phosphate. Silver phosphate samples are pyrolitically decomposed to carbon monoxide and analyzed for δ 18 O. Silver phosphate samples weighing 0.7 mg (3.5 µmol oxygen) can be analyzed routinely with an average standard deviation of about 0.3‰. There is no isotope fractionation during extraction and blanks are negligible within analytical error. Reproducibility was determined for both laboratory standards and natural samples by multiple analyses. A comparison between filtered and unfiltered natural seawater samples was also conducted and no appreciable difference was observed for the samples tested. The δ 18 O values of DIP in seawater determined using this method range from 18.6‰ to 22.3‰, suggesting small but detectable natural variability in seawater. For the San Francisco Bay estuary DIP δ 18 O is more variable, ranging from 11.4‰ near the San Joaquin River to 20.1‰ near the Golden Gate Bridge, and was well correlated with salinity, phosphate concentration, and δ 18 O of water.

California↗

Preliminary investigations on computational methods for solving the two-point seismic ray-tracing problem in a heterogeneous and isotropic medium

Two-point seismic ray tracing is an important problem in seismology. In this report, we study this problem by first deriving the differential equations governing seismic wave propagation between two end points in a 3-dimensional heterogeneous and isotropic medium. We then formulate three numerical schemes to solve the two-point ray-tracing problem. In the first method, introduced originally by R. L. Wesson, the second order ray equations are reduced by the central difference approximation to a set of simultaneous non-linear equations, which are then solved by Newton's method. In the second method, the second order ray equations are reduced to a set of first order equations which are solved by using an adaptive finite difference method introduced by M. Lentini and V. Pereyra. The third method, due originally to L. Euler uses a sum to approximate the integral for the travel time which are solved directly for the minimum time path. This formulation may be shown to be equivalent to Wesson's method. We programed Wesson's method for a 3-dimensional inhomogeneous velocity medium, and also applied Lentini and Pereyra's method for a 2-dimensional linear velocity model. Our programs are tested for constant velocity, and linear velocity in one, and two, space variable models. Both programs give satisfactory answers in comparison with known analytic solutions, and with respect to each other.

Open-File Report↗

The oilspill risk analysis model of the U. S. Geological Survey

The U.S. Geological Survey has developed an oilspill risk analysis model to aid in estimating the environmental hazards of developing oil resources in Outer Continental Shelf (OCS) lease areas. The large, computerized model analyzes the probability of spill occurrence, as well as the likely paths or trajectories of spills in relation to the locations of recreational and biological resources which may be vulnerable. The analytical methodology can easily incorporate estimates of weathering rates , slick dispersion, and possible mitigating effects of cleanup. The probability of spill occurrence is estimated from information on the anticipated level of oil production and method and route of transport. Spill movement is modeled in Monte Carlo fashion with a sample of 500 spills per season, each transported by monthly surface current vectors and wind velocities sampled from 3-hour wind transition matrices. Transition matrices are based on historic wind records grouped in 41 wind velocity classes, and are constructed seasonally for up to six wind stations. Locations and monthly vulnerabilities of up to 31 categories of environmental resources are digitized within an 800,000 square kilometer study area. Model output includes tables of conditional impact probabilities (that is, the probability of hitting a target, given that a spill has occurred), as well as probability distributions for oilspills occurring and contacting environmental resources within preselected vulnerability time horizons. (USGS)

Open-File Report↗

The oilspill risk analysis model of the U. S. Geological Survey

The U.S. Geological Survey has developed an oilspill risk analysis model to aid in estimating the environmental hazards of developing oil resources in Outer Continental Shelf (OCS) lease areas. The large, computerized model analyzes the probability of spill occurrence, as well as the likely paths or trajectories of spills in relation to the locations of recreational and biological resources which may be vulnerable. The analytical methodology can easily incorporate estimates of weathering rates , slick dispersion, and possible mitigating effects of cleanup. The probability of spill occurrence is estimated from information on the anticipated level of oil production and method of route of transport. Spill movement is modeled in Monte Carlo fashion with a sample of 500 spills per season, each transported by monthly surface current vectors and wind velocities sampled from 3-hour wind transition matrices. Transition matrices are based on historic wind records grouped in 41 wind velocity classes, and are constructed seasonally for up to six wind stations. Locations and monthly vulnerabilities of up to 31 categories of environmental resources are digitized within an 800,000 square kilometer study area. Model output includes tables of conditional impact probabilities (that is, the probability of hitting a target, given that a spill has occurred), as well as probability distributions for oilspills occurring and contacting environmental resources within preselected vulnerability time horizons. (USGS)

Professional Paper↗

Keeping it classy: Classification of live fish and ghost PIT tags detected with a mobile PIT tag interrogation system using an innovative analytical approach

The ability of passive integrated transponder (PIT) tag data to improve demographic parameter estimates has led to the rapid advancement of PIT tag systems. However, ghost tags create uncertainty about detected tag status (i.e., live fish or ghost tag) when using mobile interrogation systems. We developed a method to differentiate between live fish and ghost tags using a random forest classification model with a novel data input structure based on known fate PIT tag detections in the San Juan River (New Mexico, Colorado, and Utah, USA). We used our model to classify detected tags with an overall error rate of 6.8% (1.6% ghost tags error rate and 21.8% live fish error rate). The important variables for classification were related to distance moved and response to monsoonal flood flows; however, habitat variables did not appear to influence model accuracy. Our results and approach allow the use of mobile detection data with confidence and allow for greater accuracy in movement, distribution, and habitat use studies, potentially helping identify influential management actions that would improve our ability to conserve and recover endangered fish.

Canadian Journal of Fisheries and Aquatic Sciences↗

Automatic measurements and computations for radiochemical analyses

In natural radioactive sources the most important radioactive daughter products useful for geochemical studies are protactinium-231, the alpha-emitting thorium isotopes, and the radium isotopes. To resolve the abundances of these thorium and radium isotopes by their characteristic decay and growth patterns, a large number of repeated alpha activity measurements on the two chemically separated elements were made over extended periods of time. Alpha scintillation counting with automatic measurements and sample changing is used to obtain the basic count data. Generation of the required theoretical decay and growth functions, varying with time, and the least squares solution of the overdetermined simultaneous count rate equations are done with a digital computer. Examples of the complex count rate equations which may be solved and results of a natural sample containing four alpha-emitting isotopes of thorium are illustrated. These methods facilitate the determination of the radioactive sources on the large scale required for many geochemical investigations.

Analytical Chemistry↗

Sediment laboratory quality-assurance project: studies of methods and materials

In August 1996 the U.S. Geological Survey initiated the Sediment Laboratory Quality-Assurance project. The Sediment Laboratory Quality Assurance project is part of the National Sediment Laboratory Quality-Assurance program. This paper addresses the fmdings of the sand/fme separation analysis completed for the single-blind reference sediment-sample project and differences in reported results between two different analytical procedures. From the results it is evident that an incomplete separation of fme- and sand-size material commonly occurs resulting in the classification of some of the fme-size material as sand-size material. Electron microscopy analysis supported the hypothesis that the negative bias for fme-size material and the positive bias for sand-size material is largely due to aggregation of some of the fine-size material into sand-size particles and adherence of fine-size material to the sand-size grains. Electron microscopy analysis showed that preserved river water, which was low in dissolved solids, specific conductance, and neutral pH, showed less aggregation and adhesion than preserved river water that was higher in dissolved solids and specific conductance with a basic pH. Bacteria were also found growing in the matrix, which may enhance fme-size material aggregation through their adhesive properties. Differences between sediment-analysis methods were also investigated as pan of this study. Suspended-sediment concentration results obtained from one participating laboratory that used a total-suspended solids (TSS) method had greater variability and larger negative biases than results obtained when this laboratory used a suspended-sediment concentration method. When TSS methods were used to analyze the reference samples, the median suspended sediment concentration percent difference was -18.04 percent. When the laboratory used a suspended-sediment concentration method, the median suspended-sediment concentration percent difference was -2.74 percent. The percent difference was calculated as follows: Percent difference = (( reported mass - known mass)/known mass ) X 100.

Conference Paper↗

Calculation of 230ThU isochrons, ages, and errors

If analytical errors are responsible for the scatter of points on a 230 Th- 234 U- 238 U isochron diagram, the isochron should be fitted by a technique that 1. weights the points according to their analytical errors and error correlations, and 2. either takes into account the presence of some of the same data in two coupled XY isochrons or (equivalently) uses a single, three-dimensional XYZ isochron. A method based on maximum-likelihood estimation is presented that fulfills these requirements, and the relevant equations for errors in age and initial 234 U 238 U "> 234 U 238 U are given. Equations for estimating the necessary isotope-ratio errors and error-correlations for both alpha-spectrometric and mass-spectrometric data are also developed.

Geochimica et Cosmochimica Acta↗

Estimation of αL, velocity, Kd and confidence limits from tracer injection test data

Bromide and boron were used as tracers during an injection experiment conducted at an artificial recharge facility near Stanton, Texas. The Ogallala aquifer at the Stanton site represents a heterogeneous alluvial environment and provides the opportunity to report scale dependent dispersivities at observation distances of 2 to 15 m in this setting. Values of longitudinal dispersivities are compared with other published values. Water samples were collected at selected depths both from piezometers and from fully screened observation wells at radii of 2, 5, 10 and 15 m. An exact analytical solution is used to simulate the concentration breakthrough curves and estimate longitudinal dispersivities and velocity parameters. Greater confidence can be placed on these data because the estimated parameters are error bounded using the bootstrap method. The non-conservative behavior of boron transport in clay rich sections of the aquifer were quantified with distribution coefficients by using bromide as a conservative reference tracer.

Texas↗

Spectroscopic remote sensing for material identification, vegetation characterization, and mapping

Identifying materials by measuring and analyzing their reflectance spectra has been an important procedure in analytical chemistry for decades. Airborne and space-based imaging spectrometers allow materials to be mapped across the landscape. With many existing airborne sensors and new satellite-borne sensors planned for the future, robust methods are needed to fully exploit the information content of hyperspectral remote sensing data. A method of identifying and mapping materials using spectral feature analyses of reflectance data in an expert-system framework called MICA (Material Identification and Characterization Algorithm) is described. MICA is a module of the PRISM (Processing Routines in IDL for Spectroscopic Measurements) software, available to the public from the U.S. Geological Survey (USGS) at http://pubs.usgs.gov/of/2011/1155/. The core concepts of MICA include continuum removal and linear regression to compare key diagnostic absorption features in reference laboratory/field spectra and the spectra being analyzed. The reference spectra, diagnostic features, and threshold constraints are defined within a user-developed MICA command file (MCF). Building on several decades of experience in mineral mapping, a broadly-applicable MCF was developed to detect a set of minerals frequently occurring on the Earth's surface and applied to map minerals in the country-wide coverage of the 2007 Afghanistan HyMap data set. MICA has also been applied to detect sub-pixel oil contamination in marshes impacted by the Deepwater Horizon incident by discriminating the C-H absorption features in oil residues from background vegetation. These two recent examples demonstrate the utility of a spectroscopic approach to remote sensing for identifying and mapping the distributions of materials in imaging spectrometer data.

Conference Paper↗

Digital mapping techniques '04, workshop proceedings : May 16-19, 2004, Portland, Oregon

The Digital Mapping Techniques ‘04 (DMT’04) workshop was attended by about 100 technical experts from 40 agencies, universities, and private companies, including representatives from 22 state geological surveys (see Appendix A). This workshop was similar in nature to the previous seven meetings, held in Lawrence, Kansas (Soller, 1997), in Champaign, Illinois (Soller, 1998), in Madison, Wisconsin (Soller, 1999), in Lexington, Kentucky (Soller, 2000), in Tuscaloosa, Alabama (Soller, 2001), in Salt Lake City, Utah (Soller, 2002), and in Millersville, Pennsylvania (Soller, 2003). This year’s meeting was hosted by the Oregon Department of Geology and Mineral Industries, from May 16–19, 2004, on the Portland State University campus in Portland, Oregon. As in the previous meetings, the objective was to foster informal discussion and exchange of technical information. This objective was well met, as attendees continued to share and exchange knowledge and information, and to renew friendships and collegial work begun at past DMT workshops. All the DMT workshops have been coordinated by the Association of American State Geologists (AASG) and U.S. Geological Survey (USGS) Data Capture Working Group, which was formed in August 1996, to support the AASG and the USGS in their effort to build a National Geologic Map Database (see Soller, Berg, and Stamm, this volume, and http://ngmdb.usgs.gov/info/standards/datacapt/ ). The Working Group was formed because increased production efficiencies, standardization, and quality of digital map products were needed for the National database—and for the State and Federal geological surveys—to provide more high-quality digital maps to the public. At the 2004 meeting, oral and poster presentations and special discussion sessions emphasized: 1) methods for creating and publishing map products (here, “publishing” includes Web-based release); 2) field data capture software and techniques; 3) digital cartographic techniques; 4) migration of digital maps into ArcGIS Geodatabase format; 5) analytical GIS techniques; 6) continued development of the National Geologic Map Database; and 7) progress toward building and implementing a standard geologic map data model and standard science language for the U.S. and for North America.

Open-File Report↗

General external uncertainty models of three-plane intersection point for 3D absolute accuracy assessment of lidar point cloud

The traditional practice to assess accuracy in lidar data involves calculating RMSEz (root mean square error of the vertical component). Accuracy assessment of lidar point clouds in full 3D (dimension) is not routinely performed. The main challenge in assessing accuracy in full 3D is how to identify a conjugate point of a ground-surveyed checkpoint in the lidar point cloud with the smallest possible uncertainty value. Relatively coarse point-spacing in airborne lidar data makes it challenging to determine a conjugate point accurately. As a result, a substantial unwanted error is added to the inherent positional uncertainty of the lidar data. Unless we keep this additional error small enough, the 3D accuracy assessment result will not properly represent the inherent uncertainty. We call this added error “external uncertainty,” which is associated with conjugate point identification. This research developed a general external uncertainty model using three-plane intersections and accounts for several factors (sensor precision, feature dimension, and point density). This method can be used for lidar point cloud data from a wide range of sensor qualities, point densities, and sizes of the features of interest. The external uncertainty model was derived as a semi-analytical function that takes the number of points on a plane as an input. It is a normalized general function that can be scaled by smooth surface precision (SSP) of a lidar system. This general uncertainty model provides a quantitative guideline on the required conditions for the conjugate point based on the geometric features. Applications of external uncertainty model was demonstrated using various lidar point cloud data from US Geological Survey (USGS) 3D Elevation Program (3DEP) library to determine the valid conditions for a conjugate point from three-plane.

Remote Sensing↗

Determination of the δ 15 N and δ 18 O of nitrate in solids; RSIL lab code 2897

The purpose of the Reston Stable Isotope Laboratory (RSIL) lab code 2897 is to determine the δ 15 N and δ 18 O of nitrate (NO 3 - ) in solids. The NO 3 - fraction of the nitrogen species is dissolved by water (called leaching) and can be analyzed by the bacterial method covered in RSIL lab code 2900. After leaching, the δ 15 N and δ 18 O of the dissolved NO 3 - is analyzed by conversion of the NO 3 - to nitrous oxide (N 2 O), which serves as the analyte for mass spectrometry. A culture of denitrifying bacteria is used in the enzymatic conversion of NO 3 - to N 2 O, which follows the pathway shown in equation 1: NO 3 - → NO 2 - → NO → 1/2 N 2 O (1) Because the bacteria Pseudomonas aureofaciens lack N 2 O reductive activity, the reaction stops at N 2 O, unlike the typical denitrification reaction that goes to N 2 . After several hours, the conversion is complete, and the N 2 O is extracted from the vial, separated from volatile organic vapor and water vapor by an automated -65 °C isopropanol-slush trap, a Nafion drier, a CO 2 and water removal unit (Costech #021020 carbon dioxide absorbent with Mg(ClO 4 ) 2 ), and trapped in a small-volume trap immersed in liquid nitrogen with a modified Finnigan MAT (now Thermo Scientific) GasBench 2 introduction system. After the N 2 O is released, it is further purified by gas chromatography before introduction to the isotope-ratio mass spectrometer (IRMS). The IRMS is a Thermo Scientific Delta V Plus continuous flow IRMS (CF-IRMS). It has a universal triple collector, consisting of two wide cups with a narrow cup in the middle; it is capable of simultaneously measuring mass/charge ( m/z ) of the N 2 O molecule 44, 45, and 46. The ion beams from these m/z values are as follows: m/z = 44 = N 2 O = 14 N 14 N 16 O; m/z = 45 = N 2 O = 14 N 15 N 16 O or 14 N 14 N 17 O; m/z = 46 = N 2 O = 14 N 14 N 18 O. The 17 O contributions to the m/z 44 and m/z 45 ion beams are accounted for before δ 15 N values are reported.

Techniques and Methods↗

White-Nose Syndrome Diagnostic Laboratory Network handbook

When responding to a wildlife disease outbreak, managers depend on consistent and clear data to make decisions. However, diagnostic methods for detecting pathogens of wildlife often lack the level of procedural and interpretational standardization that occurs in the investigation of human and domestic animal diseases. This lack of standardization can hamper diagnostic reliability in two ways. First is the inappropriate application of tests to new species or in situations that are outside of the original (in other words, validated) purpose. Second is the use of laboratory-specific modifications or analytical parameters without thorough investigation of how those changes affect result comparisons across institutions or the ability to make broader conclusions about pathogen or disease. White-nose syndrome (WNS) is a disease caused by the fungal pathogen Pseudogymnoascus destructans ( Pd ), which has spread rapidly and is causing population-level declines in some species of North American bats. During the last decade, quantitative polymerase chain reaction (qPCR) has become the most common method of testing for Pd because of qPCR’s speed, accuracy, and simplicity across a wide range of invasive and noninvasive sample types. Its widespread use by many State, Federal, Provincial, and academic institutions has inevitably led to variations in methodology and interpretation among laboratories. The progressive geographic spread of fungus and disease has also led to sampling contexts and strategies that differ from those for which the qPCR assay was originally developed and validated. These factors have resulted in inconsistencies among results tested in different laboratories and, subsequently, confusion for managers and decision makers. To address these challenges, the WNS National Response Team Diagnostic Working Group launched a project congruent with increased calls for the harmonization of wildlife disease diagnostic results, and reporting standards across disparate methodologies and laboratories. Beginning in 2019, interlaboratory testing was done to better understand how variations to Pd qPCR methodology affect diagnostic consistency and to reassess the assay’s fit for purpose in new testing contexts. This information led to expanded conversations within the Diagnostic Working Group related to best practices in Pd qPCR diagnostic testing, the development of common interpretation language for classifying test results, and the incorporation of that language into an updated WNS case definition. This handbook is the resulting product and is intended to help further harmonize Pd qPCR diagnostic testing by establishing recommendations related to voluntary participation in a WNS Diagnostic Laboratory Network, documenting the currently (2022) practiced Pd qPCR methodologies, discussing general best practices for molecular diagnostics and laboratory networks, and elaborating on the epidemiologic and diagnostic basis of the agreed-upon classification language for Pd qPCR results. Through this voluntary, consensus-based approach to diagnostic harmonization, this work aims to improve the confidence of management agencies in reported Pd qPCR results and can serve as an example of national diagnostic coordination for other unregulated wildlife diseases.

Techniques and Methods↗

Evaluation of soil zone processes and a novel radiocarbon correction approach for groundwater with mixed sources

Estimates of groundwater age based on 14C is often limited by the uncertainty in geochemical processes that alter the 14C concentration measured in water and the composition (δ13C and 14C) of carbon sources needed to appropriately parametrize 14C adjustment models. Estimated ages for samples that contain a mixture of young and old groundwater will be particularly sensitive to model parametrization as relatively small additions of modern 14C from recent recharge can mask the presence and amount of old groundwater. A novel multi-model approach based on inverse geochemical modeling and lumped parameter modeling of age tracers (3H, 3Hetrit, and SF6) was used to better constrain 14C dilution caused by dissolution of carbonates in the unsaturated zone or shallow parts of the Glacial aquifer, which extends over 2000 miles across the northern contiguous United States. Calibration of 14C inverse geochemical models to LPM computed 14C concentrations in modern water indicated that 14C of soil zone and shallow aquifer carbonates were not 14C-dead (0 pmC), as is typically assumed for 14C correction models. 14C of such carbonates was on average about 53 pmC (ranged 0-110 pmC, n = 72). This information was used to correct 14C concentrations for water recharged entirely before 1950 and water that is a mixture of pre- and post-1950 water. The multi-model approach developed here was compared to an analytical 14C-adjustment model (Revised Fontes and Garnier) that assumed solid carbonates were 14C-dead. 14C corrections using the analytical adjustment model tended to over-correct final 14C concentrations by 21 pmC and underestimates mean ages by 40% for groundwater mixtures. In fact, 14C corrections based on analytical model yielded negative ages (14C > 120 pmC) in nearly 36% of mixed samples. This work presents a new approach to constraining 14C corrections and age estimates of mixtures of young and old groundwater. The new method is applied to three well networks distributed across the spatially expansive Glacial aquifer.

Glacier aquifer system↗

Sources and magnitude of bias associated with determination of polychlorinated biphenyls in environmental samples

Recently compiled data on the composition of commercial Aroclor mixtures and ECD (electron capture detector) response factors for all 209 PCB congeners are used to develop estimates of the bias associated with determination of polychlorinated blphenyis. During quantitation of multlcomponent peaks by congener-specific procedures error is introduced because of variable ECD response to isomeric PCBs. Under worst case conditions, the magnitude of this bias can range from less than 2% to as much as 600%. Multicomponent peaks containing the more highly and the lower chlorinated congeners experience the most bias. For this reason, quantitation of Σ PCB in Aroclor mixtures dominated by these species (e.g. 1016) are potentially subject to the greatest error. Comparison of response factor data for ECDs from two laboratories shows that the sign and magnitude of calibration bias for a given multicomponent peak is variable and depends, in part, on the response characteristics of individual detectors. By using the most abundant congener (of each multicomponent peak) for purposes of calibration, one can reduce the maximum bias to less than 55%. Moreover, due to cancellation of errors, the bias resulting from summation of all peak concentrations (i.e. Σ PCB) becomes vanishingly small (<1.2%). In contrast, bias associated with determination of Σ PCB as Aroclor equivalents by the traditional Aroclor method is potentially large (>200%) and highly variable in sign and magnitude. In this case, bias originates not only from the incomplete chromatographic resolution of PCB congeners but also the overlapping patterns of the Aroclor mixtures. Together these results illustrate the advantages of the congener-specific method of PCB quantitation over the traditional Aroclor Method and the extreme difficulty of estimating bias incurred by the latter procedure on a post hoc basis.

Analytical Chemistry↗

Extending the MODPATH algorithm to rectangular unstructured grids

The recent release of MODFLOW-USG, which allows model grids to have irregular, unstructured connections, requires a modification of the particle-tracking algorithm used by MODPATH. This paper describes a modification of the semi-analytical particle-tracking algorithm used by MODPATH that allows it to be extended to rectangular-based unstructured grids by dividing grid cells with multi-cell face connections into sub-cells. The new method will be incorporated in the next version of MODPATH which is currently under development.

Ground Water↗