USGS Science⌕ Search

SEARCH · USGS Science

Results for “Techniques and 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 919 records · Page 51Linked to original sources

Book review: Bird census techniques, Second edition

Conservation concerns, federal mandates to monitor birds, and citizen science programs have spawned a variety of surveys that collect information on bird populations. Unfortunately, all too frequently these surveys are poorly designed and use inappropriate counting methods. Some of the flawed approaches reflect a lack of understanding of statistical design; many ornithologists simply are not aware that many of our most entrenched counting methods (such as point counts) cannot appropriately be used in studies that compare densities of birds over space and time. It is likely that most of the readers of The Condor have participated in a bird population survey that has been criticized for poor sampling methods. For example, North American readers may be surprised to read in Bird Census Techniques that the North American Breeding Bird Survey 'is seriously flawed in its design,' and that 'Analysis of trends is impossible from points that are positioned along roads' (p. 109). Our conservation efforts are at risk if we do not acknowledge these concerns and improve our survey designs. Other surveys suffer from a lack of focus. In Bird Census Techniques, the authors emphasize that all surveys require clear statements of objectives and an understanding of appropriate survey designs to meet their objectives. Too often, we view survey design as the realm of ornithologists who know the life histories and logistical issues relevant to counting birds. This view reflects pure hubris: survey design is a collaboration between ornithologists, statisticians, and managers, in which goals based on management needs are met by applying statistical principles for design to the biological context of the species of interest. Poor survey design is often due to exclusion of some of these partners from survey development. Because ornithologists are too frequently unaware of these issues, books such as Bird Census Techniques take on added importance as manuals for educating ornithologists about the relevance of survey design and methods and the often subtle interdisciplinary nature of surveys. Review info: Bird Census Techniques, Second Edition. By Colin J. Bibby, Neil D. Burgess, David A. Hill, and Simon H. Mustoe. 2000. Academic Press, London, UK. xvii 1 302 pp. ISBN 0- 12-095831-7.

Condor↗

Evaluating freshwater mussel sampling methodologies using a simulation model

Field surveys form the basis of many research efforts and are the foundation for estimates of population size and density that inform conservation and management practices for imperiled species. As a result, evaluating the performance of different survey methods across a range of conditions that may be encountered in the field can increase understanding of the time and effort that may be required to ensure that survey results are sufficiently accurate and reliable for conservation goals. We used a spatially explicit agent-based model to simulate four commonly used freshwater mussel field survey methodologies: simple random sampling (SRS), transect random sampling (TRS), adaptive cluster sampling (ACS), and qualitative timed searches (QTS) to investigate the influence of sampling method, spatial distribution, and mussel density on the performance ( i.e., accuracy, precision, and detection rate) of survey techniques. Our analysis suggests that mussel density, spatial distribution, and sampling effort influence sampling accuracy, precision, and species detection for all sampling methods. QTS produces highly variable catch-per-unit-effort (CPUE) metrics when mussels are dense and/or clustered, indicating the technique may be unreliable as a proxy for density. Quantitative methods like SRS and TRS may be well-suited for estimating population characteristics, but a high level of effort may be needed to obtain reasonable accuracy when mussels occur at low densities. ACS may be more efficient for mussels at low densities, but it can be challenging to plan for the level of effort required to complete an ACS protocol. Designing an ecological survey requires careful consideration of research objectives and available resources. Future research may consider the performance of qualitative and quantitative surveys in combination as a means of overcoming some of the practical challenges of applying individual survey methods.

Ecological Indicators↗

Capture-recapture of white-tailed deer using DNA from fecal pellet-groups

