USGS Science⌕ Search

SEARCH · USGS Science

Results for “Developments in Water 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 1,297 records · Page 72Linked to original sources

An evaluation of biomarkers of reproductive function and potential contaminant effects in Florida largemouth bass (Micropterus salmoides floridanus) sampled from the St. Johns River

The objective of this study was to describe and compare several reproductive parameters for Florida largemouth bass ( Micropterus salmoides floridanus ) inhabiting the St. Johns River and exposed to different types and/or degrees of contamination. Welaka was selected as the reference site in this study because of its low urban and agricultural development, Palatka is in close proximity to a paper mill plant, the Green Cove site is influenced by marine shipping activities and Julington Creek site receives discharges of domestic wastewater and storm water runoff from recreational boating marinas. For this study, bass were sampled both prior to (September 1996) and during the spawning season (February 1997). In order to characterize chemical exposure, bass livers were analyzed for up to 90 trace organics and 11 trace metal contaminants. Reproductive parameters measured included gonadosomatic index (GSI), histological evaluation of gonads and plasma concentrations of vitellogenin (VTG), 17β-estradiol (E 2 ) and 11-ketotestosterone (11-KT). In general, the sum of organic chemicals was highest in livers from Palatka bass and bass from Green Cove and Julington Creek had higher hepatic concentrations of low molecular polycyclic aromatic hydrocarbons and polychlorinated biphenyls when compared to fish from Welaka. Metals were more variable across sites, with highest mean concentrations found in bass from either Julington Creek (Ag, As, Cr, Cu, Zn) or Welaka (Cd, Hg, Pb, Se, Tn). Female bass from Palatka and Green Cove had lower concentrations of E 2 , VTG and lower GSI in relation to Welaka. Males from Palatka and Green Cove showed comparable declines in 11-KT in relation to males from Julington Creek and GSI were decreased only in Palatka males. These results indicate a geographical trend in reproductive effects, with changes being most pronounced at the site closest to the paper mill (Palatka) and decreasing as the St. Johns River flows downstream. Since reproductive alterations were most evident in bass sampled from the site closest to the paper mill discharge, it is possible that exposure to these effluents might explain at least some of the results reported here. However, the presence of reproductive alterations in fish sampled at a considerable distance from the mill discharge (Green Cove, 40 km) would suggest exposure to chemicals released from sources other than the paper mill plant. It is clear that additional studies are needed to evaluate the potential impact of these reproductive changes in populations of Florida largemouth bass inhabiting the St. Johns River.

Florida↗

The Boreal-Arctic Wetland and Lake Dataset (BAWLD)

Methane emissions from boreal and arctic wetlands, lakes, and rivers are expected to increase in response to warming and associated permafrost thaw. However, the lack of appropriate land cover datasets for scaling field-measured methane emissions to circumpolar scales has contributed to a large uncertainty for our understanding of present-day and future methane emissions. Here we present the Boreal–Arctic Wetland and Lake Dataset (BAWLD), a land cover dataset based on an expert assessment, extrapolated using random forest modelling from available spatial datasets of climate, topography, soils, permafrost conditions, vegetation, wetlands, and surface water extents and dynamics. In BAWLD, we estimate the fractional coverage of five wetland, seven lake, and three river classes within 0.5 × 0.5∘ grid cells that cover the northern boreal and tundra biomes (17 % of the global land surface). Land cover classes were defined using criteria that ensured distinct methane emissions among classes, as indicated by a co-developed comprehensive dataset of methane flux observations. In BAWLD, wetlands occupied 3.2 × 106 km2 (14 % of domain) with a 95 % confidence interval between 2.8 and 3.8 × 106 km2. Bog, fen, and permafrost bog were the most abundant wetland classes, covering ∼ 28 % each of the total wetland area, while the highest-methane-emitting marsh and tundra wetland classes occupied 5 % and 12 %, respectively. Lakes, defined to include all lentic open-water ecosystems regardless of size, covered 1.4 × 106 km2 (6 % of domain). Low-methane-emitting large lakes (>10 km2) and glacial lakes jointly represented 78 % of the total lake area, while high-emitting peatland and yedoma lakes covered 18 % and 4 %, respectively. Small (<0.1 km2) glacial, peatland, and yedoma lakes combined covered 17 % of the total lake area but contributed disproportionally to the overall spatial uncertainty in lake area with a 95 % confidence interval between 0.15 and 0.38 × 106 km2. Rivers and streams were estimated to cover 0.12 × 106 km2 (0.5 % of domain), of which 8 % was associated with high-methane-emitting headwaters that drain organic-rich landscapes. Distinct combinations of spatially co-occurring wetland and lake classes were identified across the BAWLD domain, allowing for the mapping of “wetscapes” that have characteristic methane emission magnitudes and sensitivities to climate change at regional scales. With BAWLD, we provide a dataset which avoids double-accounting of wetland, lake, and river extents and which includes confidence intervals for each land cover class. As such, BAWLD will be suitable for many hydrological and biogeochemical modelling and upscaling efforts for the northern boreal and arctic region, in particular those aimed at improving assessments of current and future methane emissions. Data are freely available at https://doi.org/10.18739/A2C824F9X (Olefeldt et al., 2021).

Earth System Science Data↗

Understanding the impacts of surface-groundwater conditions on stream fishes under altered baseflow conditions

