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 1,063 records · Page 59Linked to original sources

Precise estimation of repeating earthquake moment: Example from parkfield, california

We offer a new method for estimating the relative size of repeating earthquakes using the singular value decomposition (SVD). This method takes advantage of the highly coherent waveforms of repeating earthquakes and arrives at far more precise and accurate descriptions of earthquake size than standard catalog techniques allow. We demonstrate that uncertainty in relative moment estimates is reduced from ??75% for standard coda-duration techniques employed by the network to an uncertainty of ??6.6% when the SVD method is used. This implies that a single-station estimate of moment using the SVD method has far less uncertainty than the whole-network estimates of moment based on coda duration. The SVD method offers a significant improvement in our ability to describe the size of repeating earthquakes and thus an opportunity to better understand how they accommodate slip as a function of time.

Bulletin of the Seismological Society of America↗

Sampling protocol for monitoring abiotic and biotic characteristics of mountain ponds and lakes

This document describes field techniques and procedures used for sampling mountain ponds and lakes. These techniques and procedures will be used primarily to monitor, as part of long-term programs in National Parks and other protected areas, the abiotic and biotic characteristics of naturally occurring permanent montane lentic systems up to 75 ha in surface area. However, the techniques and procedures described herein also can be used to sample temporary or ephemeral montane lentic sites. Each Standard Operating Procedure (SOP) section addresses a specific component of the limnological investigation, and describes in detail field sampling methods pertaining to parameters to be measured for each component.

Techniques and Methods↗

User's guide to SAC, a computer program for computing discharge by slope-area method

This user's guide contains information on using the slope-area program, SAC. SAC can be used to compute peak flood discharges from measurements of high-water marks along a stream reach. The Slope-area method used by the program is the U.S. Geological Survey (USGS) procedure presented in Techniques of Water Resources Investigations of the U.S. Geological Survey, beok 3, chapter A2, "Measurement of Peak Discharge by the Slope-Area Method." The program uses input files that have formats compatible with those used by the water-surface profile program (WSPRO) described in the Federal Highways Administration publication FHWA-IP-89-027. The guide briefly describes the slope-area method documents the input requirements and the output produced, and demonstrates use of SAC.

Open-File Report↗

Hydrogeophysical methods for analyzing aquifer storage and recovery systems

Hydrogeophysical methods are presented that support the siting and monitoring of aquifer storage and recovery (ASR) systems. These methods are presented as numerical simulations in the context of a proposed ASR experiment in Kuwait, although the techniques are applicable to numerous ASR projects. Bulk geophysical properties are calculated directly from ASR flow and solute transport simulations using standard petrophysical relationships and are used to simulate the dynamic geophysical response to ASR. This strategy provides a quantitative framework for determining site‐specific geophysical methods and data acquisition geometries that can provide the most useful information about the ASR implementation. An axisymmetric, coupled fluid flow and solute transport model simulates injection, storage, and withdrawal of fresh water (salinity ∼500 ppm) into the Dammam aquifer, a tertiary carbonate formation with native salinity approximately 6000 ppm. Sensitivity of the flow simulations to the correlation length of aquifer heterogeneity, aquifer dispersivity, and hydraulic permeability of the confining layer are investigated. The geophysical response using electrical resistivity, time‐domain electromagnetic (TEM), and seismic methods is computed at regular intervals during the ASR simulation to investigate the sensitivity of these different techniques to changes in subsurface properties. For the electrical and electromagnetic methods, fluid electric conductivity is derived from the modeled salinity and is combined with an assumed porosity model to compute a bulk electrical resistivity structure. The seismic response is computed from the porosity model and changes in effective stress due to fluid pressure variations during injection/recovery, while changes in fluid properties are introduced through Gassmann fluid substitution.

Groundwater↗

Artisanal mining river dredge detection using SAR: A method comparison