Traditional methods for estimating white-tailed deer population size and density are affected by behavioral biases, poor detection in densely forested areas, and invalid techniques for estimating effective trapping area. We evaluated a noninvasive method of capture—recapture for white-tailed deer (Odocoileus virginianus) density estimation using DNA extracted from fecal pellets as an individual marker and for gender determination, coupled with a spatial detection function to estimate density (spatially explicit capture—recapture, SECR). We collected pellet groups from 11 to 22 January 2010 at randomly selected sites within a 1-km2 area located on Arnold Air Force Base in Coffee and Franklin counties, Tennessee. We searched 703 10-m radius plots and collected 352 pellet-group samples from 197 plots over five two-day sampling intervals. Using only the freshest pellets we recorded 140 captures of 33 different animals (15M:18F). Male and female densities were 1.9 (SE = 0.8) and 3.8 (SE = 1.3) deer km-2, or a total density of 5.8 deer km-2 (14.9 deer mile-2). Population size was 20.8 (SE = 7.6) over a 360-ha area, and sex ratio was 1.0 M: 2.0 F (SE = 0.71). We found DNA sampling from pellet groups improved deer abundance, density and sex ratio estimates in contiguous landscapes which could be used to track responses to harvest or other management actions.

Tennessee↗

Techniques for computing discharge at four Navigation Dams on the Illinois and Des Plaines Rivers in Illinois

Techniques for computing discharge are developed for Brandon Road Dam on the Des Plaines River and for Dresden Island, Marseilles, and Starved Rock Dams on the Illinois River. At Brandon Road Dam, streamflow is regulated by the operation of Tainter gates and headgates. At Dresden Island, Marseilles, and Starved Rock Dams, only Tainter gates are operated to regulate streamflow. The locks at all dams are equipped with culvert valves that are used to fill and empty the lock. The techniques facilitate determination of discharge at locations along the upper Illinois Waterway where no streamflow-gaging stations exist. The techniques are also useful for computing low flows when the water-surface slope between control structures on the river approaches zero and traditional methods of determining discharge based on slope are unsatisfactory. Two techniques can be used to compute discharge at the dams--gate ratings and tailwater ratings . A gate ratingdescribes the relation between discharge, gate opening, tailwater stage, and headwater stage. A tailwater rating describes the relation between tailwater stage and discharge. Gate ratings for Tainter gates at Dresden Island, Marseilles, and Starved Rock Dams are based on a total of 78 measurements of discharge that range from 569 to 86,400 cubic feet per second. Flood hydrographs developed from the gate ratings and Lockmaster records of gate opening and stage compare closely with streamflow records published for nearby streamflow-gaging stations. Additional measurements are needed to verify gate ratings for Tainter gates and headgates at Brandon Road Dam after the dam rehabilitation is completed. Extensive leakage past deteriorated headgates and sluice gates contributed to uncertainty in the ratings developed for this dam. A useful tailwater rating is developed for Marseilles Dam. Tailwater ratings for Dresden Island Dam and Starved Rock Dam are of limited use because of varying downstream channel-storage conditions. A tailwater rating could not be developed for Brandon Road Dam because its tailwater pool is substantially affected by the headwater pool of Dresden Island Dam.

Water-Resources Investigations Report↗

Assessment of PDMS-water partition coefficients: implications for passive environmental sampling of hydrophobic organic compounds

Solid-phase microextraction (SPME) has shown potential as an in situ passive-sampling technique in aquatic environments. The reliability of this method depends upon accurate determination of the partition coefficient between the fiber coating and water (K f ). For some hydrophobic organic compounds (HOCs), K f values spanning 4 orders of magnitude have been reported for polydimethylsiloxane (PDMS) and water. However, 24% of the published data examined in this review did not pass the criterion for negligible depletion, resulting in questionable K f values. The range in reported K f is reduced to just over 2 orders of magnitude for some polychlorinated biphenyls (PCBs) when these questionable values are removed. Other factors that could account for the range in reported K f , such as fiber-coating thickness and fiber manufacturer, were evaluated and found to be insignificant. In addition to accurate measurement of K f , an understanding of the impact of environmental variables, such as temperature and ionic strength, on partitioning is essential for application of laboratory-measured K f values to field samples. To date, few studies have measured K f for HOCs at conditions other than at 20 degrees or 25 degrees C in distilled water. The available data indicate measurable variations in K f at different temperatures and different ionic strengths. Therefore, if the appropriate environmental variables are not taken into account, significant error will be introduced into calculated aqueous concentrations using this passive sampling technique. A multiparameter linear solvation energy relationship (LSER) was developed to estimate log K f in distilled water at 25 degrees C based on published physicochemical parameters. This method provided a good correlation (R2 = 0.94) between measured and predicted log K f values for several compound classes. Thus, an LSER approach may offer a reliable means of predicting log K f for HOCs whose experimental log K f values are presently unavailable. Future research should focus on understanding the impact of environmental variables on K f . Obtaining the data needed for an LSER approach to estimate K f for all environmentally relevant HOCs would be beneficial to the application of SPME as a passive-sampling technique.

