USGS Science⌕ Search

SEARCH · USGS Science

Results for “Risk Analysis”

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 937 records · Page 52Linked to original sources

St. Louis area earthquake hazards mapping project; seismic and liquefaction hazard maps

We present probabilistic and deterministic seismic and liquefaction hazard maps for the densely populated St. Louis metropolitan area that account for the expected effects of surficial geology on earthquake ground shaking. Hazard calculations were based on a map grid of 0.005°, or about every 500 m, and are thus higher in resolution than any earlier studies. To estimate ground motions at the surface of the model (e.g., site amplification), we used a new detailed near‐surface shear‐wave velocity model in a 1D equivalent‐linear response analysis. When compared with the 2014 U.S. Geological Survey (USGS) National Seismic Hazard Model, which uses a uniform firm‐rock‐site condition, the new probabilistic seismic‐hazard estimates document much more variability. Hazard levels for upland sites (consisting of bedrock and weathered bedrock overlain by loess‐covered till and drift deposits), show up to twice the ground‐motion values for peak ground acceleration (PGA), and similar ground‐motion values for 1.0 s spectral acceleration (SA). Probabilistic ground‐motion levels for lowland alluvial floodplain sites (generally the 20–40‐m‐thick modern Mississippi and Missouri River floodplain deposits overlying bedrock) exhibit up to twice the ground‐motion levels for PGA, and up to three times the ground‐motion levels for 1.0 s SA. Liquefaction probability curves were developed from available standard penetration test data assuming typical lowland and upland water table levels. A simplified liquefaction hazard map was created from the 5%‐in‐50‐year probabilistic ground‐shaking model. The liquefaction hazard ranges from low (<40% of area expected to liquefy) in the uplands to severe (>60% of area expected to liquefy) in the lowlands. Because many transportation routes, power and gas transmission lines, and population centers exist in or on the highly susceptible lowland alluvium, these areas in the St. Louis region are at significant potential risk from seismically induced liquefaction and associated ground deformation

Alabama, Arkansas, Illinois, Indiana, Iowa, Kentuc↗

Female persistence during toxicant treatment predicts survival probability of offspring in invasive brown treesnakes (Boiga irregularis)

Assessing the long-term efficacy of control methods is a critical component of invasive species management. For example, if traits related to control have significant heritability or are influenced by maternal effects, control methods may lose efficacy over time. The potential for these effects can be evaluated via parent/offspring survival analysis, which concomitantly recasts adaptive management as an evolutionary force for invasive species. However, difficulties can arise when the life history of an invasive is cryptic, precluding direct observations of familial relationships. Genomic pedigree reconstruction can facilitate such analyses by assigning offspring to parents in invasive species for which mating and reproduction are difficult to study. Here, we use genomic pedigree reconstruction to quantify parental longevity and probability of offspring survival for brown treesnakes ( Boiga irregularis ) on Guam in a landscape treated with toxic baits simulating application via an aerial delivery system (ADS). To do so, we used 398 single nucleotide polymorphisms (SNPs) to update an existing multi-generational genomic pedigree for a geographically-closed population of brown treesnakes. This facilitated assignment of parents to juveniles born during three consecutive years of toxic bait application under a simulated aerial treatment program (N = 72). We found that the offspring of dams that persisted until the end of the study displayed greater survival probability (cox proportional hazard model, P < 0.001), yet there was no such effect for sires. This sex-specific relationship between parental longevity and offspring survival indicates that heritability of traits contributing to resistance to ADS is unlikely, but it supports a role of maternal effects that could undermine ADS. Our study identifies potential risks associated with control efforts and also highlights the utility of parent-offspring survival analyses informed by genomic pedigree reconstruction as a tool for adaptive management.

Global Ecology and Conservation↗

Gross alpha-particle activity and high 226Ra concentrations do not correspond with high 210Po in the Atlantic and Gulf Coastal Plain aquifers of the United States

210 Po, which is of human-health concern based on lifetime ingestion cancer risk, is indirectly regulated in drinking water through the U.S. Environmental Protection Agency’s gross alpha-particle activity (GAPA) maximum contaminant level of 15 pCi/L (picocuries per liter). This regulation requires independent measurement of 226 Ra for samples exceeding the GAPA screening level of 5 pCi/L. There is no such requirement for 210 Po. Co-occurrence of 226 Ra and 210 Po, alpha-emitting 238 U-decay-series progeny, might be helpful in locating high- 210 Po waters but is unverified. Relations among 210 Po, 226 Ra, and GAPA evaluated for samples from 257 public-supply wells from Coastal Plain aquifers showed that concentrations of 226 Ra correlated with GAPA but neither correlated with 210 Po concentrations. The highest concentrations of 226 Ra and 210 Po were found under differing geochemical conditions. The highest 226 Ra occurred in low-pH oxidizing waters and in neutral-pH reducing waters, where geochemical conditions render Fe–Mn-hydroxide sorbents inefficient. 210 Po was highest (10.1 pCi/L) in reducing waters with high pH (>7.5, which results from progressive cation exchange), where 226 Ra was lowest─exchanged to clay minerals. Because 226 Ra and 210 Po did not co-occur, the GAPA screening might not be protective for 210 Po. Independent 210 Po analysis is prudent, especially where groundwater is reducing with high pH and low 226 Ra concentrations.

Environmental Science & Technology Water↗

Assessment of tsunami hazard to the U.S. Atlantic margin

Tsunami hazard is a very low-probability, but potentially high-risk natural hazard, posing unique challenges to scientists and policy makers trying to mitigate its impacts. These challenges are illustrated in this assessment of tsunami hazard to the U.S. Atlantic margin. Seismic activity along the U.S. Atlantic margin in general is low, and confirmed paleo-tsunami deposits have not yet been found, suggesting a very low rate of hazard. However, the devastating 1929 Grand Banks tsunami along the Atlantic margin of Canada shows that these events continue to occur. Densely populated areas, extensive industrial and port facilities, and the presence of ten nuclear power plants along the coast, make this region highly vulnerable to flooding by tsunamis and therefore even low-probability events need to be evaluated. We can presently draw several tentative conclusions regarding tsunami hazard to the U.S. Atlantic coast. Landslide tsunamis likely constitute the biggest tsunami hazard to the coast. Only a small number of landslides have so far been dated and they are generally older than 10,000 years. The geographical distribution of landslides along the margin is expected to be uneven and to depend on the distribution of seismic activity along the margin and on the geographical distribution of Pleistocene sediment. We do not see evidence that gas hydrate dissociation contributes to the generation of landslides along the U.S. Atlantic margin. Analysis of landslide statistics along the fluvial and glacial portions of the margin indicate that most of the landslides are translational, were probably initiated by seismic acceleration, and failed as aggregate slope failures. How tsunamis are generated from aggregate landslides remains however, unclear. Estimates of the recurrence interval of earthquakes along the continental slope may provide maximum estimates for the recurrence interval of landslide along the margin. Tsunamis caused by atmospheric disturbances and by coastal earthquakes may be more frequent than those generated by landslides, but their amplitudes are probably smaller. Among the possible far-field earthquake sources, only earthquakes located within the Gulf of Cadiz or west of the Tore-Madeira Rise are likely to affect the U.S. coast. It is questionable whether earthquakes on the Puerto Rico Trench are capable of producing a large enough tsunami that will affect the U.S. Atlantic coast. More information is needed to evaluate the seismic potential of the northern Cuba fold-and-thrust belt. The hazard from a volcano flank collapse in the Canary Islands is likely smaller than originally stated, and there is not enough information to evaluate the magnitude and frequency of flank collapse from the Azores Islands. Both deterministic and probabilistic methods to evaluate the tsunami hazard from the margin are available for application to the Atlantic margin, but their implementation requires more information than is currently available.

