Search USGS⌕ Search

SEARCH · Search USGS

Results for “Statistical Methodology”

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

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

At least 199 records · Page 11Linked to original sources
USGS2012

Estimation of baseline daily mean streamflows for ungaged locations on Pennsylvania streams, water years 1960-2008

Water-resource managers use daily mean streamflows to generate streamflow statistics and analyze streamflow conditions. An in-depth evaluation of flow regimes to promote instream ecological health often requires streamflow information obtainable only from a time series hydrograph. Historically, it has been difficult to estimate daily mean streamflow for an ungaged location. The U.S. Geological Survey (USGS), in cooperation with the Pennsylvania Department of Environmental Protection, Susquehanna River Basin Commission, and The Nature Conservancy, has developed the Ba seline S treamflow E stimator (BaSE) to estimate baseline streamflow at a daily time scale for ungaged streams in Pennsylvania using data collected during water years 1960–2008. Baseline streamflow is minimally altered by regulation, diversion, or mining, and other anthropogenic activities. Daily mean streamflow is estimated in BaSE using a methodology that equates streamflow as a percentile from a flow duration curve for a particular day at an ungaged location with streamflow as a percentile from the flow duration curve for the same day at a reference streamgage that is considered to be hydrologically similar to the ungaged location. An appropriate reference streamgage is selected using map correlation, in which variogram models are developed that correlate streamflow at one streamgage with streamflows at all other streamgages. The percentiles from a flow duration curve for the ungaged location are converted to streamflow through the use of regression equations. Regression equations used to predict 17 flow-duration exceedance probabilities were developed for Pennsylvania using geographic information system-derived basin characteristics. The standard error of prediction for the regression equations ranged from 11 percent to 92 percent with the mean of 31 percent. The map correlation method for estimating streamflow was tested at locations within two pilot basins, the Upper Delaware River Basin and the Lower Susquehanna River Basin, before being applied statewide. Reference streamgages within the pilot basins were used as ungaged locations for analyzing the map correlation method. Correlation using Spearman’s rho and centroid distance performed as well as, or better than, the method using the closest streamgage as a reference streamgage. Map correlation using the correlation metrics identified in the pilot basins was applied to 156 streamgages in and near Pennsylvania. BaSE uses the map correlation method and flow-duration exceedance probability regression equations to estimate baseline daily mean streamflow for an ungaged location. The output from BaSE is a Microsoft Excel® report file that summarizes the reference streamgage and ungaged location information, including basin characteristics, percent difference in basin characteristics between the two locations, any warning associated with the basin characteristics, mean and median streamflow for the ungaged location, and a daily hydrograph of streamflow for water years 1960–2008 for the ungaged location. The daily mean streamflow for the ungaged location can be exported as a text file to be used as input into other statistical software packages. BaSE estimates daily mean streamflow for baseline conditions only, and any alterations to streamflow from regulation, large water use, or substantial mining are not reflected in the estimated streamflow.

Pennsylvania↗
USGS2013

A review of selected inorganic surface water quality-monitoring practices: are we really measuring what we think, and if so, are we doing it right?

Successful environmental/water quality-monitoring programs usually require a balance between analytical capabilities, the collection and preservation of representative samples, and available financial/personnel resources. Due to current economic conditions, monitoring programs are under increasing pressure to do more with less. Hence, a review of current sampling and analytical methodologies, and some of the underlying assumptions that form the bases for these programs seems appropriate, to see if they are achieving their intended objectives within acceptable error limits and/or measurement uncertainty, in a cost-effective manner. That evaluation appears to indicate that several common sampling/processing/analytical procedures (e.g., dip (point) samples/measurements, nitrogen determinations, total recoverable analytical procedures) are generating biased or nonrepresentative data, and that some of the underlying assumptions relative to current programs, such as calendar-based sampling and stationarity are no longer defensible. The extensive use of statistical models as well as surrogates (e.g., turbidity) also needs to be re-examined because the hydrologic interrelationships that support their use tend to be dynamic rather than static. As a result, a number of monitoring programs may need redesigning, some sampling and analytical procedures may need to be updated, and model/surrogate interrelationships may require recalibration.

