Geology ReportsSearch

SEARCH · Geology Reports

Results for “Physical Review D”

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

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

At least 73 records · Page 4Linked to original sources

Resource subsidies between stream and terrestrial ecosystems under global change

Streams and adjacent terrestrial ecosystems are characterized by permeable boundaries that are crossed by resource subsidies. Although the importance of these subsidies for riverine ecosystems is increasingly recognized, little is known about how they may be influenced by global environmental change. Drawing from available evidence, in this review we propose a conceptual framework to evaluate the effects of global change on the quality and spatiotemporal dynamics of stream–terrestrial subsidies. We illustrate how changes to hydrological and temperature regimes, atmospheric CO 2 concentration, land use and the distribution of nonindigenous species can influence subsidy fluxes by affecting the biology and ecology of donor and recipient systems and the physical characteristics of stream–riparian boundaries. Climate-driven changes in the physiology and phenology of organisms with complex life cycles will influence their development time, body size and emergence patterns, with consequences for adjacent terrestrial consumers. Also, novel species interactions can modify subsidy dynamics via complex bottom-up and top-down effects. Given the seasonality and pulsed nature of subsidies, alterations of the temporal and spatial synchrony of resource availability to consumers across ecosystems are likely to result in ecological mismatches that can scale up from individual responses, to communities, to ecosystems. Similarly, altered hydrology, temperature, CO 2 concentration and land use will modify the recruitment and quality of riparian vegetation, the timing of leaf abscission and the establishment of invasive riparian species. Along with morphological changes to stream–terrestrial boundaries, these will alter the use and fluxes of allochthonous subsidies associated with stream ecosystems. Future research should aim to understand how subsidy dynamics will be affected by key drivers of global change, including agricultural intensification, increasing water use and biotic homogenization. Our conceptual framework based on the match–mismatch between donor and recipient organisms may facilitate understanding of the multiple effects of global change and aid in the development of future research questions.

Global Change Biology

Diverse portfolios: Investing in tributaries for restoration of large river fishes in the Anthropocene

Rehabilitation of large Anthropocene rivers requires engagement of diverse stakeholders across a broad range of sociopolitical boundaries. Competing objectives often constrain options for ecological restoration of large rivers whereas fewer competing objectives may exist in a subset of tributaries. Further, tributaries contribute toward building a “portfolio” of river ecosystem assets through physical and biological processes that may present opportunities to enhance the resilience of large river fishes. Our goal is to review roles of tributaries in enhancing mainstem large river fish populations. We present case histories from two greatly altered and distinct large-river tributary systems that highlight how tributaries contribute four portfolio assets to support large-river fish populations: 1) habitat diversity, 2) connectivity, 3) ecological asynchrony, and 4) density-dependent processes. Finally, we identify future research directions to advance our understanding of tributary roles and inform conservation actions. In the Missouri River United States, we focus on conservation efforts for the state endangered lake sturgeon, which inhabits large rivers and tributaries in the Midwest and Eastern United States. In the Colorado River, Grand Canyon United States, we focus on conservation efforts for recovery of the federally threatened humpback chub. In the Missouri River, habitat diversity focused on physical habitats such as substrate for reproduction, and deep-water habitats for refuge, whereas augmenting habitat diversity for Colorado River fishes focused on managing populations in tributaries with minimally impaired thermal and flow regimes. Connectivity enhancements in the Missouri River focused on increasing habitat accessibility that may require removal of physical structures like low-head dams; whereas in the Colorado River, the lack of connectivity may benefit native fishes as the disconnection provides refuge from non-native fish predation. Hydrologic variability among tributaries was present in both systems, likely underscoring ecological asynchrony. These case studies also described density dependent processes that could influence success of restoration actions. Although actions to restore populations varied by river system, these examples show that these four portfolio assets can help guide restoration activities across a diverse range of mainstem rivers and their tributaries. Using these assets as a guide, we suggest these can be transferable to other large river-tributary systems.

Frontiers in Environmental Science

An assessment of stream habitat and nutrients in the Elwha River basin: implications for restoration

The Elwha River was once famous for its 10 runs of anadromous salmon which included chinook that reportedly exceeded 45 kilograms. These runs either ceased to exist or were significantly depleted after the construction of the Elwha (1912) and Glines Canyon (1927) Dams, which resulted in the blockage of more than 113 kilometers of mainstem river and tributary habitat. In 1992, in response to the loss of the salmon runs in the Elwha River Basin, President George Bush signed the Elwha River Ecosystem and Fisheries Restoration Act, which authorizes the Secretary of the Interior to remove both dams for ecosystem restoration. The objective of this U.S. Geological Survey (USGS) study was to begin describing baseline conditions for assessing changes that will result from restoration. The first step was to review available physical, chemical, and biological information on the Elwha River Basin. We found that most studies have focused on anadromous fish and habitat and that little information is available on water quality, habitat classification, geomorphic processes, and riparian and aquatic biological communities. There is also a lack of sufficient data on baseline conditions for assessing future changes if restoration occurs. The second component of this study was to collect water-quality and habitat data, filling information gaps. This information will permit a better understanding of the relation between physical habitat and nutrient conditions and changes that may result from salmon restoration. We collected data in the fall of 1997 and found that the concentrations of nitrogen and phosphorous were generally low, with most samples having concentrations below detection limits. Detectable concentrations of nitrogen were associated with sites in the lower reach of the Elwha River, whereas the few detections of phosphorus were at sites throughout the basin. Nutrient data indicate that the Elwha River and its tributaries are oligotrophic. Results of the stream classification indicated that most of the habitat that would be usable by salmon is found in the mainstem of the Elwha River due to natural gradient barriers at the lower end of most tributaries. Habitat is diverse in the mainstem due to large woody debris accumulations and the existence of secondary channels. We concluded that restoring salmon runs to the Elwha River system will affect the ecosystem profoundly. Decaying carcasses of migrating salmon will be the source of large quantities of nutrients to the Elwha River. The complex instream habitat of the mainstem will enhance cycling of these nutrients because carcasses will be retained long enough to be assimilated thereby increasing primary and secondary production, size of immature salmonids, and overall higher salmon recruitment.

Washington

Holocene multidecadal and multicentennial droughts affecting Northern California and Nevada

