Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “global shape 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 433 records · Page 24
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

FFT-based solver for higher-order and multi-phase-field fracture models applied to strongly anisotropic brittle materials

This paper presents the application of a fast Fourier transform (FFT) based method to solve two phase field models designed to simulate crack growth of strongly anisotropic materials in the brittle regime. By leveraging the ability of the FFT-based solver to generate solutions with higher-order and global continuities, we design two simple algorithms to capture the complex fracture patterns (e.g. sawtooth, and curved crack growth) common in materials with strongly anisotropic surface energy via the multi-phase-field and high-order phase-field frameworks. A staggered operator-split solver is used where both the balance of linear momentum and the phase field governing equations are formulated in the periodic domain. The unit phase field of the initial failure region is prescribed by the penalty method to alleviate the sharp material contrast between the initial failure region and the base material. The discrete frequency vectors are generalized to estimate the second and fourth order gradients such that the Gibbs effect near shape interfaces or jump conditions can be suppressed. Furthermore, a preconditioner is adopted to improve the convergence rate of the iterative linear solver. Three numerical experiments are used to systematically compare the performance of the FFT-based method in the multi-phase-field and high-order phase-field frameworks.

42 ENGINEERING↗
NASA NTRS2016

Weather and Climate Change Impacts on Human Mortality in Bangladesh

Weather and climate profoundly affect human health. Several studies have demonstrated a U-, V-, or J-shaped temperature-mortality relationship with increasing death rates at the lower and particularly upper end of the temperature distribution. The objectives of this study were (1) to analyze the relationship between temperature and mortality in Bangladesh for different subpopulations and (2) to project future heat-related mortality under climate change scenarios. We used (non-)parametric Generalized Additive Models adjusted for trend, season and day of the month to analyze the effect of temperature on daily mortality. We found a decrease in mortality with increasing temperature over a wide range of values; between the 90th and 95th percentile an abrupt increase in mortality was observed which was particularly pronounced for the elderly above the age of 65 years, for males, as well as in urban areas and in areas with a high socio-economic status. Daily historical and future temperature values were obtained from the NASA Earth Exchange Global Daily Downscaled Projections (NEX-GDDP) dataset. This dataset is comprised of downscaled climate scenarios for the globe that are derived from the General Circulation Model (GCM) runs conducted under the Coupled Model Intercomparison Project Phase 5 (CMIP5). The derived dose-response functions were used to estimate the number of heat-related deaths occurring during the 1990s (1980-2005), the 2020s (2011-2040) and the 2050s (2041-2070). We estimated that excess deaths due to heat will triple from the 1990s to the 2050s, with an annual number of 0.5 million excess deaths in 1990 to and expected number of 1.5 millions in 2050.

mortality↗
NASA NTRS1996

Towards and FVE-FAC Method for Determining Thermocapillary Effects on Weld Pool Shape

Several practical materials processes, e.g., welding, float-zone purification, and Czochralski crystal growth, involve a pool of molten metal with a free surface, with strong temperature gradients along the surface. In some cases, the resulting thermocapillary flow is vigorous enough to convect heat toward the edges of the pool, increasing the driving force in a sort of positive feedback. In this work we examine this mechanism and its effect on the solid-liquid interface through a model problem: a half space of pure substance with concentrated axisymmetric surface heating, where surface tension is strong enough to keep the liquid free surface flat. The numerical method proposed for this problem utilizes a finite volume element (FVE) discretization in cylindrical coordinates. Because of the axisymmetric nature of the model problem, the control volumes used are torroidal prisms, formed by taking a polygonal cross-section in the (r, z) plane and sweeping it completely around the z-axis. Conservation of energy (in the solid), and conservation of energy, momentum, and mass (in the liquid) are enforced globally by integrating these quantities and enforcing conservation over each control volume. Judicious application of the Divergence Theorem and Stokes' Theorem, combined with a Crank-Nicolson time-stepping scheme leads to an implicit algebraic system to be solved at each time step. It is known that near the boundary of the pool, that is, near the solid-liquid interface, the full conduction-convection solution will require extremely fine length scales to resolve the physical behavior of the system. Furthermore, this boundary moves as a function of time. Accordingly, we develop the foundation of an adaptive refinement scheme based on the principles of Fast Adaptive Composite Grid methods (FAC). Implementation of the method and numerical results will appear in a later report.