Persistence of aquatic fauna depends on the conditions and connectivity of surface water and groundwater. In light of altered baseflows and both current and future predicted increases in stream temperatures, it is important to assess current thermal conditions, examine thermal responses of aquatic fauna, and evaluate water-management practices. Our study objectives were to determine (1) how changes in baseflow levels in the Kiamichi River influence hyporheic exchange, which correspondingly influences temperature at the reach scale; (2) temperature tolerances of stream fishes as a means for predicting how habitat complexity influences stream-fish populations; and (3) assess how dam releases influence the downstream temperature and dissolved oxygen regime during the low-flow period. We quantified hyporheic exchange at four reaches and, as expected, found higher groundwater exchange via transient storage occurred at the upstream sites. The net groundwater flux estimation was negative for the majority of reaches indicating that surface water is lost to groundwater during summer (i.e., losing), baseflow conditions. We determined critical thermal maximum (CTMax) for 17 stream fishes and thermal tolerances ranged 32-38°C. We determined the average thermal tolerance for two habitat fish guilds to calculate changes in thermal stress due to hypothetical reservoir release scenarios. We developed a process-based Water Quality Analysis Simulation Program model to predict downstream temperature conditions over 74-km of river in response to reservoir releases that corresponded to discharges of 0.00 (control), 0.34, 0.59, 0.76, 1.13, and 1.50 m3/s. Based on the dissolved oxygen conditions observed in 2015 and 2017 and biological oxygen demand sampling results, reservoir releases did not directly reduce dissolved oxygen concentrations in the Kiamichi River (though dissolved oxygen concentrations are limited to current water-release strategies by the managing agency). We simulated three scenarios using three water-release temperatures: 27.64°C, 26.00°C and 24.07°C that corresponded to average reservoir temperatures at gate locations on the dam. We compared the predicted temperature time series with CTMax of two fish-habitat guilds to quantify the cumulative time when stream fishes experienced severe thermal stress downstream from Sardis Reservoir. According to our simulations, reservoir releases would be capable of regulating downstream water temperature during the summer baseflow period. The 0.00 m3/s scenario resulted in 130 h of thermal stress for benthic fishes, and 73 h for mid-column fishes. As expected, thermal relief increased with increasing release magnitude and decreasing release water temperature. The 0.34 m3/s release scenario reduced thermal stress (range is simulations from the top and bottom gate) by 11-18% for mid-column fishes and 8-12% for benthic fishes with an effective distance (where the cumulative time above CTMax was reduced by half) of 1-2 km for both guilds. The 0.59 m3/s release scenario reduced thermal stress by 18-25% for mid-column fishes and 12-20% for benthic fishes with effective distances of 4-8 km and 2-7 km, respectively. Three releases representing pre-dam flow magnitudes (0.76, 1.13 and 1.50 m3/s released from top gate) reduced thermal stress up to 46% for mid-column fishes and 41% for benthic fishes with an effective distance of 13-16 km, respectively. Lastly, we quantified temperature-induced stress via whole-body cortisol concentration of six stream fishes in response to prolonged thermal exposure at two temperatures (27°C and 32°C). We found no difference in cortisol levels between temperatures for any of the six species, indicating acclimation to elevated temperatures during the test period. However, Highland Stoneroller Campostoma spadiceum expressed cortisol concentrations greater than typical basal levels at both temperatures, suggesting stress from factors other than temperature (i.e., captivity). Our results suggest different reservoir-release options could improve downstream thermal-fish habitat during the summer baseflow period.

Cooperator Science Series↗

Occurrence of cyanazine compounds in groundwater: Degradates more prevalent than the parent compound

A recently developed analytical method using liquid chromatography/mass spectrometry was used to investigate the occurrence of cyanazine and its degradates cyanazine acid (CAC), cyanazine amide (CAM), deethylcyanazine (DEC), and deethylcyanazine acid (DCAC) in groundwater. This research represents some of the earliest data on the occurrence of cyanazine degradates in groundwater. Although cyanazine was infrequently detected in the 64 wells across Iowa sampled in 1999, cyanazine degradates were commonly found during this study. The most frequently detected cyanazine compound was DCAC (32.8%) followed by CAC (29.7%), CAM (17.2%), DEC (3.1%), and cyanazine (3.1%). The frequency of detection for cyanazine or one or more of its degradates (CYTOT) was more than 12-fold over that of cyanazine alone (39.1% for CYTOT versus 3.1% for cyanazine). Of the total measured concentration of cyanazine, only 0.2% was derived from its parent compound - with DCAC (74.1%) and CAC (18.4%) comprising 92.5% of this total. Thus, although DCAC and CAC had similar frequencies of detection, DCAC was generally present in higher concentrations. No concentrations of cyanazine compounds for this study exceeded water-quality criteria for the protection of human health. Only cyanazine, however, has such a criteria established. Nevertheless, because these cyanazine degradates are still chlorinated, they may have similar toxicity as their parent compound - similar to what has been found with the chlorinated degradates of atrazine. Thus, the results of this study documented that data on the degradates for cyanazine are critical for understanding its fate and transport in the hydrologic system. Furthermore, the prevalence of the chlorinated degradates of cyanazine found in groundwater suggests that to accurately determine the overall effect on human health and the environment from cyanazine its degradates should also be considered. In addition, because CYTOT was found in 57.6% of the samples collected from alluvial aquifers, about 2-5 times more frequently than the other major aquifer types (glacial drift, bedrock/karst, bedrock/nonkarst) under investigation, this finding has long-term implications for the occurrence of CYTOT in streams. It is anticipated that low-level concentrations of CYTOT will continue to be detected in streams for years after the use of cyanazine has terminated (scheduled for the year 2000 in the United States), primarily through its movement from groundwater into streams during base-flow conditions.

Iowa↗

Application of artificial neural networks to assess pesticide contamination in shallow groundwater

In this study, a feed-forward back-propagation neural network (BPNN) was developed and applied to predict pesticide concentrations in groundwater monitoring wells. Pesticide concentration data are challenging to analyze because they tend to be highly censored. Input data to the neural network included the categorical indices of depth to aquifer material, pesticide leaching class, aquifer sensitivity to pesticide contamination, time (month) of sample collection, well depth, depth to water from land surface, and additional travel distance in the saturated zone (i.e., distance from land surface to midpoint of well screen). The output of the neural network was the total pesticide concentration detected in the well. The model prediction results produced good agreements with observed data in terms of correlation coefficient (R = 0.87) and pesticide detection efficiency (E = 89%), as well as good match between the observed and predicted "class" groups. The relative importance of input parameters to pesticide occurrence in groundwater was examined in terms of R, E, mean error (ME), root mean square error (RMSE), and pesticide occurrence "class" groups by eliminating some key input parameters to the model. Well depth and time of sample collection were the most sensitive input parameters for predicting the pesticide contamination potential of a well. This infers that wells tapping shallow aquifers are more vulnerable to pesticide contamination than those wells tapping deeper aquifers. Pesticide occurrences during post-application months (June through October) were found to be 2.5 to 3 times higher than pesticide occurrences during other months (November through April). The BPNN was used to rank the input parameters with highest potential to contaminate groundwater, including two original and five ancillary parameters. The two original parameters are depth to aquifer material and pesticide leaching class. When these two parameters were the only input parameters for the BPNN, they were not able to predict contamination potential. However, when they were used with other parameters, the predictive performance efficiency of the BPNN in terms of R, E, ME, RMSE, and pesticide occurrence "class" groups increased. Ancillary data include data collected during the study such as well depth and time of sample collection. The BPNN indicated that the ancillary data had more predictive power than the original data. The BPNN results will help researchers identify parameters to improve maps of aquifer sensitivity to pesticide contamination. ?? 2006 Elsevier B.V. All rights reserved.

Science of the Total Environment↗

Conservation planning for the Colorado River in Utah

