USGS ScienceSearch

SEARCH · USGS Science

Results for “Regional Studies in Marine Science”

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

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

At least 127 records · Page 7Linked to original sources

Submarine groundwater discharge creates cold‐water refugia that can mitigate exposure of heat stress in nearshore corals

Coral reef mortality around the world is accelerating due to human activities and rising sea temperatures that cause bleaching, which is expected to become more frequent. Our ability to predict which corals will be most resilient, however, remains limited due to insufficient information characterizing nearshore temperature and habitat conditions. In this study, we examine how submarine groundwater discharge (SGD) reduces nearshore water temperatures and exposure of corals to heat stress, complementing the understanding that SGD can adversely affect coral when it contains elevated nutrient concentrations. Data from fixed nearshore sensors and vertical depth profiles along ~100 km of the western shoreline of the Island of Hawai’i from 2003 to 2014 demonstrate that submarine groundwater discharge (SGD) can reduce nearshore water temperatures by 1 °C–5°C and create estuarine-like conditions with salinities as low as 20 PSU, where the prevalent coral species, Pocillopora meandrina, Porites lobata , and Montipora capitata , thrive. Time-series temperature records reveal that exposure to high ambient ocean temperatures, which are known to initiate bleaching events, are reduced up to 5%–46% of the time. Coral health surveys indicated coral bleaching in response to moderately high annual temperatures in 2010 and 2011, with more colonies affected farther from cold, SGD-fed waters. Synthesis of these results, along with coral response data following the more extreme marine heat wave of 2014–2015, demonstrates lower coral loss and greater coral recovery near groundwater seeps, particularly those with higher flux and influence on reducing nearshore water temperatures. Our results demonstrate that SGD may therefore provide a beneficial ecosystem service and enhance coral reef resilience, particularly where human-related nutrient additions to groundwater can be mitigated. The implications of our findings are relevant across tropical coasts where groundwater inputs can be substantial, such as the Caribbean and Indo-Pacific, and contribute to improving our understanding of coral sensitivity to gradients in temperature and nutrient stress. Improved management of groundwater resources could thus be vital to local–regional strategies for mitigating future heat stress.

Hawaii

The statistics and kinematics of transverse sand bars on an open coast

Ten years (1987-1996) of time exposure video images of the nearshore region at Duck, NC were used to study transverse sand bars, bathymetric features of intermediate length scales (10-200 m) oriented oblique or perpendicular to the shoreline. These transverse sand bars extend seaward from both the shoreline (trough transverse bars) and the shore-parallel sand bar (offshore transverse bars). Transverse bars had not previously been observed in an energetic Coastal environment such as that at Duck, and their dynamics and role in nearshore processes is unknown. Frequency of occurrence statistics and length scales of the transverse bars were calculated using the video images. Trough and offshore transverse bars appeared a mean of 39 and 73 days per year, respectively. The offshore bars were found to be much larger features than the trough bars, with mean wavelengths (alongshore spacing between consecutive crests) of 79 and 172 m for trough and offshore bars, respectively. Both the trough and offshore bars were found to persist for periods of days to months. The alongshore movement of the bars was measured and compared to estimates of surf zone longshore currents which were calculated from wave height and wave angle data. Both sets of bars were observed to move at rates up to 40 m/day. At times, both trough and offshore bars were observed shifting in the same direction as the current was flowing, and at other times, both sets of bars remained stationary, even under relatively strong longshore currents. Trough bars were also observed moving against the current. An hypothesis, proposed by Barcilon and Lau (1973) [J. Geophys. Res. 78(15): 2656-2664], that the transverse bars were created as a sea bed instability under longshore currents, was tested by comparing the magnitude of estimated surf zone longshore currents with times of formation or presence of transverse bars. There was no evidence to suggest that the bars were formed by this simple longshore current instability mechanism. Instead, it is plausible that the combined effects of waves and currents may drive the formation of these features. (C) 2000 Elsevier Science B.V.

Marine Geology

Juxtaposition of Neoproterozoic units along the Baruda - Tulu Dimtu shear-belt in the East African Orogen of western Ethiopia

Amalgamation of East and West Gondwanaland during the Neoproterozoic East African Orogen is recorded by several shear-belts or 'suture zones', some of which are associated with ultramafic and mafic complexes that have been interpreted as ophiolite fragments. The Baruda shear-belt is a major structure of this type that belongs to the N-S trending Barka - Tulu Dimtu zone. The significance of this zone has been studied within a transect in western Ethiopia which covers a variety of metasedimentary and metavolcanic sequences, ultramafic rocks and synkinematic intrusive complexes. All rocks participated in the regional D1 event as reflected in a penetrative steep foliation in supracrustal rocks and marginal parts of the intrusions. Highly strained rocks contain a stretching lineation that plunge to the east. The several-km thick Baruda shear-belt, comprising mylonitic supracrustal and plutonic rocks including mafic-ultramafic mega-lenses, is the most prominent expression of this event. Shear-sense indicators demonstrate top-to-the-west shear. Subsequent D2 deformation is recorded in 2-300 m wide, N-S striking, subvertical shear-zones with subhorizontal stretching lineation relatable to sinistral transcurrent movements. Our data indicate that rock units on either side of the Baruda shear-belt are related, rather than being exotic to each other as implied in suture zone models, since there is no major lithologic or metamorphic difference, geochemical data on metavolcanic rocks and pre-tectonic intrusions suggest a paleotectonic link, and style and extent of deformation is similar across the shear-belt. A tentative model for the transect suggests an arc and back-arc setting which experienced later continental collision and tectonic shortening. The initial setting was that of a shallow marine platform characterised by carbonates and sandstones, which covered extensive areas prior to break-up of a pre-existing supercontinent. Continental convergence is first recorded in high-K calc-alkaline volcanism characterised by pyroclastic deposits of andesitic composition, at an active continental margin at about 800 Ma. Subaerial arc volcanism was temporally and spatially overlapping with limited arc rifting, represented by submarine basalts compositionally transitional between enriched MORB and calc-alkaline magmas, and associated dyke swarms in the older carbonate-sandstone platform sequence. It is suggested that the large, mafic-ultramafic, bodies relate to this event and were originally formed as intrusions along one or more propagating rift axis within the arc complex. The regional Baruda shear-belt formed in response to contractional D1 deformation, and its location may have been largely controlled by competence contrasts between the array of rift-related intrusions and the marble-dominated lithologies. Associated shortening of the arc and back-arc region led to crustal thickening and emplacement of synkinematic, composite, batholiths at about 570-550 Ma. These are composed of moderately peraluminous granite and coeval, intermediate to mafic intrusions of shoshonitic affinity. D2 sinistral movements succeeded the contractional deformation. ?? 2001 Elsevier Science B.V.

