Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “OBSERVATORY”

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 559 records · Page 31
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Plant-Soil Relationships Influence Observed Trends Between Manganese and Carbon Across Biomes

Manganese (Mn) is an essential plant micronutrient that plays a critical role in the litter decomposition by oxidizing and degrading complex organic molecules. Previous studies report a negative correlation between Mn concentrations and carbon (C) storage in organic horizons and suggest that high Mn concentrations in leaf litter reduce soil C storage in forest ecosystems, presumably by stimulating the oxidation of lignin by fungal enzymes. Yet, the relationship between Mn and C in the litter layer and organic soil remains poorly understood and restricted to a few biomes, hampering our ability to improve mechanistic understanding of soil C accumulation. Here, to examine plant-soil interactions that underlie observed relationships between Mn and C across a wide range of biomes, we extracted biogeochemical data reported for plants and soils from the National Ecological Observatory Network (NEON) database. We found that increased C and nitrogen (N) storage in organic horizons were associated with declines in Mn concentrations across diverse ecosystems at the continental scale, and this relationship was associated with the degree of organic matter decomposition (i.e., O i , O e , and O a ). Carbon and N stocks were more strongly correlated with Mn than with climatic variables (i.e., temperature and precipitation). Foliar Mn was strongly correlated with foliar lignin, and both these parameters increased with a decrease in soil pH, indicating links between soil pH, foliar chemistry, and litter decomposability. Our observations suggest that increased Mn bioavailability and accumulation in foliage under moderately acidic soil conditions support fungal decomposition of lignin-rich litter and contributes to lower soil C stocks.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Rapid Refreezing of a Marginal Ice Zone Across a Seafloor Distributed Acoustic Sensor

Abstract Seismic instruments such as broadband seismometers and distributed acoustic sensors (DAS) have a demonstrated potential for wide‐scale and continuous in situ monitoring of near‐surface environmental and anthropogenic processes. DAS is attractive for development as a multi‐geophysical observatory due to the prevalence of existing fiber infrastructure in regions with environmental, cultural, or strategic significance. Here we present results from a multi‐seasonal acquisition of DAS data on a seafloor cable in the Beaufort Sea, Alaska. During a November 2021 data collect we captured the rapid transition of ambient noise characteristics from an “ice‐free” state to an “ice‐bound” state. A sea ice formation front was plainly visible on the DAS record and was observed to propagate 20 km seaward over a period of 8 hr. Satellite‐based instrumentation were unable to record this event due to cloud cover, low light conditions, and orbital frequency.

Baker, Michael G.↗
DOE OSTI2022

Magnetotelluric Investigations of the Kilauea Volcano, Hawaii

In 2002 and 2003 a collaborative effort was undertaken between Lawrence Berkeley National Laboratory, Sandia National Laboratories, the U.S. Geological Survey (USGS) Menlo Park, the USGS Hawaiian Volcano Observatory, and Electromagnetic Instruments Inc. to study the Kilauea volcano in Hawaii using the magnetotelluric (MT) technique. The work was motivated by a desire to improve understanding of the magma reservoirs and conduits within Kilauea and the East and Southwest Rift zones, which has implications for understanding Kilauea's plumbing system. An improved understanding of the rift zones has implications in understanding large-scale landslides that are generated in the Hilina Slump, which produce significant impacts on coastal communities. Up to eight stations operated simultaneously, with multiple remote reference sites, and data were processed using multi-station robust processing techniques. Data was acquired at 70 sites over the Southwest and East rift zones. Good to excellent quality data were obtained even in the harshest conditions, such as those encountered on the fresh lava flows of the East Rift Zone, where electrical contact resistances are on the order of 100 kΩ. A three-dimensional (3D) MT model study was done to guide interpretation of the observed MT measurements. Synthetic modeling demonstrates that conductive bodies in the upper 3 km can be spatially resolved where MT station sampling is good. Resistivity anomalies in the 3D inversions have a high degree of spatial correlation with previously published seismic velocity anomalies beneath Kilauea. Melt fractions between 0.096 and 0.117 are calculated for the Kilauea and Pu'u'ō'ō low resistivity anomalies, respectively.

58 GEOSCIENCES↗
DOE OSTI2022

