USGS Science⌕ Search

SEARCH · USGS Science

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 433 records · Page 24Linked to original sources

An updated conceptual model of Delta Smelt biology: Our evolving understanding of an estuarine fish

The main purpose of this report is to provide an up-to-date assessment and conceptual model of factors affecting Delta Smelt (Hypomesus transpacificus) throughout its primarily annual life cycle and to demonstrate how this conceptual model can be used for scientific and management purposes. The Delta Smelt is a small estuarine fish that only occurs in the San Francisco Estuary. Once abundant, it is now rare and has been protected under the federal and California Endangered Species Acts since 1993. The Delta Smelt listing was related to a step decline in the early 1980s; however, population abundance decreased even further with the onset of the “pelagic organism decline” (POD) around 2002. A substantial, albeit short-lived, increase in abundance of all life stages in 2011 showed that the Delta Smelt population can still rebound when conditions are favorable for spawning, growth, and survival. In this report, we update previous conceptual models for Delta Smelt to reflect new data and information since the release of the last synthesis report about the POD by the Interagency Ecological Program for the San Francisco Estuary (IEP) in 2010. Specific objectives include: 1. Provide decision makers with a practical tool for evaluating difficult trade-offs associated with management and policy decisions. 2. Provide scientists with a framework from which they can formulate and evaluate hypotheses using qualitative or quantitative models. 3. Provide the general public with a new way of learning about Delta Smelt and their habitat. Our updated conceptual model describes the habitat conditions and ecosystem drivers affecting each Delta Smelt life stage, across seasons and how the seasonal effects contribute to the annual success of the species. The conceptual model consists of two nested and linked levels of increasing specificity. The general life cycle conceptual model for four Delta Smelt life stages (adults, eggs and larvae, juveniles, and subadults) includes stationary ecosystem components and dynamic environmental drivers, habitat attributes, and Delta Smelt responses. The more detailed life stage transition conceptual models for each of the four Delta Smelt life stages describe relationships between environmental drivers, key habitat attributes, and the responses of Delta Smelt to habitat attributes as they transition from one life stage to the next. Our analyses and conceptual model show that good larval recruitment is essential for setting the stage for a strong year class; however, increased growth and survival through subsequent life stages are also needed to achieve and sustain higher population abundance. We used our conceptual model to generate 16 hypotheses about the factors that may have contributed to the 2011 increase in Delta Smelt relative abundance. We then evaluated these hypotheses by comparing habitat conditions and Delta Smelt responses in the wet year 2011 to those in the prior wet year 2006 and in the drier years 2005 and 2010. Larval recruitment was similarly high in both wet years and lower in the drier antecedent years, but juvenile and adult abundance increased only in 2011. In 2005 and 2006, the population was limited by very poor survival from the larval to the juvenile life stage. We found that in 2011, Delta Smelt may have benefitted from a combination of favorable habitat conditions throughout the year, including: 1. Adults and larvae benefitted from prolonged cool spring water temperatures, high 2011 winter and spring outflows which reduced entrainment risk and possibly improved other habitat conditions, and possibly enhanced food availability in late spring. 2. Juveniles benefitted from cool water temperatures in late spring and early summer as well as from improved food availability and low levels of harmful Microcystis. 3. Subadults also benefitted from improved food availability and from favorable habitat conditions in the large, low salinity zone (salinity 1-6) located more toward Suisun Bay,

California↗

Developing and implementing the use of predictive models for estimating water quality at Great Lakes beaches