Continuous, high-resolution δ 18 O records from cored sediments of Pyramid Lake, Nevada, indicate that oscillations in the hydrologic balance occurred, on average, about every 150 years (yr) during the past 7630 calendar years (cal yr). The records are not stationary; during the past 2740 yr, drought durations ranged from 20 to 100 yr and intervals between droughts ranged from 80 to 230 yr. Comparison of tree-ring-based reconstructions of climate change for the past 1200 yr from the Sierra Nevada and the El Malpais region of northwest New Mexico indicates that severe droughts associated with Anasazi withdrawal from Chaco Canyon at 820 cal yr BP (calendar years before present) and final abandonment of Chaco Canyon, Mesa Verde, and the Kayenta area at 650 cal yr BP may have impacted much of the western United States.During the middle Holocene (informally defined in this paper as extending from 8000 to 3000 cal yr BP), magnetic susceptibility values of sediments deposited in Pyramid Lake's deep basin were much larger than late–Holocene (3000–0 cal yr BP) values, indicating the presence of a shallow lake. In addition, the mean δ 18 O value of CaCO 3 precipitated between 6500 and 3430 cal yr BP was 1.6‰ less than the mean value of CaCO 3 precipitated after 2740 cal yr BP. Numerical calculations indicate that the shift in the δ 18 O baseline probably resulted from a transition to a wetter (>30%) and cooler (3–5°C) climate. The existence of a relatively dry and warm middle-Holocene climate in the Truckee River–Pyramid Lake system is generally consistent with archeological, sedimentological, chemical, physical, and biological records from various sites within the Great Basin of the western United States. Two high-resolution Holocene-climate records are now available from the Pyramid and Owens lake basins which suggest that the Holocene was characterized by five climatic intervals. TIC and δ 18 O records from Owens Lake indicate that the first interval in the early Holocene (11,600–10,000 cal yr BP) was characterized by a drying trend that was interrupted by a brief (200 yr) wet oscillation centered at 10,300 cal yr BP. This was followed by a second early-Holocene interval (10,000–8000 cal yr BP) during which relatively wet conditions prevailed. During the early part of the middle Holocene (8000–6500 cal yr BP), high-amplitude oscillations in TIC in Owens Lake and δ 18 O in Pyramid Lake indicate the presence of shallow lakes in both basins. During the latter part of the middle Holocene (6500–3800 cal yr BP), drought conditions dominated, Owens Lake desiccated, and Lake Tahoe ceased spilling to the Truckee River, causing Pyramid Lake to decline. At the beginning of the late Holocene (∼3000 cal yr BP), Lake Tahoe rose to its sill level and Pyramid Lake increased in volume.

Quaternary Science Reviews

Development of a domestic earthquake alert protocol combining the USGS pager and FEMA Hazus systems

The U.S. Geological Survey’s PAGER automated alert system provides rapid (10-20 min) loss estimates in terms of ranges of fatalities and economic impact for all significant earthquakes around the globe. In contrast, FEMA’s Hazus software, which is currently operated manually by FEMA personnel internally within several hours of any large domestic earthquake, provides more detailed loss information quantified in terms of physical damage to the building stock, as well as a broad range of social and economic consequences estimated at a much higher spatial resolution (census-tract level). Hazus was originally designed for FEMA’s mitigation planning, but has been used of late for post-earthquake situational awareness. While the quick alerts generated automatically by the PAGER system are a suitable way to initially identify consequences for earthquakes, there remains a need to provide detailed loss information for facilitating a wider array of post-earthquake decisions among emergency managers. As such, we have developed a prototype summary product that takes advantage of the benefits of both these loss models for significant domestic earthquakes. The signature product, a new twoPAGER report, will serve as a supplement to the widely deployed standard onePAGER product for all significant domestic earthquakes. Page one is the standard, automated PAGER alert content, with summary alert levels for overall fatality and economic loss estimates, as well as summary content on recent earthquakes, structure vulnerability, and historical secondary hazards. As soon as the Hazus run using the most up-to-date version of the USGS ShakeMap is completed and reviewed by both agencies, the second alert page will be generated and delivered via the USGS website as an update. This second page will contain a summary of the more comprehensive Hazus model results, including spatially distributed estimates of: affected population, economic impact, non-fatal injuries, displaced households, and damage to structures (including potential building safety evaluation requirements). This enhancement does not affect the timeliness and the importance of the initial (onePAGER) alerting content for significant earthquakes; however, the twoPAGER report will provide more detailed, damage and loss/impact content for the Nation in the critical hours following a damaging earthquake.

Conference Paper

Global significance of a sub-Moho boundary layer (SMBL) deduced from high-resolution seismic observations

We infer the fine structure of a sub-Moho boundary layer (SMBL) at the top of the lithospheric mantle from high-resolution seismic observations of Peaceful Nuclear Explosions (PNE) on superlong-range profiles in Russia. Densely recorded seismograms permit recognition of previously unknown features of teleseismic propagation of the well known Pn and Sn phases, such as a band of incoherent, scattered, high-frequency seismic energy, developing consistently from station to station, apparent velocities of sub-Moho material, and high-frequency energy to distances of more than 3000 km with a coda band, incoherent at 10 km spacing and yet consistently observed to the end of the profiles. Estimates of the other key elements of the SMBL were obtained by finite difference calculations of wave propagation in elastic 2D models from a systematic grid search through parameter space. The SMBL consists of randomly distributed, mild velocity fluctuations of 2% or schlieren of high aspect ratios (???40) with long horizontal extent (???20 km) and therefore as thin as 0.5 km only; SMBL thickness is 60-100 km. It is suggested that the SMBL is of global significance as the physical base of the platewide observed high-frequency phases Pn and Sn. It is shown that wave propagation in the SMBL waveguide is insensitive to the background velocity distribution on which its schlieren are superimposed. This explains why the Pn and Sn phases traverse geological provinces of various age, heat flow, crustal thickness, and tectonic regimes. Their propagation appears to be independent of age. temperature, pressure, and stress. Dynamic stretching of mantle material during subduction or flow, possibly combined with chemical differentiation have to be considered as scale-forming processes in the upper mantle. However, it is difficult to distinguish with the present sets of Pn/Sn array data whether (and also where) the boundary layer is a frozen-in feature of paleo-processes or whether it is a response to an on-going processes; nevertheless, the derived quantitative estimates of the SMBL properties provide important constraints for any hypothesis on scale-forming processes. Models to be tested by future numerical and field experiments are, for example, repeated subduction-convection stretching of oceanic lithosphere (marble-cake model) and schlieren formation at mid-ocean ridges. It is also proposed that the modeling of the observed blocking of Sn and Pn propagation at active plate margins offers a new tool to study the depth range of tectonics below the crust-mantle boundary. Finally, the deduced schlieren structure of the SMBL closes an important scale gap of three to four orders of magnitude between structural dimensions studied in petrological analysis of mantle samples (xenoliths or outcrop of oceanic lithosphere) and those imaged in classical seismological studies of the lithosphere.

International Geology Review

Pinyon-juniper woodlands