Strategic planning is increasingly recognized as necessary for providing the greatest possible conservation benefits for restoration efforts. Rigorous, science-based resource assessment, combined with acknowledgement of broader basin trends, provides a solid foundation for determining effective projects. It is equally important that methods used to prioritize conservation investments are simple and practical enough that they can be implemented in a timely manner and by a variety of resource managers. With the help of local and regional natural resource professionals, we have developed a broad-scale, spatially-explicit assessment of 146 miles (~20,000 acres) of the Colorado River mainstem in Grand and San Juan Counties, Utah that will function as the basis for a systematic, practical approach to conservation planning and riparian restoration prioritization. For the assessment we have: 1) acquired, modified or created spatial datasets of Colorado River bottomland conditions; 2) synthesized those datasets into habitat suitability models and estimates of natural recovery potential, fire risk and relative cost; 3) investigated and described dominant ecosystem trends and human uses, and; 4) suggested site selection and prioritization approaches. Partner organizations (The Nature Conservancy, National Park Service, Bureau of Land Management and Utah Forestry Fire and State Lands) are using the assessment and datasets to identify and prioritize a suite of restoration actions to increase ecosystem resilience and improve habitat for bottomland species. Primary datasets include maps of bottomland cover types, bottomland extent, maps of areas inundated during high and low flow events, as well as locations of campgrounds, roads, fires, invasive vegetation treatment areas and other features. Assessment of conditions and trends in the project area entailed: 1) assemblage of existing data on geology, changes in stream flow, and predictions of future conditions; 2) identification of fish and wildlife species present and grouping species into Conservation Elements (CEs) based on habitat needs, and: 3) acquisition, review and creation of spatial datasets characterizing vegetation, fluvial geomorphic and human features within the bottomland. Interpretation of aerial imagery and assimilation of pre-existing spatial data were central to our efforts in characterizing resources conditions. Detailed maps of vegetation and channel habitat features in the project area were generated from true color, high resolution (0.3m) imagery flown September 16, 2010. We also mapped channel habitat features at high flow on 1.0-m resolution, publicly available, true color imagery. We obtained additional layers such as land ownership, roads, fire history, non-native vegetation treatment areas, and recreational use features from public sources and project partners. Habitat suitability models were created for groups of terrestrial species by combining spatial datasets with the habitat needs of conservation elements, guided by literature, where available, and extensive use of expert knowledge. Conservation elements for endangered fish species life stages were identified but not modeled. Terrestrial CE’s included: • Riparian Overstory -yellow-billed cuckoo, Bullock’s oriole, black-headed grosbeak, blue grosbeak, warbling vireo, Cooper's hawk, screech owl, saw-whet owl, and bald eagle, (best: tall trees, dense canopy, diverse shrub understory, no tamarisk); • Riparian Understory - southwestern willow flycatcher, common yellowthroat, yellow warbler, yellow-breasted chat, beaver, northern river otter, black-necked garter snake, (best: dense mesic shrubs near still water, no tamarisk); • Bat Feeding - Allen's big-eared bat, Townsend's big-eared bat, fringed myotis, Yuma myotis, big free-tailed bat, spotted bat (best: diverse vegetation, close to still water); • Bat Watering - big free-tailed and spotted bats (best: still water with no tall vegetation); •

Utah↗

Using environmental DNA (eDNA) to assess the presence of cavefish and cave crayfish populations in caves of the Ozark Highlands

Many cavefishes and cave crayfishes are considered of conservation concern; however, sampling these species is inherently difficult given their occupied environments. The goal of our project was to verify the presence of select karst organisms while developing the foundation for sampling approaches that might be useful to conservation and management agencies. Our project objectives were to develop assays to amplify deoxyribonucleic acid (DNA) from several species of Ozark cavefishes and cave crayfishes and complete an initial surveillance of locations across the Ozark Highlands using environmental DNA (eDNA). Using DNA either provided by agency cooperators or that we extracted from tissue samples, we PCR amplified and then sequenced the Cytochrome Oxidase 1 (CO1) gene for cave crayfishes and the NADH Dehydrogenase Subunit 2 (ND2) gene for cavefishes. We developed species-specific primers and probes for five cave crayfishes and two cavefishes. From February 2017 to May 2017, we sampled 1–5 sampling units from 42 caves, wells, and springs (i.e., sites) using eDNA and traditional visual surveys. We measured physicochemical parameters at each sampling unit to estimate detection probability associated with both techniques. We also calculated two occupancy covariates for each site using geospatial data. We successfully amplified Troglichthys rosae DNA from the environment and detected DNA representing this species at 24 of 40 sites. At 16 of the sites where we detected T. rosae DNA, we did not visually observe the species. Although our assay for Typlichthys eigenmanni successfully amplified the target DNA from the environment, it also resulted in false absences where the species was visually confirmed. Using eDNA to detect cave crayfishes was much more difficult. The assay for Cambarus subterraneus did not work for eDNA samples and we were unable to pick up DNA from the environment, even at locations where it was visually confirmed. Alternatively, the eDNA surveys worked well for C. tartarus and we were able to amplify DNA at every site where it was visually observed. Our assay for C. aculabrum was based on a single sample obtained from GenBank, and did not amplify eDNA from field samples. Lastly, our eDNA results from samples in the known range of Orconectes stygocaneyi suggested the species may be found at an additional cave. Detection using eDNA based on our O. stygocaneyi assay was likely low because it was designed from a pseudogene; however, positive eDNA samples were sequenced to confirm species-specific DNA. Detection probability of both cavefishes and cave crayfishes varied by survey technique and was influenced by water volume, water clarity, water velocity, and substrate. Detection of cavefishes and cave crayfishes via visual surveys decreased when water volume increased, whereas detection using eDNA increased with greater water volume. Detection between taxa using either sample method was highest in habitats classified by fine substrates, except for eDNA detection of crayfishes which was greatest in coarse substrates. Detection of cavefishes increased with water clarity, but detection of cave crayfishes increased with turbidity. Detection probability of both cavefishes and crayfishes using eDNA increased slightly with water velocity, but decreased with visual surveys as water velocity increased. Occupancy by both taxa was positively related to particular geologic series. Crayfish occupancy was negatively related to fine-scale anthropogenic disturbance (i.e., 500-m buffer around the site), whereas crayfish showed no relationship with disturbance. Our results suggest possible range extensions, provide insights to factors driving detection using both sample techniques, and suggest areas where recharge zones may be shared among caves. Future efforts focused on a comprehensive evaluation of genetic diversity among cave crayfishes to improve assay design could improve detection and the applicability of eDNA as a supplemental and non-invasive sampling approach.

Cooperator Science Series↗

Integrating recent land cover mapping efforts to update the National Gap Analysis Program's species habitat map

