Search NASA⌕ Search

SEARCH · Search NASA

Results for “root 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 775 records · Page 43

Investigation of Abnormal Level Control Oscillations in a BWR Feedwater System

In the long-term operation of nuclear power plants, the aging of systems, structures, and components can lead to maintenance issues that must be dealt with to maintain cost-effective plant operations. One common issue affecting the currently operated boiling water reactors is the onset of unexpected level oscillations in feedwater heaters. This phenomenon can cause excessive cycling of drain valves and lead to premature failures. In this work, we develop a dynamic model of a set of feedwater heaters to determine the root cause of oscillations observed in an operating plant. Simulation results of various transient scenarios were used to investigate the effects of the controller parameters, boundary conditions, and possible valve and instrument issues. The analysis led to the conclusion that the most likely causes of the observed self-sustained oscillations in the system are the nonlinear behaviors of the drain valve and the level transmitter induced by degraded equipment condition. In conclusion, a partial plug of the pressure line used for level sensing in the system can account for a significant deadtime in the level transmitter, a nonlinear effect shown to induce self-sustained oscillatory behaviors.

Boiling water reactors↗

Updated Economic Model for Estimation of GDP Losses in the MACCS Offsite Consequence Analysis Code RDEIM Model Report for MACCS v4.2

This report updates the Regional Disruption Economic Impact Model (RDEIM) GDP-based model described in Bixler et al. (2020) used in the MACCS accident consequence analysis code. MACCS is the U.S. Nuclear Regulatory Commission (NRC) used to perform probabilistic health and economic consequence assessments for atmospheric releases of radionuclides. It is also used by international organizations, both reactor owners and regulators. It is intended and most commonly used for hypothetical accidents that could potentially occur in the future rather than to evaluate past accidents or to provide emergency response during an ongoing accident. It is designed to support probabilistic risk and consequence analyses and is used by the NRC, U.S. nuclear licensees, the Department of Energy, and international vendors, licensees, and regulators. The update of the RDEIM model in version 4.2 expresses the national recovery calculation explicitly, rather than implicitly as in the previous version. The calculation of the total national GDP losses remains unchanged. However, anticipated gains from recovery are now allocated across all the GDP loss types – direct, indirect, and induced – whereas in version 4.1, all recovery gains were accounted for in the indirect loss type. To achieve this, we’ve introduced new methodology to streamline and simplify the calculation of all types of losses and recovery. In addition, RDEIM includes other kinds of losses, including tangible wealth. This includes loss of tangible assets (e.g., depreciation) and accident expenditures (e.g., decontamination). This document describes the updated RDEIM economic model and provides examples of loss and recovery calculation, results analysis, and presentation. Changes to the tangible cost calculation and accident expenditures are described in section 2.2. The updates to the RDEIM input-output (I-O) model are not expected to affect the final benchmark results Bixler et al. (2020), as the RDEIM calculation for the total national GDP losses remains unchanged. The reader is referred to the MACCS revision history for other cost modelling changes since version 4.0 that may affect the benchmark. RDEIM has its roots in a code developed by Sandia National Laboratories for the Department of Homeland Security to estimate short-term losses from natural and manmade accidents, called the Regional Economic Accounting analysis tool (REAcct). This model was adapted and modified for MACCS. It is based on I-O theory, which is widely used in economic modeling. It accounts for direct losses to a disrupted region affected by an accident, indirect losses to the national economy due to disruption of the supply chain, and induced losses from reduced spending by displaced workers. RDEIM differs from REAcct in in its treatment and estimation of indirect loss multipliers, elimination of double-counting associated with inter-industry trade in the affected area, and that it is intended to be used for extended periods that can occur from a major nuclear reactor accident, such as the one that occurred at the Fukushima Daiichi site in Japan. Most input-output models do not account for economic adaptation and recovery, and in this regard RDEIM differs from its parent, REAcct, because it allows for a user-definable national recovery period. Implementation of a recovery period was one of several recommendations made by an independent peer review panel to ensure that RDEIM is state-of-practice. For this and several other reasons, RDEIM differs from REAcct.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

Effective thermal property improves phase change paint data

