USGS Science⌕ Search

SEARCH · USGS Science

Results for “Project Completion Report”

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

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

647 records · Page 36Linked to original sources

An initial abstraction and constant loss model, and methods for estimating unit hydrographs, peak streamflows, and flood volumes for urban basins in Missouri

Streamflow data, basin characteristics, and rainfall data from 39 streamflow-gaging stations for urban areas in and adjacent to Missouri were used by the U.S. Geological Survey in cooperation with the Metropolitan Sewer District of St. Louis to develop an initial abstraction and constant loss model (a time-distributed basin-loss model) and a gamma unit hydrograph (GUH) for urban areas in Missouri. Study-specific methods to determine peak streamflow and flood volume for a given rainfall event also were developed. Distinct basin characteristics were evaluated and selected for use on the basis of their theoretical relation to flow, results from previous studies, and the ability to reliably measure the basin characteristic using digital datasets and geographic information system (GIS) technology. The key basin characteristics determined or computed for each of the 39 basins upstream from the streamflow-gaging stations were drainage area, percent impervious area, main-channel slope based on the 10- and 85-percent length method, percentage of the basin area in storage (lakes, ponds, reservoirs, wetlands), the composite Natural Resources Conservation Service curve number estimated from a combination of the soil type data and land-use characteristics, and the streamflow variability index developed for the recently completed study of low-flow regression in Missouri. Characteristics of spatial and temporal rainfall distribution came from the next generation weather radar (NEXRAD) network. Procedures were developed for this study to convert the variable radar sweep rate into a 5-minute total rainfall hyetograph using data from the radar bin at the centroid of a given basin. Additional characteristics determined for each storm on the basin included the 5-day and 14-day antecedent rainfall, estimated from the mean of daily rainfall values from various rain gages in the area. The database of observed rainfall and runoff events for the 39 basins upstream from the streamflow-gaging stations was analyzed to compute the optimal storm-specific initial abstraction and constant loss values, as well as the time to peak, peak streamflow, and shape factor values of the GUH. The optimal storm-specific values were used to develop a regional regression equation for initial abstraction; conversely, the constant loss was estimated not by regression but from either a generalized or specific regional mean value. The optimal storm-specific values of GUH time to peak, GUH peak streamflow, and GUH shape factor were used to develop regression equations for the GUH. The regression equations for the GUH initially were tested alone, and then were combined with the appropriate regional regression equation for initial abstraction and both the generalized regional and specific regional mean constant loss values. For the GUH regression equations, the interquartile range was substantially smaller than the range spanned by the minimum and maximum values, which indicates most of the errors have much smaller variation, and the minimum and maximum values may be extreme outliers. The central tendency of the regressed errors for peak streamflow and runoff hydrograph volume were both approximately zero, which implies a generally unbiased estimation of these values. The mean and median of the regressed errors for time to peak streamflow were both small but greater than zero, which implies the GUH regression equations create a hydrograph that has a peak that is later in time than observed. Specifically, the regressed times indicate an offset of about 10 minutes, on average, from observed. The mean and median of the regressed errors for widths of the runoff hydrograph at 50 and 75 percent were less than zero, which implies the GUH tends to slightly underestimate these widths compared to the observed. The appropriate regional initial abstraction regression equation was combined with both the generalized and the specific regional mean constant loss values and the GUH regression equations. Both the generalized regional mean constant loss and specific regional mean constant loss forms of the basin-loss model worked equally well to model the observed runoff hydrograph based on the error analysis, and neither model seems to make a consistently better approximation. Both initial abstraction and constant loss models combined with the GUH regression equations were further validated using several storms available after the start of the project in early 2011 with similar but consistently higher error results. If these methods are used in an urban area in Missouri other than those examined in this study, advice to the user is given to consider using the generalized regional mean values. If these methods are used in an urban area that is a subbasin of one of the basins in this study, advice to the user is given to consider using the specific regional mean values. The rainfall-runoff pairs from the storm-specific GUH analysis were further analyzed against various basin and rainfall characteristics to develop equations to estimate the peak streamflow and flood volume based on a quantity of rainfall on the basin.

Missouri↗

Preliminary geologic map of the San Guillermo Mountain Quadrangle, Ventura County, California

New 1:24,000-scale geologic mapping in the Cuyama 30' x 60' quadrangle, in support of the USGS Southern California Areal Mapping Project (SCAMP), is contributing to a more complete understanding of the stratigraphy, structure, and tectonic evolution of the complex junction area between the NW-striking Coast Ranges and EW-striking western Transverse Ranges. The 1:24,000-scale geologic map of the San Guillermo Mountain quadrangle is one of six contiguous 7 1/2' quadrangle geologic maps in the eastern part of the Cuyama map area being compiled for a more detailed portrayal and reevaluation of geologic structures and rock units shown on previous geologic maps of the area (e.g., Dibblee, 1979). The following observations and interpretations are based on the new San Guillermo Mountain geologic compilation: (1) The new geologic mapping in the northern part of the San Guillermo Mountain quadrangle allows for reinterpretation of fault architecture that bears on potential seismic hazards of the region. Previous mapping had depicted the eastern Big Pine fault (BPF) as a northeast-striking, sinistral strike-slip fault that extends for 30 km northeast of the Cuyama River to its intersection with the San Andreas fault (SAF). In contrast the new mapping indicates that the eastern BPF is a thrust fault that curves from a northeast strike to an east strike, where it is continuous with the San Guillermo thrust fault, and dies out further east about 15 km south of the SAF. This redefined segment of the BPF is a south-dipping, north-directed thrust, with dominantly dip slip components (rakes > 60 deg.), that places Middle Eocene marine rocks (Juncal and Matilija Formations) over Miocene through Pliocene(?) nonmarine rocks (Caliente, Quatal, and Morales Formations). Although a broad northeast-striking fault zone, exhibiting predominantly sinistral components of slip (rakes < 45 deg.), extends to the SAF as previously mapped, the fault zone does not connect to the southwest with the BPF but instead curves into a southwest-directed thrust fault system a short distance north of the BPF. Oligocene to Pliocene(?) nonmarine sedimentary and volcanic rocks of the Plush Ranch, Caliente, and Morales(?) Formations are folded on both sides of this fault zone (informally named the Lockwood Valley fault zone [LVFZ] on the map). South-southeast of the LVFZ overturned folds have southward vergence. Several moderate-displacement (< 50 m), mainly northwest-dipping thrust and reverse faults, exhibiting mostly sinistral-oblique slip, flank and strike parallel to the overturned folds. The fold vergence and thrust direction associated with the LVFZ is opposite to that of the redefined BPF, providing further evidence that the two faults are distinct structures. These revised fault interpretations bring into question earlier estimates of net sinistral strike-slip displacement of as much as 13 km along the originally defined eastern BPF, which assumed structural connection with the LVFZ. Also, despite sparse evidence for repeated Quaternary movement on the LVFZ (e.g., Dibblee, 1982), the potential for a large earthquake involving coseismic slip on both the LVFZ and the central BPF to the southwest may not be as great as once believed. (2) Several generations of Pleistocene and younger dissected alluvial terrace and fan deposits sit at various levels above modern stream channels throughout the quadrangle. These deposits give testimony to the recent uplift and related fault deformation that has occurred in the area. (3) A vast terrane of Eocene marine sedimentary rocks (Juncal and Matilija Formations and Cozy Dell Shale) exposed south of the Big Pine fault forms the southern two-thirds of the San Guillermo Mountain quadrangle. Benthic foraminifers collected from various shale intervals within the Juncal Formation indicate a Middle Eocene age (Ulatisian) for the entire formation (K. McDougall, unpub. data, 1998) and deposition at paleodepths as great as 2,000 m (i.e., lowe

California↗

Assessment of hydrology, water quality, and trace elements in selected placer-mined creeks in the birch creek watershed near central, Alaska, 2001-05

Executive Summary The U.S. Geological Survey, in cooperation with the Bureau of Land Management, completed an assessment of hydrology, water quality, and trace-element concentrations in streambed sediment of the upper Birch Creek watershed near Central, Alaska. The assessment covered one site on upper Birch Creek and paired sites, upstream and downstream from mined areas, on Frying Pan Creek and Harrison Creek. Stream-discharge and suspended-sediment concentration data collected at other selected mined and unmined sites helped characterize conditions in the upper Birch Creek watershed. The purpose of the project was to provide the Bureau of Land Management with baseline information to evaluate watershed water quality and plan reclamation efforts. Data collection began in September 2001 and ended in September 2005. There were substantial geomorphic disturbances in the stream channel and flood plain along several miles of Harrison Creek. Placer mining has physically altered the natural stream channel morphology and removed streamside vegetation. There has been little or no effort to re-contour waste rock piles. During high-flow events, the abandoned placer-mine areas on Harrison Creek will likely contribute large quantities of sediment downstream unless the mined areas are reclaimed. During 2004 and 2005, no substantial changes in nutrient or major-ion concentrations were detected in water samples collected upstream from mined areas compared with water samples collected downstream from mined areas on Frying Pan Creek and Harrison Creek that could not be attributed to natural variation. This also was true for dissolved oxygen, pH, and specific conductance-a measure of total dissolved solids. Sample sites downstream from mined areas on Harrison Creek and Frying Pan Creek had higher median suspended-sediment concentrations, by a few milligrams per liter, than respective upstream sites. However, it is difficult to attach much importance to the small downstream increase, less than 10 milligrams per liter, in median suspended-sediment concentration for either basin. During low-flow conditions in 2004 and 2005, previously mined areas investigated on Harrison Creek and on Frying Pan Creek did not contribute substantial suspended sediments to sample sites downstream from the mined areas. No substantial mining-related water- or sediment-quality problems were detected at any of the sites investigated in the upper Birch Creek watershed during low-flow conditions. Average annual streamflow and precipitation were near normal in 2002 and 2003. Drought conditions, extreme forest fire impact, and low annual streamflow set apart the 2004 and 2005 summer seasons. Daily mean streamflow for upper Birch Creek varied throughout the period of record-from maximums of about 1,000 cubic feet per second to minimums of about 20 cubic feet per second. Streamflow increased and decreased rapidly in response to rainfall and rapid snowmelt events because the steep slopes, thin soil cover, and permafrost areas in the watershed have little capacity to retain runoff. Median suspended-sediment concentrations for the 115 paired samples from Frying Pan Creek and 101 paired samples from Harrison Creek were less than the 20 milligrams per liter total maximum daily load. The total maximum daily load was set by the U.S. Environmental Protection Agency for the upper Birch Creek basin in 1996. Suspended-sediment paired-sample data were collected using automated samplers in 2004 and 2005, primarily during low-flow conditions. Suspended-sediment concentrations in grab samples from miscellaneous sites ranged from less than 1 milligram per liter during low-flow conditions to 1,386 milligrams per liter during a high-flow event on upper Birch Creek. Streambed-sediment samples were collected at six sites on Harrison Creek, two sites on Frying Pan Creek, and one site on upper Birch Creek. Trace-element concentrations of mercury, lead, and zinc in streambed sedimen

Scientific Investigations Report↗

Comparison of simulations of land-use specific water demand and irrigation water supply by MF-FMP and IWFM