Pinyon-juniper woodlands are one of the largest ecosystems in the Southwest and in the Middle Rio Grande Basin (Fig. 1). The woodlands have been important to the region's inhabitants since prehistoric times for a variety of natural resources and amenities. The ecosystems have not been static; their distributions, stand characteristics, and site conditions have been altered by changes in climatic patterns and human use and, often, abuse. Management of these lands since European settlement has varied from light exploitation and benign neglect, to attempts to remove the trees in favor of forage for livestock, and then to a realization that these lands contain useful resources and should be managed accordingly. Land management agencies are committed to ecosystem management. While there are several definitions of ecosystem management, the goal is to use ecological approaches to create and maintain diverse, productive, and healthy ecosystems (Kaufmann et al. 1994). Ecosystem management recognizes that people are an integral part of the system and that their needs must be considered. Ecological approaches are central to the concept, but our understanding of basic woodland ecology is incomplete, and there are different opinions and interpretations of existing information (Gottfried and Severson 1993). There are many questions concerning proper ecosystem management of the pinyon-juniper woodlands and how managers can achieve these goals (Gottfried and Severson 1993). While the broad concept of ecosystem management generally is accepted, the USDA Forest Service, other public land management agencies, American Indian tribes, and private landowners may have differing definitions of what constitutes desired conditions. Key questions about the pinyon-juniper ecosystems remain unanswered. Some concern the basic dynamics of biological and physical components of the pinyon-juniper ecosystems. Others concern the distribution of woodlands prior to European settlement and changes since the introduction of livestock and fire control. This relates to whether tree densities have been increasing or whether trees are invading grasslands and, to a lesser extent, drier ponderosa pine (Pinus ponderosa) forests. In areas where woodlands were heavily used by American Indians for fuelwood prior to European contact, the advance of pinyon and juniper could represent the slow recovery from intensive use (Samuels and Betancourt 1982). There are numerous questions regarding declines in watershed condition related to changes in pinyon-juniper tree stand densities and to the density and composition of understory vegetation. There are different opinions about proper management of woodland ecosystems. Should these lands be managed for a single resource, such as forage for livestock production, or managed for sustained production of multiple resource products and amenities? Depending on site and stand conditions, the woodlands can produce variable quantities of fuelwood, pinyon nuts, wildlife habitat, forage for livestock, and cover for watershed protection. Management must also consider increasing recreational demands, threatened and endangered species, and protection of archeological sites. Many pinyon-juniper woodland watersheds in New Mexico have unsatisfactory soil and watershed conditions (USDA Forest Service 1993); managers must develop restoration procedures that recognize the value of woodland ecosystems. The concerns, questions, and conflicts surrounding management of pinyon-juniper lands, as well as the ecological foundations of ecosystem management, require that all interested parties reevaluate attitudes toward the woodlands. Ecosystem management goals and concepts recognize diversity. Pinyonjuniper woodlands are diverse, and stand characteristics and site productivities vary. Management objectives and prescriptions must evaluate the potential of each site, and decisions must be based on sound scientific information. This information is often unavailable. Therefore, this paper describes what we do know about the characteristics, distribution, and ecology of pinyon-juniper woodlands, including the effects of natural and human factors, within the southwestern United States and particularly the Middle Rio Grande Basin. It also reviews some past and present management options in this widespread and important vegetation type. The review draws on research and management information from the Rio Grande Basin and from similar areas in the Southwest and adjacent regions. It does not attempt to review all of the relevant literature; additional sources can be found within the articles cited in the References.

General Technical Report

Reference manual for generation and analysis of Habitat Time Series: version II

The selection of an instream flow requirement for water resource management often requires the review of how the physical habitat changes through time. This review is referred to as 'Time Series Analysis." The Tune Series Library (fSLIB) is a group of programs to enter, transform, analyze, and display time series data for use in stream habitat assessment. A time series may be defined as a sequence of data recorded or calculated over time. Examples might be historical monthly flow, predicted monthly weighted usable area, daily electrical power generation, annual irrigation diversion, and so forth. The time series can be analyzed, both descriptively and analytically, to understand the importance of the variation in the events over time. This is especially useful in the development of instream flow needs based on habitat availability. The TSLIB group of programs assumes that you have an adequate study plan to guide you in your analysis. You need to already have knowledge about such things as time period and time step, species and life stages to consider, and appropriate comparisons or statistics to be produced and displayed or tabulated. Knowing your destination, you must first evaluate whether TSLIB can get you there. Remember, data are not answers . This publication is a reference manual to TSLIB and is intended to be a guide to the process of using the various programs in TSLIB. This manual is essentially limited to the hands-on use of the various programs. a TSLIB use interface program (called RTSM) has been developed to provide an integrated working environment where the use has a brief on-line description of each TSLIB program with the capability to run the TSLIB program while in the user interface. For information on the RTSM program, refer to Appendix F. Before applying the computer models described herein, it is recommended that the user enroll in the short course "Problem Solving with the Instream Flow Incremental Methodology (IFIM)." This course is offered by the Aquatic Systems Branch of the National Ecology Research Center. For more information about the TSLIB software, refer to the Memorandum of Understanding. Chapter 1 provides a brief introduction to the Instream Flow Incremental Methodology and TSLIB. Other chapters in this manual provide information on the different aspects of using the models. The information contained in the other chapters includes (2) acquisition, entry, manipulation, and listing of streamflow data; (3) entry, manipulation, and listing of the habitat-versus-streamflow function; (4) transferring streamflow data; (5) water resources systems analysis; (6) generation and analysis of daily streamflow and habitat values; (7) generation of the time series of monthly habitats; (8) manipulation, analysis, and display of month time series data; and (9) generation, analysis, and display of annual time series data. Each section includes documentation for the programs therein with at least one page of information for each program, including a program description, instructions for running the program, and sample output. The Appendixes contain the following: (A) sample file formats; (B) descriptions of default filenames; (C) alphabetical summary of batch-procedure files; (D) installing and running TSLIB on a microcomputer; (E) running TSLIB on a CDC Cyber computer; (F) using the TSLIB user interface program (RTSM); and (G) running WATSTORE on the USGS Amdahl mainframe computer. The number for this version of TSLIB--Version II-- is somewhat arbitrary, as the TSLIB programs were collected into a library some time ago; but operators tended to use and manage them as individual programs. Therefore, we will consider the group of programs from the past that were only on the CDC Cyber computer as Version 0; the programs from the past that were on both the Cyber and the IBM-compatible microcomputer as Version I; and the programs contained in this reference manual as Version II.

Report

Review of selected features of the natural system model, and suggestions for applications in south Florida

A study was conducted to review selected features of the Natural System Model, version 4.3. The Natural System Model is a regional-scale model that uses recent climatic data and estimates of historic vegetation and topography to simulate pre-canal-drainage hydrologic response in south Florida. Equations used to represent the hydrologic system and the numerical solution of these equations in the model were documented and reviewed. Convergence testing was performed using 1965 input data, and selected other aspects of the model were evaluated. Some conclusions from the evaluation of the Natural System Model include the following observations. Simulations were generally insensitive to the temporal resolution used in the model. However, reduction of the computational cell size from 2-mile by 2-mile to 2/3-mile by 2/3-mile resulted in a decrease in spatial mean ponding depths for October of 0.35 foot for a 3-hour time step. Review of the computer code indicated that there is no limit on the amount of water that can be transferred from the river system to the overland flow system, on the amount of seepage from the river to the ground-water system, on evaporation from the river system, or on evapotranspiration from the overland-flow system. Oscillations of 0.2 foot or less in simulated river stage were identified and attributed to a volume limiting function which is applied in solution of the overland-flow equations. The computation of the resistance coefficient is not consistent with the computation of overland-flow velocity. Ground-water boundary conditions do not always ensure a no-flow condition at the boundary. These inconsistencies had varying degrees of effects on model simulations, and it is likely that simulations longer than 1 year are needed to fully identify effects. However, inconsistencies in model formulations should not be ignored, even if the effects of such errors on model results appear to be small or have not been clearly defined. The Natural System Model can be a very useful tool for estimating pre-drainage hydrologic response in south Florida. The model includes all of the important physical processes needed to simulate a water balance. With a few exceptions, these hydrologic processes are represented in a reasonable manner using empirical, semiempirical, and mechanistic relations. The data sets that have been assembled to represent physical features, and hydrologic and meteorological conditions are quite extensive in their scope. Some suggestions for model application were made. Simulation results from the Natural System Model need to be interpreted on a regional basis, rather than cell by cell. The available evidence suggests that simulated water levels should be interpreted with about a plus or minus 1 foot uncertainty. It is probably not appropriate to use the Natural System Model to estimate pre-drainage discharges (as opposed to hydroperiods and water levels) at a particular location or across a set of adjacent computational cells. All simulated results for computational cells within about 10 miles of the model boundaries have a higher degree of uncertainty than results for the interior of the model domain. It is most appropriate to interpret the Natural System Model simulation results in connection with other available information. Stronger linkages between hydrologic inputs to the Everglades and the ecological response of the system would enhance restoration efforts.

