Search NASA⌕ Search

SEARCH · Search NASA

Results for “process simulation modeling”

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 1,549 records · Page 86
All sourcesNASA NTRSDOE OSTI
DOE OSTI2019

A robust process-structure model for predicting the joint interface structure in impact welding

The impact welding is gaining popularity in joining of dissimilar metals as any intermetallics are often very thin or discontinous with little to no heat affected zone resulting from the process. This work focuses on a recently developed impact welding technique, known as the Vaporizing Foil Actuator Welding (VFAW), which uses the pressure created from the electrically driven rapid vaporization of a thin metallic conductor to drive the flyer. A typical characteristic of an impact-welded joint is the wavy weld interface which is attributed to the complex interfacial kinematics resulting from the high-speed oblique impact of the metal plates. A numerical simulation framework based on the finite element method (FEM) with Eulerian formalism is developed to model the high-speed impact between the metal plates. The model accounts for the thermomechanical interactions in the process and captures the complex interfacial deformations. A thorough validation of the model is achieved by comparing the wave characteristics obtained from numerical simulations with the experimental results from vaporizing foil actuator welding. The amplitude and wavelength of the interfacial waves and the resulting joint microstructure is shown to be strongly dependent on the process conditions and the specific material system, which is well captured by the model.

42 ENGINEERING↗
DOE OSTI2018

Direct observation of dislocation structure evolution in SRF cavity niobium using electron channeling contrast imaging

Two major goals of superconducting radio frequency cavity research and development are achieving higher accelerating gradient and gaining the highest quality factor, a measure of efficiency. However, the consistent improvement of these performance metrics is restricted by many factors, one of which is microstructural defects, such as dislocation substructures within the material. In this work, dislocation evolution is compared in four samples extracted from a 2.8 mm thick large-grain niobium slice, with tensile axes chosen to generate specific dislocation structures during subsequent deformation. The four samples are (1) as-extracted, (2) extracted and annealed, (3) extracted and then deformed to 40% tensile strain, and (4) extracted, annealed at 800 °C 2 h, and then deformed to 40% strain. Electron channeling contrast imaging was performed on all samples to characterize initial dislocation density and dislocation structure evolution due to annealing and deformation, and related to the mechanical behavior observed in stress-strain curves. A fundamental understanding of dislocation evolution in niobium is necessary to develop computational models to simulate cavity forming, which could enable new processing methods for cavity fabrication to be identified that could lead to improved performance.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2019

Surface‐Atmosphere Coupling Scale, the Fate of Water, and Ecophysiological Function in a Brazilian Forest

Abstract Tropical South America plays a central role in global climate. Bowen ratio teleconnects to circulation and precipitation processes far afield, and the global CO 2 growth rate is strongly influenced by carbon cycle processes in South America. However, quantification of basin‐wide seasonality of flux partitioning between latent and sensible heat, the response to anomalies around climatic norms, and understanding of the processes and mechanisms that control the carbon cycle remains elusive. Here, we investigate simulated surface‐atmosphere interaction at a single site in Brazil, using models with different representations of precipitation and cloud processes, as well as differences in scale of coupling between the surface and atmosphere. We find that the model with parameterized clouds/precipitation has a tendency toward unrealistic perpetual light precipitation, while models with explicit treatment of clouds produce more intense and less frequent rain. Models that couple the surface to the atmosphere on the scale of kilometers, as opposed to tens or hundreds of kilometers, produce even more realistic distributions of rainfall. Rainfall intensity has direct consequences for the “fate of water,” or the pathway that a hydrometeor follows once it interacts with the surface. We find that the model with explicit treatment of cloud processes, coupled to the surface at small scales, is the most realistic when compared to observations. These results have implications for simulations of global climate, as the use of models with explicit (as opposed to parameterized) cloud representations becomes more widespread.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Adaptively remeshed multiphysical modeling of resistance forge welding with experimental validation of residual stress fields and measurement processes