Predictive models have been used at beaches to improve the timeliness and accuracy of recreational water-quality assessments over the most common current approach to water-quality monitoring, which relies on culturing fecal-indicator bacteria such as Escherichia coli (E. coli.) . Beach-specific predictive models use environmental and water-quality variables that are easily and quickly measured as surrogates to estimate concentrations of fecal-indicator bacteria or to provide the probability that a State recreational water-quality standard will be exceeded. When predictive models are used for beach closure or advisory decisions, they are referred to as “nowcasts.” During the recreational seasons of 2010-12, the U.S. Geological Survey (USGS), in cooperation with 23 local and State agencies, worked to improve existing nowcasts at 4 beaches, validate predictive models at another 38 beaches, and collect data for predictive-model development at 7 beaches throughout the Great Lakes. This report summarizes efforts to collect data and develop predictive models by multiple agencies and to compile existing information on the beaches and beach-monitoring programs into one comprehensive report. Local agencies measured E. coli concentrations and variables expected to affect E. coli concentrations such as wave height, turbidity, water temperature, and numbers of birds at the time of sampling. In addition to these field measurements, equipment was installed by the USGS or local agencies at or near several beaches to collect water-quality and metrological measurements in near real time, including nearshore buoys, weather stations, and tributary staff gages and monitors. The USGS worked with local agencies to retrieve data from existing sources either manually or by use of tools designed specifically to compile and process data for predictive-model development. Predictive models were developed by use of linear regression and (or) partial least squares techniques for 42 beaches that had at least 2 years of data (2010-11 and sometimes earlier) and for 1 beach that had 1 year of data. For most models, software designed for model development by the U.S. Environmental Protection Agency (Virtual Beach) was used. The selected model for each beach was based on a combination of explanatory variables including, most commonly, turbidity, day of the year, change in lake level over 24 hours, wave height, wind direction and speed, and antecedent rainfall for various time periods. Forty-two predictive models were validated against data collected during an independent year (2012) and compared to the current method for assessing recreational water quality-using the previous day’s E. coli concentration (persistence model). Goals for good predictive-model performance were responses that were at least 5 percent greater than the persistence model and overall correct responses greater than or equal to 80 percent, sensitivities (percentage of exceedances of the bathing-water standard that were correctly predicted by the model) greater than or equal to 50 percent, and specificities (percentage of nonexceedances correctly predicted by the model) greater than or equal to 85 percent. Out of 42 predictive models, 24 models yielded over-all correct responses that were at least 5 percent greater than the use of the persistence model. Predictive-model responses met the performance goals more often than the persistence-model responses in terms of overall correctness (28 versus 17 models, respectively), sensitivity (17 versus 4 models), and specificity (34 versus 25 models). Gaining knowledge of each beach and the factors that affect E. coli concentrations is important for developing good predictive models. Collection of additional years of data with a wide range of environmental conditions may also help to improve future model performance. The USGS will continue to work with local agencies in 2013 and beyond to develop and validate predictive models at beaches and improve existing nowcasts, restructuring monitoring activities to accommodate future uncertainties in funding and resources.

Great Lakes↗

United States Geological Survey uranium and thorium resource assessment and exploration research program, fiscal year 1977

The U.S. Geological Survey Uranium and Thorium Resource Assessment and Exploration Research Program began in July 1974 (FY 1975). The first year was spent chiefly in redirecting an on-going effort and in staffing for the nearly fourfold increase in funding. A summary of the Fiscal Year 1976 program, including objectives, program model, and project descriptions, is published as Open-File Report 75-575. The overall objectives and program model remain essentially the same as in FY 1976. It should be noted that the beginning of the fiscal year was changed from July 1 to October 1 during the 1976 calendar year and that a transitional quarter was used to make this change. Some program changes were made during the transitional quarter, but these details are not recorded in the present summary.

Open-File Report↗

Mapping monkeypox transmission risk through time and space in the Congo Basin

Monkeypox is a major public health concern in the Congo Basin area, with changing patterns of human case occurrences reported in recent years. Whether this trend results from better surveillance and detection methods, reduced proportions of vaccinated vs. non-vaccinated human populations, or changing environmental conditions remains unclear. Our objective is to examine potential correlations between environment and transmission of monkeypox events in the Congo Basin. We created ecological niche models based on human cases reported in the Congo Basin by the World Health Organization at the end of the smallpox eradication campaign, in relation to remotely-sensed Normalized Difference Vegetation Index datasets from the same time period. These models predicted independent spatial subsets of monkeypox occurrences with high confidence; models were then projected onto parallel environmental datasets for the 2000s to create present-day monkeypox suitability maps. Recent trends in human monkeypox infection are associated with broad environmental changes across the Congo Basin. Our results demonstrate that ecological niche models provide useful tools for identification of areas suitable for transmission, even for poorly-known diseases like monkeypox.

PLoS ONE↗

Use of solicited bands and separation of hunting and natural mortality: a comment

We previously presented (Conroy 1985, Conroy et al. 1989) models for analysis of band recovery data stratified into those voluntarily reported and those obtained by solicitation, similar to models described by Pollock et al. (1994). We discuss differences between the 2 modeling approaches and suggest designs and analyses involving overlapping samples of standard and reward bands to enable adjustment for temporal and geographic variation in reporting and solicitation rates,

