Search NASA⌕ Search

SEARCH · Search NASA

Results for “root dynamics”

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.

486 records · Page 27
All sourcesNASA NTRSDOE OSTI
NASA NTRS2024

Effects of Water Limitation and Competition on Tree Carbon Allocation in an Earth System Modeling Framework

Earth system models (ESMs) have a limited capacity to represent plant functional diversity and shifts in trait distributions. Approaches to improving the representation of this complexity in ESMs include (i) optimality-based approaches that predict trait–environment responses and (ii) explicitly modelling coexistence and community assembly. These approaches are expected to converge only when optimality-based approaches identify competitively dominant strategies, which often differ from strategies that maximize ecosystem functioning or fitness components in monoculture. We used two models, LM3-PPA (a vegetation demographic model designed as an ESM component) and BiomeE (a computationally efficient analog for LM3-PPA), to explore how water limitation affects carbon allocation strategies of canopy trees. We compared competitive allocation strategies and those that maximize biomass or productivity in monoculture. We did not explicitly model coexistence or community assembly. Rather, we used model experiments to identify competitive and maximizing strategies in a two-dimensional trait space under different precipitation and mortality scenarios. At 10 eastern US locations, we simulated historical, wet and dry climate scenarios, novel drought and three different mortality scenarios (low, medium or high sensitivity to water deficit). For each site and scenario, we identified the competitive strategy and three maximizing strategies (maximum biomass, productivity or drought-tolerance). Root: leaf ratios tended to increase and leaf area tended to decrease with increasing water stress (increasing water limitation and its effects on mortality). However, relative to maximizing strategies, competitive strategies shifted towards greater allocation to roots and leaves with increasing water stress. Competitive overinvestments (greater allocation to roots and leaves by competitive strategies compared with maximizing strategies) were robust across different modelling contexts, including vegetation parameter sets (Acer vs. Populus), models (LM3-PPA vs. BiomeE) and uncalibrated vs. calibrated BiomeE versions. Synthesis: The theoretical prediction that competitive and maximizing allocation strategies differ under water limitation is confirmed for a demographic model designed as an ESM component. Optimality-based trait predictions can simplify representing trait diversity in ESMs but do not always correspond to competitive outcomes. Explicitly modelling coexistence and community assembly in ESMs is challenging but is likely the most general approach to representing trait diversity.

vegetation demographic model↗
NASA NTRS1947

Correlation of the Trim Limits of Stability Obtained for a PB2Y-3 Flying Boat and a 1/8-Size Powered Dynamic Model

Tests of a PB2Y-3 flying boat were made at the U.S> Naval Air Station, Patuxent River, Md., to determine its hydrodynamic trim limits of stability. Corresponding tests were also made of a 1/8-size powered dynamic model of the same flying boat in Langley tank no. 1. During the tank tests, the full-size testing procedure was reproduced as closely as possible in order to obtain data for a direct correlation of the results. As a nominal gross load of 66,000 pounds, the lower trim limits of the full-size and model were in good agreement above a speed of 80 feet per second. As the speed decreased below 80 feet per second, the difference between the model trim limits and full-scale trim limits gradually became larger. The upper trim limit of the model with flaps deflected 0 deg was higher than that of the full-size, but the difference was small over the speed range compared. At flap deflections greater than 0 deg, it was not possible to trim either the model of the airplane to the upper limit with the center of gravity at 28 percent of the mean aerodynamic chord. The decrease in the lower trim limits with increase in flap deflection showed good agreement for the airplane and model. The lower trim limits obtained at different gross loads for the full-size airplane were reduced to approximately a single curve by plotting trim against the square root of C(sub delta (sub o)) divided by C(sub V).

Garrison, Charlie C.↗
NASA NTRS2022

Detailing Cloud Property Feedbacks with a Regime-Based Decomposition

Diagnosing the root causes of cloud feedback in climate models and reasons for inter-model disagreement is a necessary first step in understanding their wide variation in climate sensitivities. Here we bring together two analysis techniques that illuminate complementary aspects of cloud feedback. The first quantifies feedbacks from changes in cloud amount, altitude, and optical depth, while the second separates feedbacks due to cloud property changes within specific cloud regimes from those due to regime occurrence frequency changes. We find that in the global mean, shortwave cloud feedback averaged across ten models comes solely from a positive within-regime cloud amount feedback countered slightly by a negative within-regime optical depth feedback. These within-regime feedbacks are highly uniform: In nearly all regimes, locations, and models, cloud amount decreases and cloud albedo increases with warming. In contrast, global-mean across-regime components vary widely across models but are very small on average. This component, however, is dominant in setting the geographic structure of the shortwave cloud feedback: Thicker, more extensive cloud types increase at the expense of thinner, less extensive cloud types in the extratropics, and vice versa at low latitudes. The prominent negative extratropical optical depth feedback has contributions from both within- and across-regime components, suggesting that thermodynamic processes affecting cloud properties as well as dynamical processes that favor thicker cloud regimes are important. The feedback breakdown presented herein may provide additional targets for observational constraints by isolating cloud property feedbacks within specific regimes without the obfuscating effects of changing dynamics that may differ across timescales.

climate sensitivity↗
NASA NTRS2000

Test Cases for a Rectangular Supercritical Wing Undergoing Pitching Oscillations

