USGS Science⌕ Search

SEARCH · USGS Science

Results for “Nature Climate Change”

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

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

At least 1,711 records · Page 95Linked to original sources

Substantial hysteresis in emergent temperature sensitivity of global wetland CH4 emissions

Wetland methane (CH 4 ) emissions ( FCH 4 ) are important in global carbon budgets and climate change assessments. Currently, FCH 4 projections rely on prescribed static temperature sensitivity that varies among biogeochemical models. Meta-analyses have proposed a consistent FCH 4 temperature dependence across spatial scales for use in models; however, site-level studies demonstrate that FCH 4 are often controlled by factors beyond temperature. Here, we evaluate the relationship between FCH 4 and temperature using observations from the FLUXNET-CH 4 database. Measurements collected across the globe show substantial seasonal hysteresis between FCH 4 and temperature, suggesting larger FCH 4 sensitivity to temperature later in the frost-free season (about 77% of site-years). Results derived from a machine-learning model and several regression models highlight the importance of representing the large spatial and temporal variability within site-years and ecosystem types. Mechanistic advancements in biogeochemical model parameterization and detailed measurements in factors modulating CH 4 production are thus needed to improve global CH 4 budget assessments.

Nature Communications↗

Quantifying the coastal hazard risk reduction benefits of coral reef restoration in the U.S. Virgin Islands

Coastal habitat restoration, especially of coral reef ecosystems, can significantly reduce the exposure of coastal communities to natural hazards and, consequently, the risk of wave-driven flooding. Likewise, reef degradation can increase coastal flood risks to people and property. In this study, the valuation of coral reefs in the United States Virgin Islands (USVI), along the coasts of St. Croix, St. John, and St. Thomas, demonstrated the social and economic benefits provided by these natural defenses. Across the territory, more than 481 people and $31.2 million of infrastructure were estimated to receive protection from coral reefs per year (2010 U.S. dollars). In 2017, Hurricanes Irma and Maria significantly damaged coral reefs throughout the archipelago. By combining engineering, ecological, geospatial, social, and economic data and tools, this study provided a rigorous valuation of where potential coral reef restoration projects could help rebuild these damaged habitats and decrease the risks from coastal hazards faced by USVI’s reef-fronted communities. Multiple restoration scenarios were considered in the analysis, two of which are detailed in this report. These include (1) ‘Ecological’ restoration, where restoration creates a structure that is 0.25 m high and 25-m-wide reef, and (2) ‘Hybrid’ restoration, where restoration creates a structure that is 1.25 m high and 5 m wide. There are many ways that such structures could be developed. In the hydrodynamic analyses, there are no assumptions about how the restoration is developed. Many practitioners of both coral (and oyster reef) restoration consider that a reef height of 0.25 m might be delivered from planting corals alone and that 1.25 m might require a combination of artificial structures and coral planting. In a third scenario, the analysis investigated the reduction of protection benefits that would occur through the reduction of 1 meter of naturally occurring reef height due to reef degradation. The reduction of protection due to the loss of reefs can also be interpreted as the protection value of the existing reefs. In all studied restoration scenarios, it was assumed that the planting of corals would enhance hydrodynamic roughness, effectively dissipating incident wave energy and reducing the potential for coastal flooding. A standardized approach was employed to strategically locate potential restoration projects along the entire linear extent of existing reefs bordering the USVI, and to identify where coral reef restoration could offer valuable benefits in flood reduction. Potential restoration projects were only located within the existing distribution of reefs across the region, even though numerous sites were positioned far offshore (2-3 km), and some were at relatively deep depths (up to 7 m). Risk-based valuation approaches were followed to delineate flood zones at a 10 m 2 resolution along the entire region's reef-lined shorelines for all the potential coral reef restoration scenarios. These were subsequently compared to flood zones without coral reef restoration. The potential reduction in coastal flood risk provided by coral reef restoration, and the protection value of existing reefs, were quantified utilizing the latest information available at the time of analysis from the U.S. Census Bureau, Federal Emergency Management Agency (FEMA), and Bureau of Economic Analysis for return-interval storm events. The change in Expected Annual Damages (EAD), a metric indicating the annual protection gained due to coral reef restoration, was calculated based on the damages associated with each storm probability. The findings suggest that the benefits of reef restoration are spatially variable within the USVI. In some areas, the analysis showed limited benefits from reef restoration, which may be attributed to the depth or offshore distances of proposed restoration sites. However, there were a number of key areas where reef restoration could have substantial benefits for flood risk reduction. The annual flood risk reduction attributed to potential ‘ecological’ coral reef restoration in the USVI was 99 people and $6.1 million (2010 U.S. dollars). The Benefit-to-Cost Ratio (BCR) for this restoration approach was found to be larger than 1 (i.e., cost-effective) along 11% of the St. Croix coastline, 4.9% of the St. John coastline, and 8.7% of the St. Thomas coastline. This analysis offers stakeholders and decision-makers a spatially explicit and rigorous evaluation that illustrates how, where, and when potential coral reef restoration efforts in St. Croix, St. John, and St. Thomas could be instrumental to reducing coastal storm-induced flooding. Understanding areas where reef management, recovery, and restoration could effectively reduce climate hazard-related risks is crucial to protect and enhance the resilience of coastal communities in USVI.

Report↗

The evolution of mapping habitat for northern spotted owls (Strix occidentalis caurina): A comparison of photo-interpreted, Landsat-based, and lidar-based habitat maps

Wildlife habitat mapping has evolved at a rapid pace over the last few decades. Beginning with simple, often subjective, hand-drawn maps, habitat mapping now involves complex species distribution models (SDMs) using mapped predictor variables derived from remotely sensed data. For species that inhabit large geographic areas, remote sensing technology is often essential for producing range wide maps. Habitat monitoring for northern spotted owls ( Strix occidentalis caurina ), whose geographic covers about 23 million ha, is based on SDMs that use Landsat Thematic Mapper imagery to create forest vegetation data layers using gradient nearest neighbor (GNN) methods. Vegetation data layers derived from GNN are modeled relationships between forest inventory plot data, climate and topographic data, and the spectral signatures acquired by the satellite. When used as predictor variables for SDMs, there is some transference of the GNN modeling error to the final habitat map. Recent increases in the use of light detection and ranging (lidar) data, coupled with the need to produce spatially accurate and detailed forest vegetation maps have spurred interest in its use for SDMs and habitat mapping. Instead of modeling predictor variables from remotely sensed spectral data, lidar provides direct measurements of vegetation height for use in SDMs. We expect a SDM habitat map produced from directly measured predictor variables to be more accurate than one produced from modeled predictors. We used maximum entropy (Maxent) SDM modeling software to compare predictive performance and estimates of habitat area between Landsat-based and lidar-based northern spotted owl SDMs and habitat maps. We explored the differences and similarities between these maps, and to a pre-existing aerial photo-interpreted habitat map produced by local wildlife biologists. The lidar-based map had the highest predictive performance based on 10 bootstrapped replicate models (AUC = 0.809 ± 0.011), but the performance of the Landsat-based map was within acceptable limits (AUC = 0.717 ± 0.021). As is common with photo-interpreted maps, there was no accuracy assessment available for comparison. The photo-interpreted map produced the highest and lowest estimates of habitat area, depending on which habitat classes were included (nesting, roosting, and foraging habitat = 9962 ha, nesting habitat only = 6036 ha). The Landsat-based map produced an estimate of habitat area that was within this range (95% CI: 6679–9592 ha), while the lidar-based map produced an area estimate similar to what was interpreted by local wildlife biologists as nesting (i.e., high quality) habitat using aerial imagery (95% CI: 5453–7216). Confidence intervals of habitat area estimates from the SDMs based on Landsat and lidar overlapped. We concluded that both Landsat- and lidar-based SDMs produced reasonable maps and area estimates for northern spotted owl habitat within the study area. The lidar-based map was more precise and spatially similar to what local wildlife biologists considered spotted owl nesting habitat. The Landsat-based map provided a less precise spatial representation of habitat within the relatively small geographic confines of the study area, but habitat area estimates were similar to both the photo-interpreted and lidar-based maps. Photo-interpreted maps are time consuming to produce, subjective in nature, and difficult to replicate. SDMs provide a framework for efficiently producing habitat maps that can be replicated as habitat conditions change over time, provided that comparable remotely sensed data are available. When the SDM uses predictor variables extracted from lidar data, it can produce a habitat map that is both accurate and useful at large and small spatial scales. In comparison, SDMs using Landsat-based data are more appropriate for large scale analyses of amounts and general spatial patterns of habitat at regional scales.

