Geology Reports⌕ Search

SEARCH · Geology Reports

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 199 records · Page 11Linked to original sources

GEOGRAPHIC ESTIMATION OF RUNOFF-MODEL PARAMETERS.

The U. S. Geological Survey is developing techniques to estimate and evaluate unit-hydrograph and loss-rate parameter values for rainfall-runoff models using Geographic Information System (GIS) procedures. The data base includes basin, soil, and climatological characteristics that will be stored in a GIS, and unit-hydrograph and loss-rate parameters obtained from calibration of a commonly used flood-hydrograph rainfall-runoff model for 616 storms in 98 gaged drainage basins. Development of unit-hydrograph and loss-rate parameter-estimation techniques includes statistical methods (exploratory data analysis, regression analysis, and categorical data analysis) to relate the model parameters to hydrologic characteristics. The estimation techniques are evaluated by use of error analysis of simulated hydrograph characteristics (peak discharge, flood volume, and time to peak discharge). The hydrographs will be simulated with parameters estimated by the techniques for (1) 102 storms occurring at 36 gaged basins; and (2) a large storm system (one which produced floods with a 50-to 100-year recurrence interval).

Conference Paper↗

KG²B, a collaborative benchmarking exercise for estimating the permeability of the Grimsel granodiorite - Part 2: modeling, microstructures and complementary data

Measuring and modelling the permeability of tight rocks remains a challenging task. In addition to the traditional sources of errors that affect more permeable formations (e.g. sample selection, non-representative specimens, disturbance introduced during sample acquisition and preparation), tight rocks can be particularly prone to solid–fluid interactions and thus more sensitive to the methods, procedures and techniques used to measure permeability. To address this problem, it is desirable to collect, for a single material, measurements obtained by different methods and pore fluids. For that purpose, a benchmarking exercise involving 24 laboratories was organized for measuring and modelling the permeability of a single low-permeability material, the Grimsel granodiorite. The objectives of the benchmark were: (i) to compare the results for a given method, (ii) to compare the results between different methods, (iii) to analyse the accuracy of each method, (iv) to study the influence of experimental conditions (especially the nature of pore fluid), (v) to discuss the relevance of indirect methods and models and finally (vi) to suggest good practice for low-permeability measurements. To complement the data set of permeability measurements presented in a companion paper, we focus here on (i) quantitative analysis of microstructures and pore size distribution, (ii) permeability modelling and (iii) complementary measurements of permeability anisotropy and poroelastic parameters. Broad ion beam—scanning electron microscopy, micro-computerized tomography, mercury injection capillary pressure (MICP) and nuclear magnetic resonance (NMR) methods were used to characterize the microstructures and provided the input parameters for permeability modelling. Several models were used: (i) basic statistical models, (ii) 3-D pore network and effective medium models, (iii) percolation model using MICP data and (iv) free-fluid model using NMR data. The models were generally successful in predicting the actual range of measured permeability. Statistical models overestimate the permeability because they do not adequately account for the heterogeneity of the crack network. Pore network and effective medium models provide additional constraints on crack parameters such as aspect ratio, aperture, density and connectivity. MICP and advanced microscopy techniques are very useful tools providing important input data for permeability estimation. Permeability measured—orthogonal to foliation is lower that—parallel to foliation. Combining the experimental and modelling results provide a unique and rich data set.

Geophysical Journal International↗

KG²B, a collaborative benchmarking exercise for estimating the permeability of the Grimsel granodiorite - Part 1: measurements, pressure dependence and pore-fluid effects

