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 181 records · Page 10Linked to original sources
USGS2015

Estimating the phenology of elk brucellosis transmission with hierarchical models of cause-specific and baseline hazards

Understanding the seasonal timing of disease transmission can lead to more effective control strategies, but the seasonality of transmission is often unknown for pathogens transmitted directly. We inserted vaginal implant transmitters (VITs) in 575 elk ( Cervus elaphus canadensis ) from 2006 to 2014 to assess when reproductive failures (i.e., abortions or still births) occur, which is the primary transmission route of Brucella abortus , the causative agent of brucellosis in the Greater Yellowstone Ecosystem. Using a survival analysis framework, we developed a Bayesian hierarchical model that simultaneously estimated the total baseline hazard of a reproductive event as well as its 2 mutually exclusive parts (abortions or live births). Approximately, 16% (95% CI = 0.10, 0.23) of the pregnant seropositive elk had reproductive failures, whereas 2% (95% CI = 0.01, 0.04) of the seronegative elk had probable abortions. Reproductive failures could have occurred as early as 13 February and as late as 10 July, peaking from March through May. Model results suggest that less than 5% of likely abortions occurred after 6 June each year and abortions were approximately 5 times more likely in March, April, or May compared to February or June. In western Wyoming, supplemental feeding of elk begins in December and ends during the peak of elk abortions and brucellosis transmission (i.e., Mar and Apr). Years with more snow may enhance elk-to-elk transmission on supplemental feeding areas because elk are artificially aggregated for the majority of the transmission season. Elk-to-cattle transmission will depend on the transmission period relative to the end of the supplemental feeding season, elk seroprevalence, population size, and the amount of commingling. Our statistical approach allowed us to estimate the probability density function of different event types over time, which may be applicable to other cause-specific survival analyses. It is often challenging to assess the cause of death, or in this case whether the reproductive event was an abortion or live birth. Accounting for uncertainty in the event type is an important future addition to our methodological approach.

Wyoming↗
USGS2018

Mean composite fire severity metrics computed with Google Earth Engine offer improved accuracy and expanded mapping potential

Landsat-based fire severity datasets are an invaluable resource for monitoring and research purposes. These gridded fire severity datasets are generally produced with pre-and post-fire imagery to estimate the degree of fire-induced ecological change. Here, we introduce methods to produce three Landsat-based fire severity metrics using the Google Earth Engine (GEE) platform: the delta normalized burn ratio (dNBR), the relativized delta normalized burn ratio (RdNBR), and the relativized burn ratio (RBR). Our methods do not rely on time-consuming a priori scene selection and instead use a mean compositing approach in which all valid pixels (e.g. cloud-free) over a pre-specified date range (pre- and post-fire) are stacked and the mean value for each pixel over each stack is used to produce the resulting fire severity datasets. This approach demonstrates that fire severity datasets can be produced with relative ease and speed compared the standard approach in which one pre-fire and post-fire scene are judiciously identified and used to produce fire severity datasets. We also validate the GEE-derived fire severity metrics using field-based fire severity plots for 18 fires in the western US. These validations are compared to Landsat-based fire severity datasets produced using only one pre- and post-fire scene, which has been the standard approach in producing such datasets since their inception. Results indicate that the GEE-derived fire severity datasets show improved validation statistics compared to parallel versions in which only one pre-fire and post-fire scene are used. We provide code and a sample geospatial fire history layer to produce dNBR, RdNBR, and RBR for the 18 fires we evaluated. Although our approach requires that a geospatial fire history layer (i.e. fire perimeters) be produced independently and prior to applying our methods, we suggest our GEE methodology can reasonably be implemented on hundreds to thousands of fires, thereby increasing opportunities for fire severity monitoring and research across the globe.

Remote Sensing↗
USGS2023

HyWaves: Hybrid downscaling of multimodal wave spectra to nearshore areas

