Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “linear scaling”

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,567 records · Page 87
All sourcesNASA NTRSDOE OSTI
NASA NTRS1993

Study of atmospheric and bidirectional effects on surface reflectance and vegetation index time series: Application to NOAA AVHRR and preparation for future space missions

The objectives of the investigation, namely 'to characterize the atmospheric and directional effects on surface reflectance and vegetation index using the First International Satellite Cloud Climatology Project (ISLCSP) Field Experiment (FIFE) data set, develop new algorithms to obtain better Advanced Very High Resolution Radiometer (AVHRR) indices, and define possible improvements for future satellite missions', were addressed in three separate, yet complementary studies. First, it was shown, from theoretical calculations, that visible and near infrared reflectances combined linearly at optimum (one or two) viewing angles relate linearly to the fraction of photosynthetically available radiation absorbed by plants, f(sub par), can be used independently of the type of foliage and substrate, eliminate the effects of sub-pixel spatial heterogeneity, and improve the accuracy of the f(sub par) estimates when compared to the Normalized Difference Vegetation Index, NDVI. Second, it was demonstrated that NDVI, even though it is not a linear combination of radiances or reflectances, can be spatially integrated without significant loss of information from scales of 300 to 1000 m. Third, AVHRR visible and near-infrared reflectances over the FIFE site, separating temporal and bidirectional components and determining the model parameters through an original iterative scheme was successfully modeled. It appears that NDVI generated from the top-of-atmosphere reflectances normalized by the bidirectional effects (as determined in the scheme) is a better vegetation index than maximum NDVI. Details about the three studies are presented.

Frouin, Robert↗
NASA NTRS1993

Wave generation by contaminant ions near a large spacecraft

Measurements from the space shuttle flights have revealed that a large spacecraft in a low earth orbit is accompanied by an extensive gas cloud which is primarily made up of water. The charge exchange between the water molecule and the ionospheric O(+) ions produces a water ion beam traversing downstream of the spacecraft. In this report we present results from a study on the generation of plasma waves by the interaction of the water ion beams with the ionospheric plasma. Since velocity distribution function is key to the understanding of the wave generation process, we have performed a test particle simulation to determine the nature of H2O(+) ions velocity distribution function. The simulations show that at the time scales shorter than the ion cyclotron period tau(sub c), the distribution function can be described by a beam. On the other hand, when the time scales are larger than tau(sub c), a ring distribution forms. A brief description of the linear instabilities driven by an ion beam streaming across a magnetic field in a plasma is presented. We have identified two types of instabilities occurring in low and high frequency bands; the low-frequency instability occurs over the frequency band from zero to about the lower hybrid frequency for a sufficiently low beam density. As the beam density increases, the linear instability occurs at decreasing frequencies below the lower-hybrid frequency. The high frequency instability occurs near the electron cyclotron frequency and its harmonics.

Singh, N.↗
NASA NTRS1993

The heating of suprathermal ions above thunderstorm cells

We estimate ion heating in the topside ionosphere directly over thunderstorm cells. The primary heating is due to lower hybrid waves excited through linear mode coupling as intense electromagnetic (EM) whistler mode radiation from lightning is scattered from small scale (2 - 20 m) magnetic-field-aligned plasma density irregularities in the topside ionosphere. For typical radiated EM fields, we find that suprathermal H+ ions in the 6 eV and greater energy range can be heated by 20 to 40 eV as a result of a single lightning discharge. We also show how the number density of 6 eV and greater H(+) ions is enhanced by preheating resulting from the absorption of proton whistlers in the 500-1000 km altitude range. For lightning discharge rates of one or more per second over a 10 exp 4 sq km area, our model predicts a total energy gain for the H(+) ions of 400 eV to 2 KeV and a perpendicular ion flux of about 10 exp 5 to 10 exp 6/sq cm sec. These fluxes should be observable on low altitude spacecraft using presently available instrumentation.

Bell, T. F.↗
NASA NTRS1994

Equatorial wave activity derived from fluctuations in observed convection