Measuring the permeability of tight rocks remains a challenging task. In addition to the traditional sources of errors that affect more permeable formations (e.g. sample selection, non-representative specimens, disturbance introduced during sample acquisition and preparation), tight rocks can be particularly prone to solid–fluid interactions and thus more sensitive to the methods, procedures and techniques used to measure permeability. To address this problem, it is desirable to collect, for a single material, measurements obtained by different methods and pore-fluids. For that purpose a collaborative benchmarking exercise involving 24 laboratories was organized for measuring the permeability of a single low permeability material, the Grimsel granodiorite, at a common effective confining pressure (5 MPa). The objectives of the benchmark were: (i) to compare the results for a given method, (ii) to compare the results between different methods, (iii) to analyze the accuracy of each method, (iv) to study the influence of experimental conditions (especially the nature of pore fluid), (v) to discuss the relevance of indirect methods and models and finally (vi) to suggest good practice for low permeability measurements. In total 39 measurements were collected that allowed us to discuss the influence of (i) pore-fluid, (ii) measurement method, (iii) sample size and (iv) pressure sensitivity. Discarding some outliers from the bulk data set (4 out of 39) an average permeability of 1.11 × 10 −18 m² with a standard deviation of 0.57 × 10 −18 m² was obtained. The most striking result was the large difference in permeability for gas measurements compared to liquid measurements. Regardless of the method used, gas permeability was higher than liquid permeability by a factor approximately 2 ( k gas = 1.28 × 10 −18 m² compared to k liquid = 0.65 × 10 −18 m²). Possible explanations are that (i) liquid permeability was underestimated due to fluid-rock interactions (ii) gas permeability was overestimated due to insufficient correction for gas slippage and/or (iii) gases and liquids do not probe exactly the same porous networks. The analysis of Knudsen numbers shows that the gas permeability measurements were performed in conditions for which the Klinkenberg correction is sufficient. Smaller samples had a larger scatter of permeability values, suggesting that their volume were below the Representative Elementary Volume. The pressure dependence of permeability was studied by some of the participating teams in the range 1–30 MPa and could be fitted to an exponential law k = k o .exp(– γP eff ) with γ = 0.093 MPa −1 . Good practice rules for measuring permeability in tight materials are also provided.

Geophysical Journal International↗

Techniques for estimating the magnitude and frequency of floods in rural basins of South Carolina, 1999

Data from 167 streamflow-gaging stations in or near South Carolina with 10 or more years of record through September 30, 1999, were used to develop two methods for estimating the magnitude and frequency of floods in South Carolina for rural ungaged basins that are not significantly affected by regulation. Floodfrequency estimates for 54 gaged sites in South Carolina were computed by fitting the water-year peak flows for each site to a log-Pearson Type III distribution. As part of the computation of flood-frequency estimates for gaged sites, new values for generalized skew coefficients were developed. Flood-frequency analyses also were made for gaging stations that drain basins from more than one physiographic province. The U.S. Geological Survey, in cooperation with the South Carolina Department of Transportation, updated these data from previous flood-frequency reports to aid officials who are active in floodplain management as well as those who design bridges, culverts, and levees, or other structures near streams where flooding is likely to occur. Regional regression analysis, using generalized least squares regression, was used to develop a set of predictive equations that can be used to estimate the 2-, 5-, 10-, 25-, 50-, 100-, 200-, and 500-year recurrence- interval flows for rural ungaged basins in the Blue Ridge, Piedmont, upper Coastal Plain, and lower Coastal Plain physiographic provinces of South Carolina. The predictive equations are all functions of drainage area. Average errors of prediction for these regression equations ranged from -16 to 19 percent for the 2-year recurrence-interval flow in the upper Coastal Plain to -34 to 52 percent for the 500-year recurrenceinterval flow in the lower Coastal Plain. A region-of-influence method also was developed that interactively estimates recurrence- interval flows for rural ungaged basins in the Blue Ridge of South Carolina. The region-of-influence method uses regression techniques to develop a unique relation between flow and basin characteristics for an individual watershed. This, then, can be used to estimate flows at ungaged sites. Because the computations required for this method are somewhat complex, a computer application was developed that performs the computations and compares the predictive errors for this method. The computer application includes the option of using the region-of-influence method, or the generalized least squares regression equations from this report to compute estimated flows and errors of prediction specific to each ungaged site. From a comparison of predictive errors using the region-of-influence method with those computed using the regional regression method, the region-of-influence method performed systematically better only in the Blue Ridge and is, therefore, not recommended for use in the other physiographic provinces. Peak-flow data for the South Carolina stations used in the regionalization study are provided in appendix A, which contains gaging station information, log- Pearson Type III statistics, information on stage-flow relations, and water-year peak stages and flows. For informational purposes, water-year peak-flow data for stations on regulated streams in South Carolina also are provided in appendix D. Other information pertaining to the regulated streams is provided in the text of the report.