Canright, David↗
NASA NTRS1995

A new model for Mars atmospheric dust based upon analysis of ultraviolet through infrared observations from Mariner 9, Viking, and Phobos

We propose key modifications to the Toon et al. (1977) model of the particle size distribution and composition of Mars atmospheric dust, based on a variety of spacecraft and wavelength observations of the dust. A much broader (r(sub eff) variance approximately 0.8 micrometers), smaller particle size (r(sub mode) approximately 0.02 micrometers) distribution coupled with a 'palagonite-like' composition is argued to fit the complete ultraviolet-to-30-micrometer absorption properties of the dust better than the montmorillonite-basalt, r(sub eff) variance = 0.4 micrometers, r(sub mode) = 0.40 dust model of Toon et al. Mariner 9 (infrared interferometer spectrometer) IRIS spectra of high atmospheric dust opacities during the 1971-1972 Mars global dust storm are analyzed in terms of the Toon et al. dust model, and a Hawaiian palagonite sample (Rousch et al., 1991) with two different size distribution models incorporating smaller dust particle sizes. Viking Infrared Thermal Mapper (IRTM) emmission-phase-function (EPF) observations at 9 micrometers are analyzed to retrieve 9-micrometer dust opacities coincident with solar band dust opacities obtained from the same EPF sequences (Clancy and Lee, 1991). These EPF dust opacities provide an independent measurement of the visible/9-micrometer extinction opacity ratio (greater than or = 2) for Mars atmospheric dust, which is consistent with a previous measurement by Martin (1986). Model values for the visible/9-micrometer opacity ratio and the ultraviolet and visible single-scattering albedos are calculated for the palagonite model with the smaller particle size distributions compared to the same properties for the Toon et al. model of dust. The montmorillonite model of the dust is found to fit the detailed shape of the dust 9-micrometer absorption well. However, it predicts structured, deep aborptions at 20 micrometers which are not observed and requires a separate ultraviolet-visible absorbing component to match the observed behavior of the dust in this wavelength region. The modeled palagonite does not match the 8-to 9-micrometer absorption presented by the dust in the IRIS spectra, probably due to its low SiO2 content (31%). However, it does provide consistent levels of ultraviolet/visible absorption, 9-to 12-micrometer absorption, and a lack of structured absorption at 20 micrometers. The ratios of dust extinction opacities at visible, 9 micrometers, and 30 micrometers are strongly affected by the dust particle size distribution. The Toon et al. dust size distribution (r(sub mode) = 0.40,r(sub eff) variance = 0.4 micrometers, r(sub cwmu) = 2.7 micrometers) predicts the correct ratio of the 9- to 30-micrometer opacity, but underpredicts the visible/9-micrometer opacity ratio considerably (1 versus greater than or = 2). A similar particle distribution width with smaller particle sizes (r(sub mode) = 0.17, r(sub eff) variance = 0.4 micrometers, r(sub cwmu) = 1.2 micrometers) will fit the observed visible/9-micrometer opacity ratio, but overpredicts the observed 9-micrometer/30-micrometer opacity ratio. A smaller and much broader particle size distribution (r(sub mode) = 0.002, r(sub eff) variance = 0.8 micrometers, r(sub cwmu) = 1.8 micrometers) can fit both dust opacity ratios. Overall, the nanocrystalline structure of palagonite coupled with a smaller, broader distribution of dust particle sizes provides a more consistent fit than the Toon et al. model of the dust to the IRIS spectra, the observed visible/9-micrometer dust opacity ratio, the Phobos occulation measurements of the dust particle sizes (Chassefiere et al., 1992), and the weakness of surface near IR absorptions expected for clay minerals (Clark, 1992; Bell and Crisp, 1993).

Clancy, R. T.↗
NASA NTRS1996

Development of Multiobjective Optimization Techniques for Sonic Boom Minimization