Marine Geology↗

Simulation of groundwater flow and brine discharge to the Dolores River in the Paradox Valley, Montrose County, Colorado

Salinity, or total dissolved solids (TDS), of the Colorado River affects agricultural, municipal, and industrial water users and is an important concern in the Western United States. In the Paradox Valley of southwestern Colorado, natural discharge of sodium-chloride brine to the Dolores River from the underlying core of a salt-valley anticline accounts for about 6 percent of the salinity load to the Colorado River. Formation of the Paradox Valley began during the Miocene, and subsequent erosion exposed the Pennsylvania Paradox Formation in the core of the anticline where a cap rock, collapse features, breccia, and sodium-chloride saturated brine developed at the top of the exposed salt diapir. The discharge of brine to the Dolores River is affected by these dissolution features, along with seasonal hydrologic conditions and density-dependent flow between older dense brine and the younger fresh groundwater in the overlying alluvial aquifer. To reduce TDS concentrations in the Dolores River through the Paradox Valley, the Bureau of Reclamation has pumped brine from a series of shallow wells adjacent to the river since July 1996. The pumped brine is collected and piped to a deep disposal well where it is injected into the Mississippian Leadville Limestone at a depth of about 4,570-meters below land surface. The pumping and injection operation is collectively known as the Paradox Valley Unit (PVU), and by 2015, the PVU had substantially reduced TDS concentrations in the Dolores River by about 70 percent. Since 2019, injection-pressure limits and related seismic activity have constrained deep-well injection and thus brine pumping at the PVU. In cooperation with the Bureau of Reclamation, the U.S. Geological Survey developed a MODFLOW-6 three-dimensional, variable-density groundwater flow and TDS transport model of the Paradox Valley to evaluate the effects of PVU pumping operations on brine discharge to the Dolores River and to guide additional research. The finite-difference model grid consists of 76 rows and 48 columns oriented from northwest to southeast in alignment with valley topography and groundwater-flow directions in the near-surface freshwater alluvial aquifer. A 7-layer hydrogeologic framework was developed from existing datasets to represent the alluvial aquifer, cap rock, collapse breccia, and groundwater flow and TDS transport from the underlying Paradox Formation salt to the Dolores River. The model represents a 33-year transient calibration period from 1987 through 2020 that includes pre-PVU conditions from 1987 through June 1996 and post-PVU conditions from July 1996 through 2020. A 1,000-year simulation of groundwater flow and coupled TDS transport computed the initial conditions for the subsequent 33-year transient simulation. Observations of precipitation, streamflow, evaporation, agricultural land use, and PVU brine pumping rates were used to specify appropriate boundary conditions to the model representing time-varying recharge, tributary streamflow, groundwater underflow, evapotranspiration (ET), and PVU pumping. Values for average monthly streamflow and TDS concentration at the upstream streamgage, the Dolores River at Bedrock (USGS streamgage 09169500), were specified as model input where the Dolores River enters Paradox Valley. Observed pumping from the PVU, water levels and TDS concentrations in groundwater, and streamflow and estimated TDS concentrations at the downstream streamgage, the Dolores River near Bedrock (USGS streamgage 09171100), were calibration targets that constrained the manual calibration of model parameters representing aquifer hydraulic conductivity, storage, streambed conductance, recharge, and (ET). Two primary model-calibration targets were the match between observed and simulated TDS mass flux from PVU pumping wells and the match between estimated and simulated TDS mass flux to the Dolores River. The simulated TDS mass withdrawn by pumping wells is calculated by the model as the product of the assigned pumping rate and simulated groundwater TDS concentrations. Because actual pumping rates were assigned as simulated values, the total simulated PVU pumping for the 33-year calibration is within 0.5 percent of the observed values. However, simulated concentrations and thus mass flux of TDS withdrawn by the PVU pumping wells were consistently about 26 percent less than observed values for all the simulated time periods (33-year simulation, pre-PVU, and post-PVU). The representation of brine inflow was explored through additional modeling to evaluate the effect of the simulated brine source on groundwater TDS concentrations. Results indicated that a saturated-salt constant-flux brine source best replicated the magnitude and transient pattern observed for TDS mass flux from PVU pumping wells. The simulated TDS mass flux to the Dolores River is compared to estimates based on observed streamflow and specific conductance (SC) data for the downstream streamgage. The calibrated model provided a close fit of simulated to measured streamflow at the downstream streamgage, and the calibrated model fit to estimated TDS concentrations at the downstream streamgage was reasonable. The greatest differences between simulated and estimated values occurred during drought periods from June 2000 to March 2003, May 2012 to June 2013, and October 2013 to October 2014, when simulated TDS concentrations in the river were greater than estimated concentrations. In general, simulated TDS mass flux to the river for the pre-PVU period is in good agreement with estimated values (2-percent difference), but the model overestimated TDS mass flux to the river by about 41 percent during the post-PVU period. The model uncertainty with respect to TDS mass flux to the river indicates other processes or model parameters not well represented by the model are affecting the system, especially during drought. During model calibration, the most sensitive parameters were identified as vertical hydraulic conductivity of the alluvial aquifer, conductance of the Dolores River streambed, ET extinction depth and rate, and recharge rate. Five 5-year scenarios of conditions for 2021–25 were simulated to assist evaluation of alternative strategies to manage the discharge of brine into the Dolores River. The first scenario simulates no PVU pumping and serves as a base case for comparison to the other scenarios. Two scenarios simulate the effects of varying withdrawal timing at an annual rate about one-third less than during 2010 through 2018. During high-flow spring snowmelt runoff periods when brine discharge is naturally minimized, PVU pumping does not substantially affect salinity in the Dolores River, and comparison of these two scenarios indicates that scheduling brine withdrawals during times of low river stage is nearly as effective at reducing TDS mass flux to the river as pumping brine year-round. Cessation of pumping during periods of high river stage may be advantageous for system maintenance, brine injection, and seismic-risk reduction. The fourth scenario tested the effect of reducing irrigation-return flow on brine discharge and predicted a slight reduction of TDS mass flux to the Dolores River, but not as great a reduction as that of using the PVU to remove brine. The fifth scenario simulated 5 years of drought conditions without PVU pumping and indicates brine discharge during drought about 15 percent greater than during average hydrologic conditions. Results from scenario 5 are consistent with the calibrated model results and indicate that aquifer properties and ET processes and parameters may be affecting simulation results during drought. The Paradox Valley groundwater model provides a reasonable overall match to observed conditions in the Dolores River. The model is useful for evaluating relative differences between brine management scenarios to inform PVU operational decisions and to identify gaps in data and process understanding. Representation of the brine source, hydraulic-conductivity parameters, and recharge and ET processes were identified as potential areas for additional field and modeling research. Additional research in the Paradox Valley might include field-data collection that provides additional information on the hydrogeologic framework, groundwater levels, groundwater TDS concentrations, stream characteristics, and aquifer properties. Additional modeling efforts could benefit from applying advanced tools for model development, calibration, and visualization including parameter-estimation and sensitivity analysis. Statistical evaluation of known model uncertainties such as hydraulic conductivity, streambed conductance, representations of the brine source, recharge, and ET could improve the match between simulated and estimated TDS mass flux from PVU pumping wells and to the Dolores River further informing model predictions and system understanding for the Paradox Valley.

