Search NASA⌕ Search

SEARCH · Search NASA

Results for “spatially explicit 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 127 records · Page 7

Magnetic field of Jupiter: A generalized inverse approach

The estimation of planetary magnetic fields from observations of the magnetic field gathered along a spacecraft flyby trajectory is examined with the aid of generalized inverse techniques, with application to the internal magnetic field of Jupiter. Model nonuniqueness resulting from the limited spatial extent of the observations and noise on the data is explored and quantitative estimates of the model parameter resolution are found. The presence of a substantial magnetic field of external origin due to the currents flowing in the Jovian magnetodisc is found to be an important source of error in estimates of the internal Jovian field, and new models explicitly incorporating these currents are proposed. New internal field models are derived using the vector helium magnetometer observations and the high field fluxgate observations of Pioneer 11, and knowledge of the external current system gained from the Pioneer 10 and Voyagers 1 and 2 encounters.

Connerney, J. E. P.↗

The magnetic field of Jupiter - A generalized inverse approach

The estimation of planetary magnetic fields from observations of the magnetic field gathered along a spacecraft flyby trajectory is examined with the aid of generalized inverse techniques, with application to the internal magnetic field of Jupiter. Model non-uniqueness resulting from the limited spatial extent of the observations and noise on the data is explored and quantitative estimates of the model parameter resolution are found. The presence of a substantial magnetic field of external origin due to the currents flowing in the Jovian magnetodisc is found to be an important source of error in estimates of the internal Jovian field, and new models explicitly incorporating these currents are proposed. New internal field models are derived using the vector helium magnetometer observations and the high field fluxgate observations of Pioneer 11, and knowledge of the external current system gained from the Pioneer 10 and Voyagers 1 and 2 encounters.

Connerney, J. E. P.↗

A Thermodynamic Paradigm for Studying Disease Vector’s Habitats & Life Cycles Using NASA’s NextGen Remote Sensing Instruments

Global public health is entering a new informational age through the use of spatial models of disease vector/host ecologies driven by the use of remotely sensed data to measure environmental and structural factors critical in determining disease vector habitats, distributions, life cycles, and host interactions. The vector habitat microclimates can be quantified in terms of the surface energy budget measured by satellites. The epidemiological equations (processes) can be adapted and modified to explicitly incorporate environmental factors and interfaces required by a specific disease and its vector/host cycle. Remote sensing can be used to measure or evaluate or estimate both environment (state functions) and interface (process functions). Remote sensing can be used to measure or evaluate or estimate both environment (state functions) and interface (process functions) defining vector habitats. The products of remote sensing can be integrated directly into the epidemiological equations to significantly enhance our understanding of disease vector’s life cycles and habitats. The next generation of NASA’s remote sensing instruments recently operational will provide a significant enhancement in our ability to study disease vector’s life cycles and habitats. Three new instruments have become operational on the International Space Station (ISS) In June 2018, ECOSTRESS (The ECOsystem Spaceborne Thermal Radiometer Experiment on Space Station), a five channel, thermal IR instrument with 70 m resolution and approximately 4-5 day repeat cycle of day/night pairs along with DESISa hyperspectral sensor system 235 channels (between 400 and 1000 nanometers), 30 m resolution, and data only being acquired on demand. In November 2018. GEDI (Global Ecosystem Dynamics Investigation). GEDI is a high-resolution laser ranging of Earth’s forests and topography. These instruments set the stage for NASA’s new global designated observables set of remote sensing measurements collecting hyperspectral and multispectral thermal data called SBG (Surface Biology Geology). The SBG is currently in the planning stage. The availability and use of level 1-4 data products generated from these ISS instruments will be discussed.

thermodynamics↗

Influence of Sea-Ice Anomalies on Antarctic Precipitation Using Source Attribution in the Community Earth System Model