The phase-change coating technique presents itself as a valuable tool in determining the heat transfer rate over the surface of small complex wind tunnel models. A numerical technique is described which shows that an effective thermophysical property - the square root of the product of thermal conductivity, density, and specific heat - may significantly improve the accuracy of the phase-change coating technique with allowance for model inhomogeneity and temperature dependency in a transient environment. Results of the measured steady-state variation of the effective thermophysical property with temperature and the effect of surface heating rate on the effective thermophysical property are plotted for a representative homogeneous model material and for an extreme nonhomogeneous model material. The use of an effective thermophysical property to reduce phase-change paint data is recommended. The analysis also confirms that the apparatus described by Corwin and Kramer (1975) can be used to measure directly this effective thermophysical property for use in wind tunnel model heat-transfer measurements.

Drummond, J. P.↗

Effect of a simulated ice accretion on the aerodynamics of a swept wing

The effect of a simulated glaze ice accretion on the aerodynamic performance of a three-dimensional swept wing is studied experimentally. A semispan wing of effective aspect ratio four was mounted from the sidewall of the UIUC subsonic wind tunnel. The model uses a NACA 0012 airfoil section on a rectangular planform with interchangeable tip and root sections to allow for 0- and 30-degree sweep. A sidewall suction system is used to minmize the tunnel boundary-layer interaction with the model. Surface pressure data from five spanwise stations are compared to earlier data from a similar tunnel. A three-component sidewall balance has been designed, built and used to measure lift, drag and pitching moment on the clean and iced model. The data compare well to the integrated pressure data and to theory on the clean model. In addition, helium-bubble flow visualization has been performed on the iced model and reveals extensive spanwise flow in the separation bubble aft of the upper surface horn. This compares well to the computational results of other researchers. Sidewall suction was found to have no effect on the aerodynamics of the swept wing.

Bragg, M. B.↗

Aerodynamic measurements on a finite wing with simulated ice

The effect of a simulated glaze ice accretion on the aerodynamic performance of a three-dimensional straight and swept wing is studied experimentally. A semispan wing of effective aspect ratio five was mounted from the sidewall of the UIUC subsonic wind tunnel. The model uses an NACA 0012 airfoil section on a rectangular planform with interchangeable tip and root sections to allow for 0- and 30-deg sweep. A sidewall suction system is used to minimize the tunnel boundary-layer interaction with the model. A three-component sidewall balance has been designed, built and used to measure lift, drag and pitching moment on the clean and iced model. Fluorescent oil flow visualization has been performed on the iced model and reveals extensive spanwise flow in the separation bubble aft of the upper surface horn. These results are compared to computational results for the surface pressures, span loads and surface oil flow.

Bragg, M. B.↗

Daily water stable isotopes, transpiration, and matrix potential data for an aspen and engelmann stand in the East River Watershed (version 2)

We provide daily stable isotope (2H & 18O) ratios in soil water and xylem (plant stem) water, as well as the sap flow (transpiration) and the soil's matric potential at a forested site near Gothic, Colorado, in the East River catchment. We measured the stable isotopic composition of the transpiration and the daily transpiration flux sum of three aspen and three engelmann spruce. In both forest stands, we installed a soil profile and measured the soil matric potential at 15, 30, and 60 cm depth as well as the stable isotopes of soil pore water at 5, 10, 30, 60, and 90 cm depths. All isotope measurements were done in situ via vapor probes connected to a cavity ring down spectrometer (Picarro L1240i).We further report the daily meteorological data observed at billy barr near our study site. We also provide for each tree the relative share of root water uptake derived from the isotope measurements via a Bayesian mixing model (MixSIAR).The daily data is provided as a time series in "Iso_MP_Sap_DataDaily_ESSDiveUpload.csv" and the units are provided in "dd.csv"; the location of the instrumented trees and soil profiles are given as latitude and longitude coordinates saved as CSV and KMZ files; and a file-level metadata (flmd.csv) file that lists each file contained in the dataset with associated metadata.The data was gathered to investigate the short-term changes of the water sources (i.e., variation of root water uptake from different soil depths) of the studied subalpine trees.Update 07/16/2025: The relative and absolute plant water uptake depths were grouped to ensure that the MixSIAR model was applied with endmembers that differed in their d2H value by at least 3 permill and at least 1 permill for d18O. Whenever the difference between observed d2H values for two or more probes at neighboring depths was less than 3 permill, we used the average value for the source water endmember. For days at which probes that were not next to each other measurements did not differ at least 3 permill, the average of all probes between these two depths was used as the water source endmember representing the depth range between these two probes.

