Search NASA⌕ Search

SEARCH · Search NASA

Results for “flow pattern”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 397 records · Page 22

Holocene thinning of Darwin and Hatherton glaciers, Antarctica, and implications for grounding-line retreat in the Ross Sea

Chronologies of glacier deposits in the Transantarctic Mountains provide important constraints on grounding-line retreat during the last deglaciation in the Ross Sea. However, between Beardmore Glacier and Ross Island – a distance of some 600 km – the existing chronologies are generally sparse and far from the modern grounding line, leaving the past dynamics of this vast region largely unconstrained. We present exposure ages of glacial deposits at three locations alongside the Darwin–Hatherton Glacier System – including within 10 km of the modern grounding line – that record several hundred meters of Late Pleistocene to Early Holocene thickening relative to present. As the ice sheet grounding line in the Ross Sea retreated, Hatherton Glacier thinned steadily from about 9 until about 3 ka. Our data are equivocal about the maximum thickness and Mid-Holocene to Early Holocene history at the mouth of Darwin Glacier, allowing for two conflicting deglaciation scenarios: (1) ~500 m of thinning from 9 to 3 ka, similar to Hatherton Glacier, or (2) ~950 m of thinning, with a rapid pulse of ~600 m thinning at around 5 ka. We test these two scenarios using a 1.5-dimensional flowband model, forced by ice thickness changes at the mouth of Darwin Glacier and evaluated by fit to the chronology of deposits at Hatherton Glacier. The constraints from Hatherton Glacier are consistent with the interpretation that the mouth of Darwin Glacier thinned steadily by ~500 m from 9 to 3 ka. Rapid pulses of thinning at the mouth of Darwin Glacier are ruled out by the data at Hatherton Glacier. This contrasts with some of the available records from the mouths of other outlet glaciers in the Transantarctic Mountains, many of which thinned by hundreds of meters over roughly a 1000-year period in the Early Holocene. The deglaciation histories of Darwin and Hatherton glaciers are best matched by a steady decrease in catchment area through the Holocene, suggesting that Byrd and/or Mulock glaciers may have captured roughly half of the catchment area of Darwin and Hatherton glaciers during the last deglaciation. An ensemble of three-dimensional ice sheet model simulations suggest that Darwin and Hatherton glaciers are strongly buttressed by convergent flow with ice from neighboring Byrd and Mulock glaciers, and by lateral drag past Minna Bluff, which could have led to a pattern of retreat distinct from other glaciers throughout the Transantarctic Mountains.

54 ENVIRONMENTAL SCIENCES↗

A GPU-based compressible combustion solver for applications exhibiting disparate space and time scales

High-speed chemically active flows pose significant computational challenges due to their disparate space and time scales, with stiff chemistry often dominating simulation time. While modern scientific computing programs achieve exascale performance by leveraging graphics processing units (GPUs), existing GPU-based compressible combustion solvers face critical limitations in memory management, load balancing, and handling the highly localized nature of chemical reactions. To this end, we present a high-performance compressible reacting flow solver built on the AMReX framework and optimized for multi-GPU settings. Here, our approach addresses three GPU performance bottlenecks: memory access patterns through column-major storage optimization, computational workload variability via a bulk-sparse integration strategy for chemical kinetics, and multi-GPU load distribution for adaptive mesh refinement applications. The solver adapts existing matrix-based chemical kinetics formulations to multi-grid contexts. Using representative combustion applications, including 2D and 3D detonations and a 3D jet-in-crossflow configuration, we demonstrate 1.4–5× performance improvements over initial implementations on an in-house cluster of NVIDIA H100 GPUs, and near-ideal weak scaling on the Frontier supercomputer (Oak Ridge Leadership Computing Facility) with up to 1024 AMD Instinct MI250X GPUs. Roofline analysis reveals substantial improvements in arithmetic intensity for both convection (∼ 10 ×) and chemistry (∼ 4 ×) routines, confirming efficient utilization of GPU memory bandwidth and computational resources.

