Geology ReportsSearch

SEARCH · Geology Reports

Results for “Statistical Science”

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 523 records · Page 29Linked to original sources

Trends in plant cover derived from vegetation plot data using ordinal zero-augmented beta regression

Questions Plant cover values in vegetation plot data are bounded between 0 and 1, and cover is typically recorded in discrete classes with non-equal intervals. Consequently, cover data are skewed and heteroskedastic, which hampers the application of conventional regression methods. Recently developed ordinal beta regression models consider these statistical difficulties. Our primary question is whether we can detect species trends in vegetation plot time series data with this modelling approach. A second question is whether trends in cover have additional value compared to trends in occurrence, which are easier to assess for practitioners. Location The Netherlands, Western Europe. Methods We used vegetation plot data collected from 10,000 fixed plots which were surveyed once every four years during 1999–2022. We used the ordinal zero-augmented beta regression (OZAB) model, a hierarchical model consisting of a logistic regression for presence and an ordinal beta regression for cover. We adapted the OZAB model for longitudinal data and produced estimates of cover and occurrence for each four-year period. Thereafter we assessed trends in cover and in occurrence across all periods. Results We found evidence of a trend in cover in 318 out of the 721 species (44%) with sufficient data. Most species showed similar directional trends in occurrence and percent cover. No trend in occurrence was detected for 64 species that had evidence of a trend in cover. Declining species had stronger relative changes in cover than in occurrence. Conclusions Our model enables researchers to detect trends in cover using longitudinal vegetation plot data. Cover trends often corroborated trends in occurrence, but we also regularly found trends in cover even in the absence of evidence for trends in occurrence. Our approach thus contributes to a more complete picture of (changes in) vegetation composition based on large monitoring data sets.

Journal of Vegetation Science

Robust projections of future fire probability for the conterminous United States

Globally increasing wildfires have been attributed to anthropogenic climate change. However, providing decision makers with a clear understanding of how future planetary warming could affect fire regimes is complicated by confounding land use factors that influence wildfire and by uncertainty associated with model simulations of climate change. We use an ensemble of statistically downscaled Global Climate Models in combination with the Physical Chemistry Fire Frequency Model (PC2FM) to project changing potential fire probabilities in the conterminous United States for two scenarios representing lower (RCP 4.5) and higher (RCP 8.5) greenhouse gas emission futures. PC2FM is a physically-based and scale-independent model that predicts mean fire return intervals from both fire reactant and reaction variables, which are largely dependent on a locale's climate. Our results overwhelmingly depict increasing potential fire probabilities across the conterminous US for both climate scenarios. The primary mechanism for the projected increases is rising temperatures, reflecting changes in the chemical reaction environment commensurate with enhanced photosynthetic rates and available thermal molecular energy. Existing high risk areas, such as the Cascade Range and the Coastal California Mountains, are projected to experience greater annual fire occurrence probabilities, with relative increases of 122% and 67%, respectively, under RCP 8.5 compared to increases of 63% and 38% under RCP 4.5. Regions not currently associated with frequently occurring wildfires, such as New England and the Great Lakes, are projected to experience a doubling of occurrence probabilities by 2100 under RCP 8.5. This high resolution, continental-scale modeling study of climate change impacts on potential fire probability accounts for shifting background environmental conditions across regions that will interact with topographic drivers to significantly alter future fire probabilities. The ensemble modeling approach presents a useful planning tool for mitigation and adaptation strategies in regions of increasing wildfire risk.

Science of the Total Environment

Forage fish species prefer habitat within designated offshore wind energy areas on the US Northeast Shelf ecosystem

As the world develops sources of renewable energy, there is an intensifying interest in offshore wind energy production. The Northeast U.S. Continental Shelf (NES) ecosystem has favorable wind dynamics, with active development of wind energy. In this study, we present species distribution models that consider both occupancy and biomass responses for a broad spectrum of fish and macroinvertebrate taxa ( n = 177). Building upon prior analyses, habitat was differentiated into overall and core habitats based on statistical distributions of habitat scores. Overall habitat was used to show each species' regional distribution based on fishery-independent survey captures between 1976 and 2019, whereas core habitat represented where the focus of the species' abundance was located as a subset of overall habitat. Wind energy developments may modify the water column in ways that impact lower-trophic-level productivity; therefore, added attention was given to the response of forage species. Over 20% of species showed preferential use of putative and potential wind development areas, including a disproportionate number of forage taxa. Principal usage varied by season, with forage species like Atlantic Menhaden Brevoortia tyrannus and Atlantic Mackerel Scomber scombrus preferentially using the lease areas in spring and Round Herring Etrumeus teres and longfin inshore squid Doryteuthis pealeii using lease areas in autumn. For species with relatively low usage of the lease areas, there was a tendency for the usage related to overall habitat to be lower than usage for core habitat; in contrast, for species with high usage of the lease areas, that usage was higher for overall habitat than for core habitat. The area of habitat tended to have positive trends across species, with these positive trends being disproportionately higher among forage taxa. These results frame the importance of wind lease areas for species in the NES, particularly forage taxa that fulfill many important ecological functions.

Marine and Coastal Fisheries: Dynamics, Management

Pilot framework for fish habitat assessments across tidal and non tidal waters in the Patuxent River Basin