Precambrian Research

Sources of terrestrially-derived organic carbon in lower Mississippi River and Louisiana shelf sediments: Implications for differential sedimentation and transport at the coastal margin

In this study, we examined the temporal and spatial variability of terrestrial organic carbon sources in lower Mississippi River and Louisiana shelf sediments (during 11 cruises over a 22-month period) to further understand the sorting dynamics and selective transport of vascular plant materials within the primary dispersal system of the river. Bulk ??13C values in lower river sediments ranged from -21.90??? to -24.64??? (mean=-23.20??1.09???), these values were generally more depleted than those found in shelf sediments (-22.5??? to -21.2???). The ??8 (??8 = sum of vanillyl, syringyl and cinnamyl phenols produced from the oxidation of 100 mg of organic carbon) values in the lower river ranged from 0.71 to 3.74 (mean = 1.78??0.23). While there was no significant relationship between ??8 and river discharge (p>0.05), the highest value occurred during peak discharge in April 1999-which corresponded to the highest observed C/N value of 17.41. The ??8 values on the shelf ranged from 0.68 to 1.36 (mean = 0.54??0.30) and were significantly lower (p <0.05) than the average value for lower river sediments. The range of S/V (syringyl/vanillyl) and C/V (cinnamyl/vanillyl) ratios on the shelf, 0.11 to 0.95 and 0.01 to 0.08, respectively, were similar to that found in the lower river. These low C/V ratios are indicative a mixture of woody and non-woody carbon sources. Recent work by Goni et al. [Nature 389 (1997) 275; Geochim. Cosmochim. Acta 62 (1998) 3055], which did not include sampling transects within the primary dispersal system of the Mississippi River, showed a non-woody vascular plant signature on the Louisiana shelf. This suggests that riverine-derived woody tissues preferentially settle out of the water column, in the lower river and inner shelf, prior to the selective dispersal of C3 versus C4 non-woody materials in other regions the shelf and slope. This works further demonstrates the importance of differential settlement of particles, sampling location within the dispersal system, and river discharge, when examining biogeochemical cycles in river-dominated margins. ?? 2002 Elsevier Science B.V. All rights reserved.

Marine Chemistry

Tectonic map of the Circum-Pacific region, Pacific basin sheet

Circum-Pacific Map Project : The Circum-Pacific Map Project was a cooperative international effort designed to show the relationship of known energy and mineral resources to the major geologic features of the Pacific basin and surrounding continental areas. Available geologic, mineral, and energy-resource data are being complemented by new, project-developed data sets such as magnetic lineations, seafloor mineral deposits, and seafloor sediment. Earth scientists representing some 180 organizations from more than 40 Pacific-region countries are involved in this work. Six overlapping equal-area regional maps at a scale of 1:10,000,000 form the cartographic base for the project: the four Circum-Pacific Quadrants (Northwest, Southwest, Southeast, and Northeast), and the Antarctic and Arctic Sheets. There is also a Pacific Basin Sheet at a scale of 1:17,000,000. The Base Map Series and the Geographic Series (published from 1977 to 1990), the Plate-Tectonic Series (published in 1981 and 1982), the Geodynamic Series (published in 1984 and 1985), and the Geologic Series (published from 1984 to 1989) all include six map sheets. Other thematic map series in preparation include Mineral-Resources, Energy-Resources and Tectonic Maps. Altogether, more than 50 map sheets are planned. The maps were prepared cooperatively by the Circum-Pacific Council for Energy and Mineral Resources and the U.S. Geological Survey and are available from the Branch of Distribution, U. S. Geological Survey, Box 25286, Federal Center, Denver, Colorado 80225, U.S.A. The Circum-Pacific Map Project is organized under six panels of geoscientists representing national earth-science organizations, universities, and natural-resource companies. The six panels correspond to the basic map areas. Current panel chairmen are Tomoyuki Moritani (Northwest Quadrant), R. Wally Johnson (Southwest Quadrant), Ian W.D. Dalziel (Antarctic Region), vacant. (Southeast Quadrant), Kenneth J. Drummond (Northeast Quadrant), and George W. Moore (Arctic Region). Project coordination and final cartography was being carried out through the cooperation of the Office of the Chief Geologist of the U.S. Geological Survey, under the direction of General Chairman, George Gryc of Menlo Park, California. Project headquarters were located at 345 Middlefield Road, MS 952, Menlo Park, California 94025, U.S.A. The framework for the Circum-Pacific Map Project was developed in 1973 by a specially convened group of 12 North American geoscientists meeting in California. The project was officially launched at the First Circum-Pacific Conference on Energy and Mineral Resources, which met in Honolulu, Hawaii, in August 1974. Sponsors of the conference were the AAPG, Pacific Science Association (PSA), and the Coordinating Committee for Offshore Prospecting for Mineral Resources in Offshore Asian Areas (CCOP). The Circum-Pacific Map Project operates as an activity of the Circum-Pacific Council for Energy and Mineral Resources, a nonprofit organization that promotes cooperation among Circum-Pacific countries in the study of energy and mineral resources of the Pacific basin. Founded by Michel T. Halbouty in 1972, the Council also sponsors conferences, topical symposia, workshops and the Earth Science Series books. Tectonic Map Series : The tectonic maps distinguish areas of oceanic and continental crust. Symbols in red mark active plate boundaries, and colored patterns show tectonic units (volcanic or magmatic arcs, arc-trench gaps, and interarc basins) associated with active plate margins. Well-documented inactive plate boundaries are shown by symbols in black. The tectonic development of oceanic crust is shown by episodes of seafloor spreading. These correlate with the rift and drift sequences at passive continental margins and episodes of tectonic activity at active plate margins. The recognized episodes of seafloor spreading seem to reflect major changes in plate kinematics. Oceanic plateaus and other prominences of greater than normal oceanic crustal thickness such as hotspot traces are also shown. Colored areas on the continents show the ages of deformation and metamorphism of basement rocks and the emplacement of igneous rocks. Transitional tectonic (molassic) and reactivation basins are shown by a colored boundary, and if they are deformed, a colored horizontal line pattern indicates the age of deformation. Colored bands along basin boundaries indicate age of inception, and isopachs indicate thickness of platform strata on continental crust and cover on oceanic crust. Colored patterns at separated continental margins show the age of inception of rift and drift (breakup) sequences. Symbols mark folds and faults, and special symbols show volcanoes and other structural features. Affiliations are as of compilation of the data. This map was created in quadrants and then compiled together. They are the Northwest land, Northwest Marine (different compilers), Northeast, Southwest and Southeast, and parts in plate-boundary sections.

