SEARCH · Geology Reports

Results for “Information Bulletin”

Search indexed USGS publications on groundwater, aquifers, geologic maps, mineral resources and earthquakes. Explore source records by subject and place.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,027 records · Page 57Linked to original sources

USGS2025

Informing adaptive management to reduce ungulate aggregations: A case study involving winter feeding of elk

In the United States, wildlife managers are entrusted with preserving culturally and economically important ungulate populations in the face of the ongoing spread of chronic wasting disease (CWD). The U.S. Fish and Wildlife Service established an adaptive management plan to reduce the reliance of elk ( Cervus canadensis ) on supplemental winter feeding on the National Elk Refuge. The end goal of reducing the unnaturally high aggregation of elk during the winter is to mitigate the threat of disease outbreaks and to meet the objectives of sustainable populations of elk on the refuge. In this case study, we evaluated 6 years (2017–2022) of data from GPS-collared elk to determine the efficacy of shortening the length of the annual feeding period on the refuge to reduce elk aggregation. We measured aggregation using proximity rates, based on pairwise interactions over time, in both raw form as an index as well as predicted proximity as a function of other abiotic influences. We created a new R package, wildagg , to help with the process of computing the metrics from our study and to increase reproducibility in the future. Aggregation declined in years with less feeding on the refuge according to raw aggregation metrics when examined in isolation and dependent on the baseline feeding year used for comparison. However, accounting for abiotic factors while modeling proximity rates suggested in some years the decision to shorten the feeding period had less influence on aggregation than predicted. Our results underscore the complexity of measuring management outcomes and the usefulness of multiple approaches to evaluation.

Wyoming↗

USGS2017

Finite‐fault Bayesian inversion of teleseismic body waves

Inverting geophysical data has provided fundamental information about the behavior of earthquake rupture. However, inferring kinematic source model parameters for finite‐fault ruptures is an intrinsically underdetermined problem (the problem of nonuniqueness), because we are restricted to finite noisy observations. Although many studies use least‐squares techniques to make the finite‐fault problem tractable, these methods generally lack the ability to apply non‐Gaussian error analysis and the imposition of nonlinear constraints. However, the Bayesian approach can be employed to find a Gaussian or non‐Gaussian distribution of all probable model parameters, while utilizing nonlinear constraints. We present case studies to quantify the resolving power and associated uncertainties using only teleseismic body waves in a Bayesian framework to infer the slip history for a synthetic case and two earthquakes: the 2011 M w 7.1 Van, east Turkey, earthquake and the 2010 M w 7.2 El Mayor–Cucapah, Baja California, earthquake. In implementing the Bayesian method, we further present two distinct solutions to investigate the uncertainties by performing the inversion with and without velocity structure perturbations. We find that the posterior ensemble becomes broader when including velocity structure variability and introduces a spatial smearing of slip. Using the Bayesian framework solely on teleseismic body waves, we find rake is poorly constrained by the observations and rise time is poorly resolved when slip amplitude is low.

Bulletin of the Seismological Society of America↗

USGS2017

Long-term afterslip of the M6.0, 2004 Parkfield, California, earthquake—Implications for forecasting amount and duration of afterslip on other major creeping faults

We present the longest record of surface afterslip on a continental strike‐slip fault for the 2004 M 6.0 Parkfield, California, earthquake, from which we can derive critical information about the duration and predictability of afterslip relevant to urban displacement hazard applications. Surface slip associated with this event occurred entirely postseismically along the interseismically creeping (0.6–1.5 cm/yr) main trace of the San Andreas fault. Using the first year of afterslip data, the program AFTER correctly predicted the cumulative surface afterslip (maximum ∼35 cm) eventually attained. By 1 yr postearthquake, observed afterslip had accumulated to only ∼74% of its modeled final value u f in units of length. The 6‐yr data suggested final slip would be reached everywhere by ∼6–12 yrs. Parkfield’s afterslip lasted much longer (∼6–12 yrs) than afterslip following a 2014 M 6.0 event in Napa, California, where no interseismic creep was known, and its afterslip neared completion (∼97% of u f ) by 1 yr. The uncertainty in u f for the Napa event fell to ≤2 cm in only three months, versus in 2 yrs for the Parkfield event, mostly because duration of the power‐law stage of afterslip at Parkfield is much longer, ∼1000 (493–1666) days versus ∼100 (35–421) days for Napa. Because the urban Hayward fault near San Francisco, California, like the Parkfield section, exhibits interseismic creep in a similar geological regime, significant afterslip might last for up to a decade following an anticipated M ≥6.7 earthquake, potentially delaying postearthquake recovery.