The spectrum of equatorial wave activity propagating vertically into the stratosphere is calculated from high-resolution imagery of the global convective pattern. Synoptic Global Cloud Imagery (GCI), constructed from six satellites simultaneously observing the earth, is used to diabatically force the linearized primitive equations. Having resolution of 0.5 deg and 3 h, that imagery captures the dominant scales of organized convection, including several harmonics of the diurnal cycle. Its global coverage with high space-time resolution allows the GCI to represent heating variability and dynamical behavior excited by it over a wide range of scales. The dynamical response above the heating is evaluated globally in terms of a space-time spectrum of Hough modes, one which includes planetary-scale Kelvin waves, Rossby waves, and gravity waves down to the resolution of the GCI. The geopotential response, which is indicative of temperature fluctuations observed by satellite, is very red in frequency. Therefore, planetary-scale waves with periods longer than two days dominate the spectrum of geopotential, while high-frequency gravity waves make a comparatively small contribution. Some 80% of the geopotential variance is accounted for by the Kelvin and gravest-symmetric Rossby modes, while the Rossby-gravity mode is comparatively weak. In horizontal eddy motion, the excited wave spectrum is still dominated by planetary-scale components. However, meridional wind fluctuations associated with the Rossby-gravity mode have variance comparable to that of zonal wind fluctuations associated with the Kelvin mode, even though the Rossby-gravity mode is nearly invisible in the geopotential response. Estimates of tropospheric heating lead to amplitudes and propagation characteristics that are broadly consistent with satellite and radiosonde observations of wave activity in the lower stratosphere. The space-time spectrum of EP flux is significantly whiter than the response in either geopotential or motion. Gravity waves of small scale and high frequency carry a large fraction of the upward flux. Although it dominates eastward variance of geopotential and motion, the Kelvin mode carries only about 50% of the eastward EP flux at phase speeds of 20-40 m/s and only 35% of the total eastward flux transmitted to the stratosphere. The remainder is carried by the gravity wave spectrum, which carries nearly all of the westward flux at phase speeds greater than 20 m/s. The gravity wave spectrum also contributes significantly at phase speeds of 10-20 m/s, where only 25% of the flux is accounted for by zonal wavenumbers less than 20. The broad nature of the gravity wave spectrum suggests its absorption at critical levels will be distributed over a deep layer of the middle atmosphere.

Bergman, John W.↗
NASA NTRS1998

The Oxidation of AlN in Dry and Wet Oxygen

The oxidation kinetics of AlN containing 3.5 wt% Y2O3 were studied by thermogravimetric analysis in dry oxygen and 10% H2O/balance oxygen at temperatures between 1000 and 1200 C for times between 48 and 100 h. The oxidation kinetics for AlN in dry oxygen were parabolic and of approximately the same magnitude and temperature dependence as other alumina forming materials. In this case, diffusion of oxygen and/or aluminum through the alumina scale is the rate limiting mechanism. The oxidation kinetics for AlN in wet oxygen were nearly linear and much more rapid than rates observed in dry oxygen. Numerous micropores were observed in the alumina formed on AIN in wet oxygen. These pores provide a fast path for oxygen transport. The linear kinetics observed in this case suggest that the interface reaction rate of AlN with wet oxygen is the oxidation rate limiting step.

Opila, Elizabeth↗
NASA NTRS2001

Structure of the Jovian Magnetodisk Current Sheet: Initial Galileo Observations

The ten-degree tilt of the Jovian magnetic dipole causes the magnetic equator to move back and forth across Jupiter's rotational equator and tile Galileo orbit that lies therein. Beyond about 24 Jovian radii, the equatorial current sheet thins and tile magnetic structure changes from quasi-dipolar into magnetodisk-like with two regions of nearly radial but antiparallel magnetic field separated by a strong current layer. The magnetic field at the center of the current sheet is very weak in this region. Herein we examine tile current sheet at radial distances from 24 55 Jovian radii. We find that the magnetic structure very much resembles tile structure seen at planetary magnetopause and tail current sheet crossings. Tile magnetic field variation is mainly linear with little rotation of the field direction, At times there is almost no small-scale structure present and the normal component of the magnetic field is almost constant through the current sheet. At other times there are strong small-scale structures present in both the southward and northward directions. This small-scale structure appears to grow with radial distance and may provide the seeds for tile explosive reconnection observed at even greater radial distances oil tile nightside. Beyond about 40 Jovian radii, the thin current sheet also appears to be almost constantly in oscillatory motion with periods of about 10 min. The amplitude of these oscillations also appears to grow with radial distance. The source of these fluctuations may be dynamical events in tile more distant magnetodisk.

Russell, C. T.↗
NASA NTRS1999

