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 1,261 records · Page 70Linked to original sources
USGS2018

Golden Eagle Monitoring Plan for the Desert Renewable Energy Conservation Plan

This report describes options for monitoring the status and population trends of the golden eagle (Aquila chrysaetos) within the Desert Renewable Energy Conservation Plan (DRECP) area of Southern California in maintaining stable or increasing population in the planning area. The report profiles the ecology of golden eagles in the region and provides a range of potential sampling options to address monitoring needs and objectives. This approach also focused on links between changes in human land-use, golden eagle nesting and foraging habitat conditions, and population dynamics. The report outlines how monitoring data from demographic, prey, and habitat studies were used to develop a predictive demographic model for golden eagles in the DRECP area. Results from the model simulations suggest increases in renewable energy development could have negative consequences for population trajectories. Results also suggest site-specific conservation actions could reduce the magnitude of negative impacts to the local population of eagles. A monitoring framework is proposed including: (1) annual assessments of site-occupancy and reproduction by territorial pairs of golden eagles (including rates at which sites become colonized or vacated over time); (2) estimates of survival, movements, and intensity of use of landscapes by breeding and non-breeding golden eagles; (3) periodic (conducted every two to four years) assessments of nesting and foraging habitats, prey populations, and associations with land-use and management activities; and (4) updating the predictive demographic model with new information obtained on eagles and associated population stressors. The results of this research were published in the Journal of Rapture Research, Wiens, David,Inman, Rich D., Esque, Todd C., Longshore, Kathleen M. and Nussear, Kenneth (2017). Spatial Demographic Models to Inform Conservation Planning of Golden Eagles in Renewable Energy Landscapes. 51(3):234-257.

Report↗
USGS2003

Use of boundary fluxes when simulating solute transport with the MODFLOW ground-water transport process

This report describes modifications to a U.S. Geological Survey (USGS) threedimensional solute-transport model (MODFLOWGWT), which is incorporated into the USGS MODFLOW ground-water model as the Ground- Water Transport (GWT) Process. The modifications improve the capability of MODFLOW-GWT to accurately simulate solute transport in simulations that represent a nonzero flux across an aquifer boundary. In such situations, the new Boundary Flux Package (BFLX) will allow the user flexibility to assign the flux to specific cell faces, although that flexibility is limited for certain types of fluxes (such as recharge and evapotranspiration, which can only be assigned to the top face if either is to be represented as a boundary flux). The approach is consistent with that used in the MODPATH model. The application of the BFLX Package was illustrated using a test case in which the Lake Package was active. The results using the BFLX Package showed noticeably higher magnitudes of velocity in the cells adjacent to the lake than previous results without the BFLX Package. Consequently, solute was transported slightly faster through the lake-aquifer system when the BFLX Package is active. However, the overall solute distributions did not differ greatly from simulations made without using the BFLX Package.

Open-File Report↗
USGS2011

Rapid estimation of the economic consequences of global earthquakes

The U.S. Geological Survey's (USGS) Prompt Assessment of Global Earthquakes for Response (PAGER) system, operational since mid 2007, rapidly estimates the most affected locations and the population exposure at different levels of shaking intensities. The PAGER system has significantly improved the way aid agencies determine the scale of response needed in the aftermath of an earthquake. For example, the PAGER exposure estimates provided reasonably accurate assessments of the scale and spatial extent of the damage and losses following the 2008 Wenchuan earthquake (Mw 7.9) in China, the 2009 L'Aquila earthquake (Mw 6.3) in Italy, the 2010 Haiti earthquake (Mw 7.0), and the 2010 Chile earthquake (Mw 8.8). Nevertheless, some engineering and seismological expertise is often required to digest PAGER's exposure estimate and turn it into estimated fatalities and economic losses. This has been the focus of PAGER's most recent development. With the new loss-estimation component of the PAGER system it is now possible to produce rapid estimation of expected fatalities for global earthquakes (Jaiswal and others, 2009). While an estimate of earthquake fatalities is a fundamental indicator of potential human consequences in developing countries (for example, Iran, Pakistan, Haiti, Peru, and many others), economic consequences often drive the responses in much of the developed world (for example, New Zealand, the United States, and Chile), where the improved structural behavior of seismically resistant buildings significantly reduces earthquake casualties. Rapid availability of estimates of both fatalities and economic losses can be a valuable resource. The total time needed to determine the actual scope of an earthquake disaster and to respond effectively varies from country to country. It can take days or sometimes weeks before the damage and consequences of a disaster can be understood both socially and economically. The objective of the U.S. Geological Survey's PAGER system is to reduce this time gap to more rapidly and effectively mobilize response. We present here a procedure to rapidly and approximately ascertain the economic impact immediately following a large earthquake anywhere in the world. In principle, the approach presented is similar to the empirical fatality estimation methodology proposed and implemented by Jaiswal and others (2009). In order to estimate economic losses, we need an assessment of the economic exposure at various levels of shaking intensity. The economic value of all the physical assets exposed at different locations in a given area is generally not known and extremely difficult to compile at a global scale. In the absence of such a dataset, we first estimate the total Gross Domestic Product (GDP) exposed at each shaking intensity by multiplying the per-capita GDP of the country by the total population exposed at that shaking intensity level. We then scale the total GDP estimated at each intensity by an exposure correction factor, which is a multiplying factor to account for the disparity between wealth and/or economic assets to the annual GDP. The economic exposure obtained using this procedure is thus a proxy estimate for the economic value of the actual inventory that is exposed to the earthquake. The economic loss ratio, defined in terms of a country-specific lognormal cumulative distribution function of shaking intensity, is derived and calibrated against the losses from past earthquakes. This report describes the development of a country or region-specific economic loss ratio model using economic loss data available for global earthquakes from 1980 to 2007. The proposed model is a potential candidate for directly estimating economic losses within the currently-operating PAGER system. PAGER's other loss models use indirect methods that require substantially more data (such as building/asset inventories, vulnerabilities, and the asset values exposed at the time of earthquake) to implement on a global basis and will thus take more time to develop and implement within the PAGER system.

