Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Reproducibility”

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 685 records · Page 38
All sourcesNASA NTRSDOE OSTI
DOE OSTI2019

Diverse Mycorrhizal Associations Enhance Terrestrial C Storage in a Global Model

Accurate projections of the terrestrial carbon (C) sink are critical to understanding the future global C cycle and setting CO 2 emission reduction goals. Current earth system models (ESMs) and dynamic global vegetation models (DGVMs) with coupled carbon-nitrogen cycles project that future terrestrial C sequestration will be limited by nitrogen (N) availability, but the magnitude of N limitation remains a critical uncertainty. Plants use multiple symbiotic nutrient acquisition strategies to mitigate N limitation, but current DGVMs omit these mechanisms. Fully coupling N-acquiring plant-microbe symbioses to soil organic matter (SOM) cycling within a DGVM for the first time, we demonstrate that increases in N acquisition via SOM decomposition and atmospheric N 2 fixation could support long-term enhancement of terrestrial C sequestration at global scales under elevated CO 2 . The model reproduced elevated CO2 responses from two experiments (Duke and Oak Ridge) representing contrasting N acquisition strategies. N release from enhanced SOM decomposition supported vegetation growth at Duke, while inorganic N depletion limited growth at Oak Ridge. Global simulations reproduced spatial patterns of N-acquiring symbioses from a novel niche-based map of mycorrhizal fungi. Under a 100-ppm increase in CO2 concentrations, shifts in N acquisition pathways facilitated 200 Pg C of terrestrial C sequestration over 100 years compared to 50 Pg C for a scenario with static N acquisition pathways. Our results imply that N acquisition strategies are important determinants of terrestrial C sequestration potential under elevated CO 2 and that nitrogen-enabled DGVMs that omit symbiotic N acquisition may underestimate future terrestrial C uptake.

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

Simulations of Van Allen Probes Plasmaspheric Electron Density Observations

Here we simulate equatorial plasmaspheric electron densities using a physics-based model (Cold PLasma, CPL; used in the ring current-atmosphere interactions model) of the source and loss processes of refilling and erosion driven by empirical inputs. The performance of CPL is evaluated against in situ measurements by the Van Allen Probes (Radiation Belt Storm Probes) for two events: the 31 May to 5 June and 15 to 20 January 2013 geomagnetic storms observed in the premidnight and postmidnight magnetic local time (MLT) sectors, respectively. Overall, CPL reproduces the radial extent of the plasmasphere to within a mean absolute difference of |ΔL avg |=0.5 ± 0.1 L. The model electric field responsible for E × B convection and the parameterization of geomagnetic conditions (under the Kp-index and solar wind properties) implemented by CPL did not account for localized enhancements in the duskward electric field during increased activity. Rather, it was found to be largely dependent on the measure of the quiet time background. Finally, this property indicates that the agreement between these simulations and observations does not account for the complete set of physical processes during extreme (strong or weak) geomagnetic conditions impacting the plasmasphere. Nevertheless, at the presented resolution of the model CPL does provide good agreement in reproducing Radiation Belt Storm Probes observations of plasmaspheric density and plasmapause location.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2019

Northern Hemisphere Blocking in ~25-km-Resolution E3SM v0.3 Atmosphere-Land Simulations

The fidelity of a prerelease Energy Exascale Earth System Model (E3SM) in reproducing atmospheric blocking is evaluated with ~100- and ~25-km horizontal resolution within ensembles of simulations with Atmospheric Modeling Intercomparison Protocol-type active atmosphere and land surface configurations. Here this evaluation is conducted via a hybrid bidimensional blocking index based on geopotential height. The lower-resolution model correctly reproduces the spatial distribution of blocking frequency maxima over the Northern Hemisphere. However, it overestimates the blocking frequency over western North America, the Pacific Northeast, and Eastern Europe regions and largely underestimates blocking frequency over the North Atlantic. The high-resolution model significantly reduces the bias over the North Pacific region, particularly over western North America, but the biases over the North Atlantic blocking sector mostly persist. A diagnosis of the mean flow reveals that the subtropical jet is displaced poleward in the high-resolution model with a jet core speed that is generally more realistic. A time scale analysis of eddies suggests that synoptic eddies associated with the blocking events are increased in the high-resolution configuration. The discrepancy between the models implies that both synoptic and low-frequency eddies play significant roles in determining North Pacific blocking, whereas North Atlantic blocking is mainly driven by low-frequency eddies. The characteristic time scale associated with blocking events is also investigated, and it shows that both model resolutions produce events that are overly persistent. However, the bias in the high-resolution model is reduced, which reflects more efficient energy dispersion with higher resolution.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2019

