Search USGS⌕ Search

SEARCH · Search USGS

Results for “One Ecosystem”

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,693 records · Page 94Linked to original sources

Home range use and movement patterns of non-native feral goats in a tropical island montane dry landscape

Advances in wildlife telemetry and remote sensing technology facilitate studies of broad-scale movements of ungulates in relation to phenological shifts in vegetation. In tropical island dry landscapes, home range use and movements of non-native feral goats ( Capra hircus ) are largely unknown, yet this information is important to help guide the conservation and restoration of some of the world’s most critically endangered ecosystems. We hypothesized that feral goats would respond to resource pulses in vegetation by traveling to areas of recent green-up. To address this hypothesis, we fitted six male and seven female feral goats with Global Positioning System (GPS) collars equipped with an Argos satellite upload link to examine goat movements in relation to the plant phenology using the Normalized Difference Vegetation Index (NDVI). Movement patterns of 50% of males and 40% of females suggested conditional movement between non-overlapping home ranges throughout the year. A shift in NDVI values corresponded with movement between primary and secondary ranges of goats that exhibited long-distance movement, suggesting that vegetation phenology as captured by NDVI is a good indicator of the habitat and movement patterns of feral goats in tropical island dry landscapes. In the context of conservation and restoration of tropical island landscapes, the results of our study identify how non-native feral goats use resources across a broad landscape to sustain their populations and facilitate invasion of native plant communities.

Hawai‘i↗

Population and habitat analyses for greater sage-grouse (Centrocercus urophasianus) in the bi-state distinct population segment—2018 update

Executive Summary The Bi-State Distinct Population Segment (Bi-State DPS) of greater sage-grouse ( Centrocercus urophasianus , hereinafter “sage-grouse”) represents a genetically distinct and geographically isolated population that straddles the border between Nevada and California. The primary threat to these sage-grouse populations is the expansion of single-leaf pinyon ( Pinus monophylla ) and Utah juniper ( Juniperus osteosperma ) into sagebrush ecosystems, which fragments and reduces population connectivity and survival. Other important threats include low water availability during brood-rearing, particularly during drought, and increased predation by common ravens ( Corvus corax ), a generalist predator often associated with anthropogenic resource subsidies. Although the Bi-State DPS occurs at high elevations relative to sage-grouse range-wide, changes in historical wildfire cycles and the conversion of native shrubs to invasive annual grasslands still threaten these populations. The Bi-State DPS has undergone multiple federal status assessments and associated litigation. For example, in October of 2013, the Bi-State DPS was proposed for listing as threatened under the Endangered Species Act of 1973 by the U.S. Fish and Wildlife Service (USFWS), then withdrawn in April 2015. The withdrawal decision was challenged, and in May 2018, a Federal district court ordered the withdrawal decision to be vacated, and USFWS was required to re-open the October 2013 listing evaluation. In response, the U.S. Geological Survey (USGS), with State and Federal collaborators, embarked on a multipronged analysis to provide current and best available science regarding population status of sage-grouse within the Bi-State DPS. Using data from a long-term monitoring program, we carried out four analytical study objectives, and here, we provide preliminary results of these analyses. First, we used integrated population modeling (IPM) to predict annual population abundance and annual finite rate of population change for the Bi-State DPS, as a whole, and for each subpopulation between 1995 and 2018. Because sage-grouse exhibit population cycles (periodic increases and decreases in abundance across approximately 6- to 10-year wavelengths), we estimated trends across three nested temporal scales that represent one (11 years), two (18 years), and three (24 years) complete population cycles. These estimates of relatively long-term averaged population change account for temporal (that is, interannual) variation. Our model predicted population abundance for the Bi-State DPS during 2018 at 3,305 individuals (2,247–4,683), with the majority occupying Bodie Hills and Long Valley. The model also predicted cyclic dynamics in abundance through time with evidence of 24-year population growth and slight trends of decline over the past 18 years. Specifically, across the Bi-State DPS as a whole, we estimated annual average at 0.99, 0.99, and 1.02 over the one, two, and three population cycles, which equated to a 10.5 percent, 16.6 percent decrease, and 60.0 percent increase in abundance over the 11-, 18-, and 24-year cycles. Estimated abundance in 2018 had not reached numbers lower than those predicted during 1995. However, we found spatial variation in population trends across the three cycles. Bodie Hills subpopulation comprised the greatest (1,521) and exhibited average annual greater than 1.0 across all periods resulting in average annual increases of 7 percent. This relatively large subpopulation has grown 5 times larger than what was predicted in 1995 while experiencing cyclical dynamics within that period. Conversely, other smaller subpopulations within the Bi-State DPS exhibited average annual equal to or less than 1.0 resulting in estimated 10-year risks of extirpation ranging from 2.0 to 76.1 percent. In general, evidence of decline among smaller subpopulations was greatest for the most recent period (2008–18) compared to a period that encompassed three full population cycles (24-year). This difference coincides with an intense period of drought that began in 2012. For comparative purposes as part of this first objective, we conducted a similar analysis for populations of sage-grouse within Nevada and California but outside the Bi-State DPS. We developed a region-wide and distance-weighted IPM using lek count from Nevada Department of Wildlife (NDOW) and California Department of Fish and Wildlife (CDFW) databases and with telemetry data collected by USGS across 12 sage-grouse subpopulations. Our models predicted similar patterns in population cycling outside the Bi-State DPS but with much stronger evidence of long-term declines across 24 years. Specifically, median averaged across each year of the 11-, 18-, and 24-year periods resulted in average annual values of 0.94, 0.97, and 0.99, respectively. These values equate to 41.0 percent, 38.5 percent, and 21.3 percent declines over the corresponding periods. Second, we used lek count data in a state-space modeling framework to compare trends in population abundance across different spatial scales (that is, leks versus Bi-State DPS). This hierarchical framework allowed us to disentangle declines associated with climate conditions as opposed to other local level factors that might signal the need for management intervention. Specifically, we identified 7 leks that were both declining and recently decoupled from larger spatial scale trends, typically governed by climatic conditions (referred to as soft or hard signals). The goal of this analysis was to provide an early warning system that might have implications for conservation actions at local scales. Third, we developed phenological (spring, summer–fall, and winter) and reproductive life stage (nesting, early brood-rearing, and late-brood rearing) based resource selection functions using various environmental covariates. We report rankings of variable importance for each season and life stage, developed maps of habitat selection indices (HSI), binned categories representing low, moderate, and high classes of quality (where any category greater than or equal to low indicated selected habitat) for each phenological season and life stage, and produced composite maps by selected phenological and reproductive stage to estimate annual habitat. Fourth, we used for each lek within the Bi-State DPS to carry out a spatial analysis that quantified substantial changes in the distribution of occupied habitat across long- (24-year) and short- (11-year) term periods. Owing to differences among available datasets, the long-term analysis primarily reflected spatial shifts among subpopulations comprising the majority of the Bi-State DPS (that is, Bodie Hills and Long Valley) while the short-term analysis also quantified changes among subpopulations along the periphery. Over long and short-term periods, the overall distribution of occupied habitat (as measured by 99 percent utilization distributions intersecting any quantified habitat) was reduced by 20,573 ha and 55,492 ha, respectively. Occupied core areas (as measured by 50 percent utilization distributions intersecting any quantified habitat) over long-term periods were solely located in Bodie Hills and Long Valley. Although nearly all subpopulations experienced contractions in occupied overall and core distribution, Bodie Hills experienced spatial expansion that occurred with concomitant spatial contraction at Long Valley over both periods. Subpopulations at the northern (Pine Nuts), central (Sagehen) and southern (White-Mountains) extents of the Bi-State DPS also experienced spatial contraction over the short-term period. These findings, coupled with those of population trends, indicate long-term patterns in redistribution of sage-grouse from Long Valley and peripheral subpopulations to Bodie Hills. That is, sage-grouse subpopulations at the periphery are declining while the largest population at the core is increasing, which could have meaningful impacts on overall metapopulation persistence. We provide evidence for loss of occupied habitat (reduced distribution) given local extirpation of subpopulations. Fifth, we calculated percentages of selected phenological, life stage, and annual habitat that each subpopulation contributed to the Bi-State DPS. We then intersected these maps with a composite estimate of occupied habitat from the fourth objective and calculated percentages of selected habitat likely occupied by sage-grouse that each subpopulation contributed to the Bi-State DPS. These values provide evidence for loss of occupied habitat and subsequent reductions in spatial distribution given reductions in abundance and, in some cases, extirpation of leks within subpopulations. Lastly, we carried out an initial in-depth analysis of selection for irrigated pastures and wet meadows during the brood-rearing stage for the Long Valley subpopulation. We chose this subpopulation because it represents a population core, representing 26.5 percent of total sage-grouse within the Bi-State DPS, and has exhibited long-term declines in abundance and distribution. This subpopulation is highly sensitive to precipitation and other factors that influence water availability. Models predicted higher use of the interior portions of irrigated pastures and wet meadows during late brood-rearing period, which represented a potentially risky use of habitat that was exacerbated during periods of low moisture (for example, drought, reduced water delivery, or both). Sage-grouse typically used edges of riparian areas and pastures, largely because the interior of these mesic areas consisted of considerably less overhead concealment cover (for example, shrubs) that likely resulted in a higher risk of mortality. We found that a lack of water delivery to pastures in the form of overwinter precipitation or diversion ditches increased the movements of sage-grouse to the interior of pastures. Although further investigation of water delivery impacts on chick survival are warrented, our initial findings regarding resource selection may explain recent declines in population growth at Long Valley.

