Geology Reports⌕ Search

SEARCH · Geology Reports

Results for “Modeling Report”

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 937 records · Page 52Linked to original sources

Hydrostratigraphic interpretation of test-hole and borehole geophysical data, Kimball, Cheyenne, and Deuel Counties, Nebraska, 2011-12

Recently (2004) adopted legislation in Nebraska requires a sustainable balance between long-term supplies and uses of surface-water and groundwater and requires Natural Resources Districts to understand the effect of groundwater use on surface-water systems when developing a groundwater-management plan. The South Platte Natural Resources District (SPNRD) is located in the southern Nebraska Panhandle and overlies the nationally important High Plains aquifer. Declines in water levels have been documented, and more stringent regulations have been enacted to ensure the supply of ground-water will be sufficient to meet the needs of future generations. Because an improved understanding of the hydrogeologic characteristics of this aquifer system is needed to ensure sustainability of groundwater withdrawals, the U.S. Geological Survey, in cooperation with the SPNRD, Conservation and Survey Division of the University of Nebraska-Lincoln, and the Nebraska Environmental Trust, began a hydrogeologic study of the SPNRD to describe the lithology and thickness of the High Plains aquifer. This report documents these characteristics at 29 new test holes, 28 of which were drilled to the base of the High Plains aquifer. Herein the High Plains aquifer is considered to include all hydrologically connected units of Tertiary and Quaternary age. The depth to the base of aquifer was interpreted to range from 37 to 610 feet in 28 of the 29 test holes. At some locations, particularly northern Kimball County, the base-of-aquifer surface was difficult to interpret from drill cutting samples and borehole geophysical logs. The depth to the base of aquifer determined for test holes drilled for this report was compared with the base-of-aquifer surface interpreted by previous researchers. In general, there were greater differences between the base-of-aquifer elevation reported herein and those in previous studies for areas north of Lodgepole Creek compared to areas south of Lodgepole Creek. The largest difference was at test hole 5-SP-11, where an Ogallala-filled paleovalley prevously had been interpreted based on relatively sparse test-hole data west of 5-SP-11. The base of aquifer near test hole 5-SP-11 reported herein is approximately 230 ft higher in elevation than previously interpreted. Among other test holes that are likely to have been drilled in Ogallala-filled paleovalleys, the greatest difference in the interpreted base of aquifer was for test hole 7-CC-11, northeast of Potter, Nebraska, where the base of aquifer is 180 feet deeper than previously interpreted. Interpretation of test-hole and borehole geophysical data for 29 additional test holes will improve resource managers’ understanding of the hydrogeologic characteristics, including aquifer thickness. Aquifer thickness, which is related to total water in storage, is not well quantified in the north and south tablelands. The additional hydrostratigraphic interpretations provided in this report will improve the hydrogeologic framework used in current (2014) and future groundwater models, which are the basis for many water-management decisions.

Nebraska↗

Age, chemistry, and correlations of Neoproterozoic–Devonian igneous rocks of the Arctic Alaska–Chukotka terrane: An overview with new U-Pb ages

The Arctic Alaska–Chukotka terrane is a microcontinent with an origin exotic to Laurentia. We used a sensitive high-resolution ion microprobe (SHRIMP) to date nine samples of Neoproterozoic rock and five samples of Devonian rock from the Brooks Range and Seward Peninsula of Alaska and from the Chukotka Peninsula of northeastern Russia. Felsic magmatism occurred at 968 Ma and 742 Ma in the Brooks Range and at 865 Ma and 670–666 Ma on Seward Peninsula. Felsic igneous rocks in Chukotka were dated at 656 Ma and 574 Ma. Devonian igneous rocks are found throughout the Arctic Alaska–Chukotka terrane, and we dated samples with ages of 391 Ma, 390 Ma, 385 Ma, 371 Ma, and 363 Ma. The felsic character of the Neoproterozoic rocks suggests formation at least in part through crustal melting. The age of the crustal source rocks that melted to form the Neoproterozoic rocks is inferred to be Mesoproterozoic based on Nd model ages ranging from 1.6 to 1.4 Ga. Rocks of this age range have been reported from the basement of Baltica but are rare in Laurentia. The 565 Ma orthogneisses on Seward Peninsula have ca. 1.1 Ga Nd model ages. Devonian igneous rocks have a wide range of model ages ranging from 1.6 to 0.8 Ga. The tectonic setting of the 968 Ma, 865 Ma, and 742 Ma rocks is unknown. The ca. 670 Ma magmatism on Seward Peninsula is interpreted to have occurred in an arc setting based on geochemistry and similarities in their ages to the Avalonian–Cadomian arc system peripheral to Gondwana. Latest Neoproterozoic magmatism is inferred to have occurred in a rift setting based on composition and the Paleozoic passive margin sequence that was deposited across the Arctic Alaska–Chukokta terrane. Devonian magmatism likely occurred in an arc and/or backarc rift setting. Significant uncertainties remain concerning the age of the Arctic Alaska–Chukotka terrane basement, particularly the age of the host rocks for Neoproterozoic intrusions.

