USGS Science⌕ Search

SEARCH · USGS Science

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 1,297 records · Page 72Linked to original sources
USGS2023

Foundations of modeling resilience of tidal saline wetlands to sea-level rise along the U.S. Pacific Coast

Context Tidal saline wetlands (TSWs) are highly threatened from climate-change effects of sea-level rise. Studies of TSWs along the East Coast U.S. and elsewhere suggest significant likely losses over coming decades but needed are analytic tools gauged to Pacific Coast U.S. wetlands. Objectives We predict the impacts of sea-level rise (SLR) on the elevation capital (vertical) and migration potential (lateral) resilience of TSWs along the Pacific Coast U.S. over the period 2020 to 2150 under a 1.5-m SLR scenario, and identified TSWs at risk of most rapid loss of resilience. Here, we define vertical resilience as the amount of elevation capital and lateral resilience as the amount of TSW displacement area relative to existing area. Methods We used Bayesian network (BN) modeling to predict changes in resilience of TSWs as probabilities which can be useful in risk analysis and risk management. We developed the model using a database sample of 26 TSWs with 147 sediment core samples, among 16 estuary drainage areas along coastal California, Oregon, and Washington. Results We found that all TSW sites would lose at least 50% of their elevation capital resilience by 2060 to just before 2100, and 100% by 2070 to 2130, depending on the site. Under a 1.5-m sea-level rise scenario, nearly all sites in California will lose most or all of their lateral migration resilience. Resilience losses generally accelerated over time. In the BN model, elevation capital resilience is most sensitive to elevation capital at time t, mean tide level at time t, and change in sea level from time 0 to time t. Conclusions All TSW sites were projected with declines in resilience. Our model can further aid decision-making such as prioritizing sites for potential management adaptation strategies. We also identified variables most influencing resilience predictions and thus those potentially prioritized for monitoring or development of strategies to prevent loss regionally.

California, Oregon, Washington↗
USGS2025

A comprehensive assessment of membrane bioreactor contaminant removal efficacy through analytical chemistry, fish exposures, and microbiome characterization

Treated municipal wastewater effluent is an important pathway for Contaminants of Emerging Concern (CEC) to enter aquatic ecosystems. As the aging wastewater infrastructure in many industrialized countries requires upgrades or replacement, assessing new treatment technologies in the context of CEC effects may provide additional support for science-based resource management. Here, we used three lines of evidence, analytical chemistry, fish exposure experiments, and fish and water microbiome analysis, to assess the effectiveness of membrane bioreactor treatment (MBR) to replace traditional activated sludge treatment. To do this, we sampled a municipal wastewater treatment plant with a split wastewater stream, a portion of which was treated with an MBR and another via an oxidation ditch (OXI). The OXI and MBR treatments substantially reduced most measured CECs compared to the primary effluent (PRI). Only pesticides and some pharmaceuticals were recalcitrant to both secondary treatment methods. While ammonia toxicity of PRI prevented its inclusion in fish exposure experiments, exposure of fish with waters from the OXI or MBR treated wastewater produced only subtle biological differences with no adverse apical outcomes. These findings were consistent with low chemically derived exposure: activity ratios for OXI and MBR. Microbiome analysis of fish and wastewater highlighted the significant reduction of microbial abundance and diversity in the MBR treatment compared to all other treatments. The comparable removal efficacy of CECs in MBR makes it an attractive alternative to traditional OXI, especially when MBR may eliminate the need for tertiary treatment for wastewater disinfection.

Environmental Toxicology and Chemistry↗
USGS2019

Trends in streamflow and concentrations and flux of nutrients and total suspended solids in the Upper White River at Muncie, near Nora, and near Centerton, Indiana