Long-term and accurate wave hindcast databases are often required in different coastal engineering projects. The assessment of the nearshore wave climate is often accomplished by using downscaling techniques to translate offshore waves to coastal areas. However, dynamical downscaling approaches may incur huge computational cost. Additionally, the common use of bulk parameterizations are often not accurate for multidimensional waves. To overcome these limitations, we present a hybrid downscaling approach that combines mathematical algorithms (statistical downscaling) and numerical modeling (dynamical downscaling) over the individual spectral partitions. Every wave partition is downscaled and aggregated afterward by using principles of wave linear theory. By assuming linearity in the propagation of the wave celerity, the application of the method is limited from offshore to intermediate water depths. In addition, the method proposed uses a technique to simplify the spectral boundary conditions in complex domains. The methodology has been applied and validated in the island states of Samoa, American Samoa, Majuro, and Kwajalein, showing good skill at reproducing the spectral hourly time series of significant wave height, peak period, and peak direction. Moreover, an accurate representation of the observed energy spectrum was achieved. This study provides insight into the numerical approximation of the combined sea-swell states while improving the quality of fast spectral forecasting and early warning systems.

Ocean Modeling↗
USGS2015

Estimation of unaltered daily mean streamflow at ungaged streams of New York, excluding Long Island, water years 1961-2010

The lakes, rivers, and streams of New York State provide an essential water resource for the State. The information provided by time series hydrologic data is essential to understanding ways to promote healthy instream ecology and to strengthen the scientific basis for sound water management decision making in New York. The U.S. Geological Survey, in cooperation with The Nature Conservancy and the New York State Energy Research and Development Authority, has developed the New York Streamflow Estimation Tool to estimate a daily mean hydrograph for the period from October 1, 1960, to September 30, 2010, at ungaged locations across the State. The New York Streamflow Estimation Tool produces a complete estimated daily mean time series from which daily flow statistics can be estimated. In addition, the New York Streamflow Estimation Tool provides a means for quantitative flow assessments at ungaged locations that can be used to address the objectives of the Clean Water Act—to restore and maintain the chemical, physical, and biological integrity of the Nation’s waters. The New York Streamflow Estimation Tool uses data from the U.S. Geological Survey streamflow network for selected streamgages in New York (excluding Long Island) and surrounding States with shared hydrologic boundaries, and physical and climate basin characteristics to estimate the natural unaltered streamflow at ungaged stream locations. The unaltered streamflow is representative of flows that are minimally altered by regulation, diversion, or mining, and other anthropogenic activities. With the streamflow network data, flow-duration exceedance probability equations were developed to estimate unaltered streamflow exceedance probabilities at an ungaged location using a methodology that equates streamflow as a percentile from a flow-duration curve for a particular day at a hydrologically similar reference streamgage with streamflow as a percentile from the flow-duration curve for the same day at an ungaged location. The reference streamgage is selected using map correlation, a geostatistical method in which variogram models are developed that correlate streamflow at one streamgage with streamflows at all other locations in the study area. Regression equations used to predict 17 flow-duration exceedance probabilities were developed to estimate the flow-duration curves at ungaged locations for New York using geographic information system-derived basin characteristics. A graphical user interface, with an integrated spreadsheet summary report, has been developed to estimate and display the daily mean streamflows and statistics and to evaluate different water management or water withdrawal scenarios with the estimated monthly data. This package of regression equations, U.S. Geological Survey streamgage data, and spreadsheet application produces an interactive tool to estimate an unaltered daily streamflow hydrograph and streamflow statistics at ungaged sites in New York. Among other uses, the New York Streamflow Estimation Tool can assist water managers with permitting water withdrawals, implementing habitat protection, estimating contaminant loads, or determining the potential affect from chemical spills.

New York↗
USGS2020

Lessons learned implementing an operational continuous United States national land change monitoring capability: The Land Change Monitoring, Assessment, and Projection (LCMAP) approach