Colorado↗

Monitoring a boreal wildfire using multi-temporal Radarsat-1 intensity and coherence images

Twenty-five C-band Radarsat-1 synthetic aperture radar (SAR) images acquired from the summer of 2002 to the summer of 2005 are used to map a 2003 boreal wildfire (B346) in the Yukon Flats National Wildlife Refuge, Alaska under conditions of near-persistent cloud cover. Our analysis is primarily based on the 15 SAR scenes acquired during arctic growing seasons. The Radarsat-1 intensity data are used to map the onset and progression of the fire, and interferometric coherence images are used to qualify burn severity and monitor post-fire recovery. We base our analysis of the fire on three test sites, two from within the fire and one unburned site. The B346 fire increased backscattered intensity values for the two burn study sites by approximately 5–6 dB and substantially reduced coherence from background levels of approximately 0.8 in unburned background forested areas to approximately 0.2 in the burned area. Using ancillary vegetation information from the National Land Cover Database (NLCD) and information on burn severity from Normalized Burn Ratio (NBR) data, we conclude that burn site 2 was more severely burned than burn site 1 and that C-band interferometric coherence data are useful for mapping landscape changes due to fire. Differences in burn severity and topography are determined to be the likely reasons for the observed differences in post-fire intensity and coherence trends between burn sites.

Geomatics, Natural Hazards and Risk↗

Research to inform passage spacing for migratory amphibians and to evaluate efficacy and designs for open elevated road segment (ERS) passages

This is a multifaceted project that includes three main areas of research targeted to inform effective crossing systems for migratory amphibians, a large group of species which are at very high risk from negative impacts from roads within their habitats (Glista et al. 2008, Hamer and McDonnell 2008, Semlitsch 2008, Brehme et al. 2018). The three projects presented in this report are: 1) Movement distances along road barrier fencing and probabilities of reaching a passage: Case study with Yosemite toads in Sierra National Forest, CA. 2) Effectiveness of a novel elevated road segment (ERS) road passage system prototype in providing connectivity for amphibians, reptiles, and small mammals: Case study in Sierra National Forest, CA. 3) Concept designs and transportation engineering evaluation for the ERS on primary roads and highways. This research began in 2018 as part of a larger U.S. Geological Survey (USGS) research program in collaboration with the U.S. Forest Service (USFS), California Department of Transportation (Caltrans), and Western Transportation Institute (WTI; Montana State University) to inform best management practices for barrier and crossing systems for sensitive amphibians and reptiles in California (Langton and Clevenger 2021, Brehme and Fisher 2020). The funding from Department of Transportation (DOT) pooled fund partners (Parks Canada / Government of Canada, Federal Highway Administration (FHWA), U.S. State Departments of Transportation (AK, AZ, CA, CO, IA, MI, MN, NM, NV, OR, WA), Ontario Ministry of Transportation) and managed by the Nevada Department of Transportation (NDOT) supported 2021 field study efforts, analyses of fence movement distances for Yosemite toads, and analysis of the efficacy of a novel ERS passage system to Yosemite toads and other amphibians, reptiles and small mammals. Finally, this pooled fund project includes an assessment by transportation engineers in consultation with USGS and Caltrans to provide insight, guidance, and concept designs for similar crossing solutions that could be implemented on improved roads. This research is meant to inform the distances required between crossings to provide permeability for migratory amphibians (i.e., to allow movements necessary for population persistence across roads) as well as to assess the permeability of a new passage design for amphibians and other small animal species that may provide greater connectivity and offer an alternative to below grade tunnels. The results of these studies add to the current body of knowledge in road ecology and increase the choices of road passage designs for amphibians and other small wildlife species.

California↗

Land use and land cover change in the Greater Yellowstone Ecosystem: 1975-1995

Shifts in the demographic and economic character of the Greater Yellowstone Ecosystem (GYE) are driving patterns of land cover and land use change in the region. Such changes may have important consequences for ecosystem functioning. The objective of this paper is to quantify the trajectories and rates of change in land cover and use across the GYE for the period 1975-1995 using satellite imagery. Spectral and geographic variables were used as inputs to classification tree regression analysis (CART) to find "rules" which defined land use and land cover classes on the landscape. The resulting CART functions were used to map land cover and land use across seven Landsat TM scenes for 1995. We then used a thresholding technique to identify locations that differed in spectral properties between the 1995 and 1985 time periods. These "changed" locations were classified using CART functions derived from spectral and geographic data from 1985. This was similarly done for the year 1975 based on Landsat MSS data. Differences between the 1975, 1985, and 1995 maps were considered change in land cover and use. We calibrated and tested the accuracy of our models using data acquired through manual interpretation of aerial photos. Elevation and vegetative indices derived from the remotely sensed satellite imagery explained the most variance in the land use and land cover classes (-i.e., defined the "rules" most often). Overall accuracies from our study were good, ranging from 94% at the coarsest level of detail to 74% at the finest. The largest changes over the study period were the increases in burned, urban, and mixed conifer-herbaceous classes and decreases in woody deciduous, mixed woody deciduous-herbaceous, and conifer habitats. These changes have important implications for ecological function and biodiversity. The expansion of mixed conifer classes may increase fuel loads and enhance risk to the growing number of rural homes. The reduction of woody deciduous cover types is likely reducing population sizes for the numerous plant and animal species that specialize on this habitat type. Some of these species are also negatively influenced by the increase of rural homes in and near woody deciduous habitats.