The U.S. Geological Survey (USGS), in cooperation with The Nature Conservancy, completed a study to estimate and assess trends in streamflow and annual mean concentrations and flux of nutrients (nitrate plus nitrite, total Kjeldahl nitrogen, and total phosphorus) and total suspended solids at three USGS streamgages (hereafter referred to as “study gages”) on the Upper White River at Muncie (USGS station 03347000), near Nora (USGS station 03351000), and near Centerton (USGS station 03354000), Indiana. Water-quality data used in the analyses were collected by several agencies between calendar years 1991 and 2017, and streamflow (discharge) data were collected by the USGS. For most of the water-quality constituents, there were suitable data to facilitate an analysis of the 26-year period extending from calendar years 1991 to 2017 (water years 1992 to 2017); however, shorter analytical periods were necessary for total Kjeldahl nitrogen for the study gages at Muncie and near Centerton and for total suspended solids for the study gage near Centerton. Temporal trends in streamflows at the study gages for the period extending from water years 1978 to 2017 were assessed using Exploration and Graphics for RivEr Trends (EGRET) and Mann-Kendall and Pettitt tests. With just one exception, the annual maximum and mean daily streamflows and the annual minimum 7-day mean streamflows at the study gages demonstrated upward trends (increasing streamflows) in the EGRET analyses. The exception was the annual 7-day minimum streamflow at the study gage near Nora, which indicated no trend. Mann-Kendall tests also indicated that the average trend for the annual maximum daily, annual mean daily, and annual 7-day minimum streamflow statistics between water years 1978 and 2017 was upward at each of the study gages; however, only the trends in the annual mean daily streamflows at the study gage at Muncie and the annual maximum daily streamflows at the study gages near Nora and near Centerton were statistically significant at a 0.05 probability level. The Pettitt tests indicated that a statistically significant step trend (abrupt change) in annual mean daily streamflows occurred at each of the study gages around water year 2001. The seasonal distributions of total suspended solids, total phosphorus, nitrate plus nitrite, and total Kjeldahl nitrogen concentrations at the study gages were evaluated to identify patterns and other distinguishing characteristics by examining boxplots of concentrations as a function of month of the year. Seasonal distributions of nitrate plus nitrite concentrations and total suspended solids concentrations differed from each other but were generally similar among the three study gages for a given constituent. Median concentrations of nitrate plus nitrite were highest during the January–June months, whereas median concentrations of total suspended solids were highest during June and July. Seasonal distributions of total phosphorus concentrations were similar at the study gages near Nora and near Centerton, but the seasonal distribution was noticeably different at the study gage at Muncie, which had monthly median concentrations that were substantially lower than at the two downstream study gages (near Nora and near Centerton). The seasonal distribution of total Kjeldahl nitrogen concentrations differed in pattern among the three study gages; however, in general, some of the higher monthly median total Kjeldahl nitrogen concentrations at each study gage were associated with the late spring and summer periods. The Weighted Regressions on Time, Discharge, and Season (WRTDS) method implemented in EGRET was used to estimate water-year annual mean daily concentrations and flux of nutrients and total suspended solids, as well as estimates of concentrations and flux that were “normalized” to remove the effect of year-to-year variation in streamflow. The approximate coefficients of determination for the WRTDS regression models ranged from a high of 0.82 for total phosphorus for the study gage near Centerton to a low of 0.19 for nitrate plus nitrite for the study gage near Nora. Loads and yields of total suspended solids, total phosphorus, nitrate plus nitrite, and total Kjeldahl nitrogen were estimated for analytical periods consisting of the longest periods of concurrent record at the three study gages. Loads of each of the constituents increased sequentially from the most upstream study gage to the most downstream study gage; however, the same was not true for yields. The highest yields of total suspended solids, total phosphorus, and total Kjeldahl nitrogen occurred at the most upstream study gage (at Muncie); however, the highest yield of nitrate plus nitrite occurred at the most downstream study gage (near Centerton). WRTDS bootstrap tests were used to assess the magnitude, direction, and likelihood of changes in annual flow-normalized mean daily concentrations and flux of total suspended solids, total phosphorus, nitrate plus nitrite, and total Kjeldahl nitrogen at the study gages between water years 1997 and 2017. Changes in flow-normalized concentrations and flux of the constituents between water years 1997 and 2017 were mostly downward (decreasing). The exceptions were likely to highly likely upward (increasing) changes in (1) flow-normalized annual mean daily concentration and annual flux for total suspended solids and total phosphorus at the study gage at Muncie, (2) flow-normalized annual mean daily total phosphorus concentration at the study gage near Centerton, (3) flow-normalized annual flux of total phosphorus at the study gage near Centerton, and (4) flow-normalized annual mean daily nitrate plus nitrite concentration at the study gage near Centerton. Although an upward change in flow-normalized nitrate plus nitrite concentrations was likely at the study gage near Centerton, flow-normalized annual flux of nitrate plus nitrite at that study gage was determined to have a highly likely downward change. EGRET and Exploration and Graphics for RivEr Trends Confidence Intervals (EGRETci) analyses can be used to improve our understanding of how concentrations and flux change as functions of time and streamflow, as well as provide information on how the relations between streamflow and constituent concentrations have changed within the calendar year between any 2 years included in the analyses. Examples of those uses, illustrating changes between calendar years 1992 and 2017, were given for total suspended solids concentrations at the study gage near Nora and for nitrate plus nitrite concentrations at the study gage near Centerton.