Circum-Pacific Map

Stony coral tissue loss disease accelerated shifts in coral composition and declines in reef accretion potential in the Florida Keys

Outbreaks of coral disease have been a dominant force shaping western Atlantic coral-reef assemblages since the late 1970s. Stony coral tissue loss disease (SCTLD) is nonetheless having an unprecedented impact in the region. Whereas numerous studies over the last decade have worked to characterize this novel pathogen and its impacts on coral populations, few have quantified its functional effects on reef ecosystems. Of particular importance is how SCTLD may be impacting the essential reef-accretion process and the myriad ecosystem services it supports. Here, we evaluated the impact of SCTLD on reef-accretion potential by estimating carbonate budgets and taxon-level carbonate production at 43 sites throughout the Florida Keys from 2016−2022. Average regional reef-accretion potential declined from an already low, but positive rate of 0.30 ± 0.16 mm y -1 (mean ± standard error) in 2016 before the disease was first observed, to a state of accretionary stasis (0.08 ± 0.12 mm y -1 ) by 2022. This 70% relative decline in reef-accretion potential was driven by the loss of reef-building corals, with significant decreases in carbonate production by massive taxa including Colpophyllia natans , Montastraea cavernosa , Pseudodiploria strigosa , Orbicella spp., and Siderastrea siderea , and increasing contributions from less susceptible, weedy taxa including Millepora spp., Agaricia spp., and Porites astreoides . In general, changes in taxon-level carbonate production following the SCTLD outbreak mirror long-term shifts in reef assemblages in response to previous stressors. One striking exception, however, is S. siderea , which had become increasingly dominant in recent decades, but declined significantly in response to SCTLD. Overall, by further decimating the already depauperate reef-building coral populations in the Florida Keys, SCTLD has caused a functionally significant shift in the composition of Florida’s coral-reef assemblages and accelerated the loss of regional reef-building capacity. The dire impacts of the disease in south Florida may serve as an early warning that the persistence of the invaluable ecological and socioeconomic functions coral reefs provide will be increasingly threatened throughout the western Atlantic in the aftermath of SCTLD.

Florida

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

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

Marine and Coastal Fisheries: Dynamics, Management

Chapter 24 Lateral variability of the estuarine turbidity maximum in a tidal strait

The behavior of the estuarine turbidity maximum (ETM) in response to freshwater flow, tidal forcing, and bed dynamics has been studied extensively by many researchers. However, the majority of investigations focus on the longitudinal position and strength of the ETM, which can vary over tidal, spring-neap, and seasonal timescales. ETMs may become longitudinally fixed due to bathymetric constraints, and thus the lateral position may vary significantly on differing timescales. Lateral dynamics of the ETM may affect contaminant uptake in biologically active regions, while local deposition patterns may be affected by the dominant lateral position. A longitudinally fixed ETM in Carquinez Strait, California, was studied to specifically investigate the dynamics of lateral ETM variability during April 2004. an abrupt topographical control on the north side restricts gravitational circulation resulting in convergence and particle trapping, creating the ETM. The cross-section was continuously monitored with two upward-looking velocity profilers and four optical backscatterance sensors. In addition, cross-sectional measurements over one tidal cycle were performed during a spring tide with boat-mounted velocity and water quality profilers. The lateral and vertical positions of the ETM center of mass varied by a maximum of 250 and 5 m, respectively (20% of width and 17% of depth) over the tidal timescale, while tidally averaged lateral and vertical positions varied substantially less (50 and 1 m, respectively). ETM position responded to tidal energy (Urms), with higher vertical position and a laterally centered position resulting from increased mixing during spring tides, and a northerly lateral position from decreased mixing during neap tides. Hydrodynamic and sediment transport modeling of this period reproduces the lateral and vertical movements of the ETM center of mass. Modeling results indicate increased gravitational circulation in the strait and enhanced particle trapping on the north side during neap tides, thus displacing the ETM center of mass to the north. The south side has no topographical control, and therefore no particle trapping mechanism exists on the south side. Secondary circulation is strengthened on spring tides, distributing near-bed sediment toward the south. The field and modeling results are in agreement with previous work in Carquinez Strait and further elucidate the strong lateral variation of the ETM, even in narrow, energetic tidal straits. ?? 2008 Elsevier B.V. All rights reserved.

Proceedings in Marine Science

An exploration of open-system uranium-series history of marine terrace corals, Perachora Peninsula, Greece and San Nicolas Island, California, USA