As part of the 2014 Chesapeake Bay Watershed Agreement, all Bay States and the District of Columbia have committed to improving the condition of the Bay, which includes a goal to achieve sustainable fisheries. One outcome under that broad goal is improved effectiveness of fish habitat conservation and preservation efforts. In support of that outcome, the U.S. Geological Survey Eastern Ecological Science Center (USGS-EESC) and the National Oceanic and Atmospheric Association’s National Centers for Coastal Ocean Science (NOAA-NCCOS) are actively developing datasets, methods, and analyses to conduct fish habitat assessments in the Chesapeake Bay watershed, guided by recommendations from a regional stakeholder workshop held by the Chesapeake Bay Program’s (CBP) Fish Habitat Action Team (FHAT) in 2018. The joint USGS and NOAA team has been collaborating on methods for conducting inland and estuarine assessments and exploring whether a seamless headwater to estuary assessment could be developed. The goals of this assessment are to benefit both State and Federal fisheries managers, help advance fisheries science, and provide beneficial information for the public. While past national and regional assessments (e.g. the National Fish Habitat Partnership National Assessment) treated inland and estuarine fish habitat conditions separately due to differences in environments, GIS data representation, and data availability, a seamless habitat assessment could be of value for a broad range of stakeholders as many fish species, several of which are invasive or under federal jurisdiction, use habitats across both inland and estuarine waters. This project developed a pilot framework, explored and tested methods necessary for a finer scale, seamless assessment across both inland and estuarine waters, and demonstrated its use. Although there was interest by the CBP FHAT for the generation of a Baywide fish habitat assessment that spanned tidal salt, tidal fresh, warm non-tidal and cold non-tidal waters, there are a myriad of implementation details and considerations around conducting a Baywide assessment across all four of these general habitat areas. Therefore, the practical need to conduct a tributary-specific pilot assessment arose. At the beginning of this pilot process, members of the FHAT were presented with a decision matrix to choose a study basin using factors such as data availability and tributary size. FHAT members chose the Patuxent River basin, which has been relatively well sampled and studied. Several spatial frameworks were considered before selection of an inclusive gridded framework for summary and analysis that represented inland drainage networks and landscape influences as well as estuarine bathymetry. A suite of landscape and in-water stressor variables were summarized into the framework and were largely generalized over time. In order to assess the viability of the framework, we chose to use species distribution modeling for each of the species to test the framework’s ability to predict habitat use of non-tidal resident, estuarine resident, and migratory species. Tessellated darter (Etheostoma olmstedi), American eel (Anguilla rostrata), and white perch (Morone americana) were chosen as illustrative fish species based on data availability, and differences in life history and habitat use. A nested modeling approach, which involved successive model runs at multiple scales (1000m, 100m, and 10m raster grids) was developed to examine differences in variable importance at different spatial scales and to enhance modeling efficiency. For white perch, a complementary modeling analysis was performed for variables available only in estuarine waters. For all testing, an ensemble modeling approach was conducted, using a suite of potential statistical techniques driven by model strength and variable predictive power. The statistical testing that we conducted was intended only to test the framework and modeling approach, and not to definitively predict all habitats where specific fish species might be present. The modeling we conducted to test the framework did have some limitations. For example, the spatial distribution of favorable habitat areas for white perch was likely influenced by the predominance of fish survey locations near the center channel of the river and the use of generalized in-water conditions. For all species, the use of juvenile and adult fish survey data limits the estimation of habitat use to those life stages. Despite such limitations of the data inputs and modeling approach, we found the framework could seamlessly predict fish habitat distribution across freshwater and tidal environments and integrate the influence of landscape stressors with local in-water factors. The developed framework presented to the Sustainable Fisheries Goal Implementation Team (GIT) and FHAT is informative and could potentially be used for other modeling applications in the Chesapeake Bay watershed and elsewhere. In particular the framework and modeling approach lend themselves to evaluating living resource distributions and underlying habitat conditions in shallow tidal waters and beyond, as recommended by the recent Comprehensive Evaluation of System Response (CESR) report from the Chesapeake Bay Program.

Maryland

Tumor prevalence and biomarkers of genotoxicity in brown bullhead (Ameiurus nebulosus) in Chesapeake Bay tributaries

We surveyed four Chesapeake Bay tributaries for skin and liver tumors in brown bullhead ( Ameiurus nebulosus ). We focused on the South River, where the highest skin tumor prevalence (53%) in the Bay watershed had been reported. The objectives were to 1) compare tumor prevalence with nearby rivers (Severn and Rhode) and a more remote river (Choptank); 2) investigate associations between tumor prevalence and polynuclear aromatic hydrocarbons (PAHs) and alkylating agents; and 3) statistically analyze Chesapeake Bay bullhead tumor data from 1992 through 2008. All four South River collections exhibited high skin tumor prevalence (19% to 58%), whereas skin tumor prevalence was 2%, 10%, and 52% in the three Severn collections; 0% and 2% in the Choptank collections; and 5.6% in the Rhode collection. Liver tumor prevalence was 0% to 6% in all but one South River collection (20%) and 0% to 6% in the three other rivers. In a subset of samples, PAH-like biliary metabolites and 32 P-DNA adducts were used as biomarkers of exposure and response to polycyclic aromatic compounds (PACs). Adducts from alkylating agents were detected as O6-methyl-2′-deoxyguanosine (O6Me-dG) and O6-ethyl-2′-deoxyguanosine (O6Et-dG) modified DNA. Bullheads from the contaminated Anacostia River were used as a positive control for DNA adducts. 32 P-DNA adduct concentrations were significantly higher in Anacostia bullhead livers compared with the other rivers. We identified alkyl DNA adducts in bullhead livers from the South and Anacostia, but not the Choptank. Neither the PAH-like bile metabolite data, sediment PAH data, nor the DNA adduct data suggest an association between liver or skin tumor prevalence and exposure to PACs or alkylating agents in the South, Choptank, Severn, or Rhode rivers. Logistic regression analysis of the Chesapeake Bay database revealed that sex and length were significant covariates for liver tumors and length was a significant covariate for skin tumors.

Chesapeake Bay watershed, South River

Surrounding land cover types as predictors of palustrine wetland vegetation quality in conterminous USA

The loss of wetland habitats and their often-unique biological communities is a major environmental concern. We examined vegetation data obtained from 380 wetlands sampled in a statistical survey of wetlands in the USA. Our goal was to identify which surrounding land cover types best predict two indices of vegetation quality in wetlands at the regional scale. We considered palustrine wetlands in four regions (Coastal Plains, North Central East, Interior Plains, and West) in which the dominant vegetation was emergent, forested, or scrub-shrub. For each wetland, we calculated weighted proportions of eight land cover types surrounding the area in which vegetation was assessed, in four zones radiating from the edge of the assessment area to 2 km. Using Akaike's Information Criterion, we determined the best 1-, 2- and 3-predictor models of the two indices, using the weighted proportions of the land cover types as potential predictors. Mean values of the two indices were generally higher in the North Central East and Coastal Plains than the other regions for forested and emergent wetlands. In nearly all cases, the best predictors of the indices were not the dominant surrounding land cover types. Overall, proportions of forest (positive effect) and agriculture (negative effect) surrounding the assessment area were the best predictors of the two indices. One or both of these variables were included as predictors in 65 of the 72 models supported by the data. Wetlands surrounding the assessment area had a positive effect on the indices, and ranked third (33%) among the predictors included in supported models. Development had a negative effect on the indices and was included in only 28% of supported models. These results can be used to develop regional management plans for wetlands, such as creating forest buffers around wetlands, or to conserve zones between wetlands to increase habitat connectivity.

Science of the Total Environment

Application of health-based screening levels to ground-water quality data in a state-scale pilot effort