Two hydrologic models, MODFLOW with the Farm Process (MF-FMP) and the Integrated Water Flow Model (IWFM), are compared with respect to each model’s capabilities of simulating land-use hydrologic processes, surface-water routing, and groundwater flow. Of major concern among the land-use processes was the consumption of water through evaporation and transpiration by plants. The comparison of MF-FMP and IWFM was conducted and completed using a realistic hypothetical case study. Both models simulate the water demand for water-accounting units resulting from evapotranspiration and inefficiency losses and, for irrigated units, the supply from surface-water deliveries and groundwater pumpage. The MF-FMP simulates reductions in evapotranspiration owing to anoxia and wilting, and separately considers land-use-related evaporation and transpiration; IWFM simulates reductions in evapotranspiration related to the depletion of soil moisture. The models simulate inefficiency losses from precipitation and irrigation water applications to runoff and deep percolation differently. MF-FMP calculates the crop irrigation requirement and total farm delivery requirement, and then subtracts inefficiency losses from runoff and deep percolation. In IWFM, inefficiency losses to surface runoff from irrigation and precipitation are computed and subtracted from the total irrigation and precipitation before the crop irrigation requirement is estimated. Inefficiency losses in terms of deep percolation are computed simultaneously with the crop irrigation requirement. The seepage from streamflow routing also is computed differently and can affect certain hydrologic settings and magnitudes ofstreamflow infiltration. MF-FMP assumes steady-state conditions in the root zone; therefore, changes in soil moisture within the root zone are not calculated. IWFM simulates changes in the root zone in both irrigated and non-irrigated natural vegetation. Changes in soil moisture are more significant for non-irrigated natural vegetation areas than in the irrigated areas. Therefore, to facilitate the comparison of models, the changes in soil moisture are only simulated by IWFM for the natural vegetation areas, and soil-moisture parameters in irrigated regions in IWFM were specified at constant values . The IWFM total simulated changes in soil moisture that are related to natural vegetation areas vary from stress period to stress period but are small over the entire two-year period of simulation. In the hypothetical case study, IWFM simulates more evapotranspiration and return flows and less streamflow infiltration than MF-FMP. This causes more simulated surface-water diversions upstream and less simulated water available to downstream farms in IWFM compared to MF-FMP. The evapotranspiration simulated by the two models is well correlated even though the quantity is different. The different approaches used to simulate soil moisture, evapotranspiration, and inefficient losses yield different results for deep percolation and pumpage. In IWFM, deep percolation is a function of soil moisture; therefore, the constant soil-moisture requirement for irrigated regions, assumed for this comparison, results in a constant deep percolation rate. This led to poor correlation with the variable deep percolation rates simulated in MF-FMP, where the deep percolation rate, a fraction of inefficiency losses from precipitation and irrigation, is a function of quasi-steady state infiltration for each soil type and a function of groundwater head. Similarly, the larger simulated evapotranspiration in IWFM is mainly responsible for larger simulated groundwater pumpage demands and related lower groundwater levels in IWFM compared to MF-FMP. Because of the differences in features between MF-FMP and IWFM, the user may find that for certain hydrologic settings one model is better suited than the other. The performance of MF-FMP and IWFM in this particular hypothetical test case, with a fixed framework composed of common initial and boundary conditions and input parameter values, does not necessarily predict the performance of MF-FMP and IWFM in a real-world situation with variable framework and parameter values. These differences may affect the evaluation of policies, projects, or water-balance analysis for some hydrologic settings. Generally, both models are powerful tools that simulate a connected system of aquifer, stream networks, land surface, root zone, and runoff processes. MF-FMP simulated the hypothetical test case in about 4 minutes compared to about 58 minutes for IWFM.

Technical Information Record↗

Assessment of adult pallid sturgeon fish condition, Lower Missouri River—Application of new information to the Missouri River Recovery Program

During spring 2015, Nebraska Game and Parks Commission (NGPC) biologists noted that pallid sturgeon (Scaphirhynchus albus) were in poor condition during sampling associated with the Pallid Sturgeon Population Assessment Project and NGPC’s annual pallid sturgeon broodstock collection effort. These observations prompted concerns that reduced fish condition could compromise reproductive health and population growth of pallid sturgeon. There was a further concern that compromised condition could possibly be linked to U.S. Army Corps of Engineers management actions and increase jeopardy to the species. An evaluation request was made to the Missouri River Recovery Program and the Effects Analysis Team was chartered to evaluate the issue. Data on all Missouri River pallid sturgeon captures were requested and received from the National Pallid Sturgeon Database. All data were examined for completeness and accuracy; 12,053 records of captures between 200 millimeters fork length (mm FL) and 1,200 mm FL were accepted. We analyzed condition using (1) the condition formula (Kn) from Shuman and others (2011); (2) a second Kn formulation derived from the 12,053 records (hereafter referred to as “Alternative Kn”); and (3) an analysis of covariance (ANCOVA) approach that did not rely on a Kn formulation. The Kn data were analyzed using group (average annual Kn) and individual (percentage in low, normal, and robust conditions) approaches. Using the Shuman Kn formulation, annual mean Kn was fairly static from 2005 to 2011 (although always higher in the upper basin), declined from 2012 to 2015, then remained either static (lower basin) or increasing (upper basin) in 2016. Under the Alternative Kn formulation, the upper basin showed no decline in Kn, whereas the lower basin displayed the same trend as the Shuman Kn formulation. Using both formulations, the individual approach revealed a more complex situation; at the same times and locations that there are fish in poor condition, there are nearby fish in normal or robust condition. The ANCOVA approach revealed that fish condition at size changed between 400 and 600 mm and that some of the apparent trend in low condition was caused by differences in sample size across the size range of the population (that is, greater catch of intermediate-sized fish compared to large fish). We examined basin, year, origin (hatchery compared to wild), segment, and size class for effects on condition and concluded that, since 2012, there has been an increase in the percentage of pallid sturgeon in low condition. There are basin, year, and segment effects; origin and size class do not seem to have an effect. The lower basin, in particular segment 9 (Platte River to Kansas River), had a high percentage of low-condition fish. Within the segment, there were bend-level effects, but the bend effect was not spatially contiguous. We concluded that existing data confirm concerns about declining fish condition, especially in the segments between Sioux City, Iowa, and Kansas City, Missouri. Although the evidence is strong that fish condition has been in decline from 2011 to 2015, additional analysis of individual fish histories may provide more confidence in this conclusion; such analysis was beyond the scope of this effort but is part of our recommendations. The most recent data in 2016 indicate that decline of condition may have leveled off; however, the length of record is insufficient to determine whether recent declines are within the background range of variation. We recommend that monitoring of fish condition should be increased and enhanced with additional health metrics. We also recommend that, should condition continue to decline, processes are deployed to bring low-condition adult fish into the hatchery to improve nutrition and condition. We could not determine the cause of declining fish condition with available data, but we compiled information on several dominant hypotheses in two main categories: inter- or intraspecific competition for resources and habitat conditions. Data are insufficient to indicate a specific causation or solution, and it is possible that multiple causes apply. We make recommendations for additional research that can be pursued to address uncertainties in trends in fish health as well as potential causes.

Missouri River↗

Hydrological information products for the Off-Project Water Program of the Klamath Basin Restoration Agreement

The Klamath Basin Restoration Agreement (KBRA) was developed by a diverse group of stakeholders, Federal and State resource management agencies, Tribal representatives, and interest groups to provide a comprehensive solution to ecological and water-supply issues in the Klamath Basin. The Off-Project Water Program (OPWP), one component of the KBRA, has as one of its purposes to permanently provide an additional 30,000 acre-feet of water per year on an average annual basis to Upper Klamath Lake through "voluntary retirement of water rights or water uses or other means as agreed to by the Klamath Tribes, to improve fisheries habitat and also provide for stability of irrigation water deliveries." The geographic area where the water rights could be retired encompasses approximately 1,900 square miles. The OPWP area is defined as including the Sprague River drainage, the Sycan River drainage downstream of Sycan Marsh, the Wood River drainage, and the Williamson River drainage from Kirk Reef at the southern end of Klamath Marsh downstream to the confluence with the Sprague River. Extensive, broad, flat, poorly drained uplands, valleys, and wetlands characterize much of the study area. Irrigation is almost entirely used for pasture. To assist parties involved with decisionmaking and implementation of the OPWP, the U.S. Geological Survey (USGS), in cooperation with the Klamath Tribes and other stakeholders, created five hydrological information products. These products include GIS digital maps and datasets containing spatial information on evapotranspiration, subirrigation indicators, water rights, subbasin streamflow statistics, and return-flow indicators. The evapotranspiration (ET) datasets were created under contract for this study by Evapotranspiration, Plus, LLC, of Twin Falls, Idaho. A high-resolution remote sensing technique known as Mapping Evapotranspiration at High Resolution and Internalized Calibration (METRIC) was used to create estimates of the spatial distribution of ET. The METRIC technique uses thermal infrared Landsat imagery to quantify actual evapotranspiration at a 30-meter resolution that can be related to individual irrigated fields. Because evaporation uses heat energy, ground surfaces with large ET rates are left cooler as a result of ET than ground surfaces that have less ET. As a consequence, irrigated fields appear in the Landsat images as cooler than nonirrigated fields. Products produced from this study include total seasonal and total monthly (April-October) actual evapotranspiration maps for 2004 (a dry year) and 2006 (a wet year). Maps showing indicators of natural subirrigation were also provided by this study. "Subirrigation" as used here is the evapotranspiration of shallow groundwater by plants with roots that penetrate to or near the water table. Subirrigation often occurs at locations where the water table is at or above the plant rooting depth. Natural consumptive use by plants diminishes the benefit of retiring water rights in subirrigated areas. Some agricultural production may be possible, however, on subirrigated lands for which water rights are retired. Because of the difficulty in precisely mapping and quantifying subirrigation, this study presents several sources of spatially mapped data that can be used as indicators of higher subirrigation probability. These include the floodplain boundaries defined by stream geomorphology, water-table depth defined in Natural Resources Conservation Service (NRCS) soil surveys, and soil rooting depth defined in NRCS soil surveys. The two water-rights mapping products created in the study were "points of diversion" (POD) and "place of use" (POU) for surface-water irrigation rights. To create these maps, all surface-water rights data, decrees, certificates, permits, and unadjudicated claims within the entire 1,900 square mile study area were aggregated into a common GIS geodatabase. Surface-water irrigation rights within a 5-mile buffer of the study area were then selected and identified. The POU area was then totaled by water right for primary and supplemental water rights. The maximum annual volume (acre-feet) allowed under each water right also was calculated using the POU area and duty (allowable annual irrigation application in feet). In cases where a water right has more than one designated POD, the total volume for the water right was equally distributed to each POD listed for the water right. Because of this, mapped distribution of diversion rates for some rights may differ from actual practice. Water-right information in the map products was from digital datasets obtained from the Oregon Water Resources Department and was, at the time acquired, the best available compilation of water-right information available. Because the completeness and accuracy of the water-right data could not be verified, users are encouraged to check directly with the Oregon Water Resources Department where specific information on individual rights or locations is essential. A dataset containing streamflow statistics for 72 subbasins in the study area was created for the study area. The statistics include annual flow durations (5-, 10-, 25-, 50-, and 95-percent exceedances) and 7-day, 10-year (7Q10) and 7-day, 2-year (7Q2) low flows, and were computed using regional regression equations based on measured streamflow records in the region. Daily streamflow records used were adjusted as needed for crop consumptive use; therefore the statistics represent streamflow under more natural conditions as though irrigation diversions did not exist. Statistics are provided for flow rates resulting from streamflow originating from within the entire drainage area upstream of the subbasin pour point (referring to the outlet of the contributing drainage basin). The statistics were computed for the purpose of providing decision makers with the ability to estimate streamflow that would be expected after water conservation techniques have been implemented or a water right has been retired. A final product from the study are datasets of indicators of the potential for subsurface return flow of irrigation water from agricultural areas to nearby streams. The datasets contain information on factors such as proximity to surface-water features, geomorphic floodplain characteristics, and depth to water. The digital data, metadata, and example illustrations for the datasets described in this report are available on-line from the USGS Water Resources National Spatial Data Infrastructure (NSDI) Node Website http://water.usgs.gov/lookup/getgislist or from the U.S. Government website DATA.gov at http://www.data.gov with links provided in a Microsoft&reg; Excel&reg; workbook in appendix A.