Idaho, Montana, Wyoming↗

Compilation and evaluation of data used to identify groundwater sources under the direct influence of surface water in Pennsylvania

A study was conducted to compile and evaluate data used to identify groundwater sources that are under the direct influence of surface water (GUDI) in Pennsylvania. In the early 1990s, the Pennsylvania Department of Environmental Protection (PADEP) implemented the Surface Water Identification Protocol (SWIP) for the identification of GUDI sources. Since the establishment of the SWIP, PADEP has classified more than 500 individual sources across Pennsylvania as GUDI, but Pennsylvania’s complex geology and physiography provide a challenge for a uniform method of GUDI determination. Components used in this study to compile and evaluate data associated with GUDI determination include: (1) a preliminary review of file information for 43 public water-supply wells, (2) quality control and addition of data to PADEP’s database for public water-supply systems to prepare data for analysis, and (3) exploratory evaluation of existing GUDI sources in the database with respect to hydrogeologic and source-construction characteristics that are currently utilized in the assessment methodology. Case files for 43 wells from PADEP’s Northcentral and Southcentral regions were reviewed to: (1) provide a better understanding of how the SWIP was applied in practice, (2) verify and compile missing data, and (3) find additional attributes not previously available that might explain a well’s categorization as GUDI. Review of file information showed that the SWIP outlined in PADEP technical guidance was usually followed, but for some sources, the GUDI determination was more complex and could not be easily summarized. Data compiled for study analyses provided by PADEP include source data derived from public water-supply system case files, a source-information database for public water-supply systems, and Microscopic Particulate Analysis (MPA) results and associated water-quality data for public water-supply system groundwater sources. Data from the Pennsylvania Drinking Water Information System (PADWIS) , which is PADEP’s database for public water-supply systems, were also used for this study. The PADWIS database originally included data for 12,147 groundwater sources (11,812 groundwater sources not under the direct influence of surface water (non-GUDI) wells and 335 GUDI wells). A subset (4,018 wells consisting of 3,842 non-GUDI wells and 175 GUDI wells) of the PADWIS database was created for an analysis and includes only community wells evaluated in accordance with the SWIP. MPA results for 631 community and noncommunity wells were compiled, along with associated water-quality data (alkalinity, chloride, Escherichia coli , fecal coliform, nitrate, pH, sodium, specific conductance, sulfate, total coliform, total dissolved solids, total residue, and turbidity) populated from the PADEP Bureau of Laboratories Sample Information System. Data compiled from sources other than PADEP include spatial data, both naturogenic (for example, average precipitation or distance to closest hydrologic feature) and anthropogenic (for example, percentage of developed or agricultural land cover within a specific vicinity of a public water-supply system well) data representing spatially derived variables. Comparison among wells in the PADWIS dataset subset using the nonparametric Kruskal-Wallis test showed that GUDI wells had significantly older median construction years, shallower depths, and static water levels closer to the land surface than non-GUDI wells and that carbonate aquifers had the highest percentages of wells designated as GUDI (12 percent; 57 wells). Further comparison of wells in the PADWIS database subset using the Spearman’s rho monotonic correlation test illustrated that public water-supply wells designated as GUDI largely occur in unconfined aquifers and have high average yield and shallow static water levels. Assessment of the MPA database subset using the Kruskal-Wallis test showed wells with MPA total risk-factor scores that exceeded zero had older median construction years and shallower casing depths than wells with MPA total risk-factor scores of zero and that carbonate aquifers had the highest percentages of wells with MPA total risk-factor scores exceeding zero (30 percent; 63 wells). Spearman’s rho correlations showed that wells completed in aquifers with depths to major water-bearing zones closer to the land-surface had higher total risk-factor scores resulting from MPA samples. Based on the results of the analyses described in this report, broad conclusions can be drawn regarding site-specific well characteristics as well as anthropogenic and naturogenic factors that could be responsible for a well being designated as GUDI, but the accuracy of these results is dependent on the quality of the data being analyzed. Ultimately, study results serve as an added resource for initial desktop screening of wells to determine if additional site-specific investigation is warranted and underscore the need for field evaluation.

Pennsylvania↗

Quality of pesticide data for groundwater analyzed for the National Water-Quality Assessment Project, 2013–18