Welding processes used in the production of pressure vessels impart residual stresses in the manufactured component. Computational modeling is critical to predicting these residual stress fields and understanding how they interact with notches and flaws to impact pressure vessel durability. Here, in this work, we present a finite element model for a resistance forge weld and validate it using laboratory measurements. Extensive microstructural changes, near-melt temperatures, and large localized deformations along the weld interface pose significant challenges to Lagrangian finite element modeling. The proposed modeling approach overcomes these roadblocks in order to provide a high-fidelity simulation that can predict the residual stress state in the manufactured pressure vessel; a rich microstructural constitutive model accounts for material recrystallization dynamics, a frictional-to-tied contact model is coordinated with the constitutive model to represent interfacial bonding, and adaptive remeshing is employed to alleviate severe mesh distortion. An interrupted-weld approach is applied to the simulation to facilitate comparison to displacement measures. Several techniques are employed for residual stress measurement in order to validate the finite element model: neutron diffraction, the contour method, and the slitting method. Model-measurement comparisons are supplemented with detailed simulations that reflect the configurations of the residual-stress measurement processes themselves. The model results show general agreement with experimental measurements, and we observe some similarities in the features around the weld region. Factors that contribute to model-measurement differences are identified. Finally, we conclude with some discussion of the model development and residual stress measurement strategies, including how to best leverage the efforts put forth here for other weld problems.

36 MATERIALS SCIENCE↗
DOE OSTI2014

Sensitivity of Global Terrestrial Gross Primary Production to Hydrologic States Simulated by the Community Land Model Using Two Runoff Parameterizations

The terrestrial water and carbon cycles interact strongly at various spatio-temporal scales. To elucidate how hydrologic processes may influence carbon cycle processes, differences in terrestrial carbon cycle simulations induced by structural differences in two runoff generation schemes were investigated using the Community Land Model 4 (CLM4). Simulations were performed with runoff generation using the default TOPMODEL-based and the Variable Infiltration Capacity (VIC) model approaches under the same experimental protocol. The comparisons showed that differences in the simulated gross primary production (GPP) are mainly attributed to differences in the simulated leaf area index (LAI) rather than soil moisture availability. More specifically, differences in runoff simulations can influence LAI through changes in soil moisture, soil temperature, and their seasonality that affect the onset of the growing season and the subsequent dynamic feedbacks between terrestrial water, energy, and carbon cycles. As a result of a relative difference of 36% in global mean total runoff between the two models and subsequent changes in soil moisture, soil temperature, and LAI, the simulated global mean GPP differs by 20.4%. However, the relative difference in the global mean net ecosystem exchange between the two models is small (2.1%) due to competing effects on total mean ecosystem respiration and other fluxes, although large regional differences can still be found. Our study highlights the significant interactions among the water, energy, and carbon cycles and the need for reducing uncertainty in the hydrologic parameterization of land surface models to better constrain carbon cycle modeling.

carbon↗
DOE OSTI2025

Semi-supervised permutation invariant particle-level anomaly detection

The development of analysis methods to distinguish potential beyond the Standard Model phenomena in a model-agnostic way can significantly enhance the discovery reach in collider experiments. However, the typical machine learning (ML) algorithms employed for this task require fixed length and ordered inputs that break the natural permutation invariance in collision events. To address this, a semi-supervised anomaly detection tool is presented that takes a variable number of particle-level inputs and leverages a signal model to encode this information into a permutation invariant, event-level representation via supervised training with a Particle Flow Network (PFN). Data events are then encoded into this representation and given as input to an autoencoder for unsupervised ANomaly deTEction on particLe flOw latent sPacE (ANTELOPE), classifying anomalous events based on a low-level and permutation invariant input modeling. Performance of the ANTELOPE architecture is evaluated on simulated samples of hadronic processes in a high energy collider experiment, showing good capability to distinguish disparate models of new physics.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
NASA NTRS2006

A Mission to Observe Ice in Clouds from Space

To date there have been multiple satellite missions to observe and retrieve cloud top properties and the liquid in, and precipitation from, clouds. There are currently a few missions that attempt to measure cloud ice properties as a byproduct of other observations. However, we do not yet quantitatively understand the processes that control the water budget of the upper troposphere where ice is the predominant phase, and how these processes are linked to precipitation processes and the radiative energy budget. The ice in clouds either melts into rain or is detrained, and persists, as cirrus clouds affecting the hydrological and energy cycle, respectively. Fully modeling the Earth's climate and improving weather and climate forecasts requires accurate satellite measurements of various cloud properties at the temporal and spatial scales of cloud processes. The uncertainty in knowledge of these ice characteristics is reflected in the large discrepancies in model simulations of the upper tropospheric water budget. Model simulations are sensitive to the partition of ice between precipitation and outflow processes, i.e., to the parameterization of ice clouds and ice processes. This presentation will describe the Submillimeter-wave InfraRed Ice Cloud Experiment (SIRICE) concept, a satellite mission designed to acquire global Earth radiance measurements in the infrared and submillimeter-wave region (183-874 GHz). If successful, this mission will bridge the measurement gap between microwave sounders and shorter-wavelength infrared and visible sensors. The brightness temperatures at submillimeter-wave frequencies are especially sensitive to cirrus ice particle sizes (because they are comparable to the wavelength). This allows for more accurate ice water path estimates when multiple channels are used to probe into the cloud layers. Further, submillimeter wavelengths offer simplicity in the retrieval algorithms because they do not probe into the liquid and near surface portions of clouds, thus requiring only one term of the radiative transfer equation (ice scattering) to relate brightness temperatures to ice. Scientific justification and the SIRICE approach to measuring ice water path and particle size that span a range encompassing both the hydrologically active and radiatively active components of cloud systems will be presented.