42 ENGINEERING↗

Staircase formation by resonant and non-resonant transport of potential vorticity

The E × B staircase is a quasi-periodic pattern of pressure profile corrugations. In this work, we present a new mechanism for E × B staircase formation that involves resonant transport versus non-resonant transport. We start from a potential vorticity evolution system and use quasi-linear theory, a model dispersion relation, and a bi-Lorentzian spectrum approximation, to construct the relation between the fluxes and the profiles. With these fluxes, we close the profile evolution equations and the extended turbulence intensity evolution equation, which together constitute a turbulence-profile evolution system. In this system, the Doppler effect from the E × B mean flow can cause resonance between trapped ion precession motion and the trapped ion mode, which drives a resonant transport contribution to the fluxes. The profiles will be flattened where the resonant transport is switched on. In contrast, for the regions of non-resonant transport, profiles are steeper. A quasi-periodic pattern of profile corrugation (the E × B staircase) spontaneously emerges in this system, which is the two states mentioned above, arranged as alternating layers in space. The feedback processes during the staircase pattern formation are identified. An estimate of the critical value of the boundary heat flux is obtained, above which the staircase formation will be triggered. An estimate scaling of the step size in the staircase pattern is obtained. The resonant turbulent transport is also a mechanism for collisionless saturation of zonal flow. Furthermore, this work is related to internal transport barrier formation and suggests some new scenarios, such as an enhanced confined L mode.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Electric current paths in a Si:P delta-doped device imaged by nitrogen-vacancy diamond magnetic microscopy

The recently-developed ability to control phosphorous-doping of silicon at an atomic level using scanning tunneling microscopy, a technique known as atomic precision advanced manufacturing (APAM), has allowed us to tailor electronic devices with atomic precision, and thus has emerged as a way to explore new possibilities in Si electronics. In these applications, critical questions include where current flow is actually occurring in or near APAM structures as well as whether leakage currents are present. In general, detection and mapping of current flow in APAM structures are valuable diagnostic tools to obtain reliable devices in digital-enhanced applications. In this report, we used nitrogen-vacancy (NV) centers in diamond for wide-field magnetic imaging (with a few-mm field of view and micron-scale resolution) of magnetic fields from surface currents flowing in an APAM test device made of a P delta-doped layer on a Si substrate, a standard APAM witness material. We integrated a diamond having a surface NV ensemble with the device (patterned in two parallel mm-sized ribbons), then mapped the magnetic field from the DC current injected in the APAM device in a home-built NV wide-field microscope. The 2D magnetic field maps were used to reconstruct the surface current densities, allowing us to obtain information on current paths, device failures such as choke points where current flow is impeded, and current leakages outside the APAM-defined P-doped regions. Analysis on the current density reconstructed map showed a projected sensitivity of ~0.03 A m -1 , corresponding to a smallest-detectable current in the 200 μm wide APAM ribbon of ~6 μA. These results demonstrate the failure analysis capability of NV wide-field magnetometry for APAM materials, opening the possibility to investigate other cutting-edge microelectronic devices.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

Experimental Investigation of a Phase Change Material Charged Finned-Tube Heat Exchanger