South Carolina↗

Aging small Canada geese by neck plumage

The neck plumage method, a new technique for separating immature from adult Canada geese ( Branta canadensis ) in the hand, was evaluated by comparison with the notched tail feather and cloacal examination methods. Two (1.4 percent) of 141 geese examined were misaged, resulting in a 6 percent error in the immature-adult ratio obtained by the neck plumage method. The neck plumage method is a rapid aging method and reasonable accuracy (94 percent) can be obtained. It can also be used to differentiate immatures from adults on the ground at distances up to 175 yards, but was almost impossible to use when geese were in flight. As yet, the neck plumage method has only been tested on the subspecies ( B. c. hutchinsii-parvipes complex) in the Tall-Grass Prairie population of small Canada geese.

South Dakota↗

Sediment accumulation and water volume in Loch Raven Reservoir, Baltimore County, Maryland

Baltimore City and its metropolitan area are supplied with water from three reservoirs, Liberty Reservoir, Prettyboy Reservoir, and Loch Raven Reservoir. Prettyboy and Loch Raven Reservoirs are located on the Gunpowder Falls (figure 1). The many uses of the reservoir system necessitate coordination and communication among resource managers. The 1996 Amendment to the Safe Drinking Water Act require States to complete source-water assessments for public drinking-water supplies. As part of an ongoing effort to provide safe drinking water and as a direct result of these laws, the City of Baltimore and the Maryland Department of the Environment (MDE), in cooperation with other State and local agencies, are studying the Gunpowder Falls Basin and its role as a source of water supply to the Baltimore area. As a part of this study, the U.S. Geological Survey (USGS), in cooperation with the Maryland Geological Survey (MGS), with funding provided by the City of Baltimore and MDE, is examining sediment accumulation in Loch Raven Reservoir. The Baltimore City Department of Public Works periodically determines the amount of water that can be stored in its reservoirs. To make this determination, field crews measure the water depth along predetermined transects or ranges. These transects provide consistent locations where water depth, or bathymetric, measurements can be made. Range surveys are repeated to provide a record of the change in storage capacity due to sediment accumulation over time. Previous bathymetric surveys of Loch Raven Reservoir were performed in 1943, 1961, 1972, and 1985. Errors in data-collection and analysis methods have been assessed and documented (Baltimore City Department of Public Works, 1989). Few comparisons can be made among survey results because of changing data-collection techniques and analysis methods.

Water-Resources Investigations Report↗

Choosing a DIVA: a comparison of emerging digital imagery vegetation analysis techniques

Question: What is the precision of five methods of measuring vegetation structure using ground-based digital imagery and processing techniques? Location: Lincoln, Nebraska, USA Methods: Vertical herbaceous cover was recorded using digital imagery techniques at two distinct locations in a mixed-grass prairie. The precision of five ground-based digital imagery vegetation analysis (DIVA) methods for measuring vegetation structure was tested using a split-split plot analysis of covariance. Variability within each DIVA technique was estimated using coefficient of variation of mean percentage cover. Results: Vertical herbaceous cover estimates differed among DIVA techniques. Additionally, environmental conditions affected the vertical vegetation obstruction estimates for certain digital imagery methods, while other techniques were more adept at handling various conditions. Overall, percentage vegetation cover values differed among techniques, but the precision of four of the five techniques was consistently high. Conclusions: DIVA procedures are sufficient for measuring various heights and densities of standing herbaceous cover. Moreover, digital imagery techniques can reduce measurement error associated with multiple observers' standing herbaceous cover estimates, allowing greater opportunity to detect patterns associated with vegetation structure.

Nebraska↗

Time-to-detection occupancy modeling: An efficient method for analyzing the occurrence of amphibians and reptiles