California↗

USGS2013

Methods to estimate distribution and range extent of grizzly bears in the Greater Yellowstone Ecosystem

The distribution of the Greater Yellowstone Ecosystem grizzly bear ( Ursus arctos ) population has expanded into areas unoccupied since the early 20th century. Up-to-date information on the area and extent of this distribution is crucial for federal, state, and tribal wildlife and land managers to make informed decisions regarding grizzly bear management. The most recent estimate of grizzly bear distribution (2004) utilized fixed-kernel density estimators to describe distribution. This method was complex and computationally time consuming and excluded observations of unmarked bears. Our objective was to develop a technique to estimate grizzly bear distribution that would allow for the use of all verified grizzly bear location data, as well as provide the simplicity to be updated more frequently. We placed all verified grizzly bear locations from all sources from 1990 to 2004 and 1990 to 2010 onto a 3-km × 3-km grid and used zonal analysis and ordinary kriging to develop a predicted surface of grizzly bear distribution. We compared the area and extent of the 2004 kriging surface with the previous 2004 effort and evaluated changes in grizzly bear distribution from 2004 to 2010. The 2004 kriging surface was 2.4% smaller than the previous fixed-kernel estimate, but more closely represented the data. Grizzly bear distribution increased 38.3% from 2004 to 2010, with most expansion in the northern and southern regions of the range. This technique can be used to provide a current estimate of grizzly bear distribution for management and conservation applications.

Wildlife Society Bulletin↗

USGS2014

Processes of arroyo filling in northern New Mexico, USA

We documented arroyo evolution at the tree, trench, and arroyo scales along the lower Rio Puerco and Chaco Wash in northern New Mexico, USA. We excavated 29 buried living woody plants and used burial signatures in their annual rings to date stratigraphy in four trenches across the arroyos. Then, we reconstructed the history of arroyo evolution by combining trench data with arroyo-scale information from aerial imagery, light detection and ranging (LiDAR), longitudinal profiles, and repeat surveys of cross sections. Burial signatures in annual rings of salt cedar and willow dated sedimentary beds greater than 30 cm thick with annual precision. Along both arroyos, incision occurred until the 1930s in association with extreme high flows, and subsequent filling involved vegetation development, channel narrowing, increased sinuosity, and finally vertical aggradation. A strongly depositional sediment transport regime interacted with floodplain shrubs to produce a characteristic narrow, trapezoidal channel. The 55 km study reach along the Rio Puerco demonstrated upstream progression of arroyo widening and filling, but not of arroyo incision, channel narrowing, or floodplain vegetation development. We conclude that the occurrence of upstream progression within large basins like the Rio Puerco makes precise synchrony across basins impossible. Arroyo wall retreat is now mostly limited to locations where meanders impinge on the arroyo wall, forming hairpin bends, for which entry to and exit from the wall are stationary. Average annual sediment storage within the Rio Puerco study reach between 1955 and 2005 was 4.8 × 10 5 t/yr, 16% of the average annual suspended sediment yield, and 24% of the long-term bedrock denudation rate. At this rate, the arroyo would fill in 310 yr.

New Mexico↗

USGS2020

Understanding sportsperson retention and reactivation through license purchasing behavior