Over the past decade, great progress has been made to develop national extent land cover mapping products to address natural resource issues. One of the core products of the GAP Program is range-wide species distribution models for nearly 2000 terrestrial vertebrate species in the U.S. We rely on deductive modeling of habitat affinities using these products to create models of habitat availability. That approach requires that we have a thematically rich and ecologically meaningful map legend to support the modeling effort. In this work, we tested the integration of the Multi-Resolution Landscape Characterization Consortium's National Land Cover Database 2011 and LANDFIRE's Disturbance Products to update the 2001 National GAP Vegetation Dataset to reflect 2011 conditions. The revised product can then be used to update the species models. We tested the update approach in three geographic areas (Northeast, Southeast, and Interior Northwest). We used the NLCD product to identify areas where the cover type mapped in 2011 was different from what was in the 2001 land cover map. We used Google Earth and ArcGIS base maps as reference imagery in order to label areas identified as "changed" to the appropriate class from our map legend. Areas mapped as urban or water in the 2011 NLCD map that were mapped differently in the 2001 GAP map were accepted without further validation and recoded to the corresponding GAP class. We used LANDFIRE's Disturbance products to identify changes that are the result of recent disturbance and to inform the reassignment of areas to their updated thematic label. We ran species habitat models for three species including Lewis's Woodpecker ( Melanerpes lewis ) and the White-tailed Jack Rabbit ( Lepus townsendii ) and Brown Headed nuthatch ( Sitta pusilla ). For each of three vertebrate species we found important differences in the amount and location of suitable habitat between the 2001 and 2011 habitat maps. Specifically, Brown headed nuthatch habitat in 2011 was −14% of the 2001 modeled habitat, whereas Lewis's Woodpecker increased by 4%. The white-tailed jack rabbit ( Lepus townsendii ) had a net change of −1% (11% decline, 10% gain). For that species we found the updates related to opening of forest due to burning and regenerating shrubs following harvest to be the locally important main transitions. In the Southeast updates related to timber management and urbanization are locally important.

The International Archives of the Photogrammetry, ↗

Development of a Kemp’s ridley sea turtle stock assessment model

We developed a Kemp’s ridley ( Lepidochelys kempii ) stock assessment model to evaluate the relative contributions of conservation efforts and other factors toward this critically endangered species’ recovery. The Kemp’s ridley demographic model developed by the Turtle Expert Working Group (TEWG) in 1998 and 2000 and updated for the binational recovery plan in 2011 was modified for use as our base model. The TEWG model uses indices of the annual reproductive population (number of nests) and hatchling recruitment to predict future annual numbers of nests on the basis of a series of assumptions regarding age and maturity, remigration interval, sex ratios, nests per female, juvenile mortality, and a putative ‘‘turtle excluder device effect’’ multiplier starting in 1990. This multiplier was necessary to fit the number of nests observed in 1990 and later. We added the effects of shrimping effort directly, modified by habitat weightings, as a proxy for all sources of anthropogenic mortality. Additional data included in our model were incremental growth of Kemp’s ridleys marked and recaptured in the Gulf of Mexico, and the length frequency of stranded Kemp’s ridleys. We also added a 2010 mortality factor that was necessary to fit the number of nests for 2010 and later (2011 and 2012). Last, we used an empirical basis for estimating natural mortality, on the basis of a Lorenzen mortality curve and growth estimates. Although our model generated reasonable estimates of annual total turtle deaths attributable to shrimp trawling, as well as additional deaths due to undetermined anthropogenic causes in 2010, we were unable to provide a clear explanation for the observed increase in the number of stranded Kemp’s ridleys in recent years, and subsequent disruption of the species’ exponential growth since the 2009 nesting season. Our consensus is that expanded data collection at the nesting beaches is needed and of high priority, and that 2015 be targeted for the next stock assessment to evaluate the 2010 event using more recent nesting and in-water data.

Gulf of Mexico Science↗

Vulnerability to climate change of managed stocks in the California Current large marine ecosystem

Introduction: Understanding how abundance, productivity and distribution of individual species may respond to climate change is a critical first step towards anticipating alterations in marine ecosystem structure and function, as well as developing strategies to adapt to the full range of potential changes. Methods: This study applies the NOAA (National Oceanic and Atmospheric Administration) Fisheries Climate Vulnerability Assessment method to 64 federally-managed species in the California Current Large Marine Ecosystem to assess their vulnerability to climate change, where vulnerability is a function of a species’ exposure to environmental change and its biological sensitivity to a set of environmental conditions, which includes components of its resiliency and adaptive capacity to respond to these new conditions. Results: Overall, two-thirds of the species were judged to have Moderate or greater vulnerability to climate change, and only one species was anticipated to have a positive response. Species classified as Highly or Very Highly vulnerable share one or more characteristics including: 1) having complex life histories that utilize a wide range of freshwater and marine habitats; 2) having habitat specialization, particularly for areas that are likely to experience increased hypoxia; 3) having long lifespans and low population growth rates; and/or 4) being of high commercial value combined with impacts from non-climate stressors such as anthropogenic habitat degradation. Species with Low or Moderate vulnerability are either habitat generalists, occupy deep-water habitats or are highly mobile and likely to shift their ranges. Discussion: As climate-related changes intensify, this work provides key information for both scientists and managers as they address the long-term sustainability of fisheries in the region. This information can inform near-term advice for prioritizing species-level data collection and research on climate impacts, help managers to determine when and where a precautionary approach might be warranted, in harvest or other management decisions, and help identify habitats or life history stages that might be especially effective to protect or restore.

California, Oregon, Washington↗

Seismic and tsunami hazard in Puerto Rico and the Virgin Islands