Ackerman, S.↗
NASA NTRS2018

Radiative Transfer Modeling of Phytoplankton Fluorescence Quenching Processes

We report the first radiative transfer model that is able to simulate phytoplankton fluorescence with both photo chemical and non-photo chemical quenching included. The fluorescence source term in the inelastic radiative transfer equation is proportional to both the quantum yield and scalar irradiance at excitation wavelengths. The photo chemical and non photo chemical quenching processes change the quantum yield based on the photosynthetic active radiation. A sensitivity study was performed to demonstrate the dependence of the fluorescence signal on chlorophyll a concentration, aerosol optical depths and solar zenith angles. This work enables us to better model the phytoplankton fluorescence, which can be used in the design of new space-based sensors that can provide sufficient sensitivity to detect the phytoplankton fluorescence signal. It could also lead to more accurate remote sensing algorithms for the study of phytoplankton physiology.

Inelastic scattering↗
NASA NTRS1998

Numerical Simulations of a Lobed Fuel Injector

Numerical modeling of the nonreactive mixing processes associated with a lobed fuel injector in a coflowing air stream is presented. The lobed fuel injector is a device which generates strong streamwise vorticity, producing locally high strain rates which can enhance the molecular mixing of reactants while delaying ignition in a controlled manner. Vortex element modeling is used to simulate flow field evolution and fuel element mixing characteristics for this lobed fuel injector. Quantitative predictions for vorticity generation and qualitative results for streamwise rollup compare well qualitatively with recent experimental investigations of this flow field [Smith et al, Phys. Fluids 9, 667 (1997)]. Parametric studies of the effects of lobe amplitude-to-wavelength ratio, lobe angle, and lobe shape for given flow conditions suggest that geometrical features may be optimized to enhance mixing and control reaction processes.

Strickland, J. H.↗
DOE OSTI2018

On the Coupling Time of the Heat-Bath Process for the Fortuin–Kasteleyn Random–Cluster Model

We consider the coupling from the past implementation of the random–cluster heat-bath process, and study its random running time, or coupling time. We focus on hypercubic lattices embedded on tori, in dimensions one to three, with cluster fugacity at least one. We make a number of conjectures regarding the asymptotic behaviour of the coupling time, motivated by rigorous results in one dimension and Monte Carlo simulations in dimensions two and three. Amongst our findings, we observe that, for generic parameter values, the distribution of the appropriately standardized coupling time converges to a Gumbel distribution, and that the standard deviation of the coupling time is asymptotic to an explicit universal constant multiple of the relaxation time. Perhaps surprisingly, we observe these results to hold both off criticality, where the coupling time closely mimics the coupon collector’s problem, and also at the critical point, provided the cluster fugacity is below the value at which the transition becomes discontinuous. Finally, we consider analogous questions for the single-spin Ising heat-bath process.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2019

Subgrid variations of the cloud water and droplet number concentration over the tropical ocean: satellite observations and implications for warm rain simulations in climate models