Understanding different thermal performance aspects of phase change material heat exchangers (PCM-HX) is critical for designing future energy-efficient thermal energy storage systems. This paper presents an experimental performance analysis of a vertically oriented annular finned-tube PCM-HX in a cylindrical container tested under six different melting and solidification test conditions, varied by heat transfer fluid’s (HTF) inlet temperatures and flow directions. The PCM-HX structure with 20 evenly spaced fins was 3D printed in aluminum, and the PCM-HX container was 285 mm in height. The PCM-HX average porosity was 89.7%, where the PCM mass was 246.4 g. Water served as the HTF, and a PCM with a nominal phase change temperature of 35 °C was selected. Instead of opaque conventional insulation sheets, a combination of air-and-argon-filled see-through insulation cylinders was assembled with the PCM-HX to visualize the phase-change phenomena, while minimizing the uncertainties due to heat loss. In total, there were 78 T-type thermocouples installed in the test section which allows a more precise assessment of thermal mass and heat loss, which are often neglected in similar studies found in the literature. We validated the reliability of the test facility as the repeated test results had negligible deviations under the same test conditions. The energy balances between the estimated theoretical and experimentally calculated PCM storage, considering heat loss and thermal mass, were within ±9% deviation for all test conditions. From the experimental investigations, we made several following observations. It was found that the melting pattern of the finned-tube PCM-HX with small gaps between the container wall and the fins was largely affected by the HTF flow directions, unlike the typical straight annular tube PCM-HX. For effective charging and discharging processes, a combination of upward melting and downward solidification is recommended to shorten the test duration. However, an upward solidification should be considered to minimize the cavity formation between the fins caused by the PCM volume contraction during the solidification process.

25 ENERGY STORAGE↗

Wireless Frequency‐Multiplexed Acoustic Array‐Based Acoustofluidics

Abstract Acoustofluidics has shown great potential in enabling on‐chip technologies for driving liquid flows and manipulating particles and cells for engineering, chemical, and biomedical applications. To introduce on‐demand liquid sample processing and micro/nano‐object manipulation functions to wearable and embeddable electronics, wireless acoustofluidic chips are highly desired. This paper presents wireless acoustofluidic chips to generate acoustic waves carrying sufficient energy and achieve key acoustofluidic functions, including arranging particles and cells, generating fluid streaming, and enriching in‐droplet particles. To enable these functions, the wireless acoustofluidic chips leverage mechanisms, including inductive coupling‐based wireless power transfer (WPT), frequency multiplexing‐based control of multiple acoustic waves, and the resultant acoustic radiation and drag forces. For validation, the wirelessly generated acoustic waves are measured using laser vibrometry when different materials (e.g., bone, tissue, and hand) are inserted between the WPT transmitter and receiver. Moreover, the wireless acoustofluidic chips successfully arrange nanoparticles into different patterns, align cells into parallel pearl chains, generate streaming, and enrich in‐droplet microparticles. This research is anticipated to facilitate the development of embeddable wireless on‐chip flow generators, wearable sensors with liquid sample processing functions, and implantable devices with flow generation and acoustic stimulation abilities for engineering, veterinary, and biomedical applications.

Li, Jiali↗

Antecedent Hydrologic Conditions Reflected in Stream Lithium Isotope Ratios During Storms

Antecedent hydrological conditions are recorded through the evolution of dissolved lithium isotope signatures (δ 7 Li) by juxtaposing two storm events in an upland watershed subject to a Mediterranean climate. Discharge and δ 7 Li are negatively correlated in both events,but mean δ 7 Li ratios and associated ranges of variation are distinct between them. We apply a previously developed reactive transport model (RTM) for the site to these event-scale flow perturbations, but observed shifts in stream δ 7 Li are not reproduced. To reconcile the stability of the subsurface solute weathering profile with our observations of dynamic stream δ 7 Li signatures, we couple the RTM to a distribution of fluid transit times that evolve based on storm hydrographs. The approach guides appropriate flux-weighting of fluid from the RTM over a range of flow path lengths, or equivalently fluid residence times. This flux-weighted RTM approach accurately reproduces dynamic storm δ 7 Li-discharge patterns distinguished by the antecedent conditions of the watershed.

58 GEOSCIENCES↗

High-pressure strengthening in ultrafine-grained metals