Executive Summary Puerto Rico and the Virgin Islands are located at an active plate boundary between the North American plate and the northeast corner of the Caribbean plate. The region was subject in historical times to large magnitude earthquakes and devastating tsunamis. A major downward tilt of the sea floor north of Puerto Rico and the Virgin Islands, large submarine rockslides, and an unusually large negative gravity anomaly are also indicative of a tectonically active region. Scientists have so far failed to explain the deformation of this region in a coherent and predictable picture, such as in California, and this has hampered their ability to assess seismic and tsunami hazards in the region. The NE corner of the Caribbean is unique among the seismically-active regions of the United States in that it is mostly covered by water. This fact presents an additional challenge for seismic and tsunami hazard assessment and mitigation. The workshop, convened in San Juan on March 23-24, 1999, was "historic" in that it brought together for the first time a broad spectrum of scientists, engineers, and public and private sector officials who deal with such diverse questions as tectonic models, probabilistic assessment of seismic hazard, prediction of tsunami runup, strong ground motion, building codes, stability of man-made structures, and the public’s preparedness for natural disasters. It was an opportunity for all the participants to find out how their own activity fit into the broad picture of science and how it aids society in hazard assessment and mitigation. In addition, the workshop was offered as a continuing education course at the Colegio de Ingenieros y Agrimensores de Puerto Rico, which assured a rapid dissemination of the results to the local community. A news conference which took place during the workshop alerted the public to the efforts of the USGS, other Federal agencies, the Commonwealth of Puerto Rico, universities and the private sector. During the first day of the workshop, participants from universities, federal institutions, and consulting firms in Puerto Rico, the Virgin Islands, the continental U.S., Dominican Republic, and Europe reviewed the present state of knowledge including a review and discussion of present plate models, recent GPS and seismic reflection data, seismicity, paleoseismology, and tsunamis. The state of earthquake/tsunami studies in Puerto Rico was presented by several faculty members from the University of Puerto Rico at Mayaguez. A preliminary seismic hazard map was presented by the USGS and previous hazard maps and economic loss assessments were considered. During the second day, the participants divided into working groups and prepared specific recommendations for future activities in the region along the six following topics below. Highlights of these recommended activities are: Marine geology and geophysics – Acquire deep-penetration seismic reflection and refraction data, deploy temporary ocean bottom seismometer arrays to record earthquakes, collect high-resolution multibeam bathymetry and side scan sonar data of the region, and in particular, the near shore region, and conduct focussed high-resolution seismic studies around faults. Determine slip rates of specific offshore faults. Assemble a GIS database for available marine geological and geophysical data. Paleoseismology and active faults - Field reconnaissance aimed at identifying Quaternary faults and determining their paleoseismic chronology and slip rates, as well as identifying and dating paleoliquefaction features from large earthquakes. Quaternary mapping of marine terraces, fluvial terraces and basins, beach ridges, etc., to establish framework for understanding neotectonic deformation of the island. Interpretation of aerial photography to identify possible Quaternary faults. Earthquake seismology – Determine an empirical seismic attenuation function using observations from local seismic networks and recently-installed broad-band stations. Evaluate existing earthquake catalogs from local networks and regional stations, complete the catalogs. Transcribe the pre-1991 network data from 9-track tape onto more stable archival media. Calibrate instruments of local networks. Use GPS measurement to constrain deformation rates used in seismic-hazard maps. Engineering – Prepare liquefaction susceptibility maps for the urban areas. Update and improve databases for types of site conditions. Collect site effect observations and near-surface geophysical measurements for future local (urban-area) hazard maps. Expand the number of instruments in the strong motion program. Develop fragility curves for Puerto Rico construction types and details, and carry out laboratory testing on selected types of mass-produced construction. Consider tsunami design in shoreline construction projects. Tsunami hazard - Extract tsunami observations from archives and develop a Caribbean historical tsunami database. Analyze prehistoric tsunami deposits. Collect accurate, up-to-date, near-shore topography and bathymetry for accurate inundation models. Prepare tsunami flooding and evacuation maps. Establish a Caribbean Tsunami Warning System for Puerto Rico and the Virgin Islands. Evaluate local, regional, national, and global seismic networks and equipment, and their role in a tsunami warning system. Societal concerns – Prepare warning messages, protocols, and evacuation routes for earthquake, tsunami, and landslide hazards for Puerto Rico and the U.S. Virgin Islands. Advocate enforcement of existing building codes. Prepare non-technical hazard assessment maps for political and educational uses. Raise the awareness of potentially affected populations by presentations at elementary schools, by the production of a tsunami video, and by distribution of earthquake preparedness manuals in newspaper supplements. Promote partnerships at state and federal level for long-term earthquake and tsunami hazard mitigation. This partnership should also include the private sector such as the insurance industry, telecommunication companies, and the engineering community. The following reports of the various working groups are the cumulative recommendations of the community of scientists, engineers, and public officials, who participated in the workshop. The list of participants and the workshop’s agenda are given in the appendix. Marine and Geology and Geophysics Working Group Paleoseismology and Active Faults Working Group Joint Working Group for Earthquake Seismology and Engineering Tsunami Working Group Societal Concerns Working Group

Puerto Rico and the Virgin Islands↗

Aeolian dust in Colorado Plateau soils: Nutrient inputs and recent change in source

Aeolian dust (windblown silt and clay) is an important component in arid-land ecosystems because it may contribute to soil formation and furnish essential nutrients. Few geologic surfaces, however, have been characterized with respect to dust-accumulation history and resultant nutrient enrichment. We have developed a combination of methods to identify the presence of aeolian dust in arid regions and to evaluate the roles of this dust in ecosystem processes. Unconsolidated sandy sediment on isolated surfaces in the Canyonlands region of the Colorado Plateau differs greatly in mineralogical and chemical composition from associated bedrock, mainly aeolian sandstone. Detrital magnetite in the surficial deposits produces moderately high values of magnetic susceptibility, but magnetite is absent in nearby bedrock. A component of the surficial deposits must be aeolian to account for the abundance of magnetite, which formed originally in far-distant igneous rocks. Particle-size analysis suggests that the aeolian dust component is typically as much as 20–30%. Dust inputs have enriched the sediments in many elements, including P, Mg, Na, K, and Mo, as well as Ca, at sites where bedrock lacks calcite cement. Soil-surface biologic crusts are effective dust traps that apparently record a change in dust sources over the past several decades. Some of the recently fallen dust may result from human disturbance of land surfaces that are far from the Canyonlands, such as the Mojave Desert. Some land-use practices in the study area have the potential to deplete soil fertility by means of wind-erosion removal of aeolian silt. Many studies have addressed the presence of aeolian dust in soils and surficial deposits in deserts to provide important geologic and ecologic information bearing on landscape dynamics (1–16). From this body of work, we have improved understanding about: ( i ) current and past sources and flux of dust, hence changing conditions of dust emission; ( ii ) the genesis of desert soils; ( iii ) the influences of aeolian silt and clay on water-infiltration rates in soil; ( iv ) the evolution of desert surfaces (such as desert pavement) relevant to surface stability, as well as the distribution of surface and subsurface water; and ( v ) interrelations among aeolian dust, distribution of plants and soil crust, rain-water runoff, and productivity. Nevertheless, we lack fundamental knowledge about the accumulation history of aeolian dust on most landscapes and about how to discriminate between contributions from parent material and aeolian dust to the biotic system. Although many different methods provide clear evidence for aeolian input into soils (4, 6, 8, 9, 11, 13, 17, 18), ecosystem studies would benefit from rapid assessment of aeolian components in soils. Here we demonstrate that magnetic methods, which characterize the type and distribution of iron oxide minerals, can be applied to detect the presence of aeolian dust in young surficial sediments and soils over large arid-land areas. The magnetic results, combined with geochemical and textural analyses, form a basis for understanding the influence of fine-grained aeolian inputs on soil fertility of the central Colorado Plateau, Utah.

Utah↗

Lateglacial and Holocene climate, disturbance and permafrost peatland dynamics on the Seward Peninsula, western Alaska