Most state and provincial fish and wildlife agencies have access to important information about patterns in sportsperson participation through their license databases. Using transaction data from Nebraska Game and Parks Commission's electronic hunting and fishing license system, we tracked license purchases of Nebraska, USA, resident license holders in 2010 through 2017. We categorized sportspersons by gender and yearly purchases as hunting only (Hunter), fishing only (Angler), a combination of hunting and fishing (Hunter–Angler), or no purchases (Inactive). The probability of movement among active sportsperson groups was limited and varied little based on initial group participation. The Angler group had the greatest probability of an individual transitioning to the Inactive group (females = 0.39; males = 0.33). The Hunter–Angler group had the greatest probability of an individual remaining within the same group (females = 0.65; males = 0.76). There was a relatively low probability of an individual in the Hunter group moving to the Angler group and vice versa (≤0.02). The sportsperson population is dynamic and understanding patterns of sportsperson participation is important for the future of fish and wildlife management in North America. Using data readily available to most fish and wildlife agencies has the potential to significantly improve our understanding of hunter and angler participation and aid management agencies and conservation organizations in the development of more effective strategies for managing sportspersons. © 2020 The Wildlife Society.

Wildlife Society Bulletin↗

USGS2013

Conditional spectrum computation incorporating multiple causal earthquakes and ground-motion prediction models

The conditional spectrum (CS) is a target spectrum (with conditional mean and conditional standard deviation) that links seismic hazard information with ground-motion selection for nonlinear dynamic analysis. Probabilistic seismic hazard analysis (PSHA) estimates the ground-motion hazard by incorporating the aleatory uncertainties in all earthquake scenarios and resulting ground motions, as well as the epistemic uncertainties in ground-motion prediction models (GMPMs) and seismic source models. Typical CS calculations to date are produced for a single earthquake scenario using a single GMPM, but more precise use requires consideration of at least multiple causal earthquakes and multiple GMPMs that are often considered in a PSHA computation. This paper presents the mathematics underlying these more precise CS calculations. Despite requiring more effort to compute than approximate calculations using a single causal earthquake and GMPM, the proposed approach produces an exact output that has a theoretical basis. To demonstrate the results of this approach and compare the exact and approximate calculations, several example calculations are performed for real sites in the western United States. The results also provide some insights regarding the circumstances under which approximate results are likely to closely match more exact results. To facilitate these more precise calculations for real applications, the exact CS calculations can now be performed for real sites in the United States using new deaggregation features in the U.S. Geological Survey hazard mapping tools. Details regarding this implementation are discussed in this paper.

Bulletin of the Seismological Society of America↗

USGS2011

Monitoring the Earthquake source process in North America

With the implementation of the USGS National Earthquake Information Center Prompt Assessment of Global Earthquakes for Response system (PAGER), rapid determination of earthquake moment magnitude is essential, especially for earthquakes that are felt within the contiguous United States. We report an implementation of moment tensor processing for application to broad, seismically active areas of North America. This effort focuses on the selection of regional crustal velocity models, codification of data quality tests, and the development of procedures for rapid computation of the seismic moment tensor. We systematically apply these techniques to earthquakes with reported magnitude greater than 3.5 in continental North America that are not associated with a tectonic plate boundary. Using the 0.02-0.10 Hz passband, we can usually determine, with few exceptions, moment tensor solutions for earthquakes with M w as small as 3.7. The threshold is significantly influenced by the density of stations, the location of the earthquake relative to the seismic stations and, of course, the signal-to-noise ratio. With the existing permanent broadband stations in North America operated for rapid earthquake response, the seismic moment tensor of most earthquakes that are M w 4 or larger can be routinely computed. As expected the nonuniform spatial pattern of these solutions reflects the seismicity pattern. However, the orientation of the direction of maximum compressive stress and the predominant style of faulting is spatially coherent across large regions of the continent.

Bulletin of the Seismological Society of America↗

USGS2006

Earthquake ground-motion prediction equations for eastern North America