54 ENVIRONMENTAL SCIENCES↗

A generalized framework for reduced-order modeling of a wind turbine wake

A reduced-order model for a wind turbine wake is sought from large eddy simulation data. Fluctuating velocity fields are combined in the correlation tensor to form the kernel of the proper orthogonal decomposition (POD). Proper orthogonal decomposition modes resulting from the decomposition represent the spatially coherent turbulence structures in the wind turbine wake; eigenvalues delineate the relative amount of turbulent kinetic energy associated with each mode. Back-projecting the POD modes onto the velocity snapshots produces dynamic coefficients that express the amplitude of each mode in time. A reduced-order model of the wind turbine wake (wakeROM) is defined through a series of polynomial parameters that quantify mode interaction and the evolution of each POD mode coefficients. The resulting system of ordinary differential equations models the wind turbine wake composed only of the large-scale turbulent dynamics identified by the POD. Tikhonov regularization is used to recalibrate the dynamical system by adding additional constraints to the minimization seeking polynomial parameters, reducing error in the modeled mode coefficients. The wakeROM is periodically reinitialized with new initial conditions found by relating the incoming turbulent velocity to the POD mode coefficients through a series of open-loop transfer functions. The wakeROM reproduces mode coefficients to within 25.2%, quantified through the normalized root-mean-square error. Furthermore, a high-level view of the modeling approach is provided as a platform to discuss promising research directions, alternate processes that could benefit stability and efficiency, and desired extensions of the wakeROM.

17 WIND ENERGY↗

Cytochalasin D does not inhibit gravitropism in roots

It is generally thought that sedimenting plastids are responsible for gravity sensing in higher plants. We directly tested the model generated by the current statolith hypothesis that the gravity sensing that leads to gravitropism results from an interaction between the plastids and actin microfilaments. We find that the primary roots of rice, corn, and cress undergo normal gravitropism and growth even when exposed to cytochalasin D, a disruptor of actin microfilaments. These results indicate that an interaction between amyloplasts and the actin cytoskeleton is not critical for gravity sensing in higher plants and weaken the current statolith hypothesis.

NASA Discipline Plant Biology↗

Grey-box and ANN-based building models for multistep-ahead prediction of indoor temperature to implement model predictive control

Model-based predictive control (MPC) strategies for heating, ventilation, and air-conditioning (HVAC) systems present an opportunity to lower building energy consumption and operational costs. Such approaches rely on the development of a model to precisely forecast building thermal dynamics, such as room air temperature or heating/cooling rate, and make control-related decisions. The control-oriented modeling of building energy systems should be accurate in predicting indoor conditions and present low computational complexity. These features are the key challenge of implementing advanced control methods such as MPC. Extant studies on building modeling for MPC have focused on step-ahead forecasting techniques to forecast building thermal dynamics, while multistep-ahead forecasting is essential. Moreover, machine learning model suitable in case of the domain-based engineering expertise are also not available. To this aim, we perform a comparative analysis of the grey-box model based on a resistance-capacitance (RC) thermal network and a machine learning model composed of an artificial neural network (ANN) for multistep-ahead prediction of building thermal dynamics using current and historical data. Actual experimental data obtained from the Flexible Research Platform (FRP) in Oak Ridge National Laboratory (US) are used for estimation and validation purposes. The average root mean squared error (RMSE) of the grey-box and ANN models are 0.89 °C and 1.02°C, respectively. Finally, the results indicate that the grey-box model outperforms the ANN model in the considered validation periods in terms of accuracy and prediction stability.

42 ENGINEERING↗

Air compressor load forecasting using artificial neural network