Challenges exist in monitoring artisanal and small-scale mining (ASM) activities, given their dynamic and often informal nature. ASM takes form through various techniques and scales, including riverine dredging, which often targets the abundant alluvial gold deposits in South America. Remote sensing offers a solution to improve data collection, regulation, and monitoring of the more mobile and elusive ASM activities and their impacts. Mapping ASM riverine dredges using Synthetic Aperture Radar (SAR) is one of the application areas least explored. Three semi-automated detection approaches using Sentinel-1 SAR are compared on their ability to identify dredges with minimal false positives. The methods are: (i) Search for Unidentified Maritime Objects (SUMO), an established method for large ocean ship detection; and two techniques specifically developed for riverine environments that are introduced in this paper: (ii) a local detection method; and (iii) a global threshold method. A visual interpretation of SAR data with the inclusion of optical high-resolution data are used to generate a validation dataset. Results show it is possible to semi-automatically detect riverine dredge using SAR and that a local detection method provides the best balance between sensitivity and precision and has the lowest risk of error. Future improvements may consider further automation, more discriminatory variables, and analyzing the methods in different environments and at higher spatial resolutions.

Remote Sensing↗

Liquid chromatographic method for determining the concentration of bisazir in water

Barrier dams, traps, and lampricides are the techniques currently used by the Great Lakes Fishery Commission to control sea lampreys {Petromyzon marinug) in the Great Lakes. To augment these control techniques, a sterile-male-release research program was initiated at the Lake Huron Biological Station. Male sea lampreys were sterilized by intraperitoneal injection of the chemical sterilant P,P-bis(1-aziridinyl)- N methylphosphinothioic amide (bisazir). An analytical method was needed to quantitate the concentration of bisazir in water and to routinely verify that bisazir (>25 μg/L) does not persist in the treated effluent discharged from the sterilization facility to Lake Huron. A rapid, accurate, and sensitive liquid chromatographic (LC) method was developed for determining bisazir in water. Bisazir was dissolved in Lake Huron water; extracted and concentrated on a C 18 solid-phase extraction column; eluted with methanol; and quantitated by reversed-phase LC using a Cis column, amobile phase of 70% water and 30% methanol (v/v), and UV detection (205 nm). Bisazir retention time was 7-8 min; total run time was about 20 min. Method detection limit for bisazir dissolved in Lake Huron water was about 15 μg/L. Recovery from Lake Huron water fortified with bisazir at 100 μg/L was 94% (95% confidence interval, 90.2-98.2%).

Journal of AOAC International↗

Evaluation of laser ablation double-focusing SC-ICPMS for “common” lead isotopic measurements in silicate glasses and mineral

An analytical method for the in situ measurement of “common” Pb isotope ratios in silicate glasses and minerals using a 193-nm excimer laser ablation (LA) system with a double-focusing single-collector (SC)-ICPMS is presented and evaluated as a possible alternative to multiple-collector (MC)-ICPMS. This LA-SC-ICPMS technique employs fast-scanning ion deflectors to sequentially place a series of flat-topped isotope peaks into a single ion-counting detector at a fixed accelerating voltage and magnetic field strength. Reference materials (including NIST, MPI-DING, and USGS glasses) are used to identify two analytical artifacts on the Pb isotope ratios (expressed here as heavier/lighter isotopes) when corrected for mass bias relative to NIST SRM610. The first artifact is characterized by anomalously low Pb isotope ratios (~0.1%/AMU) when SRM610 is analyzed in raster mode as an unknown at small spot sizes (<25 µm), which may indicate that (1) SRM610 is isotopically heterogeneous on a small length scale and/or (2) there is a non-spectral matrix effect on the Pb isotope ratios related to differences in spot size. The second artifact is characterized by anomalously high Pb isotope ratios (<0.1%/AMU) for NIST SRM612 (in raster mode) and some Fe-rich glass reference materials (BCR-2G, GOR132-G, and T1-G). These offsets are thought to be caused by one or more non-spectral matrix effects related to differences in the ablation behavior, composition, or physical properties of these reference materials compared to the bracketing SRM610 standard. The precision (±2SD) of our LA-SC-ICPMS Pb isotopic measurements is similar to ( 207 Pb/ 206 Pb and 208 Pb/ 206 Pb, or 20X Pb/ 206 Pb) or better than ( 206 Pb/ 204 Pb, 207 Pb/ 204 Pb, and 208 Pb/ 204 Pb, or 20X Pb/ 204 Pb) a series of published studies that used a different type of SC-ICPMS and obtained a factor of ~3-4 higher sensitivity for Pb. An increase in the sensitivity of our LA-SC-ICPMS would likely improve the precision of the 20X Pb/ 206 Pb and 20X Pb/ 204P b ratios for low-Pb materials (<5 ppm), possibly making the technique broadly similar to LA-MC-ICPMS (particularly compared to methods that rely upon at least one ion-counting detector). Further improvement in the precision of the 20X Pb/ 206 Pb and 20X Pb/ 204 Pb ratios for high-Pb materials (>5 ppm) by LA-SC-ICPMS is unlikely, and in this case, LA-MC-ICPMS remains the preferable analytical technique.