California;Oregon↗

The potential effects of sodium bicarbonate, a major constituent from coalbed natural gas production, on aquatic life

The production water from coalbed natural gas (CBNG) extraction contains many constituents. The U.S. Environmental Protection Agency has established aquatic life criteria for some of these constituents, and it is therefore possible to evaluate their risk to aquatic life. However, of the major ions associated with produced waters, chloride is the only one with an established aquatic life criterion (U.S. Environmental Protection Agency, 1988). The focus of this research was NaHCO 3 , a compound that is a major constituent of coalbed natural gas produced waters in the Tongue and Powder River Basins. This project included laboratory experiments, field in situ experiments, a field mixing zone study, and a fishery presence/absence assessment. Though this investigation focuses on the Tongue and Powder River Basins, the information is applicable to other watersheds where sodium bicarbonate is a principle component of product water either from CBNG or from traditional or unconventional oil and gas development. These data can also be used to separate effects of saline discharges from those potentially posed by other constituents. Finally, this research effort and the additional collaboration with USGS Water Resources and Mapping, Bureau of Land Management, US Environmental Protection Agency, State of Montana, State of Wyoming, Montana State University, University of Wyoming, and others as part of a Powder River Aquatic Task Group, can be used as a model for successful approaches to studying landscapes with energy development. The laboratory acute toxicity experiments were completed with a suite of organisms, including 7 species of fish, 5 species of invertebrates, and 1 amphibian species. Experiments performed on these multiple species resulted in LC50s that ranged from 1,120 to greater than (>) 8,000 milligrams sodium bicarbonate per liter (mg NaHCO 3 /L) (also defined as 769 to >8,000 milligrams bicarbonate per liter (mg HCO 3 -/L) or total alkalinity expressed as 608 to >4,181 milligrams calcium carbonate per liter (mg CaCO 3 /L)) that varied across species and lifestage within a species. The age at which fish were exposed to NaHCO 3 significantly affected the severity of toxic responses for some organisms. The chronic toxicity of NaHCO 3 was defined in experiments that lasted from 7&mdash;60 days post-hatch. For these experiments, sublethal effects such as growth and reproduction, in addition to significant reductions in survival were included in the final determination of effects. Chronic toxicity was observed at concentrations that ranged from 450 to 800mg NaHCO 3 /L (also defined as 430 to 657 mg HCO 3 -/L or total alkalinity expressed as 354 to 539 mg CaCO 3 /L) and the specific concentration depended on the sensitivity of the four species of invertebrates and fish exposed. Sublethal investigations during chronic studies revealed percent decrease in the activity of sodium-potassium adenosine triphosphatase (Na/K ATPase, an enzyme involved in ionoregulation) and the age of the fish at the onset of the decrease may affect the ability of fathead minnow to survive exposures to NaHCO 3 . A database of toxicity evaluations of NaHCO 3 on aquatic life has been constructed. Using these data, sample acute and chronic criteria of 459 and 381 mg NaHCO 3 /L, respectively, can be calculated for the protection of aquatic life. The final derivation and implementation of such criteria is, of course, left to the discretion of the concerned management agencies. A combination of in situ experiments, static-renewal experiments performed simultaneously with in situ experiments, and static renewal experiments performed with site water in the laboratory, demonstrated that untreated coalbed natural gas (CBNG) product water from the Tongue and Powder River Basins reduces survival of fathead minnow and pallid sturgeon. More precisely, the survival of early-lifestage fathead minnow, especially those less than 6-days post hatch (dph), likely is reduced significantly in the field when concentrations of NaHCO 3 rise above 1,500 mg/L. However, age was not a factor for pallid sturgeon and they were sensitive to product water regardless of age. Treatment with the Higgins Loop&trade; technology and dilution of untreated water increased survival in the laboratory. Both of these situations reduced ammonia in addition to the concentrations of NaHCO 3 . These experiments addressed the acute toxicity of effluent waters being added to the main stem rivers, but did not address issues related to the volumes of water that may be added to the watershed. Mixing zones of the three outfalls studied ranged from approximately 800&mdash;1,200 m below the confluence and the areas within these mixing zones with acutely lethal concentrations of NaHCO 3 (as defined by the presence of concentrated dye) are limited. The areas with concentrations of NaHCO 3 more than the concentrations likely to cause significant mortality, and more than the calculated sample water-quality criteria in the Tongue and Powder River Basins appear to be limited to tributaries and parts of mixing zones with considerable additions of untreated effluent.

Montana;Wyoming↗

Results of 2018–19 water-quality and hydraulic characterization of aquifer intervals using packer tests and preliminary geophysical-log correlations for selected boreholes at and near the former Naval Air Warfare Center Warminster, Bucks County, Pennsylvania

The U.S. Geological Survey (USGS) collected data on the vertical distribution of hydraulic head, specific capacity, and water quality using aquifer-interval-isolation tests and other vertical profiling methods in 15 boreholes completed in fractured sedimentary bedrock in Northampton, Warminster, and Warwick Townships, Bucks County, Pennsylvania during 2018–19. This work was done, in cooperation with the U.S. Navy, to support detailed investigations at and near the former Naval Air Warfare Center (NAWC) Warminster, where groundwater contamination with per- and polyfluoroalkyl substances (PFAS) had become a concern since 2014. Two PFAS compounds, perfluorooctanesulfonic acid (PFOS) and perfluorooctanoic acid (PFOA), have been measured in groundwater samples from supply and monitoring wells at or near NAWC Warminster in concentrations above U.S. Environmental Protection Agency health advisory levels for drinking water. The area is underlain by the Triassic Stockton Formation, which predominantly consists of sandstone interbedded with shale and siltstone beds and forms a layered fractured-rock aquifer used for private, industrial, and public drinking water supply. The vertical distribution of aquifer properties and water quality was assessed through hydraulic tests and sampling of aquifer intervals using a straddle-packer system (13 boreholes) or depth-discrete point sampling under known borehole-flow conditions (2 boreholes). Geophysical and video logs collected by USGS during 2017–19 were used to identify potential water-bearing fractures in 15 boreholes, which ranged in depth from 210 to 604 feet (ft) and included 6 boreholes drilled in 2018 and 9 existing wells on or near the former NAWC Warminster. Measured borehole flow was predominantly downward in most of the deepest boreholes (greater than 400 ft), which were commonly located at the highest land-surface elevations, with inflow from fractures at relatively shallow depths and outflow through fractures near or below depths of 500 ft below land surface. Hydraulic head differences measured during packer tests were up to about 60 ft between shallow and deep intervals. Borehole flow was predominantly upward in most boreholes less than 400 ft in depth and farther from, and at lower land-surface elevations than, the former NAWC Warminster. Total borehole specific capacity ranged from about 0.07 to 41 gallons per minute per foot [(gal/min)/ft]. Specific-capacity values for individual intervals ranged from 0.02 to 40.0 (gal/min)/ft, with a median of 1.14 (gal/min)/ft and a large range in values at most depths. Differences in water quality of samples as indicated by field properties (pH, dissolved oxygen, and specific conductance) and concentrations of dissolved major ions, PFOA, and PFOS were apparent among isolated intervals in the boreholes. Summed concentrations of PFOA and PFOS ranged from about 11 to 10,780 nanograms per liter (ng/L) and were greater than the 2016 U.S. Environmental Protection Agency health advisory of 70 ng/L for summed PFOA and PFOS concentrations in 62 of 104 intervals and discrete depths tested. The mass ratio of PFOS to PFOA was generally higher than 1.0 in samples with summed PFOA and PFOS concentrations greater than 70 ng/L, with ratio values as high as 8.7. In many boreholes, summed concentrations of PFOA and PFOS were positively related to chloride concentrations, which were elevated above natural-background values [less than 10 milligrams per liter] in most samples and as high as 717 milligrams per liter. Sources of the elevated chloride other than, or in addition to, common rock salt (sodium chloride) were indicated by chloride to sodium molar ratios greater than 1.0. Water-quality data indicated that sampled water from some intervals with lower hydraulic heads may be affected by water from intervals with higher hydraulic heads because of vertical flow in open boreholes; samples from these intervals with lower hydraulic heads may not be fully representative due to some component of cross contamination and should be interpreted with caution. Through a preliminary correlation of natural gamma and resistivity logs of boreholes drilled at and near the former NAWC Warminster, 11 lithologic units were identified and interpreted to strike northeast and dip to the northwest. Hydraulic heads were generally highest in isolated intervals that intercepted beds which, when projected up dip, crop out at the highest land-surface elevation on the former NAWC Warminster, indicating that the dipping-bed structure and topography are factors affecting the distribution of hydraulic head in the aquifer. The hydrogeologic framework in conjunction with the vertical distribution of hydraulic heads and water quality may assist in evaluating the locations of various PFAS sources and potential migration pathways of PFAS in groundwater at and near NAWC Warminster.

Pennsylvania↗

Modeling, mapping, and measuring the risk of freshwater invasive species across Alaska