Open-File Report↗
USGS2019

Flood-frequency estimates for Ohio streamgages based on data through water year 2015 and techniques for estimating flood-frequency characteristics of rural, unregulated Ohio streams

Estimates of the magnitudes of annual peak streamflows with annual exceedance probabilities of 0.5, 0.2, 0.1, 0.04, 0.02, 0.01, and 0.002 (equivalent to recurrence intervals of 2-, 5-, 10-, 25-, 50-, 100-, and 500-years, respectively) were computed for 391 streamgages in Ohio and adjacent states based on data collected through the 2015 water year. The flood-frequency estimates were computed following guidance outlined in Bulletin 17C, developed by the Advisory Committee on Water Information. The Bulletin 17C guidelines retain the basic statistical framework of the superseded Bulletin 17B guidelines; however, the Bulletin 17C guidelines add several enhancements including an improved method of moments approach for fitting the log-Pearson Type III (LPIII) distribution to the flood peaks (called the expected moments algorithm), a generalization of the Grubbs Beck low-outlier test (called the Multiple Grubbs Beck test) that permits identification of multiple potentially influential low floods, and new methods for estimating regional skew and uncertainty. Equations for estimating flood-frequency characteristics at ungaged sites on rural, unregulated streams in Ohio were developed with a two-step process involving ordinary least-squares and generalized least-squares regression techniques. Data from 333 streamgages with 10 or more years of unregulated record were screened for redundancy and a regression dataset was selected that was composed of flood-frequency and basin-characteristic data for 275 streamgages in Ohio and adjacent states. Two sets of equations were developed—one set, referred to as the “simple model,” uses regression region and drainage area as regressor variables, and a second set, referred to as the “full model,” uses regression region, drainage area, main-channel slope, and the percentage of the watershed covered by water and wetlands as regressor variables. The average standard errors of prediction ranged from about 40.5 to 46.5 percent for the simple-model equations and from about 37.2 to 40.3 percent for the full-model equations. For sites meeting the rural, unregulated criteria, flood-frequency estimates determined by means of LPIII analyses are reported along with weighted flood-frequency estimates, computed as a function of the LPIII estimates and the regression estimates. For sites with homogenous periods of regulation, flood-frequency estimates determined by means of LPIII analyses are reported. Ninety-five percent confidence limits are reported for all estimates. Values of regressor variables were determined from digital spatial datasets by means of a geographic information system (GIS). The GIS datasets and the new full-model equations have been incorporated into Ohio’s StreamStats application, a web-based, GIS-backed system designed to facilitate the estimation of streamflow statistics at ungaged locations on streams. Seasonal patterns in peak flows were assessed for 295 streamgages in Ohio. Annual peak flows occurred most frequently between January and April, with March having the highest frequency of occurrence. The month with the fewest number of annual peaks was October. Peak-of-record flows occurred most frequently in March, followed by January (months in which two of Ohio’s most severe widespread floods in recent history occurred). None of the peak-of-record flows occurred in October and only two occurred in November. Temporal trend in annual peak flows were assessed for 133 streamgages on unregulated streams in Ohio with 30 or more years of systematic record. Trends were assessed by computing the rank correlation (as measured with the two-sided Kendall’s tau statistic) between time and annual peak flows. Weak but statistically significant trends were indicated at 15 of the 133 streamgages. Of the 15 streamgages with significant trend in annual peak flows, 12 had an upward trend (positive tau) and 3 had a downward trend (negative tau). All 12 streamgages with positive tau values were at latitudes north of 40°33', and streamgages with negative tau values were at latitudes south of 40°33'.

Ohio↗
USGS1995

User's guide to PHREEQC, a computer program for speciation, reaction-path, advective-transport, and inverse geochemical calculations