Emergent marine terraces record past periods of relatively high sea level and are common on uplifting coasts worldwide. Such landforms are extensive around the Gulf of Corinth area of Greece and along the coast of California, USA. In the Gulf of Corinth region, marine terraces record Quaternary uplift due to ongoing basin extension on the northern side of the Africa-Eurasia-Aegean Sea plate boundary. Although dating of terrace fossils has been extensive in the Perachora Peninsula part of this region, some apparent ages are much older (∼190 ka to ∼150 ka) and some are much younger (∼91 ka to ∼69 ka) than expected for a marine terrace that has long been hypothesized to be of last-interglacial (∼120 ka) age. New uranium-series analyses of terrace corals from two previously studied localities are reported here. As with previous work, the new results show that these corals have behaved as open systems, with back-calculated, initial 234 U/ 238 U values higher than those found in modern seawater. New and previously published U-series analyses of last-interglacial corals from San Nicolas Island, California (USA) illustrate the possible mechanisms of open-system behavior, which include both “early” alpha-recoil-derived additions of 230 Th and 234 U, and “continuing” alpha-recoil-derived additions of these nuclides. Many corals from the Perachora Peninsula appear to have experienced the “continuing” alpha-recoil-derived open-system history, which explains the apparent older ages, and permits an interpretation that all analyzed corals likely date to the last interglacial period, ∼120 ka. If the estimated age of ∼120 ka reported here is correct, this age, along with considerations of glacial isostatic adjustment (GIA) modeling of past sea level, indicates that the southern side of the Perachora Peninsula has experienced uplift rates ranging from ∼0.23 m/kyr to ∼0.57 m/kyr, increasing to the southeast, in agreement with previous work. Interpretation of the geomorphology of the northern side of the peninsula has been controversial for more than four decades. The “stepped” topography has been interpreted as either a series of successively higher and older marine terraces, as on the southern side, or a marine terrace of one age that has been displaced by a series of shore-parallel, northward-dipping normal faults. In the second interpretation, topographic steps in the landscape are interpreted to be fault scarps, not ancient sea cliffs. Interpretation of the open-system U-series ages indicates that marine deposits at elevations of ∼22 m and ∼44 m are both likely ∼120 ka, supporting the scenario that steps in the landscape are likely fault scarps, not sea cliffs. Using the new U-series ages of ∼120 ka, along with GIA-modeled paleo-sea level confirms another previous study that uplift rates on the northern side of the peninsula are at least as high as 0.39–0.46 m/kyr.

California

Geologic implications of gas hydrates in the offshore of India: results of the National Gas Hydrate Program Expedition 01

The Indian National Gas Hydrate Program Expedition 01 (NGHP-01) is designed to study the occurrence of gas hydrate along the passive continental margin of the Indian Peninsula and in the Andaman convergent margin, with special emphasis on understanding the geologic and geochemical controls on the occurrence of gas hydrate in these two diverse settings. The NGHP-01 expedition established the presence of gas hydrates in the Krishna-Godavari and Mahanadi Basins, and the Andaman Sea. The expedition discovered in the Krishna-Godavari Basin one of the thickest gas hydrate accumulations ever documented, in the Andaman Sea one of the thickest and deepest gas hydrate stability zones in the world, and established the existence of a fully developed gas hydrate petroleum system in all three basins. The primary goal of NGHP-01 was to conduct scientific ocean drilling/coring, logging, and analytical activities to assess the geologic occurrence, regional context, and characteristics of gas hydrate deposits along the continental margins of India. This was done in order to meet the long-term goal of exploiting gas hydrate as a potential energy resource in a cost effective and safe manner. During its 113.5-day voyage, the D/V JOIDES Resolution cored and/or drilled 39 holes at 21 sites (1 site in Kerala-Konkan, 15 sites in Krishna-Godavari, 4 sites in Mahanadi, and 1 site in the Andaman deep offshore area), penetrated more than 9250 m of sedimentary section, and recovered nearly 2850 m of core. Twelve holes were logged with logging-while-drilling (LWD) tools and an additional 13 holes were wireline logged. The science team utilized extensive on-board laboratory facilities to examine and prepare preliminary reports on the physical properties, geochemistry, and sedimentology of all the data collected prior to the end of the expedition. Samples were also analyzed in additional post-expedition shore-based studies conducted in leading laboratories around the world. One of the specific objectives of this expedition was to test gas hydrate formation models and constrain model parameters, especially those that account for the formation of concentrated gas hydrate accumulations. The necessary data for characterizing the occurrence of in situ gas hydrate, such as interstitial water chlorinities, core-derived gas chemistry, physical and sedimentological properties, thermal images of the recovered cores, and downhole measured logging data (LWD and/or conventional wireline log data), were obtained from most of the drill sites established during NGHP-01. Almost all of the drill sites yielded evidence for the occurrence of gas hydrate; however, the inferred in situ concentration of gas hydrate varied substantially from site to site. For the most part, the interpretation of downhole logging data, core thermal images, interstitial water analyses, and pressure core images from the sites drilled during NGHP-01 indicate that the occurrence of concentrated gas hydrate is mostly associated with the presence of fractures in the sediments, and in some limited cases, by coarser grained (mostly sand-rich) sediments.

Marine and Petroleum Geology

Predicting sea-level rise vulnerability of terrestrial habitat and wildlife of the Northwestern Hawaiian Islands