Growing demands for temporally specific information on land surface change are fueling a new generation of maps and statistics that can contribute to understanding geographic and temporal patterns of change across large regions, provide input into a wide range of environmental modeling studies, clarify the drivers of change, and provide more timely information for land managers. To meet these needs, the U.S. Geological Survey has implemented a capability to monitor land surface change called the Land Change Monitoring, Assessment, and Projection (LCMAP) initiative. This paper describes the methodological foundations and lessons learned during development and testing of the LCMAP approach. Testing and evaluation of a suite of 10 annual land cover and land surface change data sets over six diverse study areas across the United States revealed good agreement with other published maps (overall agreement ranged from 73% to 87%) as well as several challenges that needed to be addressed to meet the goals of robust, repeatable, and geographically consistent monitoring results from the Continuous Change Detection and Classification (CCDC) algorithm. First, the high spatial and temporal variability of observational frequency led to differences in the number of changes identified, so CCDC was modified such that change detection is dependent on observational frequency. Second, the CCDC classification methodology was modified to improve its ability to characterize gradual land surface changes. Third, modifications were made to the classification element of CCDC to improve the representativeness of training data, which necessitated replacing the random forest algorithm with a boosted decision tree. Following these modifications, assessment of prototype Version 1 LCMAP results showed improvements in overall agreement (ranging from 85% to 90%).

Remote Sensing of Environment↗
USGS2011

Assessing first-order emulator inference for physical parameters in nonlinear mechanistic models

We present an approach for estimating physical parameters in nonlinear models that relies on an approximation to the mechanistic model itself for computational efficiency. The proposed methodology is validated and applied in two different modeling scenarios: (a) Simulation and (b) lower trophic level ocean ecosystem model. The approach we develop relies on the ability to predict right singular vectors (resulting from a decomposition of computer model experimental output) based on the computer model input and an experimental set of parameters. Critically, we model the right singular vectors in terms of the model parameters via a nonlinear statistical model. Specifically, we focus our attention on first-order models of these right singular vectors rather than the second-order (covariance) structure.

Journal of Agricultural, Biological, and Environme↗
USGS2026

Guidelines for producing integrated 210Pb and 14C age-models

Accurate reconstructions of past environmental changes are crucial in paleoecological research and require reliable chronologies of sedimentary archives. Establishing robust age-models and obtaining the most appropriate proxies for analysis is a complex scientific endeavor, requiring extensive resources and collaboration among specialists, including radiochronologists. Radiometric dating methods, such as 210 Pb and radiocarbon ( 14 C), are frequently employed to establish chronologies in aquatic sedimentary deposits and peat bogs. In this study, we review key aspects of sampling, analysis, and the principles underlying 210 Pb and 14 C age-models, focusing on methods for developing robust joint chronologies for paleoenvironmental research. Drawing largely from the authors' experiences and group discussions during and after a scientific workshop in 2022, we discuss important considerations for site selection, sampling strategies, and radiometric dating to construct integrated 210 Pb and 14 C age-models. Using expert consensus, this group – called Paleostats – aims to provide a set of best practices for other geochronologists with this methods paper. Among our conclusions, we emphasize the importance of accounting for site-specific factors such as prior information on sedimentation rates to establish appropriate sampling and analytical strategies. The use of appropriate coring devices can minimize disturbance to sediments and ensure the core surface remains intact and preserved until sectioning. Where excess 210 Pb is expected, sectioning at intervals of ≤1 cm provides an adequate sampling resolution for 210 Pb dating. Exceptions are possible, allowing for ~2–3 cm sections in areas with confirmed high sedimentation rates (e.g., > 1 cm yr −1 ). Recovering deeper core sections for 14 C dating with sufficient overlap allows for accounting errors in depth estimates made in the field. Special attention is advised during time intervals where validation proxies, such as the human-made radionuclides 137 Cs or post-bomb 14 C, are expected, and to determine the depth of secular equilibrium between 210 Pb and 226 Ra. Radiocarbon analyses are commonly performed by accelerator mass spectrometry, and age models are constructed mainly using Bayesian statistics with Markov Chain Monte Carlo techniques (e.g., Bacon ). A Bayesian approach ( Plum ) is now available for producing 210 Pb age-models, which infers the 210 Pb ex flux, eliminates the need for selecting an equilibrium depth, and allows dating cores with incomplete 210 Pb ex inventory. Plum offers improved chronologies by integrating raw 210 Pb and 14 C data, and these age-models can be enriched with other dating methodologies, such as identifying tephras and other well-recorded historical events. Harmonized reporting would contribute to making radiometric age-models reproducible, which would benefit from an international effort. Using 210 Pb and 14 C to produce integrated age-models may yield better insights into the interplay between natural and recent anthropogenic forcings on ecosystems. This can enhance our understanding of environmental processes and their impacts on climate change, ultimately supporting science-based assessments and decisions.