PHREEQC is a computer program written in the C programming language that is designed to perform a wide variety of aqueous geochemical calculations. PHREEQC is based on an ion-association aqueous model and has capabilities for (1) speciation and saturation-index calculations, (2) reaction-path and advective-transport calculations involving specified irreversible reactions, mixing of solutions, mineral and gas equilibria surface-complex-ation reactions, and ion-exchange reactions, and (3) inverse modeling, which finds sets of mineral and gas mole transfers that account for composition differences between waters, within specified compositional uncertainties. PHREEQC is derived from the Fortran program PHREEQE, but it has been completely rewritten in C with the addition many new capabilities. New features include the capabilities to use redox couples to distribute redox elements among their valence states in speciation calculations; to model ion-exchange and surface-complex-ation reactions; to model reactions with a fixed-pressure, multicomponent gas phase (that is, a gas bubble); to calculate the mass of water in the aqueous phase during reaction and transport calculations; to keep track of the moles of minerals present in the solid phases and determine automatically the thermodynamically stable phase assemblage; to simulate advective transport in combination with PHREEQC's reaction-modeling capability; and to make inverse modeling calculations that allow for uncertainties in the analytical data. The user interface is improved through the use of a simplified approach to redox reactions, which includes explicit mole-balance equations for hydrogen and oxygen; the use of a revised input that is modular and completely free format; and the use of mineral names and standard chemical symbolism rather than index numbers. The use of (2) eliminates nearly all limitations on army sizes, including numbers of elements, aqueous species, solutions, phases, and lengths of character strings. A new equation solver that optimizes a set of equalities subject to both equality and inequality constraints is used to determine the thermodynamically stable set of phases in equilibrium with a solution. A more complete Newton-Raphson formulation, master-species switching, and scaling of the algebraic equations reduce the number of failures of the numerical method in PHREEQC relative to PHREEQE. This report presents the equations that are the basis for chemical equilibrium and inverse-modeling calculations in PHREEQC, describes the input for the program, and presents twelve examples that demonstrate most of the program's capabilities.

Water-Resources Investigations Report↗
USGS2026

Working group on American Eel (WGAMEEL; outputs from 2024 meeting)

The Working Group on American Eel (WGAMEEL) met virtually three times in 2022-2024 to address the five Terms of Reference (ToRs) of its three-year term. The first two ToRs tasked WGAMEEL with listing and evaluating data on American eel landings, abundance indices, and spatial and habitat data and also to describe assessment methods used in Canada and the US. Canada-wide American eel trajectory was estimated using 12 fishery-independent datasets. Generally, longer datasets had more negative trends than shorter ones. Limiting data to the post-2000 years produced fewer negative trends that did not differ from zero, suggesting the observed declines occurred pre-2000. Spatial modelling for American eel requires knowledge and mapping that covers the breadth of habitat types occupied by the species, including freshwater, estuarine, and marine environments. In recent years there has been an expansion of online databases with data from the aquatic environment, particularly in freshwater, with estuarine and marine data less consistently documented. This report broadly compiles abiotic data series of relevance to American eel. A larger challenge for spatial modelling will be acquiring enough high quality, georeferenced biological data sets with suitable observations to assess occurrence, abundance, and trends over time in a spatial framework. The third ToR considered Indigenous Knowledge Systems for American eel. A survey reaching First Nations representatives from four Canadian provinces confirmed the cultural importance of eels in Indigenous communities, and that Indigenous knowledge possessed by the participant groups are place-based and contextual, especially regarding threats impacting eels. The final two ToRs focused on identifying stock assessment modelling approaches applicable to American and European eel, and assessing whether any of these approaches might be appropriate for American eel management moving forward. WGAMEEL evaluated the various approaches for assessing American eel or providing management advice. Two approaches that could be completed in the next few years because of their minimal data needs are index-based methods and catch-only method. A suite of approaches considered by WGAMEEL that would take more time and data were spatial or habitat models, management strategy evaluation, and spawner-per-recruit (SPR) models potentially paired with meta-population models.

ICES Scientific Reports↗
USGS2016

Geology and undiscovered resource assessment of the potash-bearing Central Asia Salt Basin, Turkmenistan, Uzbekistan, Tajikistan, and Afghanistan: Chapter AA in Global mineral resource assessment