Northern peatlands have accumulated large carbon (C) stocks, acting as a long-term atmospheric C sink since the last deglaciation. How these C-rich ecosystems will respond to future climate change, however, is still poorly understood. Furthermore, many northern peatlands exist in regions underlain by permafrost, adding to the challenge of projecting C balance under changing climate and permafrost dynamics. In this study, we used a paleoecological approach to examine the effect of past climates and local disturbances on vegetation and C accumulation at a peatland complex on the southern Seward Peninsula, Alaska over the past ∼15 ka (1 ka = 1000 cal yr BP). We analyzed two cores about 30 m apart, NL10-1 (from a permafrost peat plateau) and NL10-2 (from an adjacent thermokarst collapse-scar bog), for peat organic matter (OM), C accumulation rates, macrofossil, pollen and grain size analysis. A wet rich fen occurred during the initial stages of peatland development at the thermokarst site (NL10-2). The presence of tree pollen from Picea spp. and Larix laricinia at 13.5–12.1 ka indicates a warm regional climate, corresponding with the well-documented Bølling–Allerød warm period. A cold and dry climate interval at 12.1–11.1 ka is indicated by the disappearance of tree pollen and increase in Poaceae pollen and an increase in woody material, likely representing a local expression of the Younger Dryas (YD) event. Following the YD, the warm Holocene Thermal Maximum (HTM) is characterized by the presence of Populus pollen, while the presence of Sphagnum spp. and increased C accumulation rates suggest high peatland productivity under a warm climate. Toward the end of the HTM and throughout the mid-Holocene a wet climate-induced several major flooding disturbance events at 10 ka, 8.1 ka, 6 ka, 5.4 ka and 4.7 ka, as evidenced by decreases in OM, and increases in coarse sand abundance and aquatic fossils (algae Chara and water fleas Daphnia ). The initial peatland at permafrost site (NL10-1) is characterized by rapid C accumulation (66 g C m −2 yr −1 ), high OM content and a peak in Sphagnum spp. at 5.8–4.6 ka, suggesting the lack of permafrost. A transition to extremely low C accumulation rates of 6.3 g C m −2 yr −1 after 4.5 ka at this site suggests the onset of permafrost aggradation, likely in response to Neoglacial climate cooling as documented across the circum-Arctic region. A similar decrease in C accumulation rates also occurred at non-permafrost site NL10-2. Time-weighted C accumulation rates are 21.8 g C m −2 yr −1 for core NL10-1 during the last ∼6.5 ka and 14.8 g C m −2 yr −1 for core NL10-2 during the last ∼15 ka. Evidence from peat-core analysis and historical aerial photographs shows an abrupt increase in Sphagnum spp. and decrease in area of thermokarst lakes over the last century, suggesting major changes in hydrology and ecosystem structure, likely due to recent climate warming. Our results show that the thermokarst–permafrost complex was much more dynamic with high C accumulation rates under warmer climates in the past, while permafrost was stabilized and C accumulation slowed down following the Neoglacial cooling in the late Holocene. Furthermore, permafrost presence at local scales is controlled by both regional climate and site-specific factors, highlighting the challenge in projecting responses of permafrost peatlands and their C dynamics to future climate change.

Alaska↗

The making of the NEAM Tsunami Hazard Model 2018 (NEAMTHM18)