Occupancy models provide a reliable method of estimating species distributions while accounting for imperfect detectability. The cost of accounting for false absences is that detection and nondetection surveys typically require repeated visits to a site or multiple-observer techniques. More efficient methods of collecting data to estimate detection probabilities would allow additional sites to be surveyed for the same amount of effort, which would support more precise estimation of covariate effects to improve inference about underlying ecological processes. Time-to-detection surveys allow the estimation of detection probability based on a single site visit by one observer, and therefore might be an efficient technique for herpetological occupancy studies. We evaluated the use of time-to-detection surveys to estimate the occupancy of pond-breeding amphibians at Point Reyes National Seashore, California, USA, including variables that affected detection rates and the probability of occurrence. We found that detection times were short enough, and occupancy was high enough, to estimate reliably the probability of occurrence of three pond-breeding amphibians at Point Reyes National Seashore, and that survey and site conditions had species-specific effects on detection rates. In particular, pond characteristics affected detection times of all commonly detected species. Probability of occurrence of Sierran Treefrogs ( Hyliola sierra ) and Rough-Skinned Newts ( Taricha granulosa ) was negatively related to the detection of fish and pond area. Time-to-detection surveys can provide an efficient method for estimating detection probabilities and accounting for false absences in occupancy studies of reptiles and amphibians.

California↗

Seismic Station Installation Orientation Errors at ANSS and IRIS/USGS Stations

Many seismological studies depend on the published orientations of sensitive axes of seismic instruments relative to north (e.g., Li et al., 2011). For example, studies of the anisotropic structure of the Earth’s mantle through SKS‐splitting measurements (Long et al., 2009), constraints on core–mantle electromagnetic coupling from torsional normal‐mode measurements (Dumberry and Mound, 2008), and models of three‐dimensional (3D) velocity variations from surface waves (Ekström et al., 1997) rely on accurate sensor orientation. Unfortunately, numerous results indicate that this critical parameter is often subject to significant error (Laske, 1995; Laske and Masters, 1996; Yoshizawa et al., 1999; Schulte‐Pelkum et al., 2001; Larson and Ekström, 2002). For the Advanced National Seismic System (ANSS; ANSS Technical Integration Committee, 2002), the Global Seismographic Network (GSN; Butler et al., 2004), and many other networks, sensor orientation is typically determined by a field engineer during installation. Successful emplacement of a seismic instrument requires identifying true north, transferring a reference line, and measuring the orientation of the instrument relative to the reference line. Such an exercise is simple in theory, but there are many complications in practice. There are four commonly used methods for determining true north at the ANSS and GSN stations operated by the USGS Albuquerque Seismological Laboratory (ASL), including gyroscopic, astronomical, Global Positioning System (GPS), and magnetic field techniques. A particular method is selected based on site conditions (above ground, below ground, availability of astronomical observations, and so on) and in the case of gyroscopic methods, export restrictions. Once a north line has been determined, it must be translated to the sensor location. For installations in mines or deep vaults, this step can include tracking angles through the one or more turns in the access tunnel leading to the vault (e.g., GSN station WCI in Wyandotte Cave, Indiana). Finally, the third source of error comes from the ability of field engineers to orient the sensor relative to the reference line. In order to quantify bounds on the errors in each step in the orientation process, we conducted a series of tests at the ASL using twelve GSN and ANSS field engineers. The results from this exercise allow us to estimate upper bounds on the precision of our ability to orient instruments, as well as identify the sources of error in the procedures. We are also able to identify systematic bias of various true‐north‐finding methods relative to one another. Although we are unable to estimate the absolute accuracy of our orientation measurements due to our inability to identify true north without some error, the agreement between independent methods for finding true north provides confidence in the different approaches, assuming no systematic bias. Finally, our study neglects orientation errors that are beyond the control of the field engineer during a station visit. These additional errors can arise from deviations in the sensitive axes of the instruments relative to the case markings, processing errors (Holcomb, 2002) when comparing horizontal orientations relative to other sensors (e.g., borehole installations), and deviations of the sensitive axes of instruments from true orthogonality (e.g., instruments with separate modules such as the Streckeisen STS‐1).

Seismological Research Letters↗

Relocation of Wyoming mine production blasts using calibration explosions