We conduct sensitivity experiments using a general circulation model that has an explicit water source tagging capability forced by prescribed composites of pre-industrial sea-ice concentrations (SICs) and corresponding sea surface temperatures (SSTs) to understand the impact of sea-ice anomalies on regional evaporation, moisture transport and source–receptor relationships for Antarctic precipitation in the absence of anthropogenic forcing. Surface sensible heat fluxes, evaporation and column-integrated water vapor are larger over Southern Ocean (SO) areas with lower SICs. Changes in Antarctic precipitation and its source attribution with SICs have a strong spatial variability. Among the tagged source regions, the Southern Ocean (south of 50∘ S) contributes the most (40 %) to the Antarctic total precipitation, followed by more northerly ocean basins, most notably the South Pacific Ocean (27%), southern Indian Ocean (16 %) and South Atlantic Ocean (11 %). Comparing two experiments prescribed with high and low pre-industrial SICs, respectively, the annual mean Antarctic precipitation is about 150 Gt yr−1 (or 6 %) more in the lower SIC case than in the higher SIC case. This difference is larger than the model-simulated interannual variability in Antarctic precipitation (99 Gt yr−1). The contrast in contribution from the Southern Ocean, 102 Gt yr−1, is even more significant compared to the interannual variability of 35 Gt yr−1 in Antarctic precipitation that originates from the Southern Ocean. The horizontal transport pathways from individual vapor source regions to Antarctica are largely determined by large-scale atmospheric circulation patterns. Vapor from lower-latitude source regions takes elevated pathways to Antarctica. In contrast, vapor from the Southern Ocean moves southward within the lower troposphere to the Antarctic continent along moist isentropes that are largely shaped by local ambient conditions and coastal topography. This study also highlights the importance of atmospheric dynamics in affecting the thermodynamic impact of sea-ice anomalies associated with natural variability on Antarctic precipitation. Our analyses of the seasonal contrast in changes of basin-scale evaporation, moisture flux and precipitation suggest that the impact of SIC anomalies on regional Antarctic precipitation depends on dynamic changes that arise from SIC–SST perturbations along with internal variability. The latter appears to have a more significant effect on the moisture transport in austral winter than in summer.

Wang, Hailong↗

Nonlinear vortex trail dynamics. II - Analytic solutions

Spatially periodic large amplitude solutions of the von Karman model are obtained in the neighborhood of singularities. These singularities correspond to vortex clusters in the physical plane. The quasi-periodic and unbounded solutions found analytically confirm earlier numerical work and show qualitative agreement with experimental observations of large-scale phenomena of vortex trails. Separatrices or heteroclinic orbits were explicitly found for an integrable approximate equation, which indicate that the von Karman model itself supports chaotic solutions.

Lim, Chjan C.↗

Incorporation of spatial information in Bayesian image reconstruction - The maximum residual likelihood criterion

We have developed a new figure of merit, a 'maximum-residual-likelihood' (MRL) statistic, for the goodness of fit for Bayesian image restoration which explicitly incorporates spatial information. The MRL constraint provides a natural means of incorporating the prior knowledge that the residuals contain no spatial structure through the autocorrelation function of the residuals. We demonstrate that this statistic follows a Chi-square distribution and that forcing this statistic to have its most probable value leads to a restored image whose residuals are consistent with the noise model. Our numerical experiments suggest that image restoration using the MRL statistic alone is numerically robust and produces results which are independent of the initial guess for the restored image. However, we caution that using the MRL statistic without an image prior can result in overresolution in low SNR portions of the image.

Pina, R. K.↗

Sensitivity analysis for the dynamic response of viscoplastic shells of revolution

A computational procedure is presented for evaluating the sensitivity coefficients of the dynamic axisymmetric response of viscoplastic shells of revolution. The analytical formulation is based on Reissner's large deformation shell theory with the effects of transverse shear deformation, rotatory inertia and moments turning around the normal to the middle surface included. The material model is chosen to be isothermal viscoplasticity, and an associated flow rule is used with a von Mises effective stress. A mixed formulation is used with the fundamental unknowns consisting of six stress resultants, three generalized displacements and three velocity components. Spatial discretization is performed using finite elements, with discontinuous stress resultants across element interfaces. The temporal integration is performed by using an explicit central difference scheme (leap-frog method) with an implicit constitutive update. The sensitivity coefficients are evaluated using a direct differentiation approach. Numerical results are presented for a spherical cap subjected to step loading, and a circular plate subjected to impulsive loading. The sensitivity coefficients are generated by evaluating the derivatives of the response quantities with respect to thickness, mass density, Young's modulus, and two of the material parameters characterizing the viscoplastic response. Time histories of the response and sensitivity coefficients are presented, along with spatial distributions of these quantities at selected times.

Kulkarni, Makarand↗

Titan Polar Landscape Evolution

With the ongoing Cassini-era observations and studies of Titan it is clear that the intensity and distribution of surface processes (particularly fluvial erosion by methane and Aeolian transport) has changed through time. Currently however, alternate hypotheses substantially differ among specific scenarios with respect to the effects of atmospheric evolution, seasonal changes, and endogenic processes. We have studied the evolution of Titan's polar region through a combination of analysis of imaging, elevation data, and geomorphic mapping, spatially explicit simulations of landform evolution, and quantitative comparison of the simulated landscapes with corresponding Titan morphology. We have quantitatively evaluated alternate scenarios for the landform evolution of Titan's polar terrain. The investigations have been guided by recent geomorphic mapping and topographic characterization of the polar regions that are used to frame hypotheses of process interactions, which have been evaluated using simulation modeling. Topographic information about Titan's polar region is be based on SAR-Topography and altimetry archived on PDS, SAR-based stereo radar-grammetry, radar-sounding lake depth measurements, and superposition relationships between geomorphologic map units, which we will use to create a generalized topographic map.

