Search NASA⌕ Search

SEARCH · Search NASA

Results for “Way”

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 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Signatures of a Tidally Induced Spiral Arm at the Anticenter of the Milky Way and a Kinematically Extended Anticenter Stream Using DESI Data Release 2

Using the Dark Energy Spectroscopic Instrument (DESI) Milky Way Survey, we examine the six-dimensional space of the anticenter region of the Milky Way stellar disk (150° < Galactic longitude < 220°) using 61,883 main-sequence turnoff stars. We focus on two well-known stellar overdensities in the anticenter: the Monoceros Ring (MRi) and Anticenter Stream (ACS). We find that the MRi overdensity has kinematic signatures consistent with a tidally induced spiral arm, a type of dynamic spiral arm created by an interaction with a satellite galaxy, most likely the Sagittarius dwarf spheroidal galaxy (Sgr). We use the kinematics of the MRi to calculate the two most recent passage times of Sgr, finding 0.25 ± 0.09 Gyr and 1.10 ± 0.23 Gyr from the present day. We validate that the ACS is kinematically decoupled from the MRi because they are moving in opposite radial and vertical directions. We find that the kinematics associated with the ACS extends beyond our defined overdensity. The features we see in the ACS region are likely part of a broader distribution of stars with the same kinematic signature as detected in other places, like the vertical wave in the outer disk and phase spiral.