Steady and unsteady measured pressures for a Rectangular Supercritical Wing (RSW) undergoing pitching oscillations have been presented. From the several hundred compiled data points, 27 static and 36 pitching oscillation cases have been proposed for computational Test Cases to illustrate the trends with Mach number, reduced frequency, and angle of attack. The wing was designed to be a simple configuration for Computational Fluid Dynamics (CFD) comparisons. The wing had an unswept rectangular planform plus a tip of revolution, a panel aspect ratio of 2.0, a twelve per cent thick supercritical airfoil section, and no twist. The model was tested over a wide range of Mach numbers, from 0.27 to 0.90, corresponding to low subsonic flows up to strong transonic flows. The higher Mach numbers are well beyond the design Mach number such as might be required for flutter verification beyond cruise conditions. The pitching oscillations covered a broad range of reduced frequencies. Some early calculations for this wing are given for lifting pressure as calculated from a linear lifting surface program and from a transonic small perturbation program. The unsteady results were given primarily for a mild transonic condition at M = 0.70. For these cases the agreement with the data was only fair, possibly resulting from the omission of viscous effects. Supercritical airfoil sections are known to be sensitive to viscous effects (for example, one case cited). Calculations using a higher level code with the full potential equations have been presented for one of the same cases, and with the Euler equations. The agreement around the leading edge was improved, but overall the agreement was not completely satisfactory. Typically for low-aspect-ratio rectangular wings, transonic shock waves on the wing tend to sweep forward from root to tip such that there are strong three-dimensional effects. It might also be noted that for most of the test, the model was tested with free transition, but a few points were taken with an added transition strip for comparison. Some unpublished results of a rigid wing of the same airfoil and planform that was tested on the pitch and plunge apparatus mount system (PAPA) showed effects of the lower surface transition Strip on flutter at the lower subsonic Mach numbers. Significant effects of a transition strip were also obtained on a wing with a thicker supercritical section on the PAPA mount system. Both of these flutter tests on the PAPA resulted in very low reduced frequencies that may be a factor in this influence of the transition strip. However, these results indicate that correlation studies for RSW may require some attention to the estimation of transition location to accurately treat viscous effects. In this report several Test Cases are selected to illustrate trends for a variety of different conditions with emphasis on transonic flow effects. An overview of the model and tests is given and the standard formulary for these data is listed. Sample data points are presented in both tabular and graphical form. A complete tabulation and plotting of all the Test Cases is given. Only the static pressures and the real and imaginary parts of the first harmonic of the unsteady pressures are available. All the data for the test are available in electronic file form. The Test Cases are also available as separate electronic files.

Bennett, Robert M.↗
NASA NTRS2010

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↗
NASA NTRS1999

NASA/NOAA: Earth Science Electronic Theater 1999. Earth Science Observations, Analysis and Visualization: Roots in the 60s - Vision for the Next Millennium

The Etheater presents visualizations which span the period from the original Suomi/Hasler animations of the first ATS-1 GEO weather satellite images in 1966, to the latest 1999 NASA Earth Science Vision for the next 25 years. Hot off the SGI-Onyx Graphics-Supercomputer are NASA''s visualizations of Hurricanes Mitch, Georges, Fran and Linda. These storms have been recently featured on the covers of National Geographic, Time, Newsweek and Popular Science. Highlights will be shown from the NASA hurricane visualization resource video tape that has been used repeatedly this season on National and International network TV. Results will be presented from a new paper on automatic wind measurements in Hurricane Luis from 1-min GOES images that appeared in the November BAMS. The visualizations are produced by the NASA Goddard Visualization & Analysis Laboratory, and Scientific Visualization Studio, as well as other Goddard and NASA groups using NASA, NOAA, ESA, and NASDA Earth science datasets. Visualizations will be shown from the Earth Science ETheater 1999 recently presented in Tokyo, Paris, Munich, Sydney, Melbourne, Honolulu, Washington, New York, and Dallas. The presentation Jan 11-14 at the AMS meeting in Dallas used a 4-CPU SGI/CRAY Onyx Infinite Reality Super Graphics Workstation with 8 GB RAM and a Terabyte Disk at 3840 X 1024 resolution with triple synchronized BarcoReality 9200 projectors on a 60ft wide screen. Visualizations will also be featured from the new Earth Today Exhibit which was opened by Vice President Gore on July 2, 1998 at the Smithsonian Air & Space Museum in Washington, as well as those presented for possible use at the American Museum of Natural History (NYC), Disney EPCOT, and other venues. New methods are demonstrated for visualizing, interpreting, comparing, organizing and analyzing immense HyperImage remote sensing datasets and three dimensional numerical model results. We call the data from many new Earth sensing satellites, HyperImage datasets, because they have such high resolution in the spectral, temporal, spatial, and dynamic range domains. The traditional numerical spreadsheet paradigm has been extended to develop a scientific visualization approach for processing HyperImage datasets and 3D model results interactively. The advantages of extending the powerful spreadsheet style of computation to multiple sets of images and organizing image processing were demonstrated using the Distributed Image SpreadSheet (DISS). The DISS is being used as a high performance testbed Next Generation Internet (NGI) VisAnalysis of: 1) El Nino SSTs and NDVI response 2) Latest GOES 10 5-min rapid Scans of 26 day 5000 frame movie of March & April 198 weather and tornadic storms 3) TRMM rainfall and lightning 4)GOES 9 satellite images/winds and NOAA aircraft radar of hurricane Luis, 5) lightning detector data merged with GOES image sequences, 6) Japanese GMS, TRMM, & ADEOS data 7) Chinese FY2 data 8) Meteosat & ERS/ATSR data 9) synchronized manipulation of multiple 3D numerical model views; etc. will be illustrated. The Image SpreadSheet has been highly successful in producing Earth science visualizations for public outreach.