Cassini↗

Efficient stream distributions in radiative transfer theory

This paper discusses a new, computationally-efficient method for approximating the integro-differential equation of radiative transfer with a finite set of coupled differential equations for discrete streams. The method uses recommended spatial distributions of streams that are quite different from those typically used in that they are based on the symmetry of several regular Platonic solids. To facilitate the use of such distributions, an explicit, one-parameter relationship between the physical radiance and the abstract stream is formulated. The parameter is used to determine the minimum number of streams required in the radiative transfer model. Accuracy and computational efficiency are shown to be served best by choosing a stream distribution that is invariant to a large number of three space rotations. For various values of the above-mentioned parameter, the resulting recommended stream distribution is shown to be more computationally efficient than more conventional stream distributions. Finally, the incorporation of polarization in the stream definition is described.

Whitney, C.↗

Generation of Continental Scale Percent Tree Cover Product Using Deep-learning and Multi-scale Remote Sensing Data

Spatially explicit percent tree cover (TC) estimation is critical for mapping forest aboveground biomass and its dynamics. While various TC products have been developed, there has not been a generalized framework that can be applied to diverse terrestrial ecosystems due to underlain extreme complexities. Deep learning algorithms can learn a spatial pattern and radiometric characteristics of tree canopy as a robust approximation of physical or empirical models, and thus have emerged as promising and efficient tools for large-scale TC mapping. In this study, we synergistically use very high-resolution aerial imageries (National Agriculture Imagery Program, NAIP) and medium resolution Landsat data to map continental-scale TC (CONUS and Mexico) through a hierarchical deep learning approach (Convolutional Neural Network), i.e., NAIP TC generated from a NAIP model is utilized to train a Landsat model. The produced TC product (hereafter, NEX-TC) is able to capture the spatial pattern of TC distribution and its changes driven by natural disturbance and human land management. We further explore and analyze the reliability and potential uncertainty of the NEX-TC by comparing it to lidar- (lidar-TC), National Land Cover Database (NLCD-TC), and MODIS Vegetation Continuous Field (MODIS-TC). This evaluation practice reveals that TC products based on passive optical sensors tend to underestimate TC across all land cover types while Landsat-based TCs (i.e., NEX-TC & NLCD-TC) perform better than the coarser MODIS TC estimate. Our results show that the NEX-TC is generally comparable to NLCD-TC but it particularly outperforms NLCD-TC and MODIS-TC over the dense forests where lidar-TC indicates >80% TC. These results indicate that our hierarchical deep learning approach and TC product will be effective and useful for characterizing large-scale tree cover and possibly associated carbon dynamics.

Landsat↗

Sensitivity analysis of the non-linear dynamic viscoplastic response of 2-d structures with respect to material parameters

A computational procedure is presented for evaluating the sensitivity coefficients of the viscoplastic response of structures subjected to dynamic loading. A state of plane stress is assumed to exist in the structure, a velocity strain-Cauchy stress formulation is used, and the geometric non-linearities arising from large strains are incorporated. The Jaumann rate is used as a frame indifferent stress rate. The material model is chosen to be isothermal viscoplasticity, and an associated flow rule is used with a von Mises effective stress. The equations of motion emanating from a finite element semi-discretization are integrated using an explicit central difference scheme with an implicit stress update. The sensitivity coefficients are evaluated using a direct differentiation approach. Since the domain of integration is the current configuration, the sensitivity coefficients of the spatial derivatives of the shape functions must be included. Numerical results are presented for a thin plate with a central cutout subjected to an in-plane compressive loading. The sensitivity coefficients are generated by evaluating the derivatives of the response quantities with respect to Young's modulus, and two of the material parameters characterizing the viscoplastic response. Time histories of the response and sensitivity coefficients, and spatial distributions at selected times are presented.

Kulkarni, Makarand↗

Supplement of: The Influence of Volcanic Eruptions on the Climate of Tropical South America During the Last Millennium in an Isotope-Enabled General Circulation Model