Air compressor systems are responsible for approximately 10% of the electricity consumed in United States and European Union industry. As many researches have proven the effectiveness of using Artificial Neural Network inair compressor performance prediction, there is still a need to forecast the air compressor electrical load profile. The objective of this study is to predict compressed air systems’ electrical load profile, which is valuable to industry practitioners as well as software providers in developing better practice and tools for load managementand look-ahead scheduling programs. Two artificial neural networks, Two-Layer Feed-Forward Neural Networkand Long Short-Term Memory were used to predict an air compressors electrical load. Compressors with three different control mechanisms are evaluated with a total number of 11,874 observations. Here, the forecasts were validated using out-of-sample datasets with 5-fold cross-validation. Models produced average coefficient ofdetermination values from 0.24 to 0.94, average root-mean-square errors from 0.05 kW - 5.83 kW, and meanabsolute scaled errors from 0.20 to 1.33. The results indicate that both artificial neural networks yield goodresults for compressors using variable speed drive (average R 2 = 0.8 and no naïve forecasting), only the longshort-term memory model gives acceptable results for compressors using on/off control (average R 2 = 0.82 and no naïve forecasting), and no satisfactory results are obtained for load/unload type air compressors (models constituting active forecasting).

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Tree hydrodynamic modelling of the soil–plant–atmosphere continuum using FETCH3

Modelling the water transport along the soil–plant–atmosphere continuum is fundamental to estimating and predicting transpiration fluxes. A Finite-difference Ecosystem-scale Tree Crown Hydrodynamics model (FETCH3) for the water fluxes across the soil–plant–atmosphere continuum is presented here. The model combines the water transport pathways into one vertical dimension, and assumes that the water flow through the soil, roots, and above-ground xylem can be approximated as flow in porous media. This results in a system of three partial differential equations, resembling the Richardson–Richards equation, describing the transport of water through the plant system and with additional terms representing sinks and sources for the transfer of water from the soil to the roots and from the leaves to the atmosphere. The numerical scheme, developed in Python 3, was tested against exact analytical solutions for steady state and transient conditions using simplified but realistic model parameterizations. The model was also used to simulate a previously published case study, where observed transpiration rates were available, to evaluate model performance. With the same model setup as the published case study, FETCH3 results were in agreement with observations. Through a rigorous coupling of soil, root xylem, and stem xylem, FETCH3 can account for variable water capacitance, while conserving mass and the continuity of the water potential between these three layers. FETCH3 provides a ready-to-use open access numerical model for the simulation of water fluxes across the soil–plant–atmosphere continuum.

54 ENVIRONMENTAL SCIENCES↗

NASA's Soil Moisture Active and Passive (SMAP) Mission

The Soil Moisture Active-Passive (SMAP) Mission is one of the first Earth observation satellites being formulated by NASA in response to the 2007 National Research Council s Decadal Survey. SMAP will make global measurements of soil moisture at the Earth's land surface and its freeze-thaw state. These measurements will allow significantly improved estimates of water, energy and carbon transfers between the land and atmosphere. Soil moisture measurements are also of great importance in assessing flooding and monitoring drought. Knowledge gained from SMAP observations can help mitigate these natural hazards, resulting in potentially great economic and social benefits. SMAP observations of soil moisture and freeze/thaw timing over the boreal latitudes will also reduce a major uncertainty in quantifying the global carbon balance and help to resolve an apparent missing carbon sink over land. The SMAP mission concept will utilize an L-band radar and radiometer sharing a rotating 6-meter mesh reflector antenna flying in a 680 km polar orbit with an 8-day exact ground track repeat aboard a 3-axis stabilized spacecraft to provide high-resolution and high-accuracy global maps of soil moisture and freeze/thaw state every two to three days. In addition, the SMAP project will use these surface observations with advanced modeling and data assimilation to provide estimates of deeper root-zone soil moisture and net ecosystem exchange of carbon. SMAP recently completed its Phase A Mission Concept Study Phase for NASA and transitioned into Phase B (Formulation and Detailed Design). A number of significant accomplishments occurred during this initial phase of mission development. The SMAP project held several open meetings to solicit community feedback on possible science algorithms, prepared preliminary draft Algorithm Theoretical Basis Documents (ATBDs) for each mission science product, and established a prototype algorithm testbed to enable testing and evaluation of the performance of candidate algorithms. SMAP conducted an Applications Workshop in September 2009 to coordinate with potential application users interested in the mission data. A draft Applications Plan describing the Project s planned outreach to potential applications users has been prepared and will be updated during Phase B. SMAP made a significant evaluation of the potential terrestrial radio frequency interference (RFI) source environment and established radiometer and radar flight hardware and ground processing mitigation approaches. SMAP finalized its science orbit and orbit injection approach to optimize launch mass and prepared launch and commissioning scenarios and timeline. A science data communications approach was developed to maximize available science data volume to improve science margins while maintaining moderately short data product latencies to support many potential applications using existing ground assets and with minimum impact to the flight system. SMAP developed rigid multi-body and flexible body dynamics and control models and system designs for the 6-meter rotating instrument reflector-boom assembly (RBA) and flight system to confirm pointing and control performance, and devised strategies to efficiently implement on-orbit balancing if needed. Industry partners were selected for the spin mechanism assembly (SMA) and RBA. Preliminary designs for the radar and radiometer were initiated, including constructing breadboards of key assemblies.

