USGS ScienceSearch

SEARCH · USGS Science

Results for “Natural Resources Research”

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.

614 records · Page 3Linked to original sources

Projected sea-level rise and high tide flooding at Dry Tortugas National Park, Florida

Introduction National parks and preserves in the South Atlantic-Gulf Region contain valuable coastal habitats such as tidal wetlands and mangrove forests, as well as irreplaceable historic buildings and archeological sites located in low-lying areas. These natural and cultural resources are vulnerable to accelerated sea-level rise and escalating high tide flooding events. Through a Natural Resources Preservation Program-funded project during 2021–23, the U.S. Geological Survey, in collaboration with the National Park Service, estimated the probability of inundation at Dry Tortugas National Park, Florida, and several other parks under various sea-level rise scenarios and contemporary high tide flooding thresholds. The maps produced for this effort can be used to assess potential habitat change and explore how infrastructure and cultural resources within the park may be exposed to future flooding-related hazards.

Florida

A probabilistic assessment methodology for the evaluation of geologic energy storage capacity—Natural gas storage in depleted hydrocarbon reservoirs

The need for energy storage, particularly underground, where capacity and duration may far exceed battery storage technologies, is especially relevant given the increasing demands for reliable power alongside the development of intermittent renewable electricity sources. Geologic energy storage facilities already exist, and expanded use would enable storing gases such as methane and hydrogen. In 2018, a National Academies of Sciences, Engineering, and Medicine report, “Future Directions for the U.S. Geological Survey's Energy Resources Program,” recommended that the U.S. Geological Survey (USGS) prioritize assessing underground energy storage in geologic formations in the United States. The U.S. Geological Survey has since developed a methodology for assessing natural gas storage capacities in depleted hydrocarbon reservoirs on a national scale. The methodology introduced in this report prescribes three approaches for calculating gas storage capacity. This methodology relies on the availability of input data, including cumulative hydrocarbon production records, reservoir petrophysical properties, and reservoir pressure data. Assessment inputs can be obtained from public, State-level databases and propriety national-scale databases, although the use of analogs could be warranted for estimating input parameters. Probabilistic assessment results are aggregated to play, petroleum province, regional, and national scales. The steps defined in this report are demonstrated on the Michigan Basin Province, which includes the Mississippian Sandstone Gas Play and the Clinton Structural Play. This methodology could be used to systematically and consistently assess hydrocarbon plays and provinces for natural gas storage capacity across the United States.

Michigan

Storm impact scale for barrier islands

A new scale is proposed that categorizes impacts to natural barrier islands resulting from tropical and extra-tropical storms. The proposed scale is fundamentally different than existing storm-related scales in that the coupling between forcing processes and the geometry of the coast is explicitly included. Four regimes, representing different levels of impact, are defined. Within each regime, patterns and relative magnitudes of net erosion and accretion are argued to be unique. The borders between regimes represent thresholds defining where processes and magnitudes of impacts change dramatically. Impact level 1 is the 'swash' regime describing a storm where runup is confined to the foreshore. The foreshore typically erodes during the storm and recovers following the storm; hence, there is no net change. Impact level 2 is the 'collision' regime describing a storm where the wave runup exceeds the threshold of the base of the foredune ridge. Swash impacts the dune forcing net erosion. Impact level 3 is the 'overwash' regime describing a storm where wave runup overtops the berm or, if present, the foredune ridge. The associated net landward sand transport contributes to net migration of the barrier landward. Impact level 4 is the 'inundation' regime describing a storm where the storm surge is sufficient to completely and continuously submerge the barrier island. Sand undergoes net landward transport over the barrier island; limited evidence suggests the quantities and distance of transport are much greater than what occurs during the 'overwash' regime.

Journal of Coastal Research

Projected sea-level rise and high tide flooding at Biscayne National Park, Florida