A state-scale pilot effort was conducted to evaluate a Health-Based Screening Level (HBSL) approach developed for communicating findings from the U.S. Geological Survey (USGS) National Water-Quality Assessment Program in a human-health context. Many aquifers sampled by USGS are used as drinking-water sources, and water-quality conditions historically have been assessed by comparing measured contaminant concentrations to established drinking-water standards and guidelines. Because drinking-water standards and guidelines do not exist for many analyzed contaminants, HBSL values were developed collaboratively by the USGS, U.S. Environmental Protection Agency (USEPA), New Jersey Department of Environmental Protection, and Oregon Health & Science University, using USEPA toxicity values and USEPA Office of Water methodologies. The main objective of this report is to demonstrate the use of HBSL approach as a tool for communicating water-quality data in a human-health context by conducting a retrospective analysis of ground-water quality data from New Jersey. Another important objective is to provide guidance on the use and interpretation of HBSL values and other human-health benchmarks in the analyses of water-quality data in a human-health context. Ground-water samples collected during 1996-98 from 30 public-supply, 82 domestic, and 108 monitoring wells were analyzed for 97 pesticides and 85 volatile organic compounds (VOCs). The occurrence of individual pesticides and VOCs was evaluated in a human-health context by calculating Benchmark Quotients (BQs), defined as ratios of measured concentrations of regulated compounds (that is, compounds with Federal or state drinking-water standards) to Maximum Contaminant Level (MCL) values and ratios of measured concentrations of unregulated compounds to HBSL values. Contaminants were identified as being of potential human-health concern if maximum detected concentrations were within a factor of 10 of the associated MCL or HBSL (that is, maximum BQ value (BQmax) greater than or equal to 0.1) in any well type (public supply, domestic, monitoring). Most (57 of 77) pesticides and VOCs with human-health benchmarks were detected at concentrations well below these levels (BQmax less than 0.1) for all three well types; however, BQmax values ranged from 0.1 to 3,000 for 6 pesticides and 14 VOCs. Of these 20 contaminants, one pesticide (dieldrin) and three VOCs (1,2-dibromoethane, tetrachloroethylene, and trichloroethylene) both (1) were measured at concentrations that met or exceeded MCL or HBSL values, and (2) were detected in more than 10 percent of samples collected from raw ground water used as sources of drinking water (public-supply and (or) domestic wells) and, therefore, are particularly relevant to human health. The occurrence of multiple pesticides and VOCs in individual wells also was evaluated in a human-health context because at least 53 different contaminants were detected in each of the three well types. To assess the relative human-health importance of the occurrence of multiple contaminants in different wells, the BQ values for all contaminants in a given well were summed. The median ratio of the maximum BQ to the sum of all BQ values for each well ranged from 0.83 to 0.93 for all well types, indicating that the maximum BQ makes up the majority of the sum for most wells. Maximum and summed BQ values were statistically greater for individual public-supply wells than for individual domestic and monitoring wells. The HBSL approach is an effective tool for placing water-quality data in a human-health context. For 79 of the 182 compounds analyzed in this study, no USEPA drinking-water standards or guidelines exist, but new HBSL values were calculated for 39 of these 79 compounds. The new HBSL values increased the number of detected pesticides and VOCs with human-health benchmarks from 65 to 77 (of 97 detected compounds), thereby expanding the basis for interpreting contaminant-occu

Scientific Investigations Report

Archive of digital boomer and CHIRP seismic reflection data collected during USGS field activity 08LCA03 in Lake Panasoffkee, Florida, May 2008

In May of 2008, the U.S. Geological Survey (USGS) conducted geophysical surveys in Lake Panasoffkee, located in central Florida, as part of the USGS Lakes and Coastal Aquifers (LCA) study. This report serves as an archive of unprocessed digital boomer and Compressed High Intensity Radar Pulse (CHIRP)* seismic reflection data, trackline maps, navigation files, Field Activity Collection System (FACS) logs, Geographic Information System (GIS) files, and formal Federal Geographic Data Committee (FGDC) metadata. Filtered and gained (a relative increase in signal amplitude) digital images of the seismic profiles and geospatially corrected interactive profiles are also provided. Refer to the Acronyms page for expansions of acronyms and abbreviations used in this report. *Due to poor data acquisition conditions associated with the lake bottom sediments, only two CHIRP tracklines were collected during this field activity. The archived trace data are in standard Society of Exploration Geophysicists (SEG) SEG-Y format (Barry and others, 1975) and may be downloaded and processed with commercial or public domain software such as Seismic Unix (SU). Example SU processing scripts and USGS software for viewing the SEG-Y files (Zihlman, 1992) are provided. The USGS Florida Integrated Science Center (FISC) - St. Petersburg assigns a unique identifier to each cruise or field activity. For example, 08LCA03 tells us the data were collected in 2008 for the Lakes and Coastal Aquifers (LCA) study and the data were collected during the third field activity for that study in that calendar year. Refer to http://walrus.wr.usgs.gov/infobank/programs/html/definition/activity.html for a detailed description of the method used to assign the field activity ID. The naming convention used for each seismic line is as follows: yye##a, where 'yy' are the last two digits of the year in which the data were collected, 'e' is a 1-letter abbreviation for the equipment type (for example, b for boomer and c for CHIRP), '##' is a 2-digit number representing a specific track, and 'a' is a letter representing the section of a line if recording was prematurely terminated or rerun for quality or acquisition problems. The boomer plate is an acoustic energy source that consists of capacitors charged to a high voltage and discharged through a transducer in the water. The transducer is towed on a sled floating on the water surface and, when discharged, emits a short acoustic pulse, or shot, which propagates through the water, sediment column, or rock beneath. The acoustic energy is reflected at density boundaries (such as the seafloor, sediment, or rock layers beneath the seafloor), detected by the receiver, and recorded by a PC-based seismic acquisition system. This process is repeated at timed intervals (for example, 0.5 s) and recorded for specific intervals of time (for example, 100 ms). In this way, a two-dimensional (2-D) vertical profile of the shallow geologic structure beneath the ship track is produced. Figure 1 displays the boomer acquisition geometry. The EdgeTech SB-424 CHIRP system used for this survey has a vertical resolution of 4 - 8 cm, a penetration depth that is usually less than 2 m beneath the seafloor, and uses a signal of continuously varying frequency. The towfish is a sound source and receiver, which is typically towed 2 - 5 m above the seafloor. The acoustic energy is reflected at density boundaries (such as the seafloor or sediment layers beneath the seafloor), detected by a receiver, and recorded by a PC-based seismic acquisition system. This process is repeated at timed intervals (for example, 0.125 s) and recorded for specific intervals of time (for example, 50 ms); the resulting profile is a two-dimensional vertical image of the shallow geologic structure beneath the ship track. Figure 2 displays the acquisition geometry for the CHIRP system. Refer to table 1 for a summary of acquisition parameters and table 2 for trackline statistics.

Florida

Spatiotemporal patterns of mercury accumulation in lake sediments of western North America

For the Western North America Mercury Synthesis, we compiled mercury records from 165 dated sediment cores from 138 natural lakes across western North America. Lake sediments are accepted as faithful recorders of historical mercury accumulation rates, and regional and sub-regional temporal and spatial trends were analyzed with descriptive and inferential statistics. Mercury accumulation rates in sediments have increased, on average, four times (4×) from 1850 to 2000 and continue to increase by approximately 0.2 μg/m 2 per year. Lakes with the greatest increases were influenced by the Flin Flon smelter, followed by lakes directly affected by mining and wastewater discharges. Of lakes not directly affected by point sources, there is a clear separation in mercury accumulation rates between lakes with no/little watershed development and lakes with extensive watershed development for agricultural and/or residential purposes. Lakes in the latter group exhibited a sharp increase in mercury accumulation rates with human settlement, stabilizing after 1950 at five times (5×) 1850 rates. Mercury accumulation rates in lakes with no/little watershed development were controlled primarily by relative watershed size prior to 1850, and since have exhibited modest increases (in absolute terms and compared to that described above) associated with (regional and global) industrialization. A sub-regional analysis highlighted that in the ecoregion Northwestern Forest Mountains, <1% of mercury deposited to watersheds is delivered to lakes. Research is warranted to understand whether mountainous watersheds act as permanent sinks for mercury or if export of “legacy” mercury (deposited in years past) will delay recovery when/if emissions reductions are achieved.