If current climate change trends continue, rising sea levels may inundate low-lying islands across the globe, placing island biodiversity at risk. Recent models predict a rise of approximately one meter (1 m) in global sea level by 2100, with larger increases possible in areas of the Pacific Ocean. Pacific Islands are unique ecosystems home to many endangered endemic plant and animal species. The Northwestern Hawaiian Islands (NWHI), which extend 1,930 kilometers (km) beyond the main Hawaiian Islands, are a World Heritage Site and part of the Papahanaumokuakea Marine National Monument. These NWHI support the largest tropical seabird rookery in the world, providing breeding habitat for 21 species of seabirds, 4 endemic land bird species and essential foraging, breeding, or haul-out habitat for other resident and migratory wildlife. In recent years, concern has grown about the increasing vulnerability of the NWHI and their wildlife populations to changing climatic patterns, particularly the uncertainty associated with potential impacts from global sea-level rise (SLR) and storms. In response to the need by managers to adapt future resource protection strategies to climate change variability and dynamic island ecosystems, we have synthesized and down scaled analyses for this important region. This report describes a 2-year study of a remote northwestern Pacific atoll ecosystem and identifies wildlife and habitat vulnerable to rising sea levels and changing climate conditions. A lack of high-resolution topographic data for low-lying islands of the NWHI had previously precluded an extensive quantitative model of the potential impacts of SLR on wildlife habitat. The first chapter (chapter 1) describes the vegetation and topography of 20 islands of Papahanaumokuakea Marine National Monument, the distribution and status of wildlife populations, and the predicted impacts for a range of SLR scenarios. Furthermore, this chapter explores the potential effects of SLR on wildlife breeding habitats for each island. The subsequent chapter (chapter 2) details a study of the Laysan Island ecosystem, describing a quantitative model that incorporates SLR, storm wave, and rising groundwater inundation. Wildlife, storm, and oceanographic data allowed for an assessment of the phenological and spatial vulnerability of Laysan Island's breeding bird species to SLR and storms. Using remote sensing and geospatial techniques, we estimated topography, classified vegetation, modeled SLR, and evaluated a range of climate change scenarios. On the basis of high-resolution airborne data collected during 2010-11 (root-mean-squared error = 0.05-0.18 m), we estimated the maximum elevation of 20 individual islands extending from Kure Atoll to French Frigate Shoals (range: 1.8-39.7 m) and computed the mean elevation (1.7 m, standard deviation 1.1 m) across all low-lying islands. We also analyzed general climate models to describe rainfall and temperature scenarios expected to influence adaptation of some plants and animals for this region. Outcomes for the NWHI predicted an increase in temperature of 1.8-2.6 degrees Celsius (&deg;C) and an annual decrease in precipitation of 24.7-76.3 millimeters (mm) across the NWHI by 2100. Our models of passive SLR (excluding wave-driven effects, erosion, and accretion) showed that approximately 4 percent of the total land area in the NWHI will be lost with scenarios of +1.0 m of SLR and 26 percent will be lost with +2.0 m of SLR. Some atolls are especially vulnerable to SLR. For example, at Pearl and Hermes Atoll our analysis indicated substantial habitat losses with 43 percent of the land area inundated at +1.0 m SLR and 92 percent inundated at +2.0 m SLR. Across the NWHI, seven islands will be completely submerged with +2.0 m SLR. The limited global ranges of some tropical nesting birds make them particularly vulnerable to climate change impacts in the NWHI. Climate change scenarios and potential SLR impacts presented here emphasize the need for early climate change adaptation and mitigation planning, especially for species with limited breeding distributions and/or ranges restricted primarily to the low-lying NWHI: Cyperus pennatiformis var. bryanii , Black-footed Albatross ( Phoebastria nigripes ), Laysan Albatross ( P. immutabilis ), Bonin Petrel ( Pterodroma hypoleuca ), Gray-backed Tern ( Onychoprion lunatus ), Laysan Teal ( Anas laysanensis ), Laysan Finch ( Telespiza cantans ), and Hawaiian monk seal ( Monachus schauinslandi ). Furthermore, SLR scenarios that include the effects of wave dynamics and groundwater rise may indicate amplified vulnerability to climate change driven habitat loss on low-lying islands. In chapter 2, we incorporated the combined effects of SLR, dynamic wave-driven inundation, and rising groundwater in a quantitative study specifically for the Laysan Island ecosystem. This is the first hydrodynamic model to simulate the combined impacts of SLR and wave-driven inundation in the NWHI. We developed a high-resolution digital elevation model (mean vertical accuracy of 0.32 m) for the island. Then using recent satellite imagery, geospatial models, and historical oceanographic, storm, and biological data we estimated potential inundation extent, habitat loss, and wildlife population impacts for a range of potential SLR scenarios (0.00, +0.50, +1.00, +1.50, and +2.00 m) that may occur over the next century. Additionally, we estimated the carrying capacity of Laysan Island for five species based on the available population monitoring data and described how potential losses in nesting habitat could influence population dynamics for Black-footed Albatross, Laysan Albatross, Red-footed Booby (Sula sula), Laysan Teal, and Laysan Finch. For some other seabird populations (Masked Booby, S. dactylatra ; Brown Booby, S. leucogaster ; Great Frigatebird, Fregata minor ; and Sooty Tern, Onychoprion fuscata ), we used recent colony distribution data, land cover maps, and nesting behavior to estimate potential losses of nesting habitat from SLR and wave-driven inundation. We observed far greater potential impacts of SLR to wildlife with the dynamic wave-driven modeling approach than with the passive modeling approach. Depending on SLR scenario and coastal orientation, during storms under a +2.00 m SLR scenario, the wave-driven inundation model predicted three times more inundation than the passive model (17.2 percent of total terrestrial area versus 4.6 percent, respectively). Large-wave events generally added 1 m of water height to passive inundation surfaces, therefore our dynamic models (during storm events) forecasted comparable inundation extents earlier than passive models. Although wave-driven water levels were highest in the northwest quadrant of Laysan Island, the greatest extent of inundation occurred in the southeast where coastal dunes less than 3 m above mean sea level provide little protection from wave-driven inundation. When wave-driven inundation was included in the SLR model for Laysan Island greater nesting habitat loss and potential impacts on wildlife population dynamics were predicted. The consequences of habitat loss due to SLR may be worse for species with colonies in the wave-exposed coastal zones (for example, Black-footed Albatross) and for populations already near the island's carrying capacity (for example, Laysan Teal). Species whose peak incubation and chick-rearing periods coincide with seasonally high wave heights also will be increasingly vulnerable, especially those species nesting on the ground in areas vulnerable to inundation, such as Gray-backed Tern and Black-footed Albatross. Other species that have space for population growth, or are not restricted to a narrow range of habitat types on Laysan (for instance, Sooty Terns), may be less sensitive to habitat loss from SLR over the next century. Our assessments of inundation risk, habitat loss, and wildlife species vulnerability synthesize current knowledge about individual islands and contribute to a broader understanding of the impacts of inundation from SLR and storm-induced waves. Yet, most NWHI and their bird populations lack monitoring data to evaluate adaptations to and impacts of climate change. Exceptions include some data sets from long-term monitoring of wildlife populations, tides, or weather at French Frigate Shoals, Laysan Island, and Midway Atoll. These data sets are potentially valuable baselines, which could be informative for adaptive learning (integrating management and science) to predict, adapt, and mitigate the effects of climate change on NWHI wildlife in the future. This study provides the first quantitative vulnerability assessment for all of the low-lying NWHI, and results identify biological communities, locales, and resident endangered species of Papahanaumokuakea Marine National Monument expected to be at risk from SLR. This report is also intended as a reference for managers and conservation planners, a tool to identify and potentially reduce uncertainty, and a starting place for developing climate change monitoring priorities and future scientific studies.

Hawai'i

Assessing per- and polyfluoroalkyl substances (PFAS) in sediments and fishes in a large, urbanized estuary and the potential human health implications