Undiscovered potash resources in the Central Asia Salt Basin (CASB) of Turkmenistan, Uzbekistan, Tajikistan, and Afghanistan were assessed as part of a global mineral resource assessment led by the U.S. Geological Survey. The term “potash” refers to potassium-bearing, water-soluble salts derived from evaporite basins, where seawater dried up and precipitated various salt compounds; the word for the element “potassium” is derived from potash. Potash is produced worldwide at amounts exceeding 30 million metric tons per year, mostly for use in fertilizers. The term “potash” is used by industry to refer to potassium chloride, as well as potassium in sulfate, nitrate, and oxide forms. For the purposes of this assessment, the term “potash” refers to potassium ores and minerals and potash ore grades. Resource and production values are usually expressed by industry in terms of K 2 O (potassium oxide) or muriate of potash (KCl, potassium chloride). The CASB hosts significant discovered potash resources and originated in an inland sea during Late Jurassic time. Seawater flowed into the CASB, mostly from its extreme northwestern margin near the modern Caspian Sea, during several evaporation episodes that deposited at least five different packages of evaporites, with virtually all potash in the second and fourth packages. In this study, the CASB was subdivided into three tracts (permissive areas) for evaluation: the Amu Darya tract in the west, the Gissar tract in the center, and the Afghan-Tajik tract in the east. The Gissar and Amu Darya tracts were quantitatively assessed, whereas the Afghan‑Tajik tract was only qualitatively assessed because of the commonly extreme depth (as deep as 7 km) of the Jurassic salt, extensive deformation, and a lack of known potash deposits. Two approaches were used to estimate amounts of undiscovered potash in the CASB. Stratabound evaporite deposits in the Amu Darya tract were evaluated using an Adaptive Geometric Estimation (AGE) approach, which estimates in-place potash volumes and tonnages. The Gissar tract was evaluated by using the AGE approach for stratabound deposits and the three-part form of assessment of Singer and Menzie (2005) for discrete halokinetic deposits. In the three-part form of assessment, numbers of undiscovered deposits were estimated and combined with grade and tonnage models to probabilistically forecast the amount of undiscovered potash. The Amu Darya tract is estimated to contain 38 billion metric tons of undiscovered potash as K 2 O by using the AGE approach for stratabound deposits. The hybrid stratabound-halokinetic Gissar tract is estimated to contain between 1 and 16 billion metric tons of undiscovered potash as K 2 O. Chapter 1 of this report provides an overview of the history of the CASB and summarizes evaporite potash deposition, halokinesis, and dissolution processes that have affected the current distribution of potash-bearing salt in the CASB. Chapter 2 describes the Gissar tract, an uplifted region that contains a mix of stratabound and halokinetic potash deposits and all of the discovered and exploited potash deposits of the CASB. Chapter 3 describes the Amu Darya tract, where evaporite deposits remain flat-lying and undeformed since their original deposition. Chapter 4 describes the highly deformed and compressed Afghan-Tajik tract and what is known of the deeply-buried Jurassic salt. Chapter 5 describes the spatial databases included with this report, which contain a collection of CASB potash information. Appendixes A and B summarize descriptive models for stratabound and halokinetic potash-bearing salt deposits, respectively. Appendix C summarizes the AGE method used to evaluate the Gissar and Amu Darya tracts. Appendixes D and E contain grade and thickness data for the Gissar and Amu Darya tracts. Appendix F provides the SYSTAT script used to estimate undiscovered K2 O in a CASB tract. Appendix G provides a potash glossary, and appendix H provides biographies of assessment participants.

Scientific Investigations Report↗
USGS2026

Appendix 6: Overview of iPOM 2.0: Updates to wolf abundance estimation in Montana

Montana Fish, Wildlife and Parks (MFWP) uses the integrated Patch Occupancy Model (iPOM) to estimate statewide and regional wolf abundance. During the early years of recovery, wolf abundance was tracked using minimum counts. As the population increased, limitations of this method became apparent, motivating development of a model-based approach for estimating abundance. iPOM development began in 2006, was operationalized in 2013, and has been used in its current form (iPOM 1.0) since 2021. Since its inception, iPOM has been developed and presented as an adaptive modeling framework intended to evolve as new data and analytical methods become available1. Here, we briefly present ongoing work to produce the next iteration of iPOM, hereafter iPOM 2.0, and summarize major model updates and resulting preliminary wolf abundance estimates. iPOM uses models for occupancy, territory size, and pack size, along with monitoring data collected by MFWP wolf specialists and potential wolf sightings reported by deer and elk hunters during the 5-week general rifle season each Nov-Dec. These monitoring data, wolf sightings, and covariates for habitat and survey effort are used to model occupancy of packs on a statewide grid of 600-km 2 cells. The area used by wolves is then estimated as each cell’s occupancy estimate multiplied by grid cell area. iPOM estimates pack abundance by dividing this area by predicted territory sizes, and in-pack wolf abundance by multiplying pack abundance by estimated group size. This is combined with a lone wolf rate to estimate total wolf abundance (see Sells et al. 1 for a full description of iPOM 1.0). This technical summary provides a high-level overview of the proposed iPOM updates and their implications for wolf abundance estimates. The estimates presented here represent results as of June 2026 and are provided to inform upcoming discussions surrounding harvest management and regulations (e.g., season dates, quotas, methods). Results are preliminary and remain subject to final model refinement and peer review.

Montana↗
USGS2015

Flood-inundation maps for Grand River, Red Cedar River, and Sycamore Creek near Lansing, Michigan