Freshwater ecosystems of the Alaskan Arctic and Subarctic provide resources that are culturally, ecologically, and economically invaluable. Presently, these regions are relatively free of the impacts from invasive species compared to southern latitudes. To date, there have been relatively few verified introductions of aquatic invasive species (AIS) to freshwater ecosystems in Alaska. The expanding list and distribution of AIS has led to significant negative ecological and economic impacts (e.g., waterweed Elodea nuttalli ; E. canadensis and northern pike Esox Lucius introduced outside its native range in Alaska). Escalating human activity across Alaskan lands and waters, coupled with rapidly shifting environmental conditions, increases the potential for new species introductions and subsequent establishment. Creating a proactive framework for well-informed decision-making and action can improve the effectiveness of prevention efforts and bolster decision support tools that help resource managers direct limited resources. Prioritizing AIS that may be introduced and become established, as well as the locations at highest risk of invasion, is foundational to building a proactive invasive species management framework in Alaska. This project sought to identify and prioritize AIS known to be invasive in the contiguous United States, evaluate current and future habitat suitability for AIS in Alaska, and assess potential for AIS to be transported to habitats across Alaska, utilizing similar assessment methods as implemented for Bering Sea marine invasive species and non-native plants in Alaska. To accomplish this goal, the objectives of the project were to: 1) develop a formal ranked list of potential AIS to freshwater systems of Alaska; 2) assess the level of establishment risk for potential AIS by developing habitat suitability models for waterbodies across Alaska; and 3), identify potential pathways and specific vectors for high-risk AIS to invade Alaska and develop a framework for how vector analysis will be completed to understand transport risk. Overall, our goal is horizon scanning which is defined by Roy et al. (2019) as “a systematic examination of potential threats and opportunities, within a given context, and likely future developments, which are at the margin of current thinking and planning.” The scans include pathway analyses and risk screening of species present at pathway origin points, with a focus on identifying species at high risk of being introduced, becoming established, spreading, and causing harm. We refined a list of 28 AIS from a list of hundreds based on characterizations of species’ invasiveness and species’ proximity to Alaska (USGS 2020; GBIF 2022). Next, we evaluated the relative invasiveness of individual species to create an initial AIS ranking. We sought to characterize habitat suitability of AIS by selecting variables that were continental in scale, covering North America to include Alaska as well as the lower 48 states comparing natural discharge, sub-basin average terrain slope (degrees), average silt fraction, average organic carbon, lithological class, and human footprint in sub-basin in 2009. We estimated AIS habitat suitability across the entire state of Alaska using the physiological tolerances of the AIS (Appendix 2). We also evaluated pathways and vectors for the introduction of AIS (Appendix 2). Many pathways and vectors considered did not meet the criteria for Alaska or freshwater systems. Of the 28 ranked species that we categorized as very high, high, and moderate levels of invasiveness; all three risk groups included fish and mollusks (Appendix 2). One commonality of the very high-invasiveness-ranked species was the availability of Ecological Risk Screening Summary documents (USFWS, 2022) produced by U.S. Fish and Wildlife Service (USFWS), except for the goldfish ( Carassius auratus ) and the New Zealand mudsnail ( Potamopyrgus antipodarum ). The Ecological Risk Screening Summary is now available for New Zealand mudsnails. In general, fish species often ranked very high or high in invasiveness and included sportfish and aquarium fish, suggesting the importance of pathways such as aquarium trade, fishing industry, intentional (but illegal) introductions of sportfishes and aquarium fishes for establishment. The technique we used for habitat suitability models necessitated aquatic environmental datasets that were continental in scale, which was often interpolated from very coarse resolution source data layers, particularly in Alaska. Better spatial data representing aquatic environments would likely improve this approach. While the lack of introductions in Alaska and nearby provinces and states is encouraging, the lack of occurrence data for the focal species also created complications for habitat suitability modeling. Despite the challenges, the habitat suitability models indicated limited suitability for warmwater species while some species, such as Brook trout (S alvelinus fontinalis ), have high habitat suitability across Alaska no matter what threshold approach is taken. Some environmental predictors were more important than others. Specifically, the most important predictor variable, ‘frost free days,’ was critical for 15 out of 28 species as expected due to harsh winter conditions in Arctic and Subarctic regions. The second most important predictor was ‘subbasin land surface runoff’, a variable that indicates the amount of discharge and runoff, while the third most important predictor was ‘snow cover’ another indication of winter conditions. Overall, the ability to understand the effect of future climate scenarios on the establishment of AIS was challenging. A detailed dataset of freshwater temperatures and water chemistry (e.g., pH, calcium) would greatly improve the ability to predict invasiveness of freshwater species to Alaska’s ecosystems on a regional basis. Future studies may benefit from a more focused geographic scope examining a group of subbasins or a regional basin rather than the entire state. These drainages could be selected based upon the mostly likely locations of introduction pathways. The two most prevalent pathway risks for AIS are in-state transfer and stowaways/contaminants. Although there are examples of introductions from other pathways, the risk is somewhat mitigated by Alaska’s climate and regulations. However, variable application of protocols for inspection and cleaning of fishing gear, watercraft, and other similar items while traveling into Alaska as well as transferring from waterbody to waterbody within the state creates a substantial risk in introducing invasive species. We plot cumulative invasive vulnerability for all subbasins and for the top 10% of subbasins (Appendix 3).

Alaska↗

Ground-water resources of Riverton irrigation project area, Wyoming, with a section on chemical quality of ground water

The Riverton irrigation project area is in the northwestern part of the Wind River basin in west-central Wyoming. Because the annual precipitation is only about 9 inches, agriculture, which is the principal occupation in the area, is dependent upon irrigation. Irrigation by surface-water diversion was begum is 1906; water is now supplied to 77,716 acres and irrigation has been proposed for an additional 31,344 acres. This study of the geology and ground-water resources of the Riverton irrigation project, of adjacent irrigated land, and of nearby land proposed for irrigation was begun during the summer of 1948 and was completed in 1951. The purpose of the investigation was to evaluate the ground-water resources of the area and to study the factors that should be considered in the solution of drainage and erosional problems within the area. The Riverton irrigation project area is characterized by flat to gently sloping stream terraces, which are flanked by a combination of badlands, pediment slopes, and broad valleys. These features were formed by long-continued erosion in an arid climate of the essentially horizontal, poorly consolidated beds of the Wind River formation. The principal streams of the area flow south-eastward. Wind River and Fivemile Creek are perennial streams and the others are intermittent. Ground-water discharge and irrigation return flow have created a major problem in erosion control along Fivemile Creek. Similar conditions might develop along Muddy and lower Cottonwood Creeks when land in their drainage basins is irrigated. The bedrock exposed in the area ranges in age from Late Cretaceous to early Tertiary (middle Eocene). The Wind River formation of early and middle Eocene age forms the uppermost bedrock formation in the greater part of the area. Unconsolidated deposits of Quaternary age, which consist of terrace gravel, colluvium, eolian sand and silt. and alluvium, mantle the Wind River formation in much of the area. In the irrigated parts of the project, water from domestic use is obtained chiefly from the sandstone beds of the Wind River formation although some is obtained from the alluvium underlying the bottom land and from the unconsolidated deposits underlying the lower terraces along the Wind River. Although adequate quantities if water for domestic use are available from the Wind River formation, there quantities are not considered to be large enough to warrant pumping of ground water for irrigation. Only a few wells are in the nonirrigated part of the area. When this new land is irrigated, a body of ground water will gradually form in the terrace deposits and the alluvial and colluvial-alluvial deposits. Eventually, the terrace deposits may yield adequate quantities of water for domestic and stock use, but only locally are the alluvial and colluvial-alluvial deposits likely to become suitable aquifers. In the Riverton irrigation project area, ground water occurs under water-table conditions near the surface and under artesian conditions in certain strata at both shallow and greater depths. Irrigation is the principal source of recharge to the shallow aquifers; the water level in wells that tap these aquifers fluctuates with irrigation. The depth to water in the shallow wells ranges from less than 1 foot to about 30 feet below the land surface, depending on the season of the year and on the length of time the land has been irrigated. The water level in the wells that tap the deep confined aquifers , which receive recharge indirectly from surface sources, fluctuates only slightly because the recharge and discharge are more constant. In most places the depth to water in wells penetrating the deep confined aquifers is mush greater than that in shallow wells. but in certain low areas water from the deep aquifers flows at the surface from wells. Ground water moves from the area of recharge in the direction of the hydraulic gradient and is discharges either by evapotranspiration; by inflow into streams, drains, or lakes; by pumping or flow of wells; or by flow of springs. Waterlogging and the associated development of saline soils are common in parts of the Riverton irrigation project and adjacent irrigated land. The waterlogging is in part the result of the infiltration of irrigation water in excess of the capacity of the aquifers to store and transmit this added recharge. The solution of the drainage problems involves the consideration of a number of factors, some of which are inadequately known in some parts of the area and require further investigation before fully effective drainage measures can be designed. The results of an aquifer test to determine the hydrologic characteristics of the Wind River formation at Riverton indicate a transmissibility of 10,000 gallons per day per foot (10,000 gpd per ft) and a storage coefficient of 2 x 10 -4 . The results of the test provide a part of the necessary foundation for the solution of present and future water-supply problems at Riverton and throughout the project area. Water from shallow aquifers in irrigated tracts in the Riverton irrigation project area generally contains large amounts of dissolved solids that were leached from the soil and rocks by infiltrating irrigation water. However, wells tapping beds that receive considerable recharge from influent canal and drain seepage yield water of relatively low mineralizatoin. Dilute water is obtained also from some shallow wells in the alluvial bottom lands and on low stream terraces that border the Wind Rover. Water from deep aquifers generally is more dilute than that from shallow aquifers. However, ground water from the deep aquifers, unmixed with irrigation water, generally has a percent sodium greater than 80. Analyses of salt crusts on the ground surface in low areas that are affected by effluent seepage and a high water table show predominance of sodium sulfate salinity, and from determinations of the water-soluble and acid-soluble substances in several samples of soil and shale it is apparent that harmful concentrations of salts are being deposited in poorly drained area. Although most of the soul in the Midvale irrigation district is of the normal arid type, analyses of soil samples show that saline, nonsaline alkaline, and saline alkaline types also are present.

Wyoming↗

Geologic evolution, sedimentation, and paleoenvironments of the Angola Basin and adjacent Walvis Ridge: Synthesis of results of Deep Sea Drilling Project Leg 75

The section recovered at Site 530 (Holes 53OA and 530B) consists of eight sedimentary units and one basalt unit. The composition of the basalt recovered in Hole 53OA is distinct from typical mid-ocean ridge basalts (MORBs) but is similar to that of Hawaiian tholeiites and basalt from the central part of Walvis Ridge. Throughout most of its history, the southern Angola Basin received large volumes of redeposited material in the form of turbidites and, most recently, debris-flow deposits. Most of this material was derived from Walvis Ridge to the south, but thickness trends of acoustic units suggest that some of the sediment was derived from the African continental margin to the east. The basal sedimentary unit (Albian to Santonian) at Site 530 contains 262 beds of black shale that are interbedded with green and red claystone. Black shale makes up less than 10% of the total section, but in two cores of early Turonian age, black shale beds compose about 50% of the section. The black shales contain up to 19% organic carbon (average of about 5%) that is mainly of autochthonous marine origin but with significant contributions from terrigenous organic matter. The origin of these more- and less-reduced interbedded lithologies with varying amounts and types of organic matter, and variable amounts of pelagic, hemipelagic, and turbiditic sediment is complex and cannot be explained by any one simple process. Many factors affecting the concentration of dissolved oxygen in the bottom waters of the Angola Basin varied throughout the middle Cretaceous to produce bottom-water conditions that fluctuated between mildly oxic and oxygen-deficient, but most of the time bottom-waters and sediment-interstitial waters were sufficiently oxic to permit the accumulation of red oxidized sediment. A relatively complete sedimentary record of the Cretaceous/Tertiary boundary was recovered within a sequence of mudstone and marlstone turbidites in Hole 530A. There is a significant increase in the concentration of iridium above background levels at the boundary. High concentrations of many other elements also occur within the same stratigraphic interval as the iridium anomaly. Furthermore, there is a marked decrease in CaCO3 in the Tertiary strata above the iridium anomaly which suggests that the production of shallow-water carbonate also may have been affected by whatever caused elevated concentrations of iridium and other elements. These observations are consistent with the asteroid-impact theory proposed to explain the worldwide occurrence of an iridium anomaly at the Cretaceous/Tertiary boundary. The Cenozoic history of the Angola Basin was controlled mainly by (1) restriction of bottom-water flow from the south by Walvis Ridge; (2) development of glaciation on Antarctica; (3) opening of circulation passages in the southern oceans; (4) rapid turnover of cold, nutrient-rich waters that resulted in high productivity of diatoms; (5) influx of terrigenous sediment mainly by turbidity currents; and (6) production and preservation of carbonate sediment. The most distinctive Cenozoic event recorded in the section at Site 530 is the beginning of extensive glaciation on Anarctica and concomitant initiation of modern thermohaline bottom-water circulation that is manifested as a middle Eocene to middle Oligocene unconformity or compressed section accompanied by a drastic decrease in accumulation of CaCO3. Diatom abundances in HPC cores from Walvis Ridge (Site 532) and Angola Basin (Hole 53OB) indicate that Benguela upwelling in these areas began in the late Miocene, reached a peak in the late Pliocene to early Pleistocene, and declined thereafter. Short-term variations in sediment composition at Site 532 are manifested as cyclic variations in concentrations of clay, CaCO3, and organic carbon with average periodicities of about 30-60 k.y. The main variability that produced the cycles probably was the influx of terrigenous clastic material which diluted the CaCO3. The sediment at Site 532 also contains several percent organic carbon that is dominantly of marine origin, but with significant terrigenous components. Data from multichannel seismic, gravity, and magnetic surveys were used to define the regional stratigraphic and structural evolution of Walvis Ridge and adjacent Cape and Angola basins. Six structural provinces are recognized, four on Walvis Ridge and two additional provinces that correspond to the Cape and Angola basins. The two eastern structural provinces on Walvis Ridge are underlain by continental crust. The two western structural provinces are underlain by oceanic basement. Two main directions of faults are evident in seismic profiles, one trending N 10° and one trending N 60°. The N 60° trend corresponds to the general orientation of the northern and southern flanks of Walvis Ridge as well as to the dominant direction of fracture zones. During the first phase of separation of Africa from South America (ca. 120-130 m.y. ago), a voluminous mass of volcanics was emplaced simultaneous with the emplacement of basalt in the Parana Basin of Brazil and the Kaokoveld Region of South Africa. This period of volcanism also formed the series of seaward-dipping internal basement reflectors that are characteristic of the two structural provinces of Walvis Ridge. A system of fault blocks developed in the brittle upper part of the newly formed crust. During the second phase of rifting, which ended before late Aptian, more tilted fault blocks were created in the upper brittle stratified continental crust. Magnetic lineations in basement rocks in the Angola and Cape basins in the vicinity of Walvis Ridge are not distinct but suggest that oceanic crust began to be emplaced between 120 and 112 m.y. ago (Barremian to early Aptian). At least part of the oceanic crust of the central plateau of eastern Walvis Ridge (structural province 3) may have been emplaced before any oceanic crust formed in the adjacent basins. A ridge jump occurred during the late Aptian to early Albian in the southern part of the Angola Basin which translated the previously formed oceanic crust and its overlying evaporite deposits on the South American side. Several ridge jumps occurred on both sides of Walvis Ridge during the Late Cretaceous and early Tertiary to produce a 500-km-long segment of mid-ocean ridge.