Alaska, Alberta, Arizona, British Columbia, Califo

A hybrid machine learning model to predict and visualize nitrate concentration throughout the Central Valley aquifer, California, USA

Intense demand for water in the Central Valley of California and related increases in groundwater nitrate concentration threaten the sustainability of the groundwater resource. To assess contamination risk in the region, we developed a hybrid, non-linear, machine learning model within a statistical learning framework to predict nitrate contamination of groundwater to depths of approximately 500 m below ground surface. A database of 145 predictor variables representing well characteristics, historical and current field and landscape-scale nitrogen mass balances, historical and current land use, oxidation/reduction conditions, groundwater flow, climate, soil characteristics, depth to groundwater, and groundwater age were assigned to over 6000 private supply and public supply wells measured previously for nitrate and located throughout the study area. The boosted regression tree (BRT) method was used to screen and rank variables to predict nitrate concentration at the depths of domestic and public well supplies. The novel approach included as predictor variables outputs from existing physically based models of the Central Valley. The top five most important predictor variables included two oxidation/reduction variables (probability of manganese concentration to exceed 50 ppb and probability of dissolved oxygen concentration to be below 0.5 ppm), field-scale adjusted unsaturated zone nitrogen input for the 1975 time period, average difference between precipitation and evapotranspiration during the years 1971–2000, and 1992 total landscape nitrogen input. Twenty-five variables were selected for the final model for log-transformed nitrate. In general, increasing probability of anoxic conditions and increasing precipitation relative to potential evapotranspiration had a corresponding decrease in nitrate concentration predictions. Conversely, increasing 1975 unsaturated zone nitrogen leaching flux and 1992 total landscape nitrogen input had an increasing relative impact on nitrate predictions. Three-dimensional visualization indicates that nitrate predictions depend on the probability of anoxic conditions and other factors, and that nitrate predictions generally decreased with increasing groundwater age.

California

Diving deeper into seep distribution along the Cascadia Convergent Margin, USA

Previous margin-wide studies of methane seep distribution along the Cascadia Subduction Zone indicate peaks in seep density within the landward limit of the of gas hydrate stability zone (GHSZ; ≤500 m depth), suggesting a link between current ocean warming, acceleration of hydrate dissociated, and methane emissions. This inferred connection, however, may not account for regional geologic and/or structural complexities driving methane seepage. Expanding upon an existing seep database by adding new seeps data, we conducted statistical and spatial analyses to determine margin-wide distribution trends and offer a tectonic framework for understanding the tendency toward non-normality and spatial clustering. We then highlight the role of local-scale drivers of seep formation in addition to the first-order tectonic framework, using systematic geologic/geomorphic characterization of seep emission sites in southern Cascadia and case studies using meta-attribute analysis of seismic reflection data.. Seep distribution along the margin is non-random, but instead of clustering along the 500-m isobath, regions of high seep density occur in canyons and topographic highs. New findings from this study conclude that co-location of the outer arc high (OAH) and the landward limit of the GHSZ may explain high concentrations of seeps where deformation is the greatest and hydrates are unstable. Detailed analysis of the spatial relationships between seep sites and geologic-geomorphic features in southern Cascadia reveal a link between seeps and anticlines, with 52% of the seeps found in association with anticlines, 36% found at faults, 16% associated with canyons, and 11% at seafloor failure scarps. Given that a majority of anticlines are located along or seaward of the OAH in the actively deforming outer wedge, we suggest that the location of the OAH is a primary structural control on seep distribution. This scenario is supported by neural network analysis of multichannel seismic data revealing zones of probable fluid migration along vertical pipes, faults, and chimneys in the vicinity of active seep sites on anticlines. Determining linkages between seeps and submarine tectonic geomorphology is a crucial first step for understanding and forecasting the distribution of methane seepage, but also a necessity for evaluating causal relationships between ocean warming and gas-hydrate stability.

California, Oregon, Washington

Concentrations, loads, and associated trends of nutrients entering the Sacramento-San Joaquin Delta, California

Statistical modeling of water-quality data collected at the Sacramento River at Freeport and San Joaquin River near Vernalis, California, USA, was used to examine trends in concentrations and loads of various forms of dissolved and particulate nitrogen and phosphorus that entered the Sacramento–San Joaquin River Delta (Delta) from upstream sources between 1970 and 2019. Ammonium concentrations and loads decreased at the Sacramento River site from the mid-1970s through 1990 because of the consolidation of wastewater treatment and continuously reduced from the mid-1970s to 2019 at the San Joaquin River site. Current ammonium concentrations are mostly below 4 µM (0.056 mg N L –1 ) at both sites, a concentration above which reductions in phytoplankton productivity or changes in algal species composition may occur. The Sacramento River at Freeport site is located upstream of the Sacramento Regional County Sanitation District’s treatment facility’s discharge point; nutrient water quality there is representative of upstream sources. Inorganic nitrogen (nitrate plus ammonium) concentrations and loading differed at both sites. At the Sacramento River location, concentrations decrease in the summer agricultural season, reducing the molar ratios of nitrogen to phosphorus. In contrast, inorganic nitrogen concentrations increase in the San Joaquin River during the agricultural season as a result of irrigation runoff, increasing the molar ratio of nitrogen to phosphorus. This increase suggests a possible nitrogen limitation in the northern Delta and a phosphorus limitation in the southern Delta, as indicated by the molar ratios of bioavailable nitrogen to bioavailable phosphorus. Planned upgrades to the Sacramento Regional Wastewater Treatment Plant (SRWTP) will reduce inorganic nitrogen inputs to the northern Delta. Consequently, the supply of bioavailable nitrogen throughout the upper estuary should diminish. Source modeling of nitrogen and phosphorus identifies agriculture, atmospheric deposition, and wastewater effluent as sources of total nitrogen in the Central Valley. In contrast, geologic sources, agriculture, and wastewater discharge are the primary sources of phosphorus.

California

Analysis and review of fishery-dependent data for Hawaiian nearshore noncommercial fisheries