Kellogg, Kent↗

Initial geomagnetic field model from Magsat vector data

Magsat data from the magnetically quiet days of November 5-6, 1979, were used to derive a thirteenth degree and order spherical harmonic geomagnetic field model, MGST(6/80). The model utilized both scalar and high-accuracy vector data and fit that data with root-mean-square deviations of 8.2, 6.9, 7.6 and 7.4 nT for the scalar magnitude, B(r), B(theta), and B(phi), respectively. The model includes the three first-order coefficients of the external field. Comparison with averaged Dst indicates that zero Dst corresponds with 25 nT of horizontal field from external sources. When compared with earlier models, the earth's dipole moment continues to decrease at a rate of about 26 nT/yr. Evaluation of earlier models with Magsat data shows that the scalar field at the Magsat epoch is best predicted by the POGO(2/72) model but that the WC80, AWC/75 and IGS/75 are better for predicting vector fields.

Langel, R. A.↗

Dynamic sea surface topography, gravity and improved orbit accuracies from the direct evaluation of SEASAT altimeter data

A method for the simultaneous solution of dynamic ocean topography, gravity and orbits using satellite altimeter data is described. A GEM-T1 based gravitational model called PGS-3337 that incorporates Seasat altimetry, surface gravimetry and satellite tracking data has been determined complete to degree and order 50. The altimeter data is utilized as a dynamic observation of the satellite's height above the sea surface with a degree 10 model of dynamic topography being recovered simultaneously with the orbit parameters, gravity and tidal terms in this model. PGS-3337 has a geoid uncertainty of 60 cm root-mean-square (RMS) globally, with the uncertainty over the altimeter tracked ocean being in the 25 cm range. Doppler determined orbits for Seasat, show large improvements, with the sub-30 cm radial accuracies being achieved. When altimeter data is used in orbit determination, radial orbital accuracies of 20 cm are achieved. The RMS of fit to the altimeter data directly gives 30 cm fits for Seasat when using PGS-3337 and its geoid and dynamic topography model. This performance level is two to three times better than that achieved with earlier Goddard earth models (GEM) using the dynamic topography from long-term oceanographic averages. The recovered dynamic topography reveals the global long wavelength circulation of the oceans with a resolution of 1500 km. The power in the dynamic topography recovery is now found to be closer to that of oceanographic studies than for previous satellite solutions. This is attributed primarily to the improved modeling of the geoid which has occurred. Study of the altimeter residuals reveals regions where tidal models are poor and sea state effects are major limitations.

Marsh, J. G.↗

Comparison of Wind-Tunnel and Flight Measurements of Stability and Control Characteristics of a Douglas A-26 Airplane

Stability and control characteristics determined from tests in the Langley 19-foot pressure tunnel of a 0.2375-scale model of the Douglas XA-26 airplane are compared with those measured in flight tests of a Douglas A-26 airplane. Agreement regarding static longitudinal stability as indicated by the elevator-fixed neutral points and by the variation of elevator deflection in both straight and turning flight was found to be good except at speeds approaching the stall. At these low speeds the airplane possessed noticeably improved stability, which was attributed to pronounced stalling at the root of the production wing. The pronounced root stalling did not occur on the smooth, well-faired model wing. Elevator tab effectiveness determined from model tests agreed well with flight-test tab effectiveness, but control-force variations with speed and acceleration were not in good agreement. The use of model hinge-moment data obtained at zero sideslip appeared to be satisfactory for the determination of aileron forces in sideslip. Fairly good correlation in aileron effectiveness and control forces was obtained; fabric distortion may have been responsible to some extent for higher flight values of aileron force at high speeds. Estimation of sideslip developed in an abrupt aileron roll was fair, but determination of the rudder deflection required to maintain zero sideslip in a rapid aileron roll was not entirely satisfactory.