Florida

Crustal architecture beneath the southern Midcontinent (USA) and controls on Mesoproterozoic iron-oxide mineralization from 3D geophysical models

Several types of critical mineral-bearing ore deposits in the southern Midcontinent region of the U.S. are hosted in Mesoproterozoic igneous rocks largely concealed beneath Paleozoic cover. Discerning the architecture of igneous intrusions and volcanic centers in the crust is fundamental to understanding the geologic evolution of this vast region and its mineral resources. To advance the understanding of the geologic framework beneath the Southeast Missouri Iron Metallogenic Province, we invert continental-scale magnetic and gravity anomaly data to three-dimensional (3D) physical property models. The regional models image altered and mineralized igneous rocks near the Precambrian basement surface and underlying intrusive complexes that extend down to the Moho. At shallow crustal levels, our models confirm that iron oxide-apatite ± rare earth element (IOA±REE) deposits and iron oxide-copper-gold ± cobalt (IOCG) deposits occur within or near the edges of large low density/low susceptibility early Mesoproterozoic (ca 1.4 Ga) silicic calderas and (ca 1.3 Ga) granitic plutons. Previous isotopic and geochemical studies conclude that the iron deposits and their volcanic host rocks originated from mantle-derived and crustal melts that erupted during regional extension. Extension was associated with thermal event(s) that produced the large-scale silicic magmatism related to the ca 1.45 Ga Eastern Granite Rhyolite Province (EGRP) and the 1.35 Ga Southern Granite Rhyolite Province (SGRP). We postulate that early in the evolution of the EGRP, several trans-crustal magmatic plumbing systems developed that are evident in the 3D models. The Southeast Missouri Metallogenic Province is underlain by one such magmatic system that is expressed as a northwest-trending ~ 50 km-wide by 200 km-long elongate track of high susceptibility at deep crustal levels. The high susceptibility corridor splays upward through the crust to the Precambrian surface where the iron deposits are the epigenetic manifestation of this magmatic event. Our findings confirm that the iron deposits, with no distinct connection at the surface, are connected to one large magmatic system at depth. We propose that other similar susceptibility tracks, which are present along the top of the mantle, mark additional feeder zones that allowed magma to ascend to the main eruptive centers that produced the Granite Rhyolite Provinces. The early Mesoproterozoic extensional tectonic framework established crustal-scale pathways that controlled the distribution of subsequent magmatic activity, including the ca 1.4 Ga calderas and underlying intrusions, ca 1.3 Ga silicic plutons and Phanerozoic alkaline intrusions. If these interpretations are correct, our study has identified large areas that are prospective for critical mineral-bearing ore deposits and, importantly, suggests that the Mesoproterozoic architecture may have influenced subsequent magmatism and hydrothermal activity in the southern Midcontinent of the U.S.

Ore Geology Reviews

Coastal Circulation and Sediment Dynamics in War-in-the-Pacific National Historical Park, Guam; measurements of waves, currents, temperature, salinity, and turbidity, June 2007-January 2008

Flow in and around coral reefs affects a number of physical, chemical and biologic processes that influence the health and sustainability of coral reef ecosystems. These range from the residence time of sediment and contaminants to nutrient uptake and larval retention and dispersal. As currents approach a coast they diverge to flow around reef structures, causing high horizontal and vertical shear. This can result in either the rapid advection of material in localized jets, or the retention of material in eddies that form in the lee of bathymetric features. The high complexity and diversity both within and between reefs, in conjunction with past technical restrictions, has limited our understanding of the nature of flow and the resulting flux of physical, chemical, and biologic material in these fragile ecosystems. Sediment, nutrients, and other pollutants from a variety of land-based activities adversely impact many coral reef ecosystems in the U.S. and around the world. These pollutants are transported in surface water runoff, groundwater seepage, and atmospheric fallout into coastal waters, and there is compelling evidence that the sources have increased globally as a result of human-induced changes to watersheds. In Guam, and elsewhere on U.S. high islands in the Pacific and Caribbean, significant changes in the drainage basins due to agriculture, feral grazing, fires, and urbanization have in turn altered the character and volume of land-based pollution released to coral reefs. Terrigenous sediment run-off (and the associated nutrients and contaminants often absorbed to it) and deposition on coral reefs are recognized to potentially have significant impact on coral health by blocking light and inhibiting photosynthesis, directly smothering and abrading coral, and triggering increases in macro algae. Studies that combine information on watershed, surface water- and groundwater-flow, transport and fate of sediment and other pollutants in the reef environment, and their impact on reef health and ecology are essential for effective reef management. Two of the main anthropogenic activities along west-central Guam's coastline that may impact the region's coral reef ecosystems include pollution and coastal land use/development, as discussed in the review by Porter and others (2005). The pollution threats include point-sources, such as municipal wastewater (Northern District, Hagatna, Naval Station Guam, and Agat-Santa Rita Waster Water Treatment Plants), cooling water (Tanguisson Steam and Cabras Power Plants), and numerous storm water, ballast water, and tank bottom draw outfalls; nonpoint sources include septic systems, urban runoff, illegal dumping, and groundwater discharges. Poor land-use practices include development without the use of runoff management measures, increased areal extent of impervious surfaces and decreased extent of vegetative barriers, and recreational off-road vehicle use. Furthermore, feral ungulates and illegal wildfires remove protective vegetative cover and generally result in increased soil erosion. While anthropogenic point-sources have been reduced in many areas due to better management practices, nonpoint sources have either stayed constant or increased. Between 1975 and 1999, it is estimated that Guam lost more than a quarter of its tree cover, and more than 750 wildfires each year have resulted in a greater proportion of badlands and other erosion-prone land surfaces with high erosion rates (Forestry and Soil Resources Division, 1999). Approximately 1.8 square kilometers (km2) of Asan Bay, west-central Guam, lies within the National Park Service's (NPS) War-in-the-Pacific National Historical Park's (WAPA) Asan Unit; the bay is the sink for material coming out of the Asan watershed. Anthropogenic modifications of the watersheds adjacent to Asan Bay, which include intentionally-set wildfires, construction, and agriculture (Minton, 2005), are believed to have increased over the past 25

Open-File Report

Lithospheric buoyancy and continental intraplate stresses