Environmental Science & Technology↗
USGS2020

A modeling workflow that balances automation and human intervention to inform invasive plant management decisions at multiple spatial scales

Predictions of habitat suitability for invasive plant species can guide risk assessments at regional and national scales and inform early detection and rapid-response strategies at local scales. We present a general approach to invasive species modeling and mapping that meets objectives at multiple scales. Our methodology is designed to balance trade-offs between developing highly customized models for few species versus fitting non-specific and generic models for numerous species. We developed a national library of environmental variables known to physiologically limit plant distributions and relied on human input based on natural history knowledge to further narrow the variable set for each species before developing habitat suitability models. To ensure efficiency, we used largely automated modeling approaches and human input only at key junctures. We explore and present uncertainty by using two alternative sources of background samples, including five statistical algorithms, and constructing model ensembles. We demonstrate the use and efficiency of the Software for Assisted Habitat Modeling [SAHM 2.1.2], a package in VisTrails, which performs the majority of the modeling analyses. Our workflow includes solicitation of expert feedback on model outputs such as spatial prediction results and variable response curves, and iterative improvement based on new data availability and directed field validation of initial model results. We highlight the utility of the models for decision-making at regional and local scales with case studies of two plant species that invade natural areas: fountain grass ( Pennisetum setaceum ) and goutweed ( Aegopodium podagraria ). By balancing model automation with human intervention, we can efficiently provide land managers with mapped predicted distributions for multiple invasive species to inform decisions across spatial scales.

PLoS ONE↗
USGS2013

Salmonids, stream temperatures, and solar loading--modeling the shade provided to the Klamath River by vegetation and geomorphology

The U.S. Geological Survey is studying approaches to characterize the thermal regulation of water and the dynamics of cold water refugia. High temperatures have physiological impacts on anadromous fish species. Factors affecting the presence, variability, and quality of thermal refugia are known, such as riverine and watershed processes, hyporheic flows, deep pools and bathymetric factors, thermal stratification of reservoirs, and other broader climatic considerations. This research develops a conceptual model and methodological techniques to quantify the change in solar insolation load to the Klamath River caused by riparian and floodplain vegetation, the morphology of the river, and the orientation and topographic characteristics of its watersheds. Using multiple scales of input data from digital elevation models and airborne light detection and ranging (LiDAR) derivatives, different analysis methods yielded three different model results. These models are correlated with thermal infrared imagery for ground-truth information at the focal confluence with the Scott River. Results from nonparametric correlation tests, geostatistical cross-covariograms, and cross-correlograms indicate that statistical relationships between the insolation models and the thermal infrared imagery exist and are significant. Furthermore, the use of geostatistics provides insights to the spatial structure of the relationships that would not be apparent otherwise. To incorporate a more complete representation of the temperature dynamics in the river system, other variables including the factors mentioned above, and their influence on solar loading, are discussed. With similar datasets, these methods could be applied to any river in the United States—especially those listed as temperature impaired under Section 303(d) of the Clean Water Act—or international riverine systems. Considering the importance of thermal refugia for aquatic species, these methods can help investigate opportunities for riparian restoration, identify problematic reaches unlikely to provide good habitat, and simulate changes to solar loading estimates from alternative landscape configurations.

California↗
USGS2002

Standard errors of annual discharge and change in reservoir content data from selected stations in the lower Colorado River streamflow-gaging station network, 1995-99