Practical Aspects of Stabilized FEM Discretizations of Nonlinear Conservation Law Systems with Convex Extension

This talk considers simplified finite element discretization techniques for first-order systems of conservation laws equipped with a convex (entropy) extension. Using newly developed techniques in entropy symmetrization theory, simplified forms of the Galerkin least-squares (GLS) and the discontinuous Galerkin (DG) finite element method have been developed and analyzed. The use of symmetrization variables yields numerical schemes which inherit global entropy stability properties of the POE system. Central to the development of the simplified GLS and DG methods is the Degenerative Scaling Theorem which characterizes right symmetrizes of an arbitrary first-order hyperbolic system in terms of scaled eigenvectors of the corresponding flux Jacobean matrices. A constructive proof is provided for the Eigenvalue Scaling Theorem with detailed consideration given to the Euler, Navier-Stokes, and magnetohydrodynamic (MHD) equations. Linear and nonlinear energy stability is proven for the simplified GLS and DG methods. Spatial convergence properties of the simplified GLS and DO methods are numerical evaluated via the computation of Ringleb flow on a sequence of successively refined triangulations. Finally, we consider a posteriori error estimates for the GLS and DG demoralization assuming error functionals related to the integrated lift and drag of a body. Sample calculations in 20 are shown to validate the theory and implementation.

Barth, Timothy↗
NASA NTRS2013

On the Star Formation-AGN Connection at zeta (is) approximately greater than 0.3

Using the spectra of a sample of approximately 28,000 nearby obscured active galaxies from Data Release 7 of the Sloan Digital Sky Survey (SDSS), we probe the connection between active galactic nucleus (AGN) activity and star formation over a range of radial scales in the host galaxy. We use the extinction-corrected luminosity of the [O iii] 5007A line as a proxy of intrinsic AGN power and supermassive black hole (SMBH) accretion rate. The star formation rates (SFRs) are taken from the MPA-JHU value-added catalog and are measured through the 3 inch SDSS aperture. We construct matched samples of galaxies covering a range in redshifts. With increasing redshift, the projected aperture size encompasses increasing amounts of the host galaxy. This allows us to trace the radial distribution of star formation as a function of AGN luminosity. We find that the star formation becomes more centrally concentrated with increasing AGN luminosity and Eddington ratio. This implies that such circumnuclear star formation is associated with AGN activity, and that it increasingly dominates over omnipresent disk star formation at higher AGN luminosities, placing critical constraints on theoretical models that link host galaxy star formation and SMBH fueling. We parameterize this relationship and find that the star formation on radial scales (is) less than 1.7 kpc, when including a constant disk component, has a sub-linear dependence on SMBH accretion rate: SFR in proportion to solar mass(sup 0.36), suggesting that angular momentum transfer through the disk limits accretion efficiency rather than the supply from stellar mass loss.

LaMassa, Stephanie M.↗
NASA NTRS2015

Evaluation of Long-Term Calibrations of the AVHRR Visible Radiances

Two systematic calibrations have been compiled for the visible radiances measured by the series of AVHRR instruments flown on the NOAA operational polar weather satellites: one by the International Satellite Cloud Climatology Project (ISCCP), anchored on NASA ER-2 underflights in the 1980s and early 1990s and covering the period 1981-2009, and one by the PATMOS-x project, anchored on comparisons to the MODIS instruments on the Aqua and Terra satellites in the 2000s and covering the period 1979-2010 (this result also includes calibration for the near-IR channels). Both methods have had to extend their anchor calibrations over a long series of instruments using different vicarious approaches, so a comparison provides an opportunity to evaluate how well this extension works by cross-checking the results at the anchor points. The basic result of this comparison is that for the ''afternoon'' series of AVHRRs, the calibrations agree to within their mutual uncertainties. However, this retrospective evaluation also shows that the representation of the time variations can be simplified. The ISCCP procedure had much more difficulty extending the calibration to the ''morning'' series of AVHRRs with the calibrations for NOAA-15 and NOAA-17 exceeding the estimated uncertainties. Given the general agreement, a new calibration for all AVHRR visible radiances (except TIROS-N, NOAA-6, NOAA-19, and MetOp-A) is proposed that is based on the average of the best linear fits to the two time records. The estimated uncertainty of these calibrations is 63% absolute (scaled radiance units).

Advanced Very High Resolution Radiometer↗
NASA NTRS2017