The Hall–Petch relationship, according to which the strength of a metal increases as the grain size decreases, has been reported to break down at a critical grain size of around 10 to 15 nanometres. As the grain size decreases beyond this point, the dominant mechanism of deformation switches from a dislocation-mediated process to grain boundary sliding, leading to material softening. In one previous approach, stabilization of grain boundaries through relaxation and molybdenum segregation was used to prevent this softening effect in nickel–molybdenum alloys with grain sizes below 10 nanometres. Here we track in situ the yield stress and deformation texturing of pure nickel samples of various average grain sizes using a diamond anvil cell coupled with radial X-ray diffraction. Our high-pressure experiments reveal continuous strengthening in samples with grain sizes from 200 nanometres down to 3 nanometres, with the strengthening enhanced (rather than reduced) at grain sizes smaller than 20 nanometres. Here, we achieve a yield strength of approximately 4.2 gigapascals in our 3-nanometre-grain-size samples, ten times stronger than that of a commercial nickel material. A maximum flow stress of 10.2 gigapascals is obtained in nickel of grain size 3 nanometres for the pressure range studied here. We see similar patterns of compression strengthening in gold and palladium samples down to the smallest grain sizes. Simulations and transmission electron microscopy reveal that the high strength observed in nickel of grain size 3 nanometres is caused by the superposition of strengthening mechanisms: both partial and full dislocation hardening plus suppression of grain boundary plasticity. These insights contribute the ongoing search for ultrastrong metals via materials engineering.

36 MATERIALS SCIENCE↗

Seismic refraction tracks porosity generation and possible CO 2 production at depth under a headwater catchment

Significance To understand groundwater in weathered bedrock aquifers, we explore (bio)geochemical reactions using geophysical measurements. The result leads to a conceptual model for the three-dimensional distribution of geochemical alteration in the subsurface—well beyond the typical approach of inferring (bio)geochemical reactions in the subsurface from pits, boreholes, and outcrops. With this geophysical technique we discovered the prevalence of gas beneath the water table at this site. This geophysical viewing of the subsurface highlights patterns of water movement, reactions of clays, and microbial and abiotic gas-producing reactions recorded in the subsurface and thus illuminates storage, flow, and reactions of groundwater in bedrock systems. Protecting the groundwater resource beneath our feet will rely on such enhanced capacities to see into the subsurface.

58 GEOSCIENCES↗

How MJO Teleconnections and ENSO Interference Impacts U.S. Precipitation

A composite analysis reveals how the Madden–Julian oscillation (MJO) impacts North American rainfall through perturbations in both the upper-tropospheric flow and regional low-level moisture availability. Upper-level divergence associated with the MJO tropical convection drives a quasi-stationary Rossby wave response to the midlatitudes. This forces a midlatitude upper-level dipolar geopotential height anomaly that is accompanied by a westward retraction of the jet stream and reduced rainfall over the central-eastern North Pacific. A reverse effect is found as the MJO propagates eastward across the Maritime Continent. These large differences in the extratropical upper-level flow, combined with anomalies in the regional supply of water vapor, have a profound impact on southeastern U.S. rainfall. The low-frequency variability, including that associated with ENSO, can modify the seasonal background flow (e.g., El Niño and La Niña basic states) affecting the distribution, strength, and propagation of the intraseasonal oscillation and the extratropical teleconnection patterns. The combined effects of the ENSO and the MJO signals result in both spatial and temporal patterns of interference and modulation of North American rainfall. The results from this study show that during a particular phase of an active MJO, the extratropical response can considerably enhance or mask the interannual ENSO signal in the United States, potentially resulting in anomalies of the opposite sign than that expected during a specific ENSO phase. Analyses of specific MJO events during an El Niño or La Niña episode reveal significant contributions to extreme events via constructive and destructive interference of the MJO and ENSO signals.

54 ENVIRONMENTAL SCIENCES↗

User's Manual for the FE/NETL Onshore CO 2 EOR Cost Model, Version 1