Indiana↗
USGS1958

Determination of the oxidation state of uranium in apatite and phosphorite deposits

Geological and mineralogical evidence indicate that the uranium present in apatite may proxy for calcium in the mineral structure as U(IV). An experimental investigation was conducted and chemical evidence was obtained that establishes the presence of U(IV) in apatite. The following analytical procedure was developed for the determination of U(IV). Carbonatefluorapatite is dissolved in 1.5 M orthophosphoric acid at a temperature of 5°C or slightly below and fluorapatite is dissolved in cold 1.2 M hydrochloric acid (approximately 5°C) containing 1.5 g of hydroxylamine hydrochloride per 100 ml. Uranium(IV) is precipitated by cupferron using titanium as a carrier. The uranium in the precipitate is separated by use of the ethyl acetate extraction procedure and determined fluorimetrically. The validity and the limitations of the method have been established by spike experiments.

Geochimica et Cosmochimica Acta↗
USGS2001

Local site effects and dynamic soil behavior

Amplitudes of seismic waves increase significantly as they pass through soft soil layers near the earth's surface. This phenomenon, commonly known as site amplification, is a major factor influencing the extent of damage on structures. It is crucial that site amplification is accounted for when designing structures on soft soils. The characteristics of site amplification at a given site can be estimated by analytical models, as well as field tests. Analytical models require as inputs the geometry of all soil layers from surface to bedrock, their dynamic properties (e.g. density, wave velocity, damping), and the incident bedrock motions. Field tests involve recording and analyzing the dynamic response of sites to artificial excitations, ambient forces, and actual earthquakes. The most reliable estimates of site amplification are obtained by analyzing the recorded motions of the site during strong earthquakes. This paper presents a review of the types and the generating mechanisms of site amplification, and the models and methods that are used to characterize them from earthquake records. ?? 2001 Published by Elsevier Science Ltd.

Soil Dynamics and Earthquake Engineering↗
USGS2005

Computation of type curves for flow to partially penetrating wells in water-table aquifers

Evaluation of Neuman's analytical solution for flow to a well in a homogeneous, anisotropic, water-table aquifer commonly requires large amounts of computation time and can produce inaccurate results for selected combinations of parameters. Large computation times occur because the integrand of a semi-infinite integral involves the summation of an infinite series. Each term of the series requires evaluation of the roots of equations, and the series itself is sometimes slowly convergent. Inaccuracies can result from lack of computer precision or from the use of improper methods of numerical integration. In this paper it is proposed to use a method of numerical inversion of the Laplace transform solution, provided by Neuman, to overcome these difficulties. The solution in Laplace space is simpler in form than the real-time solution; that is, the integrand of the semi-infinite integral does not involve an infinite series or the need to evaluate roots of equations. Because the integrand is evaluated rapidly, advanced methods of numerical integration can be used to improve accuracy with an overall reduction in computation time. The proposed method of computing type curves, for which a partially documented computer program (WTAQ1) was written, was found to reduce computation time by factors of 2 to 20 over the time needed to evaluate the closed-form, real-time solution.