The Bureau of Reclamation is currently (1995–2001) testing the Lower Colorado River Accounting System as a method to estimate the consumptive use of Colorado River water by diverters from Hoover Dam to Mexico. Consumptive use is estimated in the Lower Colorado River Accounting System method, in part, on the basis of the annual discharge or annual change in reservoir contents, as well as the variance of estimate of the annual discharge or the annual change in reservoir contents at several surface-water gaging stations in the lower Colorado River stream-gaging network. The standard error and the variance of estimate were determined for the annual discharge at 14 streamflow-gaging stations and for the annual change in content at 2 reservoir-content gaging stations used in the Lower Colorado River Accounting System for calendar years 1995–99. The standard error of the annual discharge was determined by using modifications to an existing method that assumes that the uncertainty in the discharge-rating shift is the main source of uncertainty in computed discharges and that the discharge-rating shift behaves as a first-order Markovian process. The method uses Kalman filtering of a first-order Markovian process as a statistical analogy to computing streamflow with a shifted discharge rating. Temporally unbiased residuals from a discharge rating are used as a surrogate for the actual shifts used to compute discharge. The standard error of the annual discharge is determined by using Kalman-filter theory and estimates of four parameters: (1) the measurement variance of the discharge measurements used to determine the discharge-rating shift, (2) the process variance of the discharge-rating residuals, (3) the serial correlation of the discharge-rating residuals, and (4) the frequency of the discharge measurements. The existing methodology was improved by estimating the measurement variance from a semivariogram of the discharge-rating residuals, rather than on the basis of empirically derived error estimates for discharge measurements. The process variance and serial correlation of the discharge-rating residuals are estimated from the semivariogram, rather than a variogram, of the discharge-rating residuals. The empirically derived estimates are based on characteristics of the discharge measurements such as number of depth and velocity observation sections, type of current meter, and bed material composition and stability. Measurement variance determined from the semivariograms was site specific and is therefore considered a better estimate than measurement variance determined from the empirically-derived estimates. The method of estimating the standard error of the annual discharge requires the assumption of unbiased discharge-rating residuals, and for this reason, the standard errors presented in this report only represent the random error in the annual discharge data. Estimates of the standard error of the annual change in reservoir content were determined on the basis of the reservoir-surface area and the standard error of reservoir-stage readings. The standard error of the annual discharge, as a percentage, ranged from 0.11 percent for the All- American Canal near Imperial Dam in 1998 to 12.3 percent for the Colorado River below Imperial Dam in 1996. The standard error of the annual discharge was less than 2 percent for all 5 years for 11 of the 14 streamflow-gaging stations. In terms of flow volume, the standard error of the annual discharge ranged from 97 acre-feet for the Mittry Lake Diversions in 1995 to 77,000 acre-feet for the Colorado River at the northerly international boundary with Mexico in 1998. In general, the standard error of the annual discharge, as a percentage, was smallest at streamflow-gaging stations on the main stem of the Colorado River; however, the standard error of the annual discharge in acre-feet was largest at these stations because of the large annual discharge on the main stem. The standard error of the annual change in content for the two reservoirs ranged from 1,590 acre-feet for Lake Havasu in 1996 to 2,790 acre-feet for Lake Mohave in 1995. The variance of estimate of the annual discharge for a streamflow-gaging station can be reduced by making additional discharge measurements; either by increasing the number of discharge measurements made per site visit, or by increasing the frequency of site visits. Measurement error can be reduced by using the average shift for two or more discharge measurements made during a site visit. For a streamflow-gaging station where measurement error is much greater than process error and the serial correlation of the discharge-rating residuals is high, an improved gaging strategy would involve making multiple discharge measurements per site visit. In contrast, for a streamflow-gaging station where process error is much greater than measurement error and the serial correlation of discharge-rating residuals is low, the gaging strategy would consist of several single discharge-measurement site visits. For a given operating cost or for a given variance of estimate of the annual discharge at a streamflow-gaging station, the optimal site-visit and discharge-measurement strategy can be determined, providing that the travel costs as well as the measurement variance, process variance, and serial correlation of discharge-rating residuals are known.

Arizona;California↗
USGS2016