The National Water-Quality Assessment (NAWQA) Project of the U.S. Geological Survey (USGS) submitted nearly 1,900 samples collected from groundwater sites across the United States in 2013–18 for analysis of 225 pesticide compounds (pesticides and pesticide degradates, hereafter referred to as “pesticides”) by USGS National Water Quality Laboratory schedule 2437 (S2437). For the associated NAWQA study of pesticide occurrence and concentration in groundwater, and for other studies using pesticide results determined by S2437, it is necessary to assess the ability of reported results to meet data-quality requirements that will allow study objectives to be achieved. This assessment of the quality of S2437 results reported in 2013–18 examined data from field and laboratory quality-control samples, along with third-party performance assessment samples, to estimate bias and variability and to identify their potential sources, with an emphasis on implications for the interpretation of pesticide data for groundwater. Results indicate that measurements produced by the S2437 method for most pesticides have bias and variability that would be considered acceptable for many interpretative studies, which could therefore use the results without qualification or censoring. However, the reported data for a subset of pesticides have the potential for unacceptable contamination bias, high or low recovery bias, or high variability as a consequence of method performance and (or) nonlaboratory factors that could preclude their use for certain common objectives or could necessitate adjustment or qualification to meet those objectives. Based on data for laboratory blanks, censoring of some detections for a subset of pesticides reported by the laboratory in environmental samples might be necessary or desirable to avoid an unacceptably high likelihood of a false-positive result caused by laboratory contamination. The 90-percent upper confidence limit for the 95th percentile of laboratory blank concentration equals or exceeds the minimum reported groundwater concentration in at least 1 water year for 28 pesticides. During at least 1 water year, this upper confidence limit exceeds the maximum laboratory detection limit for 17 pesticides and exceeds the maximum laboratory reporting limit for 3 pesticides (ametryn, atrazine, and diazinon). The level of contamination indicated by this upper confidence limit should not substantially affect the suitability of reported environmental concentrations for any compound for comparison with corresponding human-health benchmarks. Despite being subjected to the same laboratory processes as laboratory blanks, field blanks indicated little evidence of contamination bias. This observation could largely be the consequence of data-reporting practices, which utilize detections in laboratory blanks to censor results in associated field samples (including blanks and environmental samples) when relative concentrations indicate that a result could have a substantial contribution from laboratory contamination. Laboratory censoring appears likely to reduce the risk of false-positive results in environmental samples below the level that laboratory blank results alone would imply. Whereas data available for third-party blind blank samples analyzed in 2018 indicate that only propoxur had any false-positive results, data for pesticides that were not spiked into blind spike samples analyzed in 2013–18 indicate that the false-positive rates for 31 pesticides exceeded 1 percent when considering only detections reported at concentrations greater than the maximum detection limit. Although about half of these pesticides lack substantial supporting evidence of contamination bias based on laboratory blank or field blank detections, indicating that spiking issues or degradation of parent compounds within the spiked samples might be a contributing factor to some false-positive results, these results indicate the need to closely examine detections reported for some pesticides in environmental samples analyzed during a similar period for possible contributions from contamination bias. Data for blind spike samples that were spiked at concentrations above the maximum reporting limit indicate that false-negative rates for eight pesticides exceed 10 percent; substantial low bias could affect results reported for these pesticides in environmental samples analyzed during a similar period. Data for laboratory reagent spikes, which measure recovery of pesticides in blank water, show little evidence for unacceptable recovery bias for S2437 pesticides. However, field matrix spikes, which measure recovery of pesticides in environmental matrices, indicate that degradation and (or) matrix effects could result in moderate to substantial low bias for groundwater results for several pesticides. Low bias could cause some reported concentrations to be categorized as being below a benchmark when the actual concentration in groundwater is greater than the benchmark. Occurrence and concentrations in groundwater could be substantially underrepresented for six pesticides with benchmarks (1H-1,2,4-triazole, asulam, bifenthrin, cis-permethrin, fenbutatin oxide, and naled) that have median recoveries between zero and 50 percent in field matrix spikes. Two compounds (didealkylatrazine and 2-hydroxy-6-ethylamino-4-amino-s-triazine) have median recoveries near or greater than 150 percent in field matrix spikes, indicating a substantial high bias. Plots of data for all spike types show clear changes in the typical recovery with time for some pesticides, which would require further examination for evaluation of temporal trends in environmental concentrations. Data for laboratory reagent spikes indicate that nearly all S2437 pesticides have acceptable variability resulting from random measurement error. Only two compounds (fenbutatin oxide and naled) have F-pseudosigma values greater than 30 percent for recovery, which implies the potential for relatively high variability in reported concentrations and could affect comparison of concentrations to benchmarks and determination of whether concentrations for samples collected at separate locations or times are truly different with a specified level of confidence. Data for third-party blind spike samples show relatively high variability for a greater number of pesticides, although these results likely reflect the influence of degradation and (or) differences in the magnitude and variability of concentrations used for blind spikes relative to laboratory reagent spikes. Detailed analysis of variability using field replicate data is possible for only 12 pesticides on S2437; low variability in analyte detection and concentration is indicated for most of these pesticides in groundwater.

Scientific Investigations Report↗

Nodal seismic deployment on Mauna Loa volcano, Hawaii: Dataset and preliminary insights

Mauna Loa is the largest active volcano on Earth, comprising ∼51% of the Island of Hawai‘i’s landmass and posing significant risks to the island’s communities, infrastructure, and natural environment. Historical eruptions have produced lava flows that have reached the ocean in as little as 3 hr. The timing and location of such lava flows in the past 200 yr underscore how critical determination of the location and geometry of magma storage and structure is for volcanic hazard assessment and eruption forecasting. Now, after nearly 38 yr of volcanic quiescence, Mauna Loa has erupted again. On 27 November 2022, fissures initiated within the summit caldera and then migrated to the northeast rift zone, where they generated a large lava flow that threatened a major highway. To improve our understanding of the geometry of this magma system, we deployed a temporary nodal array on Mauna Loa in the summer of 2024. This increased our seismic coverage sufficiently to image this magmatic system. This nodal array consists of 33 seismometers distributed on and around the volcano and was deployed for over three months to record seismic signals. The primary objective of this project is to resolve the high‐resolution seismic velocity structure and characterize seismic features associated with magma storage and ascent pathways. In this article, we present an overview of the deployment, evaluate the quality of the data, and show example recordings to evaluate the suitability of the data set for future seismic investigations, including earthquake relocation, seismic tomography, and receiver function analysis. Comparisons with nearby permanent broadband and short‐period seismic stations demonstrate that the nodal array recorded high‐quality waveforms, making it a valuable resource for constraining the magmatic system beneath Mauna Loa at multiple scales.

Hawaii↗

Avian influenza shedding patterns in waterfowl: implications for surveillance, environmental transmission, and disease spread

Despite the recognized importance of fecal/oral transmission of low pathogenic avian influenza (LPAI) via contaminated wetlands, little is known about the length, quantity, or route of AI virus shed by wild waterfowl. We used published laboratory challenge studies to evaluate the length and quantity of low pathogenic (LP) and highly pathogenic (HP) virus shed via oral and cloacal routes by AI-infected ducks and geese, and how these factors might influence AI epidemiology and virus detection. We used survival analysis to estimate the duration of infection (from virus inoculation to the last day virus was shed) and nonlinear models to evaluate temporal patterns in virus shedding. We found higher mean virus titer and longer median infectious period for LPAI-infected ducks (10&ndash;11.5 days in oral and cloacal swabs) than HPAI-infected ducks (5 days) and geese (7.5 days). Based on the median bird infectious dose, we found that environmental contamination is two times higher for LPAI- than HPAI-infectious ducks, which implies that susceptible birds may have a higher probability of infection during LPAI than HPAI outbreaks. Less environmental contamination during the course of infection and previously documented shorter environmental persistence for HPAI than LPAI suggest that the environment is a less favorable reservoir for HPAI. The longer infectious period, higher virus titers, and subclinical infections with LPAI viruses favor the spread of these viruses by migratory birds in comparison to HPAI. Given the lack of detection of HPAI viruses through worldwide surveillance, we suggest monitoring for AI should aim at improving our understanding of AI dynamics (in particular, the role of the environment and immunity) using long-term comprehensive live bird, serologic, and environmental sampling at targeted areas. Our findings on LPAI and HPAI shedding patterns over time provide essential information to parameterize environmental transmission and virus spread in predictive epizootiologic models of disease risks.

Journal of Wildlife Diseases↗

Building on crossvalidation for increasing the quality of geostatistical modeling