Alaska↗

Implementation of the century ecosystem model for an eroding hillslope in Mississippi

The objective of this study was to parameterize and implement the Century ecosystem model for an eroding, cultivated site near Senatobia, in Panola County, Mississippi, in order to understand the loss and replacement of soil organic carbon on an eroding cropland. The sites chosen for this study are located on highly eroded loess soils where USDA has conducted studies on rates of soil erosion. We used USDA sediment data from the study site and historical erosion estimates from the nearby area as model input for soil loss; in addition, inputs for parametization include particle-size data, climate data, and rainfall/runoff data that were collected and reported in companion papers. A cropping scenario was implemented to simulate a research site at the USDA watershed 2 at the Nelson Farm. Model output was compiled for comparison with data collected and reported in companion reports; interpretive comparisons are reported in Harden et al, in press.

Open-File Report↗

Simulation of ground-water flow in glaciofluvial aquifers in the Grand Rapids area, Minnesota

A calibrated steady-state, finite-difference, ground-waterflow model was constructed to simulate ground-water flow in three glaciofluvial aquifers, defined in this report as the upper, middle, and lower aquifers, in an area of about 114 mi2 surrounding the city of Grand Rapids in north-central Minnesota. The calibrated model will be used by Minnesota Department of Health and communities in the Grand Rapids area in the development of wellhead protection plans for their water supplies. The model was calibrated through comparison of simulated ground-water levels to measured static water levels in 351 wells, and comparison of simulated base-flow rates to estimated base-flow rates for reaches of the Mississippi and Prairie Rivers. Model statistics indicate that the model tends to overestimate ground-water levels. The root mean square errors ranged from +12.83 ft in wells completed in the upper aquifer to +19.10 ft in wells completed in the middle aquifer. Mean absolute differences between simulated and measured water levels ranged from +4.43 ft for wells completed in the upper aquifer to +9.25 ft for wells completed in the middle aquifer. Mean algebraic differences ranged from +9.35 ft for wells completed in the upper aquifer to +14.44 ft for wells completed in the middle aquifer, with the positive differences indicating that the simulated water levels were higher than the measured water levels. Percentage errors between simulated and estimated base-flow rates for the three monitored reaches all were less than 10 percent, indicating good agreement. Simulated ground-water levels were most sensitive to changes in general-head boundary conductance, indicating that this characteristic is the predominant model input variable controlling steady-state water-level conditions. Simulated groundwater flow to stream reaches was most sensitive to changes in horizontal hydraulic conductivity, indicating that this characteristic is the predominant model input variable controlling steady-state flow conditions.

Minnesota↗

Groundwater source water assessment area delineation of public supply wells in the Treasure Valley, western Idaho

The Idaho Department of Environmental Quality (IDEQ) is responsible for assessing the source water for public water supply systems including the delineation of source water assessment areas for wells, which are defined as the subsurface areas contributing water to a well during a period of time. Part of the IDEQ assessment plan is to update source water assessments over time. The U.S. Geological Survey (USGS) developed a new groundwater-flow model for part of the western Snake River Plain aquifer in southwestern Idaho. This groundwater- flow model provided a new opportunity for IDEQ to update source water assessment for public water supply wells within the model boundary. The USGS, in cooperation with IDEQ, used the USGS Treasure Valley Groundwater Flow Model along with a MODPATH particle tracking post-processing program to reassess public-water supply source water assessment areas. The USGS created an automated script that requests a well location, depth, and pumpage rate, and then incorporates that information into the Treasure Valley Groundwater Flow Model combined with the MODPATH model to create 3-, 6-, and 10-year time of travel zones that are output to a geographic information system shapefile. The TVSWA_Delineator script provides a fast, transparent tool for IDEQ to reassess or create new public-water supply source water assessment area delineations based on new and updated information and is available as a Python script and a Windows executable with the compiled script and all necessary components. TVSWA_Delineator utilizes a quadtree approach to rediscretize the Treasure Valley Groundwater Flow Model grid around the pumping well location from 1-mile cell lengths telescoped down to 1/16 mile. This rediscretization increases the source water assessment area delineation resolution near the pumping well and minimizes the effect of weak sinks in the MODPATH model. The stress inputs for each month of the final 10 years (2006 to 2015) of the Treasure Valley Groundwater Flow model historical simulation were averaged for the particle tracking analysis. The MODPATH model used backwards tracking of particles released at each stress period to create the source water assessment area delineations. This report documents the creation of the TVSWA_Delineator Python script, including adaptations made to the Treasure Valley Groundwater Flow Model, creation of a companion MODPATH model, and post-processing results to provide source water assessment areas.

Idaho↗

Approaches in highly parameterized inversion - GENIE, a general model-independent TCP/IP run manager