Lithospheric buoyancy, the product of lithospheric density and thickness, is an important physical property that influences both the long-term stability of continents and their state of stress. We have determined lithospheric buoyancy by applying the simple isostatic model of Lachenbruch and Morgan (1990). We determine the crustal portion of lithospheric buoyancy using the USGS global database of more than 1700 crustal structure determinations (Mooney et al., 2002), which demonstrates that a simple relationship between crustal thickness and surface elevation does not exist. In fact, major regions of the crust at or near sea level (0-200 m elevation) have crustal thicknesses that vary between 25 and 55 km. Predicted elevations due to the crustal component of buoyancy in the model exceed observed elevations in nearly all cases (97% of the data), consistent with the existence of a cool lithospheric mantle lid that is denser than the asthenosphere on which it floats. The difference between the observed and predicted crustal elevation is assumed to be equal to the decrease in elevation produced by the negative buoyancy of the mantle lid. Mantle lid thickness was first estimated from the mantle buoyancy and a mean lid density computed using a basal crust temperature determined from extrapolation of surface heat flow, assuming a linear thermal gradient in the mantle lid. The resulting values of total lithosphere thickness are in good agreement with thicknesses estimated from seismic data, except beneath cratonic regions where they are only 40-60% of the typical estimates (200-350 km) derived from seismic data. This inconsistency is compatible with petrologic data and tomography and geoid analyses that have suggested that cratonic mantle lids are ??? 1% less dense than mantle lids elsewhere. By lowering the thermally determined mean mantle lid density in cratons by 1%, our model reproduces the observed 200-350+ km cratonic lithospheric thickness. We then computed gravitational potential energy by taking a vertical integral over the computed lithosphere density. Our computed values suggest that the thick roots beneath cratons lead to strong negative potential energy differences relative to surrounding regions, and hence exert compressive stresses superimposed on the intraplate stresses derived from plate boundary forces. Forces related to this lithosphere structure thus may explain the dominance of reverse-faulting earthquakes in cratons. Areas of high elevation and a thin mantle lid (e.g., western U.S. Basin and Range, East African rift, and Baikal rift) are predicted to be in extension, consistent with the observed stress regime in these areas.

International Geology Review

Earthquake locations determined by the Southern Alaska seismograph network for October 1971 through May 1989

This report describes the instrumentation and evolution of the U.S. Geological Survey’s regional seismograph network in southern Alaska, provides phase and hypocenter data for seismic events from October 1971 through May 1989, reviews the location methods used, and discusses the completeness of the catalog and the accuracy of the computed hypocenters. Included are arrival time data for explosions detonated under the Trans-Alaska Crustal Transect (TACT) in 1984 and 1985. The U.S. Geological Survey (USGS) operated a regional network of seismographs in southern Alaska from 1971 to the mid 1990s. The principal purpose of this network was to record seismic data to be used to precisely locate earthquakes in the seismic zones of southern Alaska, delineate seismically active faults, assess seismic risks, document potential premonitory earthquake phenomena, investigate current tectonic deformation, and study the structure and physical properties of the crust and upper mantle. A task fundamental to all of these goals was the routine cataloging of parameters for earthquakes located within and adjacent to the seismograph network. The initial network of 10 stations, 7 around Cook Inlet and 3 near Valdez, was installed in 1971. In subsequent summers additions or modifications to the network were made. By the fall of 1973, 26 stations extended from western Cook Inlet to eastern Prince William Sound, and 4 stations were located to the east between Cordova and Yakutat. A year later 20 additional stations were installed. Thirteen of these were placed along the eastern Gulf of Alaska with support from the National Oceanic and Atmospheric Administration (NOAA) under the Outer Continental Shelf Environmental Assessment Program to investigate the seismicity of the outer continental shelf, a region of interest for oil exploration. Since then the region covered by the network remained relatively fixed while efforts were made to make the stations more reliable through improved electronic instrumentation and strengthened antenna systems. The majority of the stations installed since 1980 were operated only temporarily (from one to several years) for special studies in various areas within the network. Due to reduced funding, the network was trimmed substantially in the summer of 1985 with the closure of 15 stations, 13 of which were located in and around the Yakataga seismic gap. To further reduce costs, two telephone circuits were dropped and multiple radio relays were installed in their place. This economy reduced the reliability of these telemetry links. In addition, data collection from the areas around Cordova and Yakutat was compromised by the necessity of relying on triggered event recording using PC-based systems (Rogers, 1993) that were not fully developed and which proved to be less reliable than anticipated. The principal means of recording throughout the time period of this catalog was 20-channel oscillographs on 16-mm film (Teledyne Geotech Develocorder, Model RF400 and 4000D). Initially one Develocorder was operated at the USGS Alaskan headquarters in Anchorage, but in 1972 recording was shifted to the National Oceanic and Atmospheric Administration (NOAA) Palmer Observatory (currently the West Coast and Alaska Tsunami Warning Center). The Develocorders were turned off at the end of May 1989, and after that time recording was done in digital format at the Geophysical Institute of the University of Alaska in Fairbanks (GIUA). Thus, this catalog covers the entire period of film recording.

Alaska

Geochemistry, Comparative Analysis, and Physical and Chemical Characteristics of the Thermal Waters East of Hot Springs National Park, Arkansas, 2006-09

A study was conducted by the U.S Geological Survey in cooperation with the Arkansas State Highway and Transportation Department to characterize the source and hydrogeologic conditions responsible for thermal water in a domestic well 5.5 miles east of Hot Springs National Park, Hot Springs, Arkansas, and to determine the degree of hydraulic connectivity between the thermal water in the well and the hot springs in Hot Springs National Park. The water temperature in the well, which was completed in the Stanley Shale, measured 33.9 degrees Celsius, March 1, 2006, and dropped to 21.7 degrees Celsius after 2 hours of pumping - still more than 4 degrees above typical local groundwater temperature. A second domestic well located 3 miles from the hot springs in Hot Springs National Park was discovered to have a thermal water component during a reconnaissance of the area. This second well was completed in the Bigfork Chert and field measurement of well water revealed a maximum temperature of 26.6 degrees Celsius. Mean temperature for shallow groundwater in the area is approximately 17 degrees Celsius. The occurrence of thermal water in these wells raised questions and concerns with regard to the timing for the appearance of the thermal water, which appeared to coincide with construction (including blasting activities) of the Highway 270 bypass-Highway 70 interchange. These concerns were heightened by the planned extension of the Highway 270 bypass to the north - a corridor that takes the highway across a section of the eroded anticlinal complex responsible for recharge to the hot springs of Hot Springs National Park. Concerns regarding the possible effects of blasting associated with highway construction near the first thermal well necessitated a technical review on the effects of blasting on shallow groundwater systems. Results from available studies suggested that propagation of new fractures near blasting sites is of limited extent. Vibrations from blasting can result in rock collapse for uncased wells completed in highly fractured rock. However, the propagation of newly formed large fractures that potentially could damage well structures or result in pirating of water from production wells appears to be of limited possibility based on review of relevant studies. Characteristics of hydraulic conductivity, storage, and fracture porosity were interpreted from flow rates observed in individual wells completed in the Bigfork Chert and Stanley Shale; from hydrographs produced from continuous measurements of water levels in wells completed in the Arkansas Novaculite, the Bigfork Chert, and Stanley Shale; and from a potentiometric-surface map constructed using water levels in wells throughout the study area. Data gathered from these three separate exercises showed that fracture porosity is much greater in the Bigfork Chert relative to that in the Stanley Shale, shallow groundwater flows from elevated recharge areas with exposures of Bigfork Chert along and into streams within the valleys formed on exposures of the Stanley Shale, and there was no evidence of interbasin transfer of groundwater within the shallow flow system. Fifteen shallow wells and two cold-water springs were sampled from the various exposed formations in the study area to characterize the water quality and geochemistry for the shallow groundwater system and for comparison to the geochemistry of the hot springs in Hot Springs National Park. For the quartz formations (novaculite, chert, and sandstone formations), total dissolved solids concentrations were very low with a median concentration of 23 milligrams per liter, whereas the median concentration for groundwater from the shale formations was 184 milligrams per liter. Ten hot springs in Hot Springs National Park were sampled for the study. Several chemical constituents for the hot springs, including pH, total dissolved solids, major cations and anions, and trace metals, show similarity with the shale formations