Noncommercial, shore-based fisheries provide economic, social, and cultural services to communities throughout the Hawaiian Islands. The State of Hawai‘i Department of Land and Natural Resources (DLNR), Division of Aquatic Resources (DAR) routinely conducts surveys to monitor noncommercial fisheries such that estimates of fishing effort and catch by gear type can be generated and used to implement more sustainable management practices. DAR executes both the Hawai‘i Marine Recreational Fishery Survey (HMRFS), a nationally standardized survey that focuses on intercepting fishers at access points (i.e., boat ramps) across the main Hawaiian Islands, and a set of roving creel surveys on O‘ahu, Maui Nui, and Kaua‘i that observe and intercept fishers at locations along the shoreline outside of those targeted by HMRFS. The latter set of creel surveys were designed to complement HMRFS by expanding its geographic coverage and thus providing a more representative picture of noncommercial fishing in Hawai‘i. Sustainable management priorities set by DAR rely on the availability of statewide, fishery- dependent data. Thus, we collate information from island-based roving creel surveys into a cohesive Statewide Creel Survey Database. Further, we provide preliminary analyses and describe ways that surveys could be streamlined to improve future data collection, analysis, and utility. In so doing, we synthesize the most detailed information to-date about noncommercial shore-based fisheries of Hawai‘i. The unprecedented spatial and temporal coverage of DAR’s dataset reveals the value of their survey efforts over the last decade to address fishery management needs. Our primary objectives, results, and conclusions are summarized below: 1) Integrate DAR roving creel survey data from different islands into a single Statewide Creel Survey Dataset (Chapter II). We describe the collation of creel survey data from O‘ahu, Maui Nui, and Kaua‘i into a statewide dataset. We also offer ways in which these surveys could be streamlined to meet the needs of managers and decision makers. Briefly, these are to create a statewide strategic plan, standardize the execution of standard operating procedures, centralize the creel survey database and associated metadata, and consider using technology that improves the data pipeline, including transitioning from paper-based to electronic systems for data entry and processing. 2) Assess whether the new Statewide Creel Survey Dataset can provide inputs for length-based stock assessments (Chapter III). Only on Maui were interviews conducted with associated catch data. There was reasonably high taxonomic coverage (42 species from 186 interviews with 310 fishers), but low sample sizes for nearly all species precluded the development of length-based stock assessments. We provide summary statistics from the existing data and briefly discuss how technologies could be used to automate analysis of images of noncommercial catch. 3) Analyze the Statewide Creel Survey Dataset for spatial and temporal patterns in fishing effort (Chapter IV): a. Visualizing noncommercial fishing pressure . We found that fishing effort (mean number of fishers observed per survey event at a site) on O'ahu was over three times greater than that recorded during similar surveys conducted on Maui or Kaua'i. We create maps that display the distribution of angling and spearfishing effort around each of the three islands. b. Factors that predict fishing “hotspots” around Maui . Fishing effort on Maui was associated with areas with more wave power and less parking availability. There were half as many fishers in areas with parking lots than in areas with parking on the road shoulder only. c. Changes in fishing effort during the COVID-19 pandemic . There was no change in fishing effort on O‘ahu during the first year of the pandemic, but there was a 20% decline in year 2 and a 33% decline in year 3, both in comparison to pre-pandemic levels. Pre-pandemic creel survey data were unavailable for Maui and Kaua‘i, but fishing effort on these islands also declined as the pandemic progressed at similar or greater rates than those observed on O‘ahu. 4) Quantify potential bias in survey methods by experimentally deriving fisher detection probabilities of shore-based and drone-based surveys (Chapter V): a. Shore-based surveys . We conducted roving creel surveys for four months at three locations around Hilo Bay, designed to emulate and estimate the efficacy of DAR standard operating procedures. There was high agreement between paired observers in counting fishers, leading to near-perfect detection probabilities of both anglers (94%) and spearfishers (97%), but relatively low agreement and detection probabilities of other fishers (throw net, ‘opihi picking, etc.) (52%). b. Drone-based surveys . We used an unmanned aerial vehicle (UAV; operated by DAR staff) to collect imagery of fishers along the Hilo Bay shoreline. We used still images and video clips (with known fishing activity) to build an online survey that was distributed to DAR and HCFRU personnel, asking them to count and categorize resource users as a snorkeler, spearfisher, angler, or other fisher. Only 40.0% of the responses correctly counted and categorized resource users in the image. Anglers were correctly identified and enumerated in 90.0% of the responses, but the correct response rates of the other three user categories ranged from 67.8% – 79.4%. Snorkelers and anglers tended to be undercounted while spearfishers and other fishers were overcounted. 5) Review the potential for incorporating emerging technologies that will improve, augment, and evolve creel survey data collection, especially for spearfishing (Chapter VI). Within the context of monitoring shore-based noncommercial fishing, we review the use of electronic data entry/processing systems with geospatial and image capabilities, field cameras, drones, smart buoys, citizen science apps, data mining social media, artificial intelligence and machine learning. We highlight several of the challenges and considerations when implementing these technologies into creel surveys and provide a synthesis of options that could be used to better estimate spearfishing. The general conclusion of this assessment is that the DAR roving creel survey program is collecting valuable data that supplement the existing HMRFS efforts. However, there are a number of areas that could be improved to make these efforts a more effective tool for decision-making processes in resource management and conservation: 1) Establishment of clear statewide and island objectives for the Statewide Creel Survey Dataset. Currently, data collection efforts are focused towards addressing a very broad purpose – supplementing the HMRFS data collection efforts. However, the results of the preliminary analyses conducted as part of this project suggest that the data could be used to address other areas of need if these objectives were clearly defined. Further, the design of the creel survey would benefit from greater standardization of survey protocols between islands and an effort to define a) the acceptable margins of error associated with the estimates generated by these data and b) the minimum level of change that the surveys would need to detect to be useful to managers. 2) Centralization of data entry, data quality assessment, and data accessibility. Currently, each DAR office manages data entry, checks the data for errors, and is responsible for managing and storing the data. Instituting a centralized data entry system, particularly an online database that can receive survey data from tablets or smartphones running a standardized data collection application would improve efficiency, reduce data entry errors, and accelerate the availability of data to managers. A substantial amount of time and effort from the project described in this report was devoted to checking the dataset for errors. The development and application of data quality assurance protocols would ensure that the data are reliable and available in a timely fashion to support management decisions. 3) Address lingering questions regarding the efficacy of current survey protocols to capture and characterize the spearfishing component of the noncommercial fishery. The results presented in the report suggest that the current creel survey protocols do a good job detecting spearfishers when present but are not capturing sufficient data about their catch or total effort. There are also questions remaining as to whether the survey times and sites are sufficiently capturing the behavior of spearfishers in Hawai‘i. A more thorough assessment – whether through additional research, alteration of survey design, or review of data by representatives of the spearfishing community – would provide insight on how to use the Statewide Creel Survey Database to inform management of spearfishing. 4) Investigate the integration of technological advancements into the creel survey methods. As priorities and needs are developed and formalized, it would be valuable to consider how various technological advancements might enhance and streamline data collection or open new avenues of inquiry.

Hawaii

An approach for decomposing river water-quality trends into different flow classes