The primary source of chronic exposures to per- and polyfluoroalkyl substances (PFASs) in humans is through the ingestion of contaminated foods and drinking water, with fish and other seafood being a major contributor. Nevertheless, there is scant literature on the dietary exposure to PFASs for the general United States (U.S.) population. The Tampa Bay (Florida, USA) region has the highest population density in the State and communities and their attendant support services are arrayed in an urban to semi-rural continuum from the head of the Bay to the ocean mouth. Tampa Bay supports productive recreational and commercial fisheries, providing a diverse community of species. A variety of potential PFAS sources surround Tampa Bay including airports, industry, wastewater treatment plants, fire-fighting training areas and military installations. The objective of this study is to quantify PFASs in sediment and fishes collected from Tampa Bay to further estimate human health risks from dietary exposures. Sediment ( n = 17) and fish (24 species, n = 140) were collected throughout Tampa Bay in 2020 and 2021 and analyzed for 25 PFAS compounds. Concentrations of PFASs in sediments and edible tissues of fish ranged from 36.8 to 2,990 ng kg -1 (dry weight) and 307 to 33,600 ng kg -1 (wet weight), respectively. Generally, levels were highest in Old Tampa Bay and decreased south towards the Gulf of Mexico. Profiles in both matrices were generally dominated by perfluorooctane sulfonic acid (PFOS) with variations by location. Estimated human health risks from the consumption of contaminated fish collected in Tampa Bay exceeded concentration thresholds for minimum risk levels (MRLs) and tolerable weekly intake (TWIs) values for adults and youths. Additionally, concentrations of PFOS in edible fish tissues of several recreationally important species collected in Tampa Bay exceeded consumption guideline levels established by several governmental agencies. In the current context, the elevated levels of PFAS in Tampa Bay and the exceedances of available thresholds for potential human health risks are a cause for concern and justify a more intensive examination especially for more heavily utilized species, particularly those used in subsistence-level fishing, which, as elsewhere may be significantly under documented.

Florida

Seabird vulnerability to oil: Exposure potential, sensitivity, and uncertainty in the northern Gulf of Mexico

The northern Gulf of Mexico (nGoM) is a globally important region for oil extraction and supports a diverse assemblage of marine birds. Due to their frequent contact with surface waters, diverse foraging strategies, and the ease with which oil adheres to feathers, seabirds are particularly susceptible to hydrocarbon contamination. Given the chronic and acute exposure of seabirds to oiling and a lack of studies that focus on the exposure of seabirds to oiling in sub-tropical and tropical regions, a greater understanding of the vulnerability of seabirds to oil in the nGoM appears warranted. We present an oil vulnerability index for seabirds in the nGoM tailored to the current state of knowledge using new, spatiotemporally expensive vessel-based seabird observations. We use information on the exposure and sensitivity of seabirds to oil to rank seabird vulnerability. Exposure variables characterized the potential to encounter oil and gas (O&G). Sensitivity variables characterized the potential impact of seabirds interacting with O&G and are related to life history and productivity. We also incorporated uncertainty in each variable, identifying data gaps. We found that the percent of seabirds’ habitat defined as highly suitable within 10 km of an O&G platform ranged from 0%-65% among 24 species. Though O&G platforms only overlap with 15% of highly suitable seabird habitat, overlap occurs in areas of moderate to high vulnerability of seabirds, particularly along the shelf-slope. Productivity-associated sensitivity variables were primarily responsible for creating the gradient in vulnerability scores and had greater uncertainty than exposure variables. Highly vulnerable species (e.g., Northern gannet ( Morus bassanus )) tended to have high exposure to the water surface via foraging behaviors (e.g., plunge-diving), older age at first breeding, and an extended incubating and fledging period compared to less vulnerable species (e.g., Pomarine jaeger ( Stercorarius pomarinus )). Uncertainty related to productivity could be reduced through at-colony monitoring. Strategic seabird satellite tagging could help target monitoring efforts to colonies known to use the nGoM, and continued vessel-based observations could improve habitat characterization. As offshore energy development in the nGoM continues, managers and researchers could use these vulnerability ranks to identify information gaps to prioritize research and focal species.

northern Gulf of Mexico

Coastal groundwater dynamics off Santa Barbara, California: combining geochemical tracers, electromagnetic seepmeters, and electrical resistivity

This paper presents repeat field measurements of 222 Rn and 223,224,226,228 Ra, electromagnetic seepage meter-derived advective fluxes, and multi-electrode, stationary and continuous marine resistivity surveys collected between November 2005 and April 2007 to study coastal groundwater dynamics within a marine beach in Santa Barbara, California. The study provides insight into magnitude and dynamics of submarine groundwater discharge (SGD) and associated nutrient loadings into near-shore coastal waters, where the predominant SGD drivers can be both spatially and temporally separated. Rn-222 and 223,224,226,228 Ra were utilized to quantify the total and saline contribution, respectively, of SGD. The two short-lived 224,223 Ra isotopes provided an estimate of apparent near-shore water mass age, as well as an estimate of the Ra-derived eddy diffusion coefficient, K h ( 224 Ra = 2.86 ± 0.7 m 2 s −1 ; 223 Ra = 1.32 ± 0.5 m 2 s −1 ). Because 222 Rn ( t ½ = 3.8 day) and 224 Ra ( t ½ = 3.66 day) have comparable half-lives and production terms, they were used in concert to examine respective water column removal rates. Electromagnetic seepage meters recorded the physical, bi-directional exchange across the sediment/water interface, which ranged from −6.7 to 14.5 cm day −1 , depending on the sampling period and position relative to the low tide line. Multi-day time-series 222 Rn measurements in the near-shore water column yielded total (saline + fresh) SGD rates that ranged from 3.1 ± 2.6 to 9.2 ± 0.8 cm day −1 , depending on the sampling season. Offshore 226 Ra ( t ½ = 1600 year) and 222 Rn gradients were used with the calculated K h values to determine seabed flux estimates (dpm m −2 day −1 ), which were then converted into SGD rates (7.1 and 7.9 cm day −1 , respectively). Lastly, SGD rates were used to calculate associated nutrient loads for the near-shore coastal waters off Santa Barbara. Depending on both the season and the SGD method utilized, the following SGD-derived nutrient inputs were computed (mol per day per meter of shoreline): NH 4 + = 0.06–0.29 mol day −1 m −1 ; SiO 4 = 0.22–0.29 mol day −1 m −1 ; PO 4 3− = 0.04–0.17 mol day −1 m −1 ; [NO 2 − + NO 3 − ] = 0–0.52 mol day −1 m −1 ; dissolved inorganic nitrogen (DIN) = 0.01–0.17 mol day −1 m −1 , and dissolved organic nitrogen (DON) = 0.08–0.09 mol day −1 m −1 . Compared to the ephemeral nature of fluvial and marine inputs into this region, such SGD-derived loadings can provide a sustained source of select nutrients to the coastal waters off Santa Barbara, California that should be accounted for in mass balance estimates.