Initial Reports of the D.S.D.P.↗

Preliminary simulation of chloride transport in the Equus Beds aquifer and simulated effects of well pumping and artificial recharge on groundwater flow and chloride transport near the city of Wichita, Kansas, 1990 through 2008

The Equus Beds aquifer in south-central Kansas is a primary water-supply source for the city of Wichita. Water-level declines because of groundwater pumping for municipal and irrigation needs as well as sporadic drought conditions have caused concern about the adequacy of the Equus Beds aquifer as a future water supply for Wichita. In March 2006, the city of Wichita began construction of the Equus Beds Aquifer Storage and Recovery project, a plan to artificially recharge the aquifer with excess water from the Little Arkansas River. Artificial recharge will raise groundwater levels, increase storage volume in the aquifer, and deter or slow down a plume of chloride brine approaching the Wichita well field from the Burrton, Kansas area caused by oil production activities in the 1930s. Another source of high chloride water to the aquifer is the Arkansas River. This study was prepared in cooperation with the city of Wichita as part of the Equus Beds Aquifer Storage and Recovery project. Chloride transport in the Equus Beds aquifer was simulated between the Arkansas and Little Arkansas Rivers near the Wichita well field. Chloride transport was simulated for the Equus Beds aquifer using SEAWAT, a computer program that combines the groundwater-flow model MODFLOW-2000 and the solute-transport model MT3DMS. The chloride-transport model was used to simulate the period from 1990 through 2008 and the effects of five well pumping scenarios and one artificial recharge scenario. The chloride distribution in the aquifer for the beginning of 1990 was interpolated from groundwater samples from around that time, and the chloride concentrations in rivers for the study period were interpolated from surface water samples. Five well-pumping scenarios and one artificial-recharge scenario were assessed for their effects on simulated chloride transport and water levels in and around the Wichita well field. The scenarios were: (1) existing 1990 through 2008 pumping conditions, to serve as a baseline scenario for comparison with the hypothetical scenarios; (2) no pumping in the model area, to demonstrate the chloride movement without the influence of well pumping; (3) double municipal pumping from the Wichita well field with existing irrigation pumping; (4) existing municipal pumping with no irrigation pumping in the model area; (5) double municipal pumping in the Wichita well field and no irrigation pumping in the model area; and (6) increasing artificial recharge to the Phase 1 Artificial Storage and Recovery project sites by 2,300 acre-feet per year. The effects of the hypothetical pumping and artificial recharge scenarios on simulated chloride transport were measured by comparing the rate of movement of the 250-milligrams-per-liter-chloride front for each hypothetical scenario with the baseline scenario at the Arkansas River area near the southern part of the Wichita well field and the Burrton plume area. The scenarios that increased the rate of movement the most compared to the baseline scenario of existing pumping between the Arkansas River and the southern boundary of the well field were those that doubled the city of Wichita’s pumping from the well field (scenarios 3 and 5), increasing the rate of movement by 50 to 150 feet per year, with the highest rate increases in the shallow layer and the lowest rate increases in the deepest layer. The no pumping and no irrigation pumping scenarios (2 and 4) slowed the rate of movement in this area by 150 to 210 feet per year and 40 to 70 feet per year, respectively. In the double Wichita pumping scenario (3), the rate of movement in the shallow layer of the Burrton area decreased by about 50 feet per year. Simulated chloride rate of movement in the deeper layers of the Burrton area was decreased in the no pumping and no irrigation scenarios (2 and 4) by 80 to 120 feet per year and 50 feet per year, respectively, and increased in the scenarios that double Wichita’s pumping (3 and 5) from the well field by zero to 130 feet per year, with the largest increases in the deepest layer. In the increased Phase 1 artificial recharge scenario (6), the rate of chloride movement in the Burrton area increased in the shallow layer by about 30 feet per year, and decreased in the middle and deepest layer by about 10 and 60 feet per year, respectively. Comparisons of the rate of movement of the simulated 250-milligrams-per-liter-chloride front in the hypothetical scenarios to the baseline scenario indicated that, in general, increases to pumping in the well field area increased the rate of simulated chloride movement toward the well field area by as much as 150 feet per year. Reductions in pumping slowed the advance of chloride toward the well field by as much as 210 feet per year, although reductions did not stop the movement of chloride toward the well field, including when pumping rates were eliminated. If pumping is completely discontinued, the rate of chloride movement is about 500 to 600 feet per year in the area between the Arkansas River and the southern part of the Wichita well field, and 70 to 500 feet per year in the area near Burrton with the highest rate of movement in the shallow aquifer layer. The averages of simulated water-levels in index monitoring wells in the Wichita well field at the end of 2008 were calculated for each scenario. Compared to the baseline scenario, the average simulated water level was 5.05 feet higher for the no pumping scenario, 4.72 feet lower for the double Wichita pumping with existing irrigation scenario, 2.49 feet higher for the no irrigation pumping with existing Wichita pumping scenario, 1.53 feet lower for the double Wichita pumping with no irrigation scenario, and 0.48 feet higher for the increased Phase 1 artificial recharge scenario. The groundwater flow was simulated with a preexisting groundwater-flow model, which was not altered to calibrate the solute-transport model to observed chloride-concentration data. Therefore, some areas in the model had poor fit between simulated chloride concentrations and observed chloride concentrations, including the area between Arkansas River and the southern part of the Wichita well field, and the Hollow-Nikkel area about 6 miles north of Burrton. Compared to the interpreted location of the 250-milligrams per liter-chloride front based on data collected in 2011, in the Arkansas River area the simulated 250-milligrams per liter-chloride front moved from the river toward the well field about twice the rate of the actual 250-milligrams per liter-chloride front in the shallow layer and about four times the rate of the actual 250-milligrams per liter-chloride front in the deep layer. Future groundwater-flow and chloride-transport modeling efforts may achieve better agreement between observed and simulated chloride concentrations in these areas by taking the chloride-transport model fit into account when adjusting parameters such as hydraulic conductivity, riverbed conductance, and effective porosity during calibration. Results of the hypothetical scenarios simulated indicate that the Burrton chloride plume will continue moving toward the well field regardless of pumping in the area and that one alternative may be to increase pumping from within the plume area to reverse the groundwater-flow gradients and remove the plume. Additionally, the results of modeling these scenarios indicate that eastward movement of the Burrton plume could be slowed by the additional artificial recharge at the Phase 1 sites and that decreasing pumping along the Arkansas River or increasing water levels could retard the movement of chloride and may prevent further encroachment into the southern part of the well field area.

Kansas↗

A summary of information on the rust Puccinia psidii Winter (guava rust) with emphasis on means to prevent introduction of additional strains to Hawaii