Modeling the effects of land cover and use on landscape capability for urban ungulate populations

Expanding ungulate populations are causing concerns for wildlife professionals and residents in many urban areas worldwide. Nowhere is the phenomenon more apparent than in the eastern US, where urban white-tailed deer (Odocoileus virginianus) populations are increasing. Most habitat suitability models for deer have been developed in rural areas and across large (>1000 km2) spatial extents. Only recently have we begun to understand the factors that contribute to space use by deer over much smaller spatial extents. In this study, we explore the concepts, terminology, methodology and state-of-the-science in wildlife abundance modeling as applied to overabundant deer populations across heterogeneous urban landscapes. We used classified, high-resolution digital orthoimagery to extract landscape characteristics in several urban areas of upstate New York. In addition, we assessed deer abundance and distribution in 1-km2 blocks across each study area from either aerial surveys or ground-based distance sampling. We recorded the number of detections in each block and used binomial mixture models to explore important relationships between abundance and key landscape features. Finally, we cross-validated statistical models of abundance and compared covariate relationships across study sites. Study areas were characterized along a gradient of urbanization based on the proportions of impervious surfaces and natural vegetation which, based on the best-supported models, also distinguished blocks potentially occupied by deer. Models performed better at identifying occurrence of deer and worse at predicting abundance in cross-validation comparisons. We attribute poor predictive performance to differences in deer population trajectories over time. The proportion of impervious surfaces often yielded better predictions of abundance and occurrence than did the proportion of natural vegetation, which we attribute to a lack of certain land cover classes during cold and snowy winters. Merits and limitations of our approach to habitat suitability modeling are discussed in detail.

Book chapter↗
USGS2013

Ecological prediction with nonlinear multivariate time-frequency functional data models

Time-frequency analysis has become a fundamental component of many scientific inquiries. Due to improvements in technology, the amount of high-frequency signals that are collected for ecological and other scientific processes is increasing at a dramatic rate. In order to facilitate the use of these data in ecological prediction, we introduce a class of nonlinear multivariate time-frequency functional models that can identify important features of each signal as well as the interaction of signals corresponding to the response variable of interest. Our methodology is of independent interest and utilizes stochastic search variable selection to improve model selection and performs model averaging to enhance prediction. We illustrate the effectiveness of our approach through simulation and by application to predicting spawning success of shovelnose sturgeon in the Lower Missouri River.

Lower Missouri River↗
USGS2020

Novel molecular resources to facilitate future genetics research on freshwater mussels (Bivalvia: Unionidae)

Molecular data have been an integral tool in the resolution of the evolutionary relationships and systematics of freshwater mussels, despite the limited number of nuclear markers available for Sanger sequencing. To facilitate future studies, we evaluated the phylogenetic informativeness of loci from the recently published anchored hybrid enrichment (AHE) probe set Unioverse and developed novel Sanger primer sets to amplify two protein-coding nuclear loci with high net phylogenetic informativeness scores: fem-1 homolog C (FEM1) and UbiA prenyltransferase domain-containing protein 1 (UbiA). We report the methods used for marker development, along with the primer sequences and optimized PCR and thermal cycling conditions. To demonstrate the utility of these markers, we provide haplotype networks, DNA alignments, and summary statistics regarding the sequence variation for the two protein-coding nuclear loci (FEM1 and UbiA). Additionally, we compare the DNA sequence variation of FEM1 and UbiA to three loci commonly used in freshwater mussel genetic studies: the mitochondrial genes cytochrome c oxidase subunit 1 (CO1) and NADH dehydrogenase subunit 1 (ND1), and the nuclear internal transcribed spacer 1 (ITS1). All five loci distinguish among the three focal species (Potamilus fragilis, Potamilus inflatus, and Potamilus purpuratus), and the sequence variation was highest for ND1, followed by CO1, ITS1, UbiA, and FEM1, respectively. The newly developed Sanger PCR primers and methodologies for extracting additional loci from AHE probe sets have great potential to facilitate molecular investigations targeting supraspecific relationships in freshwater mussels, but may be of limited utility at shallow taxonomic scales.