Introduction National parks and preserves in the South Atlantic-Gulf Region contain valuable coastal habitats such as tidal wetlands and mangrove forests, as well as irreplaceable historic buildings and archeological sites located in low-lying areas. These natural and cultural resources are vulnerable to accelerated sea-level rise and escalating high tide flooding events. Through a Natural Resources Preservation Program-funded project during 2021–23, the U.S. Geological Survey, in collaboration with the National Park Service, estimated the probability of inundation at Biscayne National Park, Florida, and several other parks under various sea-level rise scenarios and contemporary high tide flooding thresholds. The maps produced for this effort can be used to assess potential habitat change and explore how infrastructure and cultural resources within the park may be exposed to future flooding-related hazards.

Florida

Projected sea-level rise and high tide flooding at San Juan National Historic Site, Puerto Rico

Introduction National parks and preserves in the South Atlantic-Gulf Region contain valuable coastal habitats such as tidal wetlands and mangrove forests, as well as irreplaceable historic buildings and archeological sites located in low-lying areas. These natural and cultural resources are vulnerable to accelerated sea-level rise and escalating high tide flooding events. Through a Natural Resources Preservation Program-funded project during 2021–23, the U.S. Geological Survey, in collaboration with the National Park Service, estimated the probability of inundation at San Juan National Historic Site, Puerto Rico, and several other parks under various sea-level rise scenarios and contemporary high tide flooding thresholds. The maps produced for this effort can be used to assess potential habitat change and explore how infrastructure and cultural resources within the park may be exposed to future flooding-related hazards.

San Juan National Historic Site

Scoping decision-maker needs and science availability to support regional natural capital accounting in the U.S. Colorado River Basin

Natural capital accounting has the potential to yield important policy insights at multiple scales, but there remains a disconnect between regional-scale natural capital accounts and their use for informing policy. In this paper, we propose a roadmap that could lead to the creation of policy-relevant regional accounts, with steps split across an initial scoping phase and a subsequent development phase. We demonstrate the scoping steps in action with an application to the Colorado River Basin (“Basin”), a large watershed in the southwestern United States (U.S.) that has faced aridification and substantial high-profile tradeoffs around the use of its water and other natural resources. Drawing on prior U.S. Geological Survey science co-production efforts, we conducted a series of eight discussion sessions with 41 scientists and science representatives whose work is relevant to Basin water, riparian and riverine ecosystems, upland ecosystems and energy and minerals. We summarise participants' thoughts on key topics and economic linkages, their insights and questions of interest and their recommendations on existing scientific data sources and gaps. We evaluate the suitability of the available data for construction of System of Environmental-Economic Accounting (SEEA) Central Framework and SEEA Ecosystem Accounting accounts, including those for land, water, forests, energy and minerals and ecosystems (covering extent, condition and ecosystem services). We present a series of lessons learned during the scoping phase, as well as lessons that could be relevant for future practitioners engaging in the development phase. The information can help guide the development of timely and relevant regional-scale environmental-economic accounts in the U.S. and beyond.

Arizona, California, Colorado, Nevada, New Mexico,

Projected sea-level rise and high tide flooding at Big Cypress National Preserve, Florida

Introduction National parks and preserves in the South Atlantic-Gulf Region contain valuable coastal habitats such as tidal wetlands and mangrove forests, as well as irreplaceable historic buildings and archeological sites located in low-lying areas. These natural and cultural resources are vulnerable to accelerated sea-level rise and escalating high tide flooding events. Through a Natural Resources Preservation Program-funded project during 2021–23, the U.S. Geological Survey, in collaboration with the National Park Service, estimated the probability of inundation at Big Cypress National Preserve, Florida, and several other parks under various sea-level rise scenarios and contemporary high tide flooding thresholds. The maps produced for this effort can be used to assess potential habitat change and explore how infrastructure and cultural resources within the park may be exposed to future flooding-related hazards.

Florida

Decoding the oxidative digestion mechanism for polystyrene nanoplastic detection in the Great Lakes using a customizable Raman spectral processing algorithm