The neotropical rust fungus Puccinia psidii ( P. psidii ) was originally described from the host common guava in its native Brazil but has been found since on hosts throughout the myrtle family (Myrtaceae), including a dramatic host jump to nonnative Eucalyptus plantations. Most rust fungi are able to live only on a very narrow range of host species. P. psidii is unusual both for having a broad host range and for the intensity of its damage to susceptible young growth. This rust first got a foothold in the United States in Florida more than three decades ago. The U.S. Department of Agriculture (USDA) has since considered it a nonactionable, nonreportable pest. Hawaii and Florida are the only two states with native species in the myrtle family. Over a period of 30 years, this rust has done little damage to any of the scattered native Myrtaceae in Florida, although the host range of the rust has gradually grown to about 30 mostly nonnative species in the family, apparently because of increasing genetic variety of the rust by repeated introductions. However, Florida&rsquo;s native Myrtaceae are among the roughly 1,100 neotropical species that are largely resistant to P. psidii . The 3,000 species of non-neotropical Myrtaceae of the Pacific, Australia, Asia, and Africa are expected to prove much more vulnerable to P. psidii . Little is known about the genetics or genetic strains of P. psidii , although existing literature shows that there are numerous strains that have differential ability to infect suites of host plants. The rust was first recorded in the state of Hawaii on Oahu in April 2005 and quickly spread throughout the Hawaiian Islands. The main concern in Hawaii became the potential threat to ohia, Metrosideros polymorpha (Myrtaceae), the endemic forest tree species overwhelmingly important in Hawaii&rsquo;s nature and culture. The potential ecological consequences of a virulent strain of rust on ohia forests are immense, due to its role as a foundation tree species and the diversity of niches it fills in Hawaii. A single genetic strain of the rust is established in Hawaii, apparently composed of a single genotype lacking sexual reproduction. P. psidii has been found statewide in Hawaii attacking Myrtaceae from near sea level to about 1,200 m elevation in areas with rainfall ranging from 750&ndash;5,000 mm. Five of eight native Myrtaceae and at least 15 nonnative species have been observed as hosts of P. psidii in Hawaii. The federally endangered Eugenia koolauensis (nioi) and the nonendangered indigenous species Eugenia reinwardtiana are severely damaged. The introduced (an Asian species) and invasive rose apple, Syzygium jambos , is severely affected at a landscape scale, with widespread crown dieback and many instances of complete tree death. In spite of billions of wind-dispersed rust spores produced from rose apple infestations during 2006 to 2008, adjacent ohia have been little affected to date by the rust strain in Hawaii. Within the elevation range of the rust, P. psidii is found on less than 5 percent of the ohia trees in the wild; on those ohia trees on which the rust is found, it is normally found on less than 5 percent of the leaves. The strain in Hawaii has not attacked many of the species known to be infected by the rust elsewhere, including common guava. On the basis of the very substantial genetic diversity of the much-studied, crop-damaging species of the genus Puccinia , there is good reason to believe that there are at minimum dozens and likely hundreds or thousands of genotypes of P. psidii , likely concentrated in the core range in Brazil but with potential for dispersal by globalization. Multiple genotypes are believed already present in the United States and certain to spread freely in the absence of restrictions. The U.S. Forest Service has initiated a major collaborative project in Brazil to investigate the genetics of susceptibility of Hawaii&rsquo;s ohia to P. psidii , but initial results will likely not be available for several years. If just one more strain reaches Hawaii, the consequences could be dire for ohia, with each new genotype arriving having an unknown likelihood of increasing damage to ohia; possibilities for mutation and (or) genetic mixing, even with asexual strains, are apparently substantial, based on what is known about other Puccinia species. Investigations are needed to clarify rust-nioi relationships. However, it is likely that keeping out new strains of P. psidii may be important for long-term survival of nioi as well as for the health of ohia forest. The source of Hawaii&rsquo;s initial invasion by P. psidii is uncertain but is strongly suspected to have been decorative foliage of species in the myrtle family from the mainland United States, most likely California, where there had been outbreaks of this rust on cultivated myrtle in 2005. In 2006&ndash;7, Maui&rsquo;s Hawaii Department of Agriculture (HDOA) inspectors intercepted several P. psidii infected shipments of foliage myrtle, shipped from several California counties. Recognizing the huge threat of the rust to Hawaii&rsquo;s one million acres of ohia forests, and consequently to Hawaii&rsquo;s watersheds and biodiversity, Hawaii&rsquo;s Board of Agriculture unanimously approved an interim rule in August 2007 banning importation of plants in the myrtle family from &ldquo;infested areas,&rdquo; specified as South America, Florida, and California. However, the interim rule has not been made permanent by HDOA, and the department has stated that it needs further information to formulate a long-term rule that imposes appropriate measures. Rust spores can survive for 2 to 3 months, and the pathogen can be transported to Hawaii on Myrtaceae from anywhere in the world through the United States mainland. There is much geographic reshuffling of flowers and foliage among the far-flung firms in the trade, especially for bouquet making. Because P. psidii is a nonactionable and nonreportable pest in the United States, foliage and flowers of the myrtle family can move freely into the country (usually but not necessarily always through the ports of Miami or Los Angeles), and from state to state. Currently, the State of Hawaii regulates incoming plant material in the family Myrtaceae by visual inspection. Inspection capacity and latent (asymptomatic) infections limit the ability to detect the rust. New molecular tests could improve detection efficiency, but the cost and the time required to process samples currently precludes their routine use in ports of entry. Interdiction, which has effectively kept coffee rust ( Hemileia vastatrix ) out of Hawaii for 120 years, offers the strongest protection for Hawaii&rsquo;s native ecosystems from P. psidii . Interdiction of Myrtaceae from the continental United States could have the important supplementary benefit of preventing establishment in Hawaii of other very significant pests of multiple species of Myrtaceae that are already in the country, including: the Eugenia psyllid Trioza eugeniae (Hemiptera: Psyllidae); Chrysophtharta m-fuscum , the Eucalyptus tortoise beetle (Coleoptera: Chrysomelidae); Leptocybe invasa , the blue gum chalcid wasp (Hymenoptera: Chalcidae); and the fungal pathogens Mycosphaerella molleriana (Ascomycota: Mycosphaerelliaceae, crinkle leaf disease of Eucalyptus spp.) and Neofusicoccum parvum (Ascomycota: Botryosphaeriaceae), currently causing serious damage to Syzygium paniculatum in south Florida nurseries. Each of these pests would be likely to cause very significant damage to native and (or) cultivated Myrtaceae in Hawaii. Each of these pests is a prime candidate for transport by the foliage and (or) nursery stock pathways from Florida and California into Hawaii. Hawaii Department of Agriculture has a clear mandate to protect Hawaii&rsquo;s natural environment, forestry and cultivated Myrtaceae. Principles of the World Trade Organization&rsquo;s Treaty on Sanitary and Phytosanitary Measures and the International Plant Protection Convention are consistent with the right of Hawaii to take action. The current threat of P. psidii and the other five serious threats to Myrtaceae are primarily posed by the importation of infected plants from the continental United States; however, that may change in the future. If Hawaii were to decide to take a stand (through State regulation) to protect its native and introduced Myrtaceae, there is a possibility that USDA would consider Federal regulation of Myrtaceae from foreign countries.

Hawaii↗

Hydrology and water quality of Shell Lake, Washburn County, Wisconsin, with special emphasis on the effects of diversion and changes in water level on the water quality of a shallow terminal lake

Shell Lake is a relatively shallow terminal lake (tributaries but no outlets) in northwestern Wisconsin that has experienced approximately 10 feet (ft) of water-level fluctuation over more than 70 years of record and extensive flooding of nearshore areas starting in the early 2000s. The City of Shell Lake (City) received a permit from the Wisconsin Department of Natural Resources in 2002 to divert water from the lake to a nearby river in order to lower water levels and reduce flooding. Previous studies suggested that water-level fluctuations were driven by long-term cycles in precipitation, evaporation, and runoff, although questions about the lake&rsquo;s connection with the groundwater system remained. The permit required that the City evaluate assumptions about lake/groundwater interactions made in previous studies and evaluate the effects of the water diversion on water levels in Shell Lake and other nearby lakes. Therefore, a cooperative study between the City and U.S. Geological Survey (USGS) was initiated to improve the understanding of the hydrogeology of the area and evaluate potential effects of the diversion on water levels in Shell Lake, the surrounding groundwater system, and nearby lakes. Concerns over deteriorating water quality in the lake, possibly associated with changes in water level, prompted an additional cooperative project between the City and the USGS to evaluate efeffects of changes in nutrient loading associated with changes in water levels on the water quality of Shell Lake. Numerical models were used to evaluate how the hydrology and water quality responded to diversion of water from the lake and historical changes in the watershed. The groundwater-flow model MODFLOW was used to simulate groundwater movement in the area around Shell Lake, including groundwater/surface-water interactions. Simulated results from the MODFLOW model indicate that groundwater flows generally northward in the area around Shell Lake, with flow locally converging toward the lake. Total groundwater inflow to Shell Lake is small (approximately 5 percent of the water budget) compared with water entering the lake from precipitation (83 percent) and surface-water runoff (13 percent). The MODFLOW model also was used to simulate average annual hydrologic conditions from 1949 to 2009, including effects of the removal of 3 billion gallons of water during 2003&ndash;5. The maximum decline in simulated average annual water levels for Shell Lake due to the diversion alone was 3.3 ft at the end of the diversion process in 2005. Model simulations also indicate that although water level continued to decline through 2009 in response to local weather patterns (local drought), the effects of the diversion decreased after the diversion ceased; that is, after 4 years of recovery (2006&ndash;9), drawdown attributable to the diversion alone decreased by about 0.6 ft because of increased groundwater inflow and decreased lake-water outflow to groundwater caused by the artificially lower lake level. A delayed response in drawdown of less than 0.5 ft was transmitted through the groundwater-flow system to upgradient lakes. This relatively small effect on upgradient lakes is attributed in part to extensive layers of shallow clay that limit lake/groundwater interaction in the area. Data collected in the lake indicated that Shell Lake is polymictic (characterized by frequent deep mixing) and that its productivity is limited by the amount of phosphorus in the lake. The lake was typically classified as oligotrophic-mesotrophic in June, mesotrophic in July, and mesotrophic-eutrophic in August. In polymictic lakes like Shell Lake, phosphorus released from the sediments is not trapped near the bottom of the lake but is intermittently released to the shallow water, resulting in deteriorating water quality as summer progresses. Because the productivity of Shell Lake is limited by phosphorus, the sources of phosphorus to the lake were quantified, and the response in water quality to changes in phosphorus inputs were evaluated by means of eutrophication models. During 2009, the total input of phosphorus to Shell Lake was 1,730 pounds (lb), of which 1,320 lb came from external sources (76 percent) and 414 lb came from internal loading from sediments in the lake (24 percent). The largest external source was from surface-water runoff, which delivered about 52 percent of the total phosphorus load compared with about 13 percent of the water input. The second largest source was from precipitation (wetfall and dryfall), which delivered 19 percent of the load compared to about 83 percent of the water input. Contributions from septic systems and groundwater accounted for about 3 and 2 percent, respectively. Increased runoff raises water levels in the lake but does not necessarily increase phosphorus loading because phosphorus concentrations in the tributaries decline during increased flow, possibly because of shorter retention times in upstream wetlands. Phosphorus loading to the lake in 2009 represented what occurred after a series of dry years; therefore, this information was combined with data from 2011, a wet year, to estimate phosphorus loading during a range of hydrologic conditions by estimating loading from each component of the phosphorus budget for each year from 1949 to 2011. Comparisons of historical water-quality records with historical water levels and applications of a hydrodynamic model (Dynamic Lake Model, DLM) and empirical eutrophication models were used to understand how changes in water level and the coinciding changes in phosphorus loading affect the water quality of Shell Lake. DLM simulations indicate that large changes in water level (approximately 10 ft) affect the persistence of stratification in the lake. During periods with low water levels, the lake is a well-mixed, polymictic system, with water quality degrading slightly as summer progresses. During periods with high water levels, the lake is more stratified, and phosphorus from internal loading is trapped in the hypolimnion and released later in summer, which results in more extreme seasonality in water quality and better clarity in early summer. Results of eutrophication model simulations using a range in external phosphorus inputs illustrate how water quality in Shell Lake (phosphorus and chlorophyll a concentrations and Secchi depths) responds to changes in external phosphorus loading. Results indicate that a 50-percent reduction in external loading from that measured in 2009 would be required to change phosphorus concentrations from 0.018 milligram per liter (mg/L) (measured in 2009) to 0.012 mg/L (estimated for the mid-1800s from analysis of diatoms in sediment cores). Such reductions in phosphorus loading cannot be accomplished by targeting septic systems or internal loading alone because septic systems contribute only about 3 percent of the phosphorus input to the lake, and internal loading from the sediments of Shell Lake contributes only about 25 percent of phosphorus input. Complete elimination of phosphorus from septic systems and internal loading would decrease the phosphorus concentrations in the lake by 0.003&ndash;0.004 mg/L. Therefore, reducing phosphorus concentration in the lake more than by 0.004 mg/L requires decreasing phosphorus loading from surface-water contributions, primarily runoff to the lake. Reconstructed changes in water quality from 1860 to 2010, based on changes in the diatom communities archived in the sediments and eutrophication model simulations, suggest that anthropogenic changes in the watershed (sawmill construction in 1881; the establishment of the village of Shell Lake; and land-use changes in the 1920s, including increased agriculture) had a much larger effect on water quality than the natural changes associated with fluctuations in water level. Although the effects of natural changes in water level on water quality appear to be small, changes in water level do have a modest effect on water quality, primarily manifested as small improvements during higher water levels. Fluctuations in water level, however, have a larger effect on the seasonality of water-quality patterns, with better water quality, especially increased Secchi depths, in early summer during years with high water levels.

Wisconsin↗

Evaluation of larval lamprey survival following salvage: A pilot study