Journal of Analytical Atomic Spectrometry↗

Borehole activation analysis by delayed and capture gamma rays using a 252Cf neutron source

Theoretical analysis and experimental comparison of the radiative capture and delayed gamma-ray activation techniques indicate the latter to be more efficient for the detection of copper, whereas the radiative capture method is preferable for nickel. A conservative lower detection limit for both copper and nickel is '-0.5%. Borehole spectra by both techniques were made in a copper- and nickel-bearing gabbro, utilizing a Ge(Li) propane-cooled detector. Al, Mn, Na, Mg, Cu, and V were readily activated and detected by the delayed method. H, Fe, Si, and Ni were not usually present in the delayed spectra but they responded well in the capture mode. It is shown that the borehole sonde can be configured to permit simultaneous delayed and capture spectra, permitting detection of all of these elements. Simulated borehole experiments indicate that, in the delayed gamma-ray mode, an infinite sample is achieved when the ore layer has a vertical thickness of ∼20 cm and a horizontal distance of '10 cm into the wallrock. The depth resolution is thus relatively good but horizontal penetration through the wallrock is limited. © 1972 Society of Economic Geologists, Inc.

Economic Geology↗

Effect of experimental technique on the determination of strontium distribution coefficients of a surficial sediment from the Idaho National Engineering Laboratory, Idaho