Using Information Theory to Evaluate Directional Precipitation Interactions Over the West Sahel Region in Observations and Models

Water availability has historically been one of the most significant threats to African regional social and economic well-being. Over the Sahel region, a megadrought during the 1960s and 1970s induced by an abrupt and substantial rainfall reduction caused widespread famine and death. The postdrought recovery, which is still ongoing, has been characterized by gradual increases in rainfall, but with dramatic fluctuations. The large negative human impacts, slow recovery, and variability raise important questions of why and how rainfall dynamics evolve and interact with other components of the regional climate system. Here we provide an observational assessment of rainfall interactions mechanisms (informed by directional transfer of information entropy) that regulate Sahel rainfall. We quantitatively demonstrate that (1) sea surface temperature over the Gulf of Guinea controls moisture advection and transport to the West Sahel region and (2) strong bidirectional interactions exist between local vegetation dynamics and rainfall patterns. Then we assess the directional interaction patterns from nine state-of-the-art Earth System Models (ESMs). We find that most ESMs are able to represent either the unidirectional control of sea surface temperature on precipitation or the bidirectional interaction between vegetation and precipitation. However, none of the ESMs represents both interactive patterns. The GFDL and IPSL-CM5A-LR models successfully reproduced observed patterns over ~50% of the West Sahel region but were not accurate in reproducing observed precipitation regional trends or interannual variation. We propose that the directional information transfer is a powerful mechanistic benchmark to assess model fidelity at the process level.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2019

Projected Changes in United States Regional Extreme Heat Days Derived From Bivariate Quantile Mapping of CMIP5 Simulations

Projections of multivariate climate extremes require methodological approaches that can maintain relationships among variables. Here we apply a piecewise multivariate quantile mapping approach to temperature and humidity projections from a subset of models from the Fifth Coupled Model Intercomparison Project and analyze the resulting climatology of extreme heat days (EHDs) with explicit consideration of the prevailing humidity. The piecewise multivariate bias correction method shows good fidelity in reproducing the frequency of different types of extreme heat events in the historical period. Projections for U.S. regions and individual cities for both the middle and end of the century are characterized by increases in the frequency of EHDs, and especially those characterized by high humidity. For many regions and individual cities, there is no overlap between the frequency of high-humidity EHDs in general circulation model ensembles from the historical and future periods, indicating that increases in extreme heat are robust. Analysis of 500-mb height, sea-level pressure, and low-level circulation composites for historical and future periods indicates that the Fifth Coupled Model Intercomparison Project models reproduce basic large-scale circulation features associated with EHDs in U.S. regions and that future changes in extreme heat are related primarily to large-scale warming rather than enhancement of regional circulation anomalies.

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

Modeling Transient Soil Moisture Limitations on Microbial Carbon Respiration