Scientific Investigations Report

Magnitude and Frequency of Rural Floods in the Southeastern United States, through 2006: Volume 2, North Carolina

Reliable estimates of the magnitude and frequency of floods are required for the economical and safe design of transportation and water-conveyance structures. A multistate approach was used to update methods for estimating the magnitude and frequency of floods in rural, ungaged basins in North Carolina, South Carolina, and Georgia that are not substantially affected by regulation, tidal fluctuations, or urban development. In North Carolina, annual peak-flow data available through September 2006 were available for 584 sites; 402 of these sites had a total of 10 or more years of systematic record that is required for at-site, flood-frequency analysis. Following data reviews and the computation of 20 physical and climatic basin characteristics for each station as well as at-site flood-frequency statistics, annual peak-flow data were identified for 363 sites in North Carolina suitable for use in this analysis. Among these 363 sites, 19 sites had records that could be divided into unregulated and regulated/ channelized annual peak discharges, which means peak-flow records were identified for a total of 382 cases in North Carolina. Considering the 382 cases, at-site flood-frequency statistics are provided for 333 unregulated cases (also used for the regression database) and 49 regulated/channelized cases. The flood-frequency statistics for the 333 unregulated sites were combined with data for sites from South Carolina, Georgia, and adjacent parts of Alabama, Florida, Tennessee, and Virginia to create a database of 943 sites considered for use in the regional regression analysis. Flood-frequency statistics were computed by fitting logarithms (base 10) of the annual peak flows to a log-Pearson Type III distribution. As part of the computation process, a new generalized skew coefficient was developed by using a Bayesian generalized least-squares regression model. Exploratory regression analyses using ordinary least-squares regression completed on the initial database of 943 sites resulted in defining five hydrologic regions for North Carolina, South Carolina, and Georgia. Stations with drainage areas less than 1 square mile were removed from the database, and a procedure to examine for basin redundancy (based on drainage area and periods of record) also resulted in the removal of some stations from the regression database. Flood-frequency estimates and basin characteristics for 828 gaged stations were combined to form the final database that was used in the regional regression analysis. Regional regression analysis, using generalized least-squares regression, was used to develop a set of predictive equations that can be used for estimating the 50-, 20-, 10-, 4-, 2-, 1-, 0.5-, and 0.2-percent chance exceedance flows for rural ungaged, basins in North Carolina, South Carolina, and Georgia. The final predictive equations are all functions of drainage area and the percentage of drainage basin within each of the five hydrologic regions. Average errors of prediction for these regression equations range from 34.0 to 47.7 percent. Discharge estimates determined from the systematic records for the current study are, on average, larger in magnitude than those from a previous study for the highest percent chance exceedances (50 and 20 percent) and tend to be smaller than those from the previous study for the lower percent chance exceedances when all sites are considered as a group. For example, mean differences for sites in the Piedmont hydrologic region range from positive 0.5 percent for the 50-percent chance exceedance flow to negative 4.6 percent for the 0.2-percent chance exceedance flow when stations are grouped by hydrologic region. Similarly for the same hydrologic region, median differences range from positive 0.9 percent for the 50-percent chance exceedance flow to negative 7.1 percent for the 0.2-percent chance exceedance flow. However, mean and median percentage differences between the estimates from the previous and curre

North Carolina

Late-kinematic timing of orogenic gold deposits and significance for computer-based exploration techniques with emphasis on the Yilgarn Block, Western Australia

Orogenic gold deposits are a widespread coherent group of epigenetic ore deposits that are sited in accretionary or collisional orogens. They formed over a large crustal-depth range from deep-seated low-salinity H 2 O–CO 2 ±CH 4 ±N 2 ore fluids and with Au transported as thio-complexes. Regional structures provide the main control on deposit distribution. In many terranes, first-order faults or shear zones appear to have controlled regional fluid flow, with greatest ore-fluid fluxes in, and adjacent to, lower-order faults, shear zones and/or large folds. Highly competent and/or chemically reactive rocks are the most common hosts to the larger deposits. Focusing of supralithostatic ore fluids into dilatant zones appears to occur late during the evolutionary history of the host terranes, normally within D 3 or D 4 in a D 1 –D 4 deformation sequence. Reactivation of suitably oriented pre-existing structures during a change in far-field stress orientation is a factor common to many deposits, and repeated reactivation may account for multiple mineralization episodes in some larger deposits. Absolute robust ages of mineralization support their late-kinematic timing, and, in general, suggest that deposits formed diachronously towards the end of the 100 to 200 m.y. long evolutionary history of hosting orogens. For example, in the Yilgarn Block, a region specifically emphasised in this study, orogenic gold deposits formed in the time interval between 40 and 90 m.y., with most about 60 to 70 m.y., after the youngest widespread basic-ultrabasic volcanism and towards the end of felsic magmatism. The late timing of orogenic gold deposits is pivotal to geologically-based exploration methodologies. This is because the present structural geometries of: (i) the deposits, (ii) the hosting goldfields, and (iii) the enclosing terranes are all essentially similar to those during gold mineralization, at least in their relative position to each other. Thus, interpretation of geological maps and cross-sections and three-dimensional models can be used to accurately simulate the physical conditions that existed at the time of ore deposition. It is particularly significant that the deposits are commonly related to repetitive and predictable geometries, such as structural heterogeneities within or adjacent to first-order structures, around rigid granitoid bodies, or in specific “locked-up” fold-thrust structures. Importantly, the two giant greenstone-hosted goldfields, Kalgoorlie and Timmins, show a remarkably similar geometry at the regional scale. Computer-based stress mapping and GIS-based prospectivity mapping are two computer-based quantitative methodologies that can utilize and take advantage of the late timing aspect of this deposit type to provide important geological aids in exploration, both in broad regions and more localized goldfields. Both require an accurate and consistent solid geology map, stress mapping requires knowledge of the far-field stresses during mineralization, and the empirical prospectivity mapping requires data from a significant number of known deposits in the terrane. The Kalgoorlie Terrane, in the Yilgarn Block, meets these criteria, and illustrates the potential of these methodologies in the exploration for orogenic gold deposits. Low minimum stress anomalies, interpreted to represent dilational zones during gold-related deformation, coincide well with the positions of known goldfields rather than individual gold deposits in the terrane, and there are additional as-yet unexplained anomalies. The prospectivity analysis confirms that predictable and repetitive factors controlling the siting of deposits are: (i) proximity to, and orientation and curvature of, granitoid-greenstone contacts, (ii) proximity to segments of crustal faults which strike in a preferred direction, (iii) proximity to specific lithological contacts which have similar preferred strike, (iv) proximity to anticlinal structures, and (v) the presence of preferred reactive host rocks (e.g., dolerite). The prospectivity map defines a series of anomalous areas, which broadly conform to those of the stress map (>78% correspondence). The most prospective category on this map covers less than 0.3% of the greenstone belts and yet hosts 16% of the known deposits, which have produced>80% of known gold. Thus, it discriminates in favour of the larger economically more-attractive deposits in the terrane. The successful application of stress mapping and prospectivity mapping to geology-based exploration for orogenic gold deposits indicates that more quantitative analysis of geological map data is a profitable line of research. The computer-based nature of these methodologies is ideal for the production of an ultimate, integrated, deposit target map, which can be compared to other, more conventional, targeting parameters such as geophysical and geochemical anomalies. Such an integrated strategy appears the way forward in the increasingly difficult task of cost-effective global exploration for orogenic gold deposits in poorly exposed terranes.