Despite the concerns surging around nanoplastics (NPs) regarding their prevalence and bioavailability in freshwater systems, robust detection of NPs in complex environmental matrices is hindered by the lack of standardized sample pretreatment and a mechanistic understanding of oxidative digestion. Here, we systematically investigate the interaction between hydrogen peroxide (H 2 O 2 ) and polystyrene (PS) NPs during digestion in deionized (DI) water and four environmental matrices from in and around the Great Lakes Basin. To facilitate high-throughput analysis, we develop Pre_peak, a customizable Raman spectral processing algorithm that achieves >99% accuracy for both NP identification and interference rejection, allowing reliable NP quantification via pixel counting and systematic decoding of the oxidative digestion mechanisms. In DI water, varying H 2 O 2 doses from 0 to 30% has negligible effects on the recovery and Raman signal intensity of PS NPs over 24 hours of digestion. However, morphological changes and aggregation of PS NPs are observed when the H 2 O 2 dose exceeds 20%. Prolonged digestion further leads to progressive NP loss. In natural waters, the optimal dosage and digestion duration depend on matrix characteristics, including dissolved organic matter (DOM) and ion composition. This study provides mechanistic insights into NP–oxidant interactions and underscores the need for matrix-tailored digestion protocols to advance standardized NP detection in freshwater environments.

Great Lakes

Peak streamflow trends in Montana and northern Wyoming and their relation to changes in climate, water years 1921–2020

Frequency analysis on annual peak streamflow (hereinafter, peak flow) is essential to water-resources management applications, including critical structure design (for example, bridges and culverts) and floodplain mapping. Nonstationarity is a statistical property of a peak-flow series such that the distributional properties (the mean, variance, or skew) change either gradually (monotonic trend) or abruptly (shift, step change or change point) through time. Not incorporating or accounting for observed nonstationarity into peak-flow frequency analysis might result in a poor representation of the true probability of large floods and thus misrepresent the actual flood risks to life and property. This report summarizes how hydroclimatic variability might affect the temporal and spatial distributions of peak-flow data in the State of Montana (and northern Wyoming) and is part of a larger study to document peak-flow nonstationarity and hydroclimatic changes across a nine-State region consisting of Illinois, Iowa, Michigan, Minnesota, Missouri, Montana, North Dakota, South Dakota, and Wisconsin. A wide range of analyses and statistical approaches are applied to document the primary mechanisms controlling floods and characterize temporal changes in hydroclimatic variables and peak flows. This study was completed in cooperation with the Montana Department of Natural Resources and Conservation. The purpose of this report is to characterize temporal and spatial patterns of nonstationarity in peak flows and hydroclimatology in Montana and northern Wyoming. In this evaluation, peak-flow, daily streamflow, and model-simulated gridded climatic data were examined for monotonic trends, change points, and other statistical properties that might indicate changing climatic and environmental conditions. This report includes background information on the study area, the history of U.S. Geological Survey peak-flow data collection and frequency analysis in Montana, and the review of research relating to hydroclimatic variability and change in Montana. This study might help provide a framework for addressing potential nonstationarity issues in peak-flow frequency updates that commonly are completed by the U.S. Geological Survey in cooperation with other agencies throughout the Nation. The analytical structure of this study includes analyses of monotonic trends and change points in numerous hydroclimatic variables in assigned 30-, 50-, 75-, and 100-year analysis periods. For Montana and part of Wyoming, the 30-, 50-, 75, and 100-year analyses included 157, 70, 48, and 12 streamgages, respectively. For those streamgages, nonstationarities were analyzed in the following variables: (1) climatic variables, including annual and seasonal (winter, spring, summer, and fall) temperature and precipitation; (2) daily streamflow variables, including the annual center of volume duration, annual center of volume median, and peaks over threshold with a mean of four events per year; and (3) annual peak-flow variables, including peak-flow timing and magnitude. A likelihood approach was used to express statistical confidence and assign the nonstationarity results as likely upward or downward (highest statistical confidence), somewhat likely upward or downward (less statistical confidence), or about as likely as not (little statistical confidence; hereinafter, neutral). For the nonstationarity analyses of the climatic, daily streamflow, and peak-flow variables, the results are presented in detail and discussed with respect to statewide patterns and geographic variability. For each of the 30-, 50-, and 75-year analyses, peak-flow change-point and monotonic trend analyses were compiled for streamgages classified with likely downward or likely upward trends. For those streamgages, the associated basin characteristics and nonstationarity results for peak-flow timing, daily streamflow, and climatic variables were investigated and statistically compared to discern associations among other variables that might contribute to the peak-flow nonstationarity results. The 50- and 75-year peak-flow nonstationarities identified in this study are mostly downward, in association with mostly upward temperature and potential evapotranspiration:precipitation monotonic trends. For the 50-, 75-, and 100-year analyses, the peak-flow change points are predominantly downward and are concentrated in the 1970s and 1980s, which indicates general consistency among the longer trend periods. These findings are in association with substantial research documenting globally rising temperature and atmospheric greenhouse gas concentrations that might be largely attributed to anthropogenic activities. Anthropogenic effects might represent long-term (on the order of several decades to more than a century) climate changes that might happen within highly variable natural climate fluctuations. Several paleo studies in the north-central United States have indicated that hydroclimatic extremes (that is, low- and high-streamflow conditions) before European settlement have been outside of extremes since the 1900s. Depending on the interactions of anthropogenic effects and natural climate variability, extreme high-streamflow conditions might occur in the future, even in the presence of long-term downward peak-flow trends.