Ground Water↗
USGS2019

Spatially-structured statistical network models for landscape genetics

A basic understanding of how the landscape impedes, or creates resistance to, the dispersal of organisms and hence gene flow is paramount for successful conservation science and management. Spatially structured ecological networks are often used to represent spatial landscape‐genetic relationships, where nodes represent individuals or populations and resistance to movement is represented using non‐binary edge weights. Weights are typically assigned or estimated by the user, rather than observed, and validating such weights is challenging. We provide a synthesis of current methods used to estimate edge weights and an overview of common model types, stressing the advantages and disadvantages of each approach and their ability to model landscape‐genetic data. We further explore a set of spatial‐statistical methods that provide ecologists with alternative approaches for modeling spatially explicit processes that may affect genetic structure. This includes an overview of spatial autoregressive models, with a particular focus on how correlation and partial correlation are used to represent neighborhood structure with the inverse of the covariance matrix (i.e., precision matrix). We then demonstrate how to model resistance by specifying an appropriate statistical model on the nodes, conditioned on the edge weights, through the precision matrix. This integration of network ecology and spatial statistics provides a practical analytical framework for landscape‐genetic studies. The results can be used to make statistical inferences about the relative importance of individual landscape characteristics, such as the vegetative cover, hillslope, or the presence of roads or rivers, on gene flow. In addition, the R code we include allows readers to explore landscape‐genetic structure in their own datasets, which will potentially provide new insights into the evolutionary processes that generated ecological networks, as well as valuable information about the optimal characteristics of conservation corridors.

Ecological Monographs↗
USGS2002

Sequestration of priority pollutant PAHs from sediment pore water employing semipermeable membrane devices

Semipermeable membrane devices (SPMDs) were employed to sample sediment pore water in static exposure studies under controlled laboratory conditions using (control pond and formulated) sediments fortified with 15 priority pollutant polycyclic aromatic hydrocarbons (PPPAHs). The sediment fortification level of 750 ng/g was selected on the basis of what might be detected in a sediment sample from a contaminated area. The sampling interval consisted of 0, 4, 7, 14, and 28 days for each study. The analytical methodologies, as well as the extraction and sample cleanup procedures used in the isolation, characterization, and quantitation of 15 PPPAHs at different fortification levels in SPMDs, water, and sediment were reported previously (Williamson, M.S. Thesis, University of Missouri - Columbia, USA; Williamson et al., Chemosphere (This issue - PII: S0045-6535(02)00394-6)) and used for this project. Average (mean) extraction recoveries for each PPPAH congener in each matrix are reported and discussed. No procedural blank extracts (controls) were found to contain any PPPAH residues above the method quantitation limit, therefore, no matrix interferences were detected. The focus of this publication is to demonstrate the ability to sequester environmental contaminants, specifically PPPAHs, from sediment pore water using SPMDs and two different types of fortified sediment.

Chemosphere↗
USGS2001

GRA prospectus: optimizing design and management of protected areas

Protected areas comprise one major type of global conservation effort that has been in the form of parks, easements, or conservation concessions. Though protected areas are increasing in number and size throughout tropical ecosystems, there is no systematic method for optimally targeting specific local areas for protection, designing the protected area, and monitoring it, or for guiding follow-up actions to manage it or its surroundings over the long run. Without such a system, conservation projects often cost more than necessary and/or risk protecting ecosystems and biodiversity less efficiently than desired. Correcting these failures requires tools and strategies for improving the placement, design, and long-term management of protected areas. The objective of this project is to develop a set of spatially based analytical tools to improve the selection, design, and management of protected areas. In this project, several conservation concessions will be compared using an economic optimization technique. The forest land use portfolio model is an integrated assessment that measures investment in different land uses in a forest. The case studies of individual tropical ecosystems are developed as forest (land) use and preservation portfolios in a geographic information system (GIS). Conservation concessions involve a private organization purchasing development and resource access rights in a certain area and retiring them. Forests are put into conservation, and those people who would otherwise have benefited from extracting resources or selling the right to do so are compensated. Concessions are legal agreements wherein the exact amount and nature of the compensation result from a negotiated agreement between an agent of the conservation community and the local community. Funds are placed in a trust fund, and annual payments are made to local communities and regional/national governments. The payments are made pending third-party verification that the forest expanse and quality have been maintained.