Western Australia

Groundwater quality of the Gulf Coast aquifer system, Houston, Texas, 2010

During March–December 2010, the U.S. Geological Survey, in cooperation with the city of Houston, collected source-water samples from 60 municipal supply wells in the Houston area. These data were collected as part of an ongoing study to determine concentrations, spatial extent, and associated geochemical conditions that might be conducive for mobility and transport of selected naturally occurring contaminants (selected trace elements and radionuclides) in the Gulf Coast aquifer system in the Houston area. In the summers of 2007 and 2008, a reconnaissance-level survey of these constituents in untreated water from 28 municipal supply wells was completed in the Houston area. Included in this report are the complete analytical results for 47 of the 60 samples collected in 2010—those results which were received from the laboratories and reviewed by the authors as of December 31, 2010. All of the wells sampled were screened in the Gulf Coast aquifer system; 22 were screened entirely in the Evangeline aquifer, and the remaining 25 wells contained screened intervals that intersected both Evangeline and Chicot aquifers. The data documented in this report were collected as part of an ongoing study to characterize source-water-quality conditions in untreated groundwater prior to drinking-water treatment. An evaluation of contaminant occurrence in source water provides background information regarding the presence of a contaminant in the environment. Because source-water samples were collected prior to any treatment or blending that potentially could alter contaminant concentrations, the water-quality results documented by this report represent the quality of the source water, not the quality of finished drinking water provided to the public. Samples were analyzed for major ions (calcium, magnesium, potassium, sodium, bromide, chloride, fluoride, silica, and sulfate), residue on evaporation (dissolved solids), trace elements (arsenic, barium, boron, chromium, iron, lithium, manganese, molybdenum, selenium, strontium, and vanadium), and selected radionuclides (gross alpha- and beta-particle activity [at 72 hours and 30 days], carbon-14, radium-226, radon-222, and uranium). Field measurements were made of selected physicochemical (relating to both physical and chemical) properties (oxidation-reduction potential, turbidity, dissolved-oxygen concentration, pH, specific conductance, water temperature, and alkalinity) and unfiltered sulfides. Similar to the results from the reconnaissance survey, physicochemical properties, major ions, and trace elements varied considerably. The ranges of selected physicochemical properties were as follows: oxidation-reduction potential ranged from -173 to 466 millivolts, dissolved oxygen ranged from less than 0.1 to 4.4 milligrams per liter, pH ranged from 7.2 to 7.8, specific conductance ranged from 439 to 724 microsiemens per centimeter at 25 degrees Celsius, and alkalinity ranged from 159 to 276 milligrams per liter as calcium carbonate. The largest ranges in concentration for filtered major ion constituents were obtained for cations sodium and calcium and for anions chloride and sulfate. Arsenic concentrations measured in samples from the 47 wells ranged from 1.6 to 23.5 micrograms per liter. The maximum concentration of arsenic (23.5 micrograms per liter) was measured in the source-water sample from well LJ-65-12-328. Quantifiable concentrations of barium, boron, lithium, molybdenum, and strontium were measured in all 47 filtered, source-water samples. Quantifiable concentrations of manganese were measured in 46 source-water samples, and an estimated concentration of manganese was measured in 1 sample. Chromium, iron, selenium, and vanadium were detected in 24 or more of the 47 source-water samples. Gross alpha-particle activities and beta-particle activities for all 47 samples were analyzed at 72 hours after sample collection and again at 30 days after sample collection, allowing for the measurement of the activity of short-lived isotopes. Gross alpha-particle activities reported in this report were not adjusted for activity contributions by radon or uranium and, therefore, are conservatively high estimates if compared to the U.S. Environmental Protection Agency National Primary Drinking Water Regulation for adjusted gross alpha-particle activity. The gross alpha-particle activities at 30 days in the samples ranged from R0.60 to 25.5 picocuries per liter and at 72 hours ranged from 2.58 to 39.7 picocuries per liter, and the "R" preceding the value of 0.60 picocuries per liter refers to a nondetected result less than the sample-specific critical level. Gross beta-particle activities measured at 30 days ranged from 1.17 to 14.4 picocuries per liter and at 72 hours ranged from 1.97 to 4.4 picocuries per liter. Filtered uranium was detected in quantifiable amounts in all of the 47 wells sampled. The uranium concentrations ranged from 0.03 to 42.7 micrograms per liter. One sample was analyzed for carbon-14, and the amount of modern atmospheric carbon was reported as 0.2 percent. Six source-water samples collected from municipal supply wells were analyzed for radium-226, and all of the concentrations were considered detectable concentrations (greater than their associated sample-specific critical level). Three source-water samples collected were analyzed for radon-222, and all of the concentrations were substantially greater than the associated sample-specific critical level.

Texas

Managed aquifer recharge in snow-fed river basins: What, why and how?