Abstract. One of the challenges in representing warm rain processes in global climate models (GCMs) is related to the representation of the subgrid variability of cloud properties, such as cloud water and cloud droplet number concentration (CDNC), and the effect thereof on individual precipitation processes such as autoconversion. This effect is conventionally treated by multiplying the resolved-scale warm rain process rates by an enhancement factor ( E q ) which is derived from integrating over an assumed subgrid cloud water distribution. The assumed subgrid cloud distribution remains highly uncertain. In this study, we derive the subgrid variations of liquid-phase cloud properties over the tropical ocean using the satellite remote sensing products from Moderate Resolution Imaging Spectroradiometer (MODIS) and investigate the corresponding enhancement factors for the GCM parameterization of autoconversion rate. We find that the conventional approach of using only subgrid variability of cloud water is insufficient and that the subgrid variability of CDNC, as well as the correlation between the two, is also important for correctly simulating the autoconversion process in GCMs. Using the MODIS data which have near-global data coverage, we find that E q shows a strong dependence on cloud regimes due to the fact that the subgrid variability of cloud water and CDNC is regime dependent. Our analysis shows a significant increase of E q from the stratocumulus (Sc) to cumulus (Cu) regions. Furthermore, the enhancement factor E N due to the subgrid variation of CDNC is derived from satellite observation for the first time, and results reveal several regions downwind of biomass burning aerosols (e.g., Gulf of Guinea, east coast of South Africa), air pollution (i.e., East China Sea), and active volcanos (e.g., Kilauea, Hawaii, and Ambae, Vanuatu), where the E N is comparable to or even larger than E q , suggesting an important role of aerosol in influencing the E N . MODIS observations suggest that the subgrid variations of cloud liquid water path (LWP) and CDNC are generally positively correlated. As a result, the combined enhancement factor, including the effect of LWP and CDNC correlation, is significantly smaller than the simple product of E q ∙ E N . Given the importance of warm rain processes in understanding the Earth's system dynamics and water cycle, we conclude that more observational studies are needed to provide a better constraint on the warm rain processes in GCMs.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2015

GPU-based simulation of the two-dimensional unstable structure of gaseous oblique detonations

In this paper, the two-dimensional structure of unstable oblique detonations induced by the wedge from a supersonic combustible gas flow is simulated using the reactive Euler equations with a one-step Arrhenius chemistry model. A wide range of activation energy of the combustible mixture is considered. Computations are performed on the Graphical Processing Unit (GPU) to reduce the simulation runtimes. A large computational domain covered by a uniform mesh with high grid resolution is used to properly capture the development of instabilities and the formation of different transverse wave structures. After the initiation point, where the oblique shock transits into a detonation, an instability begins to manifest and in all cases, the left-running transverse waves first appear, followed by the subsequent emergence of right-running transverse waves forming the dual-head triple point structure. This study shows that for low activation energies, a long computational length must be carefully considered to reveal the unstable surface due to the slow growth rate of the instability. For high activation energies, the flow behind the unstable oblique detonation features the formation of unburnt gas pockets and strong vortex-pressure wave interaction resulting in a chaotic-like vortical structure.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS2018

Improving Spacecraft Design and Operability for Europa Clipper through High-Fidelity, Mission-Level Modeling and Simulation

NASA’s planned Europa Clipper mission seeks to assess the habitability of Jupiter’s moon Europa, which exhibits strong evidence of an ocean of liquid water underneath its icy crust. The sheer number of unique instruments, all of which require quiescent environments in order to operate, compounded with Jupiter’s distance from the Sun and Earth, make this planned mission challenging and resource-constrained. High fidelity, mission-level simulations that model the spacecraft, ground, and environment from launch to end-of-mission with a given trajectory and mission plan have been employed early in the project lifecycle to better understand the interactions between various components of the mission and how design changes impact the entire system. These simulations have already resulted in tangible benefits to the project by providing vital input to key spacecraft trades, assessing impacts to operability, and quantifying how well the scientific objectives of the mission can be achieved. Improvements to simulation performance and to the process by which information defining the system is gathered and built into models used by simulations have the potential to further expand the scope of their use on Europa Clipper and future missions.

Lawler, C.R.↗
DOE OSTI2021

Developing and testing capabilities for simulating cases with heterogeneous land/water surfaces in a novel atmospheric large eddy simulation code

