Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “functional 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 217 records · Page 12
All sourcesNASA NTRSDOE OSTI
NASA NTRS1984

Statistical Time Series Models of Pilot Control with Applications to Instrument Discrimination

A general description of the methodology used in obtaining the transfer function models and verification of model fidelity, frequency domain plots of the modeled transfer functions, numerical results obtained from an analysis of poles and zeroes obtained from z plane to s-plane conversions of the transfer functions, and the results of a study on the sequential introduction of other variables, both exogenous and endogenous into the loop are contained.

Altschul, R. E.↗
DOE OSTI2022

Comparing Model Representations of Physiological Limits on Transpiration at a Semi-arid Ponderosa Pine Site

Mechanistic representations of biogeochemical processes in ecosystem models are rapidly advancing, requiring advancements in model evaluation approaches. Here we quantify multiple aspects of model functional performance to evaluate improved process representations in ecosystem models. We compare semi-empirical stomatal models with hydraulic constraints against more mechanistic representations of stomatal and hydraulic functioning at a semi-arid pine site using a suite of metrics and analytical tools. We find that models generally perform similarly under unstressed conditions, but performance diverges under atmospheric and soil drought. The more empirical models better capture synergistic information flows between soil water potential and vapor pressure deficit to transpiration, while the more mechanistic models are overly deterministic. Although models can be parameterized to yield similar functional performance, alternate parameterizations could not overcome structural model constraints that underestimate the unique information contained in soil water potential about transpiration. Additionally, both multilayer canopy and big-leaf models were unable to capture the magnitude of canopy temperature divergence from air temperature, and we demonstrate that errors in leaf temperature can propagate to considerable error in simulated transpiration. This study demonstrates the value of merging underutilized observational data streams with emerging analytical tools to characterize ecosystem function and discriminate among model process representations.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRSDate not supplied

Integrated Bosch Process System Models for In-Situ Oxygen and Carbon Production