California, Nevada↗

Management of remnant tallgrass prairie by grazing or fire: Effects on plant communities and soil properties

Tallgrass prairie is a disturbance‐dependent ecosystem that has suffered steep declines in the midwestern United States. The necessity of disturbance, typically fire or grazing, presents challenges to managers who must apply them on increasingly small and fragmented parcels. The goal of this study was to compare effects of management using cattle grazing or fire on vegetation and soil characteristics to aid managers in making decisions regarding the kind of disturbance to apply. We selected 73 sites, of which 27 were managed solely by cattle grazing and 46 solely by fire, for at least 11 yr leading up to the study. We stratified the sites by prairie type (dry, mesic, and wet) and sampled frequency of plant species on randomly placed transects, supplemented with botanist‐directed walks, and collected and composited five soil cores on a randomly selected transect within each prairie type at each site. We calculated rarefied richness and Shannon evenness from the transect data and mean coefficient of conservatism (CofC) from the total list of species. Soil samples were analyzed for texture, bulk density, total N and C, and potential net N nitrification and mineralization. A nonmetric multidimensional scaling analysis of the plant community data revealed differences in species associated with mesic and wet prairies, but no separation by management type. Similarly, none of the vegetation variables we calculated varied by management type, as determined by mixed‐effects models, but soil bulk density was 17.5% higher and total N was 22% higher on grazed sites than burned sites. Sites burned more recently had higher species richness and mean CofC, but fire was not associated with any soil variables. Sites grazed more recently had higher bulk density, total N and C, and faster N cycling rates. Overall, 28% of plant species were found exclusively in one management type or the other, but these species did not vary in mean CofC. We conclude that, at the levels of burning and grazing intensity we studied, both management approaches produce similar C storage and vegetation responses. To maintain maximum diversity across the landscape, however, both approaches are necessary.

Minnesota↗

A Review of Land-Cover Mapping Activities in Coastal Alabama and Mississippi