The effect of experimental technique on strontium distribution coefficients (K(d)'s) was determined as part of an investigation of strontium geochemical transport properties of surficial sediment from the Idaho National Engineering Laboratory, Idaho. The investigation was conducted by the U.S. Geological Survey and Idaho State University, in cooperation with the U.S. Department of Energy. Batch experiments were conducted to quantify the effect of different experimental techniques on experimentally derived strontium K(d)'s at a fixed pH of 8.0. Combinations of three variables were investigated: method of sample agitation (rotating-mixer and shaker table), ratio of the mass-of-sediment to the volume-of-reaction-solution (1:2 and 1:20), and method of sediment preparation (crushed and non-crushed). Strontium K(d)'s ranged from 11 to 23 mlg-1 among all three experimental variables examined. Strontium K(d)'s were bimodally grouped around 12 and 21 mlg-1. Among the three experimental variables examined, the mass-to-volume ratio appeared to be the only one that could account for this bimodal distribution. The bimodal distribution of the derived strontium K(d)'s may occur because the two different mass-to-volume ratios represent different natural systems. The high mass-to-volume ratio of 1:2 models a natural system, such as an aquifer, in which there is an abundance of favorable sorption sites relative to the amount of strontium in solution. The low mass-to-volume ratio of 1:20 models a natural system, such as a stream, in which the relative amount of strontium in solution exceeds the favorable surface sorption site concentration. Except for low mass-to-volume ratios of non-crushed sediment using a rotating mixer, the method of agitation and sediment preparation appears to have little influence on derived strontium K(d)'s.The effect of experimental technique on strontium distribution coefficients (Kd's) was determined as part of an investigation of strontium geochemical transport properties of surficial sediment from the Idaho National Engineering Laboratory, Idaho. The investigation was conducted by the U.S. Geological Survey and Idaho State University, in cooperation with the U.S. Department of Energy. Batch experiments were conducted to quantify the effect of different experimental techniques on experimentally derived strontium Kd's at a fixed pH of 8.0. Combinations of three variables were investigated: method of sample agitation (rotating-mixer and shaker table), ratio of the mass-of-sediment to the volume-of-reaction-solution (1:2 and 1:20), and method of sediment preparation (crushed and non-crushed). Strontium Kd's ranged from 11 to 23 mlg-1 among all three experimental variables examined. Strontium Kd's were bimodally grouped around 12 and 21 mlg-1. Among the three experimental variables examined, the mass-to-volume ratio appeared to be the only one that could account for this bimodal distribution. The bimodal distribution of the derived strontium Kd's may occur because the two different mass-to-volume ratios represent different natural systems. The high mass-to-volume ratio of 1:2 models a natural system, such as an aquifer, in which there is an abundance of favorable sorption sites relative to the amount of strontium in solution. The low mass-to-volume ratio of 1:20 models a natural system, such as a stream, in which the relative amount of strontium in solution exceeds the favorable surface sorption site concentration. Except for low mass-to-volume ratios of non-crushed sediment using a rotating mixer, the method of agitation and sediment preparation appears to have little influence on derived strontium Kd's.

Water Research↗

Integrated pest management and allocation of control efforts for vector-borne diseases

Applications of various control methods were evaluated to determine how to integrate methods so as to minimize the number of human cases of vector-borne diseases. These diseases can be controlled by lowering the number of vector-human contacts (e.g., by pesticide applications or use of repellents), or by lowering the proportion of vectors infected with pathogens (e.g., by lowering or vaccinating reservoir host populations). Control methods should be combined in such a way as to most efficiently lower the probability of human encounter with an infected vector. Simulations using a simple probabilistic model of pathogen transmission suggest that the most efficient way to integrate different control methods is to combine methods that have the same effect (e.g., combine treatments that lower the vector population; or combine treatments that lower pathogen prevalence in vectors). Combining techniques that have different effects (e.g., a technique that lowers vector populations with a technique that lowers pathogen prevalence in vectors) will be less efficient than combining two techniques that both lower vector populations or combining two techniques that both lower pathogen prevalence, costs being the same. Costs of alternative control methods generally differ, so the efficiency of various combinations at lowering human contact with infected vectors should be estimated at available funding levels. Data should be collected from initial trials to improve the effects of subsequent interventions on the number of human cases.

Journal of Vector Ecology↗

Discriminating silt-and-clay from suspended-sand in rivers using side-looking acoustic profilers

The ability to accurately monitor suspended-sediment flux in rivers is needed to support many types of studies, because the sediment that typically travels in suspension affects geomorphology and aquatic habitat in a variety of ways (e.g. bank and floodplain deposition, bar morphology, light penetration and primary productivity, tidal wetland deposition in the context of sea-level rise, sediment-associated contaminants, reservoir sedimentation and potential erosion during dam removal, among others). In addition, human-induced changes to the landscape have resulted in substantially altered suspended-sediment loads (Syvitski et al., 2005). Thus, accurate monitoring of suspended-sediment flux is necessary for informed resource management of rivers. Because of this need, a variety of techniques have been developed and applied for suspendedsediment monitoring. The traditional approach in the United States, which was developed and has been used extensively by the U.S. Geological Survey (USGS), is to collect an isokinetic, velocity-weighted sample from a river cross-section, analyze the sample in the laboratory, and use water-discharge records to compute a record of suspended-sediment flux (Guy, 1969, Guy, 1970, Edwards and Glysson, 1999, Porterfield, 1972). The labor and expense associated with this traditional approach is substantial such that the number of USGS gages reporting daily records of suspended-sediment flux decreased from 364 in 1981 to 120 in 2003 (Osterkamp et al., 2004). Also, the traditional sampling approach is limited with respect to the temporal resolution that can be achieved, thus requiring the use of approximate relations between suspended-sediment concentration and water discharge to fill gaps between samples. To address these limitations, several indirect or "surrogate" measures have been investigated (see e.g. Gray and Gartner, 2009) most notably optical backscatter (i.e. turbidity), laser-diffraction, and acoustic backscatter. These indirect techniques rely on measurements of ancillary properties that correlate with suspended-sediment concentration and particle size and thus require the collection of traditional samples for calibration. Through in situ deployments, these methods can provide the high temporal resolution that cannot be achieved through traditional sampling. Here we focus on the evaluation of acoustic profiling techniques (e.g. acoustic-Doppler sideways-looking profilers, or ADPs). One major advantage of acoustic profiling is the ability to concurrently measure water velocity (using Doppler-shift methods) and suspended-sediment concentration such that suspended-sediment flux can be directly computed using data from a single instrument. Acoustic-Doppler profilers have become popular for measuring water velocity and discharge in rivers, through both moving-boat operations and from fixed deployments such as bank-mounted sideways-looking instruments (Hirsch and Costa, 2004, Muste et al., 2007). The method presented herein is most suited to sideways-looking applications as a complement to the "index velocity" technique, whereby an index velocity from a sideways-looking instrument is related to the cross-section average velocity (determined from moving-boat discharge measurements) as a means for developing a continuous water-discharge record (Ruhl and Simpson, 2005). Topping et al. (2007) presented a method for discriminating silt-and-clay from suspended sand, using single frequency ADPs. This method takes advantage of the relations among acoustic backscatter, sediment-induced acoustic attenuation, suspended-sediment concentration (SSC), and particle size distribution (PSD). Backscatter is the amount of sound scattered back and received at the transducer while sediment-induced attenuation is the amount of sound scattered in other directions and absorbed by the sediment particles. Both of these parameters can be measured with an ADP, and their different dependencies on SSC and PSD allow for the discrimination of suspended silt-and-clay from suspended sand. Topping et al. (2007) describe application of the method at several sites along the Colorado River in Grand Canyon, and herein we present an example application of the technique for the Gunnison River, CO. However, the methods general applicability in rivers has yet to be evaluated due to a lack of concurrent acoustic and sediment data at a range of sites. To this end, the objective of the analysis presented herein is to evaluate the potential general applicability of the method, drawing from the extensive USGS database on SSC and PSD. We refer to it as "potential" general applicability because it relies on the theory underlying the previous empirical results. Use of the theoretical relations is necessary due to the lack of concurrent ADP and SSC/PSD data, but also serves the additional purpose of providing further justification of the empirical calibrations developed for the Colorado and Gunnison Rivers.

Conference Paper↗

An autocorrelation method to detect low frequency earthquakes within tremor

Recent studies have shown that deep tremor in the Nankai Trough under western Shikoku consists of a swarm of low frequency earthquakes (LFEs) that occur as slow shear slip on the down-dip extension of the primary seismogenic zone of the plate interface. The similarity of tremor in other locations suggests a similar mechanism, but the absence of cataloged low frequency earthquakes prevents a similar analysis. In this study, we develop a method for identifying LFEs within tremor. The method employs a matched-filter algorithm, similar to the technique used to infer that tremor in parts of Shikoku is comprised of LFEs; however, in this case we do not assume the origin times or locations of any LFEs a priori. We search for LFEs using the running autocorrelation of tremor waveforms for 6 Hi-Net stations in the vicinity of the tremor source. Time lags showing strong similarity in the autocorrelation represent either repeats, or near repeats, of LFEs within the tremor. We test the method on an hour of Hi-Net recordings of tremor and demonstrates that it extracts both known and previously unidentified LFEs. Once identified, we cross correlate waveforms to measure relative arrival times and locate the LFEs. The results are able to explain most of the tremor as a swarm of LFEs and the locations of newly identified events appear to fill a gap in the spatial distribution of known LFEs. This method should allow us to extend the analysis of Shelly et al. (2007a) to parts of the Nankai Trough in Shikoku that have sparse LFE coverage, and may also allow us to extend our analysis to other regions that experience deep tremor, but where LFEs have not yet been identified. Copyright 2008 by the American Geophysical Union.

Geophysical Research Letters↗

Red-shouldered hawk occupancy surveys in central Minnesota, USA

Forest-dwelling raptors are often difficult to detect because many species occur at low density or are secretive. Broadcasting conspecific vocalizations can increase the probability of detecting forest-dwelling raptors and has been shown to be an effective method for locating raptors and assessing their relative abundance. Recent advances in statistical techniques based on presence-absence data use probabilistic arguments to derive probability of detection when it is <1 and to provide a model and likelihood-based method for estimating proportion of sites occupied. We used these maximum-likelihood models with data from red-shouldered hawk (Buteo lineatus) call-broadcast surveys conducted in central Minnesota, USA, in 1994-1995 and 2004-2005. Our objectives were to obtain estimates of occupancy and detection probability 1) over multiple sampling seasons (yr), 2) incorporating within-season time-specific detection probabilities, 3) with call type and breeding stage included as covariates in models of probability of detection, and 4) with different sampling strategies. We visited individual survey locations 2-9 times per year, and estimates of both probability of detection (range = 0.28-0.54) and site occupancy (range = 0.81-0.97) varied among years. Detection probability was affected by inclusion of a within-season time-specific covariate, call type, and breeding stage. In 2004 and 2005 we used survey results to assess the effect that number of sample locations, double sampling, and discontinued sampling had on parameter estimates. We found that estimates of probability of detection and proportion of sites occupied were similar across different sampling strategies, and we suggest ways to reduce sampling effort in a monitoring program.