California

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

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

California, Oregon, Washington

Dissolved organic carbon dynamics and fluxes in Mississippi-Atchafalaya deltaic system impacted by an extreme flood event and hurricanes: A multi-satellite approach using Sentinel-2/3 and Landsat-8/9 data

Transport of riverine and wetland-derived dissolved organic carbon (DOC) spanning tidal wetlands, estuaries, and continental shelf waters functionally connects terrestrial and aquatic carbon reservoirs, yet the magnitude and ecological significance of this variable and its spatiotemporal linkage remains uncertain for coastal deltaic regions, such as Mississippi River Delta Plain, which includes Mississippi (MR) and Atchafalaya (AR) rivers and estuaries with vast expanses of wetlands and coastal forests. We examined DOC dynamics and fluxes in this large river-dominated wetland-estuarine system for the period between 2019 and 2021 that included an extreme river flood event in 2019, two major hurricanes (Barry in 2019 and Ida in 2021), and cold front passage using an improved adaptive quasi-analytical algorithm (QAA-AD) applied to multi-satellite sensors (Sentinel 3A/B OLCI, Landsat-8/OLI and Sentinel-2A/B MSI) with varying spectral and spatial (10/30/300 m) resolutions. The DOC estimates from multi-satellite sensors in combination with water fluxes were used to assess DOC fluxes from two large rivers (MR and AR) and small channels across the delta plain. Overall, this system delivered a total of 6.7 Tg C yr -1 (1 Tg = 10 12 g) into the estuarine zone and the northern Gulf of Mexico (nGoM) during 2019. High DOC fluxes from the AR (1.3 Tg C yr -1 ) and MR (4.5 Tg C yr -1 ) were associated with the extreme flood event in 2019. Hurricanes that occurred in the study period also contributed to the wetland and estuarine DOC fluxes into continental shelf waters; for example, the passage of Hurricane Barry in July 2019, delivered over a 3-day period ~1.33 ×10 9 g DOC from Barataria Basin into the nGoM. Sentinel 2-MSI land and water classification revealed that Hurricane Ida eroded a total of 1.34×10 8 m 2 of marshes in middle Barataria Basin, converting those habitats into open water with 3.0 m inundation depth and high DOC concentrations (16.4 mg L -1 ), a potentially large DOC source to the coastal waters. Overall, storms and flood events are major sources of DOC flux that facilitate transport of upstream carbon as well as transformation of carbon in the wetlands, through the conversion of vegetated wetland to open water.

Louisiana

The United States National Climate Assessment - Alaska Technical Regional Report

The Alaskan landscape is changing, both in terms of effects of human activities as a consequence of increased population, social and economic development and their effects on the local and broad landscape; and those effects that accompany naturally occurring hazards such as volcanic eruptions, earthquakes, and tsunamis. Some of the most prevalent changes, however, are those resulting from a changing climate, with both near term and potential upcoming effects expected to continue into the future. Alaska's average annual statewide temperatures have increased by nearly 4&deg;F from 1949 to 2005, with significant spatial variability due to the large latitudinal and longitudinal expanse of the State. Increases in mean annual temperature have been greatest in the interior region, and smallest in the State's southwest coastal regions. In general, however, trends point toward increases in both minimum temperatures, and in fewer extreme cold days. Trends in precipitation are somewhat similar to those in temperature, but with more variability. On the whole, Alaska saw a 10-percent increase in precipitation from 1949 to 2005, with the greatest increases recorded in winter. The National Climate Assessment has designated two well-established scenarios developed by the Intergovernmental Panel on Climate Change (Nakicenovic and others, 2001) as a minimum set that technical and author teams considered as context in preparing portions of this assessment. These two scenarios are referred to as the Special Report on Emissions Scenarios A2 and B1 scenarios, which assume either a continuation of recent trends in fossil fuel use (A2) or a vigorous global effort to reduce fossil fuel use (B1). Temperature increases from 4 to 22&deg;F are predicted (to 2070-2099) depending on which emissions scenario (A2 or B1) is used with the least warming in southeast Alaska and the greatest in the northwest. Concomitant with temperature changes, by the end of the 21st century the growing season is expected to lengthen by 15-25 days in some areas of Alaska, with much of that corresponding with earlier spring snow melt. Future projections of precipitation (30-80 years) over Alaska show an increase across the State, with the largest changes in the northwest and smallest in the southeast. Because of increasing temperatures and growing season length, however, increased precipitation may not correspond with increased water availability, due to temperature related increased evapotranspiration. The extent of snow cover in the Northern Hemisphere has decreased by about 10 percent since the late 1960s, with stronger trends noted since the late 1980s. Alaska has experienced similar trends, with a strong decrease in snow cover extent occurring in May. When averaged across the State, the disappearance of snow in the spring has occurred from 4 to 6 days earlier per decade, and snow return in fall has occurred approximately 2 days later per decade. This change appears to be driven by climate warming rather than a decrease in winter precipitation, with average winter temperatures also increasing by about 2.5&deg;F. The extent of sea ice has been declining, as has been widely published in both national and scientific media outlets, and is projected to continue to decline during this century. The observed decline in annual sea ice minimum extent (September) has occurred more rapidly than was predicted by climate models and has been accompanied by decreases in ice thickness and in the presence of multi-year ice. This decrease was first documented by satellite imagery in the late 1970s for the Bering and Chukchi Seas, and is projected to continue, with the potential for the disappearance of summer sea ice by mid- to late century. A new phenomenon that was not reported in previous assessments is ocean acidification. Uptake of carbon dioxide (CO2) by oceans has a significant effect on marine biogeochemistry by reducing seawater pH. Ocean acidification is of particular concern in Alaska, because cold sea water absorbs CO2 more rapidly than warm water, and a decrease in sea ice extent has allowed increased sea surface exposure and more uptake of CO2 into these northern waters. Ocean acidification will likely affect the ability of organisms to produce and maintain shell material, such as aragonite or calcite (calcium carbonate minerals structured from carbonate ions), required by many shelled organism, from mollusks to corals to microscopic organisms at the base of the food chain. Direct biological effects in Alaska further along the food chain have yet to be studied and may vary among organisms. Some of the potentially most significant changes to Alaska that could result from a changing climate are the effects on the terrestrial cryosphere - particularly glaciers and permafrost. Alaskan glaciers are changing at a rapid rate, the primary driver appearing to be temperature. Statewide, glaciers lost 13 cubic miles of ice annually from the 1950s to the 1990s, and that rate doubled in the 2000s. However, like temperature and precipitation, glacier ice loss is not spatially uniform; most glaciers are losing mass, yet some are growing (for example Hubbard Glacier in southeast Alaska). Alaska glaciers with the most rapid loss are those terminating in sea water or lakes. With this increasing rate of melt, the contribution of surplus fresh water entering into the oceans from Alaska's glaciers, as well as those in neighboring British Columbia, Canada, is approximately 20 percent of that contributed by the Greenland Ice Sheet. Permafrost degradation (that is, the thawing of ice-rich soils) is currently (2012) impacting infrastructure and surface-water availability in areas of both discontinuous and continuous ground ice. Over most of the State, the permafrost is warming, with increasing temperatures broadly consistent with increasing air temperatures. On the Arctic coastal plain of Alaska, permafrost temperatures showed some cooling in the 1950s and 1960s but have been followed by a roughly 5&deg;F increase since the 1980s. Many areas in the continuous permafrost zone have seen increases in temperature in the seasonally active layer and a decrease in re-freezing rates. Changes in the discontinuous permafrost zone are initially much more observable due to the resulting thermokarst terrain (land surface formed as ice rich permafrost thaws), most notable in boreal forested areas. Climate warming in Alaska has potentially broad implications for human health and food security, especially in rural areas, as well as increased risk for injury with changing winter ice conditions. Additionally, such warming poses the potential for increasing damage to existing water and sanitation facilities and challenges for development of new facilities, especially in areas underlain by permafrost. Non-infectious and infectious diseases also are becoming an increasing concern. For example, from 1999 to 2006 there was a statistically significant increase in medical claims for insectbite reactions in five of six regions of Alaska, with the largest percentage increase occurring in the most northern areas. The availability and quality of subsistence foods, normally considered to be very healthy, may change due to changing access, changing habitats, and spoilage of meat in food storage cellars. These and other trends and potential outcomes resulting from a changing climate are further described in this report. In addition, we describe new science leadership activities that have been initiated to address and provide guidance toward conducting research aimed at making available information for policy makers and land management agencies to better understand, address, and plan for changes to the local and regional environment. This report cites data in both metric and standard units due to the contributions by numerous authors and the direct reference of their data.