GENIE is a model-independent suite of programs that can be used to generally distribute, manage, and execute multiple model runs via the TCP/IP infrastructure. The suite consists of a file distribution interface, a run manage, a run executer, and a routine that can be compiled as part of a program and used to exchange model runs with the run manager. Because communication is via a standard protocol (TCP/IP), any computer connected to the Internet can serve in any of the capacities offered by this suite. Model independence is consistent with the existing template and instruction file protocols of the widely used PEST parameter estimation program. This report describes (1) the problem addressed; (2) the approach used by GENIE to queue, distribute, and retrieve model runs; and (3) user instructions, classes, and functions developed. It also includes (4) an example to illustrate the linking of GENIE with Parallel PEST using the interface routine.

Techniques and Methods↗

Reconnaissance investigation of the rough diamond resource potential and production capacity of Côte d’Ivoire

Ethnic and political conflict developed into open civil war in Côte d’Ivoire in 2002, leading to a de facto partitioning of the country into the government-controlled south and the rebel-controlled north. Côte d’Ivoire’s two main diamond mining areas, Séguéla and Tortiya, are located in the north, under what was, until recently, rebel-controlled territory. In an effort to prevent proceeds from diamond mining from funding the conflict, the United Nations (UN) placed an embargo on the export of rough diamonds from Côte d’Ivoire in 2005. That same year, the Kimberley Process (KP), the international initiative charged with stemming the flow of conflict diamonds, acted to enforce this ban by adopting the Moscow Resolution on Côte d’Ivoire, which contained measures to prevent the infiltration of Ivorian diamonds into the legitimate global rough diamond trade. Though under scrutiny by the international community, diamond mining activities continued in Côte d’Ivoire, with artisanal miners exploiting both alluvial deposits in fluvial systems and primary kimberlitic dike deposits. However, because of the embargo, there has been no official record of diamond production since the conflict began in 2002. This lack of production statistics represents a significant data gap and hinders efforts by the KP to understand how illicitly produced diamonds may be entering the legitimate trade. This study presents the results of a multiyear effort to monitor the diamond mining activities of Côte d’Ivoire’s two main diamond mining areas, Séguéla and Tortiya. An innovative approach was developed that integrates data acquired from archival reports and maps, high-resolution satellite imagery, and digital terrain modeling to assess the total diamond endowment of the Séguéla and Tortiya deposits and to calculate annual diamond production from 2006 to 2013. On the basis of currently available data, this study estimates that a total of 10,100,000 carats remain in Séguéla and a total of 1,100,000 carats remain in Tortiya. Production capacity was calculated for the two study areas for the years 2006–2010 and 2012–2013. Production capacity was found to range from between 38,000 carats and 375,000 carats in Séguéla and from 13,000 carats and 20,000 carats in Tortiya. Further, this study demonstrates that artisanal mining activities can be successfully monitored by using remote sensing and geologic modeling techniques. The production capacity estimates presented here fill a significant data gap and provide policy makers, the UN, and the KP with important information not otherwise available.

Scientific Investigations Report↗

Recent and possible future variations in the North American Monsoon

The dynamics and recent and possible future changes of the June–September rainfall associated with the North American Monsoon (NAM) are reviewed in this chapter. Our analysis as well as previous analyses of the trend in June–September precipitation from 1948 until 2010 indicate significant precipitation increases over New Mexico and the core NAM region, and significant precipitation decreases over southwest Mexico. The trends in June–September precipitation have been forced by anomalous cyclonic circulation centered at 15°N latitude over the eastern Pacific Ocean. The anomalous cyclonic circulation is responsible for changes in the flux of moisture and the divergence of moisture flux within the core NAM region. Future climate projections using the Coupled Model Intercomparison Project Phase 5 (CMIP5) models, as part of the Intergovernmental Panel on Climate Change (IPCC) Fifth Assessment Report (AR5), support the observed analyses of a later shift in the monsoon season in the presence of increased greenhouse gas concentrations in the atmosphere under the RCP8.5 scenario. The CMIP5 models under the RCP8.5 scenario predict significant NAM-related rainfall decreases during June and July and predict significant NAM-related rainfall increases during September and October.

Book chapter↗

Geologic map of Kundelan ore deposits and prospects, Zabul Province, Afghanistan; modified from the 1971 original map compilations of K.I. Litvinenko and others