Open-File Report↗
USGS2016

Montana StreamStats—A method for retrieving basin and streamflow characteristics in Montana: Chapter A in Montana StreamStats

The U.S. Geological Survey (USGS) provides streamflow characteristics and other related information needed by water-resource managers to protect people and property from floods, plan and manage water-resource activities, and protect water quality. Streamflow characteristics provided by the USGS, such as peak-flow and low-flow frequencies for streamflow-gaging stations, are frequently used by engineers, flood forecasters, land managers, biologists, and others to guide their everyday decisions. In addition to providing streamflow characteristics at streamflow-gaging stations, the USGS also develops regional regression equations and drainage area-adjustment methods for estimating streamflow characteristics at locations on ungaged streams. Regional regression equations can be complex and often require users to determine several basin characteristics, which are physical and climatic characteristics of the stream and its drainage basin. Obtaining these basin characteristics for streamflow-gaging stations and ungaged sites traditionally has been time consuming and subjective, and led to inconsistent results. StreamStats is a Web-based geographic information system application that was created by the USGS to provide users with access to an assortment of analytical tools that are useful for water-resource planning and management. StreamStats allows users to easily obtain streamflow and basin characteristics for USGS streamflow-gaging stations and user-selected locations on ungaged streams. The USGS, in cooperation with Montana Department of Transportation, Montana Department of Environmental Quality, and Montana Department of Natural Resources and Conservation, completed a study to develop a StreamStats application for Montana, compute streamflow characteristics at streamflow-gaging stations, and develop regional regression equations to estimate streamflow characteristics at ungaged sites. Chapter A of this Scientific Investigations Report describes the Montana StreamStats application and the datasets, streamflow-gaging stations, streamflow characteristics, and regression equations (as described fully in Chapters B through G of this report) that are used for development of the StreamStats application for Montana.

Montana↗
USGS1976

Analytical precision of one-sixth order semiquantitative spectrographic analysis

Over 2,700 separate analyses were made on 22 different geologic samples selected to cover wide concentration ranges for the 30 elements studied. The precision for low, medium, and high concentrations of each element determined is represented by superimposed frequency diagrams and displays the frequency of reported values occurring within one-sixth order reporting intervals about the mean.Results are derived from mixed analyst and instrument conditions with no controls enforced other than randomization of samples. The precision of the six-step (one-sixth order) semiquantitative spectrographic analysis utilized by the Denver-based mobile laboratories exceeds that which is necessary for exploration geochemistry. Disallowing results obtained near the detection levels, the repeatability of the method, in general, is shown to be within on adjoining reporting interval on each side of the mean, 83 percent of the time, and within two adjoining reporting intervals on each side of the mean 96 percent of the time.

Circular↗
USGS2002

Assessment of relative accuracy in the determination of organic matter concentrations in aquatic systems

Accurate determinations of total (TOC), dissolved (DOC) and particulate (POC) organic carbon concentrations are critical for understanding the geochemical, environmental, and ecological roles of aquatic organic matter. Of particular significance for the drinking water industry, TOC measurements are the basis for compliance with US EPA regulations. The results of an interlaboratory comparison designed to identify problems associated with the determination of organic matter concentrations in drinking water supplies are presented. The study involved 31 laboratories and a variety of commercially available analytical instruments. All participating laboratories performed well on samples of potassium hydrogen phthalate (KHP), a compound commonly used as a standard in carbon analysis. However, problems associated with the oxidation of difficult to oxidize compounds, such as dodecylbenzene sulfonic acid and caffeine, were noted. Humic substances posed fewer problems for analysts. Particulate organic matter (POM) in the form of polystyrene beads, freeze-dried bacteria and pulverized leaf material were the most difficult for all analysts, with a wide range of performances reported. The POM results indicate that the methods surveyed in this study are inappropriate for the accurate determination of POC and TOC concentration. Finally, several analysts had difficulty in efficiently separating inorganic carbon from KHP solutions, thereby biasing DOC results.