Earth-Science Reviews↗
USGS2015

Linking microbial community structure and microbial processes: An empirical and conceptual overview

A major goal of microbial ecology is to identify links between microbial community structure and microbial processes. Although this objective seems straightforward, there are conceptual and methodological challenges to designing studies that explicitly evaluate this link. Here, we analyzed literature documenting structure and process responses to manipulations to determine the frequency of structure-process links and whether experimental approaches and techniques influence link detection. We examined nine journals (published 2009–13) and retained 148 experimental studies measuring microbial community structure and processes. Many qualifying papers (112 of 148) documented structure and process responses, but few (38 of 112 papers) reported statistically testing for a link. Of these tested links, 75% were significant and typically used Spearman or Pearson's correlation analysis (68%). No particular approach for characterizing structure or processes was more likely to produce significant links. Process responses were detected earlier on average than responses in structure or both structure and process. Together, our findings suggest that few publications report statistically testing structure-process links. However, when links are tested for they often occur but share few commonalities in the processes or structures that were linked and the techniques used for measuring them.

FEMS Microbiology Ecology↗
USGS1994

Estimating breeding proportions and testing hypotheses about costs of reproduction with capture-recapture data

The proportion of animals in a population that breeds is an important determinant of population growth rate. Usual estimates of this quantity from field sampling data assume that the probability of appearing in the capture or count statistic is the same for animals that do and do not breed. A similar assumption is required by most existing methods used to test ecologically interesting hypotheses about reproductive costs using field sampling data. However, in many field sampling situations breeding and nonbreeding animals are likely to exhibit different probabilities of being seen or caught. In this paper, we propose the use of multistate capture-recapture models for these estimation and testing problems. This methodology permits a formal test of the hypothesis of equal capture/sighting probabilities for breeding and nonbreeding individuals. Two estimators of breeding proportion (and associated standard errors) are presented, one for the case of equal capture probabilities and one for the case of unequal capture probabilities. The multistate modeling framework also yields formal tests of hypotheses about reproductive costs to future reproduction or survival or both fitness components. The general methodology is illustrated using capture-recapture data on female meadow voles, Microtus pennsylvanicus. Resulting estimates of the proportion of reproductively active females showed strong seasonal variation, as expected, with low breeding proportions in midwinter. We found no evidence of reproductive costs extracted in subsequent survival or reproduction. We believe that this methodological framework has wide application to problems in animal ecology concerning breeding proportions and phenotypic reproductive costs.

Ecology↗
USGS2007

Integrating laboratory creep compaction data with numerical fault models: A Bayesian framework

[1] We developed a robust Bayesian inversion scheme to plan and analyze laboratory creep compaction experiments. We chose a simple creep law that features the main parameters of interest when trying to identify rate-controlling mechanisms from experimental data. By integrating the chosen creep law or an approximation thereof, one can use all the data, either simultaneously or in overlapping subsets, thus making more complete use of the experiment data and propagating statistical variations in the data through to the final rate constants. Despite the nonlinearity of the problem, with this technique one can retrieve accurate estimates of both the stress exponent and the activation energy, even when the porosity time series data are noisy. Whereas adding observation points and/or experiments reduces the uncertainty on all parameters, enlarging the range of temperature or effective stress significantly reduces the covariance between stress exponent and activation energy. We apply this methodology to hydrothermal creep compaction data on quartz to obtain a quantitative, semiempirical law for fault zone compaction in the interseismic period. Incorporating this law into a simple direct rupture model, we find marginal distributions of the time to failure that are robust with respect to errors in the initial fault zone porosity.

Journal of Geophysical Research B: Solid Earth↗
USGS1985

Population dynamics and interagency management of the bloater (Coregonus hoyi) in Lake Michigan, 1967-1982