Soil microorganisms are known to survive periods of aridity and to recover rapidly after wetting events, with the ability to transition between a dormant state in dry conditions and an active state in wet conditions. While this dynamic behavior has been previously incorporated into soil carbon respiration modeling frameworks, a direct comparison between this active-dormant transition mechanism and a more simplified first-order model has yet to be made. Here, we demonstrate the necessary extent of model complexity needed to reproduce transient carbon respiration rates obtained from a set of soil incubation experiments implemented over a range of soil depths and time intervals. Two approaches are tested, one uses simplified first-order kinetics, whereas the other employs a transition between active and dormant biomass. The performance of each model is evaluated using an Akaike Information Criterion (AIC) based on the accuracy with which they reproduce an experimental dataset consisting of two sets of time series soil incubations collected across a range of time and depth resolutions. Based on the AIC evaluation and model-data comparison, we conclude that a dormancy-enabled model featuring two distinct microbial strategists performs best for the majority of the soil profile (above 108 cm) for both high and low depth resolution and sampling frequency, despite the added parameters required. In contrast, the first-order model achieves better AIC scores when simulating our deepest soils (112–165 cm), where moisture fluctuations are expected to be less prevalent. These findings guide how and where we choose to apply more cost intensive models.

58 GEOSCIENCES↗
DOE OSTI2019

Evaluation of Clouds in Version 1 of E3SM Atmosphere Model with Satellite and Simulators

This study systematically evaluates clouds simulated by the Energy Exascale Earth System Model (E3SM) Atmosphere Model version one (EAMv1) against satellite and ground-based cloud observations. Both low (1°) and high (0.25°) resolution EAMv1 configurations generally underestimate clouds in low and midlatitutes and overestimate clouds in the Arctic although the error is smaller in the high-resolution model. The underestimate of clouds is due to the underestimate of optically thin to intermediate clouds. EAMv1 overestimates the optically intermediate to thick clouds. Other model errors include the largely under-predicted stratocumulus along the coasts and high clouds over the tropical deep convection regions. The underestimate of thin clouds results in too much LW radiation being emitted to space and too little SW radiation being reflected back to space while the overestimate of optically intermediate and thick clouds leads to too little LW radiation being emitted to space and too much SW radiation being reflected back to space. EAMv1 shows better skill in reproducing the observed distribution of clouds and their properties and has smaller radiatively relevant errors in the distribution of clouds than most of the CFMIP2 models. It produces more supercooled liquid cloud fraction than CAM5 and most CMIP5 models owing to a new ice nucleation scheme and a reduction of ice deposition growth rate. It simulates the diurnal variation of clouds during the summer season at SGP qualitatively well, with both the evolution of shallow clouds and the nocturnal peak in high clouds well captured, however, it largely overestimates the observed magnitude.

Zhang, Yuying↗
DOE OSTI2019

Modeling the Magnetopause Shadowing Loss During the June 2015 Dropout Event

Fast dropout of relativistic and ultrarelativistic electrons at both high and low L* regions were observed during the intense coronal mass ejection driven storm in June 2015. An improved radial diffusion model, using an event-specific Last Closed Drift Shell (LCDS) and newly-available radial diffusion coefficients ( DD LL ), is implemented to simulate the magnetopause shadowing loss of electrons. The model captures the fast shadowing loss of electrons well at high L* regions after both interplanetary shocks, and reproduces the initial adiabatic loss of the high-energy storage ring at low L* regions after the second strong shock. In this paper, we show for the first time that using the event-specific and K-dependent LCDS and improved DD LL is critical to reproduce the observed dropout features, including the timing, location, and the butterfly electron pitch angle distribution. Future inclusion of the EMIC wave scattering process is needed to model the observed further depletion of the storage ring.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2019

Assessing XMT-Measurement Variability of Air-Water Interfacial Areas in Natural Porous Media