Influence of Aerosol Embedded in Shallow Cumulus Cloud Fields on the Surface Solar Irradiance

Ubiquitous shallow cumulus clouds are associated with complex variability in surface solar irradiance (SSI). Aerosol embedded in the cloud field typically has a much smaller overall radiative effect, but can significantly perturb the shape of the SSI probability density function (PDF). These perturbations have important implications for several applications that utilize SSI, but are poorly quantified and are the subject of this study. Multiple cases of shallow cumulus cloud fields with embedded aerosol are simulated at the Southern Great Plains Atmospheric Observatory using large eddy simulation (LES). The LES-derived cloud and aerosol fields are then ingested into Monte Carlo three-dimensional (3D) radiative transfer to simulate SSI. We find a variety of perturbations to the SSI PDF that depend on aerosol presence and optical properties. The processes leading to these perturbations include extinction of the direct beam that often increases from the clear-sky region toward cloud edge due to aerosol hygroscopic growth, and scattering of radiation by aerosol into cloud shadows. The ability to predict the SSI PDF in the presence of aerosol is assessed by adding three representative aerosol optical properties into an existing machine learning framework. We show that machine learning accurately predicts the SSI PDF across a wide range of conditions with negligible computational expense. Importance metrics reveal the relatively high influence of aerosol optical properties in making the predictions. These new findings highlight the important role that aerosol plays in SSI variability for highly 3D cloud-aerosol environments and provides a computationally efficient route forward for its simulation.

3D radiative transfer↗
DOE OSTI2023

Look Up: Probing the Vertical Profile of New Particle Formation and Growth in the Planetary Boundary Layer With Models and Observations

Abstract The processes of new particle formation (NPF) and growth are important contributors to cloud condensation nuclei (CCN) concentrations, and CCN are important for climate from their impact on planetary radiative forcing. While the general ubiquity and importance of NPF is understood, the vertical extent and governing mechanisms of NPF and growth in the lower troposphere are uncertain. We present an analysis of four NPF events and two non‐NPF events during the HI‐SCALE field campaign at the Southern Great Plains observatory in Oklahoma, USA. First, we analyzed airborne and ground‐based observations of aerosol and gas‐phase properties. Second, we used a column aerosol chemistry and microphysics model to probe factors that influence the vertical profile of NPF. During HI‐SCALE, we found several instances of enhanced NPF occurring several hundred meters above the surface; however, the spatio‐temporal characteristics of the observed NPF made comparisons between airborne‐ and ground‐based observations difficult. The model represented the observed NPF (or lack of NPF) and particle growth at the surface to final diameters within 10 nm. The model predicted enhanced NPF rates in the upper mixed layer, and this enhancement is primarily due to the temperature dependence in the NPF schemes, but this was also dependent on the vertical profile of gas‐phase precursors measured during HI‐SCALE. We found vertical mixing in the model either enhanced or suppressed NPF rates, aerosol number concentrations, and particle growth rates at the surface. Finally, our analysis provides insights for future field campaigns and modeling efforts investigating the vertical profile of NPF.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

A Continental-Scale Estimate of Soil Organic Carbon Change at NEON Sites and Their Environmental and Edaphic Controls

Current carbon cycle models focus on the effects of climate and land-use change on primary productivity and microbial-mineral dependent carbon turnover in the topsoil, while less attention has been paid to vertical soil processes and soil-dependent response to land-use change along the profile. In this study, a spatial-temporal analysis was used to estimate soil organic carbon (SOC) change in topsoil/A horizon and subsoil/B horizon at National Ecological Observatory Network (NEON) sites, USA over 30 years. To separate the effects of land-use, environmental, and edaphic factors on SOC change, space-for-time substitution was used in combination with the Continuous Change Detection and Classification algorithm and Structural Equation Modeling. Results showed that (a) under natural vegetation, Spodosols and Inceptisols found in the eastern NEON sites had substantial topsoil SOC accumulation (+0.4 to +1.2 Mg C ha –1 year –1 ), while Inceptisols and Andisols in the west had a comparable magnitude of topsoil SOC loss (–0.5 to –1.8 Mg C ha –1 year –1 ); (b) Mollisols and Alfisols in the Central Plains sites were susceptible to significant SOC loss under farming and grazing; (c) Runoff/erosion and leaching potential, vertical translocation, and mineral sorption were the most important factors controlling SOC variation across the NEON sites. Our work could be used to parameterize ecosystem models simulating SOC change.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