An important requirement for a comprehensive seismic monitoring system is the capability to accurately locate small seismic events worldwide. Accurate event location can improve the probability of determining whether or not a small event, recorded predominantly by local and regional stations, is a nuclear explosion. For those portions of the earth where crustal velocities are not well established, reference event calibration techniques offer a method of increased locational accuracy and reduced locational bias. In this study, data from a set of mining events with good ground-truth data in the Powder River Basin region of eastern Wyoming are used to investigate the potential of event calibration techniques in the area. Results of this study are compared with locations published in the prototype International Data Center’s Reviewed Event Bulletin (REB). A Joint Hypocenter Determination (JHD) method was applied to a set of 23 events. Four of those events with superior ground-truth control (mining company report or Global Positioning System data) were used as JHD reference events, Nineteen (83%) of the solutions converged and the resulting set of station-phase travel-time corrections from the JHD results was then tested. When those travel-time corrections were applied individually to the four events with good ground-truth control, the average locational error reduced the original REB location error from 16.1 km to 5.7 km (65% improvement). The JHD locations indicated reduced locational bias and all of the individual error ellipses enclosed the actual known event locations. Given a set of well-recorded calibration events, it appears that the JHD methodology is a viable technique for improving locational accuracy of future small events where the location depends on arrival times from predominantly local and/or regional stations. In this specific case, the International Association of Seismology and the Physics of the Earth’s Interior (IASPEI) travel-time tables, coupled with JHDderived travel-time corrections, may obviate the need for an accurately known regional velocity structure in the Powder River Basin region.

Montana, Wyoming↗

Multi-scale and nested-intensity sampling techniques for archaeological survey

This paper discusses sampling techniques for archaeological survey that are directed toward evaluating the properties of surface artifact distributions. The sampling techniques we experimented with consist of a multi-scale sampling plot developed in plant ecology and the use of a nested-intensity survey design. We present results from the initial application of these methods. The sampling technique we borrowed from plant ecology is the Modified-Whittaker multiscale sampling plot, which gathers observations at the spatial scales of 1 sq m, 10 sq m, 100 sq m, and 1000 sq m. Nested-intensity surveys gather observations on the same sample units at multiple resolutions. We compare the results of a closely-spaced walking survey, a crawling survey, and a test excavation to a depth of 10 cm. These techniques were applied to ten 20 × 50 m survey plots distributed over an area of 418 ha near the Hudson-Meng Bison Bonebed in NW Nebraska. These approaches can significantly improve the accuracy of survey data. Our results show that high-resolution coverage techniques overlook more material than archaeologists have suspected. The combined approaches of multi-scale and nested-intensity sampling provide new tools to improve our ability to investigate the properties of surface records.

Journal of Field Archaeology↗

Uranium values from ion-filter water sampling compared with values from bulk-water sampling

Analytical values of uranium from well water sampled by means of ion-filters were compared with analytical values of the same water sampled by the standard bulk (1-liter) water-sampling technique (Brown, Skougstad and Fishman, 1970). A fluorometric extraction method of analysis was used on both sets of samples. Duplicate samples from pumped wells in Pueblo County, Colorado were used. Regression and correlation statistical analysis indicated that while analytical results of the two sets of samples had a highly significant correlation (r = 0.98), the ion-filter samples yielded only 20 percent as much uranium as the 1-liter water samples. It is possible to predict U values from the ion-filter analysis results that would be comparable with the values for the 1-liter samples, using the linear regression formula. The predicted values, however, compared with the 1-liter values, using a chi-square test, were found to be significantly different at the 5-percent level. Major deviations were found in two pairs, 10 percent of the samples, suggesting that the predicted values might be used if one were willing to accept the fact that an average of 10 percent of the samples would yield erroneous results. Sources of the absolute deviations between the two sets of samples may lie in chemical characteristics of the filters, or waters, or extraction method of the analytical technique or some combination of these. The results of the chi-square test also indicate that uranium values from at least 10 percent of the ion-filter samples that we obtained, deviate significantly from those obtained from the 1-liter samples. Sources of the deviations may lie in chemical characteristics of the filters, or waters, or extraction technique, part of the analytical technique or some combination of these.

Open-File Report↗

Measurement of in vitro leucocyte mitogenesis in fish: ELISA based detection of the thymidine analogue 5'-bromo-2'-deoxyuridine