Larval lampreys ( Entosphenus tridentatus and Lampetra spp.) are vulnerable to anthropogenic water-level fluctuations that can dewater their habitat. Dewatering events occur regularly in the Columbia River Basin for operation and management of hydropower facilities, seasonal or maintenance closures of irrigation diversions, and in-water construction projects, including for habitat restoration. Salvage efforts which can be initiated before, during, and after dewatering events are resource-intensive and are conducted based on the assumption that salvage will reduce lamprey mortality. This pilot study was the first formal assessment of the efficacy of salvage efforts, evaluating the survival and performance of larval lamprey following various salvage techniques. Lampreys were salvaged during dewatering events at three field sites under variable environmental conditions (summer and fall of 2020) and then held in the laboratory for 60 days to monitor survival, growth, and burrowing performance. Four salvage treatments were defined to represent combinations of typical salvage techniques and stressors, including multiple passes of standard electrofishing (SEF), lamprey-specific electrofishing (LEF), and modified lamprey-specific electrofishing (MLEF; probes in direct contact with dewatered, but moist substrate) as well as extended exposure on the surface and walking on sediment where lampreys were burrowed. Control groups did not experience dewatering and were collected using LEF in areas away from treatment groups. Treatments were designed to increase in intensity, from treatment 1 (walking and exposure) to treatment 4 (multiple passes of SEF, LEF and MLEF). Study sites included an earthen hatchery rearing pond (North Toutle Hatchery) dewatered in July, and two irrigation diversions (Wapato and Sunnyside diversions on Yakima River) dewatered at the end of the irrigation season in October. Treatments were executed inside circular 1 m2 enclosures that were randomly positioned in habitats expected to be dewatered. A solid, weighted ring at the bottom of the enclosure penetrated the sediment and netting extended through the water column to a floating upper ring. We deployed eight enclosures per treatment at each test site, executed the four salvage treatments, collected lamprey from within each enclosure and transported them to the laboratory, along with the control groups, for the 60-day holding period. Burrowing performance was tested in sand 1 day after the field effort and in field-collected sediment 30 days after the field effort. Mortality was documented and lamprey were measured at 1, 30, and 60 days in the laboratory and fish weights were used to calculate standard growth rate (SGR) for each site and treatment group. We collected 328 larval lampreys at our three test sites, including 71 controls and 257 larvae exposed to dewatering and salvage treatments. Overall mortality for the 60-day laboratory holding period was 11.9%. Most mortality occurred within 1-day after treatment (51.3%) and there was limited mortality past 30 days (2.6%). At the North Toutle Hatchery, we observed substantial mortality during the field tests in July, both inside and outside of our test enclosures. Mortality within our test enclosures ranged from 96.7 to 98.8% for treatment 1, 45.9 to 52.2% for treatment 3 and 6.7 to 7.1% for treatment 4. The elevated mortality at this site and logistical challenges with the execution of treatments 1 and 2 resulted in few fish (5 total for treatment 1) or no fish (treatment 2) available for testing in the laboratory. Only one larval lamprey died during field tests at the Wapato and Sunnyside irrigation diversions during testing in October. The single mortality was in treatment 1 (11.1%) and no mortalities were observed outside of the test enclosures. We used logistic regression to estimate survival of larval lampreys transported to the laboratory and held for 24 h. The Wapato and Sunnyside field sites were pooled for logistic regression and the North Toutle Hatchery site was analyzed separately due to dramatically different environmental conditions. We found that treatment 1 reduced larval survival more than any other treatment during both the summer and fall dewatering events. Trends among survival for treatments 2-4 were less clear. The unique stressor included in the first treatment, but not in other treatments, was a 2-hour exposure period during which larvae were left lying on the surface of the sediment. Treatment 1 also experienced a walking action (foot pressure on the surface of the exposed sediment). The walking action was also included in treatment 4, both before and after dewatering, along with multiple passes of various electrofishing techniques, as this treatment was designed to be a worst-case scenario for lamprey salvage. Despite what appeared to be significant stressors associated with treatment 4, the logistic regression for survival up to 24 hours in the laboratory showed that the odds of surviving treatment 4 were 16 times higher than the odds of surviving treatment 1 at Wapato and Sunnyside (combined). The same comparison at the North Toutle Hatchery showed the odds were 226 times higher for lamprey to survive treatment 4 compared to treatment 1. Lamprey from all study sites initiated burrowing activity with median times less than 10.5 seconds in both sand (day 1) and field-collected sediment (day 30). The fastest burrowing start times were less than 1.0 second and the slowest was 3.2 minutes. Lamprey behavioral responses during burrowing ability tests were variable. Some lampreys immediately moved from the release location near the surface of the water toward the sediment and began burrowing while others swam around the aquarium near the surface of the water before exploring the sediment to select a burrowing location. The median time to complete burrowing for all treatment groups and sample periods ranged from 9.9 to 48.1 seconds. No significant differences in SGR were detected between treatment and control groups at any test site. Laboratory water temperatures for the North Toutle Hatchery study groups were maintained at 15°C, giving lamprey a growth advantage compared to the Wapato and Sunnyside groups which were maintained at 10℃. SGR for lamprey collected at the North Toutle Hatchery ranged from 0.83% weight gain/day for controls to 2.04%/day for treatment 3. SGR at Wapato ranged from 0.27 to 0.67%/day and from 0.60 to 0.90 %/day at Sunnyside. Overall, SGR was consistently lower at every site for the controls compared to any of the treatment groups, although none of the differences were significant. The variability at some sites in initial lamprey size, combined with inherent variability in growth rates, limited our ability to make conclusions about how different salvage treatments influenced SGR. Treatment 1 stood out among the salvage treatments at all study sites. In this treatment, lampreys exposed on the surface of the sediment, awaiting salvage, were vulnerable to reduced survival, even under mild environmental conditions. The risk of mortality was greatest for the summer dewatering event at the North Toutle Hatchery. The remaining treatments, even with multiple passes of various electrofishing techniques, did not generally have large negative impacts on lamprey during our tests. Lamprey survival rates for these treatments were relatively high, especially at the fall dewatering sites when environmental conditions were mild. Thus, salvage efforts, despite being resource intensive, likely have limited negative outcomes for larval lamprey and make substantial contributions to lamprey conservation efforts.

Report↗

Effects of the earthquake of March 27, 1964, at Seward, Alaska

Seward, in south-central Alaska, was one of the towns most devastated by the Alaska earthquake of March 27, 1964. The greater part of Seward is built on an alluvial fan-delta near the head of Resurrection Bay on the southeast coast of the Kenai Peninsula. It is one of the few ports in south-central Alaska that is ice free all year, and the town’s economy is almost entirely dependent upon its port facilities. The Alaska earthquake of March 27, 1964, magnitude approximately 8.3–8.4, began at 6:36 p.m. Its epicenter was in the northern part of the Prince William Sound area; focal depth was 20–50 km. Strong ground motion at Seward lasted 3–4 minutes. During the shaking, a strip of land 50–400 feet wide along the Seward waterfront, together with docks and other harbor facilities, slid into Resurrection Bay as a result of large-scale submarine landsliding. Fractures ruptured the ground for'severa1 hundred feet back from the landslide scarps. Additional ground was fractured in the Forest Acres subdivision and on the alluvial floor of the Resurrection River valley; fountaining and sand boils accompanied the ground fracturing. Slide-generated wares, possibly seiche waves, and seismic sea waves crashed onto shore; ware runup was as much as 30 feet above mean lower low water and caused tremendous damage; fire from burning oil tanks added to the destruction. Damage from strong ground motion itself was comparatively minor. Tectonic subsidence of about 3.6 feet resulted in low areas being inundated at high tide. Thirteen people were killed and five were injured as a result of the earthquake. Eighty-six houses were totally destroyed and 260 were heavily damaged. The harbor facilities were almost completely destroyed, and the entire economic base of the town was wiped out. The total cost to replace the destroyed public and private facilities was estimated at $22 million. Seward lies on the axis of the Chugach Mountains geosyncline. The main structural trend in the mapped area, where the rocks consist almost entirely of graywacke and phyllite, is from near north to N. 20° E. Beds and cleavage of the rocks commonly dip 70° W. or NW. to near vertical. Locally, the rocks are complexly folded or contorted. So major faults were found in the mapped area, but small faults, shear zones, and joints are common. Surficial deposits of the area hare been divided for mapping into the following units: drift deposits, alluvial fan deposits, valley alluvium, intertidal deposits, landslide deposits, and artificial fill. Most of these units intergrade and were deposited more or less contemporaneously. The drift deposits consist chiefly of till that forms moraines along the lower flanks of the Resurrection River valley and up tributary valleys. The till is predominantly silt and sand and lesser amounts of clay-size particles, gravel, cobbles, and boulders. Glacial outwash and stratified ice-contact deposits constitute the remainder of the drift deposits. Fans and fan-deltas have been deposited at the valley mouths of tributary streams. Some, including the one upon which Seward built, project into Resurrection Bay, and deltaic-type deposits form their distal edges. The larger fans—composed chiefly of loosely compacted and poorly sorted silt, sand, and gravel—form broad aprons having low gradients. The fan deposits range in thickness from about 100 feet to possibly several hundred feet and, at least in some places, lie on a platform of compact drift. Smaller fans at the mouths of several canyons have steep gradients and considerable local relief. Valley alluvium, deposited chiefly by the Resurrection River, consists mostly of coarse sand and fine to medium gravel. In the axial part of the valley it is probably more than 100 feet thick. Near the head of Resurrection Bay, the alluvium is underlain by at least 75 feet of marine deltaic sediments, which are in turn underlain by 600 or more feet of drift in the deepest part of the bedrock valley. Beach, deltaic, and estuarine sediments, deposited on intertidal flats at the head of the bay and along far1 margins that extend into the bay, arc mapped as intertidal deposits. They consist mostly of silt, sand, and fine gravel, and lesser amounts of clay-size particles. The earthquake reactivated old slides and trigged new ones in the mountains. Rock and snow avalanches, debris flows, and creep of talus deposits characterized slide activity on the steeper slops. The Seward waterfront had been extended before the earthquake by adding artificial fill consisting of loose sand and gravel; part of the lagoon area had been filled with refuse. After the earth- quake, fill, consisting of silt and sand dredged from the head of the bay, was pumped onto part of the lagoon area and also on land at the northwest corner of the bay. Response to the disaster was immediate and decisive. City, State, and Federal agencies, as well as other organizations and individuals, gave unstintingly of their time and facilities. Within a few days, there was temporary restoration of water, sewerage, and electrical facilities. The U.S. Army Corps of Engineers was authorized to select sites and construct a new dock for the Alaska Railroad, a new small-boat basin, and related facilities. The firm of Shannon and Wilson, Inc., under contract to the Corps of Engineers, investigated subsurface soils extensively to determine the factors responsible for the sliding along the Seward waterfront and to assist in site selection for reconstruction of the destroyed harbor facilities. Borings also made along the Seward waterfront and at the head of the bay, and laboratory tests were conducted on pertinent samples. These studies were augmented by geophysical studies both on land and in the bay. In addition, the Corps of Engineers made shallow borings on the intertidal flats at the head of the bay and performed pile-driving and load tests. Borings also were drilled and test pits were dug in the subdivision of Forest Acres. Sliding along the Seward waterfront markedly deepened the water along the former shoreline. Post-earthquake slopes of the bay floor immediately offshore also are steeper in places than before the earthquake. The strong ground motion of the earthquake triggered the landsliding, but several factors may have contributed to the magnitude and characteristics of the slides. These factors are: (1) the long duration of strong ground motion, (2) the grain size and texture of the material involved in the sliding, (3) the probability that the finer grained materials liquefied and flowed seaward, and (4) the added load of manmade facilities built on the edge of the shore, Secondary effects of the slides themselves—sudden drawdown of water, followed by the weight of returning waves—also may have contributed to the destruction. Submarine sliding at the northwest corner of the bay occurred in fine-grained deltaic deposits whose frontal slopes probably were in metastable equilibrium under static conditions. Uplift pressures from aquifers under hydrostatic head, combined with the probable liquefaction characteristics of the sediments when vibrated by strong ground motion, probably caused the material to slide and flow seaward as a heavy slurry. Under static conditions, no major shoreline or submarine landsliding is expected in the Seward area; in the event of another severe earthquake, however, additional sliding is likely along the Seward waterfront and also in the deltaic deposits at the northwest corner of the bay. Fractured ground in back of the present shoreline along the Seward waterfront is an area of incipient landslides that would be unstable under strong shaking. For this reason the Scientific and Engineering Task Force placed the area in a high-risk classification and recommended no repair, rehabilitation, or new construction in this area involving use of Federal funds; it was further recommended that the area should be reserved for park or other uses that do not involve large congregations of people. The deltaic deposits at the head of the bay probably also would be susceptible to sliding during another large earthquake. This sliding would result in further landward retreat of the present shoreline toward the new railroad dock. Specifications for the new dock, whose seaward end is now approximately 1,100 feet from the back scarp of the subaqueous landslide, require design pro- visions to withstand seismic shock up to certain limits. Earthquake-induced fracturing of the ground in the subdivision of Forest Acres was confined to the lower part of a broad alluvial fan. There, sewer and water lines were ruptured and the foundations of some homes were heavily damaged. Landsliding, such as occurred along the shoreline of the bay, was not a contributing cause of the fracturing. Two hypotheses are offered to explain the fracturing: 1. Seismic energy was transformed into visible surface waves of such amplitude that the strength of surface layer was exceeded and rupturing occurred; tensional and compressional stresses alternately opened and closed the fractures and forced out water and mud. 2. Compaction by vibration of the fine-grained deposits of the fan caused ground settlement and fracturing; ground water under temporary hydrostatic head was forced to the surface as fountains and carried the finer material with it. Water waves that crashed onto shore, while shaking was still continuing, were generated chiefly by onshore and offshore landsliding. Waves that overran the shores about 25 minutes after shaking stopped and that continued to arrive for the next several hours are believed to be seismic sea waves (tsunamis) that originated in an uplifted area in the Gulf of Alaska. During the time of seismic sea-wave activity and perhaps preceding it, seiche wares also may have been generated within Resurrection Bay and complicated the wave effects along the shoreline.