This map and cross sections are redrafted modified versions of the Geological map of the Kundelan ore deposit area, scale 1:10,000 (graphical supplement no. 18) and the Geological map of the Kundelan deposits, scale 1:2,000 (graphical supplement no. 3) both contained in an unpublished Soviet report by Litvinenko and others (1971) (report no. 0540). The unpublished Soviet report was prepared in cooperation with the Ministry of Mines and Industries of the Royal Government of Afghanistan in Kabul during 1971. This redrafted map and cross sections illustrate the geology of the main Kundelan copper-gold skarn deposit, located within the Kundelan copper and gold area of interest (AOI), Zabul Province, Afghanistan. Areas of interest (AOIs) of non-fuel mineral resources within Afghanistan were first described and defined by Peters and others (2007) and later by the work of Peters and others (2011a). The location of the main Kundelan copper-gold skarn deposit (area of this map) and the Kundelan copper and gold AOI is shown on the index map provided on this map sheet. The estimated resources of the Kundelan copper-gold skarn deposit are 21,400 metric tons (t) of copper, 1.6 t of gold, and 133.4 t of molybdenum at an average grade of 1.21 weight percent (wt. %) copper (ranging from 0.66 to 4.03 wt. % copper); 0.9 grams per metric ton (g/t) gold (ranging from 0.3 to 3.1 g/t gold); 0.14 wt. % molybdenum; and as much as 10 g/t silver and 0.03 wt. % bismuth (Peters and others, 2011b). Small past production of gold and base metals is also reported by Douvgal and others (1971) from many prospects within the Kundelan copper and gold AOI. Outside the skarn areas, argillic hydrothermal alteration is present (Abdullah and others, 1977). Most copper and gold prospects in the Kundelan copper and gold AOI are reported to contain commercial-grade ores of copper and (or) gold, and many prospect areas have potential for these commodities to be discovered in commercial volumes. Future initial mine exploration and later development in many of the prospects, and specifically in the Kundelan copper-gold skarn deposit, could result in near-term small- to medium-sized gold mining operations (Peters and others, 2011b). The redrafted map and cross sections reproduce the topology of rock units, contacts, faults, and so forth, of the original Soviet map and cross sections, and they include modifications based on our examination of these documents and our observations made during a brief field visit in August of 2010. We have attempted to translate the original Russian terminology and rock classifications into modern English geologic usage as literally as possible without changing any genetic or process-oriented implications in the original descriptions. We also use the age designations from the original Soviet maps, except for the phase I and II intrusive igneous rocks. Phase I and II igneous rocks are reassigned an Early Cretaceous age (from lower Paleogene) based on a uranium-lead (U-Pb) zircon SHRIMP analysis from quartz diorite, which yielded an age of 104±1 Ma (mega-annum). The information provided in the description of map units is from both the source report by Litvinenko and others (1971) and the original 1:10,000 scale map (graphical supplement no. 18), also from Litvinenko and others (1971). Because of the poor quality of the original map, some map features could not be identified and some features may be misinterpreted. The rock unit colors used on the redrafted maps and cross sections differ from the colors shown on the original Soviet version. Colors were selected according to the color and pattern scheme of the Commission for the Geological Map of the World (CGMW) at http://www.ccgm.org. Elevations on the cross sections are derived from the original Soviet topography and may not match the Global Digital Elevation Model (GDEM) topography used on the redrafted map of this report. Most hydrography derived from the Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER) has not been included on our redrafted version of the map because of a poor fit with alluvial deposits from the unmodified original Soviet map (graphical supplement no. 18; Litvinenko and others, 1971).

Zabul Province↗

A bibliography of terrain modeling (geomorphometry), the quantitative representation of topography: Supplement 4.0

Terrain modeling, the practice of ground-surface quantification, is an amalgam of Earth science, mathematics, engineering, and computer science. The discipline is known variously as geomorphometry (or simply morphometry), terrain analysis, and quantitative geomorphology. It continues to grow through myriad applications to hydrology, geohazards mapping, tectonics, sea-floor and planetary exploration, and other fields. Dating nominally to the co-founders of academic geography, Alexander von Humboldt (1808, 1817) and Carl Ritter (1826, 1828), the field was revolutionized late in the 20th Century by the computer manipulation of spatial arrays of terrain heights, or digital elevation models (DEMs), which can quantify and portray ground-surface form over large areas (Maune, 2001). Morphometric procedures are implemented routinely by commercial geographic information systems (GIS) as well as specialized software (Harvey and Eash, 1996; Köthe and others, 1996; ESRI, 1997; Drzewiecki et al., 1999; Dikau and Saurer, 1999; Djokic and Maidment, 2000; Wilson and Gallant, 2000; Breuer, 2001; Guth, 2001; Eastman, 2002). The new Earth Surface edition of the Journal of Geophysical Research, specializing in surficial processes, is the latest of many publication venues for terrain modeling. This is the fourth update of a bibliography and introduction to terrain modeling (Pike, 1993, 1995, 1996, 1999) designed to collect the diverse, scattered literature on surface measurement as a resource for the research community. The use of DEMs in science and technology continues to accelerate and diversify (Pike, 2000a). New work appears so frequently that a sampling must suffice to represent the vast literature. This report adds 1636 entries to the 4374 in the four earlier publications1. Forty-eight additional entries correct dead Internet links and other errors found in the prior listings. Chronicling the history of terrain modeling, many entries in this report predate the 1999 supplement. Coverage is representative from about 1800 through early–mid 2002. Papers increasingly are published exclusively or in duplicate on the Internet's World Wide Web; the dates given here for Web addresses (URLs) that lack a print publication indicate a Web site's last update or my last access of it. The bibliography is arranged alphabetically and thus is not readily summarized. This introduction cites about 500 entries, a third of them grouped under 24 morphometric topics, as a guide to the listing's contents. Continuing the practice of previous bibliographies in the series to provide more information on a few applications (see summary of past topics in Pike, 2000a), this report elaborates further on topographic data, putative new parameters, tectonic geomorphology/neo-orometry, biogeography, ice-cap morphometry, results from the Mars Global DEM, landslide-hazard mapping, terrain modeling as physics, Hack's law, and broad-scale computer visualization. The literature of some of these subjects is large, and none of the summaries is intended to more than introduce the topic and comment on some of the current contributions of terrain modeling. Closing the essay is a discussion of pre-1900 papers that trace the evolution of ridge-line and watercourse quantification by descriptive geometry, as well as comments on some new books and an on-line bulletin board.