What does climate change mean for snow-fed river basins? Climate change poses unique challenges in snow-fed river basins across the western United States because the majority of water supply originates as snow (Dettinger, Udall, & Georgakakos, 2015). In the Sierra Nevada, recent observations include changes in snow accumulation and snowmelt, and shifts in peak streamflow timing (Barnhart et al., 2016; Hatchett et al., 2017; Kim & Jain, 2010; McCabe, Wolock, & Valentin, 2018; Mote, Li, Lettenmaier, Xiao, & Engel, 2018). Such changes upstream alter surface water deliveries downstream, as well as groundwater recharge utilized as both primary and supplemental water supply (Godsey et al., 2014; Harpold, 2016; Jasechko et al., 2014). basin where snowmelt runoff produces substantial water supply to meet diverse agricultural, environmental and urban water demand (Figure 1). The East and West Forks join at the confluence of the Carson River near the north end of the Carson Valley, a rich agricultural region (40,000 acres) that grows primarily alfalfa hay. The majority of irrigators rely on surface water delivered through a network of earthen ditches constructed in the mid-19th and early 20th centuries. Flow through these earthen networks and the practice of flood irrigation contribute significantly to groundwater recharge. Because no upstream surface water reservoirs exist, snowpack that accumulates through winter and melts slowly through spring has acted as a “natural” reservoir, providing ample supply through the summer irrigati agricultural, environmental and urban water demand (Figure 1). The East and West Forks join at the confluence of the Carson River near the north end of the season. Some irrigators have permitted access to supplemental groundwater that is useful during periods of drought for augmenting shortfalls in surface water delivery. Groundwater is the primary source of municipal and industrial water supply for surrounding communities (e.g., Carson City, Minden, Gardnerville, Dayton). Across the basin, water use is highly regulated through federal, tribal, state and local water-sharing agreements based on prior appropriation doctrine (Wilds, 2014). Carson River surface water allocations follow the Alpine Decree, initiated by the United States Department of Interior in 1925 and signed into law in 1980, following 55 years of litigation, to adjudicate surface water rights to individual parties (NDWP, 1999). The Alpine Decree acknowledges return flows to lower river segments, and thus each river segment is distributed autonomously. This means that the most junior water right on an upper segment can be fulfilled before considering the most senior water right on a lower segment. Ultimately, the ruling is at the discretion of the Federal Water Master to satisfy the needs of each water right Downstream of Carson Valley, surface water flows are stored in Lahontan Reservoir, the nation’s first desert reclamation project (est. 1906), where releases are managed to meet the Newland’s Project irrigation water demand and for environmental use on the Stillwater National Wildlife Refuge. Flows from the Carson River are supplemented through diversions from the Truckee River via the Truckee Canal, resulting in a trans-basin water supply system. How is the Water for the Seasons research program informing snow-fed river basin communities? In the Truckee-Carson River System, researchers and local water managers are working together to assess climate change impacts to water supply and explore how model simulations can produce useful information to support local climate adaptation. Twelve key water managers represent agricultural, environmental, urban and regulatory water-use communities, and bring to the table diverse input and perspectives on how to adapt to climate change. Hydrologists use this input to craft scenarios and simulations that meet the information needs of local water managers. Biannual workshops provide an opportunity for information exchange, where researchers and key water managers generate new knowledge of river system function. That is, researchers share results of models that examine the physical potential, and managers validate the on-the-ground potential, further informing the research process. Coincident to this research program, the region faced a prolonged drought period (2012-2016) with historically low snowpack, followed by a historic wet year (2017) that brought winter and spring flooding as a result of atmospheric river storm events (Sterle et al., 2019). For the Carson River, an important observation made by managers was that peak streamflow that had traditionally coincided with peak irrigation demands, had shifted to earlier in the spring, with summer baseflow also decreasing (Sterle & Singletary, 2017). Managers shared with researchers concerns over potential future impacts that changing snowpack will have on surface water deliveries and reliance on groundwater, as the region’s population and economy continue to grow. During workshops that occurred over this period, local water managers and researchers discussed ways to evaluate water distribution and use that honors the existing legal framework and accounts for changing snowpack regimes (amount, rain versus snow, timing). In response to managers growing interest, researchers introduced the concept of managed aquifer recharge as one potential strategy to adapt and enhance regional water sustainability. What is managed aquifer recharge? Simply stated, managed aquifer recharge is the intentional recharge of structures to spread water over agricultural lands, allowing water to naturally infiltrate into the groundwater system (Bouwer, 1999; Niswonger et al., 2017). The latter may occur during the irrigation season by applying excess water, or during the nonirrigation season when evapotranspiration losses are low. Figure 2 illustrates managed aquifer recharge in a snow-fed river basin, where streamflow generated from snowmelt runoff is diverted to agricultural lands to recharge the aquifer. Such flood irrigation practices, including water delivery through earthen ditch networks, provide incidental but significant aquifer recharge through seepage and deep drainage beneath fields (Niswonger, Allander, & Jeton, 2014). The effects of managed aquifer recharge can vary depending on the location and intensity of practice. For example, implementing managed aquifer recharge water into the groundwater system (Dillon, 2009). This differs from the incidental recharge that may occur as part of normal irrigation practices. Managed recharge may occur by injection into the aquifer through existingwells, or by using existing conveyance adjacent to/along the river’s floodplain has the potential to enhance late-season instream flows due to increased return flows, resulting in greater downstream deliveries as well as improving ecological conditions (Niswonger et al., 2017). Implementing managed aquifer recharge away from the river’s floodplain has the potential to enhance groundwater supply which is increasingly relied upon during surface water shortage (Green et al., 2011), by storing water in available aquifer space in the deep aquifer. At the basin scale, managed aquifer recharge may lead to regional groundwater sustainability. Is the Carson River Basin a candidate for managed aquifer recharge? The physical limitations to implementing managed aquifer recharge in the Carson River Basin hinges on three key factors. The first factor relates to the physical connectivity between rivers and streams, and the irrigation delivery network of canals and ditches that divert water to agricultural lands (Niswonger et al., 2017). In the Carson River Basin the mechanisms for getting water to fields is already in place. Thus, intentionally routing high flows that occur in wet years through this system during the nonirrigation season would mimic what occurs naturally during the irrigation season. The second factor relates to the occurrence of atmospheric river storm events that deliver large amounts of precipitation to the region, much greater than average (Dettinger et al., 2015). With increased frequency and intensity projected under a warmer climate, such events have the potential to produce excess water over short periods of time that could be stored through mechanisms such as managed aquifer recharge (Niswonger et al., 2017). The third factor relates to the change in snowpack accumulation and shifts in snowmelt timing observed elsewhere in the Sierra Nevada (e.g., Godsey et al., 2014; Mote et al., 2018). Having a mechanism in place to maximize use of earlier snowmelt and shifts in streamflow timing could be advantageous and enhance regional groundwater sustainability. As part of the Water for the Seasons study, a hypothetical scenario was developed to determine the feasibility of managed aquifer recharge in the Carson River Basin, assuming no legal constraints. During “wet” or above-average water years, irrigators in the Upper Carson Valley would divert high flows and spread water over agricultural lands during the nonirrigation season. Assuming flows are abundant and “early,” diversions would begin prior to the growing season, when water would otherwise flow downstream to the Lahontan Reservoir. During “dry” years or drought periods, when surface water availability is less, irrigators in the Upper Carson Valley could augment surface water shortages with groundwater, allowing available surface water flows to flow downstream. Researchers hypothesize the amount of water has the potential to boost baseflow to support environmental instream flows, for example. What concerns have local water managers expressed? The hypothetical managed aquifer recharge scenario was presented to water managers in a workshop setting. Presentations included an overview of the hydrologic and operations modeling tools used to evaluate managed aquifer recharge by simulating the timing and distribution of water in the upper watershed. Specifically, in the Upper Carson Valley, a hydrologic model (GSFLOW) simulates streamflow driven by snowmelt, and surface and groundwater interactions, while a river basin operations model (MODSIM) allocates water according to the prior appropriation doctrine in the basin (see Figure 1) (Morway, Niswonger, & Triana, 2016; Niswonger et al., 2017). Integrating these two modeling tools advances the evaluation of climate impacts on water availability in agricultural communities and the resulting impacts of alternative management strategies (Morway et al., 2016). When asked about the viability of managed aquifer recharge, the perspectives of 11 managers varied (Figure 3). Regardless of rating, all managers questioned, “How would thisreally work?” Several managers questioned whether models could simulate the connectivity between surface and groundwater to accurately quantify changes to instream flow. Others raised concerns that managed aquifer recharge violates the Alpine Decree and Nevada Water Law. Still others requested researchers consider alternatives that could work within the confines of current (2019) water law. Managers posed specific questions that should be considered when evaluating the potential for managed aquifer recharge. For example: What triggers implementation of managed aquifer recharge?How “high” or “low” must annual flows be to initiate managed aquifer recharge? When in the water year is this determined? Where exactly in the Carson Valley is managed aquiferre charge possible? For example, what areas away from the floodplain could ensure long-term storage? Can model simulations quantify potential benefits and consequences system-wide?Would this information support decision-making, such as permitting of additional supplemental groundwater rights? How are researchers going to address managers’ research questions? Managers’ perspectives help to validate the on-the-ground potential of particular strategies and further refine alternative management scenarios. For example, understanding that managers are concerned with oversaturated fields helps researchers to define conditions in the model, such as what defines a wet versus “too” wet type of year and where to focus irrigation for managed aquifer recharge. Incorporating these nuances provides more accurate quantification of the potential benefits and consequences for users across the basin. Modeling is underway to simulate managed aquifer recharge scenarios and explore basin-wide implications. Researchers and local water managers will convene to collaboratively review results and further assess whether this or other strategies could work under the confines of existing water law. Subsequent fact sheets will present these findings.

Nevada