This paper examines the population dynamics of the bloater (Coregonus hoyi) in Lake Michigan during a progressive decline in abundance from about the mid1960s through the mid1970s, and during a subsequent recovery that is still underway. The study focused on developing a data base and methodology for projecting fishable surpluses, in cooperation with a chub technical committee sponsored by the Great Lakes Fishery Commission. The Technical Committee was formed in 1974 because of depletion of bloaters and other deepwater ciscoes or 'chubs,' as they are known by Great Lakes fishermen. Subsequently the Technical Committee recommended a lakewide ban on chub fishing that was fully enacted by the states of Illinois, Michigan, and Wisconsin in 1976. With the Committee's help, commercial fishery statistics and stock assessment data were obtained from state and federal research files and used with various indirect analytical techniques to estimate relevant population parameters. The lakewide fishable stock in fall 1973, before the fishery was affected by several incomplete closures and then by the lakewide ban, was estimated as 48 to 73 million bloaters weighing 20 to 29 million pounds. Exploitation of the estimated stock varied considerably among 11 statistical districts in the several jurisdictions. Yield to the fishery exceeded production by the stock in some districts. Theoretical yields of bloaters totaling 3.59 to 3.72 million pounds were projected from 1979 for all waters combined. These projected yields were intended as guidelines for experimental quotas that the states might establish, because the population had stabilized and the potential for recruitment had improved in most areas.

Technical Report↗
USGS2014

Geologic framework for the national assessment of carbon dioxide storage resources: Denver Basin, Colorado, Wyoming, and Nebraska

This is a report about the geologic characteristics of five storage assessment units (SAUs) within the Denver Basin of Colorado, Wyoming, and Nebraska. These SAUs are Cretaceous in age and include (1) the Plainview and Lytle Formations, (2) the Muddy Sandstone, (3) the Greenhorn Limestone, (4) the Niobrara Formation and Codell Sandstone, and (5) the Terry and Hygiene Sandstone Members. The described characteristics, as specified in the methodology, affect the potential carbon dioxide storage resource in the SAUs. The specific geologic and petrophysical properties of interest include depth to the top of the storage formation, average thickness, net-porous thickness, porosity, permeability, groundwater quality, and the area of structural reservoir traps. Descriptions of the SAU boundaries and the overlying sealing units are also included. Assessment results are not contained in this report; however, the geologic information included here will be used to calculate a statistical Monte Carlo-based distribution of potential storage volume in the SAUs.

Colorado, Nebraska, Wyoming↗
USGS2012

Producing fractional rangeland component predictions in a sagebrush ecosystem, a Wyoming sensitivity analysis

Remote sensing information has been widely used to monitor vegetation condition and variations in a variety of ecosystems, including shrublands. Careful application of remotely sensed imagery can provide additional spatially explicit, continuous, and extensive data on the composition and condition of shrubland ecosystems. Historically, the most widely available remote sensing information has been collected by Landsat, which has offered large spatial coverage and moderate spatial resolution data globally for nearly three decades. Such medium-resolution satellite remote sensing information can quantify the distribution and variation of terrestrial ecosystems. Landsat imagery has been frequently used with other high-resolution remote sensing data to classify sagebrush components and quantify their spatial distributions (Ramsey and others, 2004; Seefeldt and Booth, 2004; Stow and others, 2008; Underwood and others, 2007). Modeling algorithms have been developed to use field measurements and satellite remote sensing data to quantify the extent and evaluate the quality of shrub ecosystem components in large geographic areas (Homer and others, 2009). The percent cover of sagebrush ecosystem components, including bare-ground, herbaceous, litter, sagebrush, and shrub, have been quantified for entire western states (Homer and others, 2012). Furthermore, research has demonstrated the use of current measurements with historical archives of Landsat imagery to quantify the variations of these components for the last two decades (Xian and others, 2012). The modeling method used to quantify the extent and spatial distribution of sagebrush components over a large area also has required considerable amounts of training data to meet targeted accuracy requirements. These training data have maintained product accuracy by ensuring that they are derived from good quality field measurements collected during appropriate ecosystem phenology and subsequently maximized by extrapolation on high-resolution remote sensing data (Homer and others, 2012). This method has proven its utility; however, to develop these products across even larger areas will require additional cost efficiencies to ensure that an adequate product can be developed for the lowest cost possible. Given the vast geographic extent of shrubland ecosystems in the western United States, identifying cost efficiencies with optimal training data development and subsequent application to medium resolution satellite imagery provide the most likely areas for methodological efficiency gains. The primary objective of this research was to conduct a series of sensitivity tests to evaluate the most optimal and practical way to develop Landsat scale information for estimating the extent and distribution of sagebrush ecosystem components over large areas in the conterminous United States. An existing dataset of sagebrush components developed from extensive field measurements, high-resolution satellite imagery, and medium resolution Landsat imagery in Wyoming was used as the reference database (Homer and others, 2012). Statistical analysis was performed to analyze the relation between the accuracy of sagebrush components and the amount and distribution of training data on Landsat scenes needed to obtain accurate predictions.