The random function is a mathematical model commonly used in the assessment of uncertainty associated with a spatially correlated attribute that has been partially sampled. There are multiple algorithms for modeling such random functions, all sharing the requirement of specifying various parameters that have critical influence on the results. The importance of finding ways to compare the methods and setting parameters to obtain results that better model uncertainty has increased as these algorithms have grown in number and complexity. Crossvalidation has been used in spatial statistics, mostly in kriging, for the analysis of mean square errors. An appeal of this approach is its ability to work with the same empirical sample available for running the algorithms. This paper goes beyond checking estimates by formulating a function sensitive to conditional bias. Under ideal conditions, such function turns into a straight line, which can be used as a reference for preparing measures of performance. Applied to kriging, deviations from the ideal line provide sensitivity to the semivariogram lacking in crossvalidation of kriging errors and are more sensitive to conditional bias than analyses of errors. In terms of stochastic simulation, in addition to finding better parameters, the deviations allow comparison of the realizations resulting from the applications of different methods. Examples show improvements of about 30% in the deviations and approximately 10% in the square root of mean square errors between reasonable starting modelling and the solutions according to the new criteria.

Stochastic Environmental Research and Risk Assessm↗

Stochastic model for simulating Souris River Basin precipitation, evapotranspiration, and natural streamflow

The Souris River Basin is a 61,000-square-kilometer basin in the Provinces of Saskatchewan and Manitoba and the State of North Dakota. In May and June of 2011, record-setting rains were seen in the headwater areas of the basin. Emergency spillways of major reservoirs were discharging at full or nearly full capacity, and extensive flooding was seen in numerous downstream communities. To determine the probability of future extreme floods and droughts, the U.S. Geological Survey, in cooperation with the North Dakota State Water Commission, developed a stochastic model for simulating Souris River Basin precipitation, evapotranspiration, and natural (unregulated) streamflow. Simulations from the model can be used in future studies to simulate regulated streamflow, design levees, and other structures; and to complete economic cost/benefit analyses. Long-term climatic variability was analyzed using tree-ring chronologies to hindcast precipitation to the early 1700s and compare recent wet and dry conditions to earlier extreme conditions. The extended precipitation record was consistent with findings from the Devils Lake and Red River of the North Basins (southeast of the Souris River Basin), supporting the idea that regional climatic patterns for many centuries have consisted of alternating wet and dry climate states. A stochastic climate simulation model for precipitation, temperature, and potential evapotranspiration for the Souris River Basin was developed using recorded meteorological data and extended precipitation records provided through tree-ring analysis. A significant climate transition was seen around1970, with 1912–69 representing a dry climate state and 1970–2011 representing a wet climate state. Although there were some distinct subpatterns within the basin, the predominant differences between the two states were higher spring through early fall precipitation and higher spring potential evapotranspiration for the wet compared to the dry state. A water-balance model was developed for simulating monthly natural (unregulated) mean streamflow based on precipitation, temperature, and potential evapotranspiration at select streamflow-gaging stations. The model was calibrated using streamflow data from the U.S. Geological Survey and Environment Canada, along with natural (unregulated) streamflow data from the U.S. Army Corps of Engineers. Correlation coefficients between simulated and natural (unregulated) flows generally were high (greater than 0.8), and the seasonal means and standard deviations of the simulated flows closely matched the means and standard deviations of the natural (unregulated) flows. After calibrating the model for a monthly time step, monthly streamflow for each subbasin was disaggregated into three values per month, or an approximately 10-day time step, and a separate routing model was developed for simulating 10-day streamflow for downstream gages. The stochastic climate simulation model for precipitation, temperature, and potential evapotranspiration was combined with the water-balance model to simulate potential future sequences of 10-day mean streamflow for each of the streamflow-gaging station locations. Flood risk, as determined by equilibrium flow-frequency distributions for the dry (1912–69) and wet (1970–2011) climate states, was considerably higher for the wet state compared to the dry state. Future flood risk will remain high until the wet climate state ends, and for several years after that, because there may be a long lag-time between the return of drier conditions and the onset of a lower soil-moisture storage equilibrium.

Manitoba, North Dakota, Saskatchewan↗

Shifting sands: The influence of coral reefs on shoreline erosion from short-term storm protection to long-term disequilibrium

Climate change is exacerbating shoreline erosion and flooding, posing significant risks to coastal communities. Although traditional coastal defenses such as seawalls, dykes, and breakwaters offer protection from these hazards, their high environmental and economic costs are driving interest in cost-competitive nature-based solutions. Coral reef restoration is a nature-based solution that may be particularly apt to mitigate tropical coastal flooding and shoreline erosion while providing benefits to local tourism, fisheries, and nature. However, the novelty of this field requires studies demonstrating the benefits of reefs for coastal protection. While the flood protection benefits of reefs have been well-documented, their effects on shoreline erosion are comparatively less understood. Here, we investigate the effects of coral reefs on shoreline erosion by comparing tropical beach responses at short and long timescales, as well as identifying important reef structural features influencing coastal erosion rates. Our analyses leveraged two key datasets created in this study: the first derived from a literature review on short-term shoreline erosion due to storm events, and another compiling >80 years of long-term erosion rates, bathymetry, habitat, and wave energy for the Hawaiian Islands of Kauaʻi, Oʻahu, and Maui. Our analyses reveal three key findings regarding the effects of reefs on shoreline erosion. Firstly, we find evidence for the role of reefs in mitigating shoreline erosion during storm events, with coral reef-protected beaches experiencing 97 % less beach volume loss than unprotected beaches. Secondly, a linear regression analysis demonstrates that coral reef structure and wave energy are important predictors of long-term shoreline erosion rates, explaining 34 % of the variation across the Hawaiian Islands. Consistent with prior research, we find beaches protected by coral reefs with shallow reef crests, wide reef flats, calmer offshore conditions, and positioned farther from the shore exhibit lower erosion rates than others. Finally, when comparing historical erosion rates of protected and unprotected beaches in Hawai'i, we find a seemingly incongruous pattern where coral reef-protected beaches eroded up to 2x faster than beaches without reefs. While the cause of the enhanced erosion is yet to be fully understood, a combination of coral reef structural degradation and sea-level rise is likely shifting the equilibrium profiles of reef-protected beaches inshore. These results emphasize the role of coral reefs in reducing coastal erosion during storm events while revealing contrasting erosion patterns over long timescales. Future studies would ideally broaden the scope to include various regions, utilize advanced sediment transport models, and undertake field experiments to deepen our understanding of coral reef-coupled shoreline dynamics.

Nature-Based Solutions↗

The utility of point count surveys to predict wildlife interactions with wind energy facilities: An example focused on golden eagles