A discrete, semi-analytical sensitivity analysis procedure has been developed for calculating aerodynamic design sensitivities. The sensitivities of the flow variables and the grid coordinates are numerically calculated using direct differentiation of the respective discretized governing equations. The sensitivity analysis techniques are adapted within a parabolized Navier Stokes equations solver. Aerodynamic design sensitivities for high speed wing-body configurations are calculated using the semi-analytical sensitivity analysis procedures. Representative results obtained compare well with those obtained using the finite difference approach and establish the computational efficiency and accuracy of the semi-analytical procedures. Multidisciplinary design optimization procedures have been developed for aerospace applications namely, gas turbine blades and high speed wing-body configurations. In complex applications, the coupled optimization problems are decomposed into sublevels using multilevel decomposition techniques. In cases with multiple objective functions, formal multiobjective formulation such as the Kreisselmeier-Steinhauser function approach and the modified global criteria approach have been used. Nonlinear programming techniques for continuous design variables and a hybrid optimization technique, based on a simulated annealing algorithm, for discrete design variables have been used for solving the optimization problems. The optimization procedure for gas turbine blades improves the aerodynamic and heat transfer characteristics of the blades. The two-dimensional, blade-to-blade aerodynamic analysis is performed using a panel code. The blade heat transfer analysis is performed using an in-house developed finite element procedure. The optimization procedure yields blade shapes with significantly improved velocity and temperature distributions. The multidisciplinary design optimization procedures for high speed wing-body configurations simultaneously improve the aerodynamic, the sonic boom and the structural characteristics of the aircraft. The flow solution is obtained using a comprehensive parabolized Navier Stokes solver. Sonic boom analysis is performed using an extrapolation procedure. The aircraft wing load carrying member is modeled as either an isotropic or a composite box beam. The isotropic box beam is analyzed using thin wall theory. The composite box beam is analyzed using a finite element procedure. The developed optimization procedures yield significant improvements in all the performance criteria and provide interesting design trade-offs. The semi-analytical sensitivity analysis techniques offer significant computational savings and allow the use of comprehensive analysis procedures within design optimization studies.

Chattopadhyay, Aditi↗
DOE OSTI2020

Evaluation of simulated soil carbon dynamics in Arctic-Boreal ecosystems

Given the magnitude of soil carbon stocks in northern ecosystems, and the vulnerability of these stocks to climate warming, land surface models must accurately represent soil carbon dynamics in these regions. We evaluate soil carbon stocks and turnover rates, and the relationship between soil carbon loss with soil temperature and moisture, from an ensemble of eleven global land surface models. We focus on the region of NASA's Arctic-Boreal vulnerability experiment (ABoVE) in North America to inform data collection and model development efforts. Models exhibit an order of magnitude difference in estimates of current total soil carbon stocks, generally under- or overestimating the size of current soil carbon stocks by greater than 50 PgC. We find that a model's soil carbon stock at steady-state in 1901 is the prime driver of its soil carbon stock a hundred years later—overwhelming the effect of environmental forcing factors like climate. The greatest divergence between modeled and observed soil carbon stocks is in regions dominated by peat and permafrost soils, suggesting that models are failing to capture the frozen soil carbon dynamics of permafrost regions. Using a set of functional benchmarks to test the simulated relationship of soil respiration to both soil temperature and moisture, we find that although models capture the observed shape of the soil moisture response of respiration, almost half of the models examined show temperature sensitivities, or Q10 values, that are half of observed. Significantly, models that perform better against observational constraints of respiration or carbon stock size do not necessarily perform well in terms of their functional response to key climatic factors like changing temperature. This suggests that models may be arriving at the right result, but for the wrong reason. The results of this work can help to bridge the gap between data and models by both pointing to the need to constrain initial carbon pool sizes, as well as highlighting the importance of incorporating functional benchmarks into ongoing, mechanistic modeling activities such as those included in ABoVE.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

High fusion performance in Super H-mode experiments on Alcator C-Mod and DIII-D