Large eddy simulations (LES) are the primary computational tool used to simulate high Reynolds number three-dimensional turbulent flows. In the context of earth system sciences, particularly atmospheric science, LES are uniquely able to resolve the scales of atmospheric motion that are key for building process-level understanding of boundary layer turbulence, atmosphere-surface interaction, clouds, and cloud-aerosol-chemistry interaction, and are a core limited-area modeling capability. Increasing demands are being placed on LES code bases as growing high performance computing resources allow LES to address a wider range of scientific problems. In addition, LES are emerging as a source of high-quality machine learning training data. These demands necessitate an agile and extensible code base that allows the model to quickly adapt to emergent needs. However, LES have largely relied on legacy Fortran code bases that lack flexibility. A new, Python-based LES capability called Predicting INteractions of Aerosol and Clouds in Large Eddy Simulation (PINACLES) has been developed as part of the Department of Energy’s Earth System Model Development (ESMD) program area’s Enabling Aerosol-cloud interactions at Global convection-permitting scalES (EAGLES) project. PINACLES was developed from the ground up with a philosophy of maximizing scientific throughput, by attempting to optimize for both model throughput and software extensibility. The initial development of PINACLES delivered a state-of-the-art idealized LES capability solving the non-hydrostatic anelastic equations of motion with doubly periodic boundary conditions and idealized homogenous surface boundary conditions. Here we provide a final report on the outcomes of a fiscal year 2021 Seed Laboratory Directed Research Project that extended PINACLES in two key ways. First, PINACLES was coupled to a state-of-the-art land surface model enabling it to simulate spatially inhomogeneous land-atmosphere interactions that are known to control key atmospheric processes. Second, the dynamical core of PINACLES was modified to permit non-periodic boundary conditions. This model enhancement enables simulation of realistic cases with boundary conditions prescribed from atmospheric reanalysis and enables nested simulations conducted on a hierarchy of computational domains with increasing resolution. Together, these extensions to PINACLES make it a formidable modeling capability and expand its potential application to diverse components of DOE’s atmospheric science portfolio.

42 ENGINEERING↗
DOE OSTI2021

Marine ice sheet experiments with the Community Ice Sheet Model

Ice sheet models differ in their numerical treatment of dynamical processes. Simulations of marine-based ice are sensitive to the choice of Stokes flow approximation and basal friction law and to the treatment of stresses and melt rates near the grounding line. We study the effects of these numerical choices on marine ice sheet dynamics in the Community Ice Sheet Model (CISM). In the framework of the Marine Ice Sheet Model Intercomparison Project 3d (MISMIP3d), we show that a depth-integrated, higher-order solver gives results similar to a 3D (Blatter–Pattyn) solver. We confirm that using a grounding line parameterization to approximate stresses in the grounding zone leads to accurate representation of ice sheet flow with a resolution of ~2 km, as opposed to ~0.5 km without the parameterization. In the MISMIP+ experimental framework, we compare different treatments of sub-shelf melting near the grounding line. In contrast to recent studies arguing that melting should not be applied in partly grounded cells, it is usually beneficial in CISM simulations to apply some melting in these cells. This suggests that the optimal treatment of melting near the grounding line can depend on ice sheet geometry, forcing, or model numerics. In both experimental frameworks, ice flow is sensitive to the choice of basal friction law. To study this sensitivity, we evaluate friction laws that vary the connectivity between the basal hydrological system and the ocean near the grounding line. CISM yields accurate results in steady-state and perturbation experiments at a resolution of ~2 km (arguably 4 km) when the connectivity is low or moderate and ~1 km (arguably 2 km) when the connectivity is strong.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2017

Coupling a three-dimensional subsurface flow and transport model with a land surface model to simulate stream–aquifer–land interactions (CP v1.0)

Abstract. A fully coupled three-dimensional surface and subsurface land model is developed and applied to a site along the Columbia River to simulate three-way interactions among river water, groundwater, and land surface processes. The model features the coupling of the Community Land Model version 4.5 (CLM4.5) and a massively parallel multiphysics reactive transport model (PFLOTRAN). The coupled model, named CP v1.0, is applied to a 400 m × 400 m study domain instrumented with groundwater monitoring wells along the Columbia River shoreline. CP v1.0 simulations are performed at three spatial resolutions (i.e., 2, 10, and 20 m) over a 5-year period to evaluate the impact of hydroclimatic conditions and spatial resolution on simulated variables. Results show that the coupled model is capable of simulating groundwater–river-water interactions driven by river stage variability along managed river reaches, which are of global significance as a result of over 30 000 dams constructed worldwide during the past half-century. Our numerical experiments suggest that the land-surface energy partitioning is strongly modulated by groundwater–river-water interactions through expanding the periodically inundated fraction of the riparian zone, and enhancing moisture availability in the vadose zone via capillary rise in response to the river stage change. Meanwhile, CLM4.5 fails to capture the key hydrologic process (i.e., groundwater–river-water exchange) at the site, and consequently simulates drastically different water and energy budgets. Furthermore, spatial resolution is found to significantly impact the accuracy of estimated the mass exchange rates at the boundaries of the aquifer, and it becomes critical when surface and subsurface become more tightly coupled with groundwater table within 6 to 7 meters below the surface. Inclusion of lateral subsurface flow influenced both the surface energy budget and subsurface transport processes as a result of river-water intrusion into the subsurface in response to an elevated river stage that increased soil moisture for evapotranspiration and suppressed available energy for sensible heat in the warm season. The coupled model developed in this study can be used for improving mechanistic understanding of ecosystem functioning and biogeochemical cycling along river corridors under historical and future hydroclimatic changes. The dataset presented in this study can also serve as a good benchmarking case for testing other integrated models.