Lambert, Mika [University of California, Santa Cru↗
DOE OSTI2026

The Backup Program of the Dark Energy Spectroscopic Instrument’s Milky Way Survey

The Milky Way Backup Program (MWBP), a survey currently underway with the Dark Energy Spectroscopic Instrument (DESI) on the Nicholas U. Mayall 4 m Telescope, works at the margins of the DESI Main surveys to obtain spectra of millions of additional stars from the Gaia catalog. Efficiently utilizing times between ∼12° and 18° twilight and poor weather conditions, the MWBP extends the range of stellar sources studied to both brighter magnitudes and lower Galactic latitude and declination than the stars studied in DESI’s Main Milky Way Survey. While the MWBP prioritizes candidate giant stars selected from the Gaia catalog (using color and parallax criteria), it also includes an unbiased sample of bright stars (i.e., 11.2 ≲ G < 16 mag) as well as fainter sources (to G ≲ 19 mag). As of 2025 March 1, the survey had obtained spectra of ∼7 million stars, approximately 1.2 million of which are included in the DESI Data Release 1. The DESI spectra cover the wavelength range from 3600 to 9800 Å at a resolution λ/Δλ varying from 2000 to 5000. The full survey, when completed, will cover an area of more than 21,000 deg 2 and include approximately 10 million Gaia sources, roughly equal to the number of stellar spectra obtained through the DESI Main Survey, while only utilizing ≈9% of all DESI observing time. This paper provides an overview of the MWBP, describing the target selection, observing strategy, and an introduction to the resulting data.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Applying Noether's theorem to matter in the Milky Way: evidence for external perturbations and non-steady-state effects from Gaia Data Release 2

We apply Noether's theorem to observations of main-sequence stars from the Gaia Data Release 2 archive to probe the matter distribution function of the Galaxy. That is, we examine the axial symmetry of stars at vertical heights $z$, $0.2 \le |z| \le 3$ kpc, to probe the quality of the angular momentum $L_z$ as an integral of motion. The failure of this symmetry test would speak to a Milky Way, in both its visible and dark matter, that is not isolated and/or not in steady state. The left-right symmetry-breaking pattern we have observed, north and south, reveals both effects, with a measured deviation from symmetry of typically 0.5%. We show that a prolate form of the gravitational distortion of the Milky Way by the Large Magellanic Cloud, determined from fits to the Orphan stream by Erkal et al., 2019, is compatible with the size and sign of the axial-symmetry-breaking effects we have discovered in our sample of up to 14.4 million main-sequence stars, speaking to a distortion of an emergent, rather than static, nature.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Seeking Self-regulating Simulations of Idealized Milky Way–like Galaxies

Abstract Precipitation is potentially a mechanism through which the circumgalactic medium (CGM) can regulate a galaxy’s star formation. Here, we present idealized simulations of isolated Milky Way–like galaxies intended to examine the ability of galaxies to self-regulate their star formation, in particular via precipitation. We also examine the impact of rotation in the CGM. Using six simulations, we explore variations in the initial CGM t cool / t ff ratio and rotation profile. Those variations affect the amount of gas accretion and star formation within the galactic disk. To encourage this accretion and better study its dependence on CGM structure, we gradually increase the efficiency of stellar feedback during the first half of our simulations. Yet despite this gradual increase, the resulting outflows quickly evacuate large, hot cavities within the CGM and even beyond r 200 . Some of the CGM gas avoids interacting with the cavities and is able to feed the disk along its midplane, but the cooling of feedback-heated gas far from the midplane is too slow to supply the disk with additional gas. Our simulations illustrate the importance of physical mechanisms in the outer CGM and IGM for star formation regulation in Milky Way–scale halos.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Chemical Doppelgangers in GALAH DR3: The Distinguishing Power of Neutron-capture Elements among Milky Way Disk Stars

The observed chemical diversity of Milky Way stars places important constraints on Galactic chemical evolution and the mixing processes that operate within the interstellar medium. Recent works have found that the chemical diversity of disk stars is low. For example, the Apache Point Observatory Galactic Evolution Experiment (APOGEE) "chemical doppelganger rate," or the rate at which random pairs of field stars appear as chemically similar as stars born together, is high, and the chemical distributions of APOGEE stars in some Galactic populations are well-described by two-dimensional models. However, limited attention has been paid to the heavy elements (Z > 30) in this context. In this work, we probe the potential for neutron-capture elements to enhance the chemical diversity of stars by determining their effect on the chemical doppelganger rate. We measure the doppelganger rate in GALactic Archaeology with HERMES DR3, with abundances rederived using The Cannon, and find that considering the neutron-capture elements decreases the doppelganger rate from ~2.2% to 0.4%, nearly a factor of 6, for stars with –0.1 < [Fe/H] < 0.1. While chemical similarity correlates with similarity in age and dynamics, including neutron-capture elements does not appear to select stars that are more similar in these characteristics. Our results highlight that the neutron-capture elements contain information that is distinct from that of the lighter elements and thus add at least one dimension to Milky Way abundance space. This work illustrates the importance of considering the neutron-capture elements when chemically characterizing stars and motivates ongoing work to improve their atomic data and measurements in spectroscopic surveys.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

The SAGA Survey. IV. The Star Formation Properties of 101 Satellite Systems around Milky Way–mass Galaxies

We present the star-forming properties of 378 satellite galaxies around 101 Milky Way analogs in the Satellites Around Galactic Analogs (SAGA) Survey, focusing on the environmental processes that suppress or quench star formation. In the SAGA stellar mass range of 10 6−10 M ⊙ , we present quenched fractions, star-forming rates, gas-phase metallicities, and gas content. The fraction of SAGA satellites that are quenched increases with decreasing stellar mass and shows significant system-to-system scatter. SAGA satellite quenched fractions are highest in the central 100 kpc of their hosts and decline out to the virial radius. Splitting by specific star formation rate (sSFR), the least star-forming satellite quartile follows the radial trend of the quenched population. The median sSFR of star-forming satellites increases with decreasing stellar mass and is roughly constant with projected radius. Star-forming SAGA satellites are consistent with the star formation rate–stellar mass relationship determined in the Local Volume, while the median gas-phase metallicity is higher and median H i gas mass is lower at all stellar masses. We investigate the dependence of the satellite quenched fraction on host properties. Quenched fractions are higher in systems with larger host halo mass, but this trend is only seen in the inner 100 kpc; we do not see significant trends with host color or star formation rate. Our results suggest that lower-mass satellites and satellites inside 100 kpc are more efficiently quenched in a Milky Way–like environment, with these processes acting sufficiently slowly to preserve a population of star-forming satellites at all stellar masses and projected radii.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Stellar Metallicities from DECam $u$-band Photometry: A Study of Milky Way Ultra-Faint Dwarf Galaxies

We conducted an in-depth analysis of candidate member stars located in the peripheries of three ultra-faint dwarf (UFD) galaxy satellites of the Milky Way: Boötes I (Boo1), Boötes II (Boo2), and Segue I (Seg1). Studying these peripheral stars has previously been difficult due to contamination from the Milky Way foreground. We used $u$-band photometry from the Dark Energy Camera (DECam) to derive metallicities to efficiently select UFD candidate member stars. This approach was validated on Boo1, where we identified both previously known and new candidate member stars beyond five half-light radii. We then applied a similar procedure to Boo2 and Seg1. Our findings hinted at evidence for tidal features in Boo1 and Seg1, with Boo1 having an elongation consistent with its proper motion and Seg1 showing some distant candidate stars, a few of which are along its elongation and proper motion. We find two Boo2 stars at large distances consistent with being candidate member stars. Using a foreground contamination rate derived from the \emph{Besan\c{c}on} Galaxy model, we ascribed purity estimates to each candidate member star. We recommend further spectroscopic studies on the newly identified high-purity members. Our technique offers promise for future endeavors to detect candidate member stars at large radii in other systems, leveraging metallicity-sensitive filters with the Legacy Survey of Space and Time and the new, narrow-band Ca HK filter on DECam.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Inference of Multichannel r -process Element Enrichment in the Milky Way Using Binary Neutron Star Merger Observations

Observations of GW170817 strongly suggest that binary neutron star (BNS) mergers produce rapid neutron-capture nucleosynthesis ( r -process) elements. However, it remains an open question whether these mergers can account for all the r -process element enrichment in the Milky Way’s history. Here, we constrain the contributions of the BNS channel using astrophysical neutron star observations. The rate and mass distributions are constrained by LIGO/Virgo/Kagra through the latest catalog GWTC-3, the neutron star equation of state by gravitational-wave, radio, and X-ray observations, and the delay time distribution by short gamma-ray burst (GRB) host galaxy associations. We present a Bayesian framework to consistently combine these observations with abundance information to quantify the contribution and uncertainties of single and multiple astrophysical enrichment sources, and obtain a distribution of per-event BNS r -process element yields consistent with geophysical and astrophysical abundance constraints. We then adopt a Galactic chemical evolution model assuming an instantaneous and fixed amount of Fe enrichment from core-collapse supernovae, and show that BNS-only enrichment scenarios remain inconsistent with the observed r-process abundance trend of disk stars in the Galaxy even with the uncertainties in BNS merger observations. Using stellar abundance observations instead of the short GRB constraints, we can infer a shorter BNS delay time distribution with power-law index α ≤ −2.0 and minimum delay time ${t}_{{\rm{\min }}}\leqslant 40$ Myr at 90% confidence, consistent with detailed Galactic chemical evolution models. Such delay times are in tension with those predicted by standard BNS formation models. Alternatively, we confirm that a two-channel scenario, in which the second channel tracks the star formation history without significant delay, can account for both Galactic stellar and short GRB observations. We estimate that 45%–90% of the r -process abundance in the Milky Way today would have been produced by this star formation-tracking channel, rather than BNS mergers with significant delay times.

gravitational wave astronomy↗
DOE OSTI2023

Improved Particle Heat Transfer by way of Bimodal Particle Distributions for High Temperature Solar Thermal Energy

High temperature solar thermal facilities are looking to increase operating temperatures through novel heat transfer media, one such being solid particles. These particles operating at high temperatures will require transferring their thermal energy into another working fluid like supercritical carbon dioxide which can be used in advanced power cycles. Achieving high heat transfer between the particles and supercritical carbon dioxide is essential to high efficiency and low-cost operation. Therefore, optimizing the thermal conductivity of these particles is one potential way to ensure high performance. Traditionally, unimodal particle distributions have been employed in high temperature particle solar power plants. However, ambient temperature testing of bimodal particle distributions has revealed a superior thermal conductivity when compared to its unimodal counterpart at the same temperature. This data was obtained by certified, off-the-shelf instruments that can effectively simulate the conditions a particle would be exposed to in a high temperature solar thermal system. Data obtained in this way suggests that the increased thermal conductively imputed by a bimodal particle distribution is significant at working temperatures in solar facilities. Furthermore, the thermal conductivity of these bimodal particle distributions peaks when the best combination of large and small particles is applied. At high temperatures, binary particle distributions are compared to monodispersed distributions of larger particles where heat transfer is more prolific due to the increased surface radiation. Various thermal conductivity, porosity and heat exchanger models are explored in conjunction with data acquired up to 700 C.

Stout, Dallin (ORCID:0009000294586091)↗
DOE OSTI2026

Discovery and Spectroscopic Characterization of a Distant, Compact Milky Way Satellite in Gemini

We present the discovery of a compact Milky Way satellite in the constellation of Gemini. This system was discovered by cross-matching detections from two independent search algorithms applied to Blanco/DECam data from the third data release of the DECam Local Volume Exploration survey (DELVE DR3), and confirmed with deeper imaging from Gemini/GMOS-N. Based on these data, we determine that the system is an ultra-faint ($M_V = -2.1^{+0.4}_{-0.6}$), compact ($r_{1/2} = 8.6^{+1.4}_{-1.2}$ pc) system located at a heliocentric distance of $120^{+7}_{-6}$ kpc. These physical properties place the system in the regime of ambiguous, ultra-faint compact Milky Way halo satellites that cannot be confidently classified as dwarf galaxies or star clusters from morphology alone; we therefore name the system DELVE 8/Gemini I. From medium-resolution Keck/DEIMOS spectroscopy, we securely identify four members including two blue horizontal branch stars, confirming the system as a bound satellite moving at a mean radial velocity of $v_{\rm hel} = -82.7^{+3.7}_{-3.9} {\rm km\,s}^{-1}$. We also use these spectra to place an upper limit of $\rm [Fe/H] \lesssim -2.5$ on the metallicity of DELVE 8/Gemini I's brightest star, supporting the classification of the system as either an ancient star cluster or ultra-faint dwarf galaxy. The discovery of faint, distant systems similar to DELVE 8/Gemini I is expected to become more common with upcoming surveys.

Overdeck, K. [Chicago U., Astron. Astrophys. Ctr.;↗
DOE OSTI2020

Constraining the Milky Way Mass Profile with Phase-space Distribution of Satellite Galaxies

We estimate the Milky Way (MW) halo properties using satellite kinematic data including the latest measurements from Gaia DR2. With a simulation-based 6D phase-space distribution function (DF) of satellite kinematics, we can infer halo properties efficiently and without bias, and handle the selection function and measurement errors rigorously in the Bayesian framework. Applying our DF from the EAGLE simulation to 28 satellites, we obtain an MW halo mass of $M={1.23}_{-0.18}^{+0.21}\times {10}^{12}{M}_{\odot }$ and a concentration of $c={9.4}_{-2.1}^{+2.8}$ with the prior based on the M–c relation. The inferred mass profile is consistent with previous measurements but with better precision and reliability due to the improved methodology and data. Potential improvement is illustrated by combining satellite data and stellar rotation curves. Using our EAGLE DF and best-fit MW potential, we provide much more precise estimates of the kinematics for those satellites with uncertain measurements. Compared to the EAGLE DF, which matches the observed satellite kinematics very well, the DF from the semi-analytical model based on the dark-matter-only simulation Millennium II (SAM-MII) over-represents satellites with small radii and velocities. We attribute this difference to less disruption of satellites with small pericenter distances in the SAM-MII simulation. Finally, by varying the disruption rate of such satellites in this simulation, we estimate a ~5% scatter in the inferred MW halo mass among hydrodynamics-based simulations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Two-point Correlation Function Studies for the Milky Way: Discovery of Spatial Clustering from Disk Excitations and Substructure

We introduce a two-particle correlation function (2PCF) for the Milky Way, constructed to probe spatial correlations in the orthogonal directions of the stellar disk in the Galactic cylindrical coordinates of R, ϕ, and z. We use this new tool to probe the structure and dynamics of the Galaxy using the carefully selected set of solar neighborhood stars (d ≲ 3 kpc) from Gaia Data Release 2 that we previously employed for studies of axial symmetry breaking in stellar number counts. We make additional, extensive tests, comparing to reference numerical simulations, to ensure our control over possibly confounding systematic effects. Supposing either axial or north–south symmetry, we divide this data set into two nominally symmetric sectors and construct the 2PCF, in the manner of the Landy–Szalay estimator, from the Gaia data. In so doing, working well away from the midplane region in which the spiral arms appear, we have discovered distinct symmetry-breaking patterns in the 2PCF in its orthogonal directions, thus establishing the existence of correlations in stellar number counts alone at subkiloparsec length scales for the very first time. In particular, we observe extensive wavelike structures of amplitude greatly in excess of what we would estimate if the system were in a steady state. We study the variations in these patterns across the Galactic disk, and with increasing ∣z∣, and we show how our results complement other observations of non-steady-state effects near the Sun, such as vertical asymmetries in stellar number counts and the Gaia snail.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Angular Momentum Variation of the Milky Way Thick Disk: The Dependence of Chemical Abundance and Evidence of the Inside-out Formation Scenario

We investigate the angular momentum of mono-abundance populations (MAPs) of the Milky Way thick disk by using a sample of 26,076 giant stars taken from APOGEE Data Release (DR) 17 and Gaia early DR3. The vertical and perpendicular angular momentum components, L Z and L P , of the MAPs in narrow bins have significant variations across the [α/M]–[M/H] plane. L Z and L P systematically change with [M/H] and [α/M] and can be alternatively quantified by the chemical gradients: d[M/H]/dL Z = 1.2 × 10 -3 dex kpc -1 km -1 s, d[M/H]/dL P = -5.0 × 10 -3 dec kpc -1 km -1 s, and d[α/M]/dL Z = -3.0 × 10 -4 dex kpc -1 km -1 s, d[α/M]/dL P = 1.2 × 10 -3 dec kpc -1 km -1 s. These correlations can also be explained as the chemical dependence of the spatial distribution shape of the MAPs. We also exhibit the corresponding age dependence of the angular momentum components. Under the assumption that the guiding radius (R g ) is proportional to L Z , this provides direct observational evidence of the inside-out structure formation scenario of the thick disk, with dR g /dAge = -1.9 kpc Gyr -1 . The progressive changes in the disk thickness can be explained by the upside-down formation or/and the consequent kinematical heating.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Simulating Atomic Dark Matter in Milky Way Analogs

Dark sector theories naturally lead to multicomponent scenarios for dark matter where a subcomponent can dissipate energy through self-interactions, allowing it to efficiently cool inside galaxies. We present the first cosmological hydrodynamical simulations of Milky Way analogs where the majority of dark matter is collisionless cold dark matter (CDM) but a subcomponent (6%) is strongly dissipative minimal atomic dark matter (ADM). The simulations, implemented in GIZMO and utilizing FIRE-2 galaxy formation physics to model the standard baryonic sector, demonstrate that the addition of even a small fraction of dissipative dark matter can significantly impact galactic evolution despite being consistent with current cosmological constraints. We show that ADM gas with roughly standard model–like masses and couplings can cool to form a rotating “dark disk” with angular momentum closely aligned with the visible stellar disk. The morphology of the disk depends sensitively on the parameters of the ADM model, which affect the cooling rates in the dark sector. The majority of the ADM gas gravitationally collapses into dark “clumps” (regions of black hole or mirror star formation), which form a prominent bulge and a rotating thick disk in the central galaxy. These clumps form early and quickly sink to the inner ∼kiloparsec of the galaxy, affecting the galaxy’s star formation history and present-day baryonic and CDM distributions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

The effects of SO 2 exposure on three-way catalysts for passive SCR operation

We report efficient lean-burn gasoline engines cannot employ traditional three-way catalysts (TWCs) to control NO X emissions due to excess oxygen. Passive selective catalytic reduction (pSCR) is a process designed to control NOx by generating NH 3 over the TWC during fuel-rich operation, storing it on a downstream SCR, and then converting NOx during fuel-lean operation. This work is focused on studying sulfur effects on this promising strategy. Two hydrothermally-aged commercial catalysts were studied for their effectiveness in the presence of SO 2 : a Pd-based TWC and a TWC with NO X storage. The sulfated samples were evaluated in a flow reactor under pSCR conditions and with targeted reactions. SO 2 inhibited the water gas shift and steam reforming reactions, decreasing (but not eliminating) NH 3 production; however, initial activity is recoverable after a simple thermal treatment at 650 °C. This work will prove valuable in understanding, mitigating, and even reversing the effects of SO 2 inhibition on real pSCR systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

CH 4 combustion over a commercial Pd/CeO 2 -ZrO 2 three-way catalyst: Impact of thermal aging and sulfur exposure

Thermal aging and sulfur poisoning are major problems influencing the efficiency and lifetime of natural gas three-way catalysts (TWCs) during real-world operation. In the present study, thermal aging and sulfur-induced deactivations on CH 4 combustion were investigated over commercial monolithic honeycombs, i.e., Pd/CexZr 1-x O 2 (CZ)-based TWCs that are thermally aged in the bench (HTA), followed by being operated behind the stoichiometric natural gas engine in high sulfur fuel (HTA + S). In this work, both HTA and HTA + S samples show degraded performance in CH 4 combustion relative to the fresh one, which are caused by the reconstruction of Pd species induced by thermal aging and the decreased reducibility of the support. The 18 O 2 labelling experiments show that in all cases CH 4 combustion proceeds via a MvK mechanism in which lattice O (OL) plays a key role. The sulfation aging (HTA + S sample) promotes the activation of CH 4 and O 2 by forming Pd δ+ -(SO$_4^{2-}$) δ- and Ce 3+ -VO couples at the interface as compared to HTA sample. However, such a promotion effect could be compromised by the blocking effect of sulfate on the exchange of O L from the bulk to the surface and across the surface. The fundamental understanding of thermal and S deactivation mechanisms from this study could help to predict a real-world lifetime use curve of natural gas TWCs.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Computational study of the ECN spray C via one-way coupling of internal nozzle flow and ensuing spray

Here, in this study, computational fluid dynamics simulations of the Engine Combustion Network Spray C injector were performed to investigate the effect of the internal nozzle flow on the development of the ensuing spray in automotive fuel injectors. A set of simulation best practices was first developed to achieve an accurate representation of the internal nozzle flow. Detailed, spatiotemporally defined information was extracted at the injector’s orifice exit and stored in maps that were subsequently employed to initialize a Lagrangian spray by means of a static, one-way coupling (OWC) approach. To highlight the impact of the in-nozzle flow on the spray, the OWC approach was compared with the conventional rate-of-injection (ROI) approach, for which the spray parcels are commonly initialized using a blob injection model. The results indicated that, using the ROI approach, considerable evaporation was observed within the plume in the near-nozzle region. On the other hand, the OWC approach showed evidence of evaporation farther downstream and at the plume periphery. The different behavior observed with the two approaches was mostly due to substantial differences in fuel mass distribution. The liquid fuel concentration predicted by the OWC approach matched well with X-ray tomography measurements of projected fuel mass per unit volume, in terms of both magnitude and distribution. In contrast, an unphysical distribution of fuel mass characterized by an extremely low liquid concentration was predicted by the ROI approach. In the far-field, the effect of the in-nozzle flow was found to be weaker as the spray penetrated and interacted with the ambient gas. In this region, the OWC and ROI approaches provided similar predictions of vapor penetration, spray morphology, and mixing-related characteristics.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2024

Two-Way Option Contracts That Facilitate Adaptive Water Reallocation in the Western United States

Many water markets in the western United States (U.S.) have the ability to reallocate water temporarily during drought, often as short-term water rights leases from lower value irrigated activities to higher value urban uses. Regulatory approval of water transfers, however, typically takes time and involves high transaction costs that arise from technical and legal analyses, discouraging short-term leasing. This leads municipalities to protect against drought-related shortfalls by purchasing large volumes of infrequently used permanent water rights. High transaction costs also result in municipal water rights rarely being leased back to irrigators in wet or normal years, reducing agricultural productivity. This research explores the development of a multi-year two-way option (TWO) contract that facilitates leasing from agricultural-to-urban users during drought and leasing from urban-to agricultural users during wet periods. The modeling framework developed to assess performance of the TWO contracts includes consideration of the hydrologic, engineered, and institutional systems governing the South Platte River Basin in Colorado where there is growing competition for water between municipalities (e.g., the city of Boulder) and irrigators. The modeling framework is built around StateMod, a network-based water allocation model used by state regulators to evaluate water rights allocations and potential rights transfers. Results suggest that the TWO contracts could allow municipalities to maintain supply reliability with significantly reduced rights holdings at lower cost, while increasing agricultural productivity in wet and normal years. Additionally, the TWO contracts provide irrigators with additional revenues via net payments of option fees from municipalities.

54 ENVIRONMENTAL SCIENCES↗