Wyoming↗
USGS2003

Methodology and Results for the Assessment of Oil and Gas resources, National Petroleum Reserve, Alaska

Oil and gas resources in each of the 24 plays within the National Petroleum Reserve in Alaska (NPRA) were estimated using a play analysis. Assessors specified geologic attributes, risks, and number of prospects for each play. Some specifications established distributions, while others were given as single values. From this information, sizes of oil and gas accumulations were simulated using a Monte Carlo algorithm. The number of such accumulations considered in a given simulation run was obtained from the distribution of the number of prospects. Each prospect in each successful simulation run was risked. This process yielded size-frequency distributions and summary statistics for the various petroleum categories. Estimates of remaining resources from individual plays were then aggregated, and measures of uncertainty computed. Technically recoverable, undiscovered oil beneath the Federal part of NPRA likely ranges between 5.9 and 13.2 billion barrels, with a mean (expected) value of 9.3 billion barrels. Technically recoverable, undiscovered nonassociated natural gas for the same area likely ranges between 39.1 and 83.2 trillion cubic feet, with a mean (expected) value of 59.7 trillion cubic feet. Mean values of the corresponding associated dissolved gas and natural gas liquid are 10.3 trillion cubic feet and 1.4 billion barrels respectively.

Open-File Report↗
USGS2022

Quantitative meta-analysis reveals no association between mercury contamination and body condition in birds

Mercury contamination is a major threat to the global environment, and is still increasing in some regions despite international regulations. The methylated form of mercury is hazardous to biota, yet its sublethal effects are difficult to detect in wildlife. Body condition can vary in response to stressors, but previous studies have shown mixed effects of mercury on body condition in wildlife. Using birds as study organisms, we provide the first quantitative synthesis of the effect of mercury on body condition in animals. In addition, we explored the influence of intrinsic, extrinsic and methodological factors potentially explaining cross-study heterogeneity in results. We considered experimental and correlative studies carried out in adult birds and chicks, and mercury exposure inferred from blood and feathers. Most experimental investigations (90%) showed a significant relationship between mercury concentrations and body condition. Experimental exposure to mercury disrupted nutrient (fat) metabolism, metabolic rates, and food intake, resulting in either positive or negative associations with body condition. Correlative studies also showed either positive or negative associations, of which only 14% were statistically significant. Therefore, the overall effect of mercury concentrations on body condition was null in both experimental (estimate ± SE = 0.262 ± 0.309, 20 effect sizes, five species) and correlative studies (−0.011 ± 0.020, 315 effect sizes, 145 species). The single and interactive effects of age class and tissue type were accounted for in meta-analytic models of the correlative data set, since chicks and adults, as well as blood and feathers, are known to behave differently in terms of mercury accumulation and health effects. Of the 15 moderators tested, only wintering status explained cross-study heterogeneity in the correlative data set: free-ranging wintering birds were more likely to show a negative association between mercury and body condition. However, wintering effect sizes were limited to passerines, further studies should thus confirm this trend in other taxa. Collectively, our results suggest that ( i ) effects of mercury on body condition are weak and mostly detectable under controlled conditions, and ( ii ) body condition indices are unreliable indicators of mercury sublethal effects in the wild. Food availability, feeding rates and other sources of variation that are challenging to quantify likely confound the association between mercury and body condition in natura . Future studies could explore the metabolic effects of mercury further using designs that allow for the estimation and/or manipulation of food intake in both wild and captive birds, especially in under-represented life-history stages such as migration and overwintering.

Biological Reviews↗
USGS2014

Users' guide to system dynamics model describing Coho salmon survival in Olema Creek, Point Reyes National Seashore, Marin County, California