Currently, little is known on how volcanic eruptions impact large-scale climate phenomena such as South American paleo-intertropical Convergence Zone (ITCZ) position and summer monsoon behavior. In this paper, an analysis of observations and model simulations is employed to assess the influence of large volcanic eruptions on the climate of tropical South America. This problem is first considered for historically recent volcanic episodes for which more observations are available but where fewer events exist and the confounding effects of El NioSouthern Oscillation (ENSO) lead to inconclusive interpretation of the impact of volcanic eruptions at the continental scale. Therefore, we also examine a greater number of reconstructed volcanic events for the period 850CE to present that are incorporated into the NASA GISS ModelE2-R simulation of the last millennium.An advantage of this model is its ability to explicitly track water isotopologues throughout the hydrologic cycle and simulating the isotopic imprint following a large eruption. This effectively removes a degree of uncertainty associated with error-prone conversion of isotopic signals into climate variables, and allows for a direct comparison between GISS simulations and paleoclimate proxy records.Our analysis reveals that both precipitation and oxygen isotope variability respond with a distinct seasonal and spatial structure across tropical South America following an eruption. During austral winter, the heavy oxygen isotope in precipitation is enriched, likely due to reduced moisture convergence in the ITCZ domain and reduced rainfall over northern South America. During austral summer, however, more negative values of the precipitation isotopic composition are simulated over Amazonia, despite reductions in rainfall, suggesting that the isotopic response is not a simple function of the amount effect. During the South American monsoon season, the amplitude of the temperature response to volcanic forcing is larger than the rather weak and spatially less coherent precipitation signal, complicating the isotopic response to changes in the hydrologic cycle.

Colose, Christopher↗

The Influence of Volcanic Eruptions on the Climate of Tropical South America During the Last Millennium in an Isotope-Enabled General Circulation Model

Currently, little is known on how volcanic eruptions impact large-scale climate phenomena such as South American paleo-intertropical Convergence Zone (ITCZ) position and summer monsoon behavior. In this paper, an analysis of observations and model simulations is employed to assess the influence of large volcanic eruptions on the climate of tropical South America. This problem is first considered for historically recent volcanic episodes for which more observations are available but where fewer events exist and the confounding effects of El Niño-Southern Oscillation (ENSO) lead to inconclusive interpretation of the impact of volcanic eruptions at the continental scale. Therefore, we also examine a greater number of reconstructed volcanic events for the period 850 CE to present that are incorporated into the NASA GISS ModelE2-R simulation of the last millennium. An advantage of this model is its ability to explicitly track water isotopologues throughout the hydrologic cycle and simulating the isotopic imprint following a large eruption. This effectively removes a degree of uncertainty associated with error-prone conversion of isotopic signals into climate variables, and allows for a direct comparison between GISS simulations and paleoclimate proxy records. Our analysis reveals that both precipitation and oxygen isotope variability respond with a distinct seasonal and spatial structure across tropical South America following an eruption. During austral winter, the heavy oxygen isotope in precipitation is enriched, likely due to reduced moisture convergence in the ITCZ domain and reduced rainfall over northern South America. During austral summer, however, more negative values of the precipitation isotopic composition are simulated over Amazonia, despite reductions in rainfall, suggesting that the isotopic response is not a simple function of the "amount effect". During the South American monsoon season, the amplitude of the temperature response to volcanic forcing is larger than the rather weak and spatially less coherent precipitation signal, complicating the isotopic response to changes in the hydrologic cycle.

Colose, Christopher M.↗

Stability analysis of numerical boundary conditions and implicit difference approximations for hyperbolic equations

Implicit, noniterative, finite difference schemes were recently developed by several authors for multidimensional systems of nonlinear hyperbolic partial differential equations. When applied to linear model equations with periodic boundary conditions those schemes are unconditionally stable (A-stable). As applied in practice the algorithms often face a severe time step restriction. A major source of the difficulty is the treatment of the numerical boundary conditions. One conjecture was that unconditional stability requires implicit numerical boundary conditions. An apparent counter example was the space time extrapolation considered by Gustafsson, Kreiss, and Sunstrom. Spatial (implicit) and space time (explicit) extrapolation using normal mode analysis for a finite and infinite number of spatial mesh intervals are examined. The results indicate that for unconditional stability with a finite number of spatial mesh intervals, the numerical boundary conditions must be implicit.

Beam, R. M.↗

Simple Analytic Expressions for the Magnetic Field of a Circular Current Loop