Journal of Wildlife Management↗

Magnetotelluric Data, Rainier Mesa/Shoshone Mountain, Nevada Test Site, Nevada

Introduction: The United States Department of Energy (DOE) and the National Nuclear Security Administration (NNSA) at their Nevada Site Office (NSO) are addressing ground-water contamination resulting from historical underground nuclear testing through the Environmental Management (EM) program and, in particular, the Underground Test Area (UGTA) project. During 2005, the U.S. Geological Survey (USGS), in cooperation with the DOE and NNSA-NSO, collected and processed data from twenty-six magnetotelluric (MT) and audio-magnetotelluric (AMT) sites at the Nevada Test Site. The 2005 data stations were located on and near Rainier Mesa and Shoshone Mountain to assist in characterizing the pre-Tertiary geology in those areas. These new stations extend the area of the hydrogeologic study previously conducted in Yucca Flat. The MT data presented in this report will help refine what is known about the character, thickness, and lateral extent of pre Tertiary confining units. Subsequent interpretation will include a three dimensional (3 D) character analysis and a two-dimensional (2 D) resistivity model. The purpose of this report is to release the MT sounding data. No interpretation of the data is included here.

Open-File Report↗

Model-based scenario planning to inform climate change adaptation in the Northern Great Plains—Final report

Public Summary We worked with managers in two focal areas to plan for the uncertain future by integrating quantitative climate change scenarios and simulation modeling into scenario planning exercises. In our central North Dakota focal area, centered on Knife River Indian Villages National Historic Site, managers are concerned about how changes in flood severity and growing conditions for native and invasive plants may affect archaeological resources and cultural landscapes associated with the Knife and Missouri Rivers. Climate projections and hydrological modeling based on those projections indicate plausible changes in spring and summer soil moisture ranging from a 7 percent decrease to a 13 percent increase and maximum winter snowpack (important for spring flooding) changes ranging from a 13 percent decrease to a 47 percent increase. Facilitated discussions among managers and scientists exploring the implications of these different climate scenarios for resource management revealed potential conflicts between protecting archeological sites and fostering riparian cottonwood forests. The discussions also indicated the need to prioritize archeological sites for excavation or protection and culturally important plant species for intensive management attention. In our southwestern South Dakota focal area, centered on Badlands National Park, managers are concerned about how changing climate will affect vegetation production, wildlife populations, and erosion of fossils, archeological artifacts, and roads. Climate scenarios explored by managers and scientists in this focal area ranged from a 13 percent decrease to a 33 percent increase in spring precipitation, which is critical to plant growth in the northern Great Plains region, and a slight decrease to a near doubling of intense rain events. Facilitated discussions in this focal area concluded that greater effort should be put into preparing for emergency protection, excavation, and preservation of exposed fossils or artifacts and revealed substantial opportunities for different agencies to learn from each other and cooperate on common management goals. Follow up quantitative simulation modeling of grassland dynamics helped quantify the degree of change expected in vegetation production under the wide range of climate scenarios and suggested that (a) low grazing rates could be adversely affecting vegetation composition in the national park and (b) understanding of the management practices needed to maintain desired vegetation conditions is incomplete.

Montana, Nebraska, North Dakota, South Dakota, Wyo↗

Enhanced and updated spatially referenced statistical assessment of dissolved-solids load sources and transport in streams of the Upper Colorado River Basin