Wind energy development is rapidly expanding in North America, often accompanied by requirements to survey potential facility locations for existing wildlife. Within the USA, golden eagles ( Aquila chrysaetos ) are among the most high-profile species of birds that are at risk from wind turbines. To m inimize golden eagle fatalities in areas proposed for wind development, modified point count surveys are usually conducted to estimate use by these birds. However, it is not always clear what drives variation in the relationship between on-site point count data and actual use by eagles of a wind energy project footprint. We used existing GPS-GSM telemetry data, collected at 15 min intervals from 13 golden eagles in 2012 and 2013, to explore the relationship between point count data and eagle use of an entire project footprint. To do this, we overlaid the telemetry data on hypothetical project footprints and simulated a variety of point count sampling strategies for those footprints. We compared the time an eagle was found in the sample plots with the time it was found in the project footprint using a metric we called “error due to sampling”. Error due to sampling for individual eagles appeared to be influenced by interactions between the size of the project footprint (20, 40, 90 or 180 km 2 ) and the sampling type (random, systematic or stratified) and was greatest on 90 km 2 plots. However, use of random sampling resulted in lowest error due to sampling within intermediate sized plots. In addition sampling intensity and sampling frequency both influenced the effectiveness of point count sampling. Although our work focuses on individual eagles (not the eagle populations typically surveyed in the field), our analysis shows both the utility of simulations to identify specific influences on error and also potential improvements to sampling that consider the context-specific manner that point counts are laid out on the landscape.

California↗

Per- and polyfluoroalkyl substances in groundwater from the Great Miami buried-valley aquifer, southwestern Ohio, 2019–20

Groundwater samples were collected during 2019 and 2020 from 23 wells in the Great Miami buried-valley aquifer (GM-BVA) in southwestern Ohio by the U.S. Geological Survey, in cooperation with the Miami Conservancy District, Dayton, Ohio, to determine concentrations of selected per- and polyfluoroalkyl substances (PFAS). The GM-BVA is a glacial outwash and alluvial fill aquifer that is the sole source of water supply for much of the region. Wells had total depths that ranged from 21 to 101 feet below land surface, and groundwater levels that ranged from 1.39 to 52.15 feet below land surface before sampling in 2019. Groundwater and related quality-control samples were sequentially collected from 22 of the 23 wells and analyzed for 24 different PFAS by 2 methods that used proprietary isotope-dilution based adaptations of U.S. Environmental Protection Agency (EPA) method 537.1, termed methods 1 and 2. Method 2 had smaller reporting limits (RL) for 22 of 24 PFAS analyzed and smaller detection limits (DLs) for all 24 PFAS analyzed compared with method 1, which made method 2 more sensitive to detect PFAS. Concentrations of perfluorooctanesulfonate (PFOS) in a groundwater (GW)-method 2 sample from well CL–275 of 1.9 nanograms per liter (ng/L) and perfluorooctanoate (PFOA) in a GW-method 2 sample from well BU–1106 of 2.1 ng/L were greater than their EPA interim health advisory guidances for drinking water (as of June 2022) by about 9,500 and 52,500 percent, respectively. The EPA interim health advisory guidances for PFOS (0.02 ng/L) and PFOA (0.004 ng/L) were also 65 and 215 times less, respectively, than the smallest method 2 DLs for PFOS (1.3 ng/L) and PFOA (0.86 ng/L). Other PFAS were either not detected in GM-BVA groundwater samples or were detected in concentrations less than Ohio action levels or Federal health-risk-based guidance. The most detected PFAS in groundwater was perfluorobutanesulfonate (PFBS), which had concentrations in samples from eight wells that ranged from 1.0 to 8.0 ng/L or from 0.05 to 0.4 percent of its EPA health advisory of 2,000 ng/L for drinking water. The similarity of PFBS (7.8 ng/L), perfluoropentanesulfonate (PFPeS; 8.1 ng/L), and perfluorohexanesulfonate (PFHxS; 14 ng/L) concentrations yielded from the GW-method 1 sample from well CL–275 on July 9, 2019, to those of PFBS (8.0 ng/L), PFPeS (7.8 ng/L), and PFHxS (16 ng/L) from the paired GW-method 2 sample demonstrated the capability of both methods to reproduce PFAS concentrations that were greater than their respective DLs. Non-detection of these PFAS in follow-up GW-method 1 and sequential replicate (Rep–GW-method 1) samples from CL–275 on April 21, 2020, indicated that the 2019 results represented a transient detection in groundwater. Eleven of twenty-three wells sampled in 2019 had from 1 to 4 PFAS detected in one or more groundwater samples or in a paired replicate sample: PFBS in 8 wells and 9 samples; PFHxS in 4 wells and 5 samples; and PFPeS, PFOS, perfluorobutanoate, perfluoropentanoate, PFOA, and perfluorooctanesulfonamide in 1 well and 1 sample each. More PFAS were detected in GW-method 2 samples than GW-method 1 samples because method 2 had smaller RLs and DLs. Results indicate benefits from the analysis of paired samples, sequential replicate samples, and other quality-control samples using analytical methods with sensitive RLs and DLs to verify PFAS concentrations in groundwater. Groundwater-age estimates indicate that water produced from all sampled wells had infiltrated to the water table within the 1947–present (2022) period of PFAS use or environmental presence. Eight wells with detectable PFBS in groundwater from 2019 samples also had groundwater-recharge dates that ranged from 1991 to 2016. Those ages coincided with the possible environmental presence of PFBS as a PFAS byproduct or use as an alternative to PFOS after about 2002. Two wells that had detections of PFHxS in 2019 groundwater samples also had post-2000 groundwater-recharge dates that coincided with the period of use of PFHxS as an alternative to PFOS. Six of nine wells with more than 66-percent of urban land use that was within 0.3 miles of each well, as of 2012, also had 1 to 4 PFAS detected in one of their groundwater samples. Seven of nine wells that produced groundwater in 2019 with an oxic redox category also had one or more PFAS detected in a sample.

Ohio↗

Structured decision-making workshop: Chronic wasting disease management in free-ranging cervids in Massachusetts