INTRODUCTION Land-use and land-cover (LULC) data provide important information for environmental management. Data pertaining to land-cover and land-management activities are a common requirement for spatial analyses, such as watershed modeling, climate change, and hazard assessment. In coastal areas, land development, storms, and shoreline modification amplify the need for frequent and detailed land-cover datasets. The northern Gulf of Mexico coastal area is no exception. The impact of severe storms, increases in urban area, dramatic changes in land cover, and loss of coastal-wetland habitat all indicate a vital need for reliable and comparable land-cover data. Four main attributes define a land-cover dataset: the date/time of data collection, the spatial resolution, the type of classification, and the source data. The source data are the foundation dataset used to generate LULC classification and are typically remotely sensed data, such as aerial photography or satellite imagery. These source data have a large influence on the final LULC data product, so much so that one can classify LULC datasets into two general groups: LULC data derived from aerial photography and LULC data derived from satellite imagery. The final LULC data can be converted from one format to another (for instance, vector LULC data can be converted into raster data for analysis purposes, and vice versa), but each subsequent dataset maintains the imprint of the source medium within its spatial accuracy and data features. The source data will also influence the spatial and temporal resolution, as well as the type of classification. The intended application of the LULC data typically defines the type of source data and methodology, with satellite imagery being selected for large landscapes (state-wide, national data products) and repeatability (environmental monitoring and change analysis). The coarse spatial scale and lack of refined land-use categories are typical drawbacks to satellite-based land-use classifications. Aerial photography is typically selected for smaller landscapes (watershed-basin scale), for greater definition of the land-use categories, and for increased spatial resolution. Disadvantages of using photography include time-consuming digitization, high costs for imagery collection, and lack of seasonal data. Recently, the availability of high-resolution satellite imagery has generated a new category of LULC data product. These new datasets have similar strengths to the aerial-photo-based LULC in that they possess the potential for refined definition of land-use categories and increased spatial resolution but also have the benefit of satellite-based classifications, such as repeatability for change analysis. LULC classification based on high-resolution satellite imagery is still in the early stages of development but merits greater attention because environmental-monitoring and landscape-modeling programs rely heavily on LULC data. This publication summarizes land-use and land-cover mapping activities for Alabama and Mississippi coastal areas within the U.S. Geological Survey (USGS) Northern Gulf of Mexico (NGOM) Ecosystem Change and Hazard Susceptibility Project boundaries. Existing LULC datasets will be described, as well as imagery data sources and ancillary data that may provide ground-truth or satellite training data for a forthcoming land-cover classification. Finally, potential areas for a high-resolution land-cover classification in the Alabama-Mississippi region will be identified.

Open-File Report↗

Water quality and landscape processes of four watersheds in eastern Puerto Rico