A number of statistical approaches have been developed to quantify the overall trend in river water quality, but most approaches are not intended for reporting separate trends for different flow conditions. We propose an approach called FN 2Q , which is an extension of the flow-normalization (FN) procedure of the well-established WRTDS (“Weighted Regressions on Time, Discharge, and Season”) method. The FN 2Q approach provides a daily time series of low-flow and high-flow FN flux estimates that represent the lower and upper half of daily riverflow observations that occurred on each calendar day across the period of record. These daily estimates can be summarized into any time period of interest (e.g., monthly, seasonal, or annual) for quantifying trends. The proposed approach is illustrated with an application to a record of total nitrogen concentration (632 samples) collected between 1985 and 2018 from the South Fork Shenandoah River at Front Royal, Virginia (USA). Results show that the overall FN flux of total nitrogen has declined in the period of 1985–2018, which is mainly attributable to FN flux decline in the low-flow class. Furthermore, the decline in the low-flow class was highly correlated with wastewater effluent loads, indicating that the upgrades of treatment technology at wastewater treatment facilities have likely led to water-quality improvement under low-flow conditions. The high-flow FN flux showed a spike around 2007, which was likely caused by increased delivery of particulate nitrogen associated with sediment transport. The case study demonstrates the utility of the FN 2Q approach toward not only characterizing the changes in river water quality but also guiding the direction of additional analysis for capturing the underlying drivers. The FN 2Q approach (and the published code) can easily be applied to widely available river monitoring records to quantify water-quality trends under different flow conditions to enhance understanding of river water-quality dynamics.

Virginia

Estimating black bear density in New Mexico using noninvasive genetic sampling coupled with spatially explicit capture-recapture methods

During the 2004–2005 to 2015–2016 hunting seasons, the New Mexico Department of Game and Fish (NMDGF) estimated black bear abundance (Ursus americanus) across the state by coupling density estimates with the distribution of primary habitat generated by Costello et al. (2001). These estimates have been used to set harvest limits. For example, a density of 17 bears/100 km2 for the Sangre de Cristo and Sacramento Mountains and 13.2 bears/100 km2 for the Sandia Mountains were used to set harvest levels. The advancement and widespread acceptance of non-invasive sampling and mark-recapture methods, prompted the NMDGF to collaborate with the New Mexico Cooperative Fish and Wildlife Research Unit and New Mexico State University to update their density estimates for black bear populations in select mountain ranges across the state. We established 5 study areas in 3 mountain ranges: the northern (NSC; sampled in 2012) and southern Sangre de Cristo Mountains (SSC; sampled in 2013), the Sandia Mountains (Sandias; sampled in 2014), and the northern (NSacs) and southern Sacramento Mountains (SSacs; both sampled in 2014). We collected hair samples from black bears using two concurrent non-invasive sampling methods, hair traps and bear rubs. We used a gender marker and a suite of microsatellite loci to determine the individual identification of hair samples that were suitable for genetic analysis. We used these data to generate mark-recapture encounter histories for each bear and estimated density in a spatially explicit capture-recapture framework (SECR). We constructed a suite of SECR candidate models using sex, elevation, land cover type, and time to model heterogeneity in detection probability and the spatial scale over which detection probability declines. We used Akaike’s Information Criterion corrected for small sample size (AICc) to rank and select the most supported model from which we estimated density. We set 554 hair traps, 117 bear rubs and collected 4,083 hair samples. We identified 725 (367 M, 358 F) individuals; the sex ratio for each study area was approximately equal. Our density estimates varied within and among mountain ranges with an estimated density of 21.86 bears/100 km2 (95% CI: 17.83 – 26.80) for the NSC, 19.74 bears/100 km2 (95% CI: 13.77 – 28.30) in the SSC, 25.75 bears/100 km2 (95% CI: 13.22 – 50.14) in the Sandias, 21.86 bears/100 km2 (95% CI: 17.83 – 26.80) in the NSacs, and 16.55 bears/100 km2 (95% CI: 11.64 – 23.53) in the SSacs. Overall detection probability for hair traps and bear rubs, combined, was low across all study areas and ranged from 0.00001 to 0.02. We speculate that detection probabilities were affected by failure of some hair samples to produce a complete genotype due to UV degradation of DNA, and our inability to set and check some sampling devices due to wildfires in the SSC. Ultraviolet radiation levels are particularly high in New Mexico compared to other states where NGS methods have been used because New Mexico receives substantial amounts of sunshine, is relatively high in elevation (1,200 m – 4,000 m), and is at a lower latitude. Despite these sampling difficulties, we were able to produce density estimates for New Mexico black bear populations with levels of precision comparable to estimated black bear densities made elsewhere in the U.S. Our ability to generate reliable black bear density estimates for 3 New Mexico mountain ranges is attributable to our use of a statistically robust study design and analytical method. There are multiple factors that need to be considered when developing future SECR-based density estimation projects. First, the spatial extent of the population of interest and the smallest average home range size must be determined; these will dictate size of the trapping array and spacing necessary between hair traps. The number of technicians needed and access to the study areas will also influence configuration of the trapping array. We believe shorter sampling occasions could be implemented to reduce degradation of DNA due to UV radiation; this might help increase amplification rates and thereby increase both the number of unique individuals identified and the number of recaptures, improving the precision of the density estimates. A pilot study may be useful to determine the length of time hair samples can remain in the field prior to collection. In addition, researchers may consider setting hair traps and bear rubs in more shaded areas (e.g., north facing slopes) to help reduce exposure to UV radiation. To reduce the sampling interval it will be necessary to either hire more field personnel or decrease the number of hair traps per sampling session. Both of these will enhance detection of long-range movement events by individual bears, increase initial capture and recapture rates, and improve precision of the parameter estimates. We recognize that all studies are constrained by limited resources, however, increasing field personnel would also allow a larger study area to be sampled or enable higher trap density. In conclusion, we estimated the density of black bears in 5 study areas within 3 mountains ranges of New Mexico. Our estimates will aid the NMDGF in setting sustainable harvest limits. Along with estimates of density, information on additional demographic rates (e.g., survival rates and reproduction) and the potential effects that climate change and future land use may have on the demography of black bears may also help inform management of black bears in New Mexico, and may be considered as future areas for research.

New Mexico

A Framework for Long-term Ecological Monitoring in Olympic National Park: Prototype for the Coniferous Forest Biome