Environmental Science & Technology↗

Effects of simple acid leaching of crushed and powdered geological materials on high-precision Pb isotope analyses

We present new results of simple acid leaching experiments on the Pb isotope composition of USGS standard reference material powders and on ocean island basalt whole rock splits and powders. Rock samples were leached with cold 6 N HCl in an ultrasonic bath, then on a hot plate, and washed with ultrapure H 2 O before sample digestion in HF‐HNO 3 and chromatographic purification of Pb. Lead isotope analyses were measured by Tl‐doped MC‐ICPMS. Intrasession and intersession analytical reproducibilities of repeated analyses of both synthetic Pb solutions and Pb from single digests of chemically processed natural samples were generally better than 100 ppm (2 SD). The comparison of leached and unleached samples shows that leaching consistently removes variable amounts of contaminants that differ in Pb isotopic composition for different starting materials. For repeated digests of a single sample, analyses of leached samples reproduce better than those of unleached ones, confirming that leaching effectively removes most of the heterogeneously distributed extraneous Pb. Nevertheless, the external reproducibility of leached samples is still up to an order of magnitude worse than that of Pb solution standards (∼100 ppm). More complex leaching methods employed by earlier studies yield Pb isotope ratios within error of those produced by our method and at similar levels of reproducibility, demonstrating that our simple leaching method is as effective as more complex leaching techniques. Therefore, any Pb isotope heterogeneity among multiple leached digests of samples in excess of the external reproducibility is attributed to inherent isotopic heterogeneity of the sample. The external precision of ∼100 ppm (2 SD) achieved for Pb isotope ratio determination by Tl‐doped MC‐ICPMS is thus sufficient for most rocks. The full advantage of the most precise Pb isotope analytical methods is only realized in cases where the natural isotopic heterogeneity among samples in a studied suite is substantially below 100 ppm.

Geochemistry, Geophysics, Geosystems↗

An evaluation of methods for estimating decadal stream loads

Effective management of water resources requires accurate information on the mass, or load of water-quality constituents transported from upstream watersheds to downstream receiving waters. Despite this need, no single method has been shown to consistently provide accurate load estimates among different water-quality constituents, sampling sites, and sampling regimes. We evaluate the accuracy of several load estimation methods across a broad range of sampling and environmental conditions. This analysis uses random sub-samples drawn from temporally-dense data sets of total nitrogen, total phosphorus, nitrate, and suspended-sediment concentration, and includes measurements of specific conductance which was used as a surrogate for dissolved solids concentration. Methods considered include linear interpolation and ratio estimators, regression-based methods historically employed by the U.S. Geological Survey, and newer flexible techniques including Weighted Regressions on Time, Season, and Discharge (WRTDS) and a generalized non-linear additive model. No single method is identified to have the greatest accuracy across all constituents, sites, and sampling scenarios. Most methods provide accurate estimates of specific conductance (used as a surrogate for total dissolved solids or specific major ions) and total nitrogen – lower accuracy is observed for the estimation of nitrate, total phosphorus and suspended sediment loads. Methods that allow for flexibility in the relation between concentration and flow conditions, specifically Beale’s ratio estimator and WRTDS, exhibit greater estimation accuracy and lower bias. Evaluation of methods across simulated sampling scenarios indicate that (1) high-flow sampling is necessary to produce accurate load estimates, (2) extrapolation of sample data through time or across more extreme flow conditions reduces load estimate accuracy, and (3) WRTDS and methods that use a Kalman filter or smoothing to correct for departures between individual modeled and observed values benefit most from more frequent water-quality sampling.

Journal of Hydrology↗

Methods of practice and guidelines for using survey-grade global navigation satellite systems (GNSS) to establish vertical datum in the United States Geological Survey