The Role of Non-Local Effects on Surface Sensible Heat Flux Under Different Types of Thermal Structures Over the Arctic Sea-Ice Surface

The effects of atmospheric thermal structure on the surface energy flux are poorly understood over the Arctic sea-ice surface. Here, we explore the mechanism of sensible heat exchange under different types of thermal structures over the Arctic sea-ice surface by using data collected during the Multidisciplinary drifting Observatory for the Study of Arctic Climate expedition. The quadrant analysis indicates that strong surface temperature inversions below 100 m enhance non-local effects on the positive (upward) sensible heat flux ($\overline{w'θ'}$) through entrainment of large eddies from the convective boundary layer aloft. However, strong surface inversions restrict the contributions of large eddies to the negative (downward) $\overline{w'θ'}$ due to intensified surface stability. By inspecting the existing parameterization schemes, we found that the European Center for Medium-Range Weather Forecasts Integrated Forecasting System scheme fails to predict the impacts of non-local processes on the positive $\overline{w'θ'}$, and an adjustment term to correct the bias of parameterized $\overline{w'θ'}$ is proposed.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Modeling the Winter Heat Conduction Through the Sea Ice System During MOSAiC

Abstract Models struggle to accurately simulate observed sea ice thickness changes, which could be partially due to inadequate representation of thermodynamic processes. We analyzed co‐located winter observations of the Arctic sea ice from the Multidisciplinary Drifting Observatory for the Study of the Arctic Climate for evaluating and improving thermodynamic processes in sea ice models, aiming to enable more accurate predictions of the warming climate system. We model the sea ice and snow heat conduction for observed transects forced by realistic boundary conditions to understand the impact of the non‐resolved meter‐scale snow and sea ice thickness heterogeneity on horizontal heat conduction. Neglecting horizontal processes causes underestimating the conductive heat flux of 10% or more. Furthermore, comparing model results to independent temperature observations reveals a ∼5 K surface temperature overestimation over ice thinner than 1 m, attributed to shortcomings in parameterizing surface turbulent and radiative fluxes rather than the conduction. Assessing the model deficiencies and parameterizing these unresolved processes is required for improved sea ice representation.

Geology↗
DOE OSTI2024

Weathering Incongruence in Mountainous Mediterranean Climates Recorded by Stream Lithium Isotope Ratios

Lithium isotope ratios (δ 7 Li) of rivers are increasingly serving as a diagnostic of the balance between chemical and physical weathering contributions to overall landscape denudation rates. Here, we show that intermediate weathering intensities and highly enriched stream δ 7 Li values typically associated with lowland floodplains can also describe small upland watersheds subject to cool, wet climates. This behavior is revealed by stream δ 7 Li between +22.4 and +23.5‰ within a Critical Zone observatory located in the Cévennes region of southern France, where dilute stream solute concentrations and significant atmospheric deposition otherwise mask evidence of incongruence. The water-rock reaction pathways underlying this behavior are quantified through a multicomponent, isotope-enabled reactive transport model. Using geochemical characterization of soil profiles, bedrock, and long-term stream samples as constraints, we evolve the simulation from an initially unweathered granite to a steady state weathering profile which reflects the balance between (a) fluid infiltration and drainage and (b) bedrock uplift and soil erosion. Enriched stream δ 7 Li occurs because Li is strongly incorporated into actively precipitating secondary clay phases beyond what prior laboratory experiments have suggested. Chemical weathering incongruence is maintained despite relatively slow reaction rates and moderate clay accumulation due to a combination of two factors. First, reactive primary mineral phases persist across the weathering profile and effectively “shield” the secondary clays from resolubilization due to their greater solubility. Second, the clays accumulating in the near-surface profile are relatively mature weathering byproducts. These factors promote characteristically low total dissolved solute export from the catchment despite significant input of exogenous dust.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Upscaling Soil Organic Carbon Measurements at the Continental Scale Using Multivariate Clustering Analysis and Machine Learning