This report is the result of a five-year collaboration between scientists of the U.S. Geological Survey Forest and Rangeland Ecosystem Science Center, Olympic Field Station, and the natural resources staff of Olympic National Park to develop a comprehensive strategy for monitoring natural resources of Olympic National Park. Olympic National Park is the National Park Serviceʼs prototype monitoring park, representing parks in the coniferous forest biome. Under the umbrella of the National Park Serviceʼs prototype parks program, U.S. Geological Survey and Olympic National Park staffs are obligated to: develop strategies and designs for monitoring the long-term health and integrity of national park ecosystems with a significant coniferous forest component. design exportable monitoring protocols that can be used by other parks within the coniferous forest biome (i.e., parks having similar environments), and create a demonstration area and ʻcenter of excellenceʼ for assisting other parks in developing ecological monitoring programs. Olympic National Park is part of the North Coast and Cascades Network, a network of seven Pacific Northwestern park units created recently by the National Park Serviceʼs Inventory and Monitoring Program to extend the monitoring of ʻvital signsʼ of park health to all National Park Service units. It is our intent and hope that the monitoring strategies and conceptual models described here will meet the overall purpose of the prototype parks monitoring program in proving useful not only to Olympic National Park, but also to parks within the North Coast and Cascades Network and elsewhere. Part I contains the conceptual design and sampling framework for the prototype long-term monitoring program in Olympic National Park. In this section, we explore key elements of monitoring design that help to ensure the spatial, ecological, and temporal integration of monitoring program elements and discuss approaches used to design an ecosystem-based monitoring program. Basic monitoring components include ecosystem drivers, (e.g., climate, atmospheric inputs, human pressures), indicators of ecosystem integrity (e.g., biogeochemical indicators), known threats (e.g., impacts of introduced mountain goats), and focal or ʻkeyʼ species (e.g., rare or listed species, Roosevelt elk). Monitoring system drivers and key indicators of ecosystem integrity provide the long-term baseline needed to judge what constitutes ʻunnaturalʼ variation in park resources and provide the earliest possible warning of unacceptable change. Monitoring effects of known threats and the status of focal species will provide information useful to park managers for dealing with current park issues. In Part I we describe the process of identifying potential indicators of ecological condition and present conceptual models of park ecosystems. In addition we report results from several workshops held in conjunction with Olympic National Park aimed at identifying potential indicators of change in the parkʼs ecosystem. First, we describe the responses of Olympic National Park staff to the generic question, “What is the most important resource to monitor in Olympic National Park and why?” followed by the responses from resource and land managers from areas adjoining the park. We also catalogue the responses of various expert groups that we asked to help identify the most appropriate system drivers and indicators of change in the Olympic National Park ecosystems. Results of the workshops provided the justification for selecting basic indicators of ecosystem integrity, effects of current threats to park resources, and focal resources of parks to detect both the currently evident and unforeseeable changes in park resources. We conclude Part I by exploring several generic statistical issues relevant to monitoring natural resources in Olympic National Park. Specifically we discuss trade-offs associated with sampling extensively versus sampling intensively in smaller geographic regions and describe a conceptual framework to guide development of a generic sampling frame for monitoring. We recommend partitioning Olympic National Park into three zones of decreasing accessibility to maximize monitoring efficiency. We present examples of how the generic sampling frame could be used to help ensure spatial integration of individual monitoring projects. Part II of the report is a record of the potential monitoring questions and indicators identified to date in our workshops. The presentation is organized according to the major system drivers, components, and processes identified in the intermediate-level working model of the Olympic National Park ecosystem. For each component of the park system, we develop the need and justification for monitoring, articulate park management issues, and describe key resources and ecosystem functions. We also present a pictorial conceptual model of each ecological subsystem, identify monitoring questions, and list potential indicators for each monitoring question. We conclude each section by identifying linkages of indicators to other ecological subsystems in our general ecosystem model, spatial and temporal contexts for monitoring (where and how often to monitor), and research and development needs. Part II represents the most current detailed listing of potential indicators—the material for subsequent discussions of monitoring priorities and selection of indicators for protocol development. Collectively, the sections of this report contain a comprehensive list of the important monitoring questions and potential indicators as well as recommendations for designing an integrated monitoring program. In Part I, Chapter 6 we provide recommendations on how to proceed with the important next steps in the design process: establishing priorities among the many possible monitoring questions and indicators, and beginning to research and design effective long-term monitoring protocols.

Information and Technology Report

Users and uses of Landsat 8 satellite imagery—2014 survey results

Executive Summary In 2013, Landsat 8 began adding high quality, global, moderate-resolution imagery to the more than 40-year archive of Landsat imagery. To assess the potential effects of the availability of Landsat 8 imagery on users and their work, the U.S. Geological Survey (USGS) Land Remote Sensing Program (LRS) initiated a survey of Landsat users. The objectives of the survey were to 1. Characterize various Landsat user groups, such as United States (U.S.) and international users and Landsat 8 and non-Landsat 8 users; 2. Identify any differences among user groups in uses and preferences; 3. Measure the importance of and satisfaction with Landsat 8 attributes; 4. Assess the importance to users of the frequency of usable imagery; and 5. Determine any challenges in using Landsat 8. The online survey was sent to 51,617 Landsat users registered with USGS in May 2014. Almost 13,000 people responded to the survey for a response rate of 25 percent (n = 12,966). Current Landsat users (users who had used Landsat in their work in the year prior to the survey) composed 89 percent of the sample (n = 11,549) and past Landsat users composed 11 percent (n = 1,417). The results reported here apply to current Landsat users registered with the USGS Earth Resources Observation and Science (EROS) Center. Users from 161 countries responded to the survey. Of those, 19 percent were citizens or permanent residents of the United States and 81 percent resided in other countries. More than 70 percent of current users had used Landsat 8 in the year prior to the survey. The majority of Landsat 8 users (65 percent) were established users who used Landsat imagery regularly both before and after Landsat 8 imagery became available. The average current Landsat user was male, 36 years old, and highly educated, with 9 years of experience using satellite imagery or geographic information system (GIS) software. Landsat 8 users had, on average, two more years of experience than non-Landsat 8 users. Users were employed predominantly by academic institutions (65 percent), followed by private businesses (13 percent), Federal governments (10 percent), State and local governments (6 percent), and nonprofit organizations (6 percent). Of the Landsat imagery obtained in the past year by current users, on average 31 percent came from a Landsat 8 sensor. An equivalent amount came from the Landsat 7 ETM+ sensor (33 percent); slightly less came from Landsats 4 and 5 TM sensors (27 percent). Much less came from Landsats 1 through 5 MSS sensors (5 percent). Overall, more than a third of users&rsquo; work used Landsat imagery (38 percent). Of this work, on average, 37 percent of the work was operational. Landsat 8 users considered a greater proportion of their work operational than non-Landsat 8 users (39 percent compared with 29 percent). Environmental sciences and management were the most commonly selected primary applications (selected by 42 percent of users). Land use/land cover (23 percent) was the second most commonly selected primary application, followed by education (12 percent), agriculture (9 percent), and planning and development (6 percent). Landsat 8 users were asked to rank the importance of certain attributes in determining whether to use Landsat 8 imagery in their work. The archive was ranked most important, followed by cost, spatial resolution, extent of coverage, data quality, and frequency of revisit. Users were asked how satisfied they were with these same attributes as they currently apply to Landsat 8 imagery. On average, users were most satisfied with lack of cost, extent of coverage, data quality, and the archive, but they were satisfied with all attributes. Users were asked how often they needed Landsat imagery to meet various requirements for their primary application. The survey question specifically asked how often users needed usable imagery, which differs from how often they would like the Landsat satellites to acquire an image. Users were asked to identify their needed frequency of usable imagery for the following levels: 1. Threshold level&mdash;the minimum frequency of usable imagery needed to be of any value to their primary application. 2. Breakthrough level&mdash;the frequency of usable imagery that would result in a significant improvement for their primary application of the imagery. 3. Target level&mdash;the frequency of usable imagery that would only provide a limited additional increase in the expected performance for their primary application. To meet the threshold level, three-quarters of users needed usable imagery every 17 days or less frequently. At the breakthrough level, two-thirds of users (64 percent) needed a usable image every 5&ndash;16 days. The current constellation of two satellites (Landsat 7 and 8) is capable of meeting the threshold and breakthrough needs of most users at least some of the time, but a single satellite would be highly unlikely to do so. Two-fifths of users (40 percent) felt that usable imagery provided every 4 days or more frequently would meet their target level which the current Landsat constellation cannot provide. Landsat 8 users were significantly more likely than non-Landsat 8 users to need usable imagery more frequently to meet their target levels. Additionally, U.S. Landsat 8 users were significantly more likely than other Landsat users to need usable imagery more frequently in order meet both their breakthrough and target levels. To explore the effect of the availability of Landsat 8 imagery on Landsat imagery use in general, established users (those who had consistently used Landsat imagery both before and after Landsat 8 imagery became available) using Landsat 8 imagery were asked about changes in the amount of Landsat imagery they used. The majority of established users using Landsat 8 imagery (60 percent) reported an average increase of 51 percent in the number of scenes obtained after Landsat 8 imagery became available. Landsat 8 users were asked if they had encountered challenges in using Landsat 8 whereas non-Landsat 8 users were asked if such challenges had played a role in why they were not using Landsat 8 imagery. Although many users did not encounter challenges when using or trying to use Landsat 8 data, slightly less than 30 percent did encounter issues with processing the data to a usable point. The most common issue reported was not being able to create or have access to a surface reflectance corrected product. Other challenges were related to the file sizes of images being too large to download, store, or analyze. There were no statistically significant differences between Landsat 8 and non-Landsat 8 users in terms of challenges encountered when using or trying to use the imagery, which indicates that users were not unduly discouraged by the challenges they may have encountered. When asked about potential consequences of not using Landsat 8, more than half of the non-Landsat 8 users did not report detrimental effects on their work from not using the imagery. Of those who did report detrimental effects, decreased quality of work, decreased scope of work, and increased time spent on work were the most common.