Remote Sensing of Environment↗

Regional and other general factors bearing on evaluation of earthquake and other geologic hazards to coastal communities of southeastern Alaska

The great Alaska earthquake of March 27, 1964, brought into sharp focus the need for engineering geologic studies in seismically active regions. As a result, nine communities in southeastern Alaska were selected for reconnaissance investigations as an integral part of an overall program to evaluate earthquake and other geologic hazards in most of the larger Alaska coastal communities. This report gives background information on the regional and other general factors that bear on these evaluations. Southeastern Alaska, about 525 miles long and averaging about 125 miles in width, consists of a narrow mainland strip and numerous islands. For the most part, it is a region of rugged relief with numerous glaciers capping many of the higher mountainous areas and with long linear fiords forming the inland waterways. A maritime climate prevails with mild winters and cool summers. The southeastern part of the region receives the highest precipitation in the continental United States. Ketchikan, with a population of 6,994 in 1970, is the largest city. Geology and structure of the area are complex. Igneous, metamorphic, and sedimentary rocks crop out and range in age from Paleozoic to Tertiary. Surficial deposits of Pleistocene and Holocene age mantle many areas. All of southeastern Alaska, except probably the highest peaks, was covered by glacier ice advances of late Pleistocene age. Major deglaciation was well advanced by 10,000 years ago--a time which approximately marks the end of the Pleistocene and the beginning of the Holocene. There followed a period of warm climate called the Hypsithermal, which in southeastern Alaska began 7,000-8,000 years ago and ended about 4,800-3,500 years ago. Glaciers in most places receded back of their present positions. The Hypsithermal was followed by an interval (termed Neoglaciation) of cooler climate and resurgence of glacier ice which continues to the present, although most glaciers are now rapidly receding. During the past 10,000 years worldwide sea level has risen about 100 feet, but during the past 4,000 years it has risen only about 10 feet or about 0.03 inch per year. With sea level used as a datum, the amount of sea-level rise must be added to the apparent uplift of land for the time under consideration to determine the actual amount of land uplift. The widespread presence of emergent marine deposits, several hundred feet above sea level, demonstrates that the land in southeastern Alaska has been uplifted since the last major deglaciation. The greatest known has been uplifted since the last major deglaciation. The greatest known uplift is in the vicinity of Juneau where glaciomarine deposits are present 750 feet above present sea level. Part of southeastern Alaska is presently undergoing one of the most rapid rates of uplift of any place in the world. The fastest emergence is occurring in the Glacier Bay area where the land is being uplifted relative to sea level approximately 3.9 cm per year. Most or all of the uplift appears to be due to rebound as a result of deglaciation. Southeastern Alaska lies within the circum-Pacific earthquake belt, one of the world's greatest zones of seismic activity. During historic time, there have been five earthquakes in the region with magnitudes of 8 or greater, three with magnitudes of 7 to 8, eight with magnitudes of 6 to 7, more than 15 with magnitudes of 5 to 6, and about 140 recorded earthquakes with magnitudes smaller than 5 or of unassigned magnitudes. All of the earthquakes with magnitudes 8 or greater, and a large proportion of the others, appear to be related to the active Fairweather- Queen Charlotte Islands fault system or its western extension, the Chugach-St. Elias fault. Earthquake epicenters on the Denali fault system, the other major fault system in southeastern Alaska, are few in comparison. However, because high microearthquake activity has been recorded recently on this system and earthquakes of moderate size have occurred on some of its segments, the Denali fault system probably should not be dismissed as a relict fault system of no current tectonic importance. There are numerous other known faults, as well as lineaments that may be faults of varying degrees of tectonic activity in southeastern Alaska, adjacent Canada, and eastern Alaska. One of these elements is the Totschunda fault system, which connects with the Denali fault system in eastern Alaska; it has been very active during Holocene time but few historical earthquake epicenters appear to be related to it. Both historical seismicity and geologic conditions, such as frequency and recency of faulting, must be considered together to permit an assessment of the future earthquake probability of an area. Data are too few for both factors for an accurate evaluation to be made of earthquake probability in southeastern Alaska. However, information compiled in the form of strain-release and seismic-zone maps permit some generalizations. Thus, it is tentatively concluded that most, if not all, of southeastern Alaska should be placed in seismic zone 3, a zone in which earthquakes of magnitude greater than 6 will occur from time to time and where there may be major damage to manmade structures. Inferred effects from future earthquakes in southeastern Alaska include: (1) surface displacement along faults and other tectonic land-level changes, (2) ground shaking, (3) compaction, (4) liquefaction in cohesionless materials, (5) reaction of sensitive and quick clays, (6) water-sediment ejection and associated subsidence and ground fracturing, (7) earthquake-induced sub aerial slides and slumps, (8) earthquake induced subaqueous slides, (9) effects on glaciers and related features, (10) effects on ground water and stream flow, and (11) tsunamis, seiches, and other abnormal water waves. Because of the reconnaissance nature of our studies in the coastal communities and the sparsity of laboratory data on physical properties of geologic units in each area studied, the inferred effects must be largely empirical and generalized. Therefore, the inferences are based in large part upon the effects of past major earthquakes in Alaska and elsewhere, particularly upon the well-documented effects of the Alaska earthquake of March 27, 1964. Buildings, highways, bridges, tunnels, harbor facilities, pipelines, canals, and other manmade structures may be severely damaged or destroyed by fault displacement or related tectonic land-level changes in southeastern Alaska. Direct damage from fault rupture would be restricted virtually to structures built directly athwart the fault. In California and Nevada, fault rupture almost always accompanies shocks of magnitude 6.5 or greater. The Alaska earthquake of March 27, 1964, and the Chilean earthquake of May 22, 1960, dramatically illustrated the severe adverse effects that can result from uplift or subsidence over a wide area. The variable most responsible for the degree of shaking at any epicentral distance is the type of ground. Generally, shaking is considerably greater in poorly consolidated deposits than in hard bedrock, particularly if the deposits are water saturated. Severe shaking of alluvial deposits and manmade fill, with resultant heavy damage, is well documented from the records of many past earthquakes. Damage commonly has been heavy as a result of ground settlement caused by compaction of loose sediments by shaking during an earthquake. This has been especially true where compaction was accompanied by tectonic downdrop of land, such as occurred during the Chilean earthquake of 1960 and the Alaska earthquake of 1964. Loosely emplaced manmade fill, deltaic deposits, beach deposits, and alluvial deposits may be susceptible to compaction in southeastern Alaska during a severe earthquake. Liquefaction of sand and silt is a fairly common effect of large earthquakes. It was well illustrated at Niigata, Japan, during the earthquake of June 16, 1964, and resulted in extensive damage. When part of a sloping soil mass liquefies, the entire mass can undergo catastrophic failure and can flow as a high-density liquid. In southeastern Alaska, deltaic deposits probably would be most susceptible to liquefaction. Sensitive and quick clays, which lose a considerable part of their strength when shaken, commonly fail during an earthquake and become rapid earthflows. Extensive studies were made of the sensitivity of the Bootlegger Cove Clay at Anchorage because of the marked loss of shear strength and dramatic failures of the deposits during the Alaska earthquake of 1964. If similar sensitive clays are present in some places in southeastern Alaska, they most likely are in some of the emergent fine-grained marine deposits; supporting data to confirm their presence, however, are largely lacking. Records of some 50 major earthquakes show that in at least half of the instances water and sediment have been ejected from surficial deposits Water-sediment ejection and associated subsidence and ground fracturing commonly cause extensive damage to the works of man. Ejecta may fill basements and other low-lying parts of buildings. Agricultural land can be covered with a blanket of infertile soils, and small ponds can be filled or made shallow. In southeastern Alaska these phenomena are most likely to occur on valley floors, deltas, tidal flats, alluvial fans, swamps, and lakeshores. Earthquake-induced sliding on land generally is confined to steep slopes but may take place in fine-grained deposits on moderately to nearly flat surfaces if the deposits are subject to liquefaction. A large rockslide triggered by the Lituya Bay, Alaska, earthquake of July 10, 1958, generated a wave that surged up the opposite wall of the inlet to a record height of 1,740 feet. During the Hebgen Lake, Montana, earthquake of August 17, 1959, a spectacular rockslide plunged into the Madison River canyon, buried 28 people, dammed the river, and created a large lake. Earthquake-records are replete with accounts of sliding of surficial deposits during moderate to large earthquakes. Most or all of the general factors that favor subaerial landsliding are present in southeastern Alaska. Earthquake-induced subaqueous slides can produce adverse effects both nearshore and some distance offshore. Nearshore sliding may progress shoreward and destroy harbor facilities and other structures, commonly with substantial loss of life. Disastrous large submarine slides occurred along the fronts of deltas in Seward and Valdez during the Alaska earthquake of 1964. In similar fashion, the largest submarine slides in southeastern Alaska likely will be triggered along the larger delta fronts. Sliding farther offshore can constitute a threat to navigation because of changes in water depths. Also underwater sliding can break communication cables. Glaciers were not greatly affected by the Alaska earthquake of 1964 despite the fact that about 20 percent of the area that underwent strong shaking is covered by ice. In contrast, the cataclysmic avalanche of ice and rock that fell from a high glacier-covered peak in Peru during the earthquake of May 31, 1970, produced devastating effects downvalley on man and his works in the form of mudflows. Most towns in southeastern Alaska are sufficiently distant from glaciers so as not be to directly affected. Both the Alaska earthquake of 1964 and the Hebgen Lake, Montana, earthquake of 1959 significantly affected ground- and surface-water regimens. Water levels in some wells declined whereas in others flow increased. Some springs discharged at a rate three times as much as normal; flow of others decreased or stopped. Discharge of many streams increased markedly. Most or all of the effects described above could occur in parts of southeastern Alaska during future large earthquakes. Tsunamis, seiches, and other abnormal water waves associated with large earthquakes commonly cause vast property damage and heavy loss of life. Tsunami effects can be devastating to coastal areas as far as many thousands of miles from their generation source. Seiche effects generally are confined to inland bodies of water or to relatively enclosed coastal bodies of water. Abnormal waves generated by submarine sliding or by subaerial sliding into water generally produce only local effects but may be highly devastating. Tsunami waves resulting from the Chilean earthquake of 1960 inflicted extensive damage and loss of life on coastal communities throughout a large part of southern Chile, and significant runups and damage were recorded in many places throughout the Pacific Ocean area. The tsunami waves generated by the Alaska earthquake of 1964 struck with devastating force along a broad stretch of the Alaska coast and produced heavy property damage and loss of life as far away as Crescent City, Calif. Seiche waves generated by that earthquake reached runup heights of 20-30 feet on some lakes in Alaska, and water-level fluctuations were recorded on streams, reservoirs, lakes, and swimming pools in States bordering the Gulf of Mexico. Waves generated by submarine sliding struck violently at a number of places during or immediately after the quake and were the major cause of loss of life and damage to property. Slide-generated waves probably would have a higher destructive potential in southeastern Alaska than either tsunami waves or seiche waves because of their possibly higher local runups and because they can hit the shores almost without warning during or immediately after an earthquake. Nonearthquake-related geologic hazards, although generally far less dramatic than those related to earthquakes, tend to occur so much more frequently or persistently that their aggregate effects can be significant. Three kinds of geologic hazards of this type are discussed: (1) nonearthquake-induced landsliding and subaqueous sliding, (2) flooding, and (3) land uplift. The potential for nonearthquake-triggered landsliding in southeastern Alaska ranges widely from place to place. Past sliding generally furnishes the clue in the prediction of where and in what materials future sliding will occur. Fast-moving rockslides, debris slides, and mudflows can be expected to occur from time to time on steep slopes and be highly destructive to highways, power plants, pipelines, buildings, and other facilities located on a slope or at its base. Present slow downslope movement of talus can be expected to continue at the same general rate unless conditions are changed by man or there are climatic changes. Snow and debris avalanches can be especially hazardous during winter months. Long-inactive landslides may be triggered into renewed activity or new slides may be created by man-induced modifications. Accelerated slope erosion and debris flows may follow large-scale clearing and cutting of timber. Subaqueous sliding can be expected to occur periodically along fronts of deltas and on other oversteepened underwater slopes. Floods have been common in parts of southeastern Alaska because of heavy precipitation and rapid runoff from steep slopes with resulting heavy damage to roads and other facilities. Continued damage can be expected in the future unless more remedial measures are taken. Current uplift of land in southeastern Alaska, although probably not affecting man significantly in a short period of time, may have some adverse long-term effects. These long-term effects should be borne in mind when facilities such as docks and boat harbors are constructed on or near the shore, where there is a critical relation between height of land and water.