This study investigates the accuracy and reproducibility of air-water interfacial areas measured with high-resolution synchrotron X-ray microtomography (XMT). Columns packed with one of two relatively coarse-grained monodisperse granular media, glass beads or a well-sorted quartz sand, were imaged over several years, encompassing changes in acquisition equipment, improved image quality, and enhancements to image acquisition and to processing software. For the glass beads, the specific solid surface area (SSSA-XMT) of 31.6 ± 1 cm -1 determined from direct analysis of the segmented solid-phase image data is statistically identical to the independently calculated geometric smooth-sphere specific solid surface area (GSSA, 32 ± 1 cm -1 ) and to the measured SSSA (28 ± 3 cm -1 ) obtained with the N2-Brunauer, Emmett, and Teller method. The maximum specific air-water interfacial area (Amax) is 27.4 (±2) cm -1 , which compares very well to the SSSA-XMT, GSSA, and SSSA-N2-Brunauer, Emmett, and Teller values. For the sand, the SSSA-XMT (111 ± 2 cm -1 ) and GSSA (113 ± 1 cm -1 ) are similar. The mean Amax is 96 ± 5 cm -1 , which compares well to both the SSSA and the GSSA values. The XMT-SSSA values deviated from the GSSA values by 7–16% for the first four experiments but were essentially identical for the later experiments. This indicates that enhancements in image acquisition and processing improved data accuracy. The Amax values ranged from 74 cm -1 to 101 cm -1 , with a coefficient of variation (COV) of 9%. The maximum capillary interfacial area ranged from 12 cm-1 to 19 cm -1 , for a COV of 10%. The COVs for both decreased to 5–6% for the latter five experiments. These results demonstrate that XMT imaging provides accurate and reproducible measurements of total and capillary interfacial areas.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2019

Mimicking Geologic Depositional Fabrics for Multiphase Flow Experiments

A novel laboratory technique is presented for generating two-dimensional beadpacks with reproducible natural geologic sedimentary bedforms. The heterogeneous beadpacks are deposited in a 0.6 m × 0.6 m × 0.02 m glass slab chamber using a programmable two-axis linear actuator. The method takes advantage of natural mechanical sorting processes, along with the ability to program the motion of the arm, to efficiently build cross-beds and ripple laminae similar to those seen in natural water-lain clastic deposits. By using a mixture of grain sizes (200 to 800 μm diameter), depositional fabrics with multiscale heterogeneity (mm to dm) are generated. Reproducibility tests show very high correlation among multiple fabric realizations produced with similar inputs. After introducing this first-of-a-kind capability to produce depositionally realistic fabrics at scales larger than that of rock cores, a demonstrative study is presented that illustrates gravity unstable multiphase flow which is relevant in the context of CO 2 sequestration. Preliminary results highlight the dramatic impact small scale capillary heterogeneities have on multiphase immiscible displacement and can help bridge understanding between core scale experiments and reservoir scale observations.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Reconstruction of Electron and Ion Distribution Functions in a Magnetotail Reconnection Diffusion Region

In the diffusion region of magnetotail reconnection, particle distributions are highly structured, exhibiting triangular shapes and multiple striations that deviate dramatically from the Maxwellian distribution. Fully kinetic simulations have been demonstrated to be capable of producing the essential structures of the observed distribution functions, yet are computationally not feasible for 3D global simulations. The fluid models used for large-scale simulations, on the other hand, do not have the kinetic physics necessary for describing reconnection accurately. Our study aims to bridge fully kinetic and fluid simulations by quantifying the information required to capture the non-Maxwellian features in the distributions underlying the closures used in the fluid code. Here, we compare the results of fully kinetic simulations with observed electron velocity distributions in a magnetotail reconnection diffusion region and use the maximum entropy model to reconstruct electron and ion distributions using various numbers of moments obtained from the simulation. Our results indicate that using only local moments, the maximum entropy model can reproduce many of the features of the distributions: (1) the electron outflow distribution with a tilted triangular structure is reproduced with 21 or more moments in agreement with Ng et al. and (2) counterstreaming distributions can be captured with the 35-moment model when the separation in velocity space between the populations is large.

58 GEOSCIENCES↗
DOE OSTI2021

A New Approach for Simultaneous Estimation of Entrainment and Detrainment Rates in Non- Precipitating Shallow Cumulus