Journal of Environmental Monitoring↗
USGS2016

Continuous-flow centrifugation to collect suspended sediment for chemical analysis

Recent advances in suspended-sediment monitoring tools and surrogate technologies have greatly improved the ability to quantify suspended-sediment concentrations and to estimate daily, seasonal, and annual suspended-sediment fluxes from rivers to coastal waters. However, little is known about the chemical composition of suspended sediment, and how it may vary spatially between water bodies and temporally within a single system owing to climate, seasonality, land use, and other natural and anthropogenic drivers. Many water-quality contaminants, such as organic and inorganic chemicals, nutrients, and pathogens, preferentially partition in sediment rather than water. Suspended sediment-bound chemical concentrations may be undetected during analysis of unfiltered water samples, owing to small water sample volumes and analytical limitations. Quantification of suspended sediment‑bound chemical concentrations is needed to improve estimates of total chemical concentrations, chemical fluxes, and exposure levels of aquatic organisms and humans in receiving environments. Despite these needs, few studies or monitoring programs measure the chemical composition of suspended sediment, largely owing to the difficulty in consistently obtaining samples of sufficient quality and quantity for laboratory analysis. A field protocol is described here utilizing continuous‑flow centrifugation for the collection of suspended sediment for chemical analysis. The centrifuge used for development of this method is small, lightweight, and portable for the field applications described in this protocol. Project scoping considerations, deployment of equipment and system layout options, and results from various field and laboratory quality control experiments are described. The testing confirmed the applicability of the protocol for the determination of many inorganic and organic chemicals sorbed on suspended sediment, including metals, pesticides, polycyclic aromatic hydrocarbons, and polychlorinated biphenyls. The particle-size distribution of the captured sediment changes to a more fine-grained sample during centrifugation, and the necessity to account for this change when extrapolating chemical concentrations on the centrifuged sediment sample to the environmental water system is discussed. The data produced using this method will help eliminate a data gap of suspended sediment-bound chemical concentrations, and will support management decisions, such as chemical source-control efforts or in-stream restoration activities. When coupled with streamflow and sediment flux data, it will improve estimates of riverine chemical fluxes, and will aid in assessing the importance and impacts of suspended sediment-bound chemicals to downstream freshwater and coastal marine ecosystems.

Techniques and Methods↗
USGS2004

A theory for modeling ground-water flow in heterogeneous media