Humid tropical regions occupy about a quarter of Earth's land surface, yet they contribute a substantially higher fraction of the water, solutes, and sediment discharged to the world's oceans. Nearly half of Earth's population lives in the tropics, and development stresses can potentially harm soil resources, water quality, and water supply and in addition increase landslide and flood hazards. Owing to Puerto Rico's steep topography, low water storage capacity, and dependence on trade-wind precipitation, the island's people, ecosystems, and water supply are vulnerable to extreme weather such as hurricanes, floods, and droughts. Eastern Puerto Rico offers a natural laboratory for separating geologic and land-cover influences from regional- and global-scale influences because of its various bedrock types and the changing land cover surrounding intact, mature forest of the Luquillo Experimental Forest. Accordingly, a multiyear assessment of hydrological and biogeochemical processes was designed to develop an understanding of the effects of these differences on local climate, streamflow, water quality, and ecosystems, and to form the basis for a long-term and event-based program of climate and hydrologic monitoring. Because infrequent, large storms play a major role in this landscape, we focused on high-runoff events, sampling 263 storms, including all major hurricanes from 1991 through 2005. The largest storms have profound geomorphic consequences, such as landslides, debris flows, deep gullying on deforested lands, excavation and suspension of sediment in stream channels, and delivery of a substantial fraction of annual stream sediment load. Large storms sometimes entrain ocean foam and spray causing high concentrations of seasalt-derived constituents in stream waters during the storm. Past deforestation and agricultural activities in the Cayaguás and Canóvanas watersheds accelerated erosion and soil loss, and this material continues to be remobilized during large storms. Nearly 5,000 routine and event samples were analyzed for parameters that allow determination of denudation rates based on suspended and dissolved loads; 860 of these samples were analyzed for a comprehensive suite of chemical constituents. The rivers studied are generally similar in water-quality characteristics, and windward or leeward aspect appears to exert a stronger influence on water quality than geology or land cover. Of samples analyzed for comprehensive chemistry and for sediment, 543 were collected at runoff rates greater than 1 millimeter per hour, 256 at rates exceeding 10 millimeters per hour, and 3 at rates exceeding 90 millimeters per hour. Streams have rarely been sampled during events with such high runoff rates. Rates of physical and chemical weathering are especially high, and physical denudation rates, forested watersheds included, are considerably greater than is expected for a steady-state system. The elevated physical erosion drives an increased particulate organic carbon flux, one that is large, important to the carbon cycle, and sustainable, because soil-carbon regeneration is rapid. The 15-year Water, Energy, and Biogeochemical Budget dataset, which includes discharge, field parameters, suspended sediment, major cations and anions, and nutrients, is available from the U.S. Geological Survey's National Water Information System (http://waterdata.usgs.gov/nwis). The dataset provides a baseline for characterizing future environmental change and will improve our understanding of the interdependencies of land, water, and biological resources and their responses to changes in climate and land use. Because eastern Puerto Rico resembles many tropical regions in terms of geology and patterns of development, implications from this study are transferable to other tropical regions facing deforestation, rapid land-use change, and climate change.

Professional Paper↗

Population and spatial dynamics of desert bighorn sheep in Grand Canyon during an outbreak of respiratory pneumonia

Introduction: Terrestrial species in riverine ecosystems face unique constraints leading to diverging patterns of population structure, connectivity, and disease dynamics. Desert bighorn sheep ( Ovis canadensis nelsoni ) in Grand Canyon National Park, a large native population in the southwestern USA, offer a unique opportunity to evaluate population patterns and processes in a remote riverine system with ongoing anthropogenic impacts. We integrated non-invasive, invasive, and citizen-science methods to address questions on abundance, distribution, disease status, genetic structure, and habitat fragmentation. Methods: We compiled bighorn sightings collected during river trips by park staff, commercial guides, and private citizens from 2000–2018 and captured bighorn in 2010–2016 to deploy GPS collars and test for disease. From 2011–2015, we non-invasively collected fecal samples and genotyped them at 9–16 microsatellite loci for individual identification and genetic structure. We used assignment tests to evaluate genetic structure and identify subpopulations, then estimated gene flow and recent migration to evaluate fragmentation. We used spatial capture-recapture to estimate annual population size, distribution, and trends after accounting for spatial variation in detection with a resource selection function model. Results and discussion: From 2010–2018, 3,176 sightings of bighorn were reported, with sightings of 56–145 bighorn annually on formal surveys. From 2012–2016, bighorn exhibiting signs of respiratory disease were observed along the river throughout the park. Of 25 captured individuals, 56% were infected by Mycoplasma ovipneumoniae , a key respiratory pathogen, and 81% were recently exposed. Pellet sampling for population estimation from 2011–2015 yielded 1,250 genotypes and 453 individuals. We detected 6 genetic clusters that exhibited mild to moderate genetic structure ( F ST 0.022–0.126). The river, distance, and likely topography restricted recent gene flow, but we detected cross-river movements in one section via genetic recaptures, no subpopulation appeared completely isolated, and genetic diversity was among the highest reported. Recolonization of one large stretch of currently empty habitat appears limited by the constrained topology of this system. Annual population estimates ranged 536–552 (95% CrI range 451–647), lamb:ewe ratios varied, and no significant population decline was detected. We provide a multi-method sampling framework useful for sampling other wildlife in remote riverine systems.

Arizona↗

Evaluating the role of land cover and climate uncertainties in computing gross primary production in Hawaiian Island ecosystems

Gross primary production (GPP) is the Earth’s largest carbon flux into the terrestrial biosphere and plays a critical role in regulating atmospheric chemistry and global climate. The Moderate Resolution Imaging Spectrometer (MODIS)-MOD17 data product is a widely used remote sensing-based model that provides global estimates of spatiotemporal trends in GPP. When the MOD17 algorithm is applied to regional scale heterogeneous landscapes, input data from coarse resolution land cover and climate products may increase uncertainty in GPP estimates, especially in high productivity tropical ecosystems. We examined the influence of using locally specific land cover and high-resolution local climate input data on MOD17 estimates of GPP for the State of Hawaii, a heterogeneous and discontinuous tropical landscape. Replacing the global land cover data input product (MOD12Q1) with Hawaii-specific land cover data reduced statewide GPP estimates by ~8%, primarily because the Hawaii-specific land cover map had less vegetated land area compared to the global land cover product. Replacing coarse resolution GMAO climate data with Hawaii-specific high-resolution climate data also reduced statewide GPP estimates by ~8% because of the higher spatial variability of photosynthetically active radiation (PAR) in the Hawaii-specific climate data. The combined use of both Hawaii-specific land cover and high-resolution Hawaii climate data inputs reduced statewide GPP by ~16%, suggesting equal and independent influence on MOD17 GPP estimates. Our sensitivity analyses within a heterogeneous tropical landscape suggest that refined global land cover and climate data sets may contribute to an enhanced MOD17 product at a variety of spatial scales.

Hawaii↗

Vulnerability of Pacific salmon to invasion of northern pike (Esox lucius) in Southcentral Alaska

The relentless role of invasive species in the extinction of native biota requires predictions of ecosystem vulnerability to inform proactive management strategies. The worldwide invasion and range expansion of predatory northern pike ( Esox lucius ) has been linked to the decline of native fishes and tools are needed to predict the vulnerability of habitats to invasion over broad geographic scales. To address this need, we coupled an intrinsic potential habitat modelling approach with a Bayesian network to evaluate the vulnerability of five culturally and economically vital species of Pacific salmon ( Oncorhynchus spp.) to invasion by northern pike. This study was conducted along 22,875 stream km in the Southcentral region of Alaska, USA. Pink salmon ( O . gorbuscha ) were the most vulnerable species, with 15.2% (2,458 km) of their calculated extent identified as “highly” vulnerable, followed closely by chum salmon ( O . keta , 14.8%; 2,557 km) and coho salmon ( O . kisutch , 14.7%; 2,536 km). Moreover, all five Pacific salmon species were highly vulnerable in 1,001 stream km of shared habitat. This simple to implement, adaptable, and cost-effective framework will allow prioritizing habitats for early detection and monitoring of invading northern pike.

Alaska↗

Factors affecting long-term trends in surface-water quality in the Gwynns Falls watershed, Baltimore City and County, Maryland, 1998–2016

Factors affecting water-quality trends in urban streams are not well understood, despite current regulatory requirements and considerable ongoing investments in gray and green infrastructure. To address this gap, long-term water-quality trends and factors affecting these trends were examined in the Gwynns Falls, Maryland, watershed during 1998–2016 in cooperation with Blue Water Baltimore. Data on water-quality constituents and potential factors of influence were obtained from multiple sources and compiled for analysis, with a focus on data collected as part of the National Science Foundation funded Long-Term Ecological Research project, the Baltimore Ecosystem Study. Variability in climate (specifically, precipitation) and land cover can overwhelm actions taken to improve water quality and can present challenges for meeting regulatory goals. Analysis of land cover during 2001–11 in the Gwynns Falls watershed indicated minimal change during the study time frame; therefore, land-cover change is likely not a factor affecting trends in water quality. However, a modest increase in annual precipitation and a significant increase in winter precipitation were apparent in the region. A higher proportion of runoff producing storms was observed in the winter and a lower proportion in the summer, indicating that climate change may affect water quality in the watershed. The increase in precipitation was not reflected in annual or seasonal trends of streamflow in the watershed. Nonetheless, these precipitation changes may exacerbate the inflow and infiltration of water to gray infrastructure and reduce the effectiveness of green infrastructure. For streamflow and most water-quality constituents examined, no discernable trends were noted over the timeframe examined. Despite the increases in precipitation, no trends were observed for annual or seasonal discharge at the various sites within the study area. In some locations, nitrate, phosphate, and total nitrogen show downward trends, and total phosphorus and chloride show upward trends. Sanitary sewer overflows (gray infrastructure) and best management practices (green infrastructure) were identified as factors affecting water-quality change. The duration of sanitary sewer overflows was positively correlated with annual loads of nutrients and bacteria, and the drainage area of best management practices was negatively correlated with annual loads of phosphate and sulfate. Results of the study indicate that continued investments in gray and green infrastructure are necessary for urban water-quality improvement. Although this outcome is not unexpected, long-term datasets such as the one used in this study, allow the effects of gray and green infrastructures to be quantified. Results of this study have implications for the Gwynns Falls watershed and its residents and Baltimore City and County managers. Moreover, outcomes are relevant to other watersheds in the metropolitan region that do not have the same long-term dataset. Further, this study has established a framework for ongoing statistical analysis of primary factors affecting urban water-quality trends as regulatory programs mature.

Maryland↗

Forecasting tidal marsh elevation and habitat change through fusion of Earth observations and a process model

Reducing uncertainty in data inputs at relevant spatial scales can improve tidal marsh forecasting models, and their usefulness in coastal climate change adaptation decisions. The Marsh Equilibrium Model (MEM), a one-dimensional mechanistic elevation model, incorporates feedbacks of organic and inorganic inputs to project elevations under sea-level rise scenarios. We tested the feasibility of deriving two key MEM inputs—average annual suspended sediment concentration (SSC) and aboveground peak biomass—from remote sensing data in order to apply MEM across a broader geographic region. We analyzed the precision and representativeness (spatial distribution) of these remote sensing inputs to improve understanding of our study region, a brackish tidal marsh in San Francisco Bay, and to test the applicable spatial extent for coastal modeling. We compared biomass and SSC models derived from Landsat 8, DigitalGlobe WorldView-2, and hyperspectral airborne imagery. Landsat 8-derived inputs were evaluated in a MEM sensitivity analysis. Biomass models were comparable although peak biomass from Landsat 8 best matched field-measured values. The Portable Remote Imaging Spectrometer SSC model was most accurate, although a Landsat 8 time series provided annual average SSC estimates. Landsat 8-measured peak biomass values were randomly distributed, and annual average SSC (30 mg/L) was well represented in the main channels (IQR: 29–32 mg/L), illustrating the suitability of these inputs across the model domain. Trend response surface analysis identified significant diversion between field and remote sensing-based model runs at 60 yr due to model sensitivity at the marsh edge (80–140 cm NAVD88), although at 100 yr, elevation forecasts differed less than 10 cm across 97% of the marsh surface (150–200 cm NAVD88). Results demonstrate the utility of Landsat 8 for landscape-scale tidal marsh elevation projections due to its comparable performance with the other sensors, temporal frequency, and cost. Integration of remote sensing data with MEM should advance regional projections of marsh vegetation change by better parameterizing MEM inputs spatially. Improving information for coastal modeling will support planning for ecosystem services, including habitat, carbon storage, and flood protection.

California↗

The risk of rodent introductions from shipwrecks to seabirds on Aleutian and Bering Sea islands

Accidental introductions of rodents present one of the greatest threats to indigenous island biota, especially seabirds. On uninhabited remote islands, such introductions are likely to come from shipwrecks. Here we use a comprehensive database of shipwrecks in Western Alaska to model the frequency of shipwrecks per Aleutian and Bering Sea island, taken as a proxy for the threat of rodent introductions, using physical variables, and the intensity of nearby fishing traffic and activity as predictors. Using data spanning from 1950 to 2013, we found that shipwrecks were particularly common in the 1980s to early 2000s, with a major peak in wrecks during the late 1980s. Amount of fishing activity within 5 km of an island was the strongest predictor of shipwrecks, followed by the strength of tidal currents and density of large-vessel traffic. Islands with the highest frequency of shipwrecks are all in the eastern Aleutians, including Unimak, Unalaska, and Akun Islands. By contrast, the largest seabird colonies are in the western Aleutian and Pribilof Islands, including Buldir, Kiska, and Saint George islands. Multiplying the frequency of a shipwreck by the number of seabirds breeding per island provides a measure of risk. The risk of rodent introductions from shipwrecks to seabirds was then greatest for Saint George (Bering Sea), Buldir (Western Aleutians) and Saint Matthew islands (Bering Sea). Keeping these high-risk islands rodent free would maintain their high a conservation value. Most islands with a high predicted frequency of shipwrecks already have established rodent populations and therefore few remaining seabirds. Of those islands with established rodent populations, Attu and Kiska Islands would make suitable targets for eradication, given their relatively low expected frequency of shipwrecks for their size. Further improvements in rat prevention on vessels and shipping safety would benefit the economy, human health and safety, and to the long-term conservation of island ecosystems.

Alaska↗

Evaluating establishment success of non-native fishes introduced to inland aquatic habitats of tropical Pacific islands

An information-theoretic approach was used to evaluate non-native freshwater fish species introduced to insular habitats of Hawaii and Guam comparing successful establishments vs. failures. Since the late 1800s, as many as 81 non-native freshwater fish species have been recorded as introduced to Hawaii and Guam (combined) and 50 (62%) of these are documented as having one or more established populations. We examined eleven independent variables to investigate establishment success by creating 21 a priori logistic regression models ranked using Akaike's Information Criterion adjusted for small sample size. An additional eight post-hoc models were included that comprised the best a priori model and various combinations of individual variables. The best overall model of establishment probability included effects of taxonomic affinity (family membership), prior establishment success on other tropical islands, and hypoxia tolerance. Establishment success in Hawaii and Guam was highest for those species established on many other islands, and according to our best model air-breathing fishes were more likely to become established. Six fish families, each with from three to 18 species introduced to Pacific islands, were highly successful at establishment: Cichlidae (16 established of 18 species introduced), Poeciliidae (seven of eight), Cyprinidae (four of seven), Centrarchidae (four of four), Clariidae (three of three), and Loricariidae (three of four). Those that successfully established include both small and moderately large-bodied taxa, while representing a diverse array of other morphological and life-history traits. Pathways and motives associated with fish introductions in the Pacific have been linked to desires to develop aquaculture, enhance wild stocks of food, sport, and bait fishes, for use as biological control agents, or are linked to the ornamental fish trade. We found that many established species were introduced via multiple pathways (up to eight) and our analyses suggest that the combination of prior establishment success on other tropical islands and presence of non-native fishes in multiple pathways was indicative of high propagule pressure. Our study results and conclusions on Pacific tropical island introductions are in general agreement with previous studies on non-native freshwater fishes in other regions of the world and similar to observations in continental ecosystems and temperate zones.

Hawaii↗

Lessons learned from wetlands research at the Cottonwood Lake Study Area, Stutsman County, North Dakota, 1967–2021

Depressional wetlands in the Prairie Pothole Region of North America have a long history of investigation owing to their importance in maintaining migratory-bird populations, especially waterfowl. One area of particularly intensive study is the Cottonwood Lake study area in Stutsman County, North Dakota. Studies at the Cottonwood Lake study area began in 1967 and continue through the present (2022). During this period of scientific discovery, meteorological conditions at the Cottonwood Lake study area varied greatly and included one of the most severe droughts of the 20th century and one of the wettest periods in the past 500 years. Persistent wet conditions that began in 1993 have contributed to state changes in many of the study area’s larger wetlands to lake-like conditions, whereas the smaller wetlands returned to seasonally ponded conditions during relatively dry years interspersed within the longer-term wet period. Additionally, some nonwetland areas of the study area developed wetland plant, hydrology, and soil characteristics during the 1993-to-present (2022) wet period. The persistently high stages of water in the larger wetlands since 1993 contributed to a buildup of dissolved solids and increases in salinity with time following an initial decrease in salinity caused by the dilution of dissolved solids within a larger volume of water. During 2021, drought conditions similar to the 1988 to 1992 period may develop if conditions persist. However, meteorological changes during the past 30 years have persisted long enough to be considered a change in climate conditions at the study area and, if such wet conditions continue, would represent a change from conditions that occurred in the past two millennia. During the period of study covered in this report (1967–2021), biotic communities responded in a variety of ways to subtle and marked changes in ponded-water depths, permanence, and salinity among the different wetland types in the study area. This report provides background information on the Cottonwood Lake study area and its context within the Prairie Pothole Region, documents techniques used to quantify environmental conditions and biotic communities, describes major trends that have been observed, presents significant findings as “lessons learned,” discusses recent modeling advances, and highlights key messages to managers. The Wetland Continuum concept was used as a framework to place research findings within an ecological context and to highlight the dynamic nature of prairie-pothole wetland ecosystems.

North Dakota↗

Report on the reptiles of Upland Savai’i

The reptile team conducted a 21 kilometre transect from the coast east of Asau to the uplands ending near Mauga Silisili at over 1720 m elevation. This transect covered the main habitats on Savai’i and allowed the team to determine where various reptile species and invasive species occurred across this elevational gradient. No previous reptile research had taken place on Savai’i above the elevation of A’opo Village. Limited sampling was also done around the Forestry Station in Asau. The team detected 11 species of lizards during these surveys, which is the majority of species known from Samoa. Noticeably absent was the Pacific black skink (Emoia nigra), which is a dominant element of the Samoan lizard fauna. Also no individuals of the Pacific boa (Candoia bibroni) were detected despite the concentrated effort spent looking for them. One boa was detected by the avifauna team at their Site 1, by the TV tower on a log in a marsh. The invasive house gecko (Hemidactylus frenatus) was also not detected along the main transect, but was the most abundant gecko on buildings in Asau. No reptiles were found above 1320 m elevation and most species were found significantly below there. Snakeeyed skinks (Cryptoblepharus poecilopleurus) were detected on Savai’i for the first time at Asau Getaway Resort then above the sawmill on the way to A’opo. Since western Savai’i is so poorly known for reptiles, this is the first time many of these species were recorded from this part of the island. Surveys for invasive species detected mammalian species (cats, rats, and pigs) and invertebrate species (Yellow Crazy ants and Big-headed ants). The mammals were found at various sites along the transect, including high elevations. The ants were found at lower elevations along the transect, but the Yellow Crazy ants appear to be irrupting currently on Savai’i and were swamping our traps from sea level to 500 m elevation. The low elevation lizard occurrences from sea level to 500 m appeared impaired by the invasive ants. Although habitat looked good in many places along the transect, certain species were rare or absent when the invasive ants were present, e.g. the Samoan skink (Emoia samoensis) which only occurred at elevations higher than the ants, whereas elsewhere in its range it occurs down to sea level. Currently the uplands over 500 m are free of invasive ants. We know from Hawai’i that invasive ants occur to over 2,000 m elevation and are ecologically very destructive to native flora and fauna. There is an immediate need to try to stop this upward ant invasion to protect this at-risk ecosystem, and to study the intact system now prior to an invasion.

Report↗

Potential climate-induced runoff changes and associated uncertainty in four Pacific Northwest estuaries

As part of a larger investigation into potential effects of climate change on estuarine habitats in the Pacific Northwest, we estimated changes in freshwater inputs into four estuaries: Coquille River estuary, South Slough of Coos Bay, and Yaquina Bay in Oregon, and Willapa Bay in Washington. We used the U.S. Geological Survey's Precipitation Runoff Modeling System (PRMS) to model watershed hydrological processes under current and future climatic conditions. This model allowed us to explore possible shifts in coastal hydrologic regimes at a range of spatial scales. All modeled watersheds are located in rainfall-dominated coastal areas with relatively insignificant base flow inputs, and their areas vary from 74.3 to 2,747.6 square kilometers. The watersheds also vary in mean elevation, ranging from 147 meters in the Willapa to 1,179 meters in the Coquille. The latitudes of watershed centroids range from 43.037 degrees north latitude in the Coquille River estuary to 46.629 degrees north latitude in Willapa Bay. We calibrated model parameters using historical climate grid data downscaled to one-sixteenth of a degree by the Climate Impacts Group, and historical runoff from sub-watersheds or neighboring watersheds. Nash Sutcliffe efficiency values for daily flows in calibration sub-watersheds ranged from 0.71 to 0.89. After calibration, we forced the PRMS models with four North American Regional Climate Change Assessment Program climate models: Canadian Regional Climate Model-(National Center for Atmospheric Research) Community Climate System Model version 3, Canadian Regional Climate Model-Canadian Global Climate Model version 3, Hadley Regional Model version 3-Hadley Centre Climate Model version 3, and Regional Climate Model-Canadian Global Climate Model version 3. These are global climate models (GCMs) downscaled with regional climate models that are embedded within the GCMs, and all use the A2 carbon emission scenario developed by the Intergovernmental Panel on Climate Change. With these climate-forcing outputs, we derived the mean change in flow from the period encompassing the 1980s (1971-1995) to the period encompassing the 2050s (2041-2065). Specifically, we calculated percent change in mean monthly flow rate, coefficient of variation, top 5 percent of flow, and 7-day low flow. The trends with the most agreement among climate models and among watersheds were increases in autumn mean monthly flows, especially in October and November, decreases in summer monthly mean flow, and increases in the top 5 percent of flow. We also estimated variance in PRMS outputs owing to parameter uncertainty and the selection of climate model using Latin hypercube sampling. This analysis showed that PRMS low-flow simulations are more uncertain than medium or high flow simulations, and that variation among climate models was a larger source of uncertainty than the hydrological model parameters. These results improve our understanding of how climate change may affect the saltwater-freshwater balance in Pacific Northwest estuaries, with implications for their sensitive ecosystems.

Oregon;Washington↗

Evaluating the negative effect of benthic egg predators on bloater recruitment in northern Lake Michigan

As the only extant deepwater cisco in Lake Michigan, bloater is currently at record low levels of abundance. Several mechanisms to regulate their recruitment have been proposed, including skewed sex ratios, predation on their larvae by adult alewife, and climatic factors during early life history stages, but none has unequivocal support. In this research, we evaluated an alternative mechanism of egg predation that was supported by an inverse relationship between bloater recruitment and biomass of slimy sculpin , which are known to be effective egg predators. To that end, we used a combination of field sampling, laboratory experiments, and modeling to estimate the proportion of bloater eggs consumed by sculpins each year between 1973 and 2008. Monthly field sampling between January through May 2009-2010 (when bloater eggs were incubating) offshore of Frankfort (Michigan), Sturgeon Bay (Wisconsin), Two Rivers (Wisconsin), and Muskegon (Michigan) provided benthivore diets for subsequent laboratory processing. Identification and enumeration of stomach contents and subsequent genetic analyses of eggs revealed that the mean proportion of bloater eggs in slimy sculpin diets (N = 1016) equaled 0.04. Bloater eggs also were consumed by deepwater sculpins (N = 699) at a slightly lower mean proportion (0.02), and only one round goby diet among 552 enumerated revealed a bloater egg. Based on the diet results, we developed daily ration models to estimate consumption for both deepwater and slimy sculpins . We conducted feeding experiments to estimate gastric evacuation (GEVAC) for water temperatures ranging 2-5 °C, similar to those observed during egg incubation. GEVAC rates equaled 0.0115/ h for slimy sculpin and 0.0147/h for deepwater sculpin , and did not vary between 2.7 and 5.1 °C for either species or between prey types ( Mysis relicta and fish eggs) for slimy sculpin . Index of fullness [(g prey/g fish weight)100%] was estimated from sculpins sampled in bottom trawls in the same seasons and years as the diets, and varied with fish size (averaging 1.93% and 1.85% for slimy and deepwater sculpins , respectively). Estimates of daily consumption ranged from 0.2-0.8% of sculpin body weight. Annual estimates of bloater egg consumption predicted higher values for deepwater sculpin than slimy sculpin between 1973 and 2005. This pattern was reversed in 2006, 2008, 2009, 2010 as slimy sculpin abundance increased while that of deepwater sculpin declined. The sum of sculpin consumption of bloater eggs exceeded 25% of bloater population egg production early (1975-1980) and late (2008-2010) in the time series. Despite the strong field pattern implicating egg predation by slimy sculpin , our consumption models failed to fully support this hypothesis. In particular, our results were unable to explain why bloater recruitment was relatively poor during 1995-2005 when the proportion of bloater eggs consumed was very low ( < 0.06). The results did, however, demonstrate that bloater recruitment was consistently poor when the proportion of eggs consumed was relatively high. In conclusion, consumption by native benthivores can be a contributing factor to poor recruitment of bloater, especially when slimy sculpin reach high levels of abundance. This result exemplifies the importance of ecosystem-based fishery management, given that the maintenance of healthy lake trout populations in the Great Lakes should control the abundance of slimy sculpin egg predators. In addition, future research will be required to fully understand the primary bottleneck to bloater recruitment in Lake Michigan so that efforts to stock and restore bloater in Lake Ontario have a greater probability of resulting in naturalized and sustainable populations.

Lake Michigan↗

Section 5: Remote sensing of vegetation in the riparian corridor of the Colorado River’s delta 2013-2018

This remote sensing section is based on Nagler et al. (in preparation for the journal Hydrological Processes) and is a summary of the USGS preliminary findings to date. This report documents the changes in green foliage density (greenness) as measured by satellite vegetation index (VI) data and corresponding evapotranspiration (ET) in the riparian corridor of the Colorado River delta associated with the Minutes 319 and 323 environmental water deliveries using time-series data from 2013 through 2018. The report focuses on what happened only within the riparian corridor’s seven reaches since the 2014 flows, and despite being a continuation of measuring greenness and ET after the 2017 end of Minute 319, this study continued the tracking of these two variables, greenness and ET, in these original riparian corridor focal areas. Two spatial scales are used here: (1) Landsat satellite imagery at 30 m pixels and (2) the EOS-1 satellite sensor the Moderate Resolution Imaging Spectrometer (MODIS) with a resolution of 250 m pixels. The focal period includes 2013 (prepulse flow) and the years 2014-2018, with a focus on imagery collected from the Summer growing seasons 2014 through 2018 (one-year, pre-pulse and several post-pulse years, respectively). This report re-creates the 2013-2017 Landsat-based results from Jarchow et al. (2017a, b) by using the same region of interest (ROI). The report now provides revised and re-created results using all new imagery acquisition and processing techniques, as well as extraction code, created by the Vegetation Index and Phenology (VIP) Lab of the Biosystems Engineering Department of the University of Arizona (UofA). In 2018, methods employed by the VIP lab (and not ArcGIS) were used. ArcGIS was only used in the newly processed data to display the final difference maps. The entire spatial tile data from NASA was downloaded and processed at the VIP Lab using satellite imagery at two resolutions: 250 m MODIS and 30 m Landsat using three sensors, Landsat 5, Landsat 7 ETM+ and Landsat 8 Operational Land Imager (OLI), with added scenes for each year based on new clear atmosphere requirements. The VIP lab clipped the river boundary and seven riparian reaches from the previously existing ROI used in Jarchow et al. (2017 a, b) for the analyses done under Minute 319. The NASA image datasets for this riparian corridor ROI in seven reaches were re-processed to produce additional vegetation index (VI) information for years 2013 to 2018 for this report. At the same time, the report acquired and processed imagery from 2000- 2018 (data outside the scope of this report and data not shown here). The additional VIs (NDVI, scaled NDVI, EVI, EVI2) were analyzed so that new assessments of greenness and ET could be produced from the imagery datasets following methods in Nagler et al. (2013). These VI choices were based on previous performance comparisons between biophysical ground-based data and radiometric satellite-based data collected from this riparian ecosystem (Nagler et al., 2001) as well as performance related to ET estimation (Nagler et al., 2005a, b) and current advancements in VIs such as EVI2.

Colorado River Delta↗

Detecting significant change in stream benthic macroinvertebrate communities in wilderness areas

A major challenge in the biological monitoring of stream ecosystems in protected wilderness areas is discerning whether temporal changes in community structure are significantly outside of a reference condition that represents natural or acceptable annual variation in population cycles. Otherwise sites could erroneously be classified as impaired. Long-term datasets are essential for understanding these trends, to ascertain whether any changes in community structure significantly beyond the reference condition are permanent shifts or with time move back to within previous limits. To this end, we searched for long-term (>8 years) quantitative data sets of macroinvertebrate communities in wadeable rivers collected by similar methods and time of year in protected wilderness areas with minimal anthropogenic disturbance. Four geographic areas with datasets that met these criteria in the USA were identified, namely: McLaughlin Nature Reserve in California (1 stream), Great Smoky Mountains National Park in Tennesse-North Carolina (14 streams), Wind River Wilderness Areas in Wyoming (3 streams) and Denali National Park and Preserve in Alaska (6 streams). Two statistical approaches were applied: Taxonomic Distinctness (TD) to describe changes in diversity over time and non-metric multidimensional scaling (MDS) to describe changes over time in community persistence (Jaccards Index) and community stability (Bray&ndash;Curtis Index). Control charts were used to determine if years in MDS plots were significantly outside a reference condition. For Hunting Creek, TD showed three years outside natural variation which could be attributed to severe hydrological events but years outside the natural-variation funnel at sites in other geographical areas were inconsistent and could not be explained by environmental variables. TD identified simulated severe pollutant events which caused the removal of entire invertebrate assemblages but not simulated water temperature shifts. Within a region, both MDS analyses typically identified similar years as exceeding reference condition variation, illustrating the utility of the approach for identifying wider spatial scale effects that influence more than one stream. MDS responded to both simulated water temperature stress and a pollutant event, and generally outlying years on MDS plots could be explained by environmental variables, particularly higher precipitation. Multivariate control charts successfully identified whether shifts in community structure identified by MDS were significant and whether the shift represented a press disturbance (long-term change) or a pulse disturbance. We consider a combination of TD and MDS with control charts to be a potentially powerful tool for determining years significantly outside of a reference condition variation.

Alaska, California, North Carolina, Tennesse, Wyom↗