58 GEOSCIENCES↗
DOE OSTI2017

Coupling a three-dimensional subsurface flow and transport model with a land surface model to simulate stream–aquifer–land interactions (CP v1.0)

A fully coupled three-dimensional surface and subsurface land model is developed and applied to a site along the Columbia River to simulate three-way interactions among river water, groundwater, and land surface processes. The model features the coupling of the Community Land Model version 4.5 (CLM4.5) and a massively parallel multiphysics reactive transport model (PFLOTRAN). The coupled model, named CP v1.0, is applied to a 400 m × 400 m study domain instrumented with groundwater monitoring wells along the Columbia River shoreline. CP v1.0 simulations are performed at three spatial resolutions (i.e., 2, 10, and 20 m) over a 5-year period to evaluate the impact of hydroclimatic conditions and spatial resolution on simulated variables. Results show that the coupled model is capable of simulating groundwater–river-water interactions driven by river stage variability along managed river reaches, which are of global significance as a result of over 30 000 dams constructed worldwide during the past half-century. Our numerical experiments suggest that the land-surface energy partitioning is strongly modulated by groundwater–river-water interactions through expanding the periodically inundated fraction of the riparian zone, and enhancing moisture availability in the vadose zone via capillary rise in response to the river stage change. Meanwhile, CLM4.5 fails to capture the key hydrologic process (i.e., groundwater–river-water exchange) at the site, and consequently simulates drastically different water and energy budgets. Furthermore, spatial resolution is found to significantly impact the accuracy of estimated the mass exchange rates at the boundaries of the aquifer, and it becomes critical when surface and subsurface become more tightly coupled with groundwater table within 6 to 7 meters below the surface. Inclusion of lateral subsurface flow influenced both the surface energy budget and subsurface transport processes as a result of river-water intrusion into the subsurface in response to an elevated river stage that increased soil moisture for evapotranspiration and suppressed available energy for sensible heat in the warm season. The coupled model developed in this study can be used for improving mechanistic understanding of ecosystem functioning and biogeochemical cycling along river corridors under historical and future hydroclimatic changes. The dataset presented in this study can also serve as a good benchmarking case for testing other integrated models.

58 GEOSCIENCES↗
DOE OSTI2014

Functional approach to exploring climatic and landscape controls of runoff generation: 1. Behavioral constraints on runoff volume

Abstract Inspired by the Dunne diagram, the climatic and landscape controls on the partitioning of annual runoff into its various components (Hortonian and Dunne overland flow and subsurface stormflow) are assessed quantitatively, from a purely theoretical perspective. A simple distributed hydrologic model has been built sufficient to simulate the effects of different combinations of climate, soil, and topography on the runoff generation processes. The model is driven by a sequence of simple hypothetical precipitation events, for a large combination of climate and landscape properties, and hydrologic responses at the catchment scale are obtained through aggregation of grid‐scale responses. It is found, first, that the water balance responses, including relative contributions of different runoff generation mechanisms, could be related to a small set of dimensionless similarity parameters. These capture the competition between the wetting, drying, storage, and drainage functions underlying the catchment responses, and in this way, provide a quantitative approximation of the conceptual Dunne diagram. Second, only a subset of all hypothetical catchment/climate combinations is found to be “behavioral,” in terms of falling sufficiently close to the Budyko curve, describing mean annual runoff as a function of climate aridity. Furthermore, these behavioral combinations are mostly consistent with the qualitative picture presented in the Dunne diagram, indicating clearly the commonality between the Budyko curve and the Dunne diagram. These analyses also suggest clear interrelationships amongst the “behavioral” climate, soil, and topography parameter combinations, implying these catchment properties may be constrained to be codependent in order to satisfy the Budyko curve.

Li, Hong‐Yi↗