In this study we present a method for the measurement of in vitro mitogenesis in fish leucocytes that is based on the incorporation of the thymidine analogue 5′-bromo-2′-deoxyuridine (BrdU) into the DNA of replicating cells, followed by ELISA-based detection. This technique, adapted from methods developed for mammalian cells, operates on a similar biological principle to 3 H-thymidine incorporation, but circumvents the logistical and safety issues inherent with the radioactive label. Because it directly measures DNA proliferation, the assay has advantages over other colorimetric methods that may be strongly influenced by leucocyte metabolic status. Using BrdU incorporation followed by ELISA, we evaluate the responsiveness of rainbow trout ( Oncorhynchus mykiss [Walbaum]) leucocytes to the mammalian T-cell mitogen Concanavalin A (Con A) as well as the differential response of white perch ( Morone americana [Gmelin]) leucocytes to Con A and pokeweed mitogen. Specific considerations intrinsic to the assay system are discussed, including the implications of utilising enzyme-based detection.

Fish and Shellfish Immunology↗

In my experience: Improved capture techniques for psittacines

Four methods for capturing psittacines were developed and tested in Puerto Rico from 1991-1993. Elevated mist nets at canopy height in scrub-forest or mangrove habitat captured fewer nontarget birds and possibly more parakeets than mist nets at ground level. Playback of conspecific calls may have attracted parakeets to the net area and reduced the need to erect more than six 12-m long nets. Using 6 nets instead of 12 also reduced the capture rate and captures of nontarget birds. This reduced disturbance at the netting site.....Mist nets, 7.3 m high, surrounding a tree with live decoys within the circle of nets captured parrots effectively when conspecific vocalizations of foraging parrots were played. More parrot captures occurred when the birds were frightened as they approached the decoy inside and below the circle of nets during the early morning. Nets at lower heights (5.8 m) and nets centered perpendicularly on a tree with picture and live decoys failed to capture parrots.

Wildlife Society Bulletin↗

A practical, low-noise coil system for magnetotellurics

Magnetotellurics is a geophysical technique which was developed by Cagnaird (1953) and Tikhonov (1950) and later refined by other scientists worldwide. The technique is a method of electromagnetic sounding of the Earth and is based upon the skin depth effect in conductive media. The electric and magnetic fields arising from natural sources are measured at the surface of the earth over broad frequency bands. An excellent review of the technique is provided in the paper by Vozoff (1972). The sources of the natural fields are found in two basic mechanisms. At frequencies above a few hertz, most of the energy arises from lightning in thunderstorm belts around the equatorial regions. This energy is propagated in a wave-guide formed by the earthionospheric cavity. Energy levels are higher at fundamental modes for this cavity, but sufficient energy exists over most of the audio range to be useful for sounding at these frequencies, in which case the technique is generally referred to as audio-magnetotellurics or AMT. At frequencies lower than audio, and in general below 1 Hz, the source of naturally occuring electromagnetic energy is found in ionospheric currents. Current systems flowing in the ionosphere generate EM waves which can be used in sounding of the earth. These fields generate a relatively complete spectrum of electromagnetic energy that extends from around 1 Hz to periods of one day. Figure 1 shows an amplitude spectrum characteristic of both the ionospheric and lightning sources, covering a frequency range from 0.0001 Hz to 1000 Hz. It can be seen that there is a minimum in signal levels that occurs at about 1 Hz, in the gap between the two sources, and that signal level increases with a decrease in frequency.

Open-File Report↗

A Comparison of seismic instrument noise coherence analysis techniques

The self-noise of a seismic instrument is a fundamental characteristic used to evaluate the quality of the instrument. It is important to be able to measure this self-noise robustly, to understand how differences among test configurations affect the tests, and to understand how different processing techniques and isolation methods (from nonseismic sources) can contribute to differences in results. We compare two popular coherence methods used for calculating incoherent noise, which is widely used as an estimate of instrument self-noise (incoherent noise and self-noise are not strictly identical but in observatory practice are approximately equivalent; Holcomb, 1989; Sleeman et al., 2006). Beyond directly comparing these two coherence methods on similar models of seismometers, we compare how small changes in test conditions can contribute to incoherent-noise estimates. These conditions include timing errors, signal-to-noise ratio changes (ratios between background noise and instrument incoherent noise), relative sensor locations, misalignment errors, processing techniques, and different configurations of sensor types.

Bulletin of the Seismological Society of America↗