Alaska↗

Assessing spatial variability of nutrients, phytoplankton, and related water-quality constituents in the California Sacramento–San Joaquin Delta at the landscape scale—2018 high resolution mapping surveys

Executive Summary This study examined the abundance and distribution of nutrients and phytoplankton in the tidal aquatic environments of the Sacramento–San Joaquin Delta (Delta) and Suisun Bay, comprising three spatial surveys conducted in May, July, and October of 2018 that used continuous underway high frequency sampling and measurements onboard a high-speed boat to characterize spatial variation across the extent of the Delta. The method used involves simultaneously collecting information about the concentration and spatial distribution of all major nutrient forms with analogous information about the major classes of phytoplankton and associated water-quality conditions. The results showed substantial variation across space and time, providing an unprecedented snapshot of the dynamic environmental processes that shape the ways nutrients interact with and affect aquatic habitats in the Delta. The purposes of this study were to improve our understanding of how hydrodynamics, landscape features, and aquatic primary productivity interact to drive nutrient cycling and transport in the Delta and to provide insights into the underlying processes most directly responsible for the conditions at the time of this study, and thus into the range of conditions that may be expected following the wide array of prospective future changes to the Delta. One major anticipated change at the time of this study was the planned upgrade to the Sacramento Regional Wastewater Treatment Plant, but the study also informs our understanding of potential effects from other changes to the Delta, such as those caused by other nutrient-management actions, flow actions, large-scale wetland restoration, drought, flood, levee failure, and changes to water management. Nutrient loading is the primary driver of nutrient concentrations in the Delta, but several other major drivers interact to shape their distribution and effects: geomorphology, hydrodynamics, landscape features, and aquatic productivity. Hydrodynamics affect timescales of transport and dilution of nutrient loads in the Delta. During transit through the system, channel geometry, tidal mixing, and water exports affect hydrodynamics in diverse ways that influence water-residence and transport times, thereby markedly affecting the range of times during which natural internal cycling can alter nutrient concentrations and forms. Channel geometry and location shape tidal energy and river currents into these observed dynamics. Interactions with Delta aquatic landscapes such as herbaceous tidal marsh, submerged aquatic vegetation, and large expanses of intertidal or subtidal sediments (all highly productive landscapes) exert demand on available nutrient supplies but can also simultaneously transform and generate nutrients. Finally, while phytoplankton require nutrients to sustain production and thus are a potential nutrient sink, the amount and form of nutrients also can influence the occurrence of harmful algal blooms (HABs) that adversely affect aquatic organisms as well as affect the occurrence of beneficial algal blooms that result in production of algae that are favorable for imperiled Delta pelagic aquatic food webs. The surveys revealed a complex mosaic of spatial variation, with nutrient concentrations varying from near zero to well above concentrations considered eutrophic; nutrient concentrations were more often related to the extent of hydrologic transport and mixing than to specific geographic locations or to specific landscape features. Similarly, the surveys identified phytoplankton abundance ranging from near detection to the level of large phytoplankton blooms, with large variation in phytoplankton community composition. Although the study occurred during a period of low bloom activity, phytoplankton productivity appeared to be the strongest potential sink for inorganic nutrients in the Delta, indicating that it is a larger control on nutrient concentrations and distribution than previously understood. Cycling and transformation within the water column only appeared to substantially lower total nutrient concentrations at the longest estimated transport timescales. Contrary to expectations, we did not observe substantial nutrient depletion near landscape-scale features such as open-water habitats, submerged aquatic vegetation beds, extensive wetlands, or exposed sediments, indicating that these habitat types did not act as major sinks for nutrients in the Delta during these surveys. These results indicated that nutrient reduction efforts may have the greatest effect on pelagic phytoplankton productivity in the more productive reaches of the Delta and estuary, but these effects are unlikely to be magnified by changes to nutrient loss within the Delta over conceivable changes in flow conditions, Delta water management actions, or large-scale wetland restoration activities. Nevertheless, local processes were shown to cause substantial loss, and thus integrating of nutrient effects with other indicators of aquatic habitat conditions will help inform planning future actions at specific sites. Finally, we note that the primary contribution of this study was intended to be the survey data themselves. Aside from the results highlighted in this report, the surveys are a benchmark against which future environmental change may be evaluated, including changes to nutrient management or water exports, drought, large-scale wetland restoration, and climate change. Further, although we highlight some of the main findings from the surveys in this report, the necessarily limited scope precludes examination of many topics for which these surveys may be highly informative. To facilitate the utility of these data to stakeholders, managers, and researchers, we have released the data online (Bergamaschi and others, 2020) and created an online data exploration portal ( https:​//ca.water​.usgs.gov/​bay-​delta/​2018-​delta-​wide-​mapping-​surveys.html ) where users may query the surveys in a variety of ways to test hypotheses, examine relationships, assess spatial trends, and download data. The data exploration portal is intended to be an immersive experience that allows users to gain greater understanding of the complex interactions that shape Delta aquatic environments. This report is intended as a companion to the portal, allowing the reader to challenge and further explore the highlighted findings. This study was a collaboration between the U.S. Geological Survey and the Delta Regional Monitoring Program, with additional funding provided from U.S. Geological Survey Cooperative Matching Funds Program.

California↗

Fishway Entrance Palisade

This technical report summarizes the work that was conducted by the University of Massachusetts Amherst and the United States Geological Survey (USGS), along with other project partners, on the Fishway Entrance Palisade (EP), a projected funded through the Department of Energy’s (DOE) funding opportunity titled ‘Innovative Solutions for Fish Passage at Hydropower Dams’ (DE‐FOA‐0001662). The period of performance ranged from September 1, 2018 through September 30, 2021. The EP is a novel fish passage engineering technology designed to provide more favorable entry conditions for fish and to reduce costs relative to conventional fishway auxiliary water systems (AWS). The EP project has four primary components. First, the Northeast United States Auxiliary Water Systems Database was created (Northeast Fishway Auxiliary Water Systems Database Section). The database, developed with material provided by the U.S. Fish and Wildlife Service, contains information on fishway type (e.g., lift, Denil, pool and weir) and Auxiliary Water System (AWS) details (e.g., water conveyance method, diffuser type) for 60 hydroelectric sites in the region. Findings indicate that nearly 4 out of every 10 fishway in the region is a fish lift and approximately 1 out of every 4 is a Denil ladder. The remainder are a mix of vertical slot fishways, pool and weirs, and Ice Harbor fishways. Furthermore, over half of all AWS systems use floor diffusers to discharge the auxiliary (or attraction) water into the entrance of a fishway, whereas only 14% use wall diffusers. Second, limited experiments on a conventional AWS with live, actively migrating fish were conducted at the USGS Easter Ecological Science Center (EESC) S.O. Conte Research Laboratory (Conventional Auxiliary Water System Experiments Section). This study determined how water velocity through a wall diffuser, without turning vanes or timber baffles to distribute the flow, affects the behavior and passage of adult American shad, a conservative surrogate species for migratory fish on the East Coast. Two gross diffuser velocity treatments were examined, 0.5 ft/s and 1.0 ft/s. These wall diffuser velocities represented current (0.5 ft/s) and past (1.0 ft/s) design criteria guidelines set forth by the USFWS North Atlantic-Appalachian Region (Rojas 2020; USFWS 2019). Six trials with a total of 151 American Shad were conducted in June of 2019 for the two treatments. No differences in American shad passage efficiency were discovered between the two treatments, while approximately 3 in every 4 attempts were successful at passing the diffuser. While these results may appear to indicate that the generally accepted gross wall diffuser velocity criteria for American shad of 0.5 ft/s could be safely increased to 1.0 ft/s, further analysis is warranted. Furthermore, it is unknown how other migratory and resident fish species that traverse these structures would be impacted by such a change. Studying the wall diffuser hydraulics led to an important AWS observation. Without turning vanes or timber baffles in this study, doubling the diffuser area was insufficient at producing the type of flow field change one may expect by halving the gross diffuser velocity. Instead, the flow fields throughout each treatments study area were similar, which led to similar results in shad performance. This not only highlights the importance of installing flow guidance devices like turning vanes, but also to the importance of properly maintaining them, which can be costly. Third, more expansive experiments on the novel EP were conducted in the spring of 2019 and 2021 (Fishway Entrance Palisade Experiments). The goal of this study was to determine how adult American shad responded to a variety of conditions at a full-scale EP. A total of six treatments were examined by changing the average auxiliary channel velocity between 1.0 and 5.0 ft/s in intervals of 1.0 ft/s and by inserting/removing an entrance gate at the opening of the fishway. Thirty trials with a total of 1,273 shad were conducted over the two years. In all treatments, at least ~7 out of every 10 fish successfully passed the EP diffuser and swam into the entrance channel within the 3.5-hour long trial, highlighting the general effectiveness of the novel AWS technology. In both study years, lower velocities through the EP diffuser led to increased shad performance, though performance peaked for the 2 ft/s velocity treatment. This treatment condition represents an approximate six-fold increase in gross diffuser velocity relative to conventional auxiliary water systems, which in turn presents opportunities for cost savings (e.g., reduction in diffuser size). Shad performance, in general, was worse in 2019 than in 2021, potentially due to the different run timing when our trials were conducted (2019 trials occurred near the end of the migration season, unlike in 2021). Treatments in 2019 had approximately a 20% reduction in entrance efficiency by the trial end, including a 16.7% drop for the 3 ft/s velocity treatment in 2019 relative to 2021 (the only carryover treatment between years). Lastly, adding an entrance gate caused a significant delay to entry. The time to 25% entry raised ~20 minutes from the near instantaneous 25% entry that was reported for the other treatments conducted in the same year (2021). Though by the end of the 3.5-hour trial, the overall entrance efficiency nearly matched those of the other 2021 treatments. The fourth and final component of the EP project was an economic analysis that focused on the cost of attraction and environmental flows (Modeling Power Generation Losses Due to Environmental and Fish Passage Attraction Flows at a Run-Of-River Hydroelectric Operation in the Northeast). The study assessed the economic impact of meeting environmental flow requirements at a representative hydroelectric facility and fish lift in the Northeast. An initial finding of the study was that there is a paucity of published data on the costs of meeting attraction and environmental flows. This is due, in part, to the proprietary nature of this data. To explore the costs associated with these flows, three types of environmental flows were assessed: upstream fishway attraction flows, downstream fishway attraction flows, and habitat maintenance minimum flows. A physics-based model was developed and calibrated with three years of hourly generation and flow data as inputs. Gage flow inputs were adjusted and used to calculate power generated. To address hydrologic variability, the model was executed to simulate 30 years of historical flows. Results indicate that both interannual and seasonal climatic factors impact the costs of meeting environmental flow requirements. Generation potential is most strongly curtailed during dry years in terms of maximizing the capacity factor (the percent of time a plant generates at capacity). Dry years, and especially dry summers, have the most significant costs associated with mitigation flows. Of the three types of flows, habitat flows are most costly in terms of power production, followed by upstream attraction flows. Downstream attraction flows are least costly. This finding is the likely result of differences in both flow rates and duration of the seasonal requirement for each flow. Overall, environmental flows represented a 2-12% loss in annual generation, but losses during a dry summer can reach over 20%.

Final Technical Report↗

Comparison of water consumption in two riparian vegetation communities along the central Platte River, Nebraska, 2008–09 and 2011

The Platte River is a vital natural resource for the people, plants, and animals of Nebraska. A recent study quantified water use by riparian woodlands along central reaches of the Platte River, Nebraska, finding that water use was mainly regulated below maximum predicted levels. A comparative study was launched through a cooperative partnership between the U.S. Geological Survey, the Central Platte Natural Resources District, the Nebraska Department of Natural Resources, and the Nebraska Environmental Trust to compare water use of a riparian woodland with that of a grazed riparian grassland along the central Platte River. This report describes the results of the 3-year study by the U.S. Geological Survey to measure the evapotranspiration (ET) rates in the two riparian vegetation communities. Evapotranspiration was measured during 2008–09 and 2011 using the eddy-covariance method at a riparian woodland near Odessa, hereinafter referred to as the “woodland site,” and a riparian grassland pasture near Elm Creek, hereinafter referred to as the “grassland site.” Overall, annual ET totals at the grassland site were 90 percent of the annual ET measured at the woodland site, with averages of 653 millimeters (mm) and 726 mm, respectively. Evapotranspiration rates were similar at the grassland site and the woodland site during the spring and fall seasons, but at the woodland site ET rates were higher than those of the grassland site during the peak-growth summer months of June through August. These seasonal differences and the slightly lower ET rates at the grassland site were likely the result of differing plant communities, disturbance effects related to grazing and flooding, and climatic differences between the sites. The annual water balance was calculated for each site and indicated that the predominant factors in the water balance at both sites were ET and precipitation. Annual precipitation for the study period ranged from near to above the normal precipitation of 640 mm. Substantial precipitation fell in May and October 2008 that caused flooding along the Platte River in May of this especially wet year. There was a deficit in precipitation compared to ET at both sites in 2009 and 2011, leading to a net groundwater use of greater than 140 mm per year at the woodland site and greater than 55 mm per year at the grassland site. This indicates that the net annual groundwater use or recharge depends predominately upon the relation between ET and precipitation in these riparian areas with shallow soil layers above the groundwater table. Prior research at the woodland site provided four additional annual water balances dating back to 2002 for comparison with the study period at the woodland site. Perhaps most striking in this comparison was the 25-percent increase in annual ET for 2008–09 and 2011 despite precipitation totals and potential ET rates that were within the range of those measured in 2002–05. As a result, the water balance indicates that groundwater was discharged 2 of the 3 years of the study. This likely was caused by higher groundwater levels and a healthier plant community in 2008–09 and 2011 relative to the drought-affected years of 2002–05. As a result of these changes, the crop coefficients developed for riparian woodlands during the prior research underestimated 2008–09 and 2011 annual ET rates by an average of 35 percent. Though new crop coefficients were developed by this study, the importance of soil-moisture stress and plant community successional dynamics need to be considered when applying these coefficients at other riparian sites or into the future. Nonetheless, their development and the data on which they are based may provide improved understanding of water consumption by riparian grasslands and riparian woodlands along the central Platte River.

Nebraska↗

Hydrogeologic setting and preliminary estimates of hydrologic components for Bull Run Lake and the Bull Run Lake drainage basin, Multnomah and Clackamas counties, Oregon

The hydrogeologic setting was described and preliminary estimates of hydrologic components prepared for the Bull Run Lake and for the Bull Run Lake drainage basin, in the Cascade Range of northwestern Oregon. The 0.73-square-mile lake and the 3.44-square-mile drainage basin lie within the Bull Run Watershed, the principal water supply for the Portland, Oregon, metropolitan area. During periods of high demand or low inflows to the watershed, the City of Portland, Bureau of Water Works, releases water from Bull Run Lake to augment the supply. Bull Run Lake is impounded by a natural dam formed by a landslide. Outflow of ground water from the lake through the landslide emerges as springflow at the toe of the landslide and forms the headwaters of the Bull Run River. The approximately 4,300-Mgal (million gallons) discharge of the Bull Run River measured below the springs during the 1993 water year is composed of (1) outflow of ground water from Bull Run Lake through the landslide (approximately 60 percent), (2) ground water originating from the contributing drainage area between the lake and the springs (approximately 34 percent), (3) streamflow from Bull Run Lake (approximately 5 percent), and (4) surface runoff (streamflow and overland flow) from the contributing drainage area between the lake and the springs (approximately 1 percent). Estimated ranges for inflows to the Bull Run Lake drainage basin during the 1993 water year were about 3,400 to 9,200 Mgal from precipitation from rain and snow, and about 0 to 3,300 Mgal from fog drip. Estimated ranges for outflows from the lake basin, listed from largest to smallest, were about 1,800 to 3,400 Mgal for ground-water outflow through the landslide; about 600 to 1,800 Mgal for evapotranspiration from the land surface; about 170 to 410 Mgal for lake evaporation; and about 0 to 400 Mgal for streamflow from the lake. Ground- water outflow through the consolidated rocks could not be evaluated owing to the lack of data. The lake storage increased by a range of from about 1,700 to 1,900 Mgal. Changes in ground-water storage and soil-moisture storage could not be evaluated as a result of insufficient data. Estimated inflows to Bull Run Lake from precipitation on the lake surface during the 1993 water year ranged from about 600 to 1,600 Mgal. Inflows from ground water and surface runoff could not be evaluated owing to the lack of data. Estimated ranges for outflows from the lake were about 1,800 to 3,400 Mgal from ground-wateroutflow through the landslide, about 170 to 410 Mgal from lake evaporation, and about 0 to 400 Mgal from streamflow. Outflow of ground water through the consolidated rocks could not be evaluated owing to the lack of data. Lake storage increased by a range of from about 1,700 to 1,900 Mgal. Suggestions for further study include (1) evaluation of the surface-runoff component of inflow to the lake; (2) use of a cross-sectional ground-water flow model to estimate ground-water inflow, outflow, and storage; (3) additional data collection to reduce the uncertainties of the hydrologic components that have large relative uncertainties; and (4) determination of long-term trends for a wide range of climatic and hydrologic conditions.

Oregon↗

Mars Global Surveyor Thermal Emission Spectrometer experiment: Investigation description and surface science results

The Thermal Emission Spectrometer (TES) investigation on Mars Global Surveyor (MGS) is aimed at determining (1) the composition of surface minerals, rocks, and ices; (2) the temperature and dynamics of the atmosphere; (3) the properties of the atmospheric aerosols and clouds; (4) the nature of the polar regions; and (5) the thermophysical properties of the surface materials. These objectives are met using an infrared (5.8- to 50-μm) interferometric spectrometer, along with broadband thermal (5.1- to 150-μm) and visible/near-IR (0.3- to 2.9-μm) radiometers. The MGS TES instrument weighs 14.47 kg, consumes 10.6 W when operating, and is 23.6×35.5×40.0 cm in size. The TES data are calibrated to a 1-σ precision of 2.5 −6 ×10 −8 W cm −2 sr −1 /cm −1 , 1.6×10 −6 W cm −2 sr −1 , and ∼0.5 K in the spectrometer, visible/near-IR bolometer, and IR bolometer, respectively. These instrument subsections are calibrated to an absolute accuracy of ∼4×10 −8 W cm −2 sr −1 /cm −1 (0.5 K at 280 K), 1–2%, and ∼1–2 K, respectively. Global mapping of surface mineralogy at a spatial resolution of 3 km has shown the following: (1) The mineralogic composition of dark regions varies from basaltic, primarily plagioclase feldspar and clinopyroxene, in the ancient, southern highlands to andesitic, dominated by plagioclase feldspar and volcanic glass, in the younger northern plains. (2) Aqueous mineralization has produced gray, crystalline hematite in limited regions under ambient or hydrothermal conditions; these deposits are interpreted to be in-place sedimentary rock formations and indicate that liquid water was stable near the surface for a long period of time. (3) There is no evidence for large-scale (tens of kilometers) occurrences of moderate-grained (>50-μm) carbonates exposed at the surface at a detection limit of ∼10%. (4) Unweathered volcanic minerals dominate the spectral properties of dark regions, and weathering products, such as clays, have not been observed anywhere above a detection limit of ∼10%; this lack of evidence for chemical weathering indicates a geologic history dominated by a cold, dry climate in which mechanical, rather than chemical, weathering was the significant form of erosion and sediment production. (5) There is no conclusive evidence for sulfate minerals at a detection limit of ∼15%. The polar region has been studied with the following major conclusions: (1) Condensed CO 2 has three distinct end-members, from fine-grained crystals to slab ice. (2) The growth and retreat of the polar caps observed by MGS is virtually the same as observed by Viking 12 Martian years ago. (3) Unique regions have been identified that appear to differ primarily in the grain size of CO 2 ; one south polar region appears to remain as black slab CO 2 ice throughout its sublimation. (4) Regional atmospheric dust is common in localized and regional dust storms around the margin and interior of the southern cap. Analysis of the thermophysical properties of the surface shows that (1) the spatial pattern of albedo has changed since Viking observations, (2) a unique cluster of surface materials with intermediate inertia and albedo occurs that is distinct from the previously identified low-inertia/bright and high-inertia/dark surfaces, and (3) localized patches of high-inertia material have been found in topographic lows and may have been formed by a unique set of aeolian, fluvial, or erosional processes or may be exposed bedrock.

Journal of Geophysical Research E: Planets↗

Geology and hydrogeology of the Caribbean Islands aquifer system of the Commonwealth of Puerto Rico and the U.S. Virgin Islands

Poorly lithified to unconsolidated carbonate and clastic sedimentary rocks of Tertiary (Oligocene to Pliocene) and Quaternary (Pleistocene to Holocene) age compose the South Coast aquifer and the North Coast limestone aquifer system of Puerto Rico; poorly lithified to unlithified carbonate rocks of late Tertiary (early Miocene to Pliocene) age make up the Kingshill aquifer of St. Croix, U.S. Virgin Islands. The South Coast aquifer, North Coast limestone aquifer system, and Kingshill aquifer are the most areally extensive and function as the major sources of ground water in the U.S. Caribbean Islands Regional Aquifer-System Analysis (CI-RASA) study area. In Puerto Rico's South Coast ground-water province, more than 1,000 meters of clastic and carbonate rocks of Oligocene to Pliocene age infill the South Coast Tertiary Basin. The pattern of lithofacies within this basin appears to have been controlled by changes in base level that were, at times, dominated by tectonic movement (uplift and subsidence), but were also influenced by eustasy. Deposition of the 70-kilometer long and 3- to 8-kilometer wide fan-delta plain that covers much of the South Coast ground-water province occurred largely in response to glacially-induced changes in sea level and climate during the Quaternary period. Tectonic movement played a much less important role during the Quaternary. The North Coast ground-water province of Puerto Rico is underlain by homoclinal coastal plain wedge of carbonate and siliciclastic rocks that infill the North Coast Tertiary Basin and thicken to more than 1,700 meters. A thin basal siliciclastic sequence of late Oligocene age is overlain by a thick section of mostly carbonate rocks of Oligocene to middle Miocene age. Globigerinid limestone of late Miocene to Pliocene age crops out and lies in the shallow subsurface areas of northwestern Puerto Rico. Oligocene to middle Miocene age rocks tentatively can be divided into five depositional sequences and associated systems tracts; these rocks record carbonate and minor siliciclastic deposition that occurred in response to changes in relative sea level. The Cibao Formation represents the most complex of these sequences and contains a varied facies of carbonate, mixed carbonate-siliciclastic, and siliciclastic rocks that reflect differential uplift, subsidence, and transgression of the sea. Uplift, graben formation, and gradual shallowing of the sea are reflected within the bathyal-dominated sedimentary facies of the Kingshill Limestone in St. Croix, U.S. Virgin Islands. Reef-tract limestone beds of Pliocene age were subject to exposure, resubmergence, and meteoric leaching of aragonitic skeletal debris; these beds contain patchy lenses of dolomite that are restricted to a small, structurally-controlled embayment. The South Coast aquifer, the principal water-bearing unit of Puerto Rico's South Coast ground-water province, consists of boulder- to silt-size detritus formed by large and small coalescing fan deltas of Pleistocene to Holocene age. Deep well data indicates that it is possible to vertically separate and group a highly complex and irregular-bedded detrital sequence that underlies distal parts of the fan-delta plain into discrete water-bearing units if correlated with 30- to 40-meter thick, eustatically-controlled depositional cycles. Lithofacies maps show that greatest hydraulic conductivity within the fan-delta plain is generally associated with proximal fan and midfan areas. Distal and interfan areas are least permeable. Alluvial valley aquifers located in the western part of the South Coast ground-water province are important local sources of water supply and appear to contain some of the same physical and hydraulic characteristics as the South Coast aquifer. Older sedimentary rocks within the basin are poor aquifers; conglomeratic beds are well-cemented, and carbonate beds do not contain well-developed solution features, except locally where the beds are overlain by alluvium. Ground-water occurs under unconfined conditions in proximal and midfan areas. Confined conditions within deeper parts of the system and in interfan and some midfan areas are created largely by the intercalated nature of discontinuous fine-grained beds that retard vertical ground-water movement. The development of water resources in southern Puerto Rico has modified the hydrologic system of the South Coast aquifer considerably. Under predevelopment conditions, the South Coast aquifer was recharged in the unconfined, proximal fan and some midfan areas by infrequent rainfall and seepage from streams near the fan apex. Discharge occurred as seabed seepage, baseflow discharge along the lower coastal reach of streams, seepage to coastal wetlands, or evapotranspiration in areas underlain by a shallow water table. Under development conditions, seepage from irrigation canals and areal recharge from furrow irrigation represented a principal mechanism for recharge to the aquifer. Increased ground-water withdrawals in the 1960's and 1970's resulted in declines in the water table to below sea level in some places and intrusion of salt water into the aquifer. By the middle 1980's, a reduction in ground-water withdrawals and a shift from furrow irrigation to drip-irrigation techniques resulted in the recovery of water levels. Under present-day (1986) conditions, regional ground-water flow is coastward but with local movement to some well fields. In addition to the discharge mechanisms described above, ground-water discharges also to coastal canals. The North Coast limestone aquifer system consists of limestone, lesser amounts of dolomite, and minor clastic detritus of Oligocene to Pliocene age that form an unconfined upper aquifer and a confined lower aquifer; these aquifers are separated by a clay, mudstone, and marl confining unit. Topographic relief and incision of carbonate coastal plain rocks by streams are the principal factors controlling the direction of ground-water flow. The North Coast limestone aquifer system is recharged principally by precipitation that enters the upper and lower aquifers where they crop out. Regional groundwater movement from the upper aquifer is to the major rivers, wells, coastal wetlands, coastal, nearshore, and offshore springs, or as seabed seepage. Regional discharge from the lower aquifer is to the major rivers along its unconfined parts or where the confining unit has been breached by streams. Discharge from the lower aquifer also occurs in the San Juan area where the Mucarabones Sand provides an avenue for diffuse upward ground-water flow. Transmissivity within the upper limestone aquifer appears to be largely regulated by the thickness of the freshwater lens. The lens is thickest and transmissivity is greatest in interstream areas that lie in a zone that closely corresponds to the landwardmost extent of the underlying saltwater wedge. Hydraulic conductivity of the upper aquifer generally increases in a coastward direction and reflects lithologic control, karstification in the upper 30 to 100 meters of the section, and enhanced permeability in a zone of freshwater and saltwater mixing. Transmissivity of the lower aquifer is an order of magnitude smaller than that of the upper aquifer; highest transmissivities in the lower aquifer largely correspond to a coarse grainstone-packstone and coral-patch-reef depositional facies contained within the outcropping parts of the Montebello Limestone Member and its subsurface equivalents. Porosity within the North Coast limestone aquifer system is high in grainstone-packstones and low in wackestone and marl. Dolomitized zones and moldic grainstone-packstone strata are the most porous carbonate rocks, but occur in thin beds that usually are only a few meters thick. Processes of karstification that include the development of caverous zones and large vugs, and dissolution along possible regional fracture sets has enhanced permeability within the upper part of the aquifer system. Stratigraphic and lithologic control play an important role controlling permeability within the lower part of the system. The Kingshill aquifer of St. Croix, in large part, is composed of deepwater limestone that contains only microscopic pores and is poorly permeable; however, the upper part of the aquifer, a shallow-water skeletal and reef limestone, is fairly permeable, but restricted in areal extent. Permeability within these uppermost beds of the aquifer has been enhanced by meteoric leaching, dissolution within a mixing zone of saltwater and fresh water, and dolomitization. However, most large-yield wells completed in the Kingshill aquifer are also screened in alluvium that overlies or infills incised channels. The alluvial deposits serve as a temporary storage zone for rainfall, runoff, and ground water slowly entering the Kingshill aquifer.

Puerto Rico↗

Summary appraisals of the nation's ground-water resources – Lower Colorado region

This report summarizes ground-water availability in the lower Colorado region and discusses the potential for greater ground-water development and increased efficiency of water use. The climate in the most highly developed southwestern part of the region is warm and dry and that in the northeastern part is cool and moist to dry. Although the regional average annual precipitation is only about 14 inches and most streambeds are dry during most of the year, about 1.5 billion acre-feet of ground water of moderate to good chemical quality is stored in aquifers of the region. Much of the water use is founded on pumped withdrawal of ground water. However, in most of the southwestern part of the region pumpage and consumptive use are in excess of replenishment, resulting in declining water levels. In the southwestern part of the region, water levels generally are from 200-500 feet below land surface and in large areas are less than 200 feet. Large-diameter water wells in this part of the region commonly produce 500-1 ,500 gallons per minute of water. In the northeastern part of the region, water levels generally are more than 500 feet below land surface, but in some large areas water levels are from 200-500 feet, and locally are less than 200 feet below the surface. In the northeastern part of the region, water wells generally have lower yields than those to the southwest. The yields range from a few to 2,000 gallons per minute, but most wells yield from 10 to less than 500 gallons per minute. In the southwestern part of the region about 1 billion acre-feet of ground water is recoverable from storage from the water table to a depth of 700 feet below the land surface. In the northeast about 150 million acre-feet can be recovered by dewatering a 100-foot-thick section of a typical aquifer (16 million-acre area). The estimated current annual rate of ground-water depletion, occurring almost entirely in the southwestern part of the region, is 2.4-3.2 million acre-feet per year. Almost 6 million acre-feet is pumped each year-about 90 percent for agriculture, 6 percent for public supply, and 3 percent for industrial use. The percentage of use for public supply and industry is increasing each year. The potential for greater development of ground water in the southwestern part of the region is constrained by land subsidence, earth cracks, increasing costs of pumping and transportation, and moderate to poor chemical quality of water. More ground water can be developed in the northeastern part of the region, where the major constraint is pumping cost owing to low to moderate well yields and depth to water. Some benefits can be realized everywhere in the region through changes in current use and greater efficiencies of use. Additional supplies may be made available by capture of natural evapotranspiration. Increasing the efficiency of use is possible hydrologically but, in the near term, is more expensive than increasing groundwater development. Decrease of irrigation, change to water-saving methods of irrigation, use of saline water, decrease of per capita public- supply use, and more reuse of water in almost every type of use could help extend the supply and thereby reduce the current rate of ground-water depletion. Financial problems have not yet caused an overall decrease in pumpage, but, locally, operating costs or partial dewatering of the aquifer has eliminated or decreased withdrawal. Current water laws in all States of the region, except Arizona, control or allocate the use of ground water.

Arizona, California, Nevada, New Mexico, Utah↗

City-scale geothermal energy everywhere to support renewable resilience – A transcontinental cooperation

Cities have important and varying incentives to transform their energy sector to all-electric with low carbon emissions. However, they often encounter a number of impediments when attempting to implement such a change. For example, while urban areas have the highest energy demand-density, cities often lack the space for installing additional energy generation and/or long-duration energy storage systems. Cities also have existing environmental issues from energy sources (e.g., pollution from dust, waste heat or noise) that make residents sensitive to energy infrastructure development. Utilizing power from conventional sources, such as natural gas, biomass and hydropower, which usually are distanced from the urban areas, also make cities more vulnerable to supply disruptions. One promising de-carbonizing energy option for cities focuses on their heating and cooling needs, which constitutes around 1/3 of U.S. and 1/2 of European energy consumption (including industrial processes like drying, pasteurization, etc.; Jadun and others, 2017; EU Commission 2022). If heating and cooling loads can be met by geothermal direct-use technologies, then the need for new electric sources can be greatly lessened. Despite the proven efficacy of geothermal energy as a city/community-scale heating and cooling resource, it is currently only a niche resource in the heating and cooling sector, though has significant potential for future growth. Historically, emphasis has been placed on geothermal electricity generation potential that requires higher temperature (greater than 90 °C) resources at drillable depths, but potentially viable areas are geographically limited and typically well removed from urban centers. Key drivers for investments were represented by greater political interest in renewable electricity production, higher revenues and less effort in distributing the produced energy via grids. In contrast, low-temperature (less than 90 °C) geothermal resources can be used directly for heating and cooling almost everywhere and are cost-effective in urban/suburban settings. In addition, the increased prominence of renewable electricity sources, such as wind and solar onto city-scale electric grids, has led to new urgency around questions of energy storage. Underground thermal energy storage (UTES), wherein surplus or waste heat is stored underground for later use, could present a long-duration energy storage solution. From October 2022 through September 2024, a transcontinental consortium consisting of geological surveys, geoscience organizations, industry representatives and universities aims to develop an understanding of the global potential for city-scale geothermal, proposing guidelines to aid in promoting the economic utilization of low temperature geothermal resources. Efforts will focus on providing city managers and other decision makers with the information needed to evaluate and implement suitable city/community-scale geothermal technologies. Funded by the U.S. Geological Survey’s John Wesley Powell Center for Analysis and Synthesis, this interdisciplinary consortium will showcase tools, datasets, and scientific recommendations to accelerate the broader understanding and adoption of renewable energy systems that access geothermal resources. The collaborative research activities include standardization of nomenclatures, resource description and characterization strategies globally. The results from these activities will be combined with a preliminary climate-driven, city-based energy needs related analysis to perform energy supply/demand matching analysis. The identification of city-specific applications that would benefit from the geothermal technologies provides the basis to up-scale city-specific determinations to regional and national assessments of resource estimates. The city-scale geothermal energy research initiative will ultimately provide the synergies and management analysis that can address benefits, environmental impacts, regulatory frameworks, sustainability, and suitability in retrofitted buildings or new as well as existing heating networks.

Conference Paper↗

Environmental setting, water budget, and stream assessment for the Broad Run watershed, Chester County, Pennsylvania

The Broad Run watershed lies almost entirely in West Bradford Township, Chester County, Pa., and drains 7.08 square miles to the West Branch Brandywine Creek. Because of the potential effect of encroaching development and other stresses on the Broad Run watershed, West Bradford Township, the Chester County Water Resources Authority, and the Chester County Health Department entered into a cooperative study with the U.S. Geological Survey to complete an annual water budget and stream assessment of overall conditions. The annual water budget quantified the basic parameters of the hydrologic cycle for the climatic conditions present from April 1, 2003, to March 31, 2004. These water-budget data identified immediate needs and (or) deficits that were present within the hydrologic cycle during that period, if present; however, an annual water budget encompassing a single year does not identify long-term trends. The stream assessment was conducted in two parts and assessed the overall condition of the watershed, an overall assessment of the fluvial-geomorphic conditions within the watershed and an overall assessment of the stream-quality conditions. The data collected will document present (2004) conditions and identify potential vulnerabilities to future disturbances. For the annual period from April 1, 2003, to March 31, 2004, determination of an annual water budget indicated that of the 67.8 inches of precipitation that fell on the Broad Run watershed, 38.8 inches drained by way of streamflow to the West Branch Brandywine Creek. Of this 38.8 inches of streamflow, local-minimum hydrograph separation techniques determined that 7.30 inches originated from direct runoff and 31.5 inches originated from base flow. The remaining precipitation went into ground-water storage (1.71 inches) and was lost to evapotranspiration (27.3 inches). Ground-water recharge for this period-35.2 inches-was based on these values and an estimated ground-water evapotranspiration rate of 2 inches. Assessment of fluvial-geomorphic conditions included large-scale mapping of stream classes within the Broad Run watershed and in-depth study of three representative stream reaches also within the Broad Run watershed. Based on the total distance of all stream reaches classified within the Broad Run watershed, 61 percent were classified as C-class, 14 percent as E-class, 13 percent as B-class, 5 percent as F-class, 4 percent as undifferentiated B- and F-class, 2 percent as G-class, and less than 1 percent as A-class. The map of stream classes indicates that the Broad Run watershed currently has no large-scale areas of stream impairment and that, generally, the stream is not entrenched and the main branch of the Broad Run has an available, functioning flood plain. Smaller tributary streams, however, showed signs of localized entrenchment due to site-specific influences such as natural stream-channel evolution, localized channelization, localized contraction due to road and driveway crossings, and (or) increased localized runoff. For example, one small reach along a tributary channel was observed to become entrenched due to runoff originating from a new housing development. Entrenched stream reaches are merely located by large-scale mapping and require individual assessment to determine potential causes of entrenchment and (or) future restorative actions. Three in-depth geomorphic study sites showed that the Broad Run can currently be considered graded or in a state of dynamic equilibrium. The sites did, however, identify certain vulnerabilities to future changes within the watershed. These vulnerabilities included disruption of the present sediment supply, including both increases and (or) reductions in the current sediment loads within the Broad Run; increases in both magnitude and duration of storm-water runoff; encroachment of development onto present flood-plain areas, and (or) alterations to riparian zones. Assessment of stream-quality conditions includ

Scientific Investigations Report↗

Hydrology of the Yucaipa groundwater subbasin: Characterization and integrated numerical model, San Bernardino and Riverside Counties, California

Executive Summary Water management in the Santa Ana River watershed in San Bernardino and Riverside Counties in southern California is a complex task with various water purveyors navigating geographic, geologic, hydrologic, and political challenges to provide a reliable water supply to stakeholders. As the population has increased throughout southern California, so has the demand for water. The Yucaipa groundwater subbasin (hereafter referred to as “Yucaipa subbasin”), one of nine groundwater subbasins in what the California Department of Water Resources (DWR) refers to as the Upper Santa Ana Valley groundwater basin (California Department of Water Resources, 2016; the DWR naming convention is used within this report), is no exception; steady population growth since the 1940s and changes in water use has forced local water purveyors to regularly adapt their water infrastructure to meet demand. Groundwater has historically been the dominant source of water in the Yucaipa subbasin although recently, imported water via the California State Water Project has augmented the total water supply. Despite the influx of imported water, overall demand for groundwater continues to rise, and there is concern by local water managers that groundwater levels may adversely impact water supply and (or) decline to a point where it will be uneconomical to produce water, severely limiting the ability of local agencies to meet water-supply demand. To better understand the hydrogeology and water resources in the Yucaipa subbasin, the U.S. Geological Survey (USGS) and the San Bernardino Valley Municipal Water District initiated a cooperative study to understand the hydrogeologic system of the Yucaipa subbasin and in the encompassing Yucaipa Valley watershed (YVW). A three-dimensional hydrogeologic framework model was constructed to quantify the structure and extent of hydrogeologic units. Historical and present-day groundwater conditions were characterized to evaluate the groundwater-flow system. Lastly, the Yucaipa Integrated Hydrological Model (YIHM) was developed to simulate the integrated surface-water and groundwater systems, including natural and anthropogenic (that is, human influenced) recharge and discharge throughout the study area from 1947 to 2014. The Yucaipa subbasin is an inland groundwater basin located about 12 miles (mi) southeast of the City of San Bernardino and about 75 mi east of Los Angeles, California. The subbasin encompasses about 39 square miles (mi 2 ), including the City of Yucaipa. The geographic extent of the Yucaipa subbasin was established by the California Department of Water Resources, who defined the boundaries of the subbasin based on hydrogeologic transitions between crystalline rock and basin-fill sediments, active fault strands, surface-water drainage divides, and a portion of an adjudicated groundwater management boundary. Two groundwater subbasins of the Upper Santa Ana Valley groundwater basin are adjacent to the Yucaipa subbasin, the San Bernardino groundwater subbasin to the west and the San Timoteo groundwater subbasin to the south. The Yucaipa subbasin is encompassed by the YVW, which is in turn comprised of three sub-watersheds that represent surface-water flow across and within the Yucaipa subbasin. Although the Yucaipa subbasin is the specific area of interest for this study, the entire YVW was considered for the purposes of characterizing the hydrogeology of the Yucaipa subbasin and for development of the YIHM. The purposes of this report are to (1) describe the hydrologic and hydrogeologic settings of the Yucaipa subbasin and aquifer system, (2) describe the construction and calibration of the fully coupled groundwater and surface-water flow model for the Yucaipa subbasin and the encompassing YVW, referred to as the YIHM, and (3) present numerical results, including water budgets and hydraulic heads, and the effect of pumping and climate stresses (precipitation and temperature) on water-budget components.

California↗

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

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

Palaeogeography, Palaeoclimatology, Palaeoecology↗

Global plant trait relationships extend to the climatic extremes of the tundra biome

The majority of variation in six traits critical to the growth, survival and reproduction of global plant species is thought to be organised along just two dimensions, corresponding to strategies of plant size and resource acquisition. However, it is unknown whether global plant trait relationships extend to climatic extremes, and if these interspecific relationships are confounded by trait variation within species. We test whether trait relationships extend to the cold extremes of life on Earth using the largest database of tundra plant traits yet compiled. We show that tundra plants demonstrate remarkably similar resource economic traits, but not size traits, compared to global distributions, and exhibit the same two dimensions of trait variation. Three quarters of trait variation occurs among species, mirroring global estimates of interspecific trait variation. Plant trait relationships are thus generalizable to the edge of global trait-space, informing prediction of plant community change in a warming world.

Nature Communications↗

Community estimate of global glacier mass changes from 2000 to 2023

Glaciers are indicators of ongoing anthropogenic climate change 1 . Their melting leads to increased local geohazards 2 , and impacts marine 3 and terrestrial 4 , 5 ecosystems, regional freshwater resources 6 , and both global water and energy cycles 7 , 8 . Together with the Greenland and Antarctic ice sheets, glaciers are essential drivers of present 9 , 10 and future 11 , 12 , 13 sea-level rise. Previous assessments of global glacier mass changes have been hampered by spatial and temporal limitations and the heterogeneity of existing data series 14 , 15 , 16 . Here we show in an intercomparison exercise that glaciers worldwide lost 273 ± 16 gigatonnes in mass annually from 2000 to 2023, with an increase of 36 ± 10% from the first (2000–2011) to the second (2012–2023) half of the period. Since 2000, glaciers have lost between 2% and 39% of their ice regionally and about 5% globally. Glacier mass loss is about 18% larger than the loss from the Greenland Ice Sheet and more than twice that from the Antarctic Ice Sheet 17 . Our results arise from a scientific community effort to collect, homogenize, combine and analyse glacier mass changes from in situ and remote-sensing observations. Although our estimates are in agreement with findings from previous assessments 14 , 15 , 16 at a global scale, we found some large regional deviations owing to systematic differences among observation methods. Our results provide a refined baseline for better understanding observational differences and for calibrating model ensembles 12 , 16 , 18 , which will help to narrow projection uncertainty for the twenty-first century 11 , 12 , 18 .

Nature↗

Drivers of seedling establishment success in dryland restoration efforts

Restoration of degraded drylands is urgently needed to mitigate climate change, reverse desertification and secure livelihoods for the two billion people who live in these areas. Bold global targets have been set for dryland restoration to restore millions of hectares of degraded land. These targets have been questioned as overly ambitious, but without a global evaluation of successes and failures it is impossible to gauge feasibility. Here we examine restoration seeding outcomes across 174 sites on six continents, encompassing 594,065 observations of 671 plant species. Our findings suggest reasons for optimism. Seeding had a positive impact on species presence: in almost a third of all treatments, 100% of species seeded were growing at first monitoring. However, dryland restoration is risky: 17% of projects failed, with no establishment of any seeded species, and consistent declines were found in seeded species as projects matured. Across projects, higher seeding rates and larger seed sizes resulted in a greater probability of recruitment, with further influences on species success including site aridity, taxonomic identity and species life form. Our findings suggest that investigations examining these predictive factors will yield more effective and informed restoration decision-making.

Nature Ecology and Evolution↗