This user's manual describes the conceptual and mathematical basis for the FE/NETL Onshore CO 2 EOR Cost Model (a Fortran program). The model performs a cash flow analysis to estimate the cost of implementing CO 2 EOR using supercritical CO 2 by incorporating oil field performance outputs for a pattern from the FE/NETL CO 2 Prophet Model (available on NETL's website along with its associated user's manuals under the Collection Name: FE/NETL CO 2 Prophet Model) and implementing patterns to develop an oil field for CO 2 EOR. The model calculates capital costs, operation and maintenance costs, and financing costs. The user’s manual also describes how to run the FE/NETL Onshore CO 2 EOR Cost Model, along with the model’s file structure, inputs and outputs. The FE/NETL Onshore CO 2 EOR Cost Model is available on NETL's website under the Collection Name: FE/NETL Onshore CO 2 EOR Cost Model.

54 ENVIRONMENTAL SCIENCES↗

Integration of Electromagnetic Geophysics Forward Simulation in Coupled Flow and Geomechanics for Monitoring a Gas Hydrate Deposit Located in the Ulleung Basin, East Sea, Korea

We investigate the feasibility of electromagnetic (EM) geophysics methods to detect the dissociation of gas hydrate specifically from a gas hydrate deposit located in the Ulleung Basin, East Sea, Korea via an integrated flow-geomechanics-EM geophysics simulation. To this end, coupled flow and geomechanics simulation is first performed with the multiple porosity model employed, where a mixed formulation with the finite volume (FV) and finite element (FE) methods are taken for the flow and geomechanics, respectively. From the saturation and porosity fields obtained from the coupled flow and geomechanics, the electrical conductivity model is established for the EM simulation. Solving the partial differential equation of electrical diffusion which is linearized using the 3D finite element method (FEM), the EM fields are then computed. For numerical experiments, particularly two approaches in the configuration for the EM methods are compared in this contribution: the surface-to-surface and the surface-to-borehole methods. When the surface-to-surface EM method is employed, the EM is found to be less sensitive, implying low detectability. Especially for the short term of production, the low detectability is attributed to the similarity of electrical resistivity between the dissociated gas (CH4) and hydrate as well as the specific dissociation pattern within the intercalated composites of the field. On the other hand, when the surface-to-borehole EM method is employed, its sensitivity to capture the produced gas flow is improved, confirming its detectability in monitoring gas flow. Hence, the EM geophysics simulation integrated with coupled flow and geomechanics can be a potential tool for monitoring gas hydrate deposits.

depressurization↗

3-D Geologic Controls of Hydrothermal Fluid Flow at Brady geothermal field, Nevada, USA

In many hydrothermal systems, fracture permeability along faults provides pathways for groundwater to transport heat from depth. Faulting generates a range of deformation styles that cross-cut heterogeneous geology, resulting in complex patterns of permeability, porosity, and hydraulic conductivity. Vertical connectivity (a throughgoing network of permeable areas that allows advection of heat from depth to the shallow subsurface) is rare and is confined to relatively small volumes that have highly variable spatial distribution. This local compartmentalization of connectivity represents a significant challenge to understanding hydrothermal circulation and for exploring, developing, and managing hydrothermal resources. Here, we present an evaluation of the geologic characteristics that control this compartmentalization in hydrothermal systems through 3-D analysis of the Brady geothermal field in western Nevada. A published 3-D geologic map of the Brady area is used as a basis to develop structural and geological variables that are hypothesized to control or effect permeability or connectivity. The 3-D distribution of these variables is compared to the distribution of productive and non-productive fluid flow intervals along production wells and non-productive wells via principal component analysis (PCA). This comparison elucidates which geologic and structural variables are most closely associated with productive fluid flow intervals. Results indicate that production intervals at Brady are located: (1) within or near to known and stress-loaded macro-scale faults, and (2) in areas of high fault and fracture density.

58 GEOSCIENCES↗

3-D Geologic Controls of Hydrothermal Fluid Flow at Brady Geothermal Field, Nevada using PCA