Data↗
USGS2018

ShakeMap-based prediction of earthquake-induced mass movements in Switzerland calibrated on historical observations

In Switzerland, nearly all historical Mw ~ 6 earthquakes have induced damaging landslides, rockslides and snow avalanches that, in some cases, also resulted in damage to infrastructure and loss of lives. We describe the customisation to Swiss conditions of a globally calibrated statistical approach originally developed to rapidly assess earthquake-induced landslide likelihoods worldwide. The probability of occurrence of such earthquake-induced effects is modelled through a set of geospatial susceptibility proxies and peak ground acceleration. The predictive model is tuned to capture the observations from past events and optimised for near-real-time estimates based on USGS-style ShakeMaps routinely produced by the Swiss Seismological Service. Our emphasis is on the use of high-resolution geospatial datasets along with additional local information on ground failure susceptibility. Even if calibrated on historic events with moderate magnitudes, the methodology presented in this paper yields sensible results also for low-magnitude recent events. The model is integrated in the Swiss ShakeMap framework. This study has a high practical relevance to many Swiss ShakeMap stakeholders, especially those managing lifeline systems, and to other global users interested in conducting a similar customisation for their region of interest.

Natural Hazards↗
USGS2016

Grassland and cropland net ecosystem production of the U.S. Great Plains: Regression tree model development and comparative analysis

This paper presents the methodology and results of two ecological-based net ecosystem production (NEP) regression tree models capable of up scaling measurements made at various flux tower sites throughout the U.S. Great Plains. Separate grassland and cropland NEP regression tree models were trained using various remote sensing data and other biogeophysical data, along with 15 flux towers contributing to the grassland model and 15 flux towers for the cropland model. The models yielded weekly mean daily grassland and cropland NEP maps of the U.S. Great Plains at 250 m resolution for 2000–2008. The grassland and cropland NEP maps were spatially summarized and statistically compared. The results of this study indicate that grassland and cropland ecosystems generally performed as weak net carbon (C) sinks, absorbing more C from the atmosphere than they released from 2000 to 2008. Grasslands demonstrated higher carbon sink potential (139 g C·m −2 ·year −1 ) than non-irrigated croplands. A closer look into the weekly time series reveals the C fluctuation through time and space for each land cover type.

Great Plains↗
USGS2009

Development and application of a mark-recapture model incorporating predicted sex and transitory behaviour

We developed an extension of Cormack-Jolly-Seber models to handle a complex mark-recapture problem in which (a) the sex of birds cannot be determined prior to first moult, but can be predicted on the basis of body measurements, and (b) a significant portion of captured birds appear to be transients (i.e. are captured once but leave the area or otherwise become ' untrappable'). We applied this methodology to a data set of 4184 serins (Serinus serinus) trapped in northeastern Spain during 1985-96, in order to investigate age-, sex-, and time-specific variation in survival rates. Using this approach, we were able to successfully incorporate the majority of ringings of serins. Had we eliminated birds not previously captured (as has been advocated to avoid the problem of transience) we would have reduced our sample sizes by >2000 releases. In addition, we were able to include 1610 releases of birds of unknown (but predicted) sex; these data contributed to the precision of our estimates and the power of statistical tests. We discuss problems with data structure, encoding of the algorithms to compute parameter estimates, model selection, identifiability of parameters, and goodness-of-fit, and make recommendations for the design and analysis of future studies facing similar problems.

Bird Study↗
USGS2009

The Ecology of Parasite-Host Interactions at Montezuma Well National Monument, Arizona - Appreciating the Importance of Parasites