The system dynamics model described in this report is the result of a collaboration between U.S. Geological Survey (USGS) scientists and National Park Service (NPS) San Francisco Bay Area Network (SFAN) staff, whose goal was to develop a methodology to integrate inventory and monitoring data to better understand ecosystem dynamics and trends using salmon in Olema Creek, Marin County, California, as an example case. The SFAN began monitoring multiple life stages of coho salmon (Oncorhynchus kisutch) in Olema Creek during 2003 (Carlisle and others, 2013), building on previous monitoring of spawning fish and redds. They initiated water-quality and habitat monitoring, and had access to flow and weather data from other sources. This system dynamics model of the freshwater portion of the coho salmon life cycle in Olema Creek integrated 8 years of existing monitoring data, literature values, and expert opinion to investigate potential factors limiting survival and production, identify data gaps, and improve monitoring and restoration prescriptions. A system dynamics model is particularly effective when (1) data are insufficient in time series length and/or measured parameters for a statistical or mechanistic model, and (2) the model must be easily accessible by users who are not modelers. These characteristics helped us meet the following overarching goals for this model: Summarize and synthesize NPS monitoring data with data and information from other sources to describe factors and processes affecting freshwater survival of coho salmon in Olema Creek. Provide a model that can be easily manipulated to experiment with alternative values of model parameters and novel scenarios of environmental drivers. Although the model describes the ecological dynamics of Olema Creek, these dynamics are structurally similar to numerous other coastal streams along the California coast that also contain anadromous fish populations. The model developed for Olema can be used, at least as a starting point, for other watersheds. This report describes each of the model elements with sufficient detail to guide the primary target audience, the NPS resource specialist, to run the model, interpret the results, change the input data to explore hypotheses, and ultimately modify and improve the model. Running the model and interpreting the results does not require modeling expertise on the part of the user. Additional companion publications will highlight other aspects of the model, such as its development, the rationale behind the methodological approach, scenario testing, and discussions of its use. System dynamics models consist of three basic elements: stocks , flows , and converters . Stocks are measurable quantities that can change over time, such as animal populations. Flows are any processes or conditions that change the quantity in a stock over time (Ford, 1999), are expressed in the model as a rate of change, and are diagrammed as arrows to or from stocks. Converters are processes or conditions that change the rate of flows. A converter is connected to a flow with an arrow indicating that it alters the rate of change. Anything that influences the rate of change (such as different environmental conditions, other external factors, or feedbacks from other stocks or flows) is modeled as a converter. For example, the number of fish in a population is appropriately modeled as a stock. Mortality is modeled as a flow because it is a rate of change over time used to determine the number of fish in the population. The density-dependent effect on mortality is modeled as a converter because it influences the rate of morality. Together, the flow and converter change the number, or stock, of juvenile coho. The instructions embedded in the stocks, flows, converters, and the sequence in which they are linked are processed by the simulation software with each completed sequence composing a model run. At each modeled time step within the model run, the stock counts will go up, down, or stay the same based on the modeled flows and the influence of converters on those flows. The model includes a user-friendly interface to change model parameters, which allows park staff and others to conduct sensitivity analyses, incorporate future knowledge, and implement scenarios for various future conditions. The model structure incorporates place holders for relationships that we hypothesize are significant but data are currently lacking. Future climate scenarios project stream temperatures higher than any that have ever been recorded at Olema Creek. Exploring climate change impacts on coho survival is a high priority for park staff, therefore the model provides the user with the option to experiment with hypothesized effects and to incorporate effects based on future observations.

California↗
USGS2016

Groundwater contaminant plume maps and volumes, 100-K and 100-N Areas, Hanford Site, Washington