A new approach is developed for estimating entrainment and detrainment rates in cumulus clouds based on aircraft observations. Equations relating entrainment and detrainment rates to gross entrainment and detrainment are derived. This approach is applied to the Holistic Interactions of Shallow Clouds, Aerosols, and Land-Ecosystems field campaign, supported by the U.S. Department of Energy's Atmospheric Radiation Measurement program. The results show that both entrainment and detrainment rates decrease with increasing height. Sensitivity tests with different detrained air assumptions yield similar results. The entrainment and detrainment rates can reproduce the cloud thermodynamic variables. Partial correlation analysis indicates that entrainment rate is positively correlated with environmental relative humidity (RH), and detrainment rate is negatively correlated with environmental RH and positively correlated with entrainment rate. This new approach can be applied to other cloud observations to obtain a data set of entrainment and detrainment rates in cumulus clouds.

Zhu, Lei↗
DOE OSTI2022

Annual Mean Arctic Amplification 1970–2020: Observed and Simulated by CMIP6 Climate Models

While the annual mean Arctic Amplification (AA) index varied between two and three during the 1970–2000 period, it reached values exceeding four during the first two decades of the 21st century. The AA did not change in a continuous fashion but rather in two sharp increases around 1986 and 1999. During those steps the mean global surface air temperature trend remained almost constant, while the Arctic trend increased. Although the “best” CMIP6 models reproduce the increasing trend of the AA in 1980s they do not capture the sharply increasing trend of the AA after 1999 including its rapid step-like increase. We propose that the first sharp AA increase around 1986 is due to external forcing, while the second step close to 1999 is due to internal climate variability, which models cannot reproduce in the observed time.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

WUS256: An Adjoint Waveform Tomography Model of the Crust and Upper Mantle of the Western United States for Improved Waveform Simulations