Geodetic surveys have evolved through the years to the use of survey-grade (centimeter level) global positioning to perpetuate and post-process vertical datum. The U.S. Geological Survey (USGS) uses Global Navigation Satellite Systems (GNSS) technology to monitor natural hazards, ensure geospatial control for climate and land use change, and gather data necessary for investigative studies related to water, the environment, energy, and ecosystems. Vertical datum is fundamental to a variety of these integrated earth sciences. Essentially GNSS surveys provide a three-dimensional position x, y, and z as a function of the North American Datum of 1983 ellipsoid and the most current hybrid geoid model. A GNSS survey may be approached with post-processed positioning for static observations related to a single point or network, or involve real-time corrections to provide positioning "on-the-fly." Field equipment required to facilitate GNSS surveys range from a single receiver, with a power source for static positioning, to an additional receiver or network communicated by radio or cellular for real-time positioning. A real-time approach in its most common form may be described as a roving receiver augmented by a single-base station receiver, known as a single-base real-time (RT) survey. More efficient real-time methods involving a Real-Time Network (RTN) permit the use of only one roving receiver that is augmented to a network of fixed receivers commonly known as Continually Operating Reference Stations (CORS). A post-processed approach in its most common form involves static data collection at a single point. Data are most commonly post-processed through a universally accepted utility maintained by the National Geodetic Survey (NGS), known as the Online Position User Service (OPUS). More complex post-processed methods involve static observations among a network of additional receivers collecting static data at known benchmarks. Both classifications provide users flexibility regarding efficiency and quality of data collection. Quality assurance of survey-grade global positioning is often overlooked or not understood and perceived uncertainties can be misleading. GNSS users can benefit from a blueprint of data collection standards used to ensure consistency among USGS mission areas. A classification of GNSS survey qualities provide the user with the ability to choose from the highest quality survey used to establish objective points with low uncertainties, identified as a Level I, to a GNSS survey for general topographic control without quality assurance, identified as a Level IV. A Level I survey is strictly limited to post-processed methods, whereas Level II, Level III, and Level IV surveys integrate variations of a RT approach. Among these classifications, techniques involving blunder checks and redundancy are important, and planning that involves the assessment of the overall satellite configuration, as well as terrestrial and space weather, are necessary to ensure an efficient and quality campaign. Although quality indicators and uncertainties are identified in post-processed methods using CORS, the accuracy of a GNSS survey is most effectively expressed as a comparison to a local benchmark that has a high degree of confidence. Real-time and post-processed methods should incorporate these "trusted" benchmarks as a check during any campaign. Global positioning surveys are expected to change rapidly in the future. The expansion of continuously operating reference stations, combined with newly available satellite signals, and enhancements to the conterminous geoid, are all sufficient indicators for substantial growth in real-time positioning and quality thereof.

Techniques and Methods↗

Refinements to the Graves and Pitarka (2010) Broadband Ground Motion Simulation Method

This brief article describes refinements to the Graves and Pitarka (2010) broadband ground motion simulation methodology (GP2010 hereafter) that have been implemented in version 14.3 of the SCEC Broadband Platform (BBP). The updated version of our method on the current SCEC BBP is referred to as GP14.3. Our simulation technique is a hybrid approach that combines low-­‐frequency and high-­‐frequency motions computed with different methods into a single broadband response. The separate low-­‐ and high-­‐frequency components have traditionally been called “deterministic” and “stochastic”, respectively; however, this nomenclature is an oversimplification. In reality, the low-­‐frequency approach includes many stochastic elements, and likewise, the high-­‐frequency approach includes many deterministic elements (e.g., Pulido and Kubo, 2004; Hartzell et al., 2005; Liu et al., 2006; Frankel, 2009; Graves and Pitarka, 2010; Mai et al., 2010). While the traditional terminology will likely remain in use by the broader modeling community, in this paper we will refer to these using the generic terminology “low-­‐frequency” and “high-­‐ frequency” approaches. Furthermore, one of the primary goals in refining our methodology is to provide a smoother and more consistent transition between the low-­‐ and high-­‐ frequency calculations, with the ultimate objective being the development of a single unified modeling approach that can be applied over a broad frequency band. GP2010 was validated by modeling recorded strong motions from four California earthquakes. While the method performed well overall, several issues were identified including the tendency to over-­‐predict the level of longer period (2-­‐5 sec) motions and the effects of rupture directivity. The refinements incorporated in GP14.3 are aimed at addressing these issues with application to the simulation of earthquakes in Western US (WUS). These refinements include the addition of a deep weak zone (details in following section) to the rupture characterization and allowing perturbations in the correlation of rise time and rupture speed with the specified slip distribution. Additionally, we have extended the parameterization of GP14.3 so that it is also applicable for simulating Eastern North America (ENA) earthquakes. This work has been guided by the comprehensive set of validation studies described in Goulet and Abrahamson (2014) and Dreger et al. (2014). The GP14.3 method shows improved performance relative to GP2010, and we direct the interested reader to Dreger et al. (2014) for a detailed assessment of the current methodology. In this paper, we concentrate on describing the modifications in more detail, and also discussing additional refinements that are currently being developed.