Abstract Estimates of soil organic carbon (SOC) stocks are essential for many environmental applications. However, significant inconsistencies exist in SOC stock estimates for the U.S. across current SOC maps. We propose a framework that combines unsupervised multivariate geographic clustering (MGC) and supervised Random Forests regression, improving SOC maps by capturing heterogeneous relationships with SOC drivers. We first used MGC to divide the U.S. into 20 SOC regions based on the similarity of covariates (soil biogeochemical, bioclimatic, biological, and physiographic variables). Subsequently, separate Random Forests models were trained for each SOC region, utilizing environmental covariates and SOC observations. Our estimated SOC stocks for the U.S. (52.6 ± 3.2 Pg for 0–30 cm and 108.3 ± 8.2 Pg for 0–100 cm depth) were within the range estimated by existing products like Harmonized World Soil Database, HWSD (46.7 Pg for 0–30 cm and 90.7 Pg for 0–100 cm depth) and SoilGrids 2.0 (45.7 Pg for 0–30 cm and 133.0 Pg for 0–100 cm depth). However, independent validation with soil profile data from the National Ecological Observatory Network showed that our approach ( R 2 = 0.51) outperformed the estimates obtained from Harmonized World Soil Database ( R 2 = 0.23) and SoilGrids 2.0 ( R 2 = 0.39) for the topsoil (0–30 cm). Uncertainty analysis (e.g., low representativeness and high coefficients of variation) identified regions requiring more measurements, such as Alaska and the deserts of the U.S. Southwest. Our approach effectively captures the heterogeneous relationships between widely available predictors and the current SOC baseline across regions, offering reliable SOC estimates at 1 km resolution for benchmarking Earth system models.

58 GEOSCIENCES↗
DOE OSTI2024

A Stereo Camera Simulator for Large-Eddy Simulations of Continental Shallow Cumulus Clouds Based on Three-Dimensional Path-Tracing

The complex spatial and temporal structure of cumulus clouds complicates their representation in weather and climate models. Classic meteorological instrumentation struggles to fully capture these features. Networks of multiple high-resolution hemispheric cameras are increasingly used to fill this data gap, and provide information on this missing multi-dimensional spatial information. In this study, a path-tracing algorithm is used to generate virtual camera images of resolved clouds in large-eddy simulations (LES). These images are then used as a camera network simulator, allowing reconstructions of three-dimensional cloud edges from the model output. Because the actual LES cloud field is fully known, the combined path-tracing and reconstruction method can be statistically analyzed. The method is applied to LES realizations of summertime shallow cumulus at the Jülich Observatory for Cloud Evolution (JOYCE), Germany, which also routinely operates a camera network. We find that the path-tracing method allows accurate reconstruction of up to 70% of the visible cloud edges. Additional sensitivity tests show that the method is robust for changes in its hyperparameters. The sensitivity to cloud optical thickness is also investigated, finding a cloud boundary placement error of approximately 182 m. This error can be considered typical for cloud boundary reconstruction using real stereo camera imagery. The results provide proof of principle for future use of the method for evaluating LES clouds against camera network imagery, and for further optimizing the configuration of such camera networks.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

In Situ Velocity‐Strain Sensitivity Near the San Jacinto Fault Zone Analyzed Through Train Tremors

Abstract We utilize train tremors as P‐wave seismic sources to investigate velocity‐strain sensitivity near the San Jacinto Fault Zone. A dense nodal array deployed at the Piñon Flat Observatory is used to detect and identify repeating train energy emitted from a railway in the Coachella valley. We construct P‐wave correlation functions across the fault zone and estimate the spatially averaged dt / t versus strain sensitivity to be 6.25 × 10 4 . Through numerical simulations, we explore how the sensitivity decays exponentially with depth. The optimal solution reveals a subsurface sensitivity of 1.2 × 10 5 and a depth decay rate of 0.05 km −1 . This sensitivity aligns with previous findings but is toward the higher end, likely due to the fractured fault‐zone rocks. The depth decay rate, previously unreported, is notably smaller than assumed in empirical models. This raises the necessity of further investigations of this parameter, which is crucial to study stress and velocity variations at seismogenic depth.