Montana, Wyoming

Global performance of remote sensing-based and reanalysis-driven models to estimate open water evaporation

Evaporation plays an essential role in the water cycle, influencing local and regional climates while directly impacting water availability in lakes. However, directly measuring evaporation over water bodies remains challenging due to the high costs of installing and maintaining the required in situ instrumentation. Although several remote sensing algorithms have been providing evaporation estimates, the lack of a global validation hinders our understanding of their relative uncertainties and performances across different regions. Here, we analyze the performance of a suite of models that leverage satellite data and meteorological reanalysis to estimate evaporation over lakes worldwide. We compare 3 remote sensing-based models, 1 reanalysis-driven model and 1 ensemble approach, using in situ observations from 27 lakes representing a diverse range of geographic and climatic regions. Our results demonstrate that, overall, the ensemble outperformed any individual model in terms of accuracy, with a RMSE and a bias of 1.3 and 0.3 mm day −1 , respectively. These findings highlight the benefits of using an ensemble approach to estimate open water evaporation with satellite-based models at the global scale, leveraging the unique strengths of each model. For the individual models, differences in the representation of heat storage changes and advection effects led to lower values of RMSE and bias, depending on the location and depth of the lakes. This study sets the path for future improvement of open water evaporation algorithms globally, while remote sensing techniques are proven satisfactory to monitoring of water loss in lakes globally, an essential step toward effective large-scale water resources management.

Water Resources Research

Modeling legacy nitrogen transport under instantaneous, steady-state, and transient groundwater flow conditions

In hydrologic settings where groundwater discharge contributes substantially to surface waters, legacy nitrogen in groundwater can confound surface water nitrogen loads estimated exclusively from current terrestrial sources. Additionally, legacy nitrogen in groundwater can contribute to lagged responses to nitrogen management efforts. Some methods of estimating groundwater contributions to surface water nitrogen loads account for legacy nitrogen, while others do not. The resulting differences are rarely quantified. We used a numerical modeling framework to compare three methods of estimating time-varying annual groundwater nitrogen loads to surface water receptors on eastern Long Island, New York. The instantaneous load method used steady-state contributing areas and includes no temporal groundwater lag. The second method used numerical simulations of nitrogen loads under steady-state flow, which captures groundwater transport lags but omits the annual variability in transient hydrologic stresses. The third method numerically simulated both transient groundwater flow and nitrogen transport to explicitly capture the effects of legacy nitrogen in groundwater. Depending on antecedent nitrogen and hydrologic conditions, historical nitrogen loads estimated from the numerical simulations were sometimes similar (<10% difference) and other times substantially different (±100%) from the instantaneous load estimates. Additionally, simulated future surface water nitrogen loads responded asymptotically over several decades following reductions in terrestrial nitrogen sources, further highlighting the effect of groundwater transport lag times. The comparison of the three methods, quantification of historical interannual variability, and prediction of lagged responses to nitrogen source reductions provide important context for decision makers using estimated groundwater nitrogen loads to help evaluate nitrogen management efficacy.