Hasler, A. Fritz↗
NASA NTRS2016

Comparison of Turbulent Thermal Diffusivity and Scalar Variance Models

In this study, several variable turbulent Prandtl number formulations are examined for boundary layers, pipe flow, and axisymmetric jets. The model formulations include simple algebraic relations between the thermal diffusivity and turbulent viscosity as well as more complex models that solve transport equations for the thermal variance and its dissipation rate. Results are compared with available data for wall heat transfer and profile measurements of mean temperature, the root-mean-square (RMS) fluctuating temperature, turbulent heat flux and turbulent Prandtl number. For wall-bounded problems, the algebraic models are found to best predict the rise in turbulent Prandtl number near the wall as well as the log-layer temperature profile, while the thermal variance models provide a good representation of the RMS temperature fluctuations. In jet flows, the algebraic models provide no benefit over a constant turbulent Prandtl number approach. Application of the thermal variance models finds that some significantly overpredict the temperature variance in the plume and most underpredict the thermal growth rate of the jet. The models yield very similar fluctuating temperature intensities in jets from straight pipes and smooth contraction nozzles, in contrast to data that indicate the latter should have noticeably higher values. For the particular low subsonic heated jet cases examined, changes in the turbulent Prandtl number had no effect on the centerline velocity decay.

computational fluid dynamics↗
NASA NTRS2016

Derivation and Validation of Supraglacial Lake Volumes on the Greenland Ice Sheet from High-Resolution Satellite Imagery