In many hydrothermal systems, fracture permeability along faults provides pathways for groundwater to transport heat from depth. Faulting generates a range of deformation styles that cross-cut heterogeneous geology, resulting in complex patterns of permeability, porosity, and hydraulic conductivity. Vertical connectivity (a through going network of permeable areas that allows advection of heat from depth to the shallow subsurface) is rare and is confined to relatively small volumes that have highly variable spatial distribution. This local compartmentalization of connectivity represents a significant challenge to understanding hydrothermal circulation and for exploring, developing, and managing hydrothermal resources. Here, we present an evaluation of the geologic characteristics that control this compartmentalization in hydrothermal systems through 3-D analysis of the Brady geothermal field in western Nevada. A published 3-D geologic map of the Brady area is used as a basis to develop structural and geological variables that are hypothesized to control or effect permeability or connectivity. The 3-D distribution of these variables is compared to the distribution of productive and non-productive fluid flow intervals along production wells and non-productive wells via principal component analysis (PCA). This comparison elucidates which geologic and structural variables are most closely associated with productive fluid flow intervals. Results indicate that production intervals at Brady are located: (1) within or near to known and stress-loaded macro-scale faults, and (2) in areas of high fault and fracture density. This submission includes the published journal article detailing this work, the published 3-D geologic map of the Brady Geothermal Area used as a basis to develop structural and geological variables that are hypothesized to control or effect permeability or connectivity, 3-D well data, along which geologic data were sampled for PCA analyses, and associated metadata file. This work was done using existing R programs.

15 GEOTHERMAL ENERGY↗

Impacts of Sea‐Level Rise on Coastal Groundwater Table Simulated by an Earth System Model With a Land‐Ocean Coupling Scheme

Abstract Sea‐level rise (SLR) poses a severe threat to the coastal environment through seawater intrusion into freshwater aquifers. The rising groundwater table also exacerbates the risk of pluvial, fluvial, and groundwater flooding in coastal regions. However, current Earth system models (ESMs) commonly ignore the exchanges of water at the land‐ocean interface. To address this gap, we developed a novel land‐ocean hydrologic coupling scheme in a state‐of‐the‐science ESM, the Energy Exascale Earth System Model version 2 (E3SMv2). The new scheme includes the lateral exchange between seawater and groundwater and the vertical infiltration of seawater driven by the SLR‐induced inundation. Simulations were performed with the updated E3SMv2 for the global land‐ocean interface to assess the impacts of SLR on coastal groundwater under a high CO 2 emission scenario. By the middle of this century, seawater infiltration on the inundated areas will be the dominant component in the land‐ocean coupling process, while the lateral subsurface flow exchange will be much smaller. The SLR‐induced seawater infiltration will raise the groundwater levels, enhance evapotranspiration, and increase runoff with distinct spatial patterns globally in the future. Although the coupling process is induced by SLR, we found topography and warming temperature have more control on the coupling impacts, probably due to the relatively modest magnitude of SLR during the selected future period. Overall, our study suggests significant groundwater and seawater exchange at the land‐ocean interface, which needs to be considered in ESMs.

54 ENVIRONMENTAL SCIENCES↗

Patterning Metal Grids for GaAs Solar Cells with Cracked Film Lithography: Quantifying the Cost/Performance Tradeoff

We introduce cracked film lithography (CFL) as a way to reduce the cost of III–V photovoltaics (PV). We spin-coat nanoparticle suspensions onto GaAs thin-film device stacks. The suspensions dry in seconds, forming crack networks that we use as templates through which to electroplate the solar cells’ front metal grids. For the first time, we show that heating the crack template allows it to flow and refill cracks, which decreases crack footprint and improves final grid transmittance. We demonstrate 24.7%-efficient single-junction GaAs solar cells using vacuum-free CFL grids. These devices are only 1.7% (absolute) less efficient than the baseline grids patterned by photolithography with the loss mostly resulting from the reduced transparency of the CFL pattern. Additional optimization could decrease this difference. Initial cost modeling suggests that CFL is more scalable than photolithography: In particular, CFL’s lower materials and equipment costs could greatly reduce the levelized cost of electricity of III–V PV at scale, a potential step toward terrestrial deployment.