California↗

A practical guide to understanding and validating complex models using data simulations

Biologists routinely fit novel and complex statistical models to push the limits of our understanding. Examples include, but are not limited to, flexible Bayesian approaches (e.g. BUGS, stan), frequentist and likelihood-based approaches (e.g. packages lme4 ) and machine learning methods. These software and programs afford the user greater control and flexibility in tailoring complex hierarchical models. However, this level of control and flexibility places a higher degree of responsibility on the user to evaluate the robustness of their statistical inference. To determine how often biologists are running model diagnostics on hierarchical models, we reviewed 50 recently published papers in 2021 in the journal Nature Ecology & Evolution , and we found that the majority of published papers did not report any validation of their hierarchical models, making it difficult for the reader to assess the robustness of their inference. This lack of reporting likely stems from a lack of standardized guidance for best practices and standard methods. Here, we provide a guide to understanding and validating complex models using data simulations. To determine how often biologists use data simulation techniques, we also reviewed 50 recently published papers in 2021 in the journal Methods Ecology & Evolution . We found that 78% of the papers that proposed a new estimation technique, package or model used simulations or generated data in some capacity (18 of 23 papers); but very few of those papers (5 of 23 papers) included either a demonstration that the code could recover realistic estimates for a dataset with known parameters or a demonstration of the statistical properties of the approach. To distil the variety of simulations techniques and their uses, we provide a taxonomy of simulation studies based on the intended inference. We also encourage authors to include a basic validation study whenever novel statistical models are used, which in general, is easy to implement. Simulating data helps a researcher gain a deeper understanding of the models and their assumptions and establish the reliability of their estimation approaches. Wider adoption of data simulations by biologists can improve statistical inference, reliability and open science practices.

Methods in Ecology and Evolution↗

Comparison of dialysis membrane diffusion samplers and two purging methods in bedrock wells

Collection of ground-water samples from bedrock wells using low-flow purging techniques is problematic because of the random spacing, variable hydraulic conductivity, and variable contamination of contributing fractures in each well's open interval. To test alternatives to this purging method, a field comparison of three ground-water-sampling techniques was conducted on wells in fractured bedrock at a site contaminated primarily with volatile organic compounds. Constituent concentrations in samples collected with a diffusion sampler constructed from dialysis membrane material were compared to those in samples collected from the same wells with a standard low-flow purging technique and a hybrid (high-flow/low-flow) purging technique. Concentrations of trichloroethene, cis-1,2-dichloroethene, vinyl chloride, calcium, chloride, and alkalinity agreed well among samples collected with all three techniques in 9 of the 10 wells tested. Iron concentrations varied more than those of the other parameters, but their pattern of variation was not consistent. Overall, the results of nonparametric analysis of variance testing on the nine wells sampled twice showed no statistically significant difference at the 95-percent confidence level among the concentrations of volatile organic compounds or inorganic constituents recovered by use of any of the three sampling techniques.

Conference Paper↗

The roughness of natural terrain: A planetary and remote sensing perspective

We examine the various methods and parameters in common use for quantifying and reporting surface topographic "roughness." It is shown that scale-dependent roughness parameters are almost always required, though not widely used. We suggest a method of standardizing the parameters that are computed and reported so that topographic data gathered by different workers using different field techniques can be directly and easily intercompared. We illustrate the proposed method by analyzing topographic data from 60 different surfaces gathered by five different groups and examine the information for common features. We briefly discuss the implications of our analysis for studies of planetary surface roughness, lander safety, and radar remote sensing modeling and analysis.