New earthquake ground-motion relations for hard-rock and soil sites in eastern North America (ENA), including estimates of their aleatory uncertainty (variability) have been developed based on a stochastic finite-fault model. The model incorporates new information obtained from ENA seismographic data gathered over the past 10 years, including three-component broadband data that provide new information on ENA source and path effects. Our new prediction equations are similar to the previous ground-motion prediction equations of Atkinson and Boore (1995), which were based on a stochastic point-source model. The main difference is that high-frequency amplitudes (f ??? 5 Hz) are less than previously predicted (by about a factor of 1.6 within 100 km), because of a slightly lower average stress parameter (140 bars versus 180 bars) and a steeper near-source attenuation. At frequencies less than 5 Hz, the predicted ground motions from the new equations are generally within 25% of those predicted by Atkinson and Boore (1995). The prediction equations agree well with available ENA ground-motion data as evidenced by near-zero average residuals (within a factor of 1.2) for all frequencies, and the lack of any significant residual trends with distance. However, there is a tendency to positive residuals for moderate events at high frequencies in the distance range from 30 to 100 km (by as much as a factor of 2). This indicates epistemic uncertainty in the prediction model. The positive residuals for moderate events at < 100 km could be eliminated by an increased stress parameter, at the cost of producing negative residuals in other magnitude-distance ranges; adjustment factors to the equations are provided that may be used to model this effect.

Bulletin of the Seismological Society of America↗

USGS1981

Environmental features and macrofauna of Kahana Estuary, Oahu, Hawaii

Lack of ecological information on Hawaiian estuaries prompted an intensive 2-year study of a small (5.7 ha) stream-mouth estuary on windward Oahu. Water quality and macrofauna were sampled weekly at seven stations. The water mass was strongly stratified vertically except during freshets. Average values for water column temperature and bottom salinity were 23.2°C and 12&permil; at the head to 28.3°C and 28&permil; at the mouth. Dissolved oxygen saturation in the water column varied from about 50% at night to 140% in the afternoon. Usually, bottom waters were 3–6°C warmer than surface waters and sometimes showed severe oxygen depletion. Macrofauna, collected primarily by seining, consisted mainly of decapod crustaceans (four species of crabs, seven species of shrimps) and fishes (24 species). Other typical estuarine taxons (mollusks, barnacles, polychaetes) were scarce or absent. Diversity increased seaward from 14 species near the estuary head to 29 species near the mouth. Three species of crustaceans and six of fishes were captured at all stations. Most abundant were the native prawn, Macrobrachium grandimanus , and mullet, Mugil cephalus . Perennially resident adults occurred among crustaceans and gobioid fishes; most other fishes were present as juveniles and sporadic adults. Comparisons with other data suggest that more than 50 species of native fishes may occur in Hawaiian estuaries, and that estuarine macrofaunal diversity on oceanic islands is much lower than on continents at similar latitudes.

Hawaii↗

USGS2020

Evaluation of ground‐motion models for U.S. Geological Survey seismic hazard forecasts: Hawaii tectonic earthquakes and volcanic eruptions