This document describes the results of a 2.5-day rapid decision prototype workshop that evaluated management activities for chronic wasting disease (CWD) in Massachusetts (MA) that were either proactive (i.e., actions taken prior to CWD arrival/detection) or reactive (i.e., actions taken after CWD arrival/detection). The workshop was led by members of the Wildlife Section of the MA Division of Fisheries and Wildlife (hereafter referred to as MassWildlife) and included a group of agency communications specialists and district managers. U. S. Geological Survey staff and a volunteer acted as decision facilitators and led the analysis of the decision. Chronic wasting disease is an always fatal neurological disease that has spread across much of North America and threatens the health of deer populations in locations where it occurs (reviewed by Escobar et al. 2020). CWD can spread into new areas via two general mechanisms: (1) natural spread (e.g., dispersal of CWD-infected male white-tailed deer [Odocoileus virginianus]), and (2) anthropogenic spread (e.g., CWD spread facilitated by human intervention; Leiss et al. 2017, Escobar et al. 2020). Once CWD arrives in a state, natural resources agencies spend eight times more on CWD than agencies with no known cases; to cover these new CWD-related management activities, the natural resources agencies are typically forced to reallocate money from existing conservation priorities (Chiavacci, 2022). As of May 2024, there were 34 U.S. states and five Canadian provinces that had detected CWD positive free-ranging and/or captive animals in the family Cervidae (collectively referred to as ‘cervid’ hereafter), and the number of new states/provinces that are detecting CWD for the first time continues to grow (U. S. Geological Survey, May 2024). As of February 2024, the closest CWD positive state to MA with CWD detected in free-ranging white-tailed deer is Pennsylvania. To date, there have been no detections of CWD in MA, but testing has been limited in MA since 2012. The growing number of CWD positive states suggests that there may be increasing risk of CWD entering and establishing in MA as the number of CWD cases increases across North America. According to a 2023 survey of hunters in MA conducted by MassWildlife, 68% of hunters were concerned about CWD entering MA, and 88% of respondents said that it was at least moderately important to keep CWD out of MA; these survey results indicate that most hunters may support CWD risk reduction actions (Martin Feehan, Massachusetts Division of Fisheries and Wildlife, oral communication, 12 Feb 2024). In addition, 23.1% of responding deer hunters in MA have hunted for cervids in CWD-positive states/provinces in the last five years (not including states/provinces that have been able to successfully eradicate CWD following a positive detection). Participants of the survey were also asked, “how many deer have you harvested that tested positive for CWD?”. A total of three respondents said that they had one deer test positive for CWD, which, when extended to the whole population of MA deer hunters, results in an estimated 32 CWD positive deer harvested in CWD-positive states and imported into MA in the last five years. When asked about how they transport harvested deer from out of state into MA, the three participants indicated either “already processed & packaged” or “not applicable.” Note, that in MA, it is a violation of regulation to import whole carcasses or high-risk parts (e.g., head, brain, spinal tissues, bones) of any member of the Cervidae family (wild or captive) from a state/province that has detected CWD; it is legal to bring in deboned meat, cleaned skull caps, hides without the head, or a fixed taxidermy mount (Massachusetts Division of Fisheries and Wildlife, 2024a). To date, testing for CWD has been limited in MA since 2012. However, the data collected from the 2023 MA hunter survey suggests that there is a real risk of CWD being imported by a MA resident who has hunted in a CWD positive state. Therefore, given the higher costs of CWD management post arrival, the potential natural spread of CWD from nearby states, and the risk of CWD introduction via humanmediated cervid movement, MassWildlife is motivated to take actions that minimize the risk of CWD introduction and spread in MA with the ultimate goal of managing thriving wildlife populations and maximizing hunter and general public satisfaction, which are both parts of the MassWildlife mission. A 2.5-day rapid prototyping structured decision making workshop was held with MassWildlife staff to develop a decision framework for CWD management in MA. During the workshop, we defined the context and extent of CWD management activities in MA. Next, we identified four fundamental objectives that help achieve the mission of MassWildlife and that address stakeholder concerns. The fundamental objectives included: (1) maximizing hunter satisfaction and participation, (2) maximizing public satisfaction (non-consumptive), (3) maximizing health and sustainability of cervids, and (4) maximizing the efficiency of CWD management. Then, we generated a list of five alternatives (i.e., strategies) that varied the intensity of proactive and reactive actions. The five strategies were: (1) minimal proactive and minimal reactive actions, (2) intermediate proactive and intermediate reactive actions, (3) intensive proactive and intermediate reactive actions, (4) minimal proactive and intensive reactive actions, and (5) intensive proactive and intensive reactive actions. Lastly, we estimated the performance of each strategy on the fundamental objectives and assessed the overall performance of strategies relative to one another. We did so by first estimating the consequences of each alternative strategy on fundamental objectives using expert elicitation, and then, we elicited objective weights from MassWildlife staff to incorporate the relative importance of different fundamental objectives. Given that it is unknown when CWD will arrive in MA, we evaluated the performance of alternative strategies against fundamental objectives given three distinct scenarios for time to arrival of CWD: introduction in 2.5, 7.5, or 10+ years. The preliminary results of the rapid prototype indicate that the performance of the CWD management strategies that we evaluated depends on when CWD first arrives in MA. If CWD were to arrive in 2.5 or 7.5 years from now (February, 2024), then the ‘minimal proactive and minimal reactive’ strategy performs the best on both the deer population and cost fundamental objectives (fundamental objectives 3 & 4), but the ‘intensive proactive and intensive reactive’ strategy performs best on both of the human dimensions fundamental objectives (fundamental objectives 1 & 2) as well as the minimize CWD prevalence objective (also related to fundamental objective 3). We also found that public trust is likely to remain high across all five alternative strategies if CWD arrives after year 10, but public trust decreases if CWD arrives in year 2.5 or 7.5. After incorporating objective weights, we found that in scenarios where CWD arrives in the near-term (in years 2.5 or 7.5), an intermediate strategy (e.g., ‘intermediate proactive and intermediate reactive’ or ‘intensive proactive and intermediate reactive’) performed best, and the ‘minimal proactive and intensive reactive’ strategy performed worst. Conversely, if CWD were to arrive after 10 years, then the ‘minimal proactive and minimal reactive’ and ‘minimal proactive and intensive reactive’ strategies performed best. Collectively, these results suggest that the decision on which alternative strategy to employ is sensitive to when CWD arrives in MA. Following the discussion of the preliminary results, we identified the following four next steps. First, we discussed how a more detailed communications plan is needed and would likely alter the performance estimates of the alternative strategies on fundamental objectives 1 & 2, which were hunter and public satisfaction, respectively. The development of the communication plan would likely be easier once the alternative actions have been identified along with the audience and message. Second, a surveillance plan could be a useful tool to inform CWD management. Surveillance for CWD was performed in MA annually from 2002 to 2012 (n = 4,356 wild white-tailed deer and moose [Alces alces] samples). Limited surveillance was conducted from 2013 to 2022; and in 2023, 242 wild samples were collected. It is not clear whether MA needs a robust or minimal surveillance plan (e.g., is a minimal surveillance plan enough to detect the pathogen at the threshold that would trigger action?), or what type of invasion event the surveillance plan should target (e.g., natural vs anthropogenic spread events). The use of decision trees and a formal risk assessment may help answer these questions. Third, some of the elicited estimates from experts during this rapid prototype could be replaced with empirical data. Lastly, given that the decision was sensitive to when CWD arrived in MA and a surveillance plan would rely on the mode of introduction, forecasting and predicting the CWD invasion front and/or the likelihood of different incursion events across MA would provide valuable insights.

Massachusetts↗