Journal of Geophysical Research E: Planets↗

The centenary of IAVCEI 1919–2019 and beyond: The people, places, and things of volcano geodesy

Over the first century of the International Association of Volcanology and Chemistry of the Earth’s Interior (IAVCEI), volcano geodesy grew from roots as an accidental and incidental system of measurements to an important method for monitoring volcanic activity and forecasting eruptions. The first practitioners in volcano geodesy were experts in other disciplines, and it was not until the latter half of the twentieth century that specialists in the field emerged—scientists who developed new methods, measured geodetic change at volcanoes, and quantitatively interpreted the results in terms of magmatic processes. Much of the early work in the field was restricted to a few volcanoes and involved techniques that had been adapted from other applications; relatively few methods were developed specifically for use on volcanoes. These volcanoes, however, provided the natural laboratories needed to advance the field. By the start of the twenty-first century, geodetic studies, especially using space-based techniques, contributed to the recognition of deformation and gravity change at hundreds of volcanoes on Earth. In coming years, IAVCEI researchers will focus on comprehensive exploitation of the growing volumes of geodetic data to better model, forecast, and track activity at volcanoes worldwide. Meanwhile, the field needs to become more diverse, better representing people who live in the shadows of volcanoes around the globe.

Bulletin of Volcanology↗

A semiobjective method for condensing classifications

A simple, semiobjective method is described to reduce the number of groups in a classification to an arbitrary level without losing contact with the geologic information contained in the evolving groups. The method, operated in a stepwise or cyclic manner, employs some of the commonly used numerical techniques, but avoids strict adherence to them to obtain geologically more meaningful results. The method is illustrated in a facies study of the upper Paleozoic rocks of southeastern Utah. ?? 1973 Plenum Publishing Corporation.

Journal of the International Association for Mathe↗

Regional intensity attenuation models for France and the estimation of magnitude and location of historical earthquakes

Intensity assignments for 33 calibration earthquakes were used to develop intensity attenuation models for the Alps, Armorican, Provence, Pyrenees and Rhine regions of France. Intensity decreases with ?? most rapidly in the French Alps, Provence and Pyrenees regions, and least rapidly in the Armorican and Rhine regions. The comparable Armorican and Rhine region attenuation models are aggregated into a French stable continental region model and the comparable Provence and Pyrenees region models are aggregated into a Southern France model. We analyse MSK intensity assignments using the technique of Bakun & Wentworth, which provides an objective method for estimating epicentral location and intensity magnitude MI. MI for the 1356 October 18 earthquake in the French stable continental region is 6.6 for a location near Basle, Switzerland, and moment magnitude M is 5.9-7.2 at the 95 per cent (??2??) confidence level. MI for the 1909 June 11 Trevaresse (Lambesc) earthquake near Marseilles in the Southern France region is 5.5, and M is 4.9-6.0 at the 95 per cent confidence level. Bootstrap resampling techniques are used to calculate objective, reproducible 67 per cent and 95 per cent confidence regions for the locations of historical earthquakes. These confidence regions for location provide an attractive alternative to the macroseismic epicentre and qualitative location uncertainties used heretofore. ?? 2006 The Authors Journal compilation ?? 2006 RAS.

Geophysical Journal International↗

Guidelines for using the Delphi Technique to develop habitat suitability index curves

Habitat Suitability Index (SI) curves are one method of presenting species habitat suitability criteria. The curves are often used with the Habitat Evaluation Procedures (HEP) and are necessary components of the Instream Flow Incremental Methodology (IFIM) (Armour et al. 1984). Bovee (1986) described three categories of SI curves or habitat suitability criteria based on the procedures and data used to develop the criteria. Category I curves are based on professional judgment, with 1ittle or no empirical data. Both Category II (utilization criteria) and Category III (preference criteria) curves have as their source data collected at locations where target species are observed or collected. Having Category II and Category III curves for all species of concern would be ideal. In reality, no SI curves are available for many species, and SI curves that require intensive field sampling often cannot be developed under prevailing constraints on time and costs. One alternative under these circumstances is the development and interim use of SI curves based on expert opinion. The Delphi technique (Pill 1971; Delbecq et al. 1975; Linstone and Turoff 1975) is one method used for combining the knowledge and opinions of a group of experts. The purpose of this report is to describe how the Delphi technique may be used to develop expert-opinion-based SI curves.