New York

The saline groundwater legacy of a large buried coastal paleo-estuary

Elevated groundwater salinity in coastal regions threatens the beneficial use of fresh groundwater. Coastal groundwater management typically focuses on preventing intrusion from modern sources of seawater; however, past geological processes can also leave a legacy of saline groundwater now hidden in the subsurface. Here, multiple extensive airborne electromagnetic surveys provide detailed evidence of residual salinity from a paleo-estuary filling a late Pleistocene incised valley impacting more than 10,000 km 2 that is now hidden beneath coastal Louisiana’s deltaic plain. Our results show that the three-dimensional pattern of saline groundwater beneath Louisiana mimics that of near-surface aquifers surrounding the modern Delaware Bay estuary, fingerprinting the signature of the past drowning of a large, incised valley of the Mississippi River following post-glacial sea-level rise. These findings demonstrate a new framework for understanding legacy sources of saltwater critical for managing stressed water resources along global coastlines.

Louisiana

Modeling seawater intrusion along the Alabama coastline using physical and machine learning models to evaluate the effects of multiscale natural and anthropogenic stresses

Seawater intrusion threatens groundwater resources in coastal regions, including southern Baldwin County, Alabama, where the freshwater-saltwater interface dynamics remain poorly understood. To address this gap, this study uses combined physics-based and machine-learning models to quantify seawater intrusion caused by natural (storm surges) and anthropogenic (human activities) perturbations. The long short-term memory network and wavelet analysis were used to assess vertical aquifer vulnerabilities, revealing that the shallow part of the Coastal lowlands aquifer system (CL1) in the southern Baldwin County region is more susceptible to sea level rise and groundwater extraction than deeper aquifers. Based on these findings, a cross-sectional numerical model (physics approach) for the CL1 aquifer was developed to evaluate tidal and storm surge effects, using Tropical Storm Claudette (June 2021) as a case study. Results showed that tidal fluctuations had a minimal impact on the saltwater-freshwater interface location, whereas storm surges caused substantial inland movement, with effects lasting for nine months. The steady-state version of the three-dimensional (3D) physical model predicted seawater intrusion across the entire area, and convolutional neural network-based modeling further validated the model results. The 3D physical model was also applied to a smaller area to assess human impact on the saltwater interface due to two groundwater pumping scenarios (± 50% of the baseline pumping rate). Results revealed that a 50% increase in groundwater withdrawals caused seawater to advance ~ 320 m inland, whereas a 50% reduction led to a ~ 270-meter retreat. This study highlights the vulnerability of Alabama’s shallow coastal aquifers to seawater intrusion due to storm surges and human activities, and demonstrates that combining physics-based models with machine learning approaches can improve groundwater predictions, though its accuracy depends on the availability of site-specific data.

Alabama

Toward a new framework to evaluate process-based model configurations and quantify data worth prior to calibration