In-Situ Resource Utilization (ISRU) technology is a vital component to NASA’s mission of a sustainable presence on the Moon and Mars. Local resources can be leveraged to reduce resupply frequency and mass. Elements of the Bosch process, combined with the carbothermal reduction process, can produce oxygen on the lunar surface with minimal consumables. The Bosch process can also produce oxygen on the Martian surface by using the CO 2 -rich environment. Between both systems, adsorption pump, solar thermal energy, carbon formation reactor, and water recovery subsystems are modeled and integrated to create a functional model in MATLAB software. The model is used to simulate performance of the system and reduce mass, power, and volume requirements. This integrated system model provides a tool to scale ISRU technologies for oxygen and carbon production. The MATLAB model is created by developing a system of independent subsystem models that are solved for their quasi-steady state values which can be integrated with respect to time to determine the change in current states. A flexible time stepping method is used to ensure a high level of accuracy during periods of rapid change while still making use of a simple explicit integration method. The flexible time step is calculated for each independent subsystem and the minimum value from those is used as the overall time step. A flexible time step is calculated by dividing a resolution value, or the maximum change per time step, by the variables current rate of change. The maximum value from all points in space is used for subsystem models that contain multiple values. The process is done for every variable that is being monitored in each subsystem and the global minimum is used as that iteration’s timestep. Several assumptions used in the MATLAB model for fluid flow dynamics, such as 1-D gas flow through the sorption pump, are supported by modeling in Ansys Fluent software. The Lunar oxygen production system is outlined in Fig. 1. The carbothermal reduction subsystem uses solar energy to heat a mixture of lunar regolith and carbon powder to produce carbon monoxide. To begin, the carbon monoxide feeds to the modified Bosch subsystem along with hydrogen gas. The reactants then enter the carbon formation reactor where water and carbon powder are produced. Solar thermal energy is used to add energy to the reactor, but waste heat from the carbothermal process is another potential heat source. The water is collected and electrolyzed to produce hydrogen which reenters the Bosch subsystem, and the oxygen is stored for downstream use. The carbon powder is collected and feeds back into the carbothermal subsystem. The Martian oxygen production system uses the full Bosch process and is outlined in Fig 2. A CO 2 adsorption pump thermally cycles to scrub and pressurize CO 2 from the environment. Along with an initial supply of hydrogen, the reactants enter the Reverse Water Gas Shift Reactor (RWGSR) which produces carbon monoxide and water. Carbon monoxide and unreacted hydrogen enter the carbon formation reactor to produce water and carbon powder. The water is collected from both reactors and electrolyzed to reintroduce hydrogen and store oxygen for propellant production or life support. Carbon is removed from the carbon formation reactor and stored. The adsorption pump utilizes rapid cycle temperature swings within a stack of zeolite coated surfaces. The subsystem model solves 1-D quasi-steady conservation laws of the quasi-steady form, shown in Eq. 1, for the gas stream and heat exchange liquid to predict performance parameters such as breakthrough capacity and optimum cycle time. The source term S is used to capture interactions between the fluid flows and the sorbent. A quasi-steady-state scheme is used where no time derivatives appear in the governing equations, except for those in the source terms. This results in an autonomous system, where ∂F/∂x = ƒ(F). The fluxes F are provided at the inlet, and an explicit method is used to solve for the spatial distribution of F. The heat and mass flows to the sorbent are then extracted from the source terms. These flows are numerically integrated to produce a 1-D solution for the system’s state as a function of both time and space. The body of the adsorption pump is separated into two semi-independent models: the heat exchanger fluid flow and gas flow through the zeolite coated surfaces. Both models are solved using the above-described method to find a 1-D solution as a function of space and interact only once a timestep is taken. The interaction point is the sorbent through which all heat transfer between the two models must occur. Sorbent mass adsorption is calculated using the Lagergren model, shown in Eq. 2, where the transfer coefficient, λ D , is found by solving a system of nondimensionalized equations derived by using the heat and mass transfer analogy for transport phenomena. Using Grade 544 Type 13X zeolite as the sorbent material, the equilibrium concentration, θ eq , is calculated using the k-site Langmuir isotherm and fit parameters. Additionally, the enthalpy of adsorption used in the model is computed by interpolation of available data [1]. The subsystem model was validated using the Rapid Cycle Temperature Swing Adsorption (RC-TSA) pump. The solar thermal energy subsystem focuses on a solar concentrator concept with a heat exchanger to heat the reactants before entering the carbon formation reactor. The subsystem model assumes a fixed solar flux and reflector efficiency to calculate the reactant temperature given the incoming temperature, pressure, and exchanger geometry. The receiver is a custom manufactured series of copper blocks with serpentine channels to increase its surface area and the residence time of the reactants to heat up to 550 °C. The subsystem model was validated using a heat exchanger developed at NASA Glenn Research Center (GRC). The solar thermal energy subsystem focuses on a solar concentrator concept with a heat exchanger to heat the reactants before entering the carbon formation reactor. The subsystem model assumes a fixed solar flux and reflector efficiency to calculate the reactant temperature given the incoming temperature, pressure, and exchanger geometry. The receiver is a custom manufactured series of copper blocks with serpentine channels to increase its surface area and the residence time of the reactants to heat up to 550 °C. The subsystem model was validated using a heat exchanger developed at NASA Glenn Research Center (GRC).

In situ Resource Utilization↗
NASA NTRS1993

Fusing modeling techniques to support domain analysis for reuse opportunities identification

Functional modeling techniques or object-oriented graphical representations, which are more useful to someone trying to understand the general design or high level requirements of a system? For a recent domain analysis effort, the answer was a fusion of popular modeling techniques of both types. By using both functional and object-oriented techniques, the analysts involved were able to lean on their experience in function oriented software development, while taking advantage of the descriptive power available in object oriented models. In addition, a base of familiar modeling methods permitted the group of mostly new domain analysts to learn the details of the domain analysis process while producing a quality product. This paper describes the background of this project and then provides a high level definition of domain analysis. The majority of this paper focuses on the modeling method developed and utilized during this analysis effort.

Hall, Susan Main↗
DOE OSTI2026

Fit-free optical determination of electronic thermalization time in nematic iron-based superconductors