The ‘Super H-Mode’ regime is predicted to enable pedestal height and fusion performance substantially higher than standard H-Mode operation. This regime exists due to a bifurcation of the pedestal pressure, as a function of density, that is predicted by the EPED model to occur in strongly shaped plasmas above a critical pedestal density. Experiments on Alcator C-Mod and DIII-D have achieved access to the Super H-Mode (and Near Super H) regime, and obtained very high pedestal pressure, including the highest achieved on a tokamak (pped ~ 80 kPa) in C-Mod experiments operating near the ITER magnetic field. DIII-D Super H experiments have demonstrated strong performance, including the highest stored energy in the present configuration of DIII-D (W ~ 2.2–3.2 MJ), while utilizing only about half of the available heating power (Pheat ~ 7–12 MW). These DIII-D experiments have obtained the highest value of peak fusion gain, QDT,equiv ~ 0.5, achieved on a medium scale (R < 2 m) tokamak. Sustained high performance operation (βN ~ 2.9, H98 ~ 1.6) has been achieved utilizing n = 3 magnetic perturbations for density and impurity control. Pedestal and global confinement has been maintained in the presence of deuterium and nitrogen gas puffing, which enables a more radiative divertor condition. A pair of simple performance metrics is developed to assess and compare regimes. Super H-Mode access is predicted for ITER and expected, based on both theoretical prediction and observed normalized performance, to allow ITER to achieve its goals (Q = 10) at Ip < 15 MA, and to potentially enable more compact, cost effective pilot plant and reactor designs.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1993

Porosity and the ecology of icy satellites

The case for a significant role for porosity in the structure and evolution of icy bodies in the Solar System has been difficult to establish. We present a relevant new data set and a series of structure models including a mechanical compression, not thermal creep, model for porosity that accounts satisfactorily for observed densities, moments of inertia, geologic activity, and sizes of tectonic features on icy satellites. Several types of observational data sets have been used to infer significant porosity, but until recently, alternative explanations have been preferred. Our first area of concern is the occurrence of cryovolcanism as a function of satellite radius; simple radiogenic heating models of icy satellites suggest minimum radii for melting and surface cryovolcanism to be 400 to 500 km, yet inferred melt deposits are seen on satellites half that size. One possible explanation is a deep, low conductivity regolith which lowers conductivity and raises internal temperatures, but other possibilities include tidal heating or crustal compositions of low conductivity. Our second area of concern is the occurrence and magnitude of tectonic strain; tectonic structures have been seen on icy satellites as small as Mimas and Proteus. The structures are almost exclusively extensional, with only a few possible compression Al features, and inferred global strains are on the order of 1 percent expansion. Expansions of this order in small bodies like Mimas and prevention of late compressional tectonics due to formation of ice mantles in larger bodies like Rhea are attained only in structure models including low-conductivity, and thus possibly high porosity, crusts. Thirdly, inferred moments of inertia less than 0.4 in Mimas and Tethys can be explained by high-porosity crusts, but also by differentiation of a high density core. Finally, the relatively low densities of smaller satellites like Mimas and Miranda relative to larger neighbors can be explained by deep porosity, but also by bulk compositional differences. Recent work has strengthened the case for significant porosity. Halley's nucleus was found to have a density near 0.6 g/cu cm, Janus and Epimethus were proposed to have densities near 0.7 g/cu cm, densities almost certainly due to high porosity. The irregular-spherical shape transition of icy satellites was quantitatively explained by low conductivity regoliths. A creative structure/thermal history model for Mimas simultaneously accounts quantitatively for Mimas' low density and moment of inertia by invoking initial high-porosity and subsequent compaction in the deep interior by thermal creep. The main problem with this promising model is that approximately 7 percent predicts a reduction in Mimas' radius, implying significant compressional failure and prevention of extensional tectonics, in contradiction to the observed extensional features and inferred 1 percent expansion in radius.

Croft, Steven K.↗
NASA NTRS2002

Localized Gravity/Topography Admittance and Correlation Spectra on Mars: Implications for Regional and Global Evolution

[i] From gravity and topography data collected by the Mars Global Surveyor spacecraft we calculate gravity/topography admittances and correlations in the spectral domain and compare them to those predicted from models of lithospheric flexure. On the basis of these comparisons we estimate the thickness of the Martian elastic lithosphere (T(sub e)) required to support the observed topographic load since the time of loading. We convert T(sub e) to estimates of heat flux and thermal gradient in the lithosphere through a consideration of the response of an elastic/plastic shell. In regions of high topography on Mars (e.g., the Tharsis rise and associated shield volcanoes), the mass-sheet (small-amplitude) approximation for the calculation of gravity from topography is inadequate. A correction that accounts for finite-amplitude topography tends to increase the amplitude of the predicted gravity signal at spacecraft altitudes. Proper implementation of this correction requires the use of radii from the center of mass (collectively known as the planetary shape ) in lieu of topography referenced to a gravitational equipotential. Anomalously dense surface layers or buried excess masses are not required to explain the observed admittances for the Tharsis Montes or Olympus Mons volcanoes when this correction is applied. Derived T, values generally decrease with increasing age of the lithospheric load, in a manner consistent with a rapid decline of mantle heat flux during the Noachian and more modest rates of decline during subsequent epochs.