Aerodynamic Coefficients from Aeroballistic Range Testing of Deployed- and Stowed-SIAD SFDT Models

This report documents a ballistic-range test campaign conducted in 2012 in order to estimate the aerodynamic stability characteristics of two configurations of the Supersonic Flight Dynamics Test (SFDT) vehicle prior to its initial flight in 2014. The SFDT vehicle was a test bed for demonstrating several new aerodynamic decelerator technologies then being developed under the Low-Density Supersonic Decelerator (LDSD) Project. Of particular interest here is the Supersonic Inflatable Aerodynamic Decelerator (SIAD), an inflatable attached torus used to increase the drag surface area of an entry vehicle during the supersonic portion of the entry trajectory. Two model configurations were tested in the ballistic range: one representing the SFDT vehicle prior to deployment of the SIAD, and the other representing the nominal shape with the SIAD inflated. Both models were fabricated from solid metal, and therefore, the effects of the flexibility of the inflatable decelerator were not considered. The test conditions were chosen to match, as close as possible, the Mach number, Reynolds number, and motion dynamics expected for the SFDT vehicle in flight, both with the SIAD stowed and deployed. For SFDT models with the SIAD stowed, 12 shots were performed covering a Mach number range of 3.2 to 3.7. For models representing the deployed SIAD, 37 shots were performed over a Mach number range of 2.0 to 3.8. Pitch oscillation amplitudes covered a range from 0.7 to 20.6 degrees RMS. Portions of this report (data analysis approach, aerodynamic modeling, and resulting aerodynamic coefficients) were originally published as an internal LDSD Project report [1] in 2012. In addition, this report provides a description of the test design approach, the test facility, and experimental procedures. Estimated non-linear aerodynamic coefficients, including pitch damping, for both model configurations are reported, and the shot-by-shot trajectory measurements, plotted in comparison with calculated trajectories based on the derived non-linear aerodynamic coefficients, are provided as appendices. Since the completion of these tests, two full-scale SFDT flights have been successfully conducted: one in June 2014 [2, 3], and one in June 2015 [3].

Ground Test↗
NASA NTRSDate not supplied

ISRU Pilot Excavator: Bucket Drum Scaling Experimental Results

NASA’s Space Technology Mission Directorate (STMD) is funding the development of a robotic excavator called the “ISRU Pilot Excavator” which will be a technology demonstration of excavating and transporting 10 metric tons of lunar regolith on the surface of the moon with a 30k-class robotic excavator. ISRU Pilot Excavator will be the next generation of robotic excavator to use bucket drums as excavation tools. This is an evolution of the Regolith Advanced Surface Systems Operations Robot (RASSOR) developed at NASA’s Kennedy Space Center. Bucket drums are hollow cylinders with regularly spaced scoops around the perimeter. The drums rotate in one direction to collect regolith with the scoops. The regolith slides down an internal baffling system inside the drum which prevents the regolith from falling back out of the scoops. The captured regolith can then be transported while held in the drum and then deposited by rotating the drum in the opposite direction allowing the regolith to slide back down the baffling and out of the excavation scoops. Bucket drums were developed by Lockheed Martin in 2008 and used on multiple robotic excavator prototypes ever since. However the forces on a bucket drum and considerations for scaling have not been measured in detail. Bucket drums are challenging to model using classical blade\bucket equations because of their unique geometry. Therefore this experiment was performed to measure the forces on three bucket drums of the same geometry at different scales. Small: 9.4” (239mm) dia. x 8.1” (206mm) width, Medium: 11.6” (294mm) dia. x 10” (254mm) width, and Large: 17” (432mm) dia. x 14.1”(358mm) width. The test stand consisted of an actuated gantry with controlled motion in the vertical (Z) and horizontal (X) axes and a single rotation axis (R). The bucket drums were individually mounted to the rotary axis of the test stand and translated across a prepared bed of BP-1 lunar regolith simulant at a specified linear speed and cutting depth. The test stand was outfitted with a torque sensor in line with the rotation of the drum (R) and a 3 axis (X, Y, and Z) load cell. In addition to the three sizes of bucket drums the linear excavation speed and cutting depth were test variables. The results of these experiments show the relationship between the three scales of bucket drums for factors such as: excavation force, torque due to regolith rotation inside the drum, excavation energy, time to fill, etc. and will be discussed in detail in this paper. This fundamental data will be used in the design of the ISRU Pilot Excavator and can inform the design of future bucket drum excavators.