Although parasites play important ecological roles through the direct interactions they have with their hosts, historically that fact has been underappreciated. Today, scientists have a growing appreciation of the scope of such impacts. Parasites have been reported to dominate food webs, alter predator-prey relationships, act as ecosystem engineers, and alter community structure. In spite of this growing awareness in the scientific community, parasites are still often neglected in the consideration of the management and conservation of resources and ecosystems. Given that at least half of the organisms on earth are probably parasitic, it should be evident that the ecological functions of parasites warrant greater attention. In this report, we explore different aspects of parasite-host relationships found at a desert spring pond within Montezuma Well National Monument, Arizona. In three separate but related chapters, we explore interactions between a novel amphipod host and two parasites. First, we identify how host behavior responds to this association and how this association affects interactions with both invertebrate non-host predators and a vertebrate host predator. Second, we look at the human dimension, investigating how human recreation can indirectly affect patterns of disease by altering patterns of vertebrate host space use. Finally - because parasites and diseases are of increasing importance in the management of wildlife species, especially those that are imperiled or of management concern - the third chapter argues that research would benefit from increased attention to the statistical analysis of wildlife disease studies. This report also explores issues of statistical parasitology, providing information that may better inform those designing research projects and analyzing data from studies of wildlife disease. In investigating the nature of parasite-host interactions, the role that relationships play in ecological communities, and how human activities alter these associations, scientists usually make inferences by methods of statistical hypotheses testing. This type of hypothesis testing places additional importance on the analysis and interpretation of parasite-host interactions. We address these ideas in this report, focusing on the following questions: (1) How do two parasites with complex life cycles alter the behavior of a novel amphipod host, and how do host and non-host predators respond to infected amphipod prey? (2) Does human recreation affect spatial patterns of infection in an otherwise natural ecosystem? (3) How is hypothesis-testing applied in studies of wildlife disease? (4) What conclusions can we make about the relative usefulness of these methodologies? and (5) How can the analysis and interpretation of wildlife disease studies be improved? Each chapter of this report contains its own literature-cited section, with tables included in appendixes at the end of the full report.

Open-File Report↗
USGS2022

Methods for robust estimates of tree biomass from pollen accumulation rates: Quantifying paleoecological reconstruction uncertainty

Pollen accumulation rates (PAR, grains cm –2 year –1 ) have been shown to be a reliable but methodologically complex bioproxy for quantitative reconstruction of past tree abundance. In a prior study, we found that the PARs of major tree taxa – Pseudotsuga , Pinus , Notholithocarpus , and the pollen group TC (Taxaceae and Cupressaceae families) – were robust and precise estimators of contemporary tree biomass. This paper expands our earlier work. Here, we more fully evaluate the errors associated with biomass reconstructions to identify weaknesses and recommend improvements in PAR-based reconstructions of forest biomass. We account for uncertainty in our biomass proxy in a formal, coherent fashion. The greatest error was introduced by the age models, underscoring the need for improved statistical approaches to age-depth modeling. Documenting the uncertainty in pollen vegetation models should be standard practice in paleoecology. We also share insights gained from the delineation of the relevant source area of pollen, advances in Bayesian 210 Pb modeling, the importance of site selection, and the use of independent data to corroborate biomass estimates. Lastly, we demonstrate our workflow with a new dataset of reconstructed tree biomass between 1850 and 2018 AD from lakes in the Klamath Mountains, California. Our biomass records followed a broad trend of low mean biomass in the ∼1850s followed by large contemporary increases, consistent with expectations of forest densification due to twentieth century fire suppression policies in the American West. More recent reconstructed tree biomass estimates also corresponded with silviculture treatments occurring within the relevant source area of pollen of our lake sites.

California↗
USGS2000

Methodology and application of combined watershed and ground-water models in Kansas