Kayten, Gerald G↗

Acceleration and transport of galactic cosmic rays

A recently proposed cosmic-ray acceleration mechanism is combined with a model of transport and loss from the Galaxy to construct a model for galactic cosmic rays. Various aspects of the cosmic-ray spectrum, including the secondary-to-primary ratio as a function of energy, are found to fit naturally into the overall picture. A solution to the dynamical-halo model with energy-dependent diffusion coefficient (kappa) is presented. This model is consistent with the observed spectra if kappa is proportional to the square root of the kinetic energy. The observed flattening of the secondary-to-primary ratio below about 1-2 GeV per nucleon is shown to be a natural consequence of the model. Similarly, the primary electron spectrum is predicted to flatten below about 1-2 GeV, which is consistent with deconvolution of the local nonthermal radio spectrum, although the flattening is not quite as much as indicated by the observations.

Jokipii, J. R.↗

Multimodal correlative imaging and modelling of phosphorus uptake from soil by hyphae of mycorrhizal fungi

Phosphorus (P) is essential for plant growth. Arbuscular mycorrhizal fungi (AMF) aid its uptake by acquiring P from sources distant from roots in return for carbon. Little is known about how AMF colonise soil pore-space, and models of AMF-enhanced P-uptake are poorly validated. We used synchrotron X-ray computed tomography to visualize mycorrhizas in soil and synchrotron X-ray fluorescence/X-ray absorption near edge structure (XRF/XANES) elemental mapping for P, sulphur (S) and aluminium (Al) in combination with modelling. We found that AMF inoculation had a suppressive effect on colonisation by other soil fungi and identified differences in structure and growth rate between hyphae of AMF and nonmycorrhizal fungi. Our results showed that AMF co-locate with areas of high P and low Al, and preferentially associate with organic-type P species over Al-rich inorganic P. We discovered that AMF avoid Al-rich areas as a source of P. Sulphur-rich regions were found to be correlated with higher hyphal density and an increased organic-associated P-pool, whilst oxidized S-species were found close to AMF hyphae. Increased S oxidation close to AMF suggested the observed changes were microbiome-related. Our experimentally-validated model led to an estimate of P-uptake by AMF hyphae that is an order of magnitude lower than rates previously estimated – a result with significant implications for the modelling of plant–soil–AMF interactions.

79 ASTRONOMY AND ASTROPHYSICS↗

Convergent Hydraulic Redistribution and Groundwater Access Supported Facilitative Dependency Between Trees and Grasses in a Semi‐Arid Environment

Abstract Hydraulic redistribution is the transport of water from wet to dry soil layers, upward or downward, through plant roots. Often in savanna and woodland ecosystems, deep‐rooted trees, and shallow‐rooted grasses coexist. The degree to which these different species compete for or share soil‐water derived from precipitation or groundwater, as well as how these interactions are altered by hydraulic redistribution, is unknown. We use a multilayer canopy model and field observations to examine how the presence of deep, but tree‐root accessible, groundwater impacts seasonal patterns of hydraulic redistribution, and interaction between coexisting vegetation species in a semiarid riparian woodland (US‐CMW). Based on the simulation, trees absorb moisture at the water table (∼10 m depth) and release it in the shallow soil depth (0–3 m) during the dry pre‐monsoon season. We observed the occurrence of a new convergent hydraulic redistribution pattern during the monsoon season, where moisture is transported from both the near‐surface (0–0.5 m) and the water table to intermediate soil layers (1–5 m) through tree roots. We found that hydraulic redistribution demonstrates a growth facilitation effect at this site, supporting 49% of growing season tree transpiration and 14% of the grass transpiration. Compared to a similarly structured upland savanna without accessible groundwater, the riparian site shows an increased amount of hydraulically redistributed water and more facilitative water use between coexisting grasses and trees. These results shed light on the linkage between accessible groundwater and the role of hydraulic redistribution on the interaction between deep‐rooted and shallow‐rooted vegetation.

Lee, E.↗