McGovern, Patrick J.↗
DOE OSTI2025

Improving Projections of Antarctic Ice Sheet Contribution to Sea‐Level Change Through 2300 by Capturing Gravitational, Rotational, and Deformational Effects

Ice-mass change induces regionally varying patterns of sea-level change due to gravitational, rotational, and deformational (GRD) effects, which in turn influence marine-based ice stability in Antarctica. For improved projection of the Antarctic Ice Sheet (AIS), there is a need for including GRD effects in modeling and improving understanding of basin-by-basin sensitivity of ice evolution to GRD effects under a range of climate scenarios. We couple a high-resolution, higher-order ice-sheet model with a 1D global sea-level model that fully captures GRD effects, and simulate ice evolution in Antarctica under the Ice Sheet Model Intercomparison Project for CMIP6 experiments. We perform two sets of coupled simulations incorporating 1D Maxwell solid Earth structure suitable for West and East Antarctica and show that the Amundsen Sea Embayment (ASE) in West Antarctica has the highest sensitivity to GRD effects—in high-emission scenarios, grounding-line retreat accelerates by hundreds of kilometers by 2300 without GRD effects, but GRD effects delay this retreat on a timescale of decades. However, we find that delay times do not show a clear relationship to the strength of climate forcing alone. Furthermore, GRD effects can influence ice-sheet dynamics more than the choice of climate model for a given emissions scenario. In contrast, East Antarctica exhibits minimal sensitivity to GRD effects throughout the study period. These findings underscore the critical role of GRD effects in shaping future West AIS evolution, highlighting the importance of constraining the regional 3D Earth structure and bed topography in West Antarctica, particularly the ASE.

58 GEOSCIENCES↗
NASA NTRS2022

NASA’s Global Precipitation Measurement Mission: Leveraging Stakeholder Engagement & Applications Activities to Inform Decision-making

The application of satellite precipitation estimates from NASA’s Global Precipitation Measurement (GPM) Mission for decision-making has been a focus for the mission since launch. As a result, GPM data have enabled a range of applications that address societal needs, including water resource management, crop forecasting, ecological monitoring, disaster response, public health, aviation, weather forecasting, and climate modeling, among others. GPM applications activities have continued to focus on user engagement through in person trainings and interviews, workshops, webinars, and educational outreach activities. The goals of these efforts are to synthesize community data needs in order to effectively support and enable decision-making across agencies, academia and the global community. While these efforts have helped the GPM mission establish a large stakeholder community that encompasses federal and state partners, academic institutions, nd private and nonprofit companies, there remains difficulties associated with accessing, processing, and applying the data to support or enable applications. In this article, we present GPM applications strategies and approaches used to enhance the applications value of GPM data, and most importantly, demonstrate how these efforts have and can inform different decision-making contexts. This work also provides a discussion on key lessons learned from the user community and how this information can be utilized to help better support and shape applications approaches for future NASA Earth Science missions.

Satellite precipitation↗
DOE OSTI2022

Linear simulation of magnetohydrodynamic plasma response to three-dimensional magnetic perturbations in high-β P plasmas