Increased irrigation in Kansas and other regions during the last several decades has caused serious water depletion, making the development of comprehensive strategies and tools to resolve such problems increasingly important. This paper makes the case for an intermediate complexity, quasi-distributed, comprehensive, large-watershed model, which falls between the fully distributed, physically based hydrological modeling system of the type of the SHE model and the lumped, conceptual rainfall-runoff modeling system of the type of the Stanford watershed model. This is achieved by integrating the quasi-distributed watershed model SWAT with the fully-distributed ground-water model MODFLOW. The advantage of this approach is the appreciably smaller input data requirements and the use of readily available data (compared to the fully distributed, physically based models), the statistical handling of watershed heterogeneities by employing the hydrologic-response-unit concept, and the significantly increased flexibility in handling stream-aquifer interactions, distributed well withdrawals, and multiple land uses. The mechanics of integrating the component watershed and ground-water models are outlined, and three real-world management applications of the integrated model from Kansas are briefly presented. Three different aspects of the integrated model are emphasized: (1) management applications of a Decision Support System for the integrated model (Rattlesnake Creek subbasin); (2) alternative conceptual models of spatial heterogeneity related to the presence or absence of an underlying aquifer with shallow or deep water table (Lower Republican River basin); and (3) the general nature of the integrated model linkage by employing a watershed simulator other than SWAT (Wet Walnut Creek basin). These applications demonstrate the practicality and versatility of this relatively simple and conceptually clear approach, making public acceptance of the integrated watershed modeling system much easier. This approach also enhances model calibration and thus the reliability of model results. (C) 2000 Elsevier Science B.V.Increased irrigation in Kansas and other regions during the last several decades has caused serious water depletion, making the development of comprehensive strategies and tools to resolve such problems increasingly important. This paper makes the case for an intermediate complexity, quasi-distributed, comprehensive, large-watershed model, which falls between the fully distributed, physically based hydrological modeling system of the type of the SHE model and the lumped, conceptual rainfall-runoff modeling system of the type of the Stanford watershed model. This is achieved by integrating the quasi-distributed watershed model SWAT with the fully-distributed ground-water model MODFLOW. The advantage of this approach is the appreciably smaller input data requirements and the use of readily available data (compared to the fully distributed, physically based models), the statistical handling of watershed heterogeneities by employing the hydrologic-response-unit concept, and the significantly increased flexibility in handling stream-aquifer interactions, distributed well withdrawals, and multiple land uses. The mechanics of integrating the component watershed and ground-water models are outlined, and three real-world management applications of the integrated model from Kansas are briefly presented. Three different aspects of the integrated model are emphasized: (1) management applications of a Decision Support System for the integrated model (Rattlesnake Creek subbasin); (2) alternative conceptual models of spatial heterogeneity related to the presence or absence of an underlying aquifer with shallow or deep water table (Lower Republican River basin); and (3) the general nature of the integrated model linkage by employing a watershed simulator other than SWAT (Wet Walnut Creek basin). These applications demonstrate the practicality and ve

Journal of Hydrology↗
USGS2013

Evaluating methodological assumptions of a catch-curve survival estimation of unmarked precocial shorebird chickes

Estimating productivity for precocial species can be difficult because young birds leave their nest within hours or days of hatching and detectability thereafter can be very low. Recently, a method for using a modified catch-curve to estimate precocial chick daily survival for age based count data was presented using Piping Plover ( Charadrius melodus ) data from the Missouri River. However, many of the assumptions of the catch-curve approach were not fully evaluated for precocial chicks. We developed a simulation model to mimic Piping Plovers, a fairly representative shorebird, and age-based count-data collection. Using the simulated data, we calculated daily survival estimates and compared them with the known daily survival rates from the simulation model. We conducted these comparisons under different sampling scenarios where the ecological and statistical assumptions had been violated. Overall, the daily survival estimates calculated from the simulated data corresponded well with true survival rates of the simulation. Violating the accurate aging and the independence assumptions did not result in biased daily survival estimates, whereas unequal detection for younger or older birds and violating the birth death equilibrium did result in estimator bias. Assuring that all ages are equally detectable and timing data collection to approximately meet the birth death equilibrium are key to the successful use of this method for precocial shorebirds.