Journal of Wildlife Management↗

The Uranium-trend dating method: Principles and application for southern California marine terrace deposits

Uranium-trend dating is an open-system method for age estimation of Quaternary sediments, using disequilibrium in the 238U234U230Th decay series. The technique has been applied to alluvium, colluvium, loess, till, and marine sediments. In this study we tested the U-trend dating method on calcareous marine terrace deposits from the Palos Verdes Hills and San Nicolas Island, California. Independent age estimates indicate that terraces in these areas range from ???80 ka to greater than 1.0 Ma. Two low terraces on San Nicolas Island yielded U-trend plots that have a clustered array of points and the ages of these deposits are indeterminate or highly suspect. Middle Pleistocene terraces and one early Pleistocene terrace on San Nicolas Island and all terraces on the Palos Verdes Hills gave reasonably linear U-trend plots and estimated ages that are stratigraphically consistent and in agreement with independent age estimates. We conclude that many marine terrace deposits are suitable for U-trend dating, but U-trend plots must be carefully evaluated and U-trend ages should be consistent with independent geologic control. ?? 1989.

Quaternary International↗

Magnitude and frequency of floods in North Carolina, technique for estimating the magnitude and frequency of floods on natural streams in North Carolina

Methods are provided to estimate the magnitude and frequency of floods on natural North Carolina streams with drainage areas greater than 0.5 square mile (1.3 square kilometers). For 257 gaged sites, the magnitudes of floods having recurrence intervals from 2 to 100 years are provided in tables. For ungaged sites, equations, graphs, and maps allow estimation of flood magnitudes. Multiple regression techniques were used to define the relation between flood peaks and seven basin and climatic variables. Drainage area is the most significant variable. Inclusion of the other six variables reduced the standard error of estimate less than 4 percent. Regression equations gave consistently different results for stations in the Coastal Plain than for stations in the mountains and Piedmont. Accordingly, stations were divided into two groups and estimating equations were developed for each geographic area.