Jason Michael Schuler↗
NASA NTRSDate not supplied

ISRU Pilot Excavator: Bucket Drum Scaling Experimental Results

NASA’s Space Technology Mission Directorate (STMD) is funding the development of a robotic excavator called the “ISRU Pilot Excavator” (IPEx) which will be a technology demonstration of excavating and transporting 10 metric tons of lunar regolith on the surface of the moon with a 30kg-class robotic excavator. IPEx will be the next generation of robotic excavators to use bucket drums as excavation tools. This is an evolution of the Regolith Advanced Surface Systems Operations Robot (RASSOR) developed at NASA’s Kennedy Space Center (KSC). Bucket drums are hollow cylinders with regularly spaced scoops around the perimeter. The drums rotate in one direction to collect regolith with the scoops. The regolith slides down an internal baffling system inside the drum which prevents the regolith from falling back out of the scoops (see Figure 1). The captured regolith can then be transported while held in the drum and then deposited by rotating the drum in the opposite direction allowing the regolith to slide back down the baffling and out of the excavation scoops. Bucket drums were developed by Lockheed Martin in 2008 and used on multiple robotic excavator prototypes ever since. However, the forces on a bucket drum and considerations for scaling have not been measured in detail. Bucket drums are challenging to model using classical blade\bucket equations because of their unique geometry. Therefore, this experiment was performed to measure the forces on three bucket drums of the same geometry at different scales. Small: 9.4” (239mm) dia. x 8.1” (206mm) width, Medium: 11.6” (294mm) dia. x 10” (254mm) width, and Large: 17” (432mm) dia. x 14.1” (358mm) width. The test stand consisted of an actuated gantry with controlled motion in the vertical (Z) and horizontal (X) axes and a single rotation axis (R). The bucket drums were individually mounted to the rotary axis of the test stand and translated across a prepared bed of BP-1 lunar regolith simulant at a specified linear speed and cutting depth. The test stand was outfitted with a torque sensor in line with the rotation of the drum (R) and a 3 axis (X, Y, and Z) load cell. In addition to the three sizes of bucket drums the linear excavation speed and cutting depth were test variables. The results of these experiments show the relationship between the three scales of bucket drums for factors such as: excavation force, torque due to regolith rotation inside the drum, excavation energy, time to fill, etc. and will be discussed in detail in this paper. This fundamental data will be used in the design of IPEx and can inform the design of future bucket drum excavators.

RASSOR↗
NASA NTRS2024

COLUTA: Custom 8-Channel 15-bit 40-MSPS ADC for the ATLAS Liquid Argon Calorimeter Readout

The stringent radiation-tolerance, quantization requirements, and the need for seamless integration of the on-detector readout electronics chain require the design of a full-custom analog-to-digital converter (ADC) for the 182,468 channels in the ATLAS Liquid Argon (LAr) calorimeter at the Large Hadron Collider (LHC) at CERN, as part of the High-Luminosity LHC (HL-LHC) upgrade. Each of the 8 channels in the prototype 65 nm CMOS ADC with 15-bit resolution and >68 dB signal-to-noise-and-distortion ratio (SNDR), or equivalently >11 effective number of bits (ENOB), consists of a Multiplying-DAC (MDAC) and a successive-approximation (SAR) ADC. A digital data processing unit (DDPU) calibrates and serially transmits the quantized data. Initial lab characterization of 18 chips shows a minimum SNDR of 69.5 dB at full-scale at about 5 MHz and a maximum power consumption of 1.17 W; differential non-linearity (DNL) measurements show no missing codes. Additional validation of the ADC is ongoing before approximately 70k chips can be mass-produced for the upgrade.

Rui Xu↗
NASA NTRS1971

On nonlinear effects in fracture mechanics.

Linear elastic treatment of fracture is considered applicable for net section stress up to about 0.8 the uniaxial tensile yield stress. Crack front plastic yield is still small enough to be viewed and treated as a small perturbation to the local crack front elastic stress field. Assuming these same circumstances and adopting the same point of view, an approach is presented for incorporating the nonlinear effects of small scale crack front plastic yield and slow crack extension in determination of the energy release rate and fracture toughness. Deviation from linearity of the load-displacement record in a fracture toughness test offers a quantifiable measure of these effects and is used to calculate the energy release rate. Fracture toughness values for one-eight inch thick 7075-T6 center cracked aluminum sheet are compared with uncorrected values and with values obtained by the Irwin method of plasticity correction.