Geology↗
DOE OSTI2024

The Vertical Structure of Turbulence Kinetic Energy Near the Arctic Sea‐Ice Surface

Atmospheric turbulence over the Arctic sea-ice surface has been understudied due to the lack of observational data. In this study, we focus on the turbulence kinetic energy (TKE) over sea ice and distinguish its two different vertical structures, the “Surface” type and the “Elevated” type, using observations during the Multidisciplinary drifting Observatory for the Study of Arctic Climate expedition (MOSAiC). The “Surface” type has the maximum TKE near the surface (at 2 m), while the “Elevated” type has the maximum TKE at a higher level (6 m). The TKE budget analysis indicates that the “Elevated” type is caused by the increased shear production of TKE at 6 m. In addition, spectral analysis reveals that the contribution to TKE by horizontal large eddies is enhanced in the “Elevated” type. Finally, how the vertical structure of TKE affects the parameterization of turbulent momentum flux is discussed.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Estimating Carbon Dioxide Emissions in Two California Cities Using Bayesian Inversion and Satellite Measurements

NASA's Orbiting Carbon Observatories (OCO-2 and OCO-3) provide measurements of column-averaged carbon dioxide concentrations (XCO 2 ) with sufficient spatial resolution and precision to constrain bottom-up estimates of CO 2 fluxes at regional scales. We use Bayesian inversion methods assimilating satellite retrievals to improve estimates of CO 2 fluxes in the South Coast Air Basin (SoCAB) which surrounds Los Angeles, and in the San Francisco Bay Area Air Basin (SFBA). We study 2020 to understand the impact of the COVID-19 lockdowns and an active wildfire season. Our results indicated that a 50% (30%) reduction in CO 2 emissions relative to 2015 during the COVID-19 lockdown period was consistent with OCO measurements for SFBA (SoCAB). We find that posterior wildfire emissions differed significantly from the prior at the scale of individual wildfires, though with large uncertainties, and that wildfire emissions in SFBA are significant, attributing 72% of the region's CO 2 emissions during August 2020 to wildfires.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Evaluation of the Planetary Boundary Layer Height From ERA5 Reanalysis With MOSAiC Observations Over the Arctic Ocean

The planetary boundary layer height (PBLH) is a crucial indicator reflecting the region of the atmosphere characterized by continuous turbulence. Here, we use radiosonde and surface meteorological observations (4–7 times per day, year-round measurements) during the Multidisciplinary drifting Observatory for the Study of Arctic Climate (MOSAiC) expedition to derive the PBLH (PBLH MOSAiC ), and further evaluate the PBLH from the ERA5 reanalysis (PBLH ERA5 ). Comparisons between PBLH MOSAiC and PBLH ERA5 from different perspectives reveal that: (a) The overestimation of PBLH ERA5 when the sea ice concentration is >90% is significant with the centered root mean squared error reaching up to 201 m; (b) The difference between the two products is notably pronounced in cold seasons, while it is comparatively diminished in warm seasons; (c) In neutral boundary layers, differences in PBLH ERA5 are larger compared with stable and convective boundary layers. In addition, the analysis of error sources indicates that the bias of PBLH ERA5 is sensitive to the bias of vertical thermal structure and wind speed profiles in ERA5 data sets in all conditions. Finally, we find a Random Forest model effectively reduces the bias of PBLH ERA5 with the index of agreement reaching up to 0.71 in the test data set, while a multiple linear regression demonstrates comparable performance to the Random Forest model.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Air‐Ice‐Ocean Coupling During a Strong Mid‐Winter Cyclone: Observing Coupled Dynamic Interactions Across Scales