Abstract We report a new model (WUS256) of radially anisotropic seismic wavespeeds of the crust and upper mantle of the western United States (WUS) obtained from adjoint waveform tomography for the purpose of improving synthetic waveform fits to observed data. WUS256 is based on inversion of over 94,000 waveforms from 72 earthquakes recorded by nearly 3,400 stations. We started with the SPiRaL global model (Simmons et al., 2021, https://doi.org/10.1093/gji/ggab277 ) and waveforms in the period band of 50–120 s. We followed a conservative multiscale inversion approach with eight stages and 256 total inversion iterations which enabled monotonic misfit reduction to 20‐s minimum‐period waves. WUS256 relied on time‐frequency (TF) phase misfits and a trust region limited memory Broyden–Fletcher–Goldfarb–Shanno (L‐BFGS) optimization. Hessian‐vector products were used to qualitatively assess model resolution. Results indicate that WUS256 has good coverage of the continental regions to depths of about 150 km and is able to resolve features on lateral scales of about 200 km. We quantify waveform fits by the reduction in TF and normalized amplitude difference misfits between WUS256 and the SPiRaL starting model. WUS256 significantly improves waveform fits with misfit reduction 64% for both inversion and validation data sets compared to the SPiRaL starting model and shows even better fits compared to other models. Waveform fits illustrate that WUS256 reproduces body‐waves, fundamental mode surface waves as well as late arriving dispersed and/or scattered short period surface waves. The improvement in waveform fit indicates that WUS256 can be used to reproduce path effects on regional complete waveforms and moment tensor inversions.

58 GEOSCIENCES↗
DOE OSTI2024

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↗
DOE OSTI2025

Evaluation of Marine Boundary Layer Clouds Over the Northeast Pacific During the CSET Campaign in E3SM Version 2

It is still challenging to reproduce marine boundary layer (MBL) clouds well in large-scale models despite their importance to the Earth's radiation budget and hydrological cycle. This study evaluates representation of the MBL and clouds in the Energy Exascale Earth System Model (E3SM) version 2. This study compares the E3SM simulation results with remote sensing and reanalysis data during the Cloud System Evolution in the Trades (CSET) field campaign to better understand the stratocumulus-to-cumulus cloud transition (SCT) over the northeast Pacific. E3SM results are extracted along the CSET Lagrangian trajectories. The comparison shows that the E3SM simulation applying horizontal wind nudging performs well in reproducing thermodynamic variables of the MBL and evolution trends of cloud variables along the trajectories. However, substantial overestimations of aerosol and cloud drop number (N d ) are observed, which is explained as an issue with version 2 of the model. Cloud fraction (CF) does decrease from the Californian coast to Hawaii in the E3SM simulation, but most CF values indicate either an overcast or almost clear sky, which differs from satellite and reanalysis data. The effect of N d overestimation on CF evolution is assessed via prescribed-N d simulations. Those simulations with N d modifications show negligible CF changes. A comparison of estimated inversion strength (EIS) also shows that the simulated EIS values are similar to those of reanalysis data. Our study suggests that cloud macrophysics and boundary layer processes are more important in improving the simulation to capture the SCT than refining the model's thermodynamics or cloud microphysics.

Choi, Kyoung Ock [Univ. of Washington, Seattle, WA↗
DOE OSTI2026

Modeling Large Dust Aerosols in the Community Earth System Model Version 2 (CESM2)

Dust aerosols have a wide size distribution from less than 0.1 to over 100 μm and dominate Earth's atmospheric aerosol mass. However, most Earth system models (ESMs) inadequately represent dust aerosols larger than 10 μm in diameter, limiting the accuracy of the simulated dust cycle and climate impacts. Here, we introduce a new modeling framework that captures the full observed size distribution of dust aerosols, incorporating recent advances into a mineral-resolved version of the Community ESM, while addressing known issues in previous versions. Comprehensive evaluation against diverse observations of bulk dust and component minerals demonstrates that the model reproduces the observed dust cycle across particle sizes. Incorporating the previously unrepresented large-dust fractions substantially alters dust budget estimates, highlighting potential changes in simulated climate impacts and underscoring the importance of comprehensive size-resolved dust modeling. Despite these advancements, uncertainties persist. Our results indicate that a size-dependent reduction in settling velocity is required to reproduce the observed dust size distribution downwind of source regions. Specifically, in the new model, the gravitational settling velocity of dust particles larger than 10 μm in diameter must be reduced by as much as 85% to achieve agreement with observations. This empirical reduction serves as a constraint on physics-based models of dust settling. Future developments should address misrepresented physical processes that hinder accurate modeling of the large dust aerosol transport. Expanding observational data sets covering the full-size distribution is also essential to better constrain the dust cycle and improve the representation of dust optical properties and climate effects.

Li, Longlei [Cornell Univ., Ithaca, NY (United Sta↗
DOE OSTI2026

An Integrated Modeling Framework for Sediment Dynamics During Urban Flooding: Application to Hurricane Harvey in Houston

Floodwater can mobilize and redistribute large volumes of sediment from upland to downstream urban areas, threatening infrastructure, water quality, and ecosystem health. However, existing modeling approaches often fail to capture sediment dynamics in urban floodplains due to the lack of integration between upland hydrological processes and riverine sediment transport. This study presents the first integrated modeling framework that couples the Energy Exascale Earth System Model (E3SM) land component, which simulates runoff and hillslope erosion, with TELEMAC-GAIA, a two-dimensional hydrodynamic and sediment transport model. This framework enables the fully distributed, process-based simulation of high-resolution (as fine as 30 m) sediment dynamics from hillslopes to floodplains. Applied to a highly urbanized watershed in Houston during Hurricane Harvey, this framework reproduced observed water levels at 16 USGS gauges (median R 2 = 0.83 and KGE = 0.78), key sediment dynamics such as sediment transport and deposition processes, and reproduced spatial deposition patterns consistent with LiDAR-derived data. Based on the simulation, we estimate 8.0 million m 3 of event-scale sediment deposition, including 5.7 million m 3 trapped in the flood-control reservoirs and 2.3 million m 3 deposited along major channels and floodplains. Using a representative unit removal cost, this corresponds to an estimated dredging cost of $581 million for total deposition. These results provide a first-order, physically based quantification of Harvey-scale sediment impacts. This study provides a valuable tool for the holistic analysis of sediment dynamics triggered by extreme urban flooding, supporting flood-resilience planning. More broadly, it highlights the importance of integrating physically based hydrological processes for urban flooding and sediment research.

Hurricane Harvey↗