North Carolina↗

Improved efficiency of maximum likelihood analysis of time series with temporally correlated errors

Most time series of geophysical phenomena have temporally correlated errors. From these measurements, various parameters are estimated. For instance, from geodetic measurements of positions, the rates and changes in rates are often estimated and are used to model tectonic processes. Along with the estimates of the size of the parameters, the error in these parameters needs to be assessed. If temporal correlations are not taken into account, or each observation is assumed to be independent, it is likely that any estimate of the error of these parameters will be too low and the estimated value of the parameter will be biased. Inclusion of better estimates of uncertainties is limited by several factors, including selection of the correct model for the background noise and the computational requirements to estimate the parameters of the selected noise model for cases where there are numerous observations. Here, I address the second problem of computational efficiency using maximum likelihood estimates (MLE). Most geophysical time series have background noise processes that can be represented as a combination of white and power-law noise, 1 / f &#x03B1; "> 1 / f α 1/fα with frequency, f . With missing data, standard spectral techniques involving FFTs are not appropriate. Instead, time domain techniques involving construction and inversion of large data covariance matrices are employed. Bos et al. (J Geod, 2013 . doi: 10.1007/s00190-012-0605-0 ) demonstrate one technique that substantially increases the efficiency of the MLE methods, yet is only an approximate solution for power-law indices >1.0 since they require the data covariance matrix to be Toeplitz. That restriction can be removed by simply forming a data filter that adds noise processes rather than combining them in quadrature. Consequently, the inversion of the data covariance matrix is simplified yet provides robust results for a wider range of power-law indices.