Open-File Report

Handbook of acute toxicity of chemicals to fish and aquatic invertebrates : summaries of toxicity tests conducted at Columbia National Fisheries Research Laboratory, 1965-78

Acute toxicity is a major subject of research at Columbia National Fisheries Research Laboratory for evaluating the impact of toxic chemicals on fishery resources. The Laboratory has played a leading role in developing research technology for toxicity testing and data interpretation. In 1965-78, more than 400 chemicals were tested against a variety of invertebrates and fish species representative of both cold- and warm-water climates. The use of acute toxicity tests for assessing the potential hazard of chemical contaminants to aquatic organisms is well documented (Boyd 1957; Henderson et al. 1960; Sanders and Cope 1966; Macek and McAllister 1970). Static acute toxicity tests provide rapid and (within limits) reproducible concentration-response curves for estimating toxic effects of chemicals on aquatic organisms. These tests provide a database for determining relative toxicity of a large number of chemicals to a variety of species and for estimating acute effects of chemical spills on natural aquatic systems; they also assist in determining priority and design of additional toxicity studies. Acute toxicity tests usually provide estimates of the exposure concentration causing 50% mortality (LC50) to test organisms during a specified period of time. For certain invertebrates, the effective concentration is based on immobilization, or some other identifiable endpoint, rather than on lethality. The application of the LC50 has gained acceptance among toxicologists and is generally the most highly rated test for assessing potential adverse effects of chemical contaminants to aquatic life (Brungs and Mount 1978; American Institute for Biological Sciences 1978 a ). The literature contains numerous papers dealing with the acute toxicity of chemicals to freshwater organisms. However, there is a tremendous need for a concise compendium of toxicity data covering a large variety of chemicals and test species. This Handbook is a compilation of a large volume of acute toxicity data from the Columbia Laboratory and its field laboratories. It presents definitive acute toxicity data on 271 chemicals tested against a variety of freshwater invertebrates and fishes. The chemicals represent all major groups of pesticides, as well as numerous industrial chemicals. This compilation should serve as a useful database for the many agencies and organizations dealing with research and management programs concerned with the impact of chemicals on aquatic resources. The Columbia Laboratory has played a major role in developing currently used standard methodology for static acute toxicity testing. The use of standardized methodology greatly reduces variation in results. The data presented here have been carefully scrutinized to eliminate tests that failed to follow acceptable procedures. Handling of test organisms and procedures for static toxicity tests followed those described by Lennon and Walker (1964) and Macek and McAllister (1970), and conform well with those recommended by Brauhn and Schoettger (1975) and the Committee on Methods for Toxicity Tests with Aquatic Organisms (1975). The species of fish and invertebrates that were tested are listed in phylogenetic order in Tables 1 and 2. Fish were obtained from Federal and State hatcheries as either eggs or fry. Original stocks of invertebrates were collected and cultured from wild populations with no known source of contamination; these populations were replenished regularly. The invertebrates were cultured in the Laboratory by methods similar to those described by Sanders and Cope (1966). Test chemicals usually consisted of technical or analytical grade samples of known purity. Formulations of the chemicals were also tested when available. When purity of test chemicals was known, all calculated concentrations were based on percent active ingredients. Stock solutions were prepared immediately before each test, with commercial grade acetone as the carrier solvent. Occasionally, ethanol or dimethyl-formamide was substituted. Solvent concentrations did not exceed 0.5 mL/L in final dilution water. Test water (dilution water) was reconstituted from deionized water of at least 10 6 ohms resistivity by the addition of appropriate reagent grade chemicals (Marking 1969). Water was buffered to maintain a pH of 7.2 to 7.5, an alkalinity of 30 to 35 mg/L, and a hardness of 40 to 50 mg/L as CaCO 3 . Test water was mixed thoroughly and aerated before transfer into test chambers. Fish were acclimated to dilution water by gradually changing the water in acclimated tanks from 100% well water to 100% reconstituted water over a 1- to 3-day period at the desired testing temperature. Invertebrates were acclimated from well water to dilution water over a 4- to 6-h period. Toxicity tests were conducted under static conditions without aeration, and the organisms were not fed during acclimation or testing. Temperature of test solutions was maintained within ± 1°C of that required for a given test. Toxicity tests with fish were conducted in 18.9-liter (5-gal) wide-mouthed jars containing 15 liters of test solution. Fingerling fish weighing 0.2 to 1.5 g were tested at each concentration. Caution was taken not to exceed 0.8 g of test organisms per liter of solution. Duplicate test chambers were used to accommodate larger fish. Test chambers varied in size for invertebrates, depending on the species used; volume of test solution ranged from 0.25 to 4 liters. At least 10 organisms were exposed to each concentration for all definitive tests. At least six concentrations were used per toxicity test. The tests began upon initial exposure to the toxicant and continued for 96 h. Immobilization tests with invertebrates were conducted for only 48 h. The number of dead or affected organisms in each test chamber were recorded and the dead organisms were removed every 24 h; general observations on the condition of test organisms were also recorded at these times. Toxicity data were analyzed by a statistical method described by Litchfield and Wilcoxon (1949) to determine LC50 (theoretical estimate of the concentration lethal to 50% of the test animals) and 95% confidence intervals. This method is recommended by the American Public Health Association (1971) and by Sprague (1969) for determining median lethal concentrations. The procedure is easily modified for computing a single LC50 when replicate tests are performed.

Resource Publication