The NEAM Tsunami Hazard Model 2018 (NEAMTHM18) is a probabilistic hazard model for tsunamis generated by earthquakes. It covers the coastlines of the North-eastern Atlantic, the Mediterranean, and connected seas (NEAM). NEAMTHM18 was designed as a three-phase project. The first two phases were dedicated to the model development and hazard calculations, following a formalized decision-making process based on a multiple-expert protocol. The third phase was dedicated to documentation and dissemination. The hazard assessment workflow was structured in Steps and Levels. There are four Steps: Step-1) probabilistic earthquake model; Step-2) tsunami generation and modeling in deep water; Step-3) shoaling and inundation; Step-4) hazard aggregation and uncertainty quantification. Each Step includes a different number of Levels. Level-0 always describes the input data; the other Levels describe the intermediate results needed to proceed from one Step to another. Alternative datasets and models were considered in the implementation. The epistemic hazard uncertainty was quantified through an ensemble modeling technique accounting for alternative models’ weights and yielding a distribution of hazard curves represented by the mean and various percentiles. Hazard curves were calculated at 2,343 Points of Interest (POI) distributed at an average spacing of ∼20 km. Precalculated probability maps for five maximum inundation heights (MIH) and hazard intensity maps for five average return periods (ARP) were produced from hazard curves. In the entire NEAM Region, MIHs of several meters are rare but not impossible. Considering a 2% probability of exceedance in 50 years (ARP≈2,475 years), the POIs with MIH >5 m are fewer than 1% and are all in the Mediterranean on Libya, Egypt, Cyprus, and Greece coasts. In the North-East Atlantic, POIs with MIH >3 m are on the coasts of Mauritania and Gulf of Cadiz. Overall, 30% of the POIs have MIH >1 m. NEAMTHM18 results and documentation are available through the TSUMAPS-NEAM project website (http://www.tsumaps-neam.eu/), featuring an interactive web mapper. Although the NEAMTHM18 cannot substitute in-depth analyses at local scales, it represents the first action to start local and more detailed hazard and risk assessments and contributes to designing evacuation maps for tsunami early warning.

Frontiers in Earth Science↗

Constraining controls on carbonate sequences with high-resolution chronostratigraphy: Upper Miocene, Cabo de Gata region, SE Spain

A high-resolution chronostratigraphy has been developed for Miocene shallow-water carbonate strata in the Cabo de Gata region of SE Spain for evaluation of local, regional and global factors that controlled platform architecture prior to and during the Messinian salinity crisis. Paleomagnetic data were collected from strata at three localities. Mean natural remanent magnetization (NRM) ranges between 1.53 ?? 10-8 and 5.2 ?? 10-3 Am2/kg. Incremental thermal and alternating field demagnetization isolated the characteristic remanent magnetization (ChRM). Rock magnetic studies show that the dominant magnetic mineral is magnetite, but mixtures of magnetite and hematite occur. A composite chronostratigraphy was derived from five stratigraphic sections. Regional stratigraphic data, biostratigraphic data, and an 40Ar/39Ar date of 8.5 ?? 0.1 Ma, for an interbedded volcanic flow, place the strata in geomagnetic polarity Chrons C4r to C3r. Sequence-stratigraphic and diagenetic evidence indicate a major unconformity at the base of depositional sequence (DS)3 that contains a prograding reef complex, suggesting that approximately 250 000 yr of record (Subchrons C3Br.2r to 3Br.1r) are missing near the Messinian-Tortonian boundary. Correlation to the GPTS shows that the studied strata represent five third- to fourth-order DSs. Basal units are temperate to subtropical ramps (DS1A, DS1B, DS2); these are overlain by subtropical to tropical reefal platforms (DS3), which are capped by subtropical to tropical cyclic carbonates (Terminal Carbonate Complex, TCC). Correlation of the Cabo de Gata record to the Melilla area of Morocco, and the Sorbas basin of Spain indicate that early - Late Tortonian ramp strata from these areas are partially time-equivalent. Similar strata are extensively developed in the Western Mediterranean and likely were influenced by a cool climate or influx of nutrients during an overall rise in global sea-level. After ramp deposition, a sequence boundary (SB3) in Cabo de Gata correlates with a sequence boundary in Morocco and a published third-order eustatic fall suggesting at least a partial eustatic control for the sequence boundary. Coral reefs began to develop earlier in Cabo de Gata than at Melilla or Sorbas, arguing for local factors affecting this major environmental transition. Later Messinian reefs (DS3) from all areas are time-equivalent, suggesting a regional or global control on their formation. Some Halimeda-rich horizons in the Western Mediterranean are not time-equivalent event strata as hypothesized by others. Correlation of the relative sea-level curve for the fringing reef complex (DS3) with a published eustatic curve suggests at least a partial third-order global eustatic control for the highstand part of the sequence. Downstepping DS3 reefs and initial subaerial exposure of earlier DS3 reef strata approximately correlate with initiation of a series of subaerial unconformities in the South Pacific. The longer-term relative fall in sea-level during DS3 downstepping reef progradation does not correlate with a published third-order eustatic fall. Eustatic sea-level fluctuations may have been associated with initiation of the Mediterranean Messinian salinity crisis, but the longer-term fall may have been linked to tectonic uplift in the Mediterranean region. Widespread distribution of 'TCC-style' cycles of approximately the same age suggests a regional (Western Mediterranean) or global control on sea-level change responsible for TCC cycles. In addition, four subaerial exposure-capped TCC cycles may correlate with similar subaerial unconformities in the South Pacific, suggesting at least a partial eustatic control on TCC cyclicity. The high rates of relative sea-level change needed to generate a minimum of 25-30 m sea-level changes associated with each cycle are consistent with glacio-eustacy along with rapid evaporitic drawdown in the Mediterranean. ?? 2001 Elsevier Science B.V. All rights reserved.

Palaeogeography, Palaeoclimatology, Palaeoecology↗

A genetic signature of the evolution of loss of flight in the Galapagos cormorant

INTRODUCTION Changes in the size and proportion of limbs and other structures have played a key role in the evolution of species. One common class of limb modification is recurrent wing reduction and loss of flight in birds. Indeed, Darwin used the occurrence of flightless birds as an argument in favor of his theory of natural selection. Loss of flight has evolved repeatedly and is found among 26 families of birds in 17 different orders. Despite the frequency of these modifications, we have a limited understanding of their underpinnings at the genetic and molecular levels. RATIONALE To better understand the evolution of changes in limb size, we studied a classic case of recent loss of flight in the Galapagos cormorant ( Phalacrocorax harrisi ). Cormorants are large water birds that live in coastal areas or near lakes, and P. harrisi is the only flightless cormorant among approximately 40 extant species. The entire population is distributed along the coastlines of Isabela and Fernandina islands in the Galapagos archipelago. P. harrisi has a pair of short wings, which are smaller than those of any other cormorant. The extreme reduction of the wings and pectoral skeleton observed in P. harrisi is an attractive model for studying the evolution of loss of flight because it occurred very recently; phylogenetic evidence suggests that P. harrisi diverged from its flighted relatives within the past 2 million years. We developed a comparative and predictive genomics approach that uses the genome sequences of P. harrisi and its flighted relatives to find candidate genetic variants that likely contributed to the evolution of loss of flight. RESULTS We sequenced and de novo assembled the whole genomes of P. harrisi and three closely related flighted cormorant species. We identified thousands of coding variants exclusive to P. harrisi and classified them according to their probability of altering protein function based on conservation. Variants most likely to alter protein function were significantly enriched in genes mutated in human skeletal ciliopathies, including Ofd1 , Evc , Wdr34 , and Ift122 . We carried out experiments in Caenorhabditis elegans to confirm that a missense variant present in the Galapagos cormorant IFT122 protein is sufficient to affect ciliary function. The primary cilium is essential for Hedgehog (Hh) signaling in vertebrates, and individuals affected by ciliopathies have small limbs and ribcages, mirroring the phenotype of P. harrisi . We also identified a 4–amino acid deletion in the regulatory domain of Cux1 , a highly conserved transcription factor that has been experimentally shown to regulate limb growth in chicken. The four missing amino acids are perfectly conserved in all birds and mammals sequenced to date. We tested the consequences of this deletion in a chondrogenic cell line and showed that it impairs the ability of CUX1 to transcriptionally up-regulate cilia-related genes (some of which contain function-altering variants in P. harrisi ) and to promote chondrogenic differentiation. Finally, we show that positive selection may have played a role in the fixation of the variants associated with loss of flight in P. harrisi . CONCLUSION Our results indicate that the combined effect of variants in genes necessary for the correct transcriptional regulation and function of the primary cilium likely contributed to the evolution of highly reduced wings and other skeletal adaptations associated with loss of flight in P. harrisi . Our approach may be generally useful for identification of variants underlying evolutionary novelty from genomes of closely related species.

Science↗

Guidance documents: Continued support to improve operations of fish hatcheries and field sites to reduce the impact or prevent establishment of New Zealand Mudsnails and other invasive mollusks

This project tested and revised a risk assessment/management tool authored by Moffitt and Stockton designed to provide hatchery biologists and others a structure to measure risk and provide tools to control, prevent or eliminate invasive New Zealand mudsnails (NZMS) and other invasive mollusks in fish hatcheries and hatchery operations. The document has two parts: the risk assessment tool, and an appendix that summarizes options for control or management. The framework of the guidance document for risk assessment/hatchery tool combines approaches used by the Hazard Analysis and Critical Control Points (HACCP) process with those developed by the Commission for Environmental Cooperation (CEC), of Canada, Mexico, and the United States, in the Tri-National Risk Assessment Guidelines for Aquatic Alien Invasive Species. The framework approach for this attached first document assesses risk potential with two activities: probability of infestation and consequences of infestation. Each activity is treated equally to determine the risk potential. These two activities are divided into seven basic elements that utilize scientific, technical, and other relevant information in the process of the risk assessment. To determine the probability of infestation four steps are used that have scores reported or determined and averaged. This assessment follows a familiar HACCP process to assess pathways of entry, entry potential, colonization potential, spread potential. The economic, environmental and social consequences are considered as economic impact, environmental impact, and social and cultural influences. To test this document, the Principal Investigator worked to identify interested hatchery managers through contacts at regional aquaculture meetings, fish health meetings, and through the network of invasive species managers and scientists participating in the Western Regional Panel on Aquatic Nuisance Species and the 100th Meridian Initiative's Columbia River Basin Team, and the Western New Zealand Mudsnail Conference in Seattle. Targeted hatchery workshops were conducted with staff at Dworshak National Fish Hatchery Complex (ID), Similkameen Pond, Oroville WA, and Ringold Springs State Hatchery (WA). As a result of communications with hatchery staff, invasive species managers, and on site assessments of hatchery facilities, the document was modified and enhanced. Additional resources were added to keep it up to date. The result is a more simplified tool that can lead hatchery or management personnel through the process of risk assessment and provide an introduction to the risk management and communication process. In addition to the typical HACCP processes, this tool adds steps to rate and consider uncertainty and the weight of evidence regarding options and monitoring results . Uncertainty of outcome exists in most tools that can be used to control or prevent NZMS or other invasive mollusks from infesting an area. In additional this document emphasizes that specific control tools and plans must be tailored to each specific setting to consider the economic, environmental and social influences. From the testing and evaluation process, there was a strong recognition that a number of control and prevention tools previously suggested and reported in the literature from laboratory and small scale trials may not be compatible with regional and national regulations, economic constraints, social or cultural constraints, engineering or water chemistry characteristics of each facility. The options for control are summarized in the second document, Review of Control Measures for Hatcheries Infested with NZMS (Appendix A) that provides sources for additional resources and specific tools, and guidance regarding the feasibility and success of each approach. This tool also emphasizes that management plans need to be adaptive and incorporate oversight from professionals familiar with measuring risks of fish diseases, and treatments (e.g. the fish health practitioners and water quality and effluent management teams). Finally, with such a team, the adaptive management approach must be ongoing, and become a regular component of hatchery operations. Although it was the intent that this two part document would be included as part of the revised National Management and Control Plan for the NZMS proposed by the U.S. Fish and Wildlife Service (USFWS) and others, it is provided as a stand-alone document.

Cooperator Science Series↗

Fishes of Missouri River, chute, and flood plain habitats: Chapter 4 in Initial biotic survey of Lisbon Bottom, Big Muddy National Fish and Wildlife Refuge

The Lisbon Bottom Unit of the Big Muddy National Fish and Wildlife Refuge (Refuge) is approximately 2,200 acres and is the first complete unit of the Refuge. Primary objectives of the Refuge are to create and restore diverse riverine aquatic habitats and reconnect the Missouri River to its flood plain where feasible. Management seeks to accomplish these objectives by encouraging natural processes of erosion, deposition, and succession to the greatest extent possible. One of the most salient aquatic features of the Lisbon Bottom Unit is a newly created 2-mile-long free-flowing chute, or side channel (Fig. 1). This chute began forming as a levee breech scour hole during the Great Flood of 1993. The chute continued to develop during the 1995 flood and finally cut through to a flowing side channel during the 1996 flood. Extensive erosion and bank sluffing continued during 1997 due to sustained high flows that occurred throughout most of the year. The chute has progressively become wider and deeper with a developing meander pattern and channel bars have begun to form. Lisbon Bottom also contains several seasonal and permanent wetlands and is subject to periodic flooding at high Missouri River stages. Eight studies have been completed or are ongoing to evaluate Missouri River fishes associated with various habitat components of Lisbon Bottom and adjacent Missouri River reaches (Table 1). Several are part of much larger investigations to evaluate fish use of flood-created habitat features, basinwide fish assessment, and endangered or candidate species concerns. At the Lisbon Bottom Unit or in the Missouri River adjacent to the unit 54 fish species were collected (Table 2). Eight of these species have either a protected status under State or Federal laws or biologists consider them to potentially qualify for protected status. The status of the following fish is listed in Table 2: pallid sturgeon x shovelnose sturgeon hybrid, paddlefish, northern pike, sturgeon chub, sicklefin chub, ghost shiner, western silvery minnow, plains minnow, and blue sucker. Equally important to note are fish species that were not collected. The flathead chub, currently listed as State endangered and proposed for Federal listing, were not collected in any of the studies. Common fish species, collected in all seven sampling studies, included gizzard shad, common carp, river carpsucker, and freshwater drum. Paddlefish, skipjack herring, silver carp, white sucker, redfin shiner, western silvery minnow, bullhead minnow, and black bullhead were each collected in only one of the seven studies. Pflieger (1971) developed a guild system for the fishes of Missouri based on distribution patterns and centers of abundance. Fish were assigned to one of four primary faunal groups: Ozark, Big River, Lowland, and Prairie. Two secondary groups (Ozark-prairie and Ozark-lowland) were defined for species equally abundant and distributed in two of the primary areas. Two Ozark-prairie species, shorthead redhorse and white sucker, were collected on or near the Lisbon Bottom Unit. Big River species are found primarily in the Missouri and Mississippi Rivers. The environmental factors controlling fish distribution in these rivers appear to be substrate, current velocity, and turbidity. On or near the Lisbon Bottom Unit 18 Big River species were collected (Table 2). Lowland species are intolerant of siltation and high turbidity. They inhabit standing or slow-moving water with sand, fine gravel, and organic debris substrates. Two Lowland species, bullhead minnow and mosquitofish, were collected in connected scours on the Lisbon Bottom Unit. Prairie species have broader ecological tolerances than the Ozark and Lowland species. They are largely absent in high gradient streams and cool, clear waters. Prairie species make up a significant proportion of Missouri and Mississippi River fishes. Seven Prairie species were collected on or near the Lisbon Bottom Unit (Table 2). The last group, Wide-ranging, was defined as the species more widespread than the other faunal groups with broader environmental tolerances. Most Wide-ranging species occur at least occasionally in all sections of the state. Nineteen Wide-ranging species were collected on or near the Lisbon Bottom Unit (Table 2). Kubisiak (1997) found that isolated scours in the Lower Missouri River were dominated by Wide-ranging species. Connected scours also held Wideranging species but they contained Prairie species and were dominated by Big River species. Faunal group diversity and species richness were compared between isolated (Study l) and connected (Studies 4, 5, 7, and 8) scours of the Lisbon Bottom Unit. Species richness was similar with 22 species collected in isolated scours and an average of 24.8 species collected in connected scours. Isolated scour fish catches were dominated by Wide-ranging fish (59%) with 14% Big River fishes and 14% Prairie fishes. Connected scours had a larger percentage of Big River fish (32%) and fewer Wide-ranging (36%) and Prairie (13%) species than the isolated scours. In the combined studies, Wide-ranging fish dominated the connected scours. A similar comparison was made of fish from chute habitats and the Missouri River. Chute habitats include the Lisbon Bottom Chute, the chute between the right bank and island located at river mile 219, and a shallow channel running between a large sandbar complex and the channel border adjacent to Lisbon Bottom. Species richness in the chutes (35 species) was greater than that of the adjacent river (23 species). Big River species (40%) dominated the chute samples, followed by Wide-ranging at 31% Cand Prairie Cut 14%. The river samples contained equal numbers of Wide-ranging species (35%) and Big River species (35%) with 17% Prairie species. Differences in fish guild diversity and species richness may be due to numerous factors including sampling effort and sampling season. River samples were collected only in Study 3 in October; chute samples were collected in July, August, and October. Fish may have been using the chute in July and August to spawn, feed, or escape high river levels. The Lisbon Bottom Chute is a truly unique feature of the Lower Missouri River. It is the first naturally created side channel to develop since the U.S. Army Corps of Engineers tamed the river. The Corps eliminated 89% of the islands in the Missouri River between 1879 and 1954 (Funk and Robinson 1974). The chutes between the islands and the shore were shallower with less current than the main channel. These areas provided diversity to the Missouri River fish habitat, serving as nursery and feeding areas. The Lisbon Bottom Unit of the Big Muddy National Fish and Wildlife Refuge will continue to be studied in an effort to understand the complex role these areas play in river systems.

Missouri↗