In this work, we report the numerical analyses of linear magnetohydrodynamics (MHD) plasma response to applied three-dimensional magnetic perturbations (MPs) in a joint DIII-D/EAST collaboration on high-β P (poloidal beta) plasmas, utilizing the extended-MHD code M3D-C1, with the purpose of realizing a better understanding of the existing experiment in which the n=3 MPs were applied to such high-β P plasmas attempting to control large amplitude type-I ELMs. Such high-β P plasmas obtained at the DIII-D tokamak feature an upper-biased double null configuration, a high edge safety factor q 95 ~6.4, and a stable internal transport barrier (ITB) leading to relatively high core pressures. Single-fluid simulations show that the plasma response to n=3 MPs, including both non-resonant/kinking and resonant components, is significantly weaker than that to n=1 or 2 MPs. To survey the impact of q 95 on plasma response to applied MPs, the SEGWAY (Self-consistent Equilibrium Generating Workflow for AnalYsis) module, developed in the OMFIT integrated modelling framework, is employed to generate a series of equilibria with a wide range of q 95 while other key parameters including the normalized beta, electron density at pedestal top, and plasma shape are kept fixed. Compared to the vacuum response, single-fluid M3D-C1 simulations predict a much more significant decrease of resonant plasma response to the applied n=3 MPs at the maximum penetration radii as q 95 increases. In contrast to single-fluid simulation results showing resonant penetration occurs only near the pedestal top where the E×B toroidal rotation frequency is zero, two-fluid simulations show two comparable resonant penetrations locating near the pedestal top and the ITB foot, where the perpendicular electron rotation frequency is zero. Such resonant field penetration near the ITB foot may be responsible for the observed formation of a staircase structure in both electron density and temperature profiles and thereby a considerable deterioration of global plasma performance when MPs are applied in high-β P plasmas. Motivated by this numerical work, we provide some ideas for the future research, with the purpose of realizing effective ELM control in such high-β P plasmas on the DIII-D and EAST devices.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2014

NASA Tech Briefs, March 2014

Topics include: Data Fusion for Global Estimation of Forest Characteristics From Sparse Lidar Data; Debris and Ice Mapping Analysis Tool - Database; Data Acquisition and Processing Software - DAPS; Metal-Assisted Fabrication of Biodegradable Porous Silicon Nanostructures; Post-Growth, In Situ Adhesion of Carbon Nanotubes to a Substrate for Robust CNT Cathodes; Integrated PEMFC Flow Field Design for Gravity-Independent Passive Water Removal; Thermal Mechanical Preparation of Glass Spheres; Mechanistic-Based Multiaxial-Stochastic-Strength Model for Transversely-Isotropic Brittle Materials; Methods for Mitigating Space Radiation Effects, Fault Detection and Correction, and Processing Sensor Data; Compact Ka-Band Antenna Feed with Double Circularly Polarized Capability; Dual-Leadframe Transient Liquid Phase Bonded Power Semiconductor Module Assembly and Bonding Process; Quad First Stage Processor: A Four-Channel Digitizer and Digital Beam-Forming Processor; Protective Sleeve for a Pyrotechnic Reefing Line Cutter; Metabolic Heat Regenerated Temperature Swing Adsorption; CubeSat Deployable Log Periodic Dipole Array; Re-entry Vehicle Shape for Enhanced Performance; NanoRacks-Scale MEMS Gas Chromatograph System; Variable Camber Aerodynamic Control Surfaces and Active Wing Shaping Control; Spacecraft Line-of-Sight Stabilization Using LWIR Earth Signature; Technique for Finding Retro-Reflectors in Flash LIDAR Imagery; Novel Hemispherical Dynamic Camera for EVAs; 360 deg Visual Detection and Object Tracking on an Autonomous Surface Vehicle; Simulation of Charge Carrier Mobility in Conducting Polymers; Observational Data Formatter Using CMOR for CMIP5; Propellant Loading Physics Model for Fault Detection Isolation and Recovery; Probabilistic Guidance for Swarms of Autonomous Agents; Reducing Drift in Stereo Visual Odometry; Future Air-Traffic Management Concepts Evaluation Tool; Examination and A Priori Analysis of a Direct Numerical Simulation Database for High-Pressure Turbulent Flows; and Resource-Constrained Application of Support Vector Machines to Imagery.

Source record↗
NASA NTRS2001

Single and Double ITCZ in Aqua-Planet Models with Globally Uniform Sea Surface Temperature and Solar Insolation: An Interpretation

It has been known for more than a decade that an aqua-planet model with globally uniform sea surface temperature and solar insolation angle can generate ITCZ (intertropical convergence zone). Previous studies have shown that the ITCZ under such model settings can be changed between a single ITCZ over the equator and a double ITCZ straddling the equator through one of several measures. These measures include switching to a different cumulus parameterization scheme, changes within the cumulus parameterization scheme, and changes in other aspects of the model design such as horizontal resolution. In this paper an interpretation for these findings is offered. The latitudinal location of the ITCZ is the latitude where the balance of two types of attraction on the ITCZ, both due to earth's rotation, exists. The first type is equator-ward and is directly related to the earth's rotation and thus not sensitive to model design changes. The second type is poleward and is related to the convective circulation and thus is sensitive to model design changes. Due to the shape of the attractors, the balance of the two types of attractions is reached either at the equator or more than 10 degrees away from the equator. The former case results in a single ITCZ over the equator and the latter case a double ITCZ straddling the equator.