Abstract Arctic cyclones are key drivers of sea ice and ocean variability. During the 2019–2020 Multidisciplinary drifting Observatory for the Study of Arctic Climate (MOSAiC) expedition, joint observations of the coupled air‐ice‐ocean system were collected at multiple spatial scales. Here, we present observations of a strong mid‐winter cyclone that impacted the MOSAiC site as it drifted in the central Arctic pack ice. The sea ice dynamical response showed spatial structure at the scale of the evolving and translating cyclonic wind field. Internal ice stress and ocean stress play significant roles, resulting in timing offsets between the atmospheric forcing and the ice response and post‐cyclone inertial ringing in the ice and ocean. Ice motion in response to the wind field then forces the upper ocean currents through frictional drag. The strongest impacts to the sea ice and ocean from the passing cyclone occur as a result of the surface impacts of a strong atmospheric low‐level jet (LLJ) behind the trailing cold front and changing wind directions between the warm‐sector LLJ and post cold‐frontal LLJ. Impacts of the cyclone are prolonged through the coupled ice‐ocean inertial response. Local impacts of the approximately 120 km wide LLJ occur over a 12 hr period or less and at scales of a kilometer to a few tens of kilometers, meaning that these impacts occur at combined smaller spatial scales and faster time scales than most satellite observations and coupled Earth system models can resolve.

Meteorology & Atmospheric Sciences↗
DOE OSTI2024

Use of Satellite, Surface Observations and Numerical Weather Prediction Model Data to Improve Cloud Base Height and Cloud Base Vertical Velocity Estimation

Cloud base height (CBH) and cloud base vertical velocity (CBVV) are important variables that impact the overall climate in a region as they influence the formulation, longevity, and evolution of clouds. Retrieval of both parameters have long used ground instrumentation (e.g., Doppler lidar (DL), ground base radar); however, retrieving CBH from satellites is particularly challenging given that space-based instruments only observe cloud tops. In this manuscript, CBH is retrieved using a multi-linear regression equation, while CBVV used a random forests model. Both retrievals combine satellite and numerical weather prediction data. The satellite data used are the Visible Infrared Imaging Radiometer Suite imagery, while measurements of CBH and CBVV include DL and radiosonde data at the Southern Great Plains (SGP) Atmospheric Radiation Measurement observatory. Data from 83 summer days (May-August) in 2018–2021 featuring cumulus clouds forced by solar heating were examined and used to train the models, with years 2022–2023 used for validation. Various spatial domains were defined with one large (2.4° longitude by 2.0° latitude) SGP domain being split into smaller sections (smallest being 0.99° and 0.61° longitude and latitude respectably). CBH and CBVV values obtained from the DL as compared to the models show root mean square errors between 150 and 200 m, with CBVV values between 0.45 and 1 ms -1 . Finally, it was found that the CBH formulation performs well over all domains, while the CBVV retrievals become less accurate due to more turbulence being introduced into the observations as the number of DL stations decreases in the smaller domains.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Sequential Precipitation Input Tagging (SPIT) to Estimate Water Transit Times and Hydrologic Tracer Dynamics Within Water‐Tagging Enabled Hydrologic Models

Determining the age distribution of water exiting a catchment is important for understanding groundwater storage and mixing. New water-tagging capabilities within models track precipitation events as they move through simulated storages, yet forward modeling of individual events may not systematically capture the full transit time distribution (TTD). Here, we present a “sequential precipitation input tagging” (SPIT) framework to tag all input precipitation at regular intervals during extended model simulations. Monthly tags over 7 years were applied at six National Ecological Observatory Network sites to calculate TTDs and derive mean virtual tracer age, $\overline{T_{V}}$, fractions of young water, F yw , and hydrologic tracer concentrations (water isotopes δ 18 O and δ 2 H) within a tagging enabled version of the Weather Research and Forecast hydrologic model (WRF-Hydro). Throughout seven simulation years, the fraction of simulated discharge derived from tagged events, F tag , increased each year, with the final year's F tag ranging from 66% to 100% and highlights the need to apply SPIT over many years to understand TTDs. When the F tag was >75%, simulated $\overline{T_{V}}$ ranged 179–923 days and F yw 0.6%–23.9%, with daily values exhibiting a power-law relationship with precipitation, discharge, and groundwater. Through implementation of SPIT, we find this hydrologic model configuration performs poorly in estimation of $\overline{T_{V}}$ and F yw (root mean squared error of 469 days and 14.4% respectively), suggesting it misrepresents subsurface mixing. Thus, the SPIT framework provides a reproducible approach to calculate watershed transit times within tagging enabled models and thereby assess and improve representation of hydrologic processes.

fraction of young water↗