Alaska↗

Distribution and demography of Coastal Cactus Wrens in Southern California, 2015–19

Surveys and monitoring for the coastal Cactus Wren ( Campylorhynchus brunneicapillus ) were completed in San Diego County between March 2015 and July 2019. A total of 383 plots were surveyed across 3 genetic clusters (Otay, Lake Jennings, and Sweetwater/Encanto). From 2015 to 2019, 317 plots were surveyed 8 times (twice per year in 2015, 2017–19). Additional plots were added in later years as wrens were discovered in new locations. We found differences in the proportion of plots occupied in the genetic clusters, with a lower proportion of plots occupied in the Otay cluster than in the Lake Jennings and Sweetwater/Encanto clusters in all years. Plot occupancy increased each year in the Otay and Sweetwater/Encanto clusters but not in the Lake Jennings cluster. The number of Cactus Wren territories increased from 2015 through 2018, and then decreased in 2019 in all three genetic clusters. We monitored nesting activities for two populations of Cactus Wrens in southern San Diego County. The Otay population consisted of two sites within the Otay genetic cluster, and the San Diego population consisted of two sites within the Sweetwater/Encanto and Lake Jennings genetic clusters. Nest monitoring occurred at 10–13 territories per year in the Otay population and 14–18 territories in the San Diego population from 2015 through 2019. All territories were occupied by pairs except two territories in 2015, five in 2016, and two in 2019. Between 46 and 74 Cactus Wren nests were monitored each year, which totaled 295 monitored nests from 2015 to 2019. To evaluate the direct influence of precipitation on breeding success, bio-year precipitation (“precipitation”) was calculated from July 1 of the prior year through June 30 of the breeding season year. Overall apparent nest success was positively influenced by precipitation with the lowest apparent nest success of 50 percent in 2015 and the highest apparent nest success of 72 percent in 2017, corresponding to the second lowest and the highest precipitation years, respectively. Apparent nest success also was higher in the Otay population than in the San Diego population. The number of brood nests initiated per pair and the number of renesting attempts per pair also were higher in years with more precipitation. Other metrics of Cactus Wren nesting success and productivity were positively influenced by the amount of precipitation, including clutch size and egg hatching success. The percent of hatchlings that fledged was greater in the Otay population than in the San Diego population but was not influenced by precipitation. The number of fledglings per pair was higher in years with more precipitation and was greater in the Otay population than in the San Diego population. Predation was the predominant cause of nest failure in both populations. Analysis of Cactus Wren daily nest survival rate indicated that there was a population, and possibly a precipitation effect on nest survival, with the daily survival rate for the Otay population significantly higher than for the San Diego population and weak increase in the daily survival rate with more precipitation. A total of 629 Cactus Wrens were banded during the course of the study, 360 in the San Diego population and 269 in the Otay population. Between 2015 and 2019, we resighted 301 color-banded adult birds that ranged between 1 and 8 years old. One additional color-banded bird was resighted in San Pasqual Valley (as part of a separate study); this bird originated in the San Diego population and was excluded from our analyses. Annual survival was higher for adult Cactus Wrens (ranging from 60 to 70 percent) than for first-year wrens (ranging from 20 to 28 percent) and varied by year. Annual survival was also weakly but positively correlated with precipitation. Annual survival was higher for first year and adult Cactus Wrens following years with increased precipitation. We found no evidence that survival differed by population. Banding also allowed us to examine whether there were differences in movement of adult and first-year Cactus Wrens by year or by population. We found that average dispersal distance for first-year Cactus Wrens was 1.9 kilometers in the Otay population and 1.6 kilometers in the San Diego population and did not differ by population or year. Dispersal between populations was not common. We detected five instances of movement of first-year wrens between the San Diego and Otay populations. All movements into and out of the San Diego population were from or into territories in the Sweetwater area. We detected no movement between the Lake Jennings site and either of the Sweetwater or Otay sites; however, we did detect one wren that dispersed from Lake Jennings to the San Pasqual Valley population in 2019, which was a distance of 26.4 kilometers. Adult Cactus Wrens were site-faithful, with 87 percent of adults remaining on the same territory between breeding seasons. Precipitation may be a weak driver of movement for adult Cactus Wrens, with adults more likely to remain on the same territory following years of increased precipitation. There was no difference in adult movement between populations. Arthropods were collected in pitfall traps and by vacuum in 23 Cactus Wren territories during 3 sampling periods in 2016 (early nesting, peak nesting, and late nesting). Arthropods of 19 orders and at least 128 families were collected. Analysis of 43 Cactus Wren fecal samples identified 10 arthropod orders that were present in more than 10 percent of fecal samples. The most abundant arthropod order collected was Hymenoptera; however, Cactus Wrens consumed arthropods in the order Hymenoptera significantly less than their availability, suggesting that this order was avoided. No other orders were significantly selected or avoided; however, selection indices of arthropod families identified that two families of arthropods (Isopoda Porcellionidae [woodlice] and Hymenoptera Formicidae [ants]) were avoided. After excluding the taxa that were avoided or not represented in fecal samples, 95 percent of Cactus Wren prey items were collected in pitfall traps and 5 percent were collected by vacuum. The most abundant prey orders captured were Diptera, Coleoptera, Hemiptera, Hymenoptera, and Aranea. Analysis of the abundance of Cactus Wren prey items by vegetation type and sampling period indicated that vegetation type by itself was not a significant predictor of arthropod abundance but interacted with sampling period. Seasonal availability of arthropods was highest in the peak nesting period, followed by early and late nesting periods for California sagebrush ( Artemisia californica ), lemonadeberry ( Rhus integrifolia ), non-native grass, and bare ground, whereas availability increased from early to late nesting periods for blue elderberry ( Sambucus mexicana spp. caerulea ), cactus ( Opuntia spp. and Cylindropuntia spp.), California buckwheat ( Eriogonum fasciculatum ), native bunch grasses, and black mustard ( Brassica nigra ). During the early nesting period, arthropods were most abundant in native bunch grasses and least abundant in lemonadeberry. During the peak nesting period, arthropods were most abundant in native bunch grasses and in areas of bare ground and were least abundant in cactus and blue elderberry. During late nesting, arthropods were most abundant in blue elderberry and non-native grass and least abundant in lemonadeberry and mustard. Each year from 2015 to 2019, vegetation data were collected at the same 23 territories where arthropods were sampled: 9 territories in the Otay population and 14 territories in the San Diego population. Cactus, California buckwheat, and non-native grasses were detected within at least 60 percent of sampling points in the Otay population. Cactus, California sagebrush, California buckwheat, non-native grass, and black mustard each were detected within an average of 40 percent of sampling points in the San Diego population. No native bunch grass or lemonadeberry were recorded at the Lake Jennings site within the San Diego population. The cover of shrub species was relatively stable throughout the 5 years. Cover of herbaceous species and bare ground had greater annual variation than shrub species. We found that vegetation cover varied widely among territories, with territory accounting for 69 percent of the variation in vegetation cover. Redundancy analysis allowed us to identify the vegetation types that accounted for the most variation. We used the top scores from the redundancy analysis to identify six vegetation types to be used in generalized linear mixed models analyzing the relationships between vegetation type, precipitation, and Cactus Wren breeding productivity. Three vegetation variables influenced the number of fledglings produced per pair. California sagebrush had a positive effect on the number of fledglings per pair whereas non-native grass and black mustard had a negative effect. Breeding productivity, survival, and movements of adult and first-year Cactus Wrens indicated that the Otay population behaved similarly to, if not out-performed, the San Diego population during the span of our project, suggesting that the driving forces behind low numbers of Cactus Wrens in the Otay population before 2015 were no longer in effect. The Cactus Wren populations in Otay and San Diego reached a peak in 2018, which followed a year of high productivity and survivorship, both of which were correlated with high precipitation. This peak in population size was consistent with reproductive timing and productivity in other bird populations in semi-arid ecosystems that were linked to precipitation and arthropod abundance. We did not find a strong link among arthropod abundance, vegetation composition, and Cactus Wren breeding productivity, likely in part because arthropod abundance varied by vegetation type and sampling period, suggesting that different vegetation types provided important sources of prey at different periods of the breeding season. Arthropod abundance also may not represent arthropod availability when vegetation structure discourages the ground foraging behavior of species such as Cactus Wrens. Cover of non-native grass negatively influenced breeding productivity, although arthropods were abundant in non-native grass. Other factors that could have influenced differential breeding productivity between the Otay and San Diego populations were habitat restoration, control of annual herbaceous vegetation, human disturbance, lingering effects of wildfire, and nest predation. Overall, precipitation appeared to be a driver of Cactus Wren breeding productivity and possibly survival, potentially obscuring proximate effects of arthropod or vegetation composition.

California↗