Digital flood-inundation maps for a total of 19.7 miles of the Grand River, the Red Cedar River, and Sycamore Creek were created by the U.S. Geological Survey (USGS) in cooperation with the City of Lansing, Michigan, and the U.S. Army Corps of Engineers. The flood-inundation maps, which can be accessed through the USGS Flood Inundation Mapping Science Web site at http://water.usgs.gov/osw/flood_inundation/ , show estimates of the areal extent and depth of flooding corresponding to selected water levels (stages) at three USGS streamgages: Grand River at Lansing, MI (04113000), Red Cedar River at East Lansing, MI (04112500), and Sycamore Creek at Holt Road near Holt, MI (04112850). Near-real-time stages at these streamgages can be obtained on the Internet from the USGS National Water Information System at http://waterdata.usgs.gov/ or the National Weather Service (NWS) Advanced Hydrologic Prediction Service at http:/water.weather.gov/ahps/ , which also forecasts flood hydrographs at all of these sites. Each set of flood profiles was computed by means of a one-dimensional step-backwater model. Each model was calibrated to the current stage-discharge relation at each streamgage and to water levels determined with stage sensors (pressure transducers) temporarily deployed along each stream reach. The hydraulic model was used to compute a set of water-surface profiles for flood stages from nearly Action Stage to above Major Flood stage, as reported by the National Weather Service. The computed water-surface profiles were then used in combination with a Geographic Information System digital elevation model derived from light detection and ranging (lidar) data to delineate the approximate areas flooded at each water level. These maps, used in conjunction with real-time USGS streamgage data and NWS forecasting, provide critical information to emergency management personnel and the public. This information is used to plan flood response actions, such as evacuations and road closures, as well as aid in postflood recovery efforts.

Michigan↗
USGS2000

Geophysical, archaeological and historical evidence support a solar-output model for climate change

Although the processes of climate change are not completely understood, an important causal candidate is variation in total solar output. Reported cycles in various climate-proxy data show a tendency to emulate a fundamental harmonic sequence of a basic solar-cycle length (11 years) multiplied by 2(N) (where N equals a positive or negative integer). A simple additive model for total solar-output variations was developed by superimposing a progression of fundamental harmonic cycles with slightly increasing amplitudes. The timeline of the model was calibrated to the Pleistocene/Holocene boundary at 9,000 years before present. The calibrated model was compared with geophysical, archaeological, and historical evidence of warm or cold climates during the Holocene. The evidence of periods of several centuries of cooler climates worldwide called 'little ice ages,' similar to the period anno Domini (A.D.) 1280-1860 and reoccurring approximately every 1,300 years, corresponds well with fluctuations in modeled solar output. A more detailed examination of the climate sensitive history of the last 1,000 years further supports the model. Extrapolation of the model into the future suggests a gradual cooling during the next few centuries with intermittent minor warmups and a return to near little-ice-age conditions within the next 500 years. This cool period then may be followed approximately 1,500 years from now by a return to altithermal conditions similar to the previous Holocene Maximum.

Conference Paper↗
USGS2018

Preliminary synthesis and assessment of environmental flows in the middle Verde River watershed, Arizona