Approximately 6.4 million tons of dissolved solids are discharged from the Upper Colorado River Basin (UCRB) to the Lower Colorado River Basin each year. This results in substantial economic damages, and tens of millions of dollars are spent annually on salinity control projects designed to reduce salinity loads in surface waters of the UCRB. Dissolved solids in surface water and groundwater have been studied extensively over the past century, and these studies have contributed to a conceptual understanding of sources and transport of dissolved solids. This conceptual understanding was incorporated into a Spatially Referenced Regressions on Watershed Attributes (SPARROW) model to examine sources and transport of dissolved solids in the UCRB. The results of this model were published in 2009. The present report documents the methods and data used to develop an updated dissolved-solids SPARROW model for the UCRB, and incorporates data defining current basin attributes not available in the previous model, including delineation of irrigated lands by irrigation type (sprinkler or flood irrigation), and calibration data from additional monitoring sites. Dissolved-solids loads estimated for 312 monitoring sites were used to calibrate the SPARROW model, which predicted loads for each of 10,789 stream reaches in the UCRB. The calibrated model provided a good fit to the calibration data as evidenced by R 2 and yield R 2 values of 0.96 and 0.73, respectively, and a root-mean-square error of 0.47. The model included seven geologic sources that have estimated dissolved-solids yields ranging from approximately 1 to 45 tons per square mile (tons/mi 2 ). Yields generated from irrigated agricultural lands are substantially greater than those from geologic sources, with sprinkler irrigated lands generating an average of approximately 150 tons/mi 2 and flood irrigated lands generating between 770 and 2,300 tons/mi 2 depending on underlying lithology. The coefficients estimated for six landscape transport characteristics that influence the delivery of dissolved solids from sources to streams, are consistent with the process understanding of dissolved-solids loading to streams in the UCRB. Dissolved-solids loads and the proportion of those loads among sources in the entire UCRB as well as in major tributaries in the basin are reported, as are loads generated from irrigated lands, rangelands, Bureau of Land Management (BLM) lands, and grazing allotments on BLM lands. Model-predicted loads also are compared with load estimates from 1957 and 1991 at selected locations in three divisions of the UCRB. At the basin scale, the model estimates that 32 percent of the dissolved-solids loads are from irrigated agricultural land sources that compose less than 2 percent of the land area in the UCRB. This estimate is less than previously reported estimates of 40 to 45 percent of basin-scale dissolved-solids loads from irrigated agricultural land sources. This discrepancy could be a result of the implementation of salinity control projects in the basin. Notably, results indicate that the conversion of flood irrigated agricultural lands to sprinkler irrigated agricultural lands is a likely process contributing to the temporal decrease in dissolved-solids loads from irrigated lands.

Colorado River Basin↗

Regional rainfall-runoff relations for simulation of streamflow for watersheds in Lake County, Illinois