14 SOLAR ENERGY↗

Development of inter-grid-cell lateral unsaturated and saturated flow model in the E3SM Land Model (v2.0)

Abstract. The lateral transport of water in the subsurface is important in modulating terrestrial water energy distribution. Although a few land surface models have recently included lateral saturated flow within and across grid cells, it is not a default configuration in the Climate Model Intercomparison Project version 6 experiments. In this work, we developed the lateral subsurface flow model within both unsaturated and saturated zones in the Energy Exascale Earth System Model (E3SM) Land Model version 2 (ELMv2.0). The new model, called ELMlat, was benchmarked against PFLOTRAN, a 3D subsurface flow and transport model, for three idealized hillslopes that included a convergent hillslope, divergent hillslope, and tilted V-shaped hillslope with variably saturated initial conditions. ELMlat showed comparable performance against PFLOTRAN in terms of capturing the dynamics of soil moisture and groundwater table for the three benchmark hillslope problems. Specifically, the mean absolute errors (MAEs) of the soil moisture in the top 10 layers between ELMlat and PFLOTRAN were within 1 %±3 %, and the MAEs of water table depth were within ±0.2 m. Next, ELMlat was applied to the Little Washita experimental watershed to assess its prediction of groundwater table, soil moisture, and soil temperature. The spatial pattern of simulated groundwater table depth agreed well with the global groundwater table benchmark dataset generated from a global model calibrated with long-term observations. The effects of lateral groundwater flow on the energy flux partitioning were more prominent in lowland areas with shallower groundwater tables, where the difference in simulated annual surface soil temperature could reach 0.3–0.4 ∘C between ELMv2.0 and ELMlat. Incorporating lateral subsurface flow in ELM improves the representation of the subsurface hydrology, which will provide a good basis for future large-scale applications.

58 GEOSCIENCES↗

Size and Distribution of Parr Produced from Natural‐ and Hatchery‐Origin Steelhead Spawning Naturally in a Small Pacific Northwest Coastal Stream

Abstract Recent studies suggest that steelhead Oncorhynchus mykiss produced from local hatchery‐origin (HOR) adults have lower lifetime fitness than their natural‐origin (NOR) counterparts. To increase our understanding of this pattern, we compared age‐1 parr size and distribution produced by a local integrated population of HOR and NOR adults spawning naturally across a range of environmental conditions. Across 8 years, we used genetic parentage assignments and field measurements of parr in conjunction with creek temperature and flow data to find small, nonbiologically significant differences in size between parr of different parent origins. This suggests that parr produced by HOR adults are acquiring enough resources to grow at a rate similar to that of parr produced by NOR adults. In contrast to origin, we found strong positive associations between mean size of parr and annual mean water temperature and summer flow. Additionally, the distribution of parr was similar and HOR‐produced parr were not skewed relative to the location of the hatchery. Parr occupying the full extent of accessible and suitable habitat might be a positive outcome for hatchery programs seeking to supplement an existing population or replace an extirpated population. However, these results also imply that HOR fish could be competing with NOR fish for food and space, which might be less desirable if the goal of the hatchery program is harvest and not supplementation. Our results highlight the importance of clearly articulated goals for the hatchery program, as the observed pattern could be deemed a benefit or a risk depending on the hatchery's purpose. Lastly, projections for higher water temperatures and reduced summer flows, when considered in the context of the correlations we observed for size of parr, imply that coastal steelhead populations could start experiencing negative climate impacts with respect to steelhead parr growth metrics during the summer in the next few decades.

Kennedy, Benjamen M.↗