Model criticism, discrimination, and selection methods often rely on calibrated model outputs. Because calibration can be computationally expensive, model criticism can first be undertaken by assessing model outputs obtained from limited prior parameter ensembles. However, such prior-based methods are often heuristic and do not formalize the notion of balancing model consistency with data and model complexity (i.e., model adequacy). We present a new framework to discriminate among candidate models prior to calibration that formalizes prior-to-calibration model adequacy into a metric to implicitly balance prior model output data coverage with model complexity represented by prior output (co)variance. The prior model adequacy metric “Mahalanobis distance deviation” quantifies the deviation of (a) the set of squared Mahalanobis distances of data from a prior model output distribution from (b) the set of squared Mahalanobis distances of data from their own distribution. A new data worth metric “discernment value” is also presented which quantifies the value of data for screening less-adequate models prior to calibration. Discernment value is calculated from the change in variance of a weighted average of prior model outputs from all candidate models due to less-adequate model outputs receiving lower weight. The framework is demonstrated using a one-dimensional groundwater flow model with eight possible configurations. A synthetic data network is used to test the framework. Results show the framework identifies the candidate models most similar to the true model used to create the synthetic data. Discernment values show variation in the value of different data types and locations for screening less-adequate models.

Water Resources Research

Invited perspectives: Integrating hydrologic information into the next generation of landslide early warning systems

Although rainfall-triggered landslides are initiated by subsurface hydro-mechanical processes related to the loading, weakening, and eventual failure of slope materials, most landslide early warning systems (LEWS) have relied solely on rainfall event information. In previous decades, several studies demonstrated the value of integrating proxies for subsurface hydrologic information to improve rainfall-based forecasting of shallow landslides. More recently, broader access to commercial sensors and telemetry for real-time data transmission has invigorated new research into hydrometeorological thresholds for LEWS. Given the increasing number of studies across the globe using hydrologic monitoring, mathematical modeling, or both in combination, it is now possible to make some insights into the advantages versus limitations of this approach. The extensive progress demonstrates the value of in situ hydrologic information for reducing both failed and false alarms, through the ability to characterize infiltration during, as well as the drainage and drying processes between major storm events. There are also some areas for caution surrounding the long-term sustainability of subsurface monitoring in landslide-prone terrain, as well as unresolved questions in hillslope hydrologic modeling, which relies heavily on the assumptions of diffuse flow and vertical infiltration but often ignores preferential flow and lateral drainage. Here, we share a collective perspective based on our previous collaborative work across Europe, North America, Africa, and Asia to discuss these challenges and provide some guidelines for integrating knowledge of hydrology and climate into the next generation of LEWS. We propose that the greatest opportunity for improvement is through a measure-and-model approach to develop an understanding of landslide hydro-climatology that accounts for local controls on subsurface storage dynamics. Additionally, new efforts focused on the subsurface hydrology are complementary to existing rainfall-based methods, so leveraging these with near-term precipitation forecasts is a priority for increasing lead times.

Alaska

Quantifying groundwater response and uncertainty in beaver-influenced mountainous floodplains using machine learning-based model calibration

Beavers ( Castor canadensis ) alter river corridor hydrology by creating ponds and inundating floodplains, and thereby improving surface water storage. However, the impact of inundation on groundwater, particularly in mountainous alluvial floodplains with permeable gravel/cobble layers overlain by a soil layer, remains uncertain. Numerical modeling across various floodplain structures considers topographic and sediment complexity and multidirectional flow, linking inundation to groundwater response. This study develops a model-data integration workflow to address uncertainty in groundwater response to beaver-induced inundations in a mountainous alluvial floodplain in the Upper Colorado River Basin. Uncertain factors include seasonal hydrologic dynamics, hydraulic conductivities, floodplain structures, and meteorological forcings. We employed an ensemble of groundwater models, based on geophysical and hydrologic data, with machine learning-based calibration using a neural density estimator. This allowed us to quantify the vertical flux from the soil layer to the permeable gravel bed, the down-valley underflow within the gravel bed, and their ratios. Results show a significant increase in the vertical flux relative to down-valley underflow, from 2% during dry pond periods to 20% during wet periods, serving as an analogy for conditions without and with beaver ponds. The study highlights the influence of floodplain structure on groundwater storage, water balance, and water quality impacted by beaver ponds. A thick gravel bed layer, with a large down-valley underflow, minimizes the effect of beaver-induced inundation on water quality. We emphasize the need for field-scale measurements of floodplain structure and improved characterization of evapotranspiration changes to reduce uncertainty in groundwater response.