Construction of a ground-water model for a field area is not a straightforward process. Data are virtually never complete or detailed enough to allow substitution into the model equations and direct computation of the results of interest. Formal model calibration through optimization, statistical, and geostatistical methods is being applied to an increasing extent to deal with this problem and provide for quantitative evaluation and uncertainty analysis of the model. However, these approaches are hampered by two pervasive problems: 1) nonlinearity of the solution of the model equations with respect to some of the model (or hydrogeologic) input variables (termed in this report system characteristics) and 2) detailed and generally unknown spatial variability (heterogeneity) of some of the system characteristics such as log hydraulic conductivity, specific storage, recharge and discharge, and boundary conditions. A theory is developed in this report to address these problems. The theory allows construction and analysis of a ground-water model of flow (and, by extension, transport) in heterogeneous media using a small number of lumped or smoothed system characteristics (termed parameters). The theory fully addresses both nonlinearity and heterogeneity in such a way that the parameters are not assumed to be effective values. The ground-water flow system is assumed to be adequately characterized by a set of spatially and temporally distributed discrete values, ?, of the system characteristics. This set contains both small-scale variability that cannot be described in a model and large-scale variability that can. The spatial and temporal variability in ? are accounted for by imagining ? to be generated by a stochastic process wherein ? is normally distributed, although normality is not essential. Because ? has too large a dimension to be estimated using the data normally available, for modeling purposes ? is replaced by a smoothed or lumped approximation y?. (where y is a spatial and temporal interpolation matrix). Set y?. has the same form as the expected value of ?, y 'line' ? , where 'line' ? is the set of drift parameters of the stochastic process; ?. is a best-fit vector to ?. A model function f(?), such as a computed hydraulic head or flux, is assumed to accurately represent an actual field quantity, but the same function written using y?., f(y?.), contains error from lumping or smoothing of ? using y?.. Thus, the replacement of ? by y?. yields nonzero mean model errors of the form E(f(?)-f(y?.)) throughout the model and covariances between model errors at points throughout the model. These nonzero means and covariances are evaluated through third and fifth-order accuracy, respectively, using Taylor series expansions. They can have a significant effect on construction and interpretation of a model that is calibrated by estimating ?.. Vector ?.. is estimated as 'hat' ? using weighted nonlinear least squares techniques to fit a set of model functions f(y'hat' ?) to a. corresponding set of observations of f(?), Y. These observations are assumed to be corrupted by zero-mean, normally distributed observation errors, although, as for ?, normality is not essential. An analytical approximation of the nonlinear least squares solution is obtained using Taylor series expansions and perturbation techniques that assume model and observation errors to be small. This solution is used to evaluate biases and other results to second-order accuracy in the errors. The correct weight matrix to use in the analysis is shown to be the inverse of the second-moment matrix E(Y-f(y?.))(Y-f(y?.))', but the weight matrix is assumed to be arbitrary in most developments. The best diagonal approximation is the inverse of the matrix of diagonal elements of E(Y-f(y?.))(Y-f(y?.))', and a method of estimating this diagonal matrix when it is unknown is developed using a special objective function to compute 'hat' ?. When considered to be an estimate of f

Professional Paper↗
USGS1989

A numerical solution for the diffusion equation in hydrogeologic systems

The documentation of a computer code for the numerical solution of the linear diffusion equation in one or two dimensions in Cartesian or cylindrical coordinates is presented. Applications of the program include molecular diffusion, heat conduction, and fluid flow in confined systems. The flow media may be anisotropic and heterogeneous. The model is formulated by replacing the continuous linear diffusion equation by discrete finite-difference approximations at each node in a block-centered grid. The resulting matrix equation is solved by the method of preconditioned conjugate gradients. The conjugate gradient method does not require the estimation of iteration parameters and is guaranteed convergent in the absence of rounding error. The matrixes are preconditioned to decrease the steps to convergence. The model allows the specification of any number of boundary conditions for any number of stress periods, and the output of a summary table for selected nodes showing flux and the concentration of the flux quantity for each time step. The model is written in a modular format for ease of modification. The model was verified by comparison of numerical and analytical solutions for cases of molecular diffusion, two-dimensional heat transfer, and axisymmetric radial saturated fluid flow. Application of the model to a hypothetical two-dimensional field situation of gas diffusion in the unsaturated zone is demonstrated. The input and output files are included as a check on program installation. The definition of variables, input requirements, flow chart, and program listing are included in the attachments. (USGS)

Water-Resources Investigations Report↗
USGS2014

Oxygen isotope systematics in the aragonite-CO2-H2O-NaCl system up to 0.7 mol/kg ionic strength at 25 °C