The selection and weighting of ground‐motion models (GMMs) introduces a significant source of uncertainty in U.S. Geological Survey (USGS) National Seismic Hazard Modeling Project (NSHMP) forecasts. In this study, we evaluate 18 candidate GMMs using instrumental ground‐motion observations of horizontal peak ground acceleration (PGA) and 5%‐damped pseudospectral acceleration (0.02–10 s) for tectonic earthquakes and volcanic eruptions, to inform logic‐tree weights for the update of the USGS seismic hazard model for Hawaii. GMMs are evaluated using two methods. The first is a total residual visualization approach that compares the probability density function (PDF), mean and standard deviations σ "> σ ⁠ , of the observed and predicted ground motion. The second GMM evaluation method we use is the common total residual probabilistic scoring method (log likelihood [LLH]). The LLH method provides a single score that can be used to weight GMMs in the Hawaii seismic hazard model logic trees. The total residual PDF approach provides additional information by preserving GMM over‐ and underprediction across a broad spectrum of periods that is not available from a single value LLH score. We apply these GMM evaluation methods to two different data sets: (1) a database of instrumental ground motions from historic earthquakes in Hawaii from 1973 to 2007 ( ⁠ M w "> M w 4–7.3) and (2) available ground motions from recent earthquakes ( ⁠ M w "> M w 4–6.9) associated with 2018 Kilauea eruptions. The 2018 Kilauea sequence contains both volcanic eruptions and tectonic earthquakes allowing for statistically significant GMM comparisons of the two event classes. The Kilauea ground observations provide an independent data set allowing us to evaluate the predictive power of GMMs implemented in the new USGS nshmp‐haz software system. We evaluate GMM performance as a function of earthquake depth and we demonstrate that short‐period volcanic eruption ground motions are not well predicted by any candidate GMMs. Nine of the initial 18 candidate GMMs fit the observed ground motions and meet established criteria for inclusion in the update of the Hawaii seismic hazard model. A weighted mean of four top performing GMMs in this study (NGAsubslab, NGAsubinter, ASK14, A10) is 50% lower for PGA than for GMMS used in the previous USGS seismic hazard model for Hawaii.

Hawaii↗

USGS2025

Complex staged emplacement of a basaltic lava: The example of the July 1974 flow of Kīlauea

Basaltic lava flows can be highly destructive. Forecasting the future path and/or behavior of an active lava flow is challenging because topography is often poorly constrained and lava has a complex rheology and emplacement history. Preserved lavas are an important source of information which, combined with observations of active flows, underpins conceptual models of lava flow emplacement. However, the value of preserved lavas is limited because pre-eruptive topography and, thus, syn-eruptive lava flow geometry are usually not known. Here, we use tree-mold data to constrain pre-eruptive topography and syn-eruptive lava flow geometry of the July 1974 flow of Kīlauea (USA). Tree molds, which are formed after advancing lava encloses standing trees, preserve the lava inundation height and the final preserved thickness of lava. We used data from 282 tree molds to reconstruct the temporal and spatial evolution of the ~ 2.1 km-long July 1974 flow. The tree mold dataset yields a detailed dynamic picture of staged emplacement, separated by intervals of ponding. In some ponded areas, flow depth during emplacement (~ 5 m) was twice the preserved thickness of the final lava (2–3 m). Drainage of the ponds led to episodic surges in flow advancement, decoupled from fluctuations in vent discharge rate. We infer that the final breakout occurred after the cessation of fountaining. Such complex emplacement histories may be common for pāhoehoe lavas at Kīlauea and elsewhere in situations where the terrain is of variable slope, and/or where lava is temporarily perched and stored.

Hawaii↗

USGS2020

Factors influencing the probability of hydraulic fracturing induced seismicity in Oklahoma

Injection‐induced seismicity became an important issue over the past decade, and although much of the rise in seismicity is attributed to wastewater disposal, a growing number of cases have identified hydraulic fracturing (HF) as the cause. A recent study identified regions in Oklahoma where ≥75% of seismicity from 2010 to 2016 correlated with nearly 300 HF wells. To identify factors associated with increased probability of induced seismicity, we gathered publicly available information about the HF operations in these regions including: injected volume, number of wells on a pad, injected fluid (gel vs. slickwater), vertical depth of the well, proximity of the well to basement rock, and the formation into which the injection occurred. To determine the statistical strength of the trends, we applied logistic regression, bootstrapping, and odds ratios. We see no trend with total injected volume in our Oklahoma dataset, in contrast to strong trends observed in Alberta and Texas, but we note those regions have many more multiwell pads leading to larger cumulative volumes within a localized area. We found a > ∼50% lower probability of seismicity with the use of gel compared to slickwater. We found that HF wells targeting older formations had a higher probability of seismicity; however, these wells also tend to be deeper, and we found the trend with well depth to be stronger than the trend with age of formation. When isolated to the Woodford formation, well depth produced the strongest relationship, increasing from ∼5%> to ∼50% probability from 1.5 to 5.5 km. However, no trend was seen in the proximity to basement parameter. Based on previously measured pore pressure gradients, we interpret the strong absolute depth relationship to be a result of the increasing formation overpressure measured in deeper portions of the basin that lower the stress change needed to induce seismicity.

Oklahoma↗

USGS2017

Exploring the historical earthquakes preceding the giant 1960 Chile earthquake in a time‐dependent seismogenic zone

New documentary findings and available paleoseismological evidence provide both new insights into the historical seismic sequence that ended with the giant 1960 south-central Chile earthquake and relevant information about the region’s seismogenic zone. According to the few available written records, this region was previously struck by earthquakes of varying size in 1575, 1737, and 1837. We expanded the existing compilations of the effects of the two latter using unpublished first-hand accounts found in archives in Chile, Peru, Spain, and New England. We further investigated their sources by comparing the newly unearthed historical data and available paleoseismological evidence with the effects predicted by hypothetical dislocations. The results reveal significant differences in the along-strike and depth distribution of the ruptures in 1737, 1837, and 1960. While the 1737 rupture likely occurred in the northern half of the 1960 region, on a narrow and deep portion of the megathrust, the 1837 rupture occurred mainly in the southern half and slipped over a wide range of depth. Such a wide rupture in 1837 disagrees with the narrow and shallow seismogenic zone currently inferred along this region. If in fact there is now a narrow zone where 200 years ago there was a wider one, it means that the seismogenic zone changes with time, perhaps between seismic cycles. Such change probably explains the evident variability in both size and location of the great earthquakes that have struck this region over the last centuries, as evidenced by written history, and through millennia, as inferred from paleoseismology.

Bulletin of the Seismological Society of America↗

USGS2019

Informing sea turtle outreach efforts to maximize effectiveness

Most sea turtle (Cheloniidae) species worldwide are endangered or threatened, with threats causing harm to sea turtles predominantly human‐induced. Thus, prevention of further declines to these imperiled species will require alteration of human behaviors. Regulations, incentives, and environmental education are 3 strategies that could be used to alter human behavior. Our goal was to determine how to maximize effectiveness of one of these strategies—education efforts. We investigated knowledge deficiencies and light pollution behaviors of individuals living in a region with nesting sea turtles, in an effort to determine the best approach to promote sea turtle conservation. During 2014, we mailed a survey to 3,000 property owners in 4 coastal counties in Florida, USA, to achieve 3 objectives: assess what topic areas were misunderstood; discern who had knowledge deficiencies; and determine who had adopted turtle‐friendly lighting practices. The best predictors of knowledge included geographic factors (county, proximity of residences to the beach), demographic characteristics (age), and behaviors (individual's beach visitation rates). One practice that can reduce harm to sea turtles was common: use of window treatments to reduce light pollution. However, other practices harmful to sea turtles were prevalent, including long durations of use of outdoor lighting and use of light bulbs with wavelengths that can disturb sea turtles. Our results suggest that educational efforts could be enhanced by specifically focusing on increasing awareness of the effects of human actions on sea turtles, targeting individuals who visit the beach infrequently and live far from it to foster greater connection with these ecosystems, and publicizing a variety of options that could reduce harm to sea turtles so individuals feel a sense of freedom of choice.

Florida↗

USGS2003

Predicting mercury in mallard ducklings from mercury in chorioallantoic membranes

Methylmercury has been suspected as a cause of impaired reproduction in wild birds, but the confounding effects of other environmental stressors has made it difficult to determine how much mercury in the eggs of these wild species is harmful. Even when a sample egg can be collected from the nest of a wild bird and the mercury concentration in that egg compared to the laboratory-derived thresholds for reproductive impairment, additional information on the mercury levels in other eggs from that nest would be helpful in determining whether harmful levels of mercury were present in the clutch. The measurement of mercury levels in chorioallantoic membranes offers a possible way to estimate how much mercury was in a chick that hatched from an egg, and also in the whole fresh egg itself. While an embryo is developing, wastes are collected in a sac called the chorioallantoic membranes, which often remain inside the eggshell and can be collected for contaminant analysis. We fed methylmercury to captive mallards to generate a broad range of mercury levels in eggs, allowed the eggs to hatch normally, and then compared mercury concentrations in the hatchling versus the chorioallantoic membranes left behind in the eggshell. When the data from eggs laid by mercury- treated females were expressed as common logarithms, a linear equation was created by which the concentration of mercury in a duckling could be predicted from the concentration of mercury in the chorioallantoic membranes from the same egg. Therefore, if it were not possible to collect a sample egg from a clutch of wild bird eggs, the collection of the chorioallantoic membranes could be substituted, and the mercury predicted to be in the chick or whole egg could be compared to the thresholds of mercury that have been shown to cause harm in controlled feeding studies with pheasants, chickens, and mallards.

Bulletin of Environmental Contamination and Toxico↗

USGS1999

Stratigraphic framework of Holocene volcaniclastic deposits, Akutan Volcano, east-central Aleutian Islands, Alaska

Akutan Volcano is one of the most active volcanoes in the Aleutian arc, but until recently little was known about its history and eruptive character. Following a brief but sustained period of intense seismic activity in March 1996, the Alaska Volcano Observatory began investigating the geology of the volcano and evaluating potential volcanic hazards that could affect residents of Akutan Island. During these studies new information was obtained about the Holocene eruptive history of the volcano on the basis of stratigraphic studies of volcaniclastic deposits and radiocarbon dating of associated buried soils and peat. A black, scoria-bearing, lapilli tephra, informally named the 'Akutan tephra,' is up to 2 m thick and is found over most of the island, primarily east of the volcano summit. Six radiocarbon ages on the humic fraction of soil A-horizons beneath the tephra indicate that the Akutan tephra was erupted approximately 1611 years B.P. At several locations the Akutan tephra is within a conformable stratigraphic sequence of pyroclastic-flow and lahar deposits that are all part of the same eruptive sequence. The thickness, widespread distribution, and conformable stratigraphic association with overlying pyroclastic-flow and lahar deposits indicate that the Akutan tephra likely records a major eruption of Akutan Volcano that may have formed the present summit caldera. Noncohesive lahar and pyroclastic-flow deposits that predate the Akutan tephra occur in the major valleys that head on the volcano and are evidence for six to eight earlier Holocene eruptions. These eruptions were strombolian to subplinian events that generated limited amounts of tephra and small pyroclastic flows that extended only a few kilometers from the vent. The pyroclastic flows melted snow and ice on the volcano flanks and formed lahars that traveled several kilometers down broad, formerly glaciated valleys, reaching the coast as thin, watery, hyperconcentrated flows or water floods. Slightly cohesive lahars in Hot Springs valley and Long valley could have formed from minor flank collapses of hydrothermally altered volcanic bedrock. These lahars may be unrelated to eruptive activity.

Bulletin of Volcanology↗

USGS2002

Spent shot availability and ingestion on areas managed for mourning doves

Mourning dove (Zenaida macroura) hunting is becoming increasingly popular, especially in managed shooting fields. Given the possible increase in the availability of lead (Pb) shot on these areas, our objective was to estimate availability and ingestion of spent shot at the Eagle Bluffs Conservation Area (EBCA, hunted with nontoxic shot) and the James A. Reed Memorial Wildlife Area (JARWA, hunted with Pb shot) in Missouri. During 1998, we collected soil samples one or 2 weeks prior to the hunting season (prehunt) and after 4 days of dove hunting (posthunt). We also collected information on number of doves harvested, number of shots fired, shotgun gauge, and shotshell size used. Dove carcasses were collected on both areas during 1998-99. At EBCA, 60 hunters deposited an estimated 64,775 pellets/ha of nontoxic shot on or around the managed field. At JARWA, approximately 1,086,275 pellets/ha of Pb shot were deposited by 728 hunters. Our posthunt estimates of spent-shot availability from soil sampling were 0 pellets/ha for EBCA and 6,342 pellets/ha for JARWA. Our findings suggest that existing soil sampling protocols may not provide accurate estimates of spent-shot availability in managed dove shooting fields. During 1998-99, 15 of 310 (4.8%) mourning doves collected from EBCA had ingested nontoxic shot. Of those doves, 6 (40.0%) contained a?Y7 shot pellets. In comparison, only 2 of 574 (0.3%) doves collected from JARWA had ingested Pb shot. Because a greater proportion of doves ingested multiple steel pellets compared to Pb pellets, we suggest that doves feeding in fields hunted with Pb shot may succumb to acute Pb toxicosis and thus become unavailable to harvest, resulting in an underestimate of ingestion rates. Although further research is needed to test this hypothesis, our findings may partially explain why previous studies have shown few doves with ingested Pb shot despite their feeding on areas with high Pb shot availability.

Missouri↗