Waterbirds↗
USGS2013

Capture-recapture methodology

Capture-recapture methods were initially developed to estimate human population abundance, but since that time have seen widespread use for fish and wildlife populations to estimate and model various parameters of population, metapopulation, and disease dynamics. Repeated sampling of marked animals provides information for estimating abundance and tracking the fate of individuals in the face of imperfect detection. Mark types have evolved from clipping or tagging to use of noninvasive methods such as photography of natural markings and DNA collection from feces. Survival estimation has been emphasized more recently as have transition probabilities between life history states and/or geographical locations, even where some states are unobservable or uncertain. Sophisticated software has been developed to handle highly parameterized models, including environmental and individual covariates, to conduct model selection, and to employ various estimation approaches such as maximum likelihood and Bayesian approaches. With these user-friendly tools, complex statistical models for studying population dynamics have been made available to ecologists. The future will include a continuing trend toward integrating data types, both for tagged and untagged individuals, to produce more precise and robust population models.

Book chapter↗
USGS2020

seawaveQ—An R package providing a model and utilities for analyzing trends in chemical concentrations in streams with a seasonal wave (seawave) and adjustment for streamflow (Q) and other ancillary variables, version 2.0.0

The seawaveQ R package provides functionality and help to fit a parametric regression model, SEAWAVE-Q, to pesticide concentration data from stream-water samples to assess trends. The model incorporates the strong seasonality and high degree of censoring common in pesticide data, and users can incorporate numerous ancillary variables such as streamflow anomalies. The model is fitted to pesticide data using maximum likelihood methods for censored data and is robust in terms of pesticide, stream location, and degree of censoring of the concentration data. This R package standardizes this methodology for trend analysis, documents the code, and provides help and tutorial information. In previous investigations, the SEAWAVE-Q model assumed a linear trend across the period analyzed. For short trend periods, this assumption of a linear trend is adequate. However, as the period of record analyzed becomes longer, the assumption of linearity is problematic because of changes in pesticide regulation and use, some of which can be abrupt. In this update to the model, a restricted cubic spline option was added for long trend periods. This option allows for more flexibility in the time component of the model. Bootstrap functionality is included to determine statistical significance. Model results with the new restricted cubic spline option are compared to the linear trend option for two pesticide-site combinations.

Open-File Report↗
USGS2010

Simultaneous modeling of habitat suitability, occupancy, and relative abundance: African elephants in Zimbabwe

The recent development of statistical models such as dynamic site occupancy models provides the opportunity to address fairly complex management and conservation problems with relatively simple models. However, surprisingly few empirical studies have simultaneously modeled habitat suitability and occupancy status of organisms over large landscapes for management purposes. Joint modeling of these components is particularly important in the context of management of wild populations, as it provides a more coherent framework to investigate the population dynamics of organisms in space and time for the application of management decision tools. We applied such an approach to the study of water hole use by African elephants in Hwange National Park, Zimbabwe. Here we show how such methodology may be implemented and derive estimates of annual transition probabilities among three dry-season states for water holes: (1) unsuitable state (dry water holes with no elephants); (2) suitable state (water hole with water) with low abundance of elephants; and (3) suitable state with high abundance of elephants. We found that annual rainfall and the number of neighboring water holes influenced the transition probabilities among these three states. Because of an increase in elephant densities in the park during the study period, we also found that transition probabilities from low abundance to high abundance states increased over time. The application of the joint habitat–occupancy models provides a coherent framework to examine how habitat suitability and factors that affect habitat suitability influence the distribution and abundance of organisms. We discuss how these simple models can further be used to apply structured decision-making tools in order to derive decisions that are optimal relative to specified management objectives. The modeling framework presented in this paper should be applicable to a wide range of existing data sets and should help to address important ecological, conservation, and management problems that deal with occupancy, relative abundance, and habitat suitability.

Hwange National Park↗