A 3-year study was undertaken to evaluate the suitability of the available modeling tools for characterizing environmental flows in the middle Verde River watershed of central Arizona, describe riparian vegetation throughout the watershed, and estimate sediment mobilization in the river. Existing data on fish and macroinvertebrates were analyzed in relation to basin characteristics, flow regimes, and microhabitat, and a pilot study was conducted that sampled fish and macroinvertebrates and the microhabitats in which they were found. The sampling for the pilot study took place at five different locations in the middle Verde River watershed. This report presents the results of this 3-year study. The Northern Arizona Groundwater Flow Model (NARGFM) was found to be capable of predicting long-term changes caused by alteration of regional recharge (such as may result from climate variability) and groundwater pumping in gaining, losing, and dry reaches of the major streams in the middle Verde River watershed. Over the period 1910 to 2006, the model simulated an increase in dry reaches, a small increase in reaches losing discharge to the groundwater aquifer, and a concurrent decrease in reaches gaining discharge from groundwater. Although evaluations of the suitability of using the NARGFM and Basin Characteristic Model to characterize various streamflow intervals showed that smallerscale basin monthly runoff could be estimated adequately at locations of interest, monthly stream-flow estimates were found unsatisfactory for determining environmental flows. Orthoimagery and Moderate Resolution Imaging Spectroradiometer data were used to quantify stream and riparian vegetation properties related to biotic habitat. The relative abundance of riparian vegetation varied along the main channel of the Verde River. As would be expected, more upland plant species and fewer lowland species were found in the upper-middle section compared to the lower-middle section, and vice-versa. Vegetation changes within the upper-middle and lower-middle reaches are related to differences in climate and hydrology. In general, the riparian vegetation of the middle Verde River watershed is that of a healthy ecosystem’s mixed age, mixed patch structure, likely a result of the mostly unaltered disturbance regime. The frequency of in-river hydrogeomorphic features (pool, riffle, run) varied along the middle Verde River channel. There was a greater abundance of riffle habitat in the upper-middle reach; the lower-middle reach included more pool habitat. The Oak Creek tributary was more homogenous in geomorphic stream habitat composition than West Clear Creek, where runs dominated the upper reaches and pools dominated many of the lower reaches. On the basis of the period of record and discharges recorded at 15-minute intervals, five flows were found to reach the gravel-transport threshold. Sediment mobilization computed with flows averaged over daily time steps yielded just three flows that reached the gravel-transport threshold, and monthly averaged flows yielded none. In the middle Verde River watershed, 15-minute data should be used when possible to evaluate sediment transport in the river system. Data from more than 300 fish surveys conducted from 1992 to 2011 were analyzed using two schemes, one that divided the river into five reaches based on basin characteristics, and a second that divided the river into five reaches based on degree of flow alteration (specifically, diversions). Fish community metrics and assemblage data were used to analyze patterns of species composition and abundance in the two approaches. Overall, native and non-native species were regularly interacting and probably competing for similar resources. Fish abundances were also analyzed in response to floods and other flow metrics. Although the data are limited, native fish abundances increased more rapidly than non-native fish abundances in response to large floods. The basin-characteristic reach analysis showed native fish in greater abundance in the upper-middle reaches of the Verde River watershed and generally decreasing with downstream distance. The median relative abundance of native fish decreased by 50 percent from reach 1 to reach 5. Using the reach scheme based on degree of flow alteration, nondiverted reaches were found to have a greater abundance of native fish than diverted reaches. In heavily diverted reaches, non-native species outnumbered native species. Fish metrics and stream-flow metrics for the 30, 90, and 365-day periods before collection were computed and the results analyzed statistically. Only abundance of all fish species was associated with the 30-day flow metrics. The 90-day flow metrics were generally positively associated with fish metrics, whereas the 365-day flow metrics had more negative correlations. In particular, significant relations were found between fish metrics and the magnitude and frequency of high flows, including maximum monthly flow, median annual number of high-flow events, and median annual maximum streamflow. Native sucker (Catostomidae) populations tended to decrease in periods of extended base flow, and fish in the non-native sunfish family (Centrarchidae) decreased in periods of flashy, high magnitude flows. A pilot study surveyed fish at five locations in the upper part of the middle Verde River watershed as a means to measure microhabitat availability and quantify native and non-native fish use of that available microhabitat. Results indicated that native and non-native species exhibit some clear differences in microhabitat use. Although at least some native and non-native fish were found in each velocity, depth, and substrate category, preferential microhabitat use was common. On a percentage basis, non-native species had a strong preference for slow-moving and deeper water with silt and sand substrate, with a secondary preference for faster moving and very shallow water and a coarse gravel substrate. Native species showed a general preference for somewhat faster, moderate depth water over coarse gravel and had no clear secondary preference. Macroinvertebrate-variables index period, high-flow year, and collection location (upper-middle Verde River, lowermiddle Verde River, or Verde River tributaries) were found to be important explanatory variables in differentiating among community metrics. Overall richness (number of unique taxa), Shannon’s diversity index, and the percent of the most dominant taxa were all highly correlated, but their response to each macroinvertebrate variable was different. The percentage of mayfly (order Ephemeroptera) taxa was significantly higher in Oak Creek and the upper-middle and lower-middle Verde River reaches, locations which have higher flows and more urbanization than other reaches. When community metrics were related to hydrologic metrics, caddisfly (order Trichoptera) populations appeared to increase and mayfly populations to decrease in response to less flashy and more stable streamflows. Conversely, caddisfly populations appeared to decrease and mayfly populations to increase in response to greater flow variability. Six locations along the Verde River were sampled for macroinvertebrates as part of a pilot study associated with this report—(1) below Granite Creek, (2) near Campbell Ranch, (3) at the U.S. Geological Survey Paulden gage, (4) at the Perkinsville Bridge, (5) at the USGS Clarkdale gage, and (6) near the Reitz Ranch property. A nonmetric multidimensional scaling ordination of macroinvertebrate assemblages showed that the Verde River below Granite Creek site was different from the five other sites and that the Perkinsville Bridge and near Reitz Ranch samples had similar community structure. The near Campbell Ranch and Paulden gage locations had similar microhabitat characteristics, with the exception of riparian cover, yet the assemblage structure was very different. The different community composition at Verde River below Granite Creek was likely due to it having the smallest substrate sizes, lowest velocities, shallowest depths, and most riparian cover of the six sites.

Arizona↗
USGS2009

Summary of survival data from juvenile coho salmon in the Klamath River, northern California, 2009

A study of the effects of the discharge from Iron Gate Dam on the Klamath River on juvenile coho salmon during their seaward migration began in 2005. Estimates of fish survival through various reaches of the river downstream of the dam were completed in 2006, 2007, 2008, and 2009. This report describes the estimates of survival during 2009, and is a complement to similar reports for 2006, 2007, and 2008. For each year, a series of numerical models were evaluated to determine apparent survival and recapture probabilities of radio-tagged fish in several river reaches between Iron Gate Hatchery at river kilometer 309 and a site at river kilometer 33. The evaluations indicate that the primary differences among years are in the survivals through reaches upstream of the confluence of the Scott River with the Klamath River. Data from 2009, one of two years when fish from both hatchery and wild origins were available for analysis, indicate that survival of wild and hatchery fish are similar.