Journal of Geodesy↗

Use of stable sulphur isotopes to monitor directly the behaviour of sulphur in coal during thermal desulphurization

A method has been developed using stable sulphur isotope analyses to monitor the behaviour of sulphur forms in a coal during thermal desulphurization. In this method, the natural stable isotopic composition of the pyritic and organic sulphur in coal is used as a tracer to follow their mobility during the desulphurization process. This tracer method is based on the fact that the isotopic compositions of pyritic and organic sulphur are significantly different in some coals. Isotopic results of pyrolysis experiments at temperatures ranging from 350 to 750 °C indicate that the sulphur released with the volatiles is predominantly organic sulphur. The pyritic sulphur is evolved in significant quantities only when pyrolysis temperatures exceed 500 °C. The presence of pyrite seems to have no effect on the amount of organic sulphur evolved during pyrolysis. The chemical and isotopic mass balances achieved from three different samples of the Herrin (No. 6) coal of the Illinois Basin demonstrate that this stable isotope tracer method is quantitative. The main disadvantage of this tracing technique is that not all coals contain isotopically distinct organic and pyritic sulphur.

Illinois, Indiana, Kentucky↗

Regional regression of flood characteristics employing historical information

Streamflow gauging networks provide hydrologic information for use in estimating the parameters of regional regression models. The regional regression models can be used to estimate flood statistics, such as the 100 yr peak, at ungauged sites as functions of drainage basin characteristics. A recent innovation in regional regression is the use of a generalized least squares (GLS) estimator that accounts for unequal station record lengths and sample cross correlation among the flows. However, this technique does not account for historical flood information. A method is proposed here to adjust this generalized least squares estimator to account for possible information about historical floods available at some stations in a region. The historical information is assumed to be in the form of observations of all peaks above a threshold during a long period outside the systematic record period. A Monte Carlo simulation experiment was performed to compare the GLS estimator adjusted for historical floods with the unadjusted GLS estimator and the ordinary least squares estimator. Results indicate that using the GLS estimator adjusted for historical information significantly improves the regression model.

Journal of Hydrology↗