Supraglacial meltwater lakes on the western Greenland Ice Sheet (GrIS) are critical components of its surface hydrology and surface mass balance, and they also affect its ice dynamics. Estimates of lake volume, however, are limited by the availability of in situ measurements of water depth,which in turn also limits the assessment of remotely sensed lake depths. Given the logistical difficulty of collecting physical bathymetric measurements, methods relying upon in situ data are generally restricted to small areas and thus their application to largescale studies is difficult to validate. Here, we produce and validate spaceborne estimates of supraglacial lake volumes across a relatively large area (1250 km(exp 2) of west Greenland's ablation region using data acquired by the WorldView-2 (WV-2) sensor, making use of both its stereo-imaging capability and its meter-scale resolution. We employ spectrally-derived depth retrieval models, which are either based on absolute reflectance (single-channel model) or a ratio of spectral reflectances in two bands (dual-channel model). These models are calibrated by usingWV-2multispectral imagery acquired early in the melt season and depth measurements from a high resolutionWV-2 DEM over the same lake basins when devoid of water. The calibrated models are then validated with different lakes in the area, for which we determined depths. Lake depth estimates based on measurements recorded in WV-2's blue (450-510 nm), green (510-580 nm), and red (630-690 nm) bands and dual-channel modes (blue/green, blue/red, and green/red band combinations) had near-zero bias, an average root-mean-squared deviation of 0.4 m (relative to post-drainage DEMs), and an average volumetric error of b1%. The approach outlined in this study - image-based calibration of depth-retrieval models - significantly improves spaceborne supraglacial bathymetry retrievals, which are completely independent from in situ measurements.

Spectrally-based depth-retrieval models↗
NASA NTRS2004

Update on the Status of the FLUKA Monte Carlo Transport Code

The FLUKA Monte Carlo transport code is a well-known simulation tool in High Energy Physics. FLUKA is a dynamic tool in the sense that it is being continually updated and improved by the authors. Here we review the progresses achieved in the last year on the physics models. From the point of view of hadronic physics, most of the effort is still in the field of nucleus--nucleus interactions. The currently available version of FLUKA already includes the internal capability to simulate inelastic nuclear interactions beginning with lab kinetic energies of 100 MeV/A up the the highest accessible energies by means of the DPMJET-II.5 event generator to handle the interactions for greater than 5 GeV/A and rQMD for energies below that. The new developments concern, at high energy, the embedding of the DPMJET-III generator, which represent a major change with respect to the DPMJET-II structure. This will also allow to achieve a better consistency between the nucleus-nucleus section with the original FLUKA model for hadron-nucleus collisions. Work is also in progress to implement a third event generator model based on the Master Boltzmann Equation approach, in order to extend the energy capability from 100 MeV/A down to the threshold for these reactions. In addition to these extended physics capabilities, structural changes to the programs input and scoring capabilities are continually being upgraded. In particular we want to mention the upgrades in the geometry packages, now capable of reaching higher levels of abstraction. Work is also proceeding to provide direct import into ROOT of the FLUKA output files for analysis and to deploy a user-friendly GUI input interface.

Pinsky, L.↗
NASA NTRS1999

On Heating the Sun's Corona by Magnetic Explosions: Feasibility in Active Regions and prospects for Quiet Regions and Coronal Holes

We build a case for the persistent strong coronal heating in active regions and the pervasive quasi-steady heating of the corona in quiet regions and coronal holes being driven in basically the same way as the intense transient heating in solar flares: by explosions of sheared magnetic fields in the cores of initially closed bipoles. We begin by summarizing the observational case for exploding sheared core fields being the drivers of a wide variety of flare events, with and without coronal mass ejections. We conclude that the arrangement of an event's flare heating, whether there is a coronal mass ejection, and the time and place of the ejection relative to the flare heating are all largely determined by four elements of the form and action of the magnetic field: (1) the arrangement of the impacted, interacting bipoles participating in the event, (2) which of these bipoles are active (have sheared core fields that explode) and which are passive (are heated by injection from impacted active bipoles), (3) which core field explodes first, and (4) which core-field explosions are confined within the closed field of their bipoles and which ejectively open their bipoles. We then apply this magnetic-configuration framework for flare heating to the strong coronal heating observed by the Yohkoh Soft X-ray Telescope in an active region with strongly sheared core fields observed by the MSFC vector magnetograph. All of the strong coronal heating is in continually microflaring sheared core fields or in extended loops rooted against the active core fields. Thus, the strong heating occurs in field configurations consistent with the heating being driven by frequent core-field explosions that are smaller but similar to those in confined flares and flaring arches. From analysis of the thermal and magnetic energetics of two selected core-field microflares and a bright extended loop, we find that (1) it is energetically feasible for the sheared core fields to drive all of the coronal heating in the active region via a staccato of magnetic microexplosions, (2) the microflares at the feet of the extended loop behave as the flares at the feet of flaring arches in that more coronal heating is driven within the active bipole than in the extended loop, (3) the filling factor of the X-ray plasma in the core field microflares and in the extended loop is approximately 0.1, and (4) to release enough magnetic energy for a typical microflare (10^27 - 10^28 erg), a microflaring strand of sheared core field need expand and/or untwist by only a few percent at most. Finally, we point out that (1) the field configurations for strong coronal heating in our example active region (i.e., neutral-line core fields, many embedded in the feet of extended loops) are present in abundance in the magnetic network in quiet regions and coronal holes, and (2) it is known that many network bipoles do microflare and that many produce detectable coronal heating. We therefore propose that exploding sheared core fields are the drivers of most of the heating and dynamics of the solar atmosphere, ranging from the largest and most powerful coronal mass ejections and flares, to the vigorous microflaring and coronal heating in active regions, to the multitude of fine-scale explosive events in the magnetic network. The low-lysing exploding core fields in the network drive microflares, spicules, global coronal heating, and ,consequently, the solar wind.

Moore, Ronald L.↗
NASA NTRSDate not supplied

Systematic Errors in Simulated L-Band Brightness Temperature in the SMAP Level-4 Soil Moisture Analysis

The NASA Soil Moisture Active Passive (SMAP) mission has been providing L-band (1.4 GHz) brightness temperature (Tb) observations since April 2015. By assimilating the Tb observations into the NASA Catchment land surface model using a spatially distributed ensemble Kalman filter (EnKF), the SMAP Level-4 Soil Moisture (L4_SM) product provides global, 3-hourly, 9-km resolution estimates of surface (0-5 cm) and root-zone (0-100 cm) soil moisture with ~2.5-day latency for use in research and applications. The EnKF-based L4_SM analysis assumes unbiased forecast errors. Consequently, the seasonally varying bias between the model forecast Tb and the observed values is removed prior to the assimilation of the SMAP Tb observations. The L4_SM system is thus designed to only correct errors in synoptic-scale and interannual variations from the long-term mean seasonal cycle while maintaining the model’s (possibly wrong) climatology. In this paper, we examine the Tb observation-minus-forecast (O-F) residuals from the L4_SM Version 7 product (computed after rescaling the Tb observations to the mean seasonal cycle of the simulated Tb). The long-term average of the Tb O-F residuals has a global mean of only 0.13 K and locally small values, ranging from -1 to 3 K. The model forecast Tb, however, still exhibits undesirable systematic errors relative to the (rescaled) Tb observations. At some locations, the time-average Tb O-F values strongly depend on surface soil moisture (SM). At the Yanco SMAP core validation site, for example, the Tb O-F residuals typically range from 5 to 15 K under dry soil moisture conditions (SM < 0.15 m3 m-3) yet are predominantly negative under wet soil moisture conditions (SM > 0.25 m3 m-3), with values ranging from 0 to -40 K. This results in soil moisture analysis increments that persistently make the soil drier under dry SM conditions and persistently make the soil wetter under wet SM conditions, suggesting an error in the dynamic range of the simulated Tb, soil moisture or soil temperature. In this paper, we describe the higher-order systematic Tb forecast errors in more detail, examine their impact on the L4_SM product quality, and explore potential avenues to improve the L4_SM algorithm.

Rolf Reichle↗
NASA NTRS2005

MESSENGER: Exploring Mercury's Magnetosphere

The MESSENGER mission to Mercury offers our first opportunity to explore this planet s miniature magnetosphere since the brief flybys of Mariner 10. Mercury s magnetosphere is unique in many respects. The magnetosphere of Mercury is among the smallest in the solar system; its magnetic field typically stands off the solar wind only - 1000 to 2000 km above the surface. For this reason there are no closed drift paths for energetic particles and, hence, no radiation belts. The characteristic time scales for wave propagation and convective transport are short and kinetic and fluid modes may be coupled. Magnetic reconnection at the dayside magnetopause may erode the subsolar magnetosphere allowing solar wind ions to impact directly the regolith. Inductive currents in Mercury s interior may act to modify the solar wind interaction by resisting changes due to solar wind pressure variations. Indeed, observations of these induction effects may be an important source of information on the state of Mercury s interior. In addition, Mercury s magnetosphere is the only one with its defining magnetic flux tubes rooted in a planetary regolith as opposed to an atmosphere with a conductive ionospheric layer. This lack of an ionosphere is probably the underlying reason for the brevity of the very intense, but short-lived, - 1-2 min, substorm-like energetic particle events observed by Mariner 10 during its first traversal of Mercury s magnetic tail. Because of Mercury s proximity to the sun, 0.3 - 0.5 AU, this magnetosphere experiences the most extreme driving forces in the solar system. All of these factors are expected to produce complicated interactions involving the exchange and re-cycling of neutrals and ions between the solar wind, magnetosphere, and regolith. The electrodynamics of Mercury s magnetosphere are expected to be equally complex, with strong forcing by the solar wind, magnetic reconnection at the magnetopause and in the tail, and the pick-up of planetary ions all driving field-aligned electric currents. However, these field-aligned currents do not close in an ionosphere, but in some other manner. In addition to the insights- into magnetospheric physics offered by study of the solar wind - Mercury system, quantitative specification of the "external" magnetic field generated by magnetospheric currents is necessary for accurate determination of the strength and multi-polar decomposition of Mercury s intrinsic magnetic field. MESSENGER S highly capable instrumentation and broad orbital coverage will greatly advance our understanding of both the origin of Mercury s magnetic field and the acceleration of charged particles in small magnetospheres. In. this article, we review what is known about Mercury s magnetosphere and describe the MESSENGER science team s strategy for obtaining answers to the outstanding science questions surrounding the interaction of the solar wind with Mercury and its small, but dynamic, magnetosphere.

Slavin, James A.↗
NASA NTRS2022

Aboveground biomass density models for NASA’s Global Ecosystem Dynamics Investigation (GEDI) lidar mission

NASA’s Global Ecosystem Dynamics Investigation (GEDI) is collecting spaceborne full waveform lidar data with a primary science goal of producing accurate estimates of forest aboveground biomass density (AGBD). This paper presents the development of the models used to create GEDI’s footprint-level (~25 m) AGBD (GEDI04_A) product, including a description of the datasets used and the procedure for final model selection. The data used to fit our models are from a compilation of globally distributed spatially and temporally coincident field and airborne lidar datasets, whereby we simulated GEDI-like waveforms from airborne lidar to build a calibration database. We used this database to expand the geographic extent of past waveform lidar studies, and divided the globe into four broad strata by Plant Functional Type (PFT) and six geographic regions. GEDI’s waveform-to-biomass models take the form of parametric Ordinary Least Squares (OLS) models with simulated Relative Height (RH) metrics as predictor variables. From an exhaustive set of candidate models, we selected the best input predictor variables, and data transformations for each geographic stratum in the GEDI domain to produce a set of comprehensive predictive footprint-level models. We found that model selection frequently favored combinations of RH metrics at the 98th, 90th, 50th, and 10th height above ground-level percentiles (RH98, RH90, RH50, and RH10, respectively), but that inclusion of lower RH metrics (e.g. RH10) did not markedly improve model performance. Second, forced inclusion of RH98 in all models was important and did not degrade model performance, and the best performing models were parsimonious, typically having only 1-3 predictors. Third, stratification by geographic domain (PFT, geographic region) improved model performance in comparison to global models without stratification. Fourth, for the vast majority of strata, the best performing models were fit using square root transformation of field AGBD and/or height metrics. There was considerable variability in model performance across geographic strata, and areas with sparse training data and/or high AGBD values had the poorest performance. These models are used to produce global predictions of AGBD, but will be improved in the future as more and better training data become available.

Laura Duncanson↗
NASA NTRS2002

Observer-Based Magnetic Bearing Controller Developed for Aerospace Flywheels

A prototype of a versatile, observer-based magnetic bearing controller for aerospace flywheels was successfully developed and demonstrated on a magnetic bearing test rig (see the photograph) and an actual flywheel module. The objective of this development included a fast, yet low risk, control development process, and a robust, high-performance controller for a large variety of flywheels. This required a good system model, an efficient development procedure, and a model-based controller that addressed the key problems associated with flywheel and bearing imbalance, sensor error, and vibration. The model used in this control development and tuning procedure included the flexible rotor dynamics and motor-induced vibrations. Such a model was essential for low-risk scheduling of speed-dependent control parameters and for reliable evaluation of novel control strategies. The successfully tested control prototype utilized an extended Kalman filter to estimate the true rotor principal-axis motion from the raw sensor position feedback. For control refinement, the extended Kalman filter also estimated and eliminated the combined effects of mass-imbalance and sensor runouts from the input data. A key advantage of the design based on the extended Kalman filter is its ability to accurately estimate both the rotor's principal-axis position and gyroscopic rates with the least amount of phase lag. This is important for control parameter scheduling to dampen the gyroscopic motions. Because of large uncertainties in the magnetic bearing and imbalance characteristics, this state-estimation scheme alone is insufficient for containing the rotor motion within the desired 1-mil excursion radius. A nonlinear gain adjustment based on an estimation of the principal-axis orbit size was needed to provide a coarse (nonoptimal), but robust, control of the orbit growth. Control current minimization was achieved with a (steepest gradient) search of synchronous errors in the principal-axis position input data. Actual flywheel tests of this observer-based controller (developed entirely in-house) at the NASA Glenn Research Center showed that the model correctly predicted the rotor orbit growth as a function of rotational speed, and it demonstrated the capability of gain adjustments to arrest this growth. Data from these tests on an actual flywheel module spun to 26,000 rpm proved that the controller was able to contain the shaft motion to within much less than 0.5 mils of radial excursion with axis currents less than 300 mA in root-mean-square estimate. The test speed range was limited because of thermal expansion concerns for this particular flywheel unit, not because of any deficiency in the controller. Simulations for this unit indicated that the controller should be robust up to its top operating speed of 60,000 rpm. Aside from these important achievements, and most significantly, it took less than 1 week to adapt this controller from the simple test rig to the actual flywheel and to demonstrate full five-axis levitation and control. This demonstration showed that both the controller and the model-based development and tuning framework are easily adaptable to a wide range of rotors and bearing configurations and, hence, are capable of reducing design risks and costs for many future flywheel technology developments.

Le, Dzu K.↗
NASA NTRS1995

Stormtime ring current and radiation belt ion transport: Simulations and interpretations

We use a dynamical guiding-center model to investigate the stormtime transport of ring current and radiation-belt ions. We trace the motion of representative ions' guiding centers in response to model substorm-associated impulses in the convection electric field for a range of ion energies. Our simple magnetospheric model allows us to compare our numerical results quantitatively with analytical descriptions of particle transport, (e.g., with the quasilinear theory of radial diffusion). We find that 10-145-keV ions gain access to L approximately 3, where they can form the stormtime ring current, mainly from outside the (trapping) region in which particles execute closed drift paths. Conversely, the transport of higher-energy ions (approximately greater than 145 keV at L approximately 3) turns out to resemble radial diffusion. The quasilinear diffusion coefficient calculated for our model storm does not vary smoothly with particle energy, since our impulses occur at specific (although randomly determined) times. Despite the spectral irregularity, quasilinear theory provides a surprisingly accurate description of the transport process for approximately greater than 145-keV ions, even for the case of an individual storm. For 4 different realizations of our model storm, the geometric mean discrepancies between diffusion coefficients D(sup sim, sub LL) obtained from the simulations and the quasilinear diffusion coefficient D(sup ql, sub LL) amount to factors of 2.3, 2.3, 1.5, and 3.0, respectively. We have found that these discrepancies between D(sup sim, sub LL) and D(sup ql, sub LL) can be reduced slightly by invoking drift-resonance broadening to smooth out the sharp minima and maxima in D(sup ql, sub LL). The mean of the remaining discrepancies between D(sup sim, sub LL) and D(sup ql, sub LL) for the 4 different storms then amount to factors of 1.9, 2.1, 1.5, and 2.7, respectively. We find even better agreement when we reduce the impulse amplitudes systematically in a given model storm (e.g., reduction of all the impulse amplitudes by half reduces the discrepancy factor by at least its square root) and also when we average our results over an ensemble of 20 model storms (agreement is within a factor of 1.2 without impulse-amplitude reduction). We use our simulation results also to map phase-space densities f in accordance with Liouville's theorem. We find that the stormtime transport of approximately greater than 145-keV ions produces little change in f-bar the drift-averaged phase-space density on any drift shell of interest. However, the stormtime transport produces a major enhancement from the pre-storm phase-space density at energies approximately 30-145 keV, which are representative of the stormtime ring current.

Lyons, Larry R.↗
NASA NTRS1999

Dynamics of Single Chains of Suspended Ferrofluid Particles

We present an experimental study of the dynamics of isolated chains made of super-paramagnetic particles under the influence of a magnetic field. The motivation of this work is to understand if the chain fluctuations exist and, if it does, how does the fluctuation affect chain aggregation. We find that single chains strongly fluctuate and that the characteristic frequency of their fluctuations is inversely proportional to the magnetic field strength. The higher the field the lower the characteristic frequency of the chain fluctuations. In the high magnetic field limit, chains behave like rigid rods without any internal motions. In this work, we used ferrofluid particles suspended in water. These particles do not have any intrinsic magnetization. Once a magnetic field is applied, a dipole moment is induced in each particle, proportional to the magnetic field. A dipolar magnetic interaction then occurs between particles. If dipole-dipole magnetic energy is higher than the thermal energy, the result is a structure change inside the dipolar fluid. The ratio of these two energies is expressed by a coupling constant lambda as: lambda = (pi(a(exp 3))(chi(exp 2))(mu(sub 0))(H(sub 0))(exp 2))/18kT Where a is the particle radius, mu(sub 0) is the vacuum magnetic permeability, H(sub 0) the applied magnetic field, k the Boltzmann constant and T the absolute temperature. If lambda > 1, magnetic particles form chains along the field direction. The lateral coalescence of several chains may form bigger aggregates especially if the particle volume fraction is high. While many studies and applications deal with the rheological properties and the structural changes of these dipolar fluids, this work focuses on the understanding of the chain dynamics. In order to probe the chain dynamics, we used dynamic light scattering (DLS) in self-beating mode as our experimental technique. The experimental geometry is such that the scattering plane is perpendicular to the magnetic field. Therefore, only motions in this plane are probed. A very dilute sample of a ferrofluid emulsion with a particle volume fraction of 10(exp -5) is used in this experiment. We chose such a low volume fraction to avoid multiple light scattering as well as lateral chain-chain aggregation. DLS measures the dynamic structure factor S(q,t) of the sample (q is the scattering wave vector, t is the time). In the absence of the magnetic field, identical particles of ferrofluid droplets are randomly distributed and S(q,t) reduces to exp(-q(exp 2)2D(sub 0)t). D(sub 0)=(kT/(6(pi)(eta)(a)) is the diffusion coefficient of Brownian particles (where Xi = (6(pi)(eta)(a)) is the Stokes frictional coefficient of a spherical particle in a fluid of viscosity eta). If interactions or polydispersity can not be ignored, an effective diffusion coefficient is introduced. Formally, D(sub eff) is defined as: D(sub eff) = - q(exp -2) partial derivative of (ln(S(q,t)) with respect to time, as t goes to 0. D(sub eff) reduces to D(sub 0) if no interactions and only a few particles size are present. Therefore, we can use DLS to measure particle size. The particle radius was found to be a=0.23 mu m with 7% of polydispersity. In this case, if we vary the scattering angle theta (and so q) we do not have any change in the measured diffusion coefficient: it is q-independent. When a magnetic field is applied, particles aggregate into chains if lambda > 1. We first studied the kinetics of the chain formation when lambda = 406. At a fixed scattering angle, we measured diffusion coefficient D(sub eff) as a function of time. Experimentally, we find that D(sub eff) decreases monotonously with time. Physically, this means that chains are becoming longer and longer. Since we are only sensitive to motions in the scattering plane and since chains have their main axis perpendicular to this plane, the measured diffusion coefficient is the trans-verse diffusion coefficient. We can relate D(sub eff) to the mean number of particles per chain N(t) at a given time and to the diffusion coefficient of an isolated particle D(sub 0) as D(sub eff)=f(N(t))D(sub 0). Since f(N) is known from other recent work, N can be expressed as a function of the time. We found a square root dependency: N(t) proportional to the square root of t. As expected for very low volume fraction, this behavior is characteristic of a diffusion-limited aggregation as suggested by several authors and by our previous work. In this study, we focus on the dependence of the effective diffusion coefficient on the scattering angle and the magnetic field strength. After the magnetic field is applied (lambda = 406) for a long time, typically 6 hours, kinetics of chain formation becomes very slow. Chain size does not vary much over the next hour period. Thus, we can perform different interesting experiments. First, at a fixed magnetic field, we measure the effective diffusion coefficient as a function of the scattering angle (from 5 to 130 deg). Our results show that the measured diffusion coefficient increases linearly with the scattering angle: D(sub eff) proportional to q. If we do the same experiment for different lambda values, D(sub eff) depends on lambda as D(sub eff) proportional to lambda(exp -1/2). We also find for different lambda values that the same asymptotic D(sub eff) value is obtained when q approaches zero. The angle dependency of D(sub eff) suggests that an additional motion exists besides chain drifting. Chain size is constant during experiment, which was verified by measuring the same diffusion coefficient at the beginning and at the end of the angle switching. If chains are rigid, D(sub eff) is independent of q. Therefore, we found that D(sub eff) not only measures the motion of the entire chain but also its internal fluctuations. These internal motions are the fluctuations of the particles in the chain. To understand the q dependency of D(sub eff), let us look at the probing length used. In our study, the characteristic length scale probed is l=2pi/q which is in the range of 0.9<l/a<20. When l is much larger than the particle radius (l/a going to infinity) we are mainly sensitive to the center of mass diffusion of the chain. D(sub eff) is then the diffusion of the entire chain and depends only on N but not on q and lambda. The value D(sub eff) allows us to obtain the number of particles per chain N. In the opposite limit (i.e. l/a < 1), we are sensitive to motions on the size of individual particles. The main contributions to the measured diffusion coefficient come from internal motions of the chains (i.e. particles' fluctuations inside the chain). We investigated also the effect of the chain size on D(sub eff). For the same value of the magnetic interaction lambda and different chain sizes, we found that the slope of D(sub eff) versus q decreases when N increases as a power law D(sub eff) proportional to N(exp -0.7). This is an important behavior because this means that the longer the chain the slower its internal fluctuations besides the whole chain motion.

Cutillas, S.↗
NASA NTRS1998

Diurnal Reflectance Changes in Vegetation Observed with AVIRIS

Among the most important short-term dynamic biological processes are diurnal changes in canopy water relations. Plant regulation of water transport through stomatal openings affects other gaseous transport processes, often dramatically decreasing photosynthetic fixation of carbon dioxide during periods of water stress. Water stress reduces stomatal conductance of water vapor through the leaf surface and alters the diurnal timing of stomatal opening. Under non-water stressed conditions, stomates typically open soon after dawn and transpire water vapor throughout the daylight period. During stress periods, stomates may close for part of the day, generally near mid-day. Under prolonged stress conditions, stomatal closure shifts to earlier times during the day; stomates may close by mid-morning and remain closed until the following morning - or remain closed entirely. Under these conditions the relationship between canopy greenness (e.g., measured with a vegetation index or by spectral mixture analysis) and photosynthetic fixation of carbon is lost and the remotely sensed vegetation metric is a poor predictor of gas exchange. Prediction of stomatal regulation and exchange of water and trace gases is critical for ecosystem and climate models to correctly estimate budgets of these gases and understand or predict other processes like gross and net ecosystem primary production. Plant gas exchange has been extensively studied by physiologists at the leaf and whole plant level and by biometeorologists at somewhat larger scales. While these energy driven processes follow a predictable if somewhat asymmetric diurnal cycle dependent on soil water availability and the constraints imposed by the solar energy budget, they are nonetheless difficult to measure at the tree and stand levels using conventional methods. Ecologists have long been interested in the potential of remote sensing for monitoring physiological changes using multi-temporal images. Much of this research has focused on day-to-day changes in water use, especially for agricultural applications. Ustin et al. showed seasonal changes in canopy water content in chaparral shrub could be estimated using optical methods. Vanderbilt et al. followed asymmetric diurnal changes in the reflectance of a walnut orchard, but could not attribute specific reflectance changes to specific changes in canopy architecture or physiology. Forests and shrub lands in California experience prolonged periods of drought, sometimes extending six months without precipitation. The conifer and evergreen chaparral communities common to the foothill region around the central valley of California retain their foliage throughout the summer and have low transpiration rates despite high net radiation and temperature conditions. In contrast, grasslands and drought resistant deciduous species in the same habitat are seasonally dormant in summer. Because of differences in the mechanisms of drought tolerance, rooting depth and physiology between different plant communities in the region, it is likely that they display differences in diurnal water relations. The presence of diverse plant communities provides an opportunity to investigate possible diurnal landscape patterns in water relations that could be observed by an airborne hyperspectral scanner. This investigation of AVIRIS data collected over forest and shrub land represents the continuation of a prior investigation involving spectral mixture analysis of diurnal effects in the same AVIRIS data set.

Vanderbilt, V. C.↗
NASA NTRSDate not supplied

Examination of L-Band Brightness Temperature Forecasts in the SMAP Level-4 Soil Moisture Analysis

The NASA Soil Moisture Active Passive (SMAP) mission [4] has been providing L-band (1.4 GHz) passive microwave brightness temperature (Tb) observations since April 2015. By assimilating the Tb observations into the NASA Catchment land surface model [5] using a spatially distributed ensemble Kalman filter, the NASA Global Modeling and Assimilation Office generates the SMAP Level-4 Soil Moisture (L4_SM) product, which provides global, 3-hourly, 9-km resolution estimates of surface (0-5 cm) and root-zone (0-100 cm) soil moisture with ~2.5-day latency for use in research and applications [6]. The L4_SM product also includes estimates of soil temperature, land surface fluxes, and assimilation diagnostics such as the model forecast and observed Tb values [7]. The output from the L4_SM system is routinely monitored by the L4_SM team. Such monitoring provides valuable information; instances of unusually large Tb observation-minus-forecast (O-F) residuals can indicate events during which soil moisture conditions are poorly described in the land modeling system [7,8]. For example, repeated occurrences of very large Tb O-F values in central Australia were traced back to deficiencies in the gauge-based precipitation product used in the land modeling system through L4_SM Version 5; this discovery prompted the use, beginning in Version 6, of satellite- and gauge-based precipitation observations outside of North America [9]. Recently, a statistical analysis of the Tb O-F residuals revealed systematic errors in the “tau-omega” L-band radiative transfer model that converts the land-model simulated soil moisture and temperature into the forecast Tb prior to the L4_SM analysis. The L4_SM analysis is built on the ensemble Kalman filter and assumes unbiased forecast errors. The presence of systematic Tb forecast errors could thus adversely impact the quality of the analyzed soil moisture. In this paper, we examine the Tb O-F residuals of the latest Version 7 L4_SM data (Science Version ID Vv7030 and Vv7032) [10,11,12]. In the L4_SM Version 7 algorithm, key parameters of the L-band radiative transfer model, including soil roughness, scattering albedo, and a (seasonally varying) climatology of vegetation opacity, are obtained from the SMAP Level-2 Radiometer retrieval product (Version 5) [1]. Additionally, we also use ground measurements of surface soil moisture and soil temperature from the SMAP core validation sites as in situ reference of soil conditions [2,3]. As part of the L4_SM system calibration, the seasonally-varying bias between the model forecast Tb and the observed values is removed prior to the assimilation of the SMAP Tb observations [7]. That is, the L4_SM system is designed to only correct errors in synoptic-scale and interannual variations from the long-term mean seasonal cycle while maintaining the model’s (potentially erroneous) climatology. The Tb O-F residuals examined here are thus computed using SMAP Tb observations after they are rescaled to the mean seasonal cycle of Tb from the modeling system without SMAP data assimilation. Consequently, the long-term average of the Tb O-F residuals has a global mean of only 0.13 K and locally small values, ranging from -1 to 3 K. Despite the small time-average values of the Tb O-F residuals, the model forecast Tb was nevertheless found to exhibit undesirable systematic errors. At some locations, the time-average Tb O-F values strongly depend on surface soil moisture (SM). At the Yanco SMAP core validation site, for example, the Tb O-F residuals typically range from 5 to 15 K under dry soil moisture conditions (SM < 0.15 m3 m-3) yet are predominantly negative under wet soil moisture conditions (SM > 0.25 m3 m-3), with values ranging from 0 to -40 K. This results in soil moisture analysis increments that persistently make the soil drier under dry SM conditions and persistently make the soil wetter under wet SM conditions, suggesting an error in the dynamic range of the simulated Tb, soil moisture or soil temperature. In this paper, we describe the higher-order systematic Tb forecast errors in more detail, examine their impact on the L4_SM product quality, and explore potential avenues to further improve the L4_SM algorithm.

Rolf Reichle↗