California↗
USGS2013

Generalized additive regression models of discharge and mean velocity associated with direct-runoff conditions in Texas: Utility of the U.S. Geological Survey discharge measurement database

A database containing more than 17,700 discharge values and ancillary hydraulic properties was assembled from summaries of discharge measurement records for 424 U.S. Geological Survey streamflow-gauging stations (stream gauges) in Texas. Each discharge exceeds the 90th-percentile daily mean streamflow as determined by period-of-record, stream-gauge-specific, flow-duration curves. Each discharge therefore is assumed to represent discharge measurement made during direct-runoff conditions. The hydraulic properties of each discharge measurement included concomitant cross-sectional flow area, water-surface top width, and reported mean velocity. Systematic and statewide investigation of these data in pursuit of regional models for the estimation of discharge and mean velocity has not been previously attempted. Generalized additive regression modeling is used to develop readily implemented procedures by end-users for estimation of discharge and mean velocity from select predictor variables at ungauged stream locations. The discharge model uses predictor variables of cross-sectional flow area, top width, stream location, mean annual precipitation, and a generalized terrain and climate index (OmegaEM) derived for a previous flood-frequency regionalization study. The mean velocity model uses predictor variables of discharge, top width, stream location, mean annual precipitation, and OmegaEM. The discharge model has an adjusted R-squared value of about 0.95 and a residual standard error (RSE) of about 0.22 base-10 logarithm (cubic meters per second); the mean velocity model has an adjusted R-squared value of about 0.67 and an RSE of about 0.063 fifth root (meters per second). Example applications and computations using both regression models are provided. - See more at: http://ascelibrary.org/doi/abs/10.1061/%28ASCE%29HE.1943-5584.0000635#sthash.jhGyPxgZ.dpuf

Journal of Hydrologic Engineering↗
USGS2006

GFDL's CM2 global coupled climate models. Part I: Formulation and simulation characteristics

The formulation and simulation characteristics of two new global coupled climate models developed at NOAA's Geophysical Fluid Dynamics Laboratory (GFDL) are described. The models were designed to simulate atmospheric and oceanic climate and variability from the diurnal time scale through multicentury climate change, given our computational constraints. In particular, an important goal was to use the same model for both experimental seasonal to interannual forecasting and the study of multicentury global climate change, and this goal has been achieved. Tw o versions of the coupled model are described, called CM2.0 and CM2.1. The versions differ primarily in the dynamical core used in the atmospheric component, along with the cloud tuning and some details of the land and ocean components. For both coupled models, the resolution of the land and atmospheric components is 2?? latitude ?? 2.5?? longitude; the atmospheric model has 24 vertical levels. The ocean resolution is 1?? in latitude and longitude, with meridional resolution equatorward of 30?? becoming progressively finer, such that the meridional resolution is 1/3?? at the equator. There are 50 vertical levels in the ocean, with 22 evenly spaced levels within the top 220 m. The ocean component has poles over North America and Eurasia to avoid polar filtering. Neither coupled model employs flux adjustments. The co ntrol simulations have stable, realistic climates when integrated over multiple centuries. Both models have simulations of ENSO that are substantially improved relative to previous GFDL coupled models. The CM2.0 model has been further evaluated as an ENSO forecast model and has good skill (CM2.1 has not been evaluated as an ENSO forecast model). Generally reduced temperature and salinity biases exist in CM2.1 relative to CM2.0. These reductions are associated with 1) improved simulations of surface wind stress in CM2.1 and associated changes in oceanic gyre circulations; 2) changes in cloud tuning and the land model, both of which act to increase the net surface shortwave radiation in CM2.1, thereby reducing an overall cold bias present in CM2.0; and 3) a reduction of ocean lateral viscosity in the extratropics in CM2.1, which reduces sea ice biases in the North Atlantic. Both models have be en used to conduct a suite of climate change simulations for the 2007 Intergovernmental Panel on Climate Change (IPCC) assessment report and are able to simulate the main features of the observed warming of the twentieth century. The climate sensitivities of the CM2.0 and CM2.1 models are 2.9 and 3.4 K, respectively. These sensitivities are defined by coupling the atmospheric components of CM2.0 and CM2.1 to a slab ocean model and allowing the model to come into equilibrium with a doubling of atmospheric CO2. The output from a suite of integrations conducted with these models is freely available online (see http://nomads.gfdl.noaa.gov/). ?? 2006 American Meteorological Society.

Journal of Climate↗
USGS2022

An introduction to current climate projections and their use in climate impacts research