We present a nematic response function model (NRFM) for fit-free direct extraction of the characteristic time of ultrafast electronic thermalization in iron-based superconductors, materials with electronic nematicity. By combining the NRFM for polarization-dependent pump–probe measurements of electronic nematic response with the two-temperature model (TTM) for sub-picosecond quasiparticle relaxation, we quantify the electronic thermalization timescales and their anisotropy. The nematic response function is modeled as the difference in normalized reflectivity signals, revealing a pronounced sub-picosecond extremum in signal evolution that directly yields the characteristic electronic thermalization time. This method demonstrates that the NRFM is consistent with TTM fits of transient optical response, yielding electronic thermalization time constants on the order of 110–230 fs for the FeSe 1−x Te x and Ba(Fe 0.92 Co 0.08 ) 2 As 2 thin films. Here, the proposed approach can be applied to any material that exhibits electronic nematicity, providing a powerful tool for direct mapping of the relaxation time in nematic materials, avoiding complex experimental data-fitting procedures.

Bartenev, Alexander [University of Puerto Rico, Ma↗
DOE OSTI2020

Non-Gaussianity in the weak lensing correlation function likelihood – implications for cosmological parameter biases

ABSTRACT We study the significance of non-Gaussianity in the likelihood of weak lensing shear two-point correlation functions, detecting significantly non-zero skewness and kurtosis in 1D marginal distributions of shear two-point correlation functions in simulated weak lensing data. We examine the implications in the context of future surveys, in particular LSST, with derivations of how the non-Gaussianity scales with survey area. We show that there is no significant bias in 1D posteriors of Ωm and σ8 due to the non-Gaussian likelihood distributions of shear correlations functions using the mock data (100 deg2). We also present a systematic approach to constructing approximate multivariate likelihoods with 1D parametric functions by assuming independence or more flexible non-parametric multivariate methods after decorrelating the data points using principal component analysis (PCA). While the use of PCA does not modify the non-Gaussianity of the multivariate likelihood, we find empirically that the 1D marginal sampling distributions of the PCA components exhibit less skewness and kurtosis than the original shear correlation functions. Modelling the likelihood with marginal parametric functions based on the assumption of independence between PCA components thus gives a lower limit for the biases. We further demonstrate that the difference in cosmological parameter constraints between the multivariate Gaussian likelihood model and more complex non-Gaussian likelihood models would be even smaller for an LSST-like survey. In addition, the PCA approach automatically serves as a data compression method, enabling the retention of the majority of the cosmological information while reducing the dimensionality of the data vector by a factor of ∼5.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1965

Human Tracking Performance in Uncoupled and Coupled Two-Axis Systems

This report presents tile results of an experimental and analytical study of human performance in uncoupled and coupled control systems. Human pilot performance in single and two-axis systems was mathematically modeled by linear second-order describing functions. Model parameters were determined using model matching techniques. Analysis of the models showed that the amplitude ratio and phase lead of the describing function increased with training indicating an increase in open loop bandwidth. The phase margin also decreased with training. Increasing the plant lag time constant resulted in an increase in the model lead time constant and a decrease in the zero frequency gain. No significant difference was found to exist in the normalized tracking error per axis between the two-axis tasks and the single-axis tasks. However tile model lead time constant was significantly greater in two-axis tracking. Manual tracking of two-axis systems with cross-coupling was studied experimentally and analytically. Approximate methods for modeling two-axis performance were developed and checked using a precise spectral analysis approach. Coupled and uncoupled, symmetrical and asymmetrical two-axis performance was compared. The results show that modeling of cross-coupled systems is feasible and that trained subjects are capable of decoupling the axes of some systems. A methodology study compared the identification performance of continuous, iterative, and extrapolation model matching techniques. An iterative technique employing sensitivity equations for the generation of influence coefficients was found to be the best technique due to its excellent identification accuracy and ease of implementation. Convergence in iterative techniques can be improved substantially by equalizing the parameter adjustment rates and limiting the maximum correction per iteration.

Todosiev, E. P.↗
NASA NTRS1991

Relation between the quasi-periodic oscillations and the low-frequency noise of GX 5-1 in the horizontal branch

Ginga observations of quasi-periodic oscillations (QPOs) and the low-frequency noise (LFN) from GX 5-1 in its horizontal-branch spectral state are presented. Power spectral fits were attempted using model functions based on simple oscillating shot models. A clear second-harmonic peak of QPO was detected. Variations in the powers of QPO and LFN on timescales of 8-256 s were also studied. These variations were significant for all of the timescales studied, and were uncorrelated with each other on timescales shorter than a few tens of seconds, and correlated on longer timescales. From simulations based on a simple shot model, it was found that the variation amplitude and the lack of correlation on short timescales are not inconsistent with the oscillating shot models. A more complex model is necessary to fully explain the observed properties.

Mitsuda, Kazuhisa↗
NASA NTRS2014

Aerosol Remote Sensing in Polar Regions

Multi-year sets of ground-based sun-photometer measurements conducted at 12 Arctic sites and 9 Antarctic sites were examined to determine daily mean values of aerosol optical thickness tau(lambda) at visible and near-infrared wavelengths, from which best-fit values of Ångström's exponent alpha were calculated. Analyzing these data, the monthly mean values of tau(0.50 micrometers) and alpha and the relative frequency histograms of the daily mean values of both parameters were determined for winter-spring and summer-autumn in the Arctic and for austral summer in Antarctica. The Arctic and Antarctic covariance plots of the seasonal median values of alpha versus tau(0.50 micrometers) showed: (i) a considerable increase in tau(0.50 micrometers) for the Arctic aerosol from summer to winter-spring, without marked changes in alpha; and (ii) a marked increase in tau(0.50 micrometer) passing from the Antarctic Plateau to coastal sites, whereas alpha decreased considerably due to the larger fraction of sea-salt aerosol. Good agreement was found when comparing ground-based sun-photometer measurements of tau(lambda) and alpha at Arctic and Antarctic coastal sites with Microtops measurements conducted during numerous AERONET/MAN cruises from 2006 to 2013 in three Arctic Ocean sectors and in coastal and off-shore regions of the Southern Atlantic, Pacific, and Indian Oceans, and the Antarctic Peninsula. Lidar measurements were also examined to characterize vertical profiles of the aerosol backscattering coefficient measured throughout the year at Ny-Ålesund. Satellite-based MODIS, MISR, and AATSR retrievals of tau(lambda) over large parts of the oceanic polar regions during spring and summer were in close agreement with ship-borne and coastal ground-based sun-photometer measurements. An overview of the chemical composition of mode particles is also presented, based on in-situ measurements at Arctic and Antarctic sites. Fourteen log-normal aerosol number size-distributions were defined to represent the average features of nuclei, accumulation and coarse mode particles for Arctic haze, summer background aerosol, Asian dust and boreal forest fire smoke, and for various background austral summer aerosol types at coastal and high-altitude Antarctic sites. The main columnar aerosol optical characteristics were determined for all 14 particle modes, based on in-situ measurements of the scattering and absorption coefficients. Diurnally averaged direct aerosol-induced radiative forcing and efficiency were calculated for a set of multimodal aerosol extinction models, using various Bidirectional Reflectance Distribution Function models over vegetation-covered, oceanic and snow-covered surfaces. These gave a reliable measure of the pronounced effects of aerosols on the radiation balance of the surface-atmosphere system over polar regions.

Sun-photometer measurements↗
NASA NTRS2003

Oxidation of Carbon Fibers in a Cracked Ceramic Matrix Composite Modeled as a Function of Temperature

The oxidation model simulates the oxidation of the reinforcing carbon fibers within a ceramic matrix composite material containing as-fabricated microcracks. The physics-based oxidation model uses theoretically and experimentally determined variables as input for the model. The model simulates the ingress of oxygen through microcracks into a two-dimensional plane within the composite material. Model input includes temperature, oxygen concentration, the reaction rate constant, the diffusion coefficient, and the crack opening width as a function of the mechanical and thermal loads. The model is run in an iterative process for a two-dimensional grid system in which oxygen diffuses through the porous and cracked regions of the material and reacts with carbon in short time steps. The model allows the local oxygen concentrations and carbon volumes from the edge to the interior of the composite to be determined over time. Oxidation damage predicted by the model was compared with that observed from microstructural analysis of experimentally tested composite material to validate the model for two temperatures of interest. When the model is run for low-temperature conditions, the kinetics are reaction controlled. Carbon and oxygen reactions occur relatively slowly. Therefore, oxygen can bypass the carbon near the outer edge and diffuse into the interior so that it saturates the entire composite at relatively high concentrations. The kinetics are limited by the reaction rate between carbon and oxygen. This results in an interior that has high local concentrations of oxygen and a similar amount of consumed carbon throughout the cross section. When the model is run for high-temperature conditions, the kinetics are diffusion controlled. Carbon and oxygen reactions occur very quickly. The carbon consumes oxygen as soon as it is supplied. The kinetics are limited by the relatively slow rate at which oxygen is supplied in comparison to the relatively fast rate at which carbon and oxygen reactions occur. This results in a sharp gradient in oxygen concentration from the edge where it is supplied to the nearest source of carbon, which is where the oxygen is quickly consumed. A moving reaction front is seen in which the outlaying carbon is consumed before the next inner layer of carbon begins to react.

Halbig, Michael C.↗
NASA NTRS1984

Thematic Mapper Protoflight Model Line Spread Function

The Thematic Mapper (TM) Protoflight Model Spatial Line Spread Function (LSF) was not measured before launch. Therefore, methodology are developed to characterize LSF with protoflight model optics and electronics measurements that were made before launch. Direct prelaunch LSF measurements that were made from the flight model TM verified the protoflight TM LSF simulation. Results for two selected protoflight TM channels are presented here. It is shown that LSF data for the other ninety-four channels could be generated in the same fashion.

Schueler, C.↗
NASA NTRS1974

A representation of Jacchia's thermospheric models in spherical harmonic functions

The Jacchia models are represented in terms of spherical harmonic functions. This representation has the advantage of ease of comparison with other global theoretical and empirical models that use this mathematical form. Furthermore, it is analytic, continuous, and has continuous derivatives all over the globe. The representation of the exospheric temperatures shows clearly the amplitudes of the various periodic terms and uses relatively few constants. An example of a similar representation for the total mass density at a particular height and level of solar activity is given as well.

Blum, P.↗
DOE OSTI2025

Impacts of benchmarking choices on inferred model skill of the Arctic–Boreal terrestrial carbon cycle

Abstract Land surface models require continuous validation against observations to improve and reduce simulation uncertainty. However, inferred model performance can be heavily influenced by subjective choices made in the selection and application of observational data products. A key area often misrepresented by models is the Arctic–Boreal region, which is a potential tipping point region in Earth’s climate system due to large permafrost carbon stocks that are vulnerable to release with climate warming. We use the International Land Model Benchmarking (ILAMB) framework to evaluate how the model skill of TRENDY-v9 models varies based on the choice of observational-based benchmark and how benchmarks are applied in model evaluation. This analysis uses global datasets integrated into ILAMB and new, regionally-specific observational products from the Arctic–Boreal Vulnerability Experiment. Our results cover the overall time period of 1979–2019 and show that model scores can vary substantially depending on the data product applied, with higher model scores indicating better model performance against observations. The lowest model scores occur when benchmarked against regional, compared to global, datasets. We also evaluate observed and modeled functional relationships between ecosystem respiration and air temperature and between gross primary production and precipitation. Here, we find that the magnitude and shape of the responses are strongly impacted by the choice of observational dataset and the approach used to construct the functional relationship benchmark. These results suggest that model evaluation studies could conclude a false sense of model skill if only using a single benchmark data product or if not applying regional data products when performing a regional model analysis. Collectively, our findings highlight the influence of benchmarking choices on model evaluation and point to the need for benchmarking guidelines when assessing model skill.

Poe, Jeralyn (ORCID:0000000318495278)↗
DOE OSTI2023

Practical considerations in using a binary pseudorandom array for instrument transfer function calibration

Binary pseudo-random array (BPRA) artifacts are useful devices for calibrating the instrument transfer function (ITF) of interferometric microscopes and other optical and non-optical surface and wavefront measurement instruments. The intrinsic white noise character of the power spectral density function of the artifact simplifies the deconvolution of the ITF from the measured power spectral density (PSD). However, resampling of the BPRA intrinsic artifact features with the measurement tool's specific sampling pattern modifies the white noise character of the intrinsic spectrum and needs to be accounted for in the ITF-based data deconvolution process. We have developed an analytic solution to the spectrum of a resampled one- and two- dimensional BPRA. The resultant nominal PSD function is a simple twoparameter cosine function with a period equal to the resampled pixel width. A transfer function model for interferometric microscopes that incorporates this function, along with an ITF that includes aliasing effects and variable numerical aperture (NA), wavelength, and obscuration factor, is used to fit to the BPRA PSDs measured by an interference microscope for a range of objective and zoom lens magnification combinations.

Takacs, PZ↗
NASA NTRS1989

Use of artificial intelligence in supervisory control

Viewgraphs describing the design and testing of an intelligent decision support system called OFMspert are presented. In this expert system, knowledge about the human operator is represented through an operator/system model referred to as the OFM (Operator Function Model). OFMspert uses the blackboard model of problem solving to maintain a dynamic representation of operator goals, plans, tasks, and actions given previous operator actions and current system state. Results of an experiment to assess OFMspert's intent inferencing capability are outlined. Finally, the overall design philosophy for an intelligent tutoring system (OFMTutor) for operators of complex dynamic systems is summarized.

Cohen, Aaron↗
NASA NTRS1987

Investigation of turbulent flow in highly curved ducts with application to turbomachinery components

Numerical predictions were performed using a semi-elliptic calculation procedure for the case of turbulent flow in passage through a 90 deg bend of square cross section. Two versions of the isotropic turbulent viscosity two equation k-epsilon model were used. The first, the wall function model (WFM), employs the logarithmic law-of-the-wall relation and the notion of equilibrium flow to set all the necessary boundary conditions at the first grid point adjacent to a solid wall. The second model, the Van Driest model (VDM), employs Prandtl's original mixing length formulation, in conjunction with Van Driest's semi-empirical relation for the mixing length, to calculate the turbulent viscosity in the near wall regions of the flow. In this case, boundary conditions for k and epsilon, required to calculate these quantities in the core of the flow, are obtained by matching the mixing length and Reynolds number model formulations in an overlapping region of the flow near the walls. In both cases the results obtained show an improvement over earlier calculations using an elliptic numerical procedure. This is attributed to the finer grids possible in the present work. Of the two models, the VDM formulation shows better overall conformity with the mean flow measurements. Neither model reproduces well the details of the stress distribution as a result of the implied isotropic turbulent viscosity.

Humphrey, J. A. C.↗
NASA NTRS1992

Adaptation of multidimensional group particle tracking and particle wall-boundary condition model to the FDNS code

A particulate two-phase flow CFD model was developed based on the FDNS code which is a pressure based predictor plus multi-corrector Navier-Stokes flow solver. Turbulence models with compressibility correction and the wall function models were employed as submodels. A finite-rate chemistry model was used for reacting flow simulation. For particulate two-phase flow simulations, a Eulerian-Lagrangian solution method using an efficient implicit particle trajectory integration scheme was developed in this study. Effects of particle-gas reaction and particle size change to agglomeration or fragmentation were not considered in this investigation. At the onset of the present study, a two-dimensional version of FDNS which had been modified to treat Lagrangian tracking of particles (FDNS-2DEL) had already been written and was operational. The FDNS-2DEL code was too slow for practical use, mainly because it had not been written in a form amenable to vectorization on the Cray, nor was the full three-dimensional form of FDNS utilized. The specific objective of this study was to reorder to calculations into long single arrays for automatic vectorization on the Cray and to implement the full three-dimensional version of FDNS to produce the FDNS-3DEL code. Since the FDNS-2DEL code was slow, a very limited number of test cases had been run with it. This study was also intended to increase the number of cases simulated to verify and improve, as necessary, the particle tracking methodology coded in FDNS.

Chen, Y. S.↗
NASA NTRS1993

Evapotranspiration and runoff from large land areas: Land surface hydrology for atmospheric general circulation models

A land surface hydrology parameterization for use in atmospheric GCM's is presented. The parameterization incorporates subgrid scale variability in topography, soils, soil moisture and precipitation. The framework of the model is the statistical distribution of a topography-soils index, which controls the local water balance fluxes, and is therefore taken to represent the large land area. Spatially variable water balance fluxes are integrated with respect to the topography-soils index to yield our large topography-soils distribution, and interval responses are weighted by the probability of occurrence of the interval. Grid square averaged land surface fluxes result. The model functions independently as a macroscale water balance model. Runoff ratio and evapotranspiration efficiency parameterizations are derived and are shown to depend on the spatial variability of the above mentioned properties and processes, as well as the dynamics of land surface-atmosphere interactions.

Famiglietti, J. S.↗