Alaska

Prediction of the probability of elevated nitrate concentrations at groundwater depths used for drinking-water supply in the Puget Sound basin, Washington, 2004–19

The Puget Sound basin encompasses the 13,700-square-mile area that drains to the Puget Sound and the adjacent marine waters of Washington State. Well more than 4 million people live within the basin, with numbers continuing to increase, who rely on the basin’s natural resources including groundwater. The Puget Sound Partnership was created by a Washington State statute to implement a science-based recovery of the Puget Sound to help address impacts to these resources. As part of the recovery, the partnership developed the Puget Sound Vital Signs as measures of ecosystem health that guide the assessment of progress toward Puget Sound recovery goals. The Puget Sound Partnership Leadership Council adopted a Drinking Water Vital Sign associated with human health and quality of life, recognizing certain indicators as integral to the sustainability of Puget Sound recovery efforts. One such Vital Sign indicator was the vulnerability of groundwater throughout the aquifers of the Puget Sound basin to elevated nitrate concentrations as defined by the probability of exceeding 2 milligrams/liter (mg/L) at a specific location and well depth. The U.S. Geological Survey (USGS) led the effort to characterize groundwater vulnerability. For this study, groundwater vulnerability refers to a probability with which a contaminant applied at or near the land surface can migrate to the aquifer of interest for a given set of land-use practices. Nitrate concentration data were selected for evaluation because elevated nitrate concentrations are typically caused by anthropogenic activities and have been associated with deleterious impacts on human health. To identify groundwater vulnerability to elevated nitrate concentrations, logistic regression was used to relate anthropogenic (human associated) and natural variables to the occurrence of elevated nitrate concentrations in untreated groundwater from large public water supply system wells found within the Washington State Department of Health Sentry database. Variables that were analyzed included well depth, soil hydraulic conductivity, precipitation, population density, fertilizer application amounts, and land-use types. Statistically significant models that predicted the probabilities of groundwater nitrate concentrations greater than 2 mg/L based on the predictor variables were created for the time periods 2000–04, 2005–09, 2010–14, and 2015–19. For all time periods, well depth and a measure of the abundance of urban and agricultural land over or near the well consistently helped explain the vulnerability of the well to elevated nitrate concentrations defined as a probability of exceeding 2 mg/L of nitrate. Precipitation and (or) soil hydraulic conductivity were also important predictor variables in the models. The models for each time period were used to create maps of groundwater vulnerability at 150- and 300-foot depths throughout the Puget Sound basin. As expected, the most vulnerable locations were associated with shallower well depths and increased agriculture and urban land cover. Across all four time periods, groundwater vulnerability throughout the Puget Sound was low, with probabilities of exceeding 2 mg/L concentrations of nitrate at depths at 150 and 300 feet typically less than 50 percent. Results also found a slight decrease in probabilities of elevated nitrate concentrations throughout the basin over time. More specifically, additional statistical tests found that groundwater with probabilities of less than about 60 percent declined from 2000 to 2019 and represented more than 75 percent of the modeled Puget Sound basin aquifer. Wells with greater than 60 percent probability increased over the same time period but represented only about 25 percent of the aquifer. The maps and statistical analysis presented in the study provide valuable and informative evaluation of the vulnerability of groundwater in the Puget Sound basin to elevated nitrate concentrations. The probability maps do not represent measured nitrate concentrations in groundwater, but rather they present the probability that nitrate concentrations exceed 2 mg/L. The models and predictions from this study are a viable indicator for the Puget Sound Partnership’s Healthy Human Population—Drinking Water Vital Sign. The logistic regression modeling approach presented here benefits water managers by allowing them to assess temporal trends in a range of probabilities, explore vulnerability changes as new regional land cover and anthropogenic data are generated, and distinguish vulnerabilities at different depths within the aquifer.

Washington