Colorado

Modeling future groundwater depletion to evaluate sustainability goals set under the Sustainable Groundwater Management Act in the critically overdrafted basins of the Central Valley, California, USA (2020–2070)

In 2014, California's Sustainable Groundwater Management Act (SGMA) mandated local agencies to devise and implement groundwater sustainability plans to address critically overdrafted conditions throughout the state's aquifers. However, the feasibility of these agencies' sustainability goals has not previously been assessed through a regional-scale, integrative lens. Here, we develop and analyze a novel, basin-wide database of 936 sustainability indicator wells located within Central Valley subbasins designated as critically overdrafted, most of which lie in the San Joaquin Valley. Our database shows 2040 groundwater elevation goals vary widely from 60 m above to 80 m below 2020 levels, with variability within and between adjacent subbasins. To evaluate the feasibility of achieving these goals, we coupled the database with a regional hydrologic model (Central Valley Hydrologic Model version 2) and simulated multiple future pumping scenarios. Results show that under increased groundwater demand, 60%–70% of indicator wells may fail to meet their 2040 goals. Even a 50% reduction from 2020 demand levels leaves nearly 40% of wells failing to meet their sustainability thresholds by 2040. Baseline models show that by 2070, up to 70% of wells could fail to meet their goals due to large-scale, spatially connected regions of groundwater depletion. This integrated framework, linking the first region-wide compilation of SGMA indicator wells with a regional groundwater model, demonstrates that many local sustainability goals may be unattainable with substantial (up to 50%) reductions in pumping. Additional management interventions, such as expanded recharge or coordinated demand reductions, may help achieve sustainability goals.

California

Effects of a low-lipid diet on the gut microbiome and head kidney transcriptome of juvenile Chinook Salmon

Objective Pacific salmon Oncorhynchus spp. reared in production hatcheries are typically fed high-lipid, energy-dense diets to achieve large size and high body condition prior to release. In contrast, juveniles in natural environments tend to consume low-lipid, high-protein diets, and fish reared for research or conservation purposes are sometimes fed diets that are formulated to mimic natural diets and promote wild-like phenotypes. Understanding how these alternative diets affect fish health beyond growth and body condition could ultimately contribute to improving hatchery fish fitness. Methods In this work, we evaluated changes in the fecal microbiome and gene expression of juvenile Chinook Salmon O. tshawytscha on a standard high-lipid hatchery diet versus a low-lipid diet formulated to mimic the nutrition profile of natural-origin fish. To evaluate the time scale at which diet alters the fecal microbiome, we collected longitudinal samples over a 12-week period and switched the diets of a subset of fish twice during the experiment. We used 16S ribosomal RNA gene amplicon sequencing to characterize fecal microbiome differences between fish on the two diets as well as hatchery-reared fish at a production hatchery, hatchery fish that had been captured after release into a stream, and natural-origin, stream-reared fish of similar ages. Additionally, we conducted RNA sequencing on head kidney samples from laboratory-reared fish to evaluate changes in gene expression in this important immune organ. Results We found that the low-lipid diet and the hatchery diet resulted in microbiomes that differed from the microbiome of natural-origin fish and from each other and that diet-driven changes to the microbiome could occur in under 14 d. The low-lipid diet did not result in a microbiome that resembled the microbiome of naturally produced fish. Instead, the low-lipid diet resulted in a microbiome community that was distinct from those of fish reared on the hatchery diet and fish sampled from the wild. The RNA sequencing results indicated differential enrichment of pathways related to immunity, metabolism, and hormone synthesis between fish that were fed the two experimental diets. Conclusions The results suggest that additional environmental factors influence the microbiome more strongly than diet formulation or that the low-lipid diet has a smaller effect on the microbiome than a natural, ­invertebrate-based diet. Given that the gut microbiome and systemic immune function contribute significantly to disease resistance, our findings highlight the importance of understanding how diets fed to fish in captivity may affect fish health beyond growth and body condition metrics.

Journal of Aquatic Animal Health