Genetic diversity of Wolbachia endosymbionts in Culex quinquefasciatus from Hawai`i, Midway Atoll, and Samoa

Incompatible insect techniques are potential methods for controlling Culex quinquefasciatus and avian disease transmission in Hawai‘i without the use of pesticides or genetically modified organisms. The approach is based on naturally occurring sperm-egg incompatibilities within the Culex pipiens complex that are controlled by different strains of the bacterial endosymbiont Wolbachia pipientis (wPip). Incompatibilities can be unidirectional (crosses between males infected with strain A and females infected with strain B are fertile, while reciprocal crosses are not) or bidirectional (reciprocal crosses between sexes with different wPip strains are infertile). The technique depends on release of sufficient numbers of male mosquitoes infected with an incompatible wPip strain to suppress mosquito populations and reduce transmission of introduced avian malaria ( Plasmodium relictum ) and Avipoxvirus in native forest bird habitats. Both diseases are difficult to manage using more traditional methods based on removal and treatment of larval habitats and coordination of multiple approaches may be needed to control this vector. We characterized the diversity of Wolbachia strains in C. quinquefasciatus from Hawai‘i, Kaua‘i, Midway Atoll, and American Samoa with a variety of genetic markers to identify compatibility groups and their distribution within and between islands. We confirmed the presence of wPip with multilocus sequence typing, tested for local genetic variability using 16 WO prophage genes, and identified similarities to strains from other parts of the world with a transposable element (tr1). We also tested for genetic differences in ankyrin motifs (ank2 and pk1) which have been used to classify wPip strains into five worldwide groups (wPip1–wPip5) that vary in compatibility with each other based on experimental crosses. We found a mixture of both widely distributed and site specific genotypes based on presence or absence of WO prophage and transposable element markers on Hawai‘i Island (Volcano, Pu‘u Wa‘awa‘a, Laupāhoehoe, Kaumana, Kahuku, Nīnole, and Maulua Gulch), Kaua‘i Island (Kawaikōī, Mōhihi, Kalāheo, Lāwa‘i and Hanapepe) and Midway Atoll. Genotypes from American Samoa were unique and formed their own clade. Based on analysis of ankyrin motifs, wPip strains from Hawai‘i, Kaua‘i, and Midway Atoll were most similar to wPip5 strains of Australasian origin. By contrast, Wolbachia strains from Culex quinquefasciatus collected in American Samoa were most similar to wPip3 strains of American origin. We detected a single Culex mosquito from Pu‘u Wa‘awa‘a on Hawai‘i Island that was infected with a unique wPip3 genotype. This discovery, plus a rarefaction analysis of genotypes from Kaua‘i and Hawai‘i Islands suggests that limited sampling may have underestimated diversity of wPip in our study. Mosquitoes infected with wPip5 and wPip3 are bidirectionally compatible with each other based on prior studies, which would support their ability to coexist within the same population on Hawai‘i Island. Available evidence from prior studies suggests that genotype wPip4 from Africa, the Middle East, Europe, and Asia is bidirectionally incompatible with genotype wPip5 and varies in compatibility with genotype wPip3 depending on geographic origin. Since wPip5 appears to be the most common compatibility group in Hawai‘i based on limited sampling, logical next steps are to 1) expand the current survey to include additional islands and localities, 2) infect a laboratory colony of Hawaiian Culex with wPip4 through tetracycline treatment of Hawaiian mosquitoes and backcross with Culex from Europe, North Africa, and the Middle East that are naturally infected with wPip4, 3) conduct cage trials to confirm bidirectional incompatibilities between Hawaiian Culex infected with wPip4 and wPip5, and 4) conduct field trials to evaluate whether release of incompatible males can be applied at small scales to suppress local populations.

American Samoa, Hawaii↗

Graphic and algebraic solutions of the discordant lead-uranium age problem

Uranium-bearing minerals that give lead-uranium and lead—lead ages that are essentially in agreement, i.e. concordant, generally are considered to have had a relatively simple geologic history and to have been unaltered since their deposition. The concordant ages obtained on such materials are, therefore, assumed to approach closely the actual age of the minerals. Many uranium-bearing samples, particularly uranium ores, give the following discordant age sequences; Pb 206 U 238 &lt; Pb 207 U 235 &#x2AA1; Pb 207 Pb 206 "> Pb 206 U 238 <Pb 207 U 235 ⪡Pb 207 Pb 206 or, less frequently, Pb 207 Pb 206 &#x2AA1; Pb 207 U 235 &lt; Pb 206 U 238 "> Pb 207 Pb 206 ⪡Pb 207 U 235 <Pb 206 U 238 . These discordant age sequences have been attributed most often to uncertainties in the common lead correction, selective loss of radio-active daughter products, loss or gain of lead or uranium, or contamination by an older generation of radiogenic lead. The evaluation of discordant lead isotope age data may be separated into two operations. The first operation, with which this paper is concerned, is mechanical in nature and involves the calculation of the different possible concordant ages corresponding to the various processes assumed to have produced the discordant ages. The second operation is more difficult to define and requires, in part, some personal judgement. It includes a synthesis of the possible concordant age solutions with other independent geologic and isotopic evidence. The concordant age ultimately chosen as most acceptable should be consistent not only with the known events in the geologic history of the area, the age relations of the enclosing rocks, and the mineralogic and paragenetic evidence, but also with other independent age measurements and the isotopic data obtained on the lead in related or associated non-radioactive minerals. The calculation of the possible concordant ages from discordant age data has been greatly simplified by Wetherill's graphical method of plotting the mole ratios of radiogenic Pb 206 U 238 "> Pb 206 U 238 ( N 206 N 238 "> N 206 N 238 ) vs. radiogenic Pb 207 U 235 "> Pb 207 U 235 ( N 207 N 235 "> N 207 N 235 ) after correcting for the contaminating common Pb 206 and Pb 207 . The linear relationships noted in this graphical procedure have been extended to plots of the mole ratios of total Pb 206 U 238 "> Pb 206 U 238 ( t N 206 N 238 "> t N 206 N 238 ) vs. total Pb 207 U 235 "> Pb 207 U 235 ( t N 207 N 235 "> t N 207 N 235 ). This modification permits the calculation of concordant ages for unaltered samples using only the Pb 207 Pb 206 "> Pb 207 Pb 206 ratio of the contaminating common lead. If isotopic data are available for two samples of the same age, x and y , from the same or related deposits or outcrops, graphs of the normalized difference ratios [ ( N 206 N 204 )x &#x2212; ( N 206 N 204 )y ( N 238 N 204 )x &#x2212;( N 238 N 204 )y ] vs. [ ( N 207 N 204 )x &#x2212; ( N 207 N 204 )y ( N 235 N 204 )x &#x2212;( N 235 N 204 )y ] "> [(N 206 N 204 )x − (N 206 N 204 )y(N 238 N 204 )x −(N 238 N 204 )y] vs. [(N 207 N 204 )x − (N 207 N 204 )y(N 235 N 204 )x −(N 235 N 204 )y] can give concordant ages corrected for unknown amounts of a common lead with an unknown Pb 207 / Pb 206 ratio. (If thorium is absent the difference ratios may be normalized with the more abundant index isotope, Pb 208 .) Similar plots of tho normalized, difference ratios for three genetically related samples ( x − y ) and( x − z ), will give concordant ages corrected, in addition, for either one unknown period of past alteration or initial contamination by an older generation of radiogenic lead of unknown Pb 207 /Pb 206 ratio. Practical numerical solutions for many of tho concordant age calculations are not currently available. However, the algebraic equivalents of these new graphical methods give equations which may be programmed for computing machines. For geologically probable parameters the equations of higher order have two positive real roots that rapidly converge on the exact concordant ages corrected for original radiogenic lead and for loss or gain of lead or uranium. Modifications of these general age equations expanded only to the second degree have been derived for use with desk calculators. These graphical and algebraic methods clearly suggest both the type and minimum number of samples necessary for adequate mathematical analysis of discordant lead isotope age data. This mathematical treatment also makes it clear that discordant lead isotope data alone cannot provide the basis for the choice of one of the possible concordant age solutions. The new equations, in particular, provide an incentive to improve our physical constants, analytical techniques and sampling methods in order that we may derive all of the useful geologic information that is available in a comprehensive lead isotope age study.

Geochimica et Cosmochimica Acta↗