This study provides an independent estimate of the areal and volumetric extent of groundwater contaminant plumes which are affected by waste disposal in the 100-K and 100-N Areas (study area) along the Columbia River Corridor of the Hanford Site. The Hanford Natural Resource Trustee Council requested that the U.S. Geological Survey perform this interpolation to assess the accuracy of delineations previously conducted by the U.S. Department of Energy and its contractors, in order to assure that the Natural Resource Damage Assessment could rely on these analyses. This study is based on previously existing chemical (or radionuclide) sampling and analysis data downloaded from publicly available Hanford Site Internet sources, geostatistically selected and interpreted as representative of current (from 2009 through part of 2012) but average conditions for groundwater contamination in the study area. The study is limited in scope to five contaminants—hexavalent chromium, tritium, nitrate, strontium-90, and carbon-14, all detected at concentrations greater than regulatory limits in the past. All recent analytical concentrations (or activities) for each contaminant, adjusted for radioactive decay, non-detections, and co-located wells, were converted to log-normal distributions and these transformed values were averaged for each well location. The log-normally linearized well averages were spatially interpolated on a 50 × 50-meter (m) grid extending across the combined 100-N and 100-K Areas study area but limited to avoid unrepresentative extrapolation, using the minimum curvature geostatistical interpolation method provided by SURFER ® data analysis software. Plume extents were interpreted by interpolating the log-normally transformed data, again using SURFER ® , along lines of equal contaminant concentration at an appropriate established regulatory concentration . Total areas for each plume were calculated as an indicator of relative environmental damage. These plume extents are shown graphically and in tabular form for comparison to previous estimates. Plume data also were interpolated to a finer grid (10 × 10 m) for some processing, particularly to estimate volumes of contaminated groundwater. However, hydrogeologic transport modeling was not considered for the interpolation. The compilation of plume extents for each contaminant also allowed estimates of overlap of the plumes or areas with more than one contaminant above regulatory standards. A mapping of saturated aquifer thickness also was derived across the 100-K and 100–N study area, based on the vertical difference between the groundwater level (water table) at the top and the altitude of the top of the Ringold Upper Mud geologic unit, considered the bottom of the uppermost unconfined aquifer. Saturated thickness was calculated for each cell in the finer (10 × 10 m) grid. The summation of the cells’ saturated thickness values within each polygon of plume regulatory exceedance provided an estimate of the total volume of contaminated aquifer, and the results also were checked using a SURFER ® volumetric integration procedure. The total volume of contaminated groundwater in each plume was derived by multiplying the aquifer saturated thickness volume by a locally representative value of porosity (0.3). Estimates of the uncertainty of the plume delineation also are presented. “Upper limit” plume delineations were calculated for each contaminant using the same procedure as the “average” plume extent except with values at each well that are set at a 95-percent upper confidence limit around the log-normally transformed mean concentrations, based on the standard error for the distribution of the mean value in that well; “lower limit” plumes are calculated at a 5-percent confidence limit around the geometric mean. These upper- and lower-limit estimates are considered unrealistic because the statistics were increased or decreased at each well simultaneously and were not adjusted for correlation among the well distributions (i.e., it is not realistic that all wells would be high simultaneously). Sources of the variability in the distributions used in the upper- and lower-extent maps include time varying concentrations and analytical errors. The plume delineations developed in this study are similar to the previous plume descriptions developed by U.S. Department of Energy and its contractors. The differences are primarily due to data selection and interpolation methodology. The differences in delineated plumes are not sufficient to result in the Hanford Natural Resource Trustee Council adjusting its understandings of contaminant impact or remediation.

Washington↗
USGS1992

Sampling vadose-zone water for a volatile organic compound at Picatinny Arsenal, New Jersey

A new method of collecting samples of unsaturated-zone water for quantitative analysis for a volatile organic compound, trichloroethene (TCE), was compared to three other, previously described sampling methodologies in the laboratory and in the field. In the laboratory, prepared water samples containing TCE in a known concentration (20 µg/L) were sampled repeatedly by using each of the four methods to quantify method precision and accuracy. To compare the four methods in the field, unsaturated-zone water above a TCE-contaminated water-table aquifer was transferred from a depth of 2 m to land surface with 0.15-m-long suction lysimeters attached to 1.85-m lengths of stainless-steel tubing. Statistical analyses of the laboratory and field data indicate that the new method, which involves collecting the water samples in gas-tight glass syringes, is superior to the other three methods for the quantitative sampling and analysis of TCE on the basis of its high precision and accuracy and ease of use. This method was used to collect additional samples from the field site to quantify the spatial variability of TCE concentrations in the unsaturated-zone water. Results of analysis of variance of the data indicate that the spatial concentration variability is important, and that differences in TCE concentration are statistically significant for horizontal distances less than 3.6 m.

New Jersey↗