FWS/OBS↗

A comparison of solute-transport solution techniques based on inverse modelling results

Five common numerical techniques (finite difference, predictor-corrector, total-variation-diminishing, method-of-characteristics, and modified-method-of-characteristics) were tested using simulations of a controlled conservative tracer-test experiment through a heterogeneous, two-dimensional sand tank. The experimental facility was constructed using randomly distributed homogeneous blocks of five sand types. This experimental model provides an outstanding opportunity to compare the solution techniques because of the heterogeneous hydraulic conductivity distribution of known structure, and the availability of detailed measurements with which to compare simulated concentrations. The present work uses this opportunity to investigate how three common types of results-simulated breakthrough curves, sensitivity analysis, and calibrated parameter values-change in this heterogeneous situation, given the different methods of simulating solute transport. The results show that simulated peak concentrations, even at very fine grid spacings, varied because of different amounts of numerical dispersion. Sensitivity analysis results were robust in that they were independent of the solution technique. They revealed extreme correlation between hydraulic conductivity and porosity, and that the breakthrough curve data did not provide enough information about the dispersivities to estimate individual values for the five sands. However, estimated hydraulic conductivity values are significantly influenced by both the large possible variations in model dispersion and the amount of numerical dispersion present in the solution technique.Five common numerical techniques (finite difference, predictor-corrector, total-variation-diminishing, method-of-characteristics, and modified-method-of-characteristics) were tested using simulations of a controlled conservative tracer-test experiment through a heterogeneous, two-dimensional sand tank. The experimental facility was constructed using randomly distributed homogeneous blocks of five sand types. This experimental model provides an outstanding opportunity to compare the solution techniques because of the heterogeneous hydraulic conductivity distribution of known structure, and the availability of detailed measurements with which to compare simulated concentrations. The present work uses this opportunity to investigate how three common types of results - simulated breakthrough curves, sensitivity analysis, and calibrated parameter values - change in this heterogeneous situation, given the different methods of simulating solute transport. The results show that simulated peak concentrations, even at very fine grid spacings, varied because of different amounts of numerical dispersion. Sensitivity analysis results were robust in that they were independent of the solution technique. They revealed extreme correlation between hydraulic conductivity and porosity, and that the breakthrough curve data did not provide enough information about the dispersivities to estimate individual values for the five sands. However, estimated hydraulic conductivity values are significantly influenced by both the large possible variations in model dispersion and the amount of numerical dispersion present in the solution technique.

Conference Paper↗

Toward an efficient framework for remote sensing of river bathymetry: Comparing sensors and algorithms on an inaccessible proglacial river in Alaska

Remote sensing can provide reliable information on river depths and this approach might be particularly valuable in areas that are difficult to survey via conventional field methods. In this study, we assessed the potential to map the bathymetry of an inaccessible proglacial river in Alaska from both aerial orthophotos and a multispectral satellite image. In addition, we evaluated a variety of depth retrieval algorithms with different input data requirements, including some methods that require field measurements of water depth for calibration and other techniques that can be applied even when such field data are not available. These approaches might enable more efficient use of remote sensing methods by resource management agencies. Our results suggest that bathymetric mapping along the turquoise-colored river we examined was not only feasible but highly accurate ( R 2 up to 0.94) for both types of image data. Algorithms that use paired observations of depth and reflectance to train depth retrieval models were the most accurate, with errors on the order of 15%–20% and little or no bias. Alternative techniques based on hydraulic and statistical concepts also led to strong agreement between predicted and observed depths but were more susceptible to systematic biases toward under- or over-estimation of depth. In contrast to clear-flowing streams, bathymetric mapping in this environment was enabled by a direct relationship between the depth and brightness of the water due to scattering by suspended sediment. In selecting an appropriate depth retrieval method, a compromise might need to be reached between the level of field effort invested and the accuracy of the resulting image-derived bathymetry. Standalone software for implementing these techniques is freely available.

Alaska↗