Chao, Winston C.↗
NASA NTRS1999

Data Structure and Parallel Decomposition Considerations on a Fibonacci Grid

The Fibonacci grid, proposed by Swinbank and Purser (see companion abstract), provides attractive properties for global numerical atmospheric prediction by offering an optimally homogeneous, geometrically regular, and approximately isotropic discretization, with only the polar regions requiring special numerical treatment. It is a mathematical idealization, applied to the sphere, of the multi-spiral patterns often found in botanical structures, such as in pine cones and sunflower heads. Computationally, it is natural to organize the domain, into zones, in each of which the same pair, or triple, of "Fibonacci spirals" dominate. But the further subdivision of such zones into "tiles" of a shape and size suitable for distribution to the processors of a massively parallel computer requires very careful consideration if the subsequent spatial computations along the respective spirals, especially those computations (such as compact differencing schemes) that involve recursion, can be implemented in an efficient "load-balanced "manner without requiring excessive amounts of inter-processor communications. In this paper we show how certain "number theoretic" properties of the Fibonacci sequence (whose numbers prescribe the multiplicity of successive spirals) may be exploited in the decomposition of grid zones into tidy arrangements of triangular grid tiles, each tile possessing one side approximately parallel to the constant-latitude zone boundary. We also describe how the spatially recursive processes may be decomposed across such a tiling, and the directionality of the recursions reversed on alternate grid lines, to ensure a very high degree of load balancing throughout the execution of the computations required for one time step of a global model.

Michalakes, John↗
NASA NTRS2014

Range-Specific High-Resolution Mesoscale Model Setup: Data Assimilation

Mesoscale weather conditions can have an adverse effect on space launch, landing, and ground processing at the Eastern Range (ER) in Florida and Wallops Flight Facility (WFF) in Virginia. During summer, land-sea interactions across Kennedy Space Center (KSC) and Cape Canaveral Air Force Station (CCAFS) lead to sea breeze front formation, which can spawn deep convection that can hinder operations and endanger personnel and resources. Many other weak locally driven low-level boundaries and their interactions with the sea breeze front and each other can also initiate deep convection in the KSC/CCAFS area. Some of these other boundaries include the Indian River breeze front, Banana River breeze front, outflows from previous convection, horizontal convective rolls, convergence lines from other inland bodies of water such as Lake Okeechobee, the trailing convergence line from convergence of sea breeze fronts due to the shape of Cape Canaveral, frictional convergence lines from the islands in the Bahamas, convergence lines from soil moisture differences, convergence lines from cloud shading, and others. All these subtle weak boundary interactions often make forecasting of operationally important weather very difficult at KSC/CCAFS during the convective season (May-Oct). These convective processes often build quickly, last a short time (60 minutes or less), and occur over small distances, all of which also poses a significant challenge to the local forecasters who are responsible for issuing weather advisories, watches, and warnings. Surface winds during the transition seasons of spring and fall pose the most difficulties for the forecasters at WFF. They also encounter problems forecasting convective activity and temperature during those seasons. Therefore, accurate mesoscale model forecasts are needed to aid in their decision making. Both the ER and WFF would benefit greatly from high-resolution mesoscale model output to better forecast a variety of unique weather phenomena. Global and national scale models cannot properly resolve important local-scale weather features at each location due to their horizontal resolutions being much too coarse. Therefore, a properly tuned model at a high resolution is needed to provide improved capability. This task is a multi-year effort in which the Applied Meteorology Unit (AMU) will tune the Weather Research and Forecasting (WRF) model individually for each range. The goal of the first year, the results of which are in this report, was to tune the WRF model based on the best model resolution and run time while using reasonable computing capabilities. To accomplish this, the ER and WFF supported the tasking of the AMU to perform a number of sensitivity tests in order to determine the best model configuration for operational use at each of the ranges to best predict winds, precipitation, and temperature (Watson 2013). This task is a continuation of that work and will provide a recommended local data assimilation (DA) and numerical forecast model design optimized for the ER and WFF to support space launch activities. The model will be optimized for local weather challenges at both ranges.

Global and national scale models↗
NASA NTRS2014