Rainfall and streamflow data collected in Lake County, Ill., from March 1990 through September 1993 were used to (1) calibrate a rainfall-runoff model for an area encompassing three watersheds (individual areas of 17.2, 35.7, and 37.0 mi 2 (square miles) and (2) verify the regional model parameter set obtained from the calibration by applying the parameter set to rainfall-runoff models for an additional small (6.3 mi 2 ) watershed and a large (59.6 mi 2 ) watershed. In addition, rainfall and streamflow data collected from April 1991 through September 1993 were used to calibrate the rainfall-runoff model for three single land-use watersheds (38.2-305 acres), called hydrologic response units (HRU's). Significant differences were found between the best parameters used in the HRU models and in the larger watershed models. The main channels in the HRU's are intermittent streams; thus, the parameters in the HRU models were selected such that a fluctuating water table could be simulated; runoff from the larger watersheds is not as sensitive to the effects of a fluctuating water table. Classification of land cover into two pervious subareas (forest and grass) and one impervious subarea (including parking lots, streets, and rooftops, among others) was sufficient to simulate the rainfall-runoff relations for all watersheds accurately. The model parameters presented in this report, which were refined through regional calibration and verified for watersheds not considered in the calibration, allow simulation of runoff in watersheds in Lake County, Ill., with approximately 93-percent accuracy in the total water balance, an average absolute error in the annual- flow estimates of 10.9 percent (and an error rarely exceeding 25 percent for annual flow), and monthly water balances with correlation coefficients of 93 percent and coefficients of model-fit efficiency of 86 percent. The models closely reproduced the partial-duration series of runoff and storm-runoff frequencies for the modeled watersheds.

Illinois↗

Effects of hypothetical management scenarios on simulated water temperatures in the Tualatin River, Oregon

In 1994, the U.S. Geological Survey (USGS) and the Unified Sewerage Agency of Washington County, Oregon (USA) began a cooperative study to better understand water-temperature variations in the Tualatin River and to assess mitigative water-management solutions. Continuous water-temperature data were collected at locations along the main stem of the river and along the major tributaries during the lowflow periods of 1994 and 1995. The 1994 data were used to develop and calibrate flow and water-temperature models characterizing conditions in the main stem. The models were used to simulate 10 hypotheti3 cal water-management scenarios, which would enable water managers to understand the effects of various human activities on water temperatures. Modeling results from the study are presented in Risley (1997); the data collected are presented in Risley and Doyle (1997). This report presents the water-temperature model simulation results of 16 additional hypothetical water-management scenarios using the 1994 and 1995 data. The additional modeling was funded by the USGS and the USA under a cooperative agreement. For a comprehensive description of the water-temperature models and their underlying assumptions, refer to Risley (1997).

Oregon↗

Estimating the effects of wetland conservation practices in croplands: Approaches for modeling in CEAP–Cropland Assessment

Quantifying the current and potential benefits of conservation practices can be a valuable tool for encouraging greater practice adoption on agricultural lands. A goal of the CEAP-Cropland Assessment is to estimate the environmental effects of conservation practices that reduce losses (exports) of soil, nutrients, and pesticides from farmlands to streams and rivers. The assessment approach combines empirical data on reported cropland practices with simulation modeling that compares field-level exports for scenarios “with practices” and “without practices.” Conserved, restored, and created wetlands collectively represent conservation practices that can influence sediment and nutrient exports from croplands. However, modeling the role of wetlands within croplands presents some challenges, including the potential for negative impacts of sediment and nutrient inputs on wetland functions. This Science Note outlines some preliminary solutions for incorporating wetlands and wetland practices into the CEAP-Cropland modeling framework. First, modeling the effects of wetland practices requires identifying wetland hydrogeomorphic type and accounting for the condition of both the wetland and an adjacent upland zone. Second, modeling is facilitated by classifying wetland-related practices into two functional categories (wetland and upland buffer). Third, simulating practice effects requires alternative field configurations to account for hydrological differences among wetland types. These ideas are illustrated for two contrasting wetland types (riparian and depressional).

CEAP-Wetlands Science Note↗

Fault slip and seismic moment of the 1700 Cascadia earthquake inferred from Japanese tsunami descriptions

The 1700 Cascadia earthquake attained moment magnitude 9 according to new estimates based on effects of its tsunami in Japan, computed coseismic seafloor deformation for hypothetical ruptures in Cascadia, and tsunami modeling in the Pacific Ocean. Reports of damage and flooding show that the 1700 Casscadia tsunami reached 1-5 m heights at seven shoreline sites in Japan. Three sets of estimated heights express uncertainty about location and depth of reported flooding, landward decline in tsunami heights from shorelines, and post-1700 land-level changes. We compare each set with tsunami heights computed from six Cascadia sources. Each source is vertical seafloor displacement calculated with a three-dimensional elastic dislocation model, for three sources the rupture extends the 1100 km length of the subduction zone and differs in width and shallow dip; for the other sources, ruptures of ordinary width extend 360-670 km. To compute tsunami waveforms, we use a linear long-wave approximation with a finite difference method, and we employ modern bathymetry with nearshore grid spacing as small as 0.4 km. The various combinations of Japanese tsunami heights and Cascadia sources give seismic moment of 1-9 ?? 1022 N m, equivalent to moment magnitude 8.7-9.2. This range excludes several unquantified uncertainties. The most likely earthquake, of moment magnitude 9.0, has 19 m of coseismic slip on an offshore, full-slip zone 1100 km long with linearly decreasing slip on a downdip partial-slip zone. The shorter rupture models require up to 40 m offshore slip and predict land-level changes inconsistent with coastal paleoseismological evidence. Copyright 2003 by the American Geophysical Union.

Journal of Geophysical Research B: Solid Earth↗

Report on the final completion of the Unified Lunar Control Network 2005 and Lunar Topographic Model

In order to highlight this project to the extraterrestrial mapping community, we repeat here our earlier abstract [1], with a corrected Figure 2. A report describing the Unified Lunar Control Network 2005 and the files associated with that network is now available as an on-line USGS Open-File Report [2] at the location http://pubs.usgs.gov/of/2006/1367/. A “Readme” file describes the available files, including the report text, the original photogrammetric solution input and output files, derived files such as information on all the control point positions, the expected vertical precision (EVP) of the points and the solution residuals, and a directory containing global DEMs derived from the point positions. While a paper providing further details about the solution is in preparation, we provide some useful additional information about the solution and possible uses for it here.

Conference Paper↗

Hydrogeology of well-field areas near Tampa, Florida, phase 2 — Development and documentation of a quasi-three-dimensional finite-difference model for simulation of steady-state ground-water flow

This report describes a quasi-three-dimensional finite-difference model for simulation of steady-state ground-water flow in the Floridan aquifer over a 932-square-mile area that contains 10 municipal well fields. The over-lying surficial aquifer contains a water table and is coupled to the Floridan aquifer by leakage term that represents flow through a confining layer separating the two aquifers. Under the steady-state condition, all storage terms are set to zero. Use of the head-controlled flux condition allows simulated head and flow changes to occur in the Floridan aquifer at the model boundaries. Procedures used to calibrate the model, test its sensitivity to input-parameter errors, and validate its accuracy for predictive purposes are described. Also included are attachments that describe setting up and running the model. Example model-interrogation runs show anticipated drawdowns under high, average, and low recharge conditions with 10 well fields pumping simultaneously at the maximum annual permitted rates totaling 186.9 million gallons per day.

Florida↗

Modeled inundation limits of potential lahars from Mount Adams in the White Salmon River Valley, Washington

Lahars large enough to reach populated areas are a hazard at Mount Adams, a massive volcano in the southern Cascade Range of Washington State (fig. 1). It is considered to be still active and has the potential to erupt again. By definition, lahars are gravity-driven flows of water-saturated mixtures of mud and rock (plus or minus ice, wood, and other debris), which originate from volcanoes and have a variety of potential triggering mechanisms (Vallance, 2000; Vallance and Iverson, 2015). Flowing mixtures can range in fluid consistency from something like a milkshake to something more like wet concrete, and they behave like flash floods, in that they can appear suddenly in river channels with little warning and commonly have boulder- or log-choked flow fronts. Lahars are hazardous because they can flow rapidly in confined valleys (commonly 20–35 miles per hour [mph] or 9–16 meters per second [m/s]), can travel more than 100 miles (mi) (161 kilometers [km]) from a source volcano, and can move with incredible destructive force, carrying multi-ton boulders and logs that can act as battering rams (Pierson, 1998). The biggest threats from lahars to downstream communities are present during eruptive activity, and impacts to communities can be dire. For example, a very large eruption-triggered lahar in Colombia in 1985 surprised and killed more than 20,000 people in a large town located about 45 mi (72 km) downstream and out of sight of the volcano that produced it (Pierson and others, 1990). Mount Adams, one of the largest volcanoes in the Cascade Range, is a composite stratocone composed primarily of andesite lava flows. It has been the most continuously active volcano within the 480-mi 2 Mount Adams volcanic field—a region covering parts of Klickitat, Skamania, Yakima, andLewis Counties and part of the Yakama Nation Reservation in Washington State (Hildreth and Fierstein,1995, 1997). About 500,000 years in age, Mount Adams reached its present size by about 15,000 years ago, primarily through the episodic effusion of lava flows; it has not had a history of major explosive eruptions like Mount St. Helens, its neighbor to the west. Timing of the most recent eruptive activity (recorded by four thin tephra layers) is on the order of 1,000 years ago; the tephras are bracketed by 2,500-year-old and 500-year-old ash layers from Mount St. Helens (Hildreth and Fierstein, 1995, 1997). Mount Adams currently shows no signs of renewed unrest. Eruptive history does not tell us everything we need to know about hazards at Mount Adams, however, which are fully addressed in the volcano hazard assessment for Mount Adams (W.E. Scott and others, 1995). This volcano has had a long-active hydrothermal system that circulated acidic hydrothermal fluids, formed by the solution of volcanic gases in heated groundwater, through fractures and permeable zones into upper parts of the volcanic cone. Acid sulfate leaching of rocks in the summit area may still be occurring, but chemical and thermal evidence suggests that the main hydrothermal system is no longer active at Mount Adams (Nathenson and Mariner, 2013). However, these rock-weakening chemical reactions have operated long enough to change about 0.4 cubic miles (mi 3 ) (1.7 cubic kilometers [km 3 ]) of the hard lava rock in the volcano’s upper cone to a much weaker clay-rich rock, thus significantly reducing rock strength and thereby slope stability in parts of the cone (Finn and others, 2007). The two largest previous lahars from Mount Adams were triggered by landslides of hydrothermally altered rock from the upper southwestern flank of the cone, and any future large lahars are likely to be triggered by the same mechanism. Mount Rainier also has had extensive hydrothermal alteration of rock in its upper edifice, and it also has a history of large landslides that transform into lahars (K.M. Scott and others, 1995; Vallance and Scott, 1997; Reid and others, 2001). The spatial depiction of modeled lahar inundation zones accompanying this report, shown in two different map perspectives, is intended to augment (not replace) the existing hazard maps for Mount Adams (W.E. Scott and others, 1995; Vallance, 1999). The maps in this report show potential areas of inundation by lahars of different initial volumes, which are determined by a computer model, LAHARZ (Iverson and others, 1998; Schilling, 1998). One map sheet presents LAHARZ-determined inundation areas on a normal plan-view shaded-relief map of the study area; the other gives an oblique perspective of the landscape with raised topography, as if one were viewing the landscape at an angle from an aircraft (Jenny and Patterson, 2007). LAHARZ was developed after the original hazard maps (based only on mapping of geologic deposits) were made. Predicted inundation zones on these maps provide an alternative approach to estimation of areas that could be inundated as lahars of different volumes pass through the valley. However, there is considerable uncertainty in the exact location of the hazard-zone boundaries shown on these maps, as well as on earlier maps.

Washington↗

A preliminary deposit model for lithium-cesium-tantalum (LCT) pegmatites

This report is part of an effort by the U.S. Geological Survey to update existing mineral deposit models and to develop new ones. We emphasize practical aspects of pegmatite geology that might directly or indirectly help in exploration for lithium-cesium-tantalum (LCT) pegmatites, or for assessing regions for pegmatite-related mineral resource potential. These deposits are an important link in the world’s supply chain of rare and strategic elements, accounting for about one-third of world lithium production, most of the tantalum, and all of the cesium.

Open-File Report↗

Comparison of two regression-based approaches for determining nutrient and sediment fluxes and trends in the Chesapeake Bay watershed

Nutrient and sediment fluxes and changes in fluxes over time are key indicators that water resource managers can use to assess the progress being made in improving the structure and function of the Chesapeake Bay ecosystem. The U.S. Geological Survey collects annual nutrient (nitrogen and phosphorus) and sediment flux data and computes trends that describe the extent to which water-quality conditions are changing within the major Chesapeake Bay tributaries. Two regression-based approaches were compared for estimating annual nutrient and sediment fluxes and for characterizing how these annual fluxes are changing over time. The two regression models compared are the traditionally used ESTIMATOR and the newly developed Weighted Regression on Time, Discharge, and Season (WRTDS). The model comparison focused on answering three questions: (1) What are the differences between the functional form and construction of each model? (2) Which model produces estimates of flux with the greatest accuracy and least amount of bias? (3) How different would the historical estimates of annual flux be if WRTDS had been used instead of ESTIMATOR? One additional point of comparison between the two models is how each model determines trends in annual flux once the year-to-year variations in discharge have been determined. All comparisons were made using total nitrogen, nitrate, total phosphorus, orthophosphorus, and suspended-sediment concentration data collected at the nine U.S. Geological Survey River Input Monitoring stations located on the Susquehanna, Potomac, James, Rappahannock, Appomattox, Pamunkey, Mattaponi, Patuxent, and Choptank Rivers in the Chesapeake Bay watershed. Two model characteristics that uniquely distinguish ESTIMATOR and WRTDS are the fundamental model form and the determination of model coefficients. ESTIMATOR and WRTDS both predict water-quality constituent concentration by developing a linear relation between the natural logarithm of observed constituent concentration and three explanatory variables—the natural log of discharge, time, and season. ESTIMATOR uses two additional explanatory variables—the square of the log of discharge and time-squared. Both models determine coefficients for variables for a series of estimation windows. ESTIMATOR establishes variable coefficients for a series of 9-year moving windows; all observed constituent concentration data within the 9-year window are used to establish each coefficient. Conversely, WRTDS establishes variable coefficients for each combination of discharge and time using only observed concentration data that are similar in time, season, and discharge to the day being estimated. As a result of these distinguishing characteristics, ESTIMATOR reproduces concentration-discharge relations that are closely approximated by a quadratic or linear function with respect to both the log of discharge and time. Conversely, the linear model form of WRTDS coupled with extensive model windowing for each combination of discharge and time allows WRTDS to reproduce observed concentration-discharge relations that are more sinuous in form. Another distinction between ESTIMATOR and WRTDS is the reporting of uncertainty associated with the model estimates of flux and trend. ESTIMATOR quantifies the standard error of prediction associated with the determination of flux and trends. The standard error of prediction enables the determination of the 95-percent confidence intervals for flux and trend as well as the ability to test whether the reported trend is significantly different from zero (where zero equals no trend). Conversely, WRTDS is unable to propagate error through the many (over 5,000) models for unique combinations of flow and time to determine a total standard error. As a result, WRTDS flux estimates are not reported with confidence intervals and a level of significance is not determined for flow-normalized fluxes. The differences between ESTIMATOR and WRTDS, with regard to model form and determination of model coefficients, have an influence on the determination of nutrient and sediment fluxes and associated changes in flux over time as a result of management activities. The comparison between the model estimates of flux and trend was made for combinations of five water-quality constituents at nine River Input Monitoring stations. The major findings with regard to nutrient and sediment fluxes are as follows: (1)WRTDS produced estimates of flux for all combinations that were more accurate, based on reduction in root mean squared error, than flux estimates from ESTIMATOR; (2) for 67 percent of the combinations, WRTDS and ESTIMATOR both produced estimates of flux that were minimally biased compared to observed fluxes(flux bias = tendency to over or underpredict flux observations); however, for 33 percent of the combinations, WRTDS produced estimates of flux that were considerably less biased (by at least 10 percent) than flux estimates from ESTIMATOR; (3) the average percent difference in annual fluxes generated by ESTIMATOR and WRTDS was less than 10 percent at 80 percent of the combinations; and (4) the greatest differences related to flux bias and annual fluxes all occurred for combinations where the pattern in observed concentration-discharge relation was sinuous (two points of inflection) rather than linear or quadratic (zero or one point of inflection). The major findings with regard to trends are as follows: (1) both models produce water-quality trends that have factored in the year-to-year variations in flow; (2) trends in water-quality condition are represented by ESTIMATOR as a trend in flow-adjusted concentration and by WRTDS as a flow normalized flux; (3) for 67 percent of the combinations with trend estimates, the WRTDS trends in flow-normalized flux are in the same direction and magnitude to the ESTIMATOR trends in flow-adjusted concentration, and at the remaining 33 percent the differences in trend magnitude and direction are related to fundamental differences between concentration and flux; and (4) the majority (85 percent) of the total nitrogen, nitrate, and orthophosphorus combinations exhibited long-term (1985 to 2010) trends in WRTDS flow-normalized flux that indicate improvement or reduction in associated flux and the majority (83 percent) of the total phosphorus (from 1985 to 2010) and suspended sediment (from 2001 to 2010) combinations exhibited trends in WRTDS flow-normalized flux that indicate degradation or increases in the flux delivered.

Delaware, Maryland, New York, Pennsylvania, Virgin↗

Model selection for the North American Breeding Bird Survey

The North American Breeding Bird Survey (BBS) provides data that can be used in complex, multiscale analyses of population change, while controlling for scale‐specific nuisance factors. Many alternative models can be fit to the data, but most model selection procedures are not appropriate for hierarchical models. Leave‐one‐out cross‐validation (LOOCV), in which relative model fit is assessed by omitting an observation and assessing the prediction of a model fit using the remainder of the data, provides a reasonable approach for assessing models, but is time consuming and not feasible to apply for all observations in large data sets. We report the first large‐scale formal model selection for BBS data, applying LOOCV to stratified random samples of observations from BBS data. Our results are for 548 species of North American birds, comparing the fit of four alternative models that differ in year effect structures and in descriptions of extra‐Poisson overdispersion. We use a hierarchical model among species to evaluate posterior probabilities that models are best for individual species. Models in which differences in year effects are conditionally independent (D models) were generally favored over models in which year effects are modeled by a slope parameter and a random year effect (S models), and models in which extra‐Poisson overdispersion effects are independent and t ‐distributed (H models) tended to be favored over models where overdispersion was independent and normally distributed. Our conclusions lead us to recommend a change from the conventional S model to D and H models for the vast majority of species (544/548). Comparison of estimated population trends based on the favored model relative to the S model currently used for BBS summaries indicates no consistent differences in estimated trends. Of the 18 species that showed large differences in estimated trends between models, estimated trends from the default S model were more extreme, reflecting the influence of the slope parameter in that model for species that are undergoing large population changes. WAIC, a computationally simpler alternative to LOOCV, does not appear to be a reliable alternative to LOOCV.

Ecological Applications↗