To investigate the oxygen isotope systematics in the aragonite-CO2-H2O-NaCl system, witherite (BaCO3) was precipitated quasi-instantaneously and quantitatively from Na-Cl-Ba-CO2 solutions of seawater-like ionic strength (I = 0.7 mol/kg) at two pH values (~7.9 and ~10.6) at 25 °C. The oxygen isotope composition of the witherite and the dissolved inorganic carbon speciation in the starting solution were used to estimate the oxygen isotope fractionations between HCO3¯ and H2O as well as between CO3 2 and H2O. Given the analytical error on the oxygen isotope composition of the witherite and uncertainties of the parent solution pH and speciation, oxygen isotope fractionation between NaHCO3° and HCO3¯, as well as between NaCO3¯ and CO3 2, is negligible under the experimental conditions investigated. The influence of dissolved NaCl concentration on the oxygen isotope fractionation in the aragonite-CO2-H2O-NaCl system also was investigated at 25 °C. Aragonite was precipitated from Na-Cl-Ca-Mg-(B)-CO2 solutions of seawater-like ionic strength using passive CO2 degassing or constant addition methods. Based upon our new experimental observations and published experimental data from lower ionic strength solutions by Kim et al. (2007b), the equilibrium aragonite-water oxygen isotope fractionation factor is independent of the ionic strength of the parent solution up to 0.7 mol/kg. Hence, our study also suggests that the aragonite precipitation mechanism is not affected by the presence of sodium and chloride ions in the parent solution over the range of concentrations investigated.

Geochimica et Cosmochimica Acta↗
USGS2023

Wide-area debris field and seabed characterization of a deep ocean dump site surveyed by autonomous underwater vehicles

Disposal of industrial and hazardous waste in the deep ocean was a pervasive global practice near coastlines in the 20th century. Uncertainty in the quantity, location, and contents of dumped materials at historical disposal sites underscores ongoing risks to marine ecosystems and human health. This study presents analysis of a 150-km2 wide-area sidescan sonar survey conducted in March 2021 with two autonomous underwater vehicles (AUVs) at an offshore dumpsite in San Pedro Basin, California. Previous camera surveys located 60 barrels and other debris; sediment analysis in the region showed varying concentrations of the insecticidal chemical Dichlorodiphenyl-trichloroethane (DDT), for which an estimated 350-700 metric tons was discarded in the San Pedro Basin between 1947 and 1961. A lack of primary historical documents from Montrose Chemical or California Salvage Company that specify the DDT acid waste disposal method has contributed to the ambiguity surrounding whether the dumping method was via bulk discharge or containerized units. Barrels and debris visually detected during previous surveys by Remotely Operated Vehicles (ROVs) Jason and SuBastian and AUV Sentry served as the basis for ground truth testing and object classification using size and acoustic intensity characteristics. Target detection and classification algorithms were developed using image and signal processing techniques that resulted in the identification of over 74,000 debris targets within the survey region. Statistical, spectral, and machine learning techniques are used to characterize seabed variability and for bottom-type classification. These analytical techniques combined with the survey capabilities of AUVs provide a framework for efficient mapping and characterization of uncharted deep-water disposal sites.

California↗
USGS2018

Rethinking the use of seabed sediment temperature profiles to trace submarine groundwater flow

Submarine groundwater fluxes across the seafloor facilitate important hydrological and biogeochemical exchanges between oceans and seabed sediment, yet few studies have investigated spatially distributed groundwater fluxes in deep‐ocean environments such as continental slopes. Heat has been previously applied as a submarine groundwater tracer using an analytical solution to a heat flow equation assuming steady state conditions and homogeneous thermal conductivity. These assumptions are often violated in shallow seabeds due to ocean bottom temperature changes or sediment property variations. Here heat tracing analysis techniques recently developed for terrestrial settings are applied in concert to examine the influences of groundwater flow, ocean temperature changes, and seabed thermal conductivity variations on deep‐ocean sediment temperature profiles. Temperature observations from the sediment and bottom ocean water on the Scotian Slope off eastern Canada are used to demonstrate how simple thermal methods for tracing groundwater can be employed if more comprehensive techniques indicate that the simplifying assumptions are valid. The spatial distribution of the inferred groundwater fluxes on the slope suggests a downward groundwater flow system with recharge occurring over the upper‐middle slope and discharge on the lower slope. We speculate that the downward groundwater flow inferred on the Scotian Slope is due to density‐driven processes arising from underlying salt domes, in contrast with upward slope systems driven by geothermal convection. Improvements in the design of future submarine hydrogeological studies are proposed for thermal data collection and groundwater flow analysis, including new equations that quantify the minimum detectable flux magnitude for a given sensor accuracy and profile length.

Water Resources Research↗