An Approach for the Long-Term 30-m Land Surface Snow-Free Albedo Retrieval from Historic Landsat Surface Reflectance and MODIS-based A Priori Anisotropy Knowledge

Land surface albedo has been recognized by the Global Terrestrial Observing System (GTOS) as an essential climate variable crucial for accurate modeling and monitoring of the Earth's radiative budget. While global climate studies can leverage albedo datasets from MODIS, VIIRS, and other coarse-resolution sensors, many applications in heterogeneous environments can benefit from higher-resolution albedo products derived from Landsat. We previously developed a "MODIS-concurrent" approach for the 30-meter albedo estimation which relied on combining post-2000 Landsat data with MODIS Bidirectional Reflectance Distribution Function (BRDF) information. Here we present a "pre-MODIS era" approach to extend 30-m surface albedo generation in time back to the 1980s, through an a priori anisotropy Look-Up Table (LUT) built up from the high quality MCD43A BRDF estimates over representative homogenous regions. Each entry in the LUT reflects a unique combination of land cover, seasonality, terrain information, disturbance age and type, and Landsat optical spectral bands. An initial conceptual LUT was created for the Pacific Northwest (PNW) of the United States and provides BRDF shapes estimated from MODIS observations for undisturbed and disturbed surface types (including recovery trajectories of burned areas and non-fire disturbances). By accepting the assumption of a generally invariant BRDF shape for similar land surface structures as a priori information, spectral white-sky and black-sky albedos are derived through albedo-to-nadir reflectance ratios as a bridge between the Landsat and MODIS scale. A further narrow-to-broadband conversion based on radiative transfer simulations is adopted to produce broadband albedos at visible, near infrared, and shortwave regimes.We evaluate the accuracy of resultant Landsat albedo using available field measurements at forested AmeriFlux stations in the PNW region, and examine the consistency of the surface albedo generated by this approach respectively with that from the "concurrent" approach and the coincident MODIS operational surface albedo products. Using the tower measurements as reference, the derived Landsat 30-m snow-free shortwave broadband albedo yields an absolute accuracy of 0.02 with a root mean square error less than 0.016 and a bias of no more than 0.007. A further cross-comparison over individual scenes shows that the retrieved white sky shortwave albedo from the "pre-MODIS era" LUT approach is highly consistent (R(exp 2) = 0.988, the scene-averaged low RMSE = 0.009 and bias = −0.005) with that generated by the earlier "concurrent" approach. The Landsat albedo also exhibits more detailed landscape texture and a wider dynamic range of albedo values than the coincident 500-m MODIS operational products (MCD43A3), especially in the heterogeneous regions. Collectively, the "pre-MODIS" LUT and "concurrent" approaches provide a practical way to retrieve long-term Landsat albedo from the historic Landsat archives as far back as the 1980s, as well as the current Landsat-8 mission, and thus support investigations into the evolution of the albedo of terrestrial biomes at fine resolution.

Albedo algorithm↗
DOE OSTI2026

Evolutionary constraints and climate variability jointly shape starch–sugar balance in woody plants

Nonstructural carbohydrates (NSC) buffer plants against carbon imbalances, yet their partitioning between storage and soluble pools remains elusive at global scales. Here, we compiled a dataset of starch to soluble sugar ratio (St : Su) for 308 woody species across 220 sites world-wide and introduce a dimensionless index that integrates storage and demand while minimizing methodological artifacts. St : Su was strongly associated with growth, identifying it as a key axis of carbon allocation. Foliage consistently exhibited lower St : Su than lignified organs, reflecting a division between transient and conservative pools. Conifers accumulated more starch in foliage but less in stems relative to angiosperms, while leaf habits and mycorrhizal associations further modulated organ-specific strategies. Contrary to expectation, foliar and root St : Su varied little among biomes, but stems exhibited higher ratios in tropical rainforests than in boreal or arid regions, reflecting differences in species composition and adaptive storage under disturbance. Phylogeny constrained stem storage, whereas climatic variability, rather than mean conditions, dominated allocation in leaves and roots. These findings establish St : Su as a robust functional trait linking allocation strategies, growth, and resilience, which can be used to improve vegetation model prediction of forest productivity and mortality under climate variability.

Li, Weibin [Lanzhou Univ. (China)] (ORCID:00000001↗