Using climate projections to evaluate future climate impacts and their associated risks requires a background knowledge of the nature of climate change, use of climate models to develop future projections, and knowledge of how to address climate scenario uncertainty. This chapter provides an overview of climate and climate change, some of the foundational climate science that underlies current climate change assessments, and a brief introduction to climate models and climate scenario uncertainty. Global projections of temperature and precipitation changes from the recent Intergovernmental Panel on Climate Change (IPCC) Sixth Assessment Report (AR6) and a brief comparison to the prior assessment (AR5) are provided. The main sources of uncertainty in these projections include climate variability, climate model differences and treatment of scientific knowledge gaps, and greenhouse gas (GHG) emissions. When projections are downscaled to local resolution, downscaling is an additional source of uncertainty. These uncertainties can be incorporated in assessments of climate impacts by choosing a range of scenarios that directly address the sources of uncertainty. Evaluating the likelihood of a given climate impact on animal health or management strategies requires consideration of the main sources of climate projection uncertainties. Adaptation requires consideration of global-to-regional contexts of climate changes and impacts, but also adaptive capacity.

Book chapter↗
USGS2014

Modeling marbled murrelet (Brachyramphus marmoratus) habitat using LiDAR-derived canopy data

LiDAR (Light Detection And Ranging) is an emerging remote-sensing tool that can provide fine-scale data describing vertical complexity of vegetation relevant to species that are responsive to forest structure. We used LiDAR data to estimate occupancy probability for the federally threatened marbled murrelet ( Brachyramphus marmoratus ) in the Oregon Coast Range of the United States. Our goal was to address the need identified in the Recovery Plan for a more accurate estimate of the availability of nesting habitat by developing occupancy maps based on refined measures of nest-strand structure. We used murrelet occupancy data collected by the Bureau of Land Management Coos Bay District, and canopy metrics calculated from discrete return airborne LiDAR data, to fit a logistic regression model predicting the probability of occupancy. Our final model for stand-level occupancy included distance to coast, and 5 LiDAR-derived variables describing canopy structure. With an area under the curve value (AUC) of 0.74, this model had acceptable discrimination and fair agreement (Cohen's κ = 0.24), especially considering that all sites in our sample were regarded by managers as potential habitat. The LiDAR model provided better discrimination between occupied and unoccupied sites than did a model using variables derived from Gradient Nearest Neighbor maps that were previously reported as important predictors of murrelet occupancy (AUC = 0.64, κ = 0.12). We also evaluated LiDAR metrics at 11 known murrelet nest sites. Two LiDAR-derived variables accurately discriminated nest sites from random sites (average AUC = 0.91). LiDAR provided a means of quantifying 3-dimensional canopy structure with variables that are ecologically relevant to murrelet nesting habitat, and have not been as accurately quantified by other mensuration methods.

Oregon↗
USGS2006

The genetic basis of hair whorl, handedness, and other phenotypes

Evidence is presented that RHD , RHCE , and other RH genes, may be interesting candidates to consider when searching for the genetic basis of hair whorl rotation (i.e., clockwise or counterclockwise), handedness (i.e., right handed, left handed or ambidextrous), speech laterality (i.e., right brained or left brained), speech dyslexia (e.g., stuttering), sexual orientation (i.e., heterosexual, homosexual, bisexual, or transsexual), schizophrenia, bipolar disorder, and autism spectrum disorder. Such evidence involves the need for a genetic model that includes maternal immunization to explain some of the empirical results reported in the literature. The complex polymorphisms present among the maternally immunizing RH genes can then be used to explain other empirical results. Easily tested hypotheses are suggested, based upon genotypic (but not phenotypic) frequencies of the RH genes. In particular, homozygous dominant individuals are expected to be less common or lacking entirely among the alternative phenotypes. If it is proven that RH genes are involved in brain architecture, it will have a profound effect upon our understanding of the development and organization of the asymmetrical vertebrate brain and may eventually lead to a better understanding of the developmental processes which occur to produce the various alternative phenotypes discussed here. In addition, if RH genes are shown to be involved in the production of these phenotypes, then the evolutionary studies can be performed to demonstrate the beneficial effect of the recessive alleles of RHD and RHCE , and why human evolution appears to be selecting for the recessive alleles even though an increase in the frequency of such alleles may imply lower average fecundity among some individuals possessing them.

Medical Hypotheses↗
USGS1985

Franciscan complex calera limestones: Accreted remnants of farallon plate oceanic plateaus

The Calera Limestone, part of the Franciscan Complex of northern California, may have formed in a palaeoenvironment similar to Hess and Shatsky Rises of the present north-west Pacific1. We report here new palaeomagnetic results, palaeontological data and recent plate-motion models that reinforce this assertion. The Calera Limestone may have formed on Farallon Plate plateaus, north of the Pacific-Farallon spreading centre as a counterpart to Hess or Shatsky Rises. In one model2, the plateaus were formed by hotspots close to the Farallon_Pacific ridge axis. On accretion to North America, plateau dissection in the late Cretaceous to Eocene (50-70 Myr) could explain the occurrence of large volumes of pillow basalt and exotic blocks of limestone in the Franciscan Complex. Partial subduction of the plateaus could have contributed to Laramide (70-40 Myr) compressional events3. ?? 1985 Nature Publishing Group.

Nature↗