Open-File Report↗

Coexistence and community structure of tropical trees in a Hawaiian montane rain forest

We measured the diameter at breast height of all trees and shrubs > 5 meters in height, including standing dead trees, on 68 0.04-hectare study plots in a montane, subtropical rain forest on Mauna Loa, Hawai`i. The canopy species consisted of 88 percent Metrosideros polymorpha (ohia) and 12 percent Acacia koa (koa). Negative associations were found between the densities of koa and ohia, the density of koa and the total basal area of ohia, and the total basal areas of koa and ohia. The two-species lottery competition model, a stochastic model in which the coexistence of two species in a space-limited community results from temporal variation in recruitment and death rates, predicts a quadratic-beta distribution for the proportion of space occupied by each species. A discrete version of the quadratic-beta distribution, the quadratic-beta binomial distribution, was fit to the live koa and ohia densities and assessed with goodness-of-fit tests. Likelihood ratio tests provided evidence that the mean adult death rates of the two species were equal but that the relative competitive abilities of the two species favored ohia. These tests were corroborated by a contingency table analysis of death rates based on standing dead trees and growth rate studies which report that koa grows much faster than ohia. The lottery model predicts a positive covariance between death rates and ohia recruitment when mean death rates are equal and koa has a higher growth rate than ohia. We argue that the competitive advantage of ohia is due to its superior dispersal ability into large gaps, which would yield the positive covariance described above, and it is this positive covariance term that skews the occupation of space in favor of ohia.

Hawaii↗

Basin, climatic, and irrigation factors associated with median summer water yields for streams in Southwestern Michigan, 1945-2015

Median summer water yields and resultant flows for streams are used in Michigan to regulate large water withdrawals to help prevent negative effects on characteristic fish populations. Large water withdrawals commonly are associated with irrigation in rural areas. In an earlier statewide report, an index-flow statistic for the period of record, defined as the median flow during the summer month of lowest flow, was used to characterize median summer flows and associated water yields. In this report, the annual series of median summer water yields for the period July 1 through September 30 within the period of record is used to characterize median summer water yields. For 27 streamgages included in both reports, the average index water yield was at the 37th percentile of the distribution of median summer water yields. In contrast to an index statistic, an annual time series provides a basis for detecting trends in median summer water yields and for determining basin, climatic, and irrigation factors affecting spatial and temporal variations in summer water yields. Daily flow data from 40 selected U.S. Geological Survey streamgages in southwestern Michigan were used in this analysis. Two mixed models were identified to estimate median summer water yields based on fixed basin characteristics and temporally varying climatic factors for 1945–2015. No irrigation data were available prior to 1970, so no irrigation variables were included in the mixed models for 1945–2015. Then, two mixed models were developed for 1970–2015, a period in which a partial annual series of county-level irrigation data also were available. One of the 1970–2015 mixed models provides a basis for estimating median summer water yields at sites in southwestern Michigan using an estimated trend component, and selected basin, climatic, and irrigation factors. Re-estimation of model parameters in this mixed model with more spatially precise information on irrigation withdrawals may improve model accuracy.

Michigan↗

Encounters of aircraft with volcanic ash clouds: A compilation of known incidents, 1953-2009

Information about reported encounters of aircraft with volcanic ash clouds from 1953 through 2009 has been compiled to document the nature and scope of risks to aviation from volcanic activity. The information, gleaned from a variety of published and other sources, is presented in database and spreadsheet formats; the compilation will be updated as additional encounters occur and as new data and corrections come to light. The effects observed by flight crews and extent of aircraft damage vary greatly among incidents, and each incident in the compilation is rated according to a severity index. Of the 129 reported incidents, 94 incidents are confirmed ash encounters, with 79 of those having various degrees of airframe or engine damage; 20 are low-severity events that involve suspected ash or gas clouds; and 15 have data that are insufficient to assess severity. Twenty-six of the damaging encounters involved significant to very severe damage to engines and (or) airframes, including nine encounters with engine shutdown during flight. The average annual rate of damaging encounters since 1976, when reporting picked up, has been approximately 2 per year. Most of the damaging encounters occurred within 24 hours of the onset of ash production or at distances less than 1,000 kilometers from the source volcanoes. The compilation covers only events of relatively short duration for which aircraft were checked for damage soon thereafter; documenting instances of long-term repeated exposure to ash (or sulfate aerosols) will require further investigation. Of 38 source volcanoes, 8 have caused 5 or more encounters, of which the majority were damaging: Augustine (United States), Chaiten (Chile), Mount St. Helens (United States), Pacaya (Guatemala), Pinatubo (Philippines), Redoubt (United States), Sakura-jima (Japan), and Soufriere Hills (Montserrat, Lesser Antilles, United Kingdom). Aircraft have been damaged by eruptions ranging from small, recurring episodes to very large, infrequent events. Moderate-size (Volcanic Explosivity Index 3) eruptions are responsible for nearly half of the damaging encounters. Vigilance is required during the early phases of eruptive activity when data about ash emission may be the most limited and warning capabilities the most strained, yet the risk the greatest. The risk-mitigation strategy for minimizing damaging encounters continues to rely on the combination of real-time volcano monitoring and rapid eruption reporting, detection and tracking of ash clouds in the atmosphere using satellite-based sensors, dispersion modeling to forecast expected ash-cloud movement, and global dissemination of specialized warning messages. To obtain the entire Data Series 545 report, download the text file and appendixes 1-4, which are available as separate files. Click on the links at right. Please Send Updates We hope that publication of this compilation will encourage more reporting of encounters by the aviation industry and civil aviation authorities. We actively seek corrections and additions to the information presented here. Persons who have corrections or additional data pertaining to incidents already in the database or who have data about previously unreported incidents are urged to contact the authors.

Data Series↗

Evaluating management alternatives for Wyoming elk feedgrounds in consideration of chronic wasting disease

Executive Summary The authors used decision and modeling analyses to evaluate management alternatives for a decision on whether to permit Cervus canadensis (elk) feeding on two sites on Bridger-Teton National Forest, Dell Creek and Forest Park. Supplemental feeding of elk could increase the transmission of chronic wasting disease (CWD) locally and disease spread regionally, potentially impacting elk populations over time with wider implications for Odocoileus hemionus (mule deer) and Odocoileus virginianus (white-tailed deer) populations and hunting, tourism, and regional revenue. Supplemental feeding is thought to improve overwinter elk survival and reduce the commingling of elk with cattle during months when brucellosis transmission risk is highest. We worked with the U.S. Department of Agriculture Forest Service to identify their fundamental objectives and associated performance metrics related to this feedground decision. We then developed disease and habitat selection models to quantify the effect of four management alternatives on select performance metrics. The four alternatives were to continue to permit feeding, phaseout permits to feed in three years, permit feeding on an emergency basis, or stop permitting feeding. In this report, we present methods and summarized results on disease and habitat selection models and summaries of other performance metrics analyzed by BIO-WEST, Inc. and Cirrus Ecological Solutions as part of an Environmental Impact Statement. Data from Wyoming Game and Fish Department (WGFD) supported the assumption that supplemental elk feeding allows for larger elk populations in a region. We documented that herd units (HU) without feedgrounds had 23 percent lower densities of elk per area of winter range when compared against HUs with feedgrounds, after accounting for differences in sightability of elk during counts on and off feedgrounds. Thus, throughout our analyses, we assumed feedground closures would reduce elk carrying capacity resulting in an average decline of previously fed elk population segments by 23 percent (5th and 95th percentiles = [11 percent, 35 percent]) by year 20. Most of that decline occurred within the first few years after a feedground ceases to operate. We used a panel of CWD experts to help estimate CWD trans-mission in fed and unfed elk population segments. In aggregate, the expert panel estimated that median values of direct and indirect transmission of CWD are expected to be 1.9 and 4 times higher, respectively, in fed elk populations compared to unfed elk. We used these disease transmission estimates in combination with local elk demographic rates and carrying capacity estimates to project disease and population dynamics. In year 20, we predicted CWD prevalence would increase to 42 percent (5th and 95th percentiles = [29 percent, 55 percent]), and 13 percent (5th and 95th percentiles = [4 percent, 26 percent]) on average for fed and unfed elk population segments, respectively, given a starting prevalence of 1.6 percent. The prevalence estimates for the unfed elk population segments are in the range of previous observations of CWD in elk in the western United States. The average CWD prevalence from 2016 to 2018 in the unfed elk population of Wind Cave National Park in South Dakota was 18 percent overall but up to 30 percent in some regions (Sargeant and others, 2021). Meanwhile, CWD prevalence in the Iron Mountain and Laramie Peak elk herds in Wyoming from 2016 to 2018 was 14 percent and 7 percent, respectively, despite being present since at least 2002 (Wyoming Game and Fish Department, 2020b). From 2016 to 2020, elk that were fed at Dell Creek and Forest Park constituted on average 12–20 percent of the total elk on their respective HUs. As a result, the differences between management alternatives are modest when considering the closure of only one feedground on a HU. The no feeding alternative for Forest Park resulted in a CWD prevalence of 17 percent (SD = 7 percent) in the Afton HU compared to 20 percent (SD = 7 percent) with continued feeding by year 20. In the Upper Green River HU, no feeding on Dell Creek resulted in a CWD prevalence of 27 percent (SD = 6 percent) compared to 30 percent (SD = 5 percent) with continued feeding. In terms of disease-associated mortality, we predicted the closure of Forest Park and Dell Creek feedgrounds would reduce the total number of CWD mortalities by 9 percent in the Upper Green River HU and 26 percent in the Afton HU during the 20-year timespan. Our spatial analyses predicted that management alternative effects vary by HU as a function of private property and other wildlife winter ranges proximity relative to feedground location. The predicted number of elk abortions on private land, as a proxy for brucellosis risk to cattle, may increase by 8–21 percent in the absence of feeding at Dell Creek and Forest Park. Eight feedgrounds are located on Bridger-Teton National Forest, all of which have permits that have expired or will expire prior to 2028. In addition, WGFD could change their management of feedgrounds given new information; therefore, we also assessed the cumulative effects of continued feeding, phaseout, and no feeding management alternatives across five HUs south of Jackson, Wyoming (Afton HU, Fall Creek HU, Piney HU, Pinedale HU, and Upper Green River HU). These five HUs ranged from about 41 to 85 percent of the elk herd using feedgrounds, which corresponded to a CWD prevalence at year 20 of 23–34 percent if all feedgrounds in those five HUs remained open relative to 12 to 14 percent if all feedgrounds were closed. We predicted feedground closures may result in immediate reductions in population size relative to alternatives that continue feeding (for example, continued feeding and emergency feeding alternatives); however, over longer periods of time, CWD-associated mortality leads to larger population reductions. The no feeding alternative resulted in higher elk population sizes compared to the continued feeding alternative after about 10 years of implementation. Delayed action under a phaseout alternative resulted in increasing the CWD prevalence to 20 percent relative to 12 to 14 percent, on average, without feeding on HUs with a large population of fed elk such as the Upper Green River HU. Summarizing our cumulative results across all five of the analyzed HUs, we predicted continued feeding will lead to fewer elk by year 20 (mean = 8,300, standard deviation [SD] = 740) compared to no feeding at U.S. Department of Agri-culture Forest Service sites (10,700, SD = 890). The closure of all feedgrounds was projected to result in the largest elk populations at year 20 (12,500, SD = 980). No feeding at all sites also resulted in the largest cumulative harvest of 57,700 (SD = 2,600) compared to 51,100 (SD = 3,800) for continued feeding at all current feedground sites on the five HUs. Continued feeding also resulted in the lowest brucellosis costs to producers ($194,600, SD = $11,500) compared to no feeding on all feedgrounds ($243,000, SD = $13,700). Assuming moderate reductions in hunter interest because of increasing CWD prevalence in elk, we predicted that no feeding resulted in regional revenues generated by hunting activities of $190 million (SD = $10 million) compared to $173 million (SD = $10 million) for continued feeding over the 20-year timeframe. Recent CWD detections in mule deer and elk in Grand Teton National Park has elevated the importance of the cur-rent decision on whether, and how, to permit elk feeding on Dell Creek and Forest Park and the management of the other feedgrounds. Aggressive male harvest has slowed, but not stopped, the increasing prevalence of CWD in mule deer (Conner and others, 2021). It is unclear whether harvest management can be an effective tool to slow the spread of CWD in elk. There are also no effective treatments or vaccines for CWD, and it is unlikely that any will be developed that can be easily deployed in the near future. Thus, reducing artificial aggregations is one of the few management approaches suggested by the Western Association of Fish and Wildlife Agencies (Almberg and others, 2017). Future surveillance and monitoring can be designed to resolve uncertainties that can improve future decision-making. If feedgrounds close, research could quantify elk population reductions in the absence of feeding, the redistribution of fed elk to other places, or the consequences of elk movement on private property. If feedgrounds remain open, research could assess how rapidly CWD spreads in artificial aggregations of elk; however, surveillance programs would need to be designed with sufficient power to detect initial changes of CWD prevalence. Delaying action on feedground management was projected to be costly. Results of the phaseout alternative relative to the no feeding alternative suggested a 3-year delay was enough for substantial long-term changes in CWD prevalence. The long-term persistence of infectious CWD prions in the environment suggests that feedground management decisions may have long-lasting consequences. Our results indicated tradeoffs in the ability of a management agency to achieve all their objectives, and all management alternatives resulted in significant reductions in elk population size. This report contains the foundational elements for formal decision analysis methods, which can be implemented to help decision makers transparently evaluate the consequences of decision alternatives and identify the set of actions that best achieve agency and stakeholder priorities.

Wyoming↗

Users manual for an expert system (HSPEXP) for calibration of the hydrological simulation program; Fortran

Expert system software was developed to assist less experienced modelers with calibration of a watershed model and to facilitate the interaction between the modeler and the modeling process not provided by mathematical optimization. A prototype was developed with artificial intelligence software tools, a knowledge engineer, and two domain experts. The manual procedures used by the domain experts were identified and the prototype was then coded by the knowledge engineer. The expert system consists of a set of hierarchical rules designed to guide the calibration of the model through a systematic evaluation of model parameters. When the prototype was completed and tested, it was rewritten for portability and operational use and was named HSPEXP. The watershed model Hydrological Simulation Program--Fortran (HSPF) is used in the expert system. This report is the users manual for HSPEXP and contains a discussion of the concepts and detailed steps and examples for using the software. The system has been tested on watersheds in the States of Washington and Maryland, and the system correctly identified the model parameters to be adjusted and the adjustments led to improved calibration.

Water-Resources Investigations Report↗

Multispecies modeling for adaptive management of horseshoe crabs and red knots in the Delaware Bay

Adaptive management requires that predictive models be explicit and transparent to improve decisions by comparing management actions, directing further research and monitoring, and facilitating learning. The rufa subspecies of red knots ( Calidris canutus rufa ), which has recently exhibited steep population declines, relies on horseshoe crab ( Limulus polyphemus ) eggs as their primary food source during stopover in Delaware Bay during spring migration. We present a model with two different parameterizations for use in the adaptive management of horseshoe crab harvests in the Delaware Bay that links red knot mass gain, annual survival, and fecundity to horseshoe crab dynamics. The models reflect prevailing hypotheses regarding ecological links between these two species. When reported crab harvest from 1998 to 2008 was applied, projections corresponded to the observed red knot population abundances depending on strengths of the demographic relationship between these species. We compared different simulated horseshoe crab harvest strategies to evaluate whether, given this model, horseshoe crab harvest management can affect red knot conservation and found that restricting harvest can benefit red knot populations. Our model is the first to explicitly and quantitatively link these two species and will be used within an adaptive management framework to manage the Delaware Bay system and learn more about the specific nature of the linkage between the two species.

Delaware, New Jersey↗

N-mix for fish: estimating riverine salmonid habitat selection via N-mixture models

Models that formulate mathematical linkages between fish use and habitat characteristics are applied for many purposes. For riverine fish, these linkages are often cast as resource selection functions with variables including depth and velocity of water and distance to nearest cover. Ecologists are now recognizing the role that detection plays in observing organisms, and failure to account for imperfect detection can lead to spurious inference. Herein, we present a flexible N-mixture model to associate habitat characteristics with the abundance of riverine salmonids that simultaneously estimates detection probability. Our formulation has the added benefits of accounting for demographics variation and can generate probabilistic statements regarding intensity of habitat use. In addition to the conceptual benefits, model application to data from the Trinity River, California, yields interesting results. Detection was estimated to vary among surveyors, but there was little spatial or temporal variation. Additionally, a weaker effect of water depth on resource selection is estimated than that reported by previous studies not accounting for detection probability. N-mixture models show great promise for applications to riverine resource selection.

Canadian Journal of Fisheries and Aquatic Sciences↗

Automated construction of Streamflow-Routing networks for MODFLOW—Application in the Mississippi Embayment region

In humid regions with dense stream networks, surface water exerts a fundamental control on the water levels and flow directions of shallow groundwater. Understanding interactions between groundwater and surface water is critical for managing groundwater resources and groundwater-dependent ecosystems. Representing streams in groundwater models has historically been arduous and error prone. In recent years, however, all the information needed to numerically describe stream boundary conditions for a model area has become readily available online, as have robust open-source software tools for translating that information to a model grid. The SFRmaker Python package leverages geospatial capabilities in the scientific Python ecosystem to robustly automate the production of input to the Streamflow-Routing (SFR) Package of MODFLOW from the National Hydrography Dataset Plus or other hydrography data. This report documents an application of SFRmaker to automate production of SFR Package input for groundwater models within the Mississippi Embayment Regional Aquifer Study area. SFR Package input was developed in three steps: (1) preprocessing to develop a single set of grid-independent flowlines from National Hydrography Dataset Plus version 2 data; (2) setting up the SFR package from the preprocessed flowlines, and (3) correcting streambed top elevations after an initial model run. Separating the hydrography preprocessing from the construction of SFR Package input was advantageous in that it minimized the need to repeat computationally expensive geoprocessing (thereby speeding model construction) and also allowed for the curation of a single set of grid-independent SFR input data that can be used for any MODFLOW model within the study area.

Mississippi Embayment Region↗