Analytic expressions for the magnetic induction (magnetic flux density, B) of a simple planar circular current loop have been published in Cartesian and cylindrical coordinates [1,2], and are also known implicitly in spherical coordinates [3]. In this paper, we present explicit analytic expressions for B and its spatial derivatives in Cartesian, cylindrical, and spherical coordinates for a filamentary current loop. These results were obtained with extensive use of Mathematica "TM" and are exact throughout all space outside of the conductor. The field expressions reduce to the well-known limiting cases and satisfy V · B = 0 and V x B = 0 outside the conductor. These results are general and applicable to any model using filamentary circular current loops. Solenoids of arbitrary size may be easily modeled by approximating the total magnetic induction as the sum of those for the individual loops. The inclusion of the spatial derivatives expands their utility to magnetohydrodynamics where the derivatives are required. The equations can be coded into any high-level programming language. It is necessary to numerically evaluate complete elliptic integrals of the first and second kind, but this capability is now available with most programming packages.

Simpson, James C.↗

WRLES: Wave Resolving Large-Eddy Simulation Code, Theory and Usage

A computational fluid dynamics code has been developed for large-eddy simulations (LES) of turbulent flow. The code uses high-order of accuracy and high-resolution numerical methods to minimize solution error and maximize the resolution of the turbulent structures. Spatial discretization is performed using explicit central differencing. The central differencing schemes in the code include 2nd- to 12th-order standard central difference methods as well as 7-, 9-, 11- and 13-point dispersion relation preserving schemes. Solution filtering and high-order shock capturing are included for stability. Time discretization is performed using multistage Runge-Kutta methods that are up to 4th order accurate. Several options are available to model turbulence including: Baldwin-Lomax and Spalart-Allmaras Reynolds-averaged Navier-Stokes turbulence models, and Smagorinsky, Dynamic Smagorinsky and Vreman sub-grid scale models for LES. This report presents the theory behind the numerical and physical models used in the code and provides a user's manual to the operation of the code.

High-Order↗

Evidence for Tropopause Layer Moistening by Convection During CRYSTAL-FACE

Measurements and analysis of the impact of deep convection on tropopause layer moisture are easily confounded by difficulties making precise observations with sufficient spatial coverage before and after convective events and difficulties distinguishing between changes due to local convection versus large-scale advection. The interactions between cloud microphysics and dynamics in the convective transport of moisture into the tropopause layer also result in a sufficiently complex and poorly characterized system to allow for considerable freedom in theoretical models of stratosphere-troposphere exchange. In this work we perform detailed large-eddy simulations with an explicit cloud microphysics model to study the impact of deep convection on tropopause layer moisture profiles observed over southern Florida during CRYSTALFACE. For four days during the campaign (July 11, 16, 28, and 29) we initialize a 100-km square domain with temperature and moisture profiles measured prior to convection at the PARSL ground site, and initiate convection with a warm bubble that produces an anvil at peak elevations in agreement with lidar and radar observations on that day. Comparing the moisture field after the anvils decay with the initial state, we find that convection predominantly moistens the tropopause layer (as defined by minimum temperature and minimum potential temperature lapse rate), although some drying is also predicted in localized layers. We will also present results of sensitivity tests designed to separate the roles of cloud microphysics and dynamics.

Ackerman, A.↗

A Comparison Between Gravity Wave Momentum Fluxes in Observations and Climate Models

For the first time, a formal comparison is made between gravity wave momentum fluxes in models and those derived from observations. Although gravity waves occur over a wide range of spatial and temporal scales, the focus of this paper is on scales that are being parameterized in present climate models, sub-1000-km scales. Only observational methods that permit derivation of gravity wave momentum fluxes over large geographical areas are discussed, and these are from satellite temperature measurements, constant-density long-duration balloons, and high-vertical-resolution radiosonde data. The models discussed include two high-resolution models in which gravity waves are explicitly modeled, Kanto and the Community Atmosphere Model, version 5 (CAM5), and three climate models containing gravity wave parameterizations,MAECHAM5, Hadley Centre Global Environmental Model 3 (HadGEM3), and the Goddard Institute for Space Studies (GISS) model. Measurements generally show similar flux magnitudes as in models, except that the fluxes derived from satellite measurements fall off more rapidly with height. This is likely due to limitations on the observable range of wavelengths, although other factors may contribute. When one accounts for this more rapid fall off, the geographical distribution of the fluxes from observations and models compare reasonably well, except for certain features that depend on the specification of the nonorographic gravity wave source functions in the climate models. For instance, both the observed fluxes and those in the high-resolution models are very small at summer high latitudes, but this is not the case for some of the climate models. This comparison between gravity wave fluxes from climate models, high-resolution models, and fluxes derived from observations indicates that such efforts offer a promising path toward improving specifications of gravity wave sources in climate models.

mathematical models↗