Liebowitz, H.↗
DOE OSTI2024

Boson sampling with Gaussian input states: Toward efficient scaling and certification

A universal quantum computer of large scale is not available yet, however, intermediate models of quantum computation would still permit demonstrations of a quantum computational advantage over classical computing and could challenge the Extended Church-Turing Thesis. One of these models based on single photons interacting via linear optics is called Boson Sampling. Although Boson Sampling was demonstrated and the threshold to claim quantum computational advantage was achieved, the question of how to scale up Boson Sampling experiments remains. To make progress with this problem, here we present a practically achievable pathway to scale Boson Sampling experiments by combining continuous-variable quantum information and temporal encoding. Here, we propose the combination of switchable dual-homodyne and single-photon detections, the temporal loop technique, and scattershot-based Boson Sampling. We detail the required assumptions for concluding computational hardness for this configuration. Furthermore, this particular combination of techniques moves towards an efficient scaling and certification of Boson Sampling, all in a single experimental setup.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Early Exploration of a Flexible Framework for Efficient Quantum Linear Solvers in Power Systems

The rapid integration of renewable energy resources presents formidable challenges in managing power grids. While advanced computing and machine learning techniques offer some solutions for accelerating grid modeling and simulation, there remain complex problems that classical computers cannot effectively address. Quantum computing, a promising technology, has the potential to fundamentally transform how we manage power systems, especially in scenarios with a higher proportion of renewable energy sources. One critical aspect is solving linear systems of equations, crucial for power system applications like power flow analysis, for which the Harrow-Hassidim-Lloyd (HHL) algorithm is a well-known quantum solution. However, HHL quantum circuits often exhibit excessive depth, making them impractical for current Noisy-Intermediate-Scale-Quantum (NISQ) devices. In this paper, we introduce a versatile framework, powered by NWQSim, that bridges the gap between power system applications and quantum linear solvers available in Qiskit. This framework empowers researchers to efficiently explore power system applications using quantum linear solvers. Through innovative gate fusion strategies, reduced circuit depth, and GPU acceleration, our simulator significantly enhances resource efficiency. Power flow case studies have demonstrated up to a eight-fold speedup compared to Qiskit Aer, all while maintaining comparable levels of accuracy.

quantum computing, Harrow-Hassidim-Lloyd, high-per↗
DOE OSTI2020

Small-Scale Recycling of Irradiated Uranium and Transuranic Elements (Np, Pu) Using 3D-Printed Centrifugal Contactors

Niowave’s technology employs super-conducting linear accelerators to induce fission on low-enriched uranium targets followed by subsequent radiochemical processing to purify 99 Mo as well as a number of other fission products. Argonne’s Radiochemistry group (CFCT), in partnership with Niowave and UNLV, will be managing the TCF project to develop a closed-cycle loop for Niowave’s uranium targets. The roles of the DOE national laboratory partner will be to 1) develop a basic chemical understanding of the separations and purifications required to meet industry standards, 2) initiate the additive manufacturing (AM) of the centrifugal contactors, 3) develop a process flowsheet using Argonne computer codes.

43 PARTICLE ACCELERATORS↗
DOE OSTI2021

Non-linear hydrologic organization

We revisit three variants of the well-known Stommel diagrams that have been used to summarize knowledge of characteristic scales in time and space of some important hydrologic phenomena and modified these diagrams focusing on spatiotemporal scaling analyses of the underlying hydrologic processes. In the present paper we focus on soil formation, vegetation growth, and drainage network organization. We use existing scaling relationships for vegetation growth and soil formation, both of which refer to the same fundamental length and timescales defining flow rates at the pore scale but different powers of the power law relating time and space. The principle of a hierarchical organization of optimal subsurface flow paths could underlie both root lateral spread (RLS) of vegetation and drainage basin organization. To assess the applicability of scaling, and to extend the Stommel diagrams, data for soil depth, vegetation root lateral spread, and drainage basin length have been accessed. The new data considered here include timescales out to 150 Myr that correspond to depths of up to 240 m and horizontal length scales up to 6400 km and probe the limits of drainage basin